Thomas E. Markland

Ring-Polymer Molecular Dynamics: Quantum Effects in Chemical Dynamics from Classical Trajectories in an Extended Phase Space

This article reviews the ring-polymer molecular dynamics model for condensed-phase quantum dynamics. This model, which involves classical evolution in an extended ring-polymer phase space, provides a practical approach to approximating the effects of quantum fluctuations on the dynamics of condensed-phase systems. The review covers the theory, implementation, applications, and limitations of the approximation.

Competing quantum effects in the dynamics of a flexible water model

Numerous studies have identified large quantum mechanical effects in the dynamics of liquid water. In this paper, we suggest that these effects may have been overestimated due to the use of rigid water models and flexible models in which the intramolecular interactions were described using simple harmonic functions. To demonstrate this, we introduce a new simple point charge model for liquid water, q-TIP4P/F, in which the O-H stretches are described by Morse-type functions. We have parametrized this model to give the correct liquid structure, diffusion coefficient, and infrared absorption frequencies in quantum (path integral-based) simulations. The model also reproduces the experimental temperature variation of the liquid density and affords reasonable agreement with the experimental melting temperature of hexagonal ice at atmospheric pressure. By comparing classical and quantum simulations of the liquid, we find that quantum mechanical fluctuations increase the rates of translational diffusion and orientational relaxation in our model by a factor of around 1.15. This effect is much smaller than that observed in all previous simulations of empirical water models, which have found a quantum effect of at least 1.4 regardless of the quantum simulation method or the water model employed. The small quantum effect in our model is a result of two competing phenomena. Intermolecular zero point energy and tunneling effects destabilize the hydrogen-bonding network, leading to a less viscous liquid with a larger diffusion coefficient. However, this is offset by intramolecular zero point motion, which changes the average water monomer geometry resulting in a larger dipole moment, stronger intermolecular interactions, and a slower diffusion. We end by suggesting, on the basis of simulations of other potential energy models, that the small quantum effect we find in the diffusion coefficient is associated with the ability of our model to produce a single broad O-H stretching band in the infrared absorption spectrum.

Efficient stochastic thermostatting of path integral molecular dynamics

The path integral molecular dynamics (PIMD) method provides a convenient way to compute the quantum mechanical structural and thermodynamic properties of condensed phase systems at the expense of introducing an additional set of high frequency normal modes on top of the physical vibrations of the system. Efficiently sampling such a wide range of frequencies provides a considerable thermostatting challenge. Here we introduce a simple stochastic path integral Langevin equation (PILE) thermostat which exploits an analytic knowledge of the free path integral normal mode frequencies. We also apply a recently developed colored noise thermostat based on a generalized Langevin equation (GLE), which automatically achieves a similar, frequency-optimized sampling. The sampling efficiencies of these thermostats are compared with that of the more conventional Nosé-Hoover chain (NHC) thermostat for a number of physically relevant properties of the liquid water and hydrogen-in-palladium systems. In nearly every case, the new PILE thermostat is found to perform just as well as the NHC thermostat while allowing for a computationally more efficient implementation. The GLE thermostat also proves to be very robust delivering a near-optimum sampling efficiency in all of the cases considered. We suspect that these simple stochastic thermostats will therefore find useful application in many future PIMD simulations.

Nuclear Quantum Effects in Water and Aqueous Systems: Experiment, Theory, and Current Challenges

Nuclear quantum effects influence the structure and dynamics of hydrogen-bonded systems, such as water, which impacts their observed properties with widely varying magnitudes. This review highlights the recent significant developments in the experiment, theory, and simulation of nuclear quantum effects in water. Novel experimental techniques, such as deep inelastic neutron scattering, now provide a detailed view of the role of nuclear quantum effects in waters properties. These have been combined with theoretical developments such as the introduction of the principle of competing quantum effects that allows the subtle interplay of waters quantum effects and their manifestation in experimental observables to be explained. We discuss how this principle has recently been used to explain the apparent dichotomy in waters isotope effects, which can range from very large to almost nonexistent depending on the property and conditions. We then review the latest major developments in simulation algorithms and theory that have enabled the efficient inclusion of nuclear quantum effects in molecular simulations, permitting their combination with on-the-fly evaluation of the potential energy surface using electronic structure theory. Finally, we identify current challenges and future opportunities in this area of research.

An efficient ring polymer contraction scheme for imaginary time path integral simulations

A quantum simulation of an imaginary time path integral typically requires around n times more computational effort than the corresponding classical simulation, where n is the number of ring polymer beads (or imaginary time slices) used in the calculation. However, this estimate neglects the fact that the potential energies of many systems can be decomposed into a sum of rapidly varying short-range and slowly varying long-range contributions. For such systems, the computational effort of the path integral simulation can be reduced considerably by evaluating the long-range forces on a contracted ring polymer with fewer beads than are needed to evaluate the short-range forces. This idea is developed and then illustrated with an application to a flexible model of liquid water in which the intramolecular forces are evaluated with 32 beads, the oxygen-oxygen Lennard-Jones forces with seven, and the intermolecular electrostatic forces with just five. The resulting static and dynamic properties are within a few percent of those of a full 32-bead calculation, and yet they are obtained with a computational effort less than six times (rather than 32 times) that of a classical simulation. We hope that this development will encourage future studies of quantum mechanical fluctuations in liquid water and aqueous solutions and in many other systems with similar interaction potentials.

Unraveling quantum mechanical effects in water using isotopic fractionation

When two phases of water are at equilibrium, the ratio of hydrogen isotopes in each is slightly altered because of their different phase affinities. This isotopic fractionation process can be utilized to analyze water’s movement in the world’s climate. Here we show that equilibrium fractionation ratios, an entirely quantum mechanical property, also provide a sensitive probe to assess the magnitude of nuclear quantum fluctuations in water. By comparing the predictions of a series of water models, we show that those describing the OH chemical bond as rigid or harmonic greatly overpredict the magnitude of isotope fractionation. Models that account for anharmonicity in this coordinate are shown to provide much more accurate results because of their ability to give partial cancellation between inter- and intramolecular quantum effects. These results give evidence of the existence of competing quantum effects in water and allow us to identify how this cancellation varies across a wide-range of temperatures. In addition, this work demonstrates that simulation can provide accurate predictions and insights into hydrogen fractionation.

Quantum diffusion of hydrogen and muonium atoms in liquid water and hexagonal ice.

We have used the ring polymer molecular dynamics method to study the diffusion of muonium, hydrogen, and deuterium atoms in liquid water and hexagonal ice over a wide temperature range (8-361 K). Quantum effects are found to dramatically reduce the diffusion of muonium in water relative to that predicted by classical simulation. This leads to a simple explanation for the lack of any significant isotope effect in the observed diffusion coefficients of these species in the room temperature liquid. Our results indicate that the mechanism of the diffusion in liquid water is similar to the intercavity hopping mechanism observed in ice, supplemented by the diffusion of the cavities in the liquid. Within the same model, we have also been able to simulate the observed crossover in the c-axis diffusion coefficients of hydrogen and deuterium in hexagonal ice. Finally, we have been able to obtain good agreement with experimental data on the diffusion of muonium in hexagonal ice at 8 K, where the process is entirely quantum mechanical.

Efficient methods and practical guidelines for simulating isotope effects.

The shift in chemical equilibria due to isotope substitution is frequently exploited to obtain insight into a wide variety of chemical and physical processes. It is a purely quantum mechanical effect, which can be computed exactly using simulations based on the path integral formalism. Here we discuss how these techniques can be made dramatically more efficient, and how they ultimately outperform quasi-harmonic approximations to treat quantum liquids not only in terms of accuracy, but also in terms of computational cost. To achieve this goal we introduce path integral quantum mechanics estimators based on free energy perturbation, which enable the evaluation of isotope effects using only a single path integral molecular dynamics trajectory of the naturally abundant isotope. We use as an example the calculation of the free energy change associated with H/D and (16)O/(18)O substitutions in liquid water, and of the fractionation of those isotopes between the liquid and the vapor phase. In doing so, we demonstrate and discuss quantitatively the relative benefits of each approach, thereby providing a set of guidelines that should facilitate the choice of the most appropriate method in different, commonly encountered scenarios. The efficiency of the estimators we introduce and the analysis that we perform should in particular facilitate accurate ab initio calculation of isotope effects in condensed phase systems.

Quantum delocalization of protons in the hydrogen-bond network of an enzyme active site

Significance Because of the low mass of the proton, nuclear quantum effects can dramatically alter the properties of hydrogen-bond networks, especially when short and strong hydrogen bonds occur. Here, we combine experiments and state-of-the-art simulations that include the quantum nature of both the electrons and nuclei to show that the enzyme ketosteroid isomerase contains a hydrogen-bond network in its active site that facilitates extensive quantum proton delocalization. This leads to a 10,000-fold increase in the acidity of an active-site residue compared with the limit where the nuclei are classical particles. This work opens up new avenues for understanding the interplay between quantum effects and hydrogen bonding in biological systems containing strong hydrogen bonds. Enzymes use protein architectures to create highly specialized structural motifs that can greatly enhance the rates of complex chemical transformations. Here, we use experiments, combined with ab initio simulations that exactly include nuclear quantum effects, to show that a triad of strongly hydrogen-bonded tyrosine residues within the active site of the enzyme ketosteroid isomerase (KSI) facilitates quantum proton delocalization. This delocalization dramatically stabilizes the deprotonation of an active-site tyrosine residue, resulting in a very large isotope effect on its acidity. When an intermediate analog is docked, it is incorporated into the hydrogen-bond network, giving rise to extended quantum proton delocalization in the active site. These results shed light on the role of nuclear quantum effects in the hydrogen-bond network that stabilizes the reactive intermediate of KSI, and the behavior of protons in biological systems containing strong hydrogen bonds.

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.