Search

The observation of long-lived electronic coherence in photosynthetic excitation energy transfer (EET) by Engel et al. [Nature (London) 446, 782 (2007)] raises questions about the role of the protein environment in protecting this coherence and the significance of the quantum coherence in light harvesting efficiency. In this paper we explore the applicability of the Redfield equation in its full form, in the secular approximation and with neglect of the imaginary part of the relaxation terms for the study of these phenomena. We find that none of the methods can give a reliable picture of the role of the environment in photosynthetic EET. In particular the popular secular approximation (or the corresponding Lindblad equation) produces anomalous behavior in the incoherent transfer region leading to overestimation of the contribution of environment-assisted transfer. The full Redfield expression on the other hand produces environment-independent dynamics in the large reorganization energy region. A companion paper presents an improved approach, which corrects these deficiencies [A. Ishizaki and G. R. Fleming, J. Chem. Phys. 130, 234111 (2009)]. [ABSTRACT FROM AUTHOR]

In this paper, we examine how and when quantum evolution can be approximated in terms of (generalized) classical dynamics in calculations of correlation functions, with a focus on the symmetrized time correlation function introduced by Schofield. To that end, this function is expressed as a path integral in complex time and written in terms of sum and difference path variables. Taylor series expansion of the path integral's exponent to first and second order in the difference variables leads to two original developments. The first order expansion is used to obtain a simple, path integral based, derivation of the so-called Schofield's quantum correction factor. The second order result is employed to show how quantum mechanical delocalization manifests itself in the approximation of the correlation function and hinders, even in the semiclassical limit, the interpretation of the propagators in terms of sets of guiding classical trajectories dressed with appropriate weights. [ABSTRACT FROM AUTHOR]

A method for constructing diabatic potential-energy matrices from ab initio quantum chemistry data is described and tested for use in exact quantum reactive scattering. The method is a refinement of that presented in a previous paper, in that it accounts for the presence of the nonremovable derivative coupling. The accuracy of quantum dynamics on this type of diabatic potential is tested by comparison with an analytic model and for an ab initio description of the two lowest-energy states of H3. [ABSTRACT FROM AUTHOR]

Expressions for reaction rate constants in multidimensional chemical systems are derived by applying semiclassical approximation to the quantum path integrals of the ImF formulation of reaction rate theory. First, the transverse degrees of freedom orthogonal to the reaction coordinate are treated within the steepest descent approximation, after which the semiclassical approximation is applied to the remaining reaction coordinate. Thus derived, the semiclassical expressions account for the multidimensional nature of quantum effects and accurately incorporate nuclear quantum effects such as multidimensional tunneling and zero point energies. The obtained expressions are applicable in the broad temperature range from the deep tunneling to high-temperature regimes. The present paper provides derivation of the semiclassical instanton expressions proposed by Kryvohuz [J. Chem. Phys. 134, 114103 (2011)]. [ABSTRACT FROM AUTHOR]

Most of today's molecular-dynamics simulations of materials are based on the Born-Oppenheimer approximation. There are many cases, however, in which the coupling of the electrons and nuclei is important and it is necessary to go beyond the Born-Oppenheimer approximation. In these methods, the non-adiabatic coupling vectors are fundamental since they represent the link between the classical atomic motion of the nuclei and the time evolution of the quantum electronic state. In this paper we analyze the calculation of non-adiabatic coupling vectors in a basis set of local orbitals and derive an expression to calculate them in a practical and computationally efficient way. Some examples of the application of this expression using a local-orbital density functional theory approach are presented for a few simple molecules: H3, formaldimine, and azobenzene. These results show that the approach presented here, using the Slater transition-state density, is a very promising way for the practical calculation of non-adiabatic coupling vectors for large systems. [ABSTRACT FROM AUTHOR]

This paper shows how a compact finite difference Hessian approximation scheme can be proficiently implemented into semiclassical initial value representation molecular dynamics. Effects of the approximation on the monodromy matrix calculation are tested by propagating initial sampling distributions to determine power spectra for analytic potential energy surfaces and for 'on the fly' carbon dioxide direct dynamics. With the approximation scheme the computational cost is significantly reduced, making ab initio direct semiclassical dynamics computationally more feasible and, at the same time, properly reproducing important quantum effects inherent in the monodromy matrix and the pre-exponential factor of the semiclassical propagator. [ABSTRACT FROM AUTHOR]

Mixed quantum/classical (MQC) simulations treat the majority of a system classically and reserve quantum mechanics only for a few degrees of freedom that actively participate in the chemical process(es) of interest. In MQC calculations, the quantum and classical degrees of freedom are coupled together using pseudopotentials. Although most pseudopotentials are developed empirically, there are methods for deriving pseudopotentials using the results of quantum chemistry calculations, which guarantee that the explicitly-treated valence electron wave functions remain orthogonal to the implicitly-treated core electron orbitals. Whether empirical or analytically derived in nature, to date all such pseudopotentials have been subject to the frozen core approximation (FCA) that ignores how changes in the nuclear coordinates alter the core orbitals, which in turn affects the wave function of the valence electrons. In this paper, we present a way to go beyond the FCA by developing pseudopotentials that respond to these changes. In other words, we show how to derive an analytic expression for a pseudopotential that is an explicit function of nuclear coordinates, thus accounting for the polarization effects experienced by atomic cores in different chemical environments. We then use this formalism to develop a coordinate-dependent pseudopotential for the bonding electron of the sodium dimer cation molecule and we show how the analytic representation of this potential can be used in one-electron MQC simulations that provide the accuracy of a fully quantum mechanical Hartree-Fock (HF) calculation at all internuclear separations. We also show that one-electron MQC simulations of Na2+ using our coordinate-dependent pseudopotential provide a significant advantage in accuracy compared to frozen core potentials with no additional computational expense. This is because use of a frozen core potential produces a charge density for the bonding electron of Na2+ that is too localized on the molecule, leading to significant overbinding of the valence electron. This means that FCA calculations are subject to inaccuracies of order ∼10% in the calculated bond length and vibrational frequency of the molecule relative to a full HF calculation; these errors are fully corrected by using our coordinate-dependent pseudopotential. Overall, our findings indicate that even for molecules like Na2+, which have a simple electronic structure that might be expected to be well-treated within the FCA, the importance of including the effects of the changing core molecular orbitals on the bonding electrons cannot be overlooked. [ABSTRACT FROM AUTHOR]

In this paper we further explore and develop the quantum continuum mechanics (QCM) of Tao et al. [Phys. Rev. Lett. 103, 086401 (2009)] with the aim of making it simpler to use in practice. Our simplifications relate to the non-interacting part of the QCM equations, and primarily refer to practical implementations in which the groundstate stress tensor is approximated by its Kohn-Sham (KS) version. We use the simplified approach to directly prove the exactness of QCM for one-electron systems via an orthonormal formulation. This proof sheds light on certain physical considerations contained in the QCM theory and their implication on QCM-based approximations. The one-electron proof then motivates an approximation to the QCM (exact under certain conditions) expanded on the wavefunctions of the KS equations. Particular attention is paid to the relationships between transitions from occupied to unoccupied KS orbitals and their approximations under the QCM. We also demonstrate the simplified QCM semianalytically on an example system. [ABSTRACT FROM AUTHOR]

The dynamical properties of liquid water play an important role in many processes in nature. In this paper, we focus on the infrared (IR) absorption spectrum of liquid water based on the linearized semiclassical initial value representation (LSC-IVR) with the local Gaussian approximation (LGA) [J. Liu and W. H. Miller, J. Chem. Phys. 131, 074113 (2009)] and an ab initio based, flexible, polarizable Thole-type model (TTM3-F) [G. S. Fanourgakis and S. S. Xantheas, J. Chem. Phys. 128, 074506 (2008)]. Although the LSC-IVR (LGA) gives the exact result for the isolated three-dimensional shifted harmonic stretching model, it yields a blueshifted peak position for the more realistic anharmonic stretching potential. By using the short-time information of the LSC-IVR correlation function; however, it is shown how one can obtain more accurate results for the position of the stretching peak. Due to the physical decay in the condensed phase system, the LSC-IVR (LGA) is a good and practical approximate quantum approach for the IR spectrum of liquid water. The present results offer valuable insight into future attempts to improve the accuracy of the TTM3-F potential or other ab initio-based models in reproducing the IR spectrum of liquid water. [ABSTRACT FROM AUTHOR]

An efficient time-dependent importance sampling method is developed for the Monte Carlo calculation of time correlation functions via the initial value representation (IVR) of semiclassical (SC) theory. A prefactor-free time-dependent sampling function weights the importance of a trajectory based on the magnitude of its contribution to the time correlation function, and global trial moves are used to facilitate the efficient sampling the phase space of initial conditions. The method can be generally applied to sampling rare events efficiently while avoiding being trapped in a local region of the phase space. Results presented in the paper for two system-bath models demonstrate the efficiency of this new importance sampling method for full SC-IVR calculations. [ABSTRACT FROM AUTHOR]

Fewest-switches surface hopping (FSSH) is a popular nonadiabatic dynamics method which treats nuclei with classical mechanics and electrons with quantum mechanics. In order to simulate the motion of a wave packet as accurately as possible, standard FSSH requires a stochastic sampling of the trajectories over a distribution of initial conditions corresponding, e.g., to the Wigner distribution of the initial quantum wave packet. Although it is well-known that FSSH does not properly account for decoherence effects, there is some confusion in the literature about whether or not this averaging over a distribution of initial conditions can approximate some of the effects of decoherence. In this paper, we not only show that averaging over initial conditions does not generally account for decoherence, but also why it fails to do so. We also show how an apparent improvement in accuracy can be obtained for a fortuitous choice of model problems, even though this improvement is not possible, in general. For a basic set of one-dimensional and two-dimensional examples, we find significantly improved results using our recently introduced augmented FSSH algorithm. [ABSTRACT FROM AUTHOR]

Internal conversion is an inherently quantum mechanical process. To date, 'on the fly' computation of internal conversion rates is limited to harmonic approximations, which would seem to be especially unsuitable, given that the typical transition to the ground electronic state occurs at energies which are far from the harmonic limit. It is thus of interest to study the applicability of the semiclassial initial value representation (SCIVR) approach which is in principle amenable to on the fly studies even with 'many' degrees of freedom. In this paper we study the applicability of the Herman-Kluk (HK) SCIVR to a model system with two coupled and anharmonic degrees of freedom. We find that (a) the HK SCIVR is a good approximation to the exact quantum dynamics; (b) computation of the first order correction to the HK-SCIVR approximation corroborates the accuracy; (c) by studying a large parameter range, we find that the harmonic approximation is mostly unsatisfactory; and (d) for the specific model used, the coupling between the modes was found to be relatively unimportant. These results imply that the HK-SCIVR methodology is a good candidate for on the fly studies of internal conversion processes of 'large' molecules. [ABSTRACT FROM AUTHOR]

The frozen Gaussian approximation to the quantum propagator may be a viable method for obtaining 'on the fly' quantum dynamical information on systems with many degrees of freedom. However, it has two severe limitations, it rapidly loses normalization and one needs to know the Gaussian averaged potential, hence it is not a purely local theory in the force field. These limitations are in principle remedied by using the Herman-Kluk (HK) form for the semiclassical propagator. The HK propagator approximately conserves unitarity for relatively long times and depends only locally on the bare potential and its second derivatives. However, the HK propagator involves a much more expensive computation due to the need for evaluating the monodromy matrix elements. In this paper, we (a) derive a new formula for the normalization integral based on a prefactor free HK propagator which is amenable to 'on the fly' computations; (b) show that a frozen Gaussian version of the normalization integral is not readily computable 'on the fly'; (c) provide a new insight into how the HK prefactor leads to approximate unitarity; and (d) how one may construct a prefactor free approximation which combines the advantages of the frozen Gaussian and the HK propagators. The theoretical developments are backed by numerical examples on a Morse oscillator and a quartic double well potential. [ABSTRACT FROM AUTHOR]

In this paper we have calculated escape rate from a meta stable state in the presence of both colored internal thermal and external nonthermal noises. For the internal noise we have considered usual Gaussian distribution but the external noise may be Gaussian or non-Gaussian in characteristic. The calculated rate is valid for low noise strength of non-Gaussian noise such that an effective Gaussian approximation of non-Gaussian noise wherein the higher order even cumulants of order '4' and higher are neglected. The rate expression we derived here reduces to the known results of the literature, as well as for purely external noise driven activated rate process. The latter exhibits how the rate changes if one switches from non-Gaussian to Gaussian character of the external noise. [ABSTRACT FROM AUTHOR]

Initial value representations (IVRs) of semiclassical (SC) theory provide a general approach for adding quantum mechanical effects to classical molecular dynamics simulations of large molecular systems. Of the various versions of SC-IVR methodology for evaluating time correlation functions, the Fourier transform forward-backward (FB) approach is the simplest one that is able to describe true quantum coherence effects, so it is of considerable importance to find efficient and systematic ways for implementing it. It is shown in this paper that a Gaussian approximation for the “structure function”—the dependence of the correlation function on the (typically) momentum jump parameter—provides an efficient and accurate way for doing so. The approach is illustrated by an application to the time-dependent radial distribution function of I2 (after photoexcitation) in a cluster of (up to 16) argon atoms. [ABSTRACT FROM AUTHOR]

We present a significant improvement to a complex time-dependent WKB (CWKB) formulation developed by Boiron and Lombardi [J. Chem. Phys. 108, 3431 (1998)] in which the time-dependent WKB equations are solved along classical trajectories that propagate in complex space. Boiron and Lombardi showed that the method gives very good agreement with the exact quantum mechanical result as long as the wavefunction does not exhibit interference effects such as oscillations and nodes. In this paper, we show that this limitation can be overcome by superposing the contributions of crossing trajectories. Secondly, we demonstrate that the approximation improves when incorporating higher order terms in the expansion. Thirdly, equations of motion for caustics and Stokes lines are implemented to help overcome Stokes discontinuities. These improvements could make the CWKB formulation a competitive alternative to current time-dependent semiclassical methods. [ABSTRACT FROM AUTHOR]

The de Broglie-Bohm formulation of the Schrödinger equation implies conservation of the wave function probability density associated with each quantum trajectory in closed systems. This conservation property greatly simplifies numerical implementations of the quantum trajectory dynamics and increases its accuracy. The reconstruction of a wave function, however, becomes expensive or inaccurate as it requires fitting or interpolation procedures. In this paper we present a method of computing wave packet correlation functions and wave function projections, which typically contain all the desired information about dynamics, without the full knowledge of the wave function by making quadratic expansions of the wave function phase and amplitude near each trajectory similar to expansions used in semiclassical methods. Computation of the quantities of interest in this procedure is linear with respect to the number of trajectories. The introduced approximations are consistent with approximate quantum potential dynamics method. The projection technique is applied to model chemical systems and to the H+H2 exchange reaction in three dimensions. [ABSTRACT FROM AUTHOR]

Various versions of quantum transition state theory (QTST) for the calculation of chemical reaction rate coefficients have been developed. The Hansen-Andersen formulation of QTST approximates the reactive flux autocorrelation function of Miller, Schwartz, and Tromp for the problem of interest by the correlation function for a one-dimensional system with a suitably chosen parabolic reference potential. In this paper, we present a systematic method for correcting this approximation in order to obtain more accurate rate constants at low temperatures. Values of the correlation function for real or imaginary times, such as might be obtained from path integral Monte Carlo calculations, are used to construct an improved approximation for the flux autocorrelation function. Our method uses a Bayesian analysis that incorporates a modified form of the maximum entropy method. In contrast to other applications of the maximum entropy method to this problem, the improved approximation is consistent with the known analyticity properties of the correlation function. The method is applied to several test problems, including one based on the H[SUB2] + Cl reaction. Our method yields significant corrections using a remarkably small amount of information. [ABSTRACT FROM AUTHOR]

Polyexponential kinetical behavior typical for condensed phase reactions in highly viscous media is studied on a simple example of one-dimensional diffusion equation with a sink modeling a chemical conversion of reactants. The corresponding polyexponential regime is demonstrated to have a thorough analogy with the quasiclassical approximation of one-dimensional quantum mechanics and a relevant approximation for the Green’s function is developed. The asymptotic short- and long-time kinetics are examined at the analytical level. Contrary to the frozen medium approximation according to which the slow diffusion motion of the medium is entirely ignored, the present quasiclassical model is fit for a qualitative description of the total time interval covering the reaction events from the initial moment up to the ultimate steady-state monoexponential evolution. The range of validity of the quasiclassical approach is discussed. Numerical tests expose some peculiarities of the present treatment for equilibrium and nonequilibrium initial distributions. The work presents a qualitative development of the theory of nonexponential kinetics pioneered by papers of Agmon and Hopfield, Sumi and Marcus, and Nadler and Marcus. © 1995 American Institute of Physics. [ABSTRACT FROM AUTHOR]

The recently proposed multilayer multiconfiguration time-dependent Hartree (ML-MCTDH) approach to evaluating reactive quantum dynamics is applied to two model condensed-phase proton transfer reactions. The models consist of a one-dimensional double-well “system” that is bilinearly coupled to a “bath” of harmonic oscillators parameterized to represent a condensed-phase environment. Numerically exact quantum-mechanical flux correlation functions and thermal rate constants are obtained for a broad range of temperatures and system-bath coupling strengths, thus demonstrating the efficacy of the ML-MCTDH approach. Particular attention is focused on the regime where low temperatures are combined with weak system-bath coupling. Under such conditions it is found that long propagation times are often required and that quantum coherence effects may prevent a rigorous determination of the rate constant. [ABSTRACT FROM AUTHOR]

Mixed quantum/classical (MQC) simulations treat the majority of a system classically and reserve quantum mechanics only for a few degrees of freedom that actively participate in the chemical process(es) of interest. In MQC calculations, the quantum and classical degrees of freedom are coupled together using pseudopotentials. Although most pseudopotentials are developed empirically, there are methods for deriving pseudopotentials using the results of quantum chemistry calculations, which guarantee that the explicitly-treated valence electron wave functions remain orthogonal to the implicitly-treated core electron orbitals. Whether empirical or analytically derived in nature, to date all such pseudopotentials have been subject to the frozen core approximation (FCA) that ignores how changes in the nuclear coordinates alter the core orbitals, which in turn affects the wave function of the valence electrons. In this paper, we present a way to go beyond the FCA by developing pseudopotentials that respond to these changes. In other words, we show how to derive an analytic expression for a pseudopotential that is an explicit function of nuclear coordinates, thus accounting for the polarization effects experienced by atomic cores in different chemical environments. We then use this formalism to develop a coordinate-dependent pseudopotential for the bonding electron of the sodium dimer cation molecule and we show how the analytic representation of this potential can be used in one-electron MQC simulations that provide the accuracy of a fully quantum mechanical Hartree-Fock (HF) calculation at all internuclear separations. We also show that one-electron MQC simulations of Na2+ using our coordinate-dependent pseudopotential provide a significant advantage in accuracy compared to frozen core potentials with no additional computational expense. This is because use of a frozen core potential produces a charge density for the bonding electron of Na2+ that is too localized on the molecule, leading to significant overbinding of the valence electron. This means that FCA calculations are subject to inaccuracies of order ∼10% in the calculated bond length and vibrational frequency of the molecule relative to a full HF calculation; these errors are fully corrected by using our coordinate-dependent pseudopotential. Overall, our findings indicate that even for molecules like Na2+, which have a simple electronic structure that might be expected to be well-treated within the FCA, the importance of including the effects of the changing core molecular orbitals on the bonding electrons cannot be overlooked. [ABSTRACT FROM AUTHOR]

The quantum instanton approximation for thermal rate constants of chemical reactions [Miller, Zhao, Ceotto, and Yang, J. Chem. Phys. 119, 1329 (2003)], which is modeled after the earlier semiclassical instanton approach, is applied to the hydrogen abstraction reaction from methane by a hydrogen atom, H+CH[sub 4]→H[sub 2]+CH[sub 3], using a modified and recalibrated version of the Jordan–Gilbert potential surface. The quantum instanton rate is evaluated using path integral Monte Carlo approaches based on the recently proposed implementation schemes [Yamamoto and Miller, J. Chem. Phys. 120, 3086 (2004)]. The calculations were carried out using the Cartesian coordinates of all the atoms (thus involving 18 degrees of freedom), thereby taking explicit account of rotational effects of the whole system and also allowing the equivalent treatment of the four methane hydrogens. To achieve such a treatment, we present extended forms of the path integral estimators for relevant quantities that may be used for general N-atom systems with any generalized reaction coordinates. The quantum instanton rates thus obtained for the temperature range T=200–2000 K show good agreement with available experimental data, which gives support to the accuracy of the underlying potential surface used. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]

In this paper, we investigate coherent quantum dynamics in a nonequilibrium environment. We focus on a two-state quantum system strongly coupled to a single classical environmental oscillator, and explore the effect of nonstationary statistical properties of the oscillator on the quantum evolution. A simple nonequilibrium model, consisting of an oscillator with a well-defined initial phase which undergoes subsequent diffusion, is introduced and studied. Approximate but accurate analytic expressions for the evolution of the off-diagonal density matrix element of the quantum system are derived in the second-order cumulant approximation. The effect of the initial phase choice on the subsequent quantum evolution is quantified. It is observed that the initial phase can have a significant effect on the preservation of coherence on short time scales, suggesting this variable as a control parameter for optimizing coherence in many-body quantum systems. [ABSTRACT FROM AUTHOR]

In this paper we suggest an approximation to the Wigner Path integral using quantum trajectories. We discuss the results in terms of Bohm’s theory and Takabayasi’s study of the relation between Bohm and Wigner theory. We discuss how this new approach is an approximation to Bohm’s theory and discuss its validity. We also suggest a method of evaluation and where such an approach would be useful. [ABSTRACT FROM AUTHOR]

Quantum molecular-dynamics simulations are an important way of gaining information on the molecular level about chemical systems. In this paper, a new method for the approximate solution of the time-dependent Schro¨dinger equation is presented. This is a reformulation of the multiconfiguration time-dependent Hartree (MCTDH) wave packet propagation method, which is transformed so that the evolution of the wave function can be represented by pseudo-particle trajectories. In this way, the poor scaling of computational resources with system size attending all exact solutions of the time-dependent Schro¨dinger equation are circumvented. The equations of motion for the trajectories and the wave function expansion coefficients (importance of each trajectory for the representation) are derived using a variational principle. Other than the MCTDH ansatz, no major approximations have been introduced, and the method converges on the numerically exact solution. Importantly, the trajectories are not classical trajectories, and are coupled by nonlocal effects. A strategy for the practical solution of the equations of motion is then detailed. © 2001 American Institute of Physics. [ABSTRACT FROM AUTHOR]

We apply the matrix product state (MPS) method to propagate the hierarchical equations of motion (HEOM). It is shown that the MPS approximation works well in different type of problems, including boson and fermion baths. The MPS method based on the time-dependent variational principle is also found to be applicable to HEOM with over one thousand effective modes. Combining the flexibility of the HEOM in defining the effective modes and the efficiency of the MPS method thus may provide a promising tool in simulating quantum dynamics in condensed phases. [ABSTRACT FROM AUTHOR]

In our earlier paper we argued that for the decomposition (1) to be valid, the operator V, which we denoted as Helec, must be written as a direct integral over clamped-nuclei electronic Hamiltonians, the integral being taken over all nuclear positions; we further suggested that after this correction an expansion analogous to (5) would be problematic. The operator Helec has purely continuous spectrum extending from some minimum value to on the real-axis; it has no true eigenvalues, and no normalizable eigenvectors in the Hilbert space L2(x, X). We emphasize that this description is entirely in agreement with the mathematical physics literature but not with the conventional Born-Huang discussion summarized above, since the distinction between Ho (Xf ) and its direct integral Helec is obvious and fundamental. The recent work of Jecko7 gives a critical mathematical summary of the Born-Oppenheimer approximation including the expansion approach described by (3).(7), and emphasizes the fundamental role of spectral projection for operators with continuous spectra. Thus, it is very much to the point that there is no complete set of eigenfunctions in Hilbert space for the clamped-nuclei Hamiltonian Ho(Xf ) required for the expansion (5) (and likewise for Helec). This is crucial since the Born-Huang method effectively writes the clamped nuclei Hamiltonian in a spectral form using the geigenfunctions to provide a resolution of the identity [ABSTRACT FROM AUTHOR]

Closed-form analytic solutions to non-secular Bloch-Redfield master equations for quantum dynamics of a V-type system driven by weak coupling to a thermal bath, relevant to light harvesting processes, are obtained and discussed. We focus on noise-induced Fano coherences among the excited states induced by incoherent driving of the V-system initially in the ground state. For suddenly turned-on incoherent driving, the time evolution of the coherences is determined by the damping parameter ζ = 1/2(γ1 + γ2)/Δp, where γi are the radiative decay rates of the excited levels i = 1, 2, and Δp =Δ² + (1 - p²)γ1γ2 depends on the excited-state level splitting Δ > 0 and the angle between the transition dipole moments in the energy basis. The coherences oscillate as a function of time in the underdamped limit (ζ ≫ 1), approach a long-lived quasi-steady state in the overdamped limit (ζ ≪ 1), and display an intermediate behavior at critical damping (ζ = 1). The sudden incoherent turn-on is shown to generate a mixture of excited eigenstates |e1> and |e2> and their in-phase coherent superposition ..., which is remarkably long-lived in the overdamped limit (where r¹ and r² are the incoherent pumping rates). Formation of this coherent superposition enhances the decay rate from the excited states to the ground state. In the strongly asymmetric V-system where the coupling strengths between the ground state and the excited states differ significantly, additional asymptotic quasistationary coherences are identified, which arise due to slow equilibration of one of the excited states. Finally, we demonstrate that noise-induced Fano coherences are maximized with respect to populations when r1 = r2 and the transition dipole moments are fully aligned. [ABSTRACT FROM AUTHOR]

We present a formalism that explicitly unifies the commonly used Nakajima-Zwanzig approach for reduced density matrix dynamics with the more versatile Mori theory in the context of nonequilibrium dynamics. Employing a Dyson-type expansion to circumvent the difficulty of projected dynamics, we obtain a self-consistent equation for the memory kernel which requires only knowledge of normally evolved auxiliary kernels. To illustrate the properties of the current approach, we focus on the spin-boson model and limit our attention to the use of a simple and inexpensive quasi-classical dynamics, given by the Ehrenfest method, for the calculation of the auxiliary kernels. For the first time, we provide a detailed analysis of the dependence of the properties of the memory kernels obtained via different projection operators, namely, the thermal (Redfield-type) and population based (NIBA-type) projection operators. We further elucidate the conditions that lead to short-lived memory kernels and the regions of parameter space to which this program is best suited. Via a thorough analysis of the different closures available for the auxiliary kernels and the convergence properties of the self-consistently extracted memory kernel, we identify the mechanisms whereby the current approach leads to a significant improvement over the direct usage of standard semi- and quasi-classical dynamics. [ABSTRACT FROM AUTHOR]

Generalized quantum master equations (GQMEs) are an important tool in modeling chemical and physical processes. For a large number of problems, it has been shown that exact and approximate quantum dynamics methods can be made dramatically more efficient, and in the latter case more accurate, by proceeding via the GQME formalism. However, there are many situations where utilizing the GQME approach with an approximate method has been observed to return the same dynamics as using that method directly. Here, for systems both in and out of equilibrium, we provide a more detailed understanding of the conditions under which using an approximate method can yield benefits when combined with the GQME formalism. In particular, we demonstrate the necessary manipulations, which are satisfied by exact quantum dynamics, that are required to recast the memory kernel in a form that can be analytically shown to yield the same result as a direct application of the dynamics regardless of the approximation used. By considering the connections between these forms of the kernel, we derive the conditions that approximate methods must satisfy if they are to offer different results when used in conjunction with the GQME formalism. These analytical results thus provide new insights as to when proceeding via the GQME approach can be used to improve the accuracy of simulations. [ABSTRACT FROM AUTHOR]

Well-defined criteria are proposed for assessing the accuracy of quantum master equations whose memory functions are approximated by Padé resummation of the first two moments in the electronic coupling. These criteria partition the parameter space into distinct levels of expected accuracy, ranging from quantitatively accurate regimes to regions of parameter space where the approach is not expected to be applicable. Extensive comparison of Padé-resummed master equations with numerically exact results in the context of the spin-boson model demonstrates that the proposed criteria correctly demarcate the regions of parameter space where the Padé approximation is reliable. The applicability analysis we present is not confined to the specifics of the Hamiltonian under consideration and should provide guidelines for other classes of resummation techniques. [ABSTRACT FROM AUTHOR]

We show that a single change in the derivation of the linearized semiclassical-initial value representation (LSC-IVR or "classical Wigner approximation") results in a classical dynamics which conserves the quantum Boltzmann distribution. We rederive the (standard) LSC-IVR approach by writing the (exact) quantum time-correlation function in terms of the normal modes of a free ring-polymer (i.e., a discrete imaginary-time Feynman path), taking the limit that the number of polymer beads N ↑ ∞, such that the lowest normal-mode frequencies take their "Matsubara" values. The change we propose is to truncate the quantum Liouvillian, not explicitly in powers of ħ² at ħ0 (which gives back the standard LSC-IVR approximation), but in the normal-mode derivatives corresponding to the lowest Matsubara frequencies. The resulting "Matsubara" dynamics is inherently classical (since all terms O(ħ²) disappear from the Matsubara Liouvillian in the limit N ↑ ∞) and conserves the quantum Boltzmann distribution because the Matsubara Hamiltonian is symmetric with respect to imaginary-time translation. Numerical tests show that the Matsubara approximation to the quantum time-correlation function converges with respect to the number of modes and gives better agreement than LSC-IVR with the exact quantum result. Matsubara dynamics is too computationally expensive to be applied to complex systems, but its further approximation may lead to practical methods. [ABSTRACT FROM AUTHOR]

An analytical solution for the electrostatic energy and potential for a system of charged, polarizable spheres in a dielectric medium is developed from a multiple scattering expansion that is equivalent to a formal solution to Poisson's equation for the system. The leading contributions emerge in the form of effective two-, three-, and four-body interactions that are explicit analytical functions of the sphere positions, charges, and internal dielectric constants and the external dielectric constant, thereby also enabling analytical computation of the electrostatic forces on the ions. Tests of successive terms demonstrate their rapid convergence. Similar methods can be used to evaluate higher order contributions and the expansion for the electrostatic field. The results will prove far more efficient for MD and MC simulations with spherical particles than current approximate methods that require the computation of surface polarization charge distributions but that apply also for systems with complex geometries. [ABSTRACT FROM AUTHOR]

A frozen-rotor approximation is formulated for the mixed quantum/classical theory of collisional energy transfer and ro-vibrational energy flow [M. Ivanov and D. Babikov, J. Chem. Phys. 134, 144107 (2011)]. Numerical tests are conducted to assess its efficiency and accuracy, compared to the original version of the method, where rotation of the molecule in space is treated explicitly and adiabatically. New approach is considerably faster and helps blocking the artificial ro-vibrational transitions at the pre- and post-collisional stages of the process. Although molecular orientation in space is fixed, the energy exchange between rotational, vibrational, and translational digresses of freedom still occurs, allowing to compute ro-vibrational excitation and quenching. Behavior of the energy transfer function through eight orders of magnitude range of values and in a broad range of ΔE is reproduced well. In the range of moderate -500 <= ΔE <= +500 cm-1 the approximate method is rather accurate. The absolute values of stabilization cross sections for scattering resonances trapped behind the centrifugal threshold are a factor 2-to-3 smaller (compared to the explicit-rotation approach). This performance is acceptable and similar to the well-known sudden-rotation approximation in the time-independent inelastic scattering methods. [ABSTRACT FROM AUTHOR]

The semiclassical instanton approach discussed by Kryvohuz [J. Chem. Phys. 134, 114103 (2011)] is applied to calculate kinetic H/D isotope effect (KIE) of intramolecular hydrogen transfer in cis-1,3-pentadiene. All 33 vibrational degrees of freedom are treated quantum mechanically with semiclassical approximation. Nuclear quantum effects such as tunneling under the barrier and zero-point energy are automatically incorporated in the theory, and are shown to be responsible for the observed appreciable kinetic isotope effect in cis-1,3-pentadiene. Over the barrier passage is also automatically included. Numerical calculations are performed on an empirical valence bond potential energy surface and compared with the previous experimental and theoretical studies. An estimation of heavy-atom 12C/13C KIE in the same system is also provided and the factors contributing to it are discussed. [ABSTRACT FROM AUTHOR]

A time-averaging approximation method developed to efficiently calculate the short-time dynamics of coupled vibrational chromophores using mixed quantum/classical theories is extended in order to be applicable to the study of vibrational dynamics at longer time scales. A quantum mechanical time propagator for long times is decomposed into the product of short-time propagators, and a time-averaging approximation is then applied to each of the latter. Using the extended time-averaging approximation, we calculate the anisotropy decay of the data obtained from impulsive vibrational pump-probe experiments on the OH stretching modes of water, which is in excellent agreement with numerically exact results. [ABSTRACT FROM AUTHOR]

The collision-induced process He+H2+(v=0-2;j=0-3)→He+H+H+ has been investigated using a time-dependent quantum mechanical wave packet approach, within the centrifugal sudden approximation. The exchange reaction He + H2+ → HeH+ + H, which has a lower threshold, dominates over the dissociation process over the entire energy range considered in this study. The reaction cross section for both the exchange and dissociation channels and the branching ratio between the two channels have been computed on the McLaughlin-Thompson-Joseph-Sathyamurthy potential-energy surface and compared with the available experimental and quasiclassical trajectory results. [ABSTRACT FROM AUTHOR]

In the study of open quantum systems, the polaron transformation has recently attracted a renewed interest as it offers the possibility to explore the strong system-bath coupling regime. Despite this interest, a clear and unambiguous analysis of the regimes of validity of the polaron transformation is still lacking. Here we provide such a benchmark, comparing second order perturbation theory results in the original untransformed frame, the polaron frame, and the variational extension with numerically exact path integral calculations of the equilibrium reduced density matrix. Equilibrium quantities allow a direct comparison of the three methods without invoking any further approximations as is usually required in deriving master equations. It is found that the second order results in the original frame are accurate for weak system-bath coupling; the results deteriorate when the bath cut-off frequency decreases. The full polaron results are accurate for the entire range of coupling for a fast bath but only in the strong coupling regime for a slow bath. The variational method is capable of interpolating between these two methods and is valid over a much broader range of parameters. [ABSTRACT FROM AUTHOR]

The exchange processes of D + H2O and D + HOD reactions are studied using initial state-selected time-dependent wave packet approach in full dimension. The total reaction probabilities for different partial waves, together with the integral cross sections, are obtained both by the centrifugal sudden (CS) approximation and exact coupled-channel (CC) calculations, for the H2O(HOD) reactant initially in the ground rovibrational state. In the CC calculations, small resonance peaks in the reaction probabilities and quick diminishing of the resonance peaks with the increase of total angular momenta J do not lead to clear step-like features just above the threshold in the cross sections for the title reactions, which are different in other isotopically substituted reactions where the hydrogen atom was included as the reactant instead of the deuterium atom [B. Fu, Y. Zhou, and D. H. Zhang, Chem. Sci. 3, 270 (2012); B. Fu and D. H. Zhang, J. Phys. Chem. A 116, 820 (2012)]. It is interesting that the shape resonance-induced features resulting from the reaction tunneling are significantly diminished accordingly in the reactions of the deuterium atom and H2O or HOD, owing to the weaker tunneling capability of the reagent deuterium atom in the title reactions than the reagent hydrogen atom in other reactions. In the CS calculations, the resonance peaks persist in many partial waves but cannot survive the partial-wave summations. The cross sections for the D′ + H2O → D′OH + H and D′ + HOD → D′OD + H reactions are substantially larger than those for the D′ + HOD → HOD′ + D reaction, indicating that the D′/H exchange reactions are much more favored than the D′/D exchange. [ABSTRACT FROM AUTHOR]

We propose to measure nonadiabaticity of molecular quantum dynamics rigorously with the quantum fidelity between the Born-Oppenheimer and fully nonadiabatic dynamics. It is shown that this measure of nonadiabaticity applies in situations where other criteria, such as the energy gap criterion or the extent of population transfer, fail. We further propose to estimate this quantum fidelity efficiently with a generalization of the dephasing representation to multiple surfaces. Two variants of the multiple-surface dephasing representation (MSDR) are introduced, in which the nuclei are propagated either with the fewest-switches surface hopping or with the locally mean field dynamics (LMFD). The LMFD can be interpreted as the Ehrenfest dynamics of an ensemble of nuclear trajectories, and has been used previously in the nonadiabatic semiclassical initial value representation. In addition to propagating an ensemble of classical trajectories, the MSDR requires evaluating nonadiabatic couplings and solving the Schrödinger (or more generally, the quantum Liouville-von Neumann) equation for a single discrete degree of freedom. The MSDR can be also used in the diabatic basis to measure the importance of the diabatic couplings. The method is tested on three model problems introduced by Tully and on a two-surface model of dissociation of NaI. [ABSTRACT FROM AUTHOR]

Cross sections and rate coefficients for low lying rotational transitions in H2O colliding with para-hydrogen pH2 are computed using an adiabatic approximation which reduces the dimensional dynamics from a 5D to a 3D problem. Calculations have been performed at the close-coupling level using the recent potential of Valiron et al. [J. Chem. Phys. 129, 134306 (2008)]. A good agreement is found between the reduced adiabatic calculations and the 5D exact calculations, with an impressive time saving and memory gain. This adiabatic reduction of dimensionality seems very promising for scattering studies involving the excitation of a heavy target molecule by a light molecular projectile. [ABSTRACT FROM AUTHOR]

We make the first application of semiclassical (SC) techniques to the plane-wavepacket formulation of time-domain (T-domain) scattering. The angular scattering of the state-to-state reaction, H + D2(vi = 0, ji = 0) → HD(vf = 3, jf = 0) + D, is analysed, where v and j are vibrational and rotational quantum numbers, respectively. It is proved that the forward-angle scattering in the T-domain, which arises from a delayed mechanism, is an example of a glory. The SC techniques used in the T-domain are: An integral transitional approximation, a semiclassical transitional approximation, a uniform semiclassical approximation (USA), a primitive semiclassical approximation and a classical semiclassical approximation. Nearside-farside (NF) scattering theory is also employed, both partial wave and SC, since a NF analysis provides valuable insights into oscillatory structures present in the full scattering pattern. In addition, we incorporate techniques into the SC theory called 'one linear fit' and 'two linear fits', which allow the derivative of the quantum deflection function, Θ′(J), to be estimated when ΘJ exhibits undulations as a function of J, the total angular momentum variable. The input to our SC analyses is numerical scattering (S) matrix data, calculated from accurate quantum collisional calculations for the Boothroyd-Keogh-Martin-Peterson potential energy surface No. 2, in the energy domain (E-domain), from which accurate S matrix elements in the T-domain are generated. In the E-domain, we introduce a new technique, called 'T-to-E domain SC analysis.' It half-Fourier transforms the E-domain accurate quantum scattering amplitude to the T-domain, where we carry out a SC analysis; this is followed by an inverse half-Fourier transform of the T-domain SC scattering amplitude back to the E-domain. We demonstrate that T-to-E USA differential cross sections (DCSs) agree well with exact quantum DCSs at forward angles, for energies where a direct USA analysis in the E-domain fails. [ABSTRACT FROM AUTHOR]

A robust interpolation between the weak- and strong-correlation regimes of quantum systems is presented. It affords approximants to the function E(ω) describing the dependence of the total energy (or other observable) on the coupling parameter ω that measures the correlation strength. The approximants conform to truncations of the asymptotic expansions of E(ω) at the ω → 0 and ω → ∞ limits with arbitrary (but given) numbers of terms. In addition, depending on the number of fitted parameters, they either reproduce or optimally (in the least-square or maximum-error sense) approximate the exact E(ω) at any given number of values of the coupling strength. Numerical tests demonstrate the high accuracy of even the low-order approximate expression for E(ω). The approximants, which do not suffer from spurious poles, possess a wide range of applicability that stems from their capability of accurately reproducing not only E(ω) but also its derivatives with respect to ω. They are equally useful for interpolation between the low- and high-temperature limits of energy and other quantities associated with various models of statistical thermodynamics. The new interpolation scheme is not applicable to the cases where the weak- and strong-correlation asymptotics involve non-analytic functions of ω or expressions dependent on logarithm of the coupling strength. Excluded are also the cases where the weak- and strong-correlation asymptotics pertain to de facto different states, e.g., the ground state of a homogeneous electron gas in three dimensions. [ABSTRACT FROM AUTHOR]

We report fully-quantum time-independent calculations of cross sections and rate coefficients for the collisional excitation and dissociation of D2 by H, two astrophysically relevant processes. Our calculations are based on the recent H3 global potential energy surface of Mielke et al. [J. Chem. Phys. 116, 4142 (2002)]. Results of exact three-dimensional calculations, i.e., including the reactive channels, are compared to pure inelastic two-dimensional calculations based on the rigid rotor approximation. A reasonable agreement is found between the two sets of inelastic cross sections over the whole energy range 10-9000 cm-1. At the highest collisional energies, where the reactive channels are significant, the rigid rotor approach slightly overestimates the cross sections, as expected. At moderate collisional energies, however, the opposite behaviour is observed. The rigid rotor approach is found to be reliable at temperatures below ∼500 K, with a significant but moderate contribution from reactive channels [ABSTRACT FROM AUTHOR]

We present a new approach to calculate real-time quantum dynamics in complex systems. The formalism is based on the partitioning of a system's environment into 'core' and 'reservoir' modes with the former to be treated quantum mechanically and the latter classically. The presented method only requires the calculation of the system's reduced density matrix averaged over the quantum core degrees of freedom which is then coupled to a classically evolved reservoir to treat the remaining modes. We demonstrate our approach by applying it to the spin-boson problem using the noninteracting blip approximation to treat the system and core, and Ehrenfest dynamics to treat the reservoir. The resulting hybrid methodology is accurate for both fast and slow baths, since it naturally reduces to its composite methods in their respective regimes of validity. In addition, our combined method is shown to yield good results in intermediate regimes where neither approximation alone is accurate and to perform equally well for both strong and weak system-bath coupling. Our approach therefore provides an accurate and efficient methodology for calculating quantum dynamics in complex systems. [ABSTRACT FROM AUTHOR]

The linearized semiclassical initial value representation is employed to describe ultrafast electron transfer processes coupled to a phonon bath and weakly coupled to a proton mode. The goal of our theoretical investigation is to understand the influence of the proton on the electronic dynamics in various bath relaxation regimes. More specifically, we study the impact of the proton on coherences and analyze if the coupling to the proton is revealed in the form of an isotope effect. This will be important in distinguishing reactions in which the proton does not undergo significant rearrangement from those in which the electron transfer is accompanied by proton transfer. Unlike other methodologies widely employed to describe nonadiabatic electron transfer, this approach treats the electronic and nuclear degrees of freedom consistently. However, due to the linearized approximation, quantum interference effects are not captured accurately. Our study shows that at small phonon bath reorganization energies, coherent oscillations and isotope effect are observed in both slow and fast bath regimes. The coherences are more substantially damped by deuterium in comparison to the proton. Further, in contrast to the dynamics of the spin-boson model, the coherences are not long-lived. At large bath reorganization energies, the decay is incoherent in the slow and fast bath regimes. In this case, the extent of the isotope effect depends on the relative relaxation timescales of the proton mode and the phonon bath. The isotope effect is magnified for baths that relax on picosecond timescales in contrast to baths that relax in femtoseconds. [ABSTRACT FROM AUTHOR]

We focus on the non-equilibrium two-bath spin-boson model, a toy model for examining quantum thermal transport in many-body open systems. Describing the dynamics within the noninteracting-blip approximation equations, applicable, e.g., in the strong system-bath coupling limit and/or at high temperatures, we derive expressions for the cumulant generating function in both the Markovian and non-Markovian limits by energy-resolving the quantum master equation of the subsystem. For a Markovian bath, we readily demonstrate the validity of a steady-state heat exchange fluctuation theorem. In the non-Markovian limit a 'weaker' symmetry relation generally holds, a general outcome of microreversibility. We discuss the reduction of this symmetry relation to the universal steady-state fluctuation theorem. Using the cumulant generating function, an analytic expression for the heat current is obtained. Our results establish the validity of the steady-state heat exchange fluctuation theorem in quantum systems with strong system-bath interactions. From the practical point of view, this study provides tools for exploring transport characteristics of the two-bath spin-boson model, a prototype for a nonlinear thermal conductor. [ABSTRACT FROM AUTHOR]

Given the exponentially scaling cost of full quantum calculations, approximations need to be employed for the simulation of the time evolution of chemical systems. We present a modified version of surface hopping that has the potential to treat larger systems. This is accomplished through an Ehrenfest-like treatment of the excited states, thereby reducing the dynamics to transitions between the ground state and a mean-field excited state. A simplified description of the excited states is achieved, while still allowing for an accurate description of disparate reaction channels. We test our mean-field approximation for the excited states on a series of model problems. Results are compared to the standard surface hopping procedure, with its explicit treatment of all excited states, and the traditional Ehrenfest approach, with its averaging together of all states. [ABSTRACT FROM AUTHOR]

We present a method which uses the results of a molecular Kohn-Sham calculation at a reference geometry to approximate the energy at many different geometries. The Kohn-Sham electron density of the reference geometry is decomposed into atomic fragments, which move with the nuclei to approximate the density at a new geometry and the energy is evaluated with the Harris-Foulkes functional. Preliminary results for a biological quantum-mechanics/molecular-mechanics trajectory are promising: the errors of reference-geometry Harris-Foulkes (compared to full self-consistent Kohn-Sham) for the PBE exchange-correlation functional have the same magnitude as the difference between the energies of PBE and BLYP. [ABSTRACT FROM AUTHOR]

For a rigorous quantum simulation of nonadiabatic dynamics of electrons and nuclei, knowledge of not only the first-order but also the second-order nonadiabatic couplings (NACs) is required. Here, we propose a method to efficiently calculate the second-order NAC from time-dependent density functional theory (TDDFT), on the basis of the Casida ansatz adapted for the computation of first-order NAC, which has been justified in our previous work and can be shown to be valid for calculating second-order NAC between ground state and singly excited states within the Tamm-Dancoff approximation. Test calculations of the second-order NAC in the immediate vicinity of Jahn-Teller and Renner-Teller intersections show that calculation results from TDDFT, combined with modified linear response theory, agree well with the prediction from the Jahn-Teller/Renner-Teller models. Contrary to the diverging behavior of the first-order NAC near all types of intersection points, the Cartesian components of the second-order NAC are shown to be negligibly small near Renner-Teller glancing intersections, while they are significantly large near the Jahn-Teller conical intersections. Nevertheless, the components of the second-order NAC can cancel each other to a large extent in Jahn-Teller systems, indicating the background of neglecting the second-order NAC in practical dynamics simulations. On the other hand, it is shown that such a cancellation becomes less effective in an elliptic Jahn-Teller system and thus the role of second-order NAC needs to be evaluated in the rigorous framework. Our study shows that TDDFT is promising to provide accurate data of NAC for full quantum mechanical simulation of nonadiabatic processes. [ABSTRACT FROM AUTHOR]