Hybrid quantum/classical molecular dynamics for a proton transfer reaction coupled to a dissipative bath

A hybrid quantum/classical molecular dynamics approach is applied to a proton transfer reaction represented by a symmetric double well system coupled to a dissipative bath. In this approach, the proton is treated quantum mechanically and all bath modes are treated classically. The transition state theory rate constant is obtained from the potential of mean force, which is generated along a collective reaction coordinate with umbrella sampling techniques. The transmission coefficient, which accounts for dynamical recrossings of the dividing surface, is calculated with a reactive flux approach combined with the molecular dynamics with quantum transitions surface hopping method. The hybrid quantum/classical results agree well with numerically exact results in the spatial-diffusion-controlled regime, which is most relevant for proton transfer in proteins. This hybrid quantum/classical approach has already been shown to be computationally practical for studying proton transfer in large biological systems. These results have important implications for future applications to hydrogen transfer reactions in solution and proteins.     

Hybrid integral equation/simulation model for enhancing free energy computations

Integral equation theory is used for extrapolating free energy data from molecular simulations of a reference state with respect to a modification of the interaction potential. The methodology is applied to the correction of artefacts arising from potential shifting and truncation. Corrective contributions for the hydration free energy with respect to the full potential are analysed for the case that both the solute-solvent as well as the solvent-solvent potentials are truncated and modified by a shifted-force term, reaching beyond the range of the dielectric continuum approximation and simple long-range correction expressions. The model systems argon in water and pure water are used as examples for apolar and polar solutes, revealing significant correction contributions even for the short-ranged dispersive interactions and the magnitude of solute-solvent and solvent-solvent components. In comparison with simulation-based extrapolation techniques the integral equation method is shown to be capable of quantitatively predicting truncation artefacts at negligible computational overhead.     

An adiabatic linearized path integral approach for quantum time correlation functions: electronic transport in metal-molten salt solutions

We generalize the linearized path integral approach to evaluate quantum time correlation functions for systems best described by a set of nuclear and electronic degrees of freedom, restricting ourselves to the adiabatic approximation. If the operators in the correlation function are nondiagonal in the electronic states, then this adiabatic linearized path integral approximation for the thermal averaged quantum dynamics presents interesting and distinctive features, which we derive and explore in this paper. The capability of these approximations to accurately reproduce the behavior of physical systems is demonstrated by calculating the diffusion constant for an excess electron in a metal-molten salt solution.     

Hybrid classical quantum force field for modeling very large molecules

A coherent computational scheme on a very large molecule in which the subsystem that undergoes the most important electronic changes is treated by a semiempirical quantum chemical method, though the rest of the molecule is described by a classical force field, has been proposed recently. The continuity between the two subsystems is obtained by a strictly localized bond orbital, which is assumed to have transferable properties determined on model molecules. The computation of the forces acting on the atoms is now operating, giving rise to a hybrid classical quantum force field (CQFF ) which allows full energy minimization and modeling chemical changes in large biomolecules. As an illustrative example, we study the short hydrogen bonds and the proton-exchange process in the histidine-aspartic acid system of the catalytic triad of human neutrophil elastase. The CQFF approach reproduces the crystallographic data quite well, in opposition to a classical force field. The method also offers the possibility of switching off the electrostatic interaction between the quantum and the classical subsystems, allowing us to analyze the various components of the perturbation exerted by the macromolecule in the reactive part. Molecular dynamics confirm a fast proton exchange between the three possible energy wells. The method appears to be quite powerful and applicable to other cases of chemical interest such as surface reactivity of nonmetallic solids. © 1996 John Wiley & Sons, Inc.     

An adiabatic linearized path integral approach for quantum time-correlation functions II: a cumulant expansion method for improving convergence

Linearized mixed quantum-classical simulations are a promising approach for calculating time-correlation functions. At the moment, however, they suffer from some numerical problems that may compromise their efficiency and reliability in applications to realistic condensed-phase systems. In this paper, we present a method that improves upon the convergence properties of the standard algorithm for linearized calculations by implementing a cumulant expansion of the relevant averages. The effectiveness of the new approach is tested by applying it to the challenging computation of the diffusion of an excess electron in a metal-molten salt solution.     

Ab initio quantum mechanical/molecular mechanical simulation of electron transfer process: fractional electron approach

Electron transfer (ET) reactions are one of the most important processes in chemistry and biology. Because of the quantum nature of the processes and the complicated roles of the solvent, theoretical study of ET processes is challenging. To simulate ET processes at the electronic level, we have developed an efficient density functional theory (DFT) quantum mechanical (QM)/molecular mechanical (MM) approach that uses the fractional number of electrons as the order parameter to calculate the redox free energy of ET reactions in solution. We applied this method to study the ET reactions of the aqueous metal complexes Fe(H(2)O)(6)(2+/3+) and Ru(H(2)O)(6)(2+/3+). The calculated oxidation potentials, 5.82 eV for Fe(II/III) and 5.14 eV for Ru(II/III), agree well with the experimental data, 5.50 and 4.96 eV, for iron and ruthenium, respectively. Furthermore, we have constructed the diabatic free energy surfaces from histogram analysis based on the molecular dynamics trajectories. The resulting reorganization energy and the diabatic activation energy also show good agreement with experimental data. Our calculations show that using the fractional number of electrons (FNE) as the order parameter in the thermodynamic integration process leads to efficient sampling and validate the ab initio QM/MM approach in the calculation of redox free energies.     

OH hydrogen abstraction reactions from alanine and glycine: A quantum mechanical approach

Density functional theory (B3LYP and BHandHLYP) and unrestricted second‐order Møller–Plesset (MP2) calculations have been performed using 3‐21G, 6‐31G(d,p), and 6‐311 G(2d,2p) basis sets, to study the OH hydrogen abstraction reaction from alanine and glycine. The structures of the different stationary points are discussed. Ring‐like structures are found for all the transition states. Reaction profiles are modeled including the formation of prereactive complexes, and very low or negative net energy barriers are obtained depending on the method and on the reacting site. ZPE and thermal corrections to the energy for all the species, and BSSE corrections for B3LYP activation energies are included. A complex mechanism involving the formation of a prereactive complex is proposed, and the rate coefficients for the overall reactions are calculated using classical transition state theory. The predicted values of the rate coefficients are 3.54×10⁸ L?mol^?1?s^?1 for glycine and 1.38×10⁹ L?mol^?1?s^?1 for alanine. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1138–1153, 2001     

Direct assessment of quantum nuclear effects on hydrogen bond strength by constrained-centroid ab initio path integral molecular dynamics

The impact of quantum nuclear effects on hydrogen (H-) bond strength has been inferred in earlier work from bond lengths obtained from path integral molecular dynamics (PIMD) simulations. To obtain a direct quantitative assessment of such effects, we use constrained-centroid PIMD simulations to calculate the free energy changes upon breaking the H-bonds in dimers of HF and water. Comparing ab initio simulations performed using PIMD and classical nucleus molecular dynamics (MD), we find smaller dissociation free energies with the PIMD method. Specifically, at 50 K, the H-bond in (HF)(2) is about 30% weaker when quantum nuclear effects are included, while that in (H(2)O)(2) is about 15% weaker. In a complementary set of simulations, we compare unconstrained PIMD and classical nucleus MD simulations to assess the influence of quantum nuclei on the structures of these systems. We find increased heavy atom distances, indicating weakening of the H-bond consistent with that observed by direct calculation of the free energies of dissociation.     

Nonequilibrium transport in quantum impurity models: exact path integral simulations

We simulate the nonequilibrium dynamics of two generic many-body quantum impurity models by employing the recently developed iterative influence-functional path integral method [Phys. Rev. B: Condens. Matter, 2010, 82, 205323]. This general approach is presented here in the context of quantum transport in molecular electronic junctions. Models of particular interest include the single impurity Anderson model and the related spinless two-state Anderson dot. In both cases we study the time evolution of the dot occupation and the current characteristics at finite temperature. A comparison to mean-field results is presented, when applicable.     

Hybrid Monte Carlo implementation of the Fourier path integral algorithm

This paper formulates a hybrid Monte Carlo implementation of the Fourier path integral (FPI-HMC) approach with partial averaging. Such a hybrid Monte Carlo approach allows one to generate collective moves through configuration space using molecular dynamics while retaining the computational advantages associated with the Fourier path integral Monte Carlo method. In comparison with the earlier Metropolis Monte Carlo implementations of the FPI algorithm, the present HMC method is shown to be significantly more efficient for quantum Lennard-Jones solids and suggests that such algorithms may prove useful for efficient simulations of a range of atomic and molecular systems.     

A comparative study of imaginary time path integral based methods for quantum dynamics

The recently introduced approximate many-body quantum simulation method, ring polymer molecular dynamics (RPMD), is compared to the centroid molecular dynamics method (CMD). Comparisons of simulation results for liquid para-hydrogen at two state points and liquid ortho-deuterium at one state point are presented. The calculated quantum correlation functions for the two methods are shown to be in good agreement with one another for a large portion of the time spectrum. However, as the quantum mechanical nature of the system increases, RPMD is less accurate in predicting the kinetic energy of the system than is CMD. A simplified and highly efficient algorithm is proposed which largely corrects this deficiency.     

Stereospecificity for the hydrogen transfer of pyridoxal enzyme reactions

We have studied the stereospecificities of various pyridoxal 5'-phosphate dependent enzymes for the hydrogen transfer between the C-4' of a bound coenzyme and the C-2 of a substrate in the transamination catalyzed by the enzymes. Prior to our studies, pyridoxal enzymes so far studied were reported to catalyze the hydrogen transfer only on the si-face of the planar imine intermediate formed from substrate and coenzyme. This finding had been considered as the evidence that pyridoxal enzymes have evolved divergently from a common ancestral protein, because identity in the stereospecificity reflects the similarity in the active-site structure, in particular in the geometrical relationship between the coenzyme and the active site base participating in the hydrogen transfer. However, we found that D-amino acid aminotransferase, branched-chain L-amino acid aminotransferase, and 4-amino-4-deoxychorismate lyase catalyze the re-face specific hydrogen transfer, and that amino acid racemases catalyze the nonstereospecific hydrogen transfer. These findings suggest the convergent evolution of pyridoxal enzymes. Crystallographical studies have shown that the stereospecificity reflects the active-site structure of the enzymes, and that the enzymes with the same fold exhibit the same stereospecificity. The active site structure with the catalytic base being situated on the specific face of the cofactor has been conserved during the evolution among the pyridoxal enzymes of the same family.     

Computer simulation of quantum dynamics in a classical spin environment

In this paper, a formalism for studying the dynamics of quantum systems coupled to classical spin environments is reviewed. The theory is based on generalized antisymmetric brackets and naturally predicts open-path off-diagonal geometric phases in the evolution of the density matrix. It is shown that such geometric phases must also be considered in the quantum–classical Liouville equation for a classical bath with canonical phase space coordinates; this occurs whenever the adiabatics basis is complex (as in the case of a magnetic field coupled to the quantum subsystem). When the quantum subsystem is weakly coupled to the spin environment, non-adiabatic transitions can be neglected and one can construct an effective non-Markovian computer simulation scheme for open quantum system dynamics in classical spin environments. In order to tackle this case, integration algorithms based on the symmetric Trotter factorization of the classical-like spin propagator are derived. Such algorithms are applied to a model comprising a quantum two-level system coupled to a single classical spin in an external magnetic field. Starting from an excited state, the population difference and the coherences of this two-state model are simulated in time while the dynamics of the classical spin is monitored in detail. It is the author’s opinion that the numerical evidence provided in this paper is a first step toward developing the simulation of quantum dynamics in classical spin environments into an effective tool. In turn, the ability to simulate such a dynamics can have a positive impact on various fields, among which, for example, nanoscience.     

Computer simulation of quantum melting in hydrogen clusters.

We introduce a new criterion, based on multipole dynamical correlations calculated within reptation quantum Monte Carlo, to discriminate between a melting versus freezing behavior in quantum clusters. This criterion is applied to small clusters of para-hydrogen molecules (both pristine and doped with a CO chromophore), for cluster sizes of around twelve molecules. This is a magic size at which para-hydrogen clusters display an icosahedral structure and a large stability. Despite the similar geometric structure of CO@(pH2)12 and (pH2)13, the first system has a rigid, crystalline, behavior; the second behaves more like a superfluid (or, possibly, a supersolid).     

A quantum generalization of intrinsic reaction coordinate using path integral centroid coordinates

We propose a generalization of the intrinsic reaction coordinate (IRC) for quantum many-body systems described in terms of the mass-weighted ring polymer centroids in the imaginary-time path integral theory. This novel kind of reaction coordinate, which may be called the "centroid IRC," corresponds to the minimum free energy path connecting reactant and product states with a least amount of reversible work applied to the center of masses of the quantum nuclei, i.e., the centroids. We provide a numerical procedure to obtain the centroid IRC based on first principles by combining ab initio path integral simulation with the string method. This approach is applied to NH(3) molecule and N(2)H(5) (-) ion as well as their deuterated isotopomers to study the importance of nuclear quantum effects in the intramolecular and intermolecular proton transfer reactions. We find that, in the intramolecular proton transfer (inversion) of NH(3), the free energy barrier for the centroid variables decreases with an amount of about 20% compared to the classical one at the room temperature. In the intermolecular proton transfer of N(2)H(5) (-), the centroid IRC is largely deviated from the "classical" IRC, and the free energy barrier is reduced by the quantum effects even more drastically.     

A classical analysis of quantum resonances in isotopic collinear H + H2 reactions

Isotopic substitution of hydrogen by muonium in the collinear H + H₂ reaction causes dramatic changes in the resonance patterns of the quantum reaction probabilities. These changes are explained by a classical model. Using the systems' resonant orbits, we show that resonances appear close to energetic thresholds of new vibrational channels.