Joseph A. Morrone

Nuclear quantum effects in water.

A path-integral Car-Parrinello molecular dynamics simulation of liquid water and ice is performed. It is found that the inclusion of nuclear quantum effects systematically improves the agreement of first-principles simulations of liquid water with experiment. In addition, the proton momentum distribution is computed utilizing a recently developed open path-integral molecular dynamics methodology. It is shown that these results are in good agreement with experimental data.

Tunneling and delocalization effects in hydrogen bonded systems: A study in position and momentum space

Novel experimental and computational studies have uncovered the proton momentum distribution in hydrogen bonded systems. In this work, we utilize recently developed open path integral Car-Parrinello molecular dynamics methodology in order to study the momentum distribution in phases of high pressure ice. Some of these phases exhibit symmetric hydrogen bonds and quantum tunneling. We find that the symmetric hydrogen bonded phase possesses a narrowed momentum distribution as compared with a covalently bonded phase, in agreement with recent experimental findings. The signatures of tunneling that we observe are a narrowed distribution in the low-to-intermediate momentum region, with a tail that extends to match the result of the covalently bonded state. The transition to tunneling behavior shows similarity to features observed in recent experiments performed on confined water. We corroborate our ice simulations with a study of a particle in a model one-dimensional double well potential that mimics some of the effects observed in bulk simulations. The temperature dependence of the momentum distribution in the one-dimensional model allows for the differentiation between ground state and mixed state tunneling effects.

Quantum fluctuations can promote or inhibit glass formation

Intuition suggests that the occurrence of large quantum fluctuations should prevent a material from forming a glass by enabling its atoms to rearrange into a lower-energy ordered state. But new simulations suggest the opposite could be true, with fluctuations sometimes enhancing glass formation.

Interplay between hydrodynamics and the free energy surface in the assembly of nanoscale hydrophobes.

Solvent plays an important role in the relative motion of nanoscopic bodies, and the study of such phenomena can help elucidate the mechanism of hydrophobic assembly, as well as the influence of solvent-mediated effects on in vivo motion in crowded cellular environments. Here we study important aspects of this problem within the framework of Brownian dynamics. We compute the free energy surface that the Brownian particles experience and their hydrodynamic interactions from molecular dynamics simulations in explicit solvent. We find that molecular scale effects dominate at short distances, thus giving rise to deviations from the predictions of continuum hydrodynamic theory. Drying phenomena, solvent layering, and fluctuations engender distinct signatures of the molecular scale. The rate of assembly in the diffusion-controlled limit is found to decrease from molecular scale hydrodynamic interactions, in opposition to the free energy driving force for hydrophobic assembly, and act to reinforce the influence of the free energy surface on the association of more hydrophilic bodies.

Molecular Dynamics with Multiple Time Scales: How to Avoid Pitfalls.

Multiple time scale methodologies have gained widespread use in molecular dynamics simulations and are implemented in a variety of ways across numerous packages. However, performance of the algorithms depends upon the details of the implementation. This is particularly important in the way in which the nonbonded interactions are partitioned. In this work, we show why some previous implementations give rise to energy drifts, and how this can be corrected. We also provide a recipe for using multiple time step methods to generate stable trajectories in large scale biomolecular simulations, where long trajectories are needed.

Displaced path integral formulation for the momentum distribution of quantum particles

The proton momentum distribution, accessible by deep inelastic neutron scattering, is a very sensitive probe of the potential of mean force experienced by the protons in hydrogen-bonded systems. In this work we introduce a novel estimator for the end-to-end distribution of the Feynman paths, i.e., the Fourier transform of the momentum distribution. In this formulation, free particle and environmental contributions factorize. Moreover, the environmental contribution has a natural analogy to a free energy surface in statistical mechanics, facilitating the interpretation of experiments. The new formulation is not only conceptually but also computationally advantageous. We illustrate the method with applications to an empirical water model, ab initio ice, and one dimensional model systems.

Are Hydrodynamic Interactions Important in the Kinetics of Hydrophobic Collapse

We study the kinetics of assembly of two plates of varying hydrophobicity, including cases where drying occurs and water strongly solvates the plate surfaces. The potential of mean force and molecular-scale hydrodynamics are computed from molecular dynamics simulations in explicit solvent as a function of particle separation. In agreement with our recent work on nanospheres [J. Phys. Chem. B 2012,116, 378-389], regions of high friction are found to be engendered by large and slow solvent fluctuations. These slow fluctuations can be due to either drying or confinement. The mean first passage times for assembly are computed by means of molecular dynamics simulations in explicit solvent and by Brownian dynamics simulations along the reaction path. Brownian dynamics makes use of the potential of mean force and hydrodynamic profile that we determined. Surprisingly, we find reasonable agreement between full-scale molecular dynamics and Brownian dynamics, despite the role of slow solvent relaxation in the assembly process. We found that molecular-scale hydrodynamic interactions are essential in describing the kinetics of assembly.

Efficient multiple time scale molecular dynamics: Using colored noise thermostats to stabilize resonances

Multiple time scale molecular dynamics enhances computational efficiency by updating slow motions less frequently than fast motions. However, in practice, the largest outer time step possible is limited not by the physical forces but by resonances between the fast and slow modes. In this paper we show that this problem can be alleviated by using a simple colored noise thermostatting scheme which selectively targets the high frequency modes in the system. For two sample problems, flexible water and solvated alanine dipeptide, we demonstrate that this allows the use of large outer time steps while still obtaining accurate sampling and minimizing the perturbation of the dynamics. Furthermore, this approach is shown to be comparable to constraining fast motions, thus providing an alternative to molecular dynamics with constraints.

How hydrophobic drying forces impact the kinetics of molecular recognition

A model of protein–ligand binding kinetics, in which slow solvent dynamics results from hydrophobic drying transitions, is investigated. Molecular dynamics simulations show that solvent in the receptor pocket can fluctuate between wet and dry states with lifetimes in each state that are long enough for the extraction of a separable potential of mean force and wet-to-dry transitions. We present a diffusive surface hopping model that is represented by a 2D Markovian master equation. One dimension is the standard reaction coordinate, the ligand–pocket separation, and the other is the solvent state in the region between ligand and binding pocket which specifies whether it is wet or dry. In our model, the ligand diffuses on a dynamic free-energy surface which undergoes kinetic transitions between the wet and dry states. The model yields good agreement with results from explicit solvent molecular dynamics simulation and an improved description of the kinetics of hydrophobic assembly. Furthermore, it is consistent with a “non-Markovian Brownian theory” for the ligand–pocket separation coordinate alone.

Momentum distribution, vibrational dynamics, and the potential of mean force in ice

By analyzing the momentum distribution obtained from path integral and phonon calculations we find that the protons in hexagonal ice experience an anisotropic quasiharmonic effective potential with three distinct principal frequencies that reflect molecular orientation. Due to the importance of anisotropy, anharmonic features of the environment cannot be extracted from existing experimental distributions that involve the spherical average. The full directional distribution is required, and we give a theoretical prediction for this quantity that could be verified in future experiments. Within the quasiharmonic context, anharmonicity in the ground-state dynamics of the proton is substantial and has quantal origin, a finding that impacts the interpretation of several spectroscopies.