Effect of coarse graining in water models for the study of kinetics and mechanisms of clathrate hydrates nucleation and growth.

Clathrate hydrates are crystalline inclusion compounds wherein a water framework encages small guest atoms/molecules within its cavities. Among the others, methane clathrates are the largest fossil fuel resource still available. They can also be used to safely transport gases and can also form spontaneously under suitable conditions plugging pipelines.

Jarzynski equality on work and free energy: Crystal indentation as a case study.

Mathematical relations concerning particle systems require knowledge of the applicability conditions to become physically relevant and not merely formal. We illustrate this fact through the analysis of the Jarzynski equality (JE), whose derivation for Hamiltonian systems suggests that the equilibrium free-energy variations can be computational or experimentally determined in almost any kind of non-equilibrium processes. This apparent generality is surprising in a mechanical theory.

Dynamical nonequilibrium molecular dynamics reveals the structural basis for allostery and signal propagation in biomolecular systems.

Abstract: A dynamical approach to nonequilibrium molecular dynamics (D-NEMD), proposed in the 1970s by Ciccotti et al., is undergoing a renaissance and is having increasing impact in the study of biological macromolecules. This D-NEMD approach, combining MD simulations in stationary (in particular, equilibrium) and nonequilibrium conditions, allows for the determination of the time-dependent structural response of a system using the Kubo-Onsager relation.

Quantum Trajectories for the Dynamics in the Exact Factorization Framework: A Proof-of-Principle Test.

In the framework of the exact factorization of the time-dependent electron-nuclear wave function, we investigate the possibility of solving the nuclear time-dependent Schrödinger equation based on trajectories. The nuclear equation is separated in a Hamilton-Jacobi equation for the phase of the wave function, and a continuity equation for its (squared) modulus. For illustrative adiabatic and nonadiabatic one-dimensional models, we implement a procedure to follow the evolution of the nuclear density along the characteristics of the Hamilton-Jacobi equation.

Adiabatic motion and statistical mechanics via mass-zero constrained dynamics.

In recent work [Coretti et al., J. Chem.

A General Mechanism for Signal Propagation in the Nicotinic Acetylcholine Receptor Family.

Nicotinic acetylcholine receptors (nAChRs) modulate synaptic activity in the central nervous system. The α7 subtype, in particular, has attracted considerable interest in drug discovery as a target for several conditions, including Alzheimer's disease and schizophrenia. Identifying agonist-induced structural changes underlying nAChR activation is fundamentally important for understanding biological function and rational drug design.

Ab initio accelerated molecular dynamics study of the hydride ligands in the ruthenium complex: Ru(H)H(P(CH)).

The dihydrogen complex Ru(H)H(P(CH)) has been investigated, via ab initio accelerated molecular dynamics, to elucidate the H ligands dynamics and possible reaction paths for H/H exchange. We have characterized the free energy landscape associated with the H atoms positional exchange around the Ru centre. From the free energy landscape, we have been able to estimate a barrier of 6 kcal mol for the H/H exchange process.

On computing the solubility of molecular systems subject to constraints using the extended Einstein crystal method.

A method to compute solubilities for molecular systems using atomistic simulations, based on an extension of the Einstein crystal method, has recently been presented [Li et al., J. Chem.

Solubility of paracetamol in ethanol by molecular dynamics using the extended Einstein crystal method and experiments.

Li and co-workers [Li et al., J. Chem.

The physics of open systems for the simulation of complex molecular environments in soft matter.

Molecular dynamics (MD) has become one of the most powerful tools of investigation in soft matter. Despite such success, simulations of large molecular environments are mostly run using the approximation of closed systems without the possibility of exchange of matter. Due to the molecular complexity of soft matter systems, an optimal simulation strategy would require the application of concurrent multiscale resolution approaches such that each part of a large system can be considered as an open subsystem at a high resolution embedded in a large coarser reservoir of energy and particles.

Communication: Constrained molecular dynamics for polarizable models.

A new algorithm to solve numerically the evolution of empirical shell models of polarizable systems is presented. It employs constrained molecular dynamics to satisfy exactly, at each time step, the crucial condition that the gradient of the potential with respect to the shell degrees of freedom is null. The algorithm is efficient, stable, and, contrary to the available alternatives, it is symplectic and time reversible.

Time-reversal symmetry for systems in a constant external magnetic field.

The time-reversal properties of charged systems in a constant external magnetic field are reconsidered in this paper. We show that the evolution equations of the system are invariant under a new symmetry operation that implies a new signature property for time-correlation functions under time reversal. We then show how these findings can be combined with a previously identified symmetry to determine, for example, null components of the correlation functions of velocities and currents and of the associated transport coefficients.

Nucleation of Molecular Crystals Driven by Relative Information Entropy.

Simulating nucleation of molecular crystals is extremely challenging for all but the simplest cases. The challenge lies in formulating effective order parameters that are capable of driving the transition process. In recent years, order parameters based on molecular pair-functions have been successfully used in combination with enhanced sampling techniques to simulate nucleation of simple molecular crystals.

Thermal Diffusion in Binary Mixtures: Transient Behavior and Transport Coefficients from Equilibrium and Nonequilibrium Molecular Dynamics.

Equilibrium and nonequilibrium molecular dynamics simulations are combined to compute the full set of coefficients that appear in the phenomenological equations describing thermal transport in a binary mixture subject to a constant thermal gradient. The Dynamical Non-Equilibrium Molecular Dynamics approach (D-NEMD) is employed to obtain the microscopic time evolution of the density and temperature fields, together with that of the mass and energy fluxes. D-NEMD enables one to study not only the steady state, but also the evolution of the fields during the transient that follows the onset of the thermal gradient, up to the establishment of the steady state.

On metastability and Markov state models for non-stationary molecular dynamics.

Unlike for systems in equilibrium, a straightforward definition of a metastable set in the non-stationary, non-equilibrium case may only be given case-by-case-and therefore it is not directly useful any more, in particular in cases where the slowest relaxation time scales are comparable to the time scales at which the external field driving the system varies. We generalize the concept of metastability by relying on the theory of coherent sets. A pair of sets A and B is called coherent with respect to the time interval [t, t] if (a) most of the trajectories starting in A at t end up in B at t and (b) most of the trajectories arriving in B at t actually started from A at t.

Particle-Based Modeling of Living Actin Filaments in an Optical Trap.

We report a coarse-grained molecular dynamics simulation study of a bundle of parallel actin filaments under supercritical conditions pressing against a loaded mobile wall using a particle-based approach where each particle represents an actin unit. The filaments are grafted to a fixed wall at one end and are reactive at the other end, where they can perform single monomer (de)polymerization steps and push on a mobile obstacle. We simulate a reactive grand canonical ensemble in a box of fixed transverse area , with a fixed number of grafted filaments N f , at temperature and monomer chemical potential μ 1 .

On the properties of a bundle of flexible actin filaments in an optical trap.

We establish the statistical mechanics framework for a bundle of Nf living and uncrosslinked actin filaments in a supercritical solution of free monomers pressing against a mobile wall. The filaments are anchored normally to a fixed planar surface at one of their ends and, because of their limited flexibility, they grow almost parallel to each other. Their growing ends hit a moving obstacle, depicted as a second planar wall, parallel to the previous one and subjected to a harmonic compressive force.

Spontaneous chiral symmetry breaking in collective active motion.

Chiral symmetry breaking is ubiquitous in biological systems, from DNA to bacterial suspensions. A key unresolved problem is how chiral structures may spontaneously emerge from achiral interactions. We study a simple model of active swimmers in three dimensions that effectively incorporates hydrodynamic interactions.

Probabilistic Derivation of Spatiotemporal Correlation Functions in the Hydrodynamic Limit.

In this paper, we use probability theory to prove in suitable conditions the equivalence of equilibrium time correlation functions of microscopic density fields with the time correlation functions of local macroscopic density fields evolved by hydrodynamics in (approximate) phenomenological continuum theories of matter. We further discuss a useful and rigorous numerical algorithm, derived from this framework, to compute macroscopic space- and time-dependent behaviors (such as the hydrodynamical one) via molecular dynamics simulations.

A semi-flexible model prediction for the polymerization force exerted by a living F-actin filament on a fixed wall.

We consider a single living semi-flexible filament with persistence length ℓp in chemical equilibrium with a solution of free monomers at fixed monomer chemical potential μ1 and fixed temperature T. While one end of the filament is chemically active with single monomer (de)polymerization steps, the other end is grafted normally to a rigid wall to mimic a rigid network from which the filament under consideration emerges. A second rigid wall, parallel to the grafting wall, is fixed at distance L < < ℓp from the filament seed.

A Structural Model of the Human α7 Nicotinic Receptor in an Open Conformation.

Nicotinic acetylcholine receptors (nAchRs) are ligand-gated ion channels that regulate chemical transmission at the neuromuscular junction. Structural information is available at low resolution from open and closed forms of an eukaryotic receptor, and at high resolution from other members of the same structural family, two prokaryotic orthologs and an eukaryotic GluCl channel. Structures of human channels however are still lacking.

Clathrate structure-type recognition: Application to hydrate nucleation and crystallisation.

For clathrate-hydrate polymorphic structure-type (sI versus sII), geometric recognition criteria have been developed and validated. These are applied to the study of the rich interplay and development of both sI and sII motifs in a variety of hydrate-nucleation events for methane and H2S hydrate studied by direct and enhanced-sampling molecular dynamics (MD) simulations. In the case of nucleation of methane hydrate from enhanced-sampling simulation, we notice that already at the transition state, ∼80% of the enclathrated CH4 molecules are contained in a well-structured (sII) clathrate-like crystallite.

Maximum probability domains for the analysis of the microscopic structure of liquids.

We introduce the concept of maximum probability domains (MPDs), developed in the context of the analysis of electronic densities, in the study of the microscopic spatial structures of liquids. The idea of locating a particle in a three dimensional region, by determining the domain where the probability of finding that, and only that, particle is maximum, gives an interesting characterization of the local structure of the liquid. The optimization procedure, required for the search of the domain of maximum probability, is carried out by the implementation of the level set method.

Diagnostic ultrasound-guided excisional testicular biopsy for small (< 1 cm) incidental nodules. A single institution experience.

Introduction And Objectives: The widespread use of scrotal ultrasound (SUS) has led to a marked increase in the number of incidentally detected testicular lesions. A small incidental nodule (STN) has defined as a non palpable (< 10 mm), asymptomatic solid lesion with normal levels of oncological testicular markers. Nowadays the lack of agreement on the topic causes managing problems to andrologists.