Quantum-classical Liouville dynamics of nonadiabatic proton transfer

A proton transfer reaction in a linear hydrogen-bonded complex dissolved in a polar solvent is studied using mixed quantum-classical Liouville dynamics. In this system, the proton is treated quantum mechanically and the remainder of the degrees of freedom is treated classically. The rates and mechanisms of the reaction are investigated using both adiabatic and nonadiabatic molecular dynamics. We use a nonadiabatic dynamics algorithm which allows the system to evolve on single adiabatic surfaces and on coherently coupled pairs of adiabatic surfaces. Reactive-flux correlation function expressions are used to compute the rate coefficients and the role of the dynamics on the coherently coupled surfaces is elucidated.     

Nonadiabatic quantum-classical reaction rates with quantum equilibrium structure

Time correlation function expressions for quantum reaction-rate coefficients are computed in a quantum-classical limit. This form for the correlation function retains the full quantum equilibrium structure of the system in the spectral density function but approximates the time evolution of the operator by quantum-classical Liouville dynamics. Approximate analytical expressions for the spectral density function, which incorporate quantum effects in the many-body environment and reaction coordinate, are derived. The results of numerical simulations of the reaction rate are presented for a reaction model in which a two-level system is coupled to a bistable oscillator which is, in turn, coupled to a bath of harmonic oscillators. The nonadiabatic quantum-classical dynamics is simulated in terms of an ensemble of surface-hopping trajectories and the effects of the quantum equilibrium structure on the reaction rate are discussed.     

Analysis of kinetic isotope effects for nonadiabatic reactions

Factors influencing the rates of quantum mechanical particle transfer reactions in many-body systems are discussed. The investigations are carried out on a simple model for a proton transfer reaction that captures generic features seen in more realistic models of condensed phase systems. The model involves a bistable quantum oscillator coupled to a one-dimensional double-well reaction coordinate, which is in turn coupled to a bath of harmonic oscillators. Reactive-flux correlation functions that involve quantum-classical Liouville dynamics for chemical species operators and quantum equilibrium sampling are used to estimate the reaction rates. Approximate analytical expressions for the quantum equilibrium structure are derived. Reaction rates are shown to be influenced significantly by both the quantum equilibrium structure and nonadiabatic dynamics. Nonadiabatic dynamical effects are found to play the major role in determining the magnitude of the kinetic isotope effect for the model transfer reaction.     

Transport properties of quantum-classical systems

Correlation function expressions for calculating transport coefficients for quantum-classical systems are derived. The results are obtained by starting with quantum transport coefficient expressions and replacing the quantum time evolution with quantum-classical Liouville evolution, while retaining the full quantum equilibrium structure through the spectral density function. The method provides a variety of routes for simulating transport coefficients of mixed quantum-classical systems, composed of a quantum subsystem and a classical bath, by selecting different but equivalent time evolution schemes of any operator or the spectral density. The structure of the spectral density is examined for a single harmonic oscillator where exact analytical results can be obtained. The utility of the formulation is illustrated by considering the rate constant of an activated quantum transfer process that can be described by a many-body bath reaction coordinate.     

Deep tunneling dominates the biologically important hydride transfer reaction from NADH to FMN in morphinone reductase

The temperature dependence of the primary kinetic isotope effect (KIE), combined temperature-pressure studies of the primary KIE, and studies of the alpha-secondary KIE previously led us to infer that hydride transfer from nicotinamide adenine dinucleotide to flavin mononucleotide in morphinone reductase proceeds via environmentally coupled hydride tunneling. We present here a computational analysis of this hydride transfer reaction using QM/MM molecular dynamics simulations and variational transition-state theory calculations. Our calculated primary and secondary KIEs are in good agreement with the corresponding experimental values. Although the experimentally observed KIE lies below the semiclassical limit, our calculations suggest that approximately 99% of the reaction proceeds via tunneling: this is the first "deep tunneling" reaction observed for hydride transfer. We also show that the dominant tunneling mechanism is controlled by the isotope at the primary rather than the secondary position: with protium in the primary position, large-curvature tunneling dominates, whereas with deuterium in this position, small-curvature tunneling dominates. Also, our study is consistent with tunneling being preceded by reorganization: in the reactant, the rings of the nicotinamide and isoalloxazine moieties are stacked roughly parallel to each other, and as the system moves toward a "tunneling-ready" configuration, the nicotinamide ring rotates to become almost perpendicular to the isoalloxazine ring.     

Proton transfer in nanoconfined polar solvents. II. Adiabatic proton transfer dynamics

The reaction dynamics for a model phenol-amine proton transfer system in a confined methyl chloride solvent have been simulated by mixed quantum-classical molecular dynamics. In this approach, the proton vibration is treated quantum mechanically (and adiabatically), while the rest of the system is described classically. Nonequilibrium trajectories are used to determine the proton transfer reaction rate constant. The reaction complex and methyl chloride solvent are confined in a smooth, hydrophobic spherical cavity, and radii of 10, 12, and 15 A have been considered. The effects of the cavity radius and the heavy atom (hydrogen bond) distance on the reaction dynamics are considered, and the mechanism of the proton transfer is examined in detail by analysis of the trajectories.     

Theories and Applications of Mixed Quantum-Classical Non-adiabatic Dynamics

Electronically non-adiabatic processes are essential parts of photochemical process, collisions of excited species, electron transfer processes, and quantum information processing. Various non-adiabatic dynamics methods and their numerical implementation have been developed in the last decades. This review summarizes the most significant development of mixed quantum-classical methods and their applications which mainly include the Liouville equation, Ehrenfest mean-field, trajectory surface hopping, and multiple spawning methods. The recently developed quantum trajectory mean-field method that accounts for the decoherence corrections in a parameter-free fashion is discussed in more detail.     

Proton-transfer reactions between nitroalkanes and hydroxide ion under non-steady-state conditions. Apparent and real kinetic isotope effects.

The kinetics of the proton-transfer reactions between 1-nitro-1-(4-nitrophenyl)ethane (NNPE(H(D))) and hydroxide ion in water/acetonitrile (50/50 vol %) were studied at temperatures ranging from 289 to 319 K. The equilibrium constants for the reactions are large under these conditions, ensuring that the back reaction is not significant. The extent of reaction/time profiles during the first half-lives are compared with theoretical data for the simple single-step mechanism and a 2-step mechanism involving initial donor/acceptor complex formation followed by unimolecular proton transfer and dissociation of ions. In all cases, the profiles for the reactions of both NNPE(H) and NNPE(D) deviate significantly from those expected for the simple single-step mechanism. Excellent fits of experimental data with theoretical data for the complex mechanism, in the pre-steady-state time period, were observed in all cases. At all base concentrations (0.5 to 5.0 mM) and at all temperatures the apparent kinetic isotope effects (KIE(app)) were observed to increase with increasing extent of reaction. Resolution of the kinetics into microscopic rate constants at 298 K resulted in a real kinetic isotope effect (KIE(real)) for the proton-transfer step equal to 22. Significant proton tunneling was further indicated by the temperature dependence of the rate constants for proton and deuteron transfers: KIE(real) ranging from 17 to 26, E(a)(D) -- E(a)(H) equal 2.8 kcal/mol, and A(D)/A(H) equal to 4.95.     

Evidence for environmentally coupled hydrogen tunneling during dihydrofolate reductase catalysis

Hydride transfer during catalysis by dihydrofolate reductase from Thermotoga maritima has been studied by stopped flow spectroscopy. The reduction of dihydrofolate by NADPH showed a biphasic temperature dependence of the deuterium kinetic isotope effect. At temperatures above 25 degrees C the KIE was temperature independent, while the reaction rates were strongly temperature dependent. Below 25 degrees C the KIE becomes dependent on temperature, and the ratio of the preexponential factors is inverse, suggesting a greater role for active dynamics that modulate the tunneling distance.     

Pressure effects on enzyme-catalyzed quantum tunneling events arise from protein-specific structural and dynamic changes

The rate and kinetic isotope effect (KIE) on proton transfer during the aromatic amine dehydrogenase-catalyzed reaction with phenylethylamine shows complex pressure and temperature dependences. We are able to rationalize these effects within an environmentally coupled tunneling model based on constant pressure molecular dynamics (MD) simulations. As pressure appears to act anisotropically on the enzyme, perturbation of the reaction coordinate (donor-acceptor compression) is, in this case, marginal. Therefore, while we have previously demonstrated that pressure and temperature dependences can be used to infer H-tunneling and the involvement of promoting vibrations, these effects should not be used in the absence of atomistic insight, as they can vary greatly for different enzymes. We show that a pressure-dependent KIE is not a definitive hallmark of quantum mechanical H-tunneling during an enzyme-catalyzed reaction and that pressure-independent KIEs cannot be used to exclude tunneling contributions or a role for promoting vibrations in the enzyme-catalyzed reaction. We conclude that coupling of MD calculations with experimental rate and KIE studies is required to provide atomistic understanding of pressure effects in enzyme-catalyzed reactions.     

混合量子-经典方法计算电荷转移速率及其在实际体系中的应用(英文)

混合量子-经典方法在复杂分子体系动力学过程的模拟方面有重要应用.我们采用Ehrenfest方法、surfacehopping方法和混合量子经典Liouville方程计算了在非绝热极限下的电荷转移速率.然后将这三种方法应用于有机半导体材料电荷转移速率的计算.研究结果发现,Ehrenfest方法和surface hopping方法可能严重偏离正确的结果.偏离的原因是这两种方法没有正确处理相干项的运动,而且这种偏离在涉及到高频模式时显得更加严重.     

Fully adaptive propagation of the quantum-classical Liouville equation

In mixed quantum-classical molecular dynamics few but important degrees of freedom of a dynamical system are modeled quantum-mechanically while the remaining ones are treated within the classical approximation. Rothe methods established in the theory of partial differential equations are used to control both temporal and spatial discretization errors on grounds of a global tolerance criterion. The TRAIL (trapezoidal rule for adaptive integration of Liouville dynamics) scheme [I. Horenko and M. Weiser, J. Comput. Chem. 24, 1921 (2003)] has been extended to account for nonadiabatic effects in molecular dynamics described by the quantum-classical Liouville equation. In the context of particle methods, the quality of the spatial approximation of the phase-space distributions is maximized while the numerical condition of the least-squares problem for the parameters of particles is minimized. The resulting dynamical scheme is based on a simultaneous propagation of moving particles (Gaussian and Dirac deltalike trajectories) in phase space employing a fully adaptive strategy to upgrade Dirac to Gaussian particles and, vice versa, downgrading Gaussians to Dirac-type trajectories. This allows for the combination of Monte-Carlo-based strategies for the sampling of densities and coherences in multidimensional problems with deterministic treatment of nonadiabatic effects. Numerical examples demonstrate the application of the method to spin-boson systems in different dimensionality. Nonadiabatic effects occurring at conical intersections are treated in the diabatic representation. By decreasing the global tolerance, the numerical solution obtained from the TRAIL scheme are shown to converge towards exact results.     

Hydrogen isotope effects and mechanism of aqueous ozone and peroxone decompositions

Hydrogen peroxide exalts the reactivity of aqueous ozone by reasons that remain obscure. Should H2O2 enhance free radical production, as it is generally believed, a chain mechanism propagated by (.OH/.O2-) species would account for O3 decomposition rates in neat H2O, HR-O3, and in peroxone (O3 + H2O2) solutions, HPR-O3. We found, however, that: (1) the radical mechanism correctly predicts HR-O3 but vastly overestimates HPR-O3, (2) solvent deuteration experiments preclude radical products from the (O3 + HO2-) reaction. The modest kinetic isotope effect (KIE) we measure in H2O/D2O: HR-O3/DR-O3 = 1.5 +/- 0.3, is compatible with a chain process driven by electron- and/or O-atom transfer processes. But the large KIE found in peroxone: HPR-O3/DPR-O3 = 19.6 +/- 4.0, is due to an elementary (O3 + HO2-) reaction involving H-O2- bond cleavage. Since the KIE for the hypothetical H-atom transfer: O3 + HO2- HO3. +.O2-, would emerge as a KIE1/2 factor in the rates of the ensuing radical chain, the magnitude of the observed KIE must be associated with the hydride transfer reaction that yields a diamagnetic species: O3 + HO2- HO3- + O2. HO3-/H2O3 may be the bactericidal trioxide recently identified in the antibody-catalyzed addition of O2(1Deltag) to H2O.     

A derivation of the mixed quantum-classical Liouville equation from the influence functional formalism

We show that the mixed quantum-classical Liouville equation is equivalent to linearizing the forward-backward action in the influence functional. Derivations are provided in terms of either the diabatic or adiabatic basis sets. An application of the mixed quantum-classical Liouville equation for calculating the memory kernel of the generalized quantum master equation is also presented. The accuracy and computational feasibility of such an approach is demonstrated in the case of a two-level system nonlinearly coupled to an anharmonic bath.