Optical linear response function with linear and diagonal quadratic electron-vibration coupling in mixed quantum-classical systems

Optical linear response function of linearly and quadratically coupled mixed quantum-classical condensed phase systems is derived. The linear response function is derived using Kapral's formalism of statistical mechanics in mixed quantum-classical systems. Our mixed quantum-classical linear dipole moment correlation function J(t) is compared with the full quantum J(t) [Y. J. Yan and S. Mukamel, J. Chem. Phys. 85, 5908 (1986)] in the high temperature limit. Model calculations and discussion of our results are presented. Various formulas of Franck-Condon factors for both linear and quadratic coupling are discussed.     

Mixed quantum-classical description of spectroscopy of dissipative systems

Mixed quantum-classical statistical mechanics is employed to calculate dipole moment correlation function and linear absorption spectra. A quantum two-level subsystem interacting with quantum vibrations (primary oscillators) which in turn are coupled to a classical bath composed of infinite set of harmonic oscillators is used as a dissipative system. Starting with mixed quantum-classical Liouville equation for the evaluation of the mixed quantum-classical dipole moment correlation function and using coherent states and the inverse of Baker-Campbell-Hausdorf formula to evaluate the trace over the primary oscillators, whereby, a closed analytical expression for the electronic dipole moment correlation function is obtained. Illustrations of several absorption spectra at different temperatures are provided. An approximate optical four-point correlation is obtained in the high temperature limit. A strategy for deriving an exact optical four-point correlation is suggested.     

Anharmonic nuclear dynamics in the mixed quantum-classical limit

This study employs mixed quantum-classical dynamics (MQCD) formalism to evaluate the linear electronic dipole moment time correlation function (DMTCF) in which a Morse oscillator serves to model the associated vibrations in a mixed quantum-classical (MQC) environment. While the main purpose of this work is to study the applicability of MQCD formalism to anharmonic systems in condensed phase, approximate schemes to physically evaluate the mathematically divergent integrals have been developed in order to deal with the essential singularities that arise while evaluating the Morse oscillator canonical partition function and the DMTCF in MQC systems in the classical limit. The motivation for numerically and analytically evaluating these divergent integrals is that a partition function of any system should lead to a finite value at any temperature and therefore this divergence is unphysical. Additionally, since a partition function is to signify the number of accessible states to the system at hand, divergent results are not physically acceptable. As such, straightforward approximate analytic expressions, at different levels of rigor, for both the classical Morse oscillator partition function and the DMTCF in MQC systems are derived, for the first time. Calculations of Morse oscillator partition function values using different approaches at various temperatures for CO, HCl, and I(2) molecules, showing good results, are presented to test the expressions derived herein. It is found that this divergence, due to singularity, diminishes upon lowering the temperature and only arises at high temperatures. The gradual diminishing of the singularity upon lowering the temperature is sensible since the Morse potential fits the parabolic potential at low temperatures. Model calculations and discussion of the DMTCF and linear absorption spectra in MQC systems using the molecular constants of CO molecule are provided. The linear absorption lineshape is derived by two methods, one of which is asymptotic expansion.     

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.     

Mixed quantum-classical Redfield master equation

Redfield master equation is derived from mixed quantum-classical Liouville equation using product initial conditions. Simple two-level system example is given and comparison with Fermi golden rule is made.     

Statistical mechanics of quantum-classical systems with holonomic constraints

The statistical mechanics of quantum-classical systems with holonomic constraints is formulated rigorously by unifying the classical Dirac bracket and the quantum-classical bracket in matrix form. The resulting Dirac quantum-classical theory, which conserves the holonomic constraints exactly, is then used to formulate time evolution and statistical mechanics. The correct momentum-jump approximation for constrained systems arises naturally from this formalism. Finally, in analogy with what was found in the classical case, it is shown that the rigorous linear-response function of constrained quantum-classical systems contains nontrivial additional terms which are absent in the response of unconstrained systems.     

Calculations of scattering of N2 molecules from Ru(0001)

Recent experimental measurements of state resolved scattering of nitrogen molecules from a Ru(0001) surface are discussed in comparison with a mixed quantum-classical theory that has given reasonable explanations for similar data on other systems. Acceptable agreement between data and calculations is obtained, but only upon assuming an effective mass of the surface equal to 2.3 times the mass of a single Ru atom.     

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.     

Extended hydrodynamic approach to quantum-classical nonequilibrium evolution. I. Theory

A mixed quantum-classical formulation is developed for a quantum subsystem in strong interaction with an N-particle environment, to be treated as classical in the framework of a hydrodynamic representation. Starting from the quantum Liouville equation for the N-particle distribution and the corresponding reduced single-particle distribution, exact quantum hydrodynamic equations are obtained for the momentum moments of the single-particle distribution coupled to a discretized quantum subsystem. The quantum-classical limit is subsequently taken and the resulting hierarchy of equations is further approximated by various closure schemes. These include, in particular, (i) a Grad-Hermite-type closure, (ii) a Gaussian closure at the level of a quantum-classical local Maxwellian distribution, and (iii) a dynamical density functional theory approximation by which the hydrodynamic pressure term is replaced by a free energy functional derivative. The latter limit yields a mixed quantum-classical formulation which has previously been introduced by I. Burghardt and B. Bagchi, Chem. Phys. 134, 343 (2006).     

多原子体系非绝热动力学的近似理论方法

非绝热动力学普遍存在于光物理和光化学过程中。描述非绝热跃迁需要处理电子-原子核间的相互耦合运动。由于计算量随体系尺度增大剧烈增长,准确的量子动力学计算目前只适用于描述小分子体系。为了研究多原子分子体系的非绝热过程,近年来发展了一些基于量子-经典动力学近似方法。本文将对典型的这类方法包括经典Ehrenfest方法、面跳跃方法、基于Wigner表示的混合量子-经典方法进行简要的介绍,并讨论如何将量子-经典动力学方法与电子结构从头算手段相结合,模拟非绝热过程。重点阐明各种方法的基本思想和优缺点,并对该领域的发展进行展望。     

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.     

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.     

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

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