Zero-point Energy Leakage Research Articles

A Strategy for Modeling Nonstatistical Reactivity Effects: Combining Chemical Activation Estimates with a Vibrational Relaxation Model.

The kinetics of many chemical reactions can be readily explained with a statistical approach, for example, using a form of transition state theory and comparing calculated Gibbs energies along the reaction coordinate(s). However, there are cases where this approach fails, notably when the vibrational relaxation of the molecule to its statistical equilibrium occurs on the same time scale as the reaction dynamics, whether it is caused by slow relaxation, a fast reaction, or both. These nonstatistical phenomena are then often explored computationally using (quasi)classical ab initio molecular dynamics by calculating a large number of trajectories while being prone to issues such as zero-point energy leakage. On the other side of the field, we see resource-intensive quantum dynamics simulations, which significantly limit the size of explorable systems. We find that using a Fermi's golden rule type of model for vibrational relaxation, based on anharmonic coupling constants, we can extract the same qualitative information while giving insights into how to enhance (or destroy) the bottlenecks causing the phenomena. We present this model as a middle ground for exploring complex nonstatistical behavior, capable of treating medium-sized organic molecules or biologically relevant fragments. We also cover the challenges involved, in particular quantifying the excess energy in terms of vibrational modes. Relying on readily available electronic structure methods and providing results in a simple master equation form, this model shows promise as a screening tool for opportunities in mode-selective chemistry without external control.

A Hessian-Free Method to Prevent Zero-Point Energy Leakage in Classical Trajectories.

The problem associated with the zero-point energy (ZPE) leak in classical trajectory calculations is well known. Since ZPE is a manifestation of the quantum uncertainty principle, there are no restrictions on energy during the classical propagation of nuclei. This phenomenon can lead to unphysical results, such as forming products without the ZPE in the internal vibrational degrees of freedom (DOFs). The ZPE leakage also permits reactions below the quantum threshold for the reaction. We have developed a new Hessian-free method, inspired by the Lowe-Andersen thermostat model, to prevent energy dipping below a threshold in the local-pair (LP) vibrational DOFs. The idea is to pump the leaked energy to the corresponding local vibrational mode taken from the other vibrational DOFs. We have applied the new correction protocol on the ab-initio ground-state molecular dynamics simulation of the water dimer (H2O)2, which dissociates due to unphysical ZPE spilling from high-frequency OH modes. The LP-ZPE method has been able to prevent the ZPE spilling of the OH stretching modes by pumping back the leaked energy into the corresponding modes, while this energy is taken from the other modes of the dimer itself, keeping the system as a microcanonical ensemble.

Quasiclassical Trajectory Study of the Si + SH Reaction on an Accurate Double Many-Body Expansion Potential Energy Surface

An accurate potential energy surface (PES) for the HSiS system based on MRCI+Q calculations extrapolated to the complete basis set limit is presented. Modeled with the double many-body expansion (DMBE) method, the PES provides an accurate description of the long-range interactions, including electrostatic and dispersion terms decaying as R-4, R-5, R-6, R-8, R-10 that are predicted from dipole moments, quadrupole moments, and dipolar polarizabilities, which are also calculated at the MRCI+Q level. The novel PES is then used in quasiclassical trajectory calculations to predict the rate coefficients of the Si + SH → SiS + H reaction, which has been shown to be a major source of the SiS in certain regions of the interstellar medium. An account of the zero-point energy leakage based on various nonactive models is also given. It is shown that the reaction is dominated by long-range forces, with the mechanism Si + SH → SiSH → SSiH → SiS + H being the most important one for all temperatures studied. Although SSiH corresponds to the global minimum of the PES, the contribution from the direct reaction Si + SH → SSiH → SiS + H is less than 0.5% for temperatures higher than 500 K. The rovibrational distributions of the products are calculated using the momentum Gaussian binning method and show that as the temperature is increased the average vibrational quantum number decreases while the rotational distribution spreads up to larger values.

Simulation of Nuclear Quantum Effects in Condensed Matter Systems via Quantum Baths

This paper reviews methods that aim at simulating nuclear quantum effects (NQEs) using generalized thermal baths. Generalized (or quantum) baths simulate statistical quantum features, and in particular zero-point energy effects, through non-Markovian stochastic dynamics. They make use of generalized Langevin Equations (GLEs), in which the quantum Bose–Einstein energy distribution is enforced by tuning the random and friction forces, while the system degrees of freedom remain classical. Although these baths have been formally justified only for harmonic oscillators, they perform well for several systems, while keeping the cost of the simulations comparable to the classical ones. We review the formal properties and main characteristics of classical and quantum GLEs, in relation with the fluctuation–dissipation theorems. Then, we describe the quantum thermostat and quantum thermal bath, the two generalized baths currently most used, providing several examples of applications for condensed matter systems, including the calculation of vibrational spectra. The most important drawback of these methods, zero-point energy leakage, is discussed in detail with the help of model systems, and a recently proposed scheme to monitor and mitigate or eliminate it—the adaptive quantum thermal bath—is summarised. This approach considerably extends the domain of application of generalized baths, leading, for instance, to the successful simulation of liquid water, where a subtle interplay of NQEs is at play. The paper concludes by overviewing further development opportunities and open challenges of generalized baths.

Ring Polymer Molecular Dynamics in Gas-Surface Reactions: Inclusion of Quantum Effects Made Simple.

Accurately modeling gas-surface collision dynamics presents a great challenge for theory, especially in the low-energy (or temperature) regime where quantum effects are important. Here, a path integral-based nonequilibrium ring polymer molecular dynamics (NE-RPMD) approach is adapted to calculate dissociative initial sticking probabilities (S0) of H2 on Cu(111) and D2O on Ni(111), revealing the distinct quantum nature in the two benchmark surface reactions. NE-RPMD successfully captures quantum tunneling in H2 dissociation at very low energies, where the quasi-classical trajectory (QCT) method suddenly fails. Additionally, QCT substantially overestimates S0 of D2O because of severe zero point energy (ZPE) leakage, even at collision energies greater than the ZPE-corrected barrier. Instead, NE-RPMD predicts S0 values of D2O in much improved agreement with reference results obtained by the quantum wavepacket method with reasonable corrections of the thermal contribution. Our results suggest NE-RPMD as a promising approach to model quantum effects in gas-surface reactions.

Quasi-Classical Trajectory Study of NH(3∑-) + NH(3∑-) Reactive Collisions.

A full-dimension quasi-classical trajectory study of collisions between two NH radicals is presented. Interatomic interactions are represented by a previously reported global six-dimensional potential energy surface for singlet electronic state of the N2H2 system. This study suggests that the formation of N2 from the collision of two NH radicals may occur via a one-step (NH + NH → N2 + H + H) or two-step (NH + NH → N2H + H → N2 + H + H) microscopic reaction mechanism. A fast vibrational energy redistribution is observed in the four-body complex in the latter mechanism. Excitation functions are presented and discussed. A variant of the vibrational energy quantum mechanical threshold method was used to correct the zero-point energy leakage in the classical calculations. The influence of reactant's rotational and vibrational energy on reactivity was investigated by state-specific calculations with one or both reactants excited. Reaction rate constants for the ground and some rotationally excited states are presented using an Arrhenius-Kooij-like functional form.

State dependent ring polymer molecular dynamics for investigating excited nonadiabatic dynamics

A recently proposed nonadiabatic ring polymer molecular dynamics (NRPMD) approach has shown to provide accurate quantum dynamics by incorporating explicit state descriptions and nuclear quantizations. Here, we present a rigorous derivation of the NRPMD Hamiltonian and investigate its performance on simulating excited state nonadiabatic dynamics. Our derivation is based on the Meyer-Miller-Stock-Thoss mapping representation for electronic states and the ring-polymer path-integral description for nuclei, resulting in the same Hamiltonian proposed in the original NRPMD approach. In addition, we investigate the accuracy of using NRPMD to simulate the photoinduced nonadiabatic dynamics in simple model systems. These model calculations suggest that NRPMD can alleviate the zero-point energy leakage problem that is commonly encountered in the classical Wigner dynamics and provide accurate excited state nonadiabatic dynamics. This work provides a solid theoretical foundation of the promising NRPMD Hamiltonian and demonstrates the possibility of using the state-dependent RPMD approach to accurately simulate electronic nonadiabatic dynamics while explicitly quantizing nuclei.

The Fluctuation-Dissipation Theorem as a Diagnosis and Cure for Zero-Point Energy Leakage in Quantum Thermal Bath Simulations.

Quantum thermal bath (QTB) simulations reproduce statistical nuclear quantum effects via a Langevin equation with a colored random force. Although this approach has proven efficient for a variety of chemical and condensed-matter problems, the QTB, as many other semiclassical methods, suffers from zero-point energy leakage (ZPEL). The absence of a reliable criterion to quantify the ZPEL without resorting to demanding comparisons with path integral-based calculations has so far hindered the use of the QTB for the simulation of real systems. In this work, we establish a quantitative connection between ZPEL in the QTB framework and deviations from the quantum fluctuation-dissipation theorem (FDT) that can be monitored along the simulation. This provides a rigorous general criterion to detect and quantify the ZPEL without any a priori knowledge of the system under study. We then use this criterion to build an adaptive QTB method that strictly enforces the quantum FDT at all frequencies via an on-the-fly, spectrally resolved fine-tuning of the system-bath coupling coefficients. The validity of the adaptive approach is first demonstrated on a simple two-oscillator model. It is then applied to two more realistic problems: the description of the vibrational properties of a model aluminum crystal at low temperature and the simulation of the liquid-solid phase transition in a 13-atom neon cluster. In both systems, the standard QTB results are strongly altered by the ZPEL, which can be essentially eliminated using the adaptive approach.

Nuclear quantum effects in molecular dynamics simulations

To take into account nuclear quantum effects on the dynamics of atoms, the path integral molecular dynamics (PIMD) method used since 1980s is based on the formalism developed by R. P. Feynman. However, the huge computation time required for the PIMD reduces its range of applicability. Another drawback is the requirement of additional techniques to access time correlation functions (ring polymer MD or centroid MD). We developed an alternative technique based on a quantum thermal bath (QTB) which reduces the computation time by a factor of ∼20. The QTB approach consists in a classical Langevin dynamics in which the white noise random force is replaced by a Gaussian random force having the power spectral density given by the quantum fluctuation-dissipation theorem. The method has yielded satisfactory results for weakly anharmonic systems: the quantum harmonic oscillator, the heat capacity of a MgO crystal, and isotope effects in 7LiH and 7LiD. Unfortunately, the QTB is subject to the problem of zero-point energy leakage (ZPEL) in highly anharmonic systems, which is inherent in the use of classical mechanics. Indeed, a part of the energy of the high-frequency modes is transferred to the low-frequency modes leading to a wrong energy distribution. We have shown that in order to reduce or even eliminate ZPEL, it is sufficient to increase the value of the frictional coefficient. Another way to solve the ZPEL problem is to combine the QTB and PIMD techniques. It requires the modification of the power spectral density of the random force within the QTB. This combination can also be seen as a way to speed up the PIMD.

A theoretical study of energy transfer in Ar(1S) + SO2( ') collisions: Cross sections and rate coefficients for vibrational transitions.

Vibrational transitions, induced by collisions between rare-gas atoms and molecules, play a key role in many problems of interest in physics and chemistry. A theoretical investigation of the translation-to-vibration (T-V) energy transfer process in argon atom and sulfur dioxide molecule collisions is presented here. For such a purpose, the framework of the quasi-classical trajectory (QCT) methodology was followed over the range of translational energies 2 ≤ Etr/kcal mol-1 ≤ 100. A new realistic potential energy surface (PES) for the ArSO2 system was developed using pairwise addition for the four-body energy term within the double many-body expansion. The topological features of the obtained function are compared with a previous one reported by Hippler et al. [J. Phys. Chem. 90, 6158 (1986)]. To test the accuracy of the PES, additional coupled cluster singles and doubles method with a perturbative contribution of connected triples calculations were carried out for the global minimum configuration. From dynamical calculations, the cross sections for the T-V excitation process indicate a barrier-type mechanism due to strong repulsive interactions between SO2 molecules and the Ar atom. Corrections to zero-point energy leakage in QCT were carried out using vibrational energy quantum mechanical threshold of the complex and variations. Rate coefficients and cross sections are calculated for some vibrational transitions using pseudo-quantization approaches of the vibrational energy of products. Main attributes of the title molecular collision are discussed and compared with available information in the literature.

Mode-specific and bond-selective dissociative chemisorption of CHD3 and CH2D2 on Ni(111) revisited using a new potential energy surface

Dissociative chemisorption of methane on a nickel surface is a prototypical system for studying mode-specific chemistry in gas-surface reactions. We recently developed a fifteen-dimensional potential energy surface for this system which has proven to be chemically accurate in reproducing the measured absolute dissociative sticking probabilities of CHD3 in thermal conditions and with vibrational excitation on Ni(111) at high incident energies. Here, using this new potential energy surface, we explored mode specificity and bond selectivity for CHD3 and CH2D2 dissociative chemisorption at low incidence energies down to ~50 kJ/mol via a quasi-classical trajectory method. Our calculated dissociation probabilities are consistent with previous theoretical and experimental ones with an average shift in translational energy of ~8 kJ/mol. Our results very well reproduce the C–H/C–D branching ratio upon the C–H local mode excitation, which can be rationalized by the sudden vector projection model. Quantitatively, however, the calculated dissociative sticking probabilities are systematically larger than experimental ones, due presumably to the artificial zero point energy leakage into reaction coordinate. Further high-dimensional quantum dynamics calculations are necessary for acquiring a chemically accurate description of methane dissociative chemisorption at low incident energies.

Herman-Kluk propagator is free from zero-point energy leakage

Quasiclassical techniques constitute a promising route to approximate quantum dynamics based on classical trajectories starting from a quantum-mechanically correct distribution. One of their main drawbacks is the so-called zero-point energy (ZPE) leakage, that is artificial redistribution of energy from the modes with high frequency and thus high ZPE to those with low frequency and ZPE due to classical equipartition. Here, we show that the elaborate semiclassical formalism based on the Herman-Kluk propagator is free from the ZPE leakage despite utilizing purely classical propagation. We demonstrate this with example applications for two- and three-dimensional anharmonically coupled oscillators. This finding opens the road to correct dynamical simulations of systems with a multitude of degrees of freedom that cannot be treated fully quantum-mechanically due to the exponential increase of the numerical effort.

A Ring Polymer Molecular Dynamics Approach to Study the Transition between Statistical and Direct Mechanisms in the H2 + H3+ → H3+ + H2 Reaction.

Because of its fundamental importance in astrochemistry, the H2 + H3+ → H3+ + H2 reaction has been studied experimentally in a wide temperature range. Theoretical studies of the title reaction significantly lag primarily because of the challenges associated with the proper treatment of the zero-point energy (ZPE). As a result, all previous theoretical estimates for the ratio between a direct proton-hop and indirect exchange (via the H5+ complex) channels deviate from the experiment, in particular, at lower temperatures where the quantum effects dominate. In this work, the ring polymer molecular dynamics (RPMD) method is applied to study this reaction, providing very good agreement with the experiment. RPMD is immune to the shortcomings associated with the ZPE leakage and is able to describe the transition from direct to indirect mechanisms below room temperature. We argue that RPMD represents a useful tool for further studies of numerous ZPE-sensitive chemical reactions that are of high interest in astrochemistry.

Zero-Point Energy Leakage in Quantum Thermal Bath Molecular Dynamics Simulations

The quantum thermal bath (QTB) has been presented as an alternative to path-integral-based methods to introduce nuclear quantum effects in molecular dynamics simulations. The method has proved to be efficient, yielding accurate results for various systems. However, the QTB method is prone to zero-point energy leakage (ZPEL) in highly anharmonic systems. This is a well-known problem in methods based on classical trajectories where part of the energy of the high-frequency modes is transferred to the low-frequency modes leading to a wrong energy distribution. In some cases, the ZPEL can have dramatic consequences on the properties of the system. Thus, we investigate the ZPEL by testing the QTB method on selected systems with increasing complexity in order to study the conditions and the parameters that influence the leakage. We also analyze the consequences of the ZPEL on the structural and vibrational properties of the system. We find that the leakage is particularly dependent on the damping coefficient and that increasing its value can reduce and, in some cases, completely remove the ZPEL. When using sufficiently high values for the damping coefficient, the expected energy distribution among the vibrational modes is ensured. In this case, the QTB method gives very encouraging results. In particular, the structural properties are well-reproduced. The dynamical properties should be regarded with caution although valuable information can still be extracted from the vibrational spectrum, even for large values of the damping term.

Applicability of Quantum Thermal Baths to Complex Many-Body Systems with Various Degrees of Anharmonicity.

The semiclassical method of quantum thermal baths by colored noise thermostats has been used to simulate various atomic systems in the molecular and bulk limits, at finite temperature and in moderately to strongly anharmonic regimes. In all cases, the method performs relatively well against alternative approaches in predicting correct energetic properties, including in the presence of phase changes, provided that vibrational delocalization is not too strong-neon appearing already as an upper limiting case. In contrast, the dynamical behavior inferred from global indicators such as the root-mean-square bond length fluctuation index or the vibrational spectrum reveals more marked differences caused by zero-point energy leakage, except in the case of isolated molecules with well separated vibrational modes. To correct for such deficiencies and reduce the undesired transfer among modes, empirical modifications of the noise power spectral density were attempted to better describe thermal equilibrium but still failed when used as semiclassical preparation for microcanonical trajectories.

Capture and dissociation in the complex-forming CH+H2→ CH2+H, CH+H2 reactions

The rate coefficients for the capture process CH + H(2)→ CH(3) and the reactions CH + H(2)→ CH(2) + H (abstraction), CH + H(2) (exchange) have been calculated in the 200-800 K temperature range, using the quasiclassical trajectory (QCT) method and the most recent global potential energy surface. The reactions, which are of interest in combustion and in astrochemistry, proceed via the formation of long-lived CH(3) collision complexes, and the three H atoms become equivalent. QCT rate coefficients for capture are in quite good agreement with experiments. However, an important zero point energy (ZPE) leakage problem occurs in the QCT calculations for the abstraction, exchange and inelastic exit channels. To account for this issue, a pragmatic but accurate approach has been applied, leading to a good agreement with experimental abstraction rate coefficients. Exchange rate coefficients have also been calculated using this approach. Finally, calculations employing QCT capture/phase space theory (PST) models have been carried out, leading to similar values for the abstraction rate coefficients as the QCT and previous quantum mechanical capture/PST methods. This suggests that QCT capture/PST models are a good alternative to the QCT method for this and similar systems.

Zero point energy leakage in condensed phase dynamics: An assessment of quantum simulation methods for liquid water

The approximate quantum mechanical ring polymer molecular dynamics (RPMD) and linearized semiclassical initial value representation (LSC-IVR) methods are compared and contrasted in a study of the dynamics of the flexible q-TIP4P/F water model at room temperature. For this water model, a RPMD simulation gives a diffusion coefficient that is only a few percent larger than the classical diffusion coefficient, whereas a LSC-IVR simulation gives a diffusion coefficient that is three times larger. We attribute this discrepancy to the unphysical leakage of initially quantized zero point energy (ZPE) from the intramolecular to the intermolecular modes of the liquid as the LSC-IVR simulation progresses. In spite of this problem, which is avoided by construction in RPMD, the LSC-IVR may still provide a useful approximation to certain short-time dynamical properties which are not so strongly affected by the ZPE leakage. We illustrate this with an application to the liquid water dipole absorption spectrum, for which the RPMD approximation breaks down at frequencies in the O-H stretching region owing to contamination from the internal modes of the ring polymer. The LSC-IVR does not suffer from this difficulty and it appears to provide quite a promising way to calculate condensed phase vibrational spectra.

Dynamics and kinetics of the S + HO2reaction: A theoretical study

We report a quasi-classical trajectory study of the S + HO2 reaction using a previously reported global potential energy surface for the ground electronic state of HSO2. Zero-point energy leakage is approximately accounted for by using the vibrational energy quantum mechanical threshold method. Calculations are carried out both for specific ro-vibrational states of the reactants and thermalized ones, with rate constants being reported as a function of temperature. The results suggest that the title reaction is capture type, with OH and SO showing as the most favorable products. The internal energy distribution of such products and the reaction mechanism are also investigated. © 2008 Wiley Periodicals, Inc. Int J Chem Kinet 40: 533–540, 2008

State-to-State Dynamics Analysis of the F + CHD3 Reaction: A Quasiclassical Trajectory Study

An exhaustive state-to-state dynamics study was performed to analyze the F + CHD3 --> FD(nu', j') + CHD2(nu) gas-phase abstraction reaction. Quasiclassical trajectory (QCT) calculations, including corrections to avoid zero-point energy leakage along the trajectories, were performed at different collision energies on an analytical potential energy surface (PES-2006) recently developed by our group. Whereas the CHD2 coproduct appears vibrationally and rotationally cold, most of the available energy appears as FD(nu') product vibrational energy, peaking at nu' = 2 and nu' = 3, with the population in the latter level growing as the energy increases. The excitation function rises from the threshold of the reaction and then levels off at higher energies, with the maximum contribution from the FD(nu' = 3) level. The state-specific FD(nu') scattering distributions correlated with the coproduct CHD2 in the nu4 = 2 and nu3 = 1 states, at different collision energies, show a steady change from backward to forward scattering as the energy increases. This similar behavior for the two coproduct vibrational states, nu4 = 2 and nu3 = 1, agrees qualitatively with the experimental measurements. Comparison with theoretical and experimental results for the isotopic analogues, F + CH4 and F + CD4, shows that the title reaction presents a direct mechanism, similar to the perdeuterated reaction, but contrasts with that of the F + CH4 reaction. These results for the dynamics of different isotopic variants, always in qualitative and sometimes in quantitative agreement with experiment, show the capacity of the PES-2006 surface to correctly describe the title reaction, even though there are differences that could be due to deficiencies of the PES but also to the known limitations of the classical treatment in the QCT method.

Role of the C−H Stretch Mode Excitation in the Dynamics of the Cl + CHD3 Reaction: A Quasi-classical Trajectory Calculation

To analyze the effect of the C-H stretch mode excitation on the dynamics of the Cl + CHD3 gas-phase abstraction reaction, an exhaustive state-to-state dynamics study was performed. This reaction can evolve along two channels: H-abstraction, CD3 + ClH, and D-abstraction, CHD2 + ClD. On an analytical potential energy surface constructed previously by our group, named PES-2005, quasi-classical trajectory calculations were performed at a collision energy of 0.18 eV, including corrections to avoid zero-point energy leakage along the trajectories. First, strong coupling between different vibrational modes in the entry valley was observed; i.e., the reaction is vibrationally nonadiabatic. Second, for the ground-state CHD3(nu=0) reaction, the diatomic fragments appeared in their ground states, and the H- and D-abstraction reactions showed similar reactivities. However, when the reactivity per atom is considered, the H is three times more reactive than the D atom. Third, when the C-H stretch mode is excited by one quantum, CHD3(nu1=1), the H-abstraction is strongly favored, and the C-H stretch excitation is maintained in the product CHD2(nu1=1) + ClD channel; i.e., the reaction shows mode selectivity, reproducing the experimental evidence, and also the reactivity of the vibrational ground state is increased, in agreement with experiment. Fourth, the state-to-state angular distributions of the CD3 and CHD2 products showed the products to be practically sideways for the reactant ground state, while the C-H excitation yielded a more forward scattering, reproducing the experimental data. The role of the zero-point energy correction was also analyzed, and we find that the dynamics results are very sensitive on how the ZPE issue is treated. Finally, a comparison is made with the similar H + CHD3(nu1=0,1) and Cl + CH4(nu1=0,1) reactions.