Semiclassical description of quantum coherence effects and their quenching: A forward–backward initial value representation study

Thermo-Field Dynamics Approach to Photo-induced Electronic Transitions Driven by Incoherent Thermal Radiation

The effects of thermal light-matter interaction on the dynamics of photo-induced electronic transitions in molecules are investigated using a novel first principles approach based on the thermo-field dynamics description of both the molecular vibrational modes and of the radiation field. The developed approach permits numerically accurate simulations of quantum dynamics of electronic/excitonic systems coupled to nuclear and photonic baths kept at different temperatures. The baths can be described by arbitrary spectral densities and can have any system-bath coupling strengths. In agreement with the results obtained previously by less rigorous methods, we show that the excitation process obtained by the continuous interaction with the suddenly turned-on thermal radiation field creates a mixed ensemble having a nonnegligible component consisting of a superposition of vibronic eigenstates which can sustain coherent oscillations for relatively long times. The results become especially relevant for the dynamics of electronic transitions upon sunlight excitation. Analytical results based on time-dependent perturbation theory support the numerical simulations and provide a simple interpretation of the time evolution of quantum observables.

Mixed Quantum-Classical Liouville Equation Treatment of Electronic Spectroscopy of Condensed Systems: Harmonic and Anharmonic Electron-Phonon Coupling

This Review integrates the use of electronic optical response function theory and the mixed quantum-classical (MQC) Liouville equation (MQCLE), thereby leading to electronic spectroscopy in MQC media. It further sheds light on the applicability, utility, and efficiency of the mixed quantum-classical dynamics (MQCD) formalism, which starts off with the MQCLE, in probing spectroscopy and dynamics of condensed systems, whereby quantum mechanics and classical mechanics are combined systematically. The author has been exploring and implementing MQCD to investigate electron-phonon coupling effects on electronic dephasing in harmonic and anharmonic systems by calculating linear and nonlinear optical transition analytically and numerically dipole moment time correlation functions in an MQC environment, thereby presenting an in depth spectral profile analysis and their shape and symmetry. The distinctive capability of the MQC time correlation functions is that ergodicity and stationarity properties are inherently satisfied as part of the mixed quantum-classical dynamics (MQCD) framework, unlike classical correlation functions. While some research groups have applied MQCLE to calculate vibrational spectra to study hydrogen-bonded complexes in a MQC environment and other groups calculated Optical response function to probe electron transfer dynamics using the basis mapping technique, the approach, purpose, rigor, applications, and path to the end results reported herein are different. Finally, the same framework is employed to study dissipative systems in the MQC limit, whereby the zero-phonon line adopts the correct width and eliminates its asymmetry. While the full quantum mechanical model, like the multimode Brownian oscillator (MBO) model, yields the correct width and inaccurate shape in the low-temperature limit, the MQCD formalism seems to produce an accurate zero-phonon profile. Nonlinear optical signals are also reviewed in MQC media to show the applicability and utility of this approach. The vibronic optical response functions developed here will account for geometry change, frequency change, and anharmonicity upon electronic excitation to accurately probe electronic dephasing, electron-phonon coupling, shape, and symmetry of profiles and present differences and similarities to the MBO model on pure electronic dephasing. Frequency change and anharmonicity are vitally crucial for accurately assessing electron-phonon coupling upon electronic excitation. This is an additional unique result obtained by the author to further demonstrate the applicability and utility of this approach over other approximation schemes in probing electronic dephasing including that of the MBO model.

Electronic dephasing of polyatomic molecules interacting with mixed quantum-classical media

This paper offers an expedient, efficient, and unique treatment of multimode quantum subsystems (polyatomic molecules) interacting with a classical environment in which the time evolution of the coupling term is governed by the algebraic rules of statistical mechanics in mixed quantum-classical systems developed by Kapral and Nielsen [S. Nielsen, R. Kapral, and G. Ciccotti, J. Chem. Phys., 2001, 115, 5805]. This unique time evolution of the coupling term is neither quantal nor classical but rather something different that relies heavily on Wigner transform, thereby leading to non-Newtonian mechanics. As such, an argument is presented that the approach provided herein for treating polyatomic molecular systems in a mixed quantum-classical environment is new and different as opposed to the many other schemes of semiclassical dynamics that are normally employed to study such systems. The merits of expediency and efficiency of the herein mixed quantum-classical dynamics calculations emanate from avoiding using integrals for time evolutions, and, instead, employing matrix mechanics whereby LU decomposition and singular value decomposition (SVD) numerical techniques are utilized for diagonalization. An electronic 2-level subsystem interacting with a classical bath through the spin-boson model to render accurate pure electronic dephasing in multimode molecular systems by eliminating the unphysical asymmetry in the line shape of the zero-phonon line (ZPL) exhibited by other models is exploited. This work has a superior advantage over the single-mode spin-boson model, published previously, whereby a multitude of types of vibrational modes (slow, fast, or both) of the quantum subsystem may readily be handled using different spectral densities. The spin-boson model used here is a composite system made up of a quantum subsystem, i.e., a subsystem bilinearly coupled to a multidimensional harmonic oscillator (representing the intermediate quantum vibrational modes between the electronic subsystem and the bath), interacting with a classical bath, where the coupling term is governed by the mixed quantum-classical Liouville equation. A multidimensional coherent-state approach is employed to deal with the time evolution of the quantum subsystem. A closed-form expression of linear and nonlinear optical electronic transition dipole moment time correlation functions in mixed quantum-classical dissipative media is derived. Pure electronic dephasing is probed using the aforementioned approach. Linear absorption spectra and 4-wave mixing signals (e.g., photon echo and pump-probe) are calculated showing a reasonable thermal broadening, temporal decay, and accurate pure dephasing.

Excited state distribution function for probing Herzberg-Teller vibronic coupling using linear optical response theory: Application to glassy pheophytin a

The goal of the present work is to develop an excited-state distribution function that can be used to calculate electronic transition dipole moment time correlation functions at a considerably low computational cost. An additional merit of the distribution function is its capability to probe the Hertzberg-Teller vibronic coupling effect in terms of the previously reported Condon correlation functions in the literature without having to start from the equilibrium density operator to probe spectral non-Condon effects, thereby exploring Hertzberg-Teller vibronic coupling by building on the Condon regime. It is easily extendable to anharmonic systems. Model calculations are reported to show the high degree of accuracy and computational efficiency of the presented approach. Application to a photosynthetic system such as pheophytin a in triethylamine is provided.

Caldeira-Leggett model vs ab initio potential: A vibrational spectroscopy test of water solvation

We present a semiclassically approximate quantum treatment of solvation with the purpose of investigating the accuracy of the Caldeira-Leggett model. We do that by simulating the vibrational features of water solvation by means of two different approaches. One is entirely based on the adoption of an accurate ab initio potential to describe water clusters of increasing dimensionality. The other one consists of a model made of a central water molecule coupled to a high-dimensional Caldeira-Leggett harmonic bath. We demonstrate the role of quantum effects in the detection of water solvation and show that the computationally cheap approach based on the Caldeira-Leggett bath is only partially effective. The main conclusion of the study is that quantum methods associated with high-level potential energy surfaces are necessary to correctly study solvation features, while simplified models, even if attractive owing to their reduced computational cost, can provide some useful insights but are not able to come up with a comprehensive description of the solvation phenomenon.

Quantum Trajectories for the Dynamics in the Exact Factorization Framework: A Proof-of-Principle Test

In the framework of the exact factorization of the time-dependent electron-nuclear wave function, we investigate the possibility of solving the nuclear time-dependent Schrödinger equation based on trajectories. The nuclear equation is separated in a Hamilton-Jacobi equation for the phase of the wave function, and a continuity equation for its (squared) modulus. For illustrative adiabatic and nonadiabatic one-dimensional models, we implement a procedure to follow the evolution of the nuclear density along the characteristics of the Hamilton-Jacobi equation. Those characteristics are referred to as quantum trajectories, since they are generated via ordinary differential equations similar to Hamilton's equations, but including the so-called quantum potential, and they can be used to reconstruct exactly the quantum-mechanical nuclear wave function, provided infinite initial conditions are propagated in time.

On the Husimi Version of the Classical Limit of Quantum Correlation Functions

We extend the Husimi (coherent state) based version of linearized semiclassical theories for the calculation of correlation functions to the case of survival probabilities. This is a case that could be dealt with before only by use of the Wigner version of linearized semiclassical theory. Numerical comparisons of the Husimi and the Wigner case with full quantum results as well as with full semiclassical ones will be given for the revival dynamics in a Morse oscillator with and without coupling to an additional harmonic degree of freedom.

Semiclassical dynamics in the mixed quantum-classical limit

The semiclassical double Herman-Kluk initial value representation is an accurate approach to computing quantum real time correlation functions, but its applications are limited by the need to evaluate an oscillatory integral. In previous work, we have shown that this "sign problem" can be mitigated using the modified Filinov filtration technique to control the extent to which individual modes of the system contribute to the overall phase of the integrand. Here, we follow this idea to a logical conclusion: we analytically derive a general expression for the mixed quantum-classical limit of the semiclassical correlation function-analytical mixed quantum-classical-initial value representation (AMQC-IVR), where the phase contributions from the "classical" modes of the system are filtered while the "quantum" modes are treated in the full semiclassical limit. We numerically demonstrate the accuracy and efficiency of the AMQC-IVR formulation in calculations of quantum correlation functions and reaction rates using three model systems with varied coupling strengths between the classical and quantum subsystems. We also introduce a separable prefactor approximation that further reduces computational cost but is only accurate in the limit of weak coupling between the quantum and classical subsystems.

The importance of the pre-exponential factor in semiclassical molecular dynamics

This paper deals with the critical issue of approximating the pre-exponential factor in semiclassical molecular dynamics. The pre-exponential factor is important because it accounts for the quantum contribution to the semiclassical propagator of the classical Feynman path fluctuations. Pre-exponential factor approximations are necessary when chaotic or complex systems are simulated. We introduced pre-exponential factor approximations based either on analytical considerations or numerical regularization. The approximations are tested for power spectrum calculations of more and more chaotic model systems and on several molecules, for which exact quantum mechanical values are available. The results show that the pre-exponential factor approximations introduced are accurate enough to be safely employed for semiclassical simulations of complex systems.

Improving the efficiency of Monte Carlo surface hopping calculations

A surface hopping method with a Monte Carlo procedure for deciding whether to hop at each step along the classical trajectories used in the semiclassical calculation is discussed. It is shown for a simple one-dimensional model problem that the numerical efficiency of the method can be improved by averaging over several copies of the sections of each trajectory that span the interaction regions. The use of Sobol sequences in the selection of the initial momentum for the trajectories is also explored. It is found that accurate results can be obtained with relatively small trajectory samples.

Effective Markovian description of decoherence in bound systems

Effective descriptions accounting for the evolution of quantum systems that are acted on by a bath are desirable. As the number of bath degrees of freedom increases and full quantum simulations turn out computationally prohibitive, simpler models become essential to understand and gain an insight into the main physical mechanisms involved in the system dynamics. In this regard, vibrational decoherence of an I2 diatomics is tackled here within the framework of Markovian quantum state diffusion. The I2 dynamics are analyzed in terms of an effective decoherence rate, Λ, and the specific choice of the initial state, in particular, Gaussian wave packets and two-state superpositions. It is found that, for Markovian baths, the relevant quantity regarding decoherence is the product of friction (η) and temperature (T); there is no distinction between varying one or the other. It is also observed that decoherence becomes faster as the energy levels involved in the system state correspond to higher eigenvalues. This effect is due to a population redistribution during the dynamical process and an eventual irreversible loss of the initial coherence. These results have been compared with those available in the literature from more detailed semiclassical IVR simulations, finding a good agreement.

Semiclassical initial-value representation of the transfer operator

The ability of semiclassical initial-value representation (IVR) methods to determine approximate energy levels for bound systems is limited due to problems associated with long classical trajectories. These difficulties become especially severe for large or classically chaotic systems. This work attempts to overcome such problems by developing an IVR expression that is classically equivalent to Bogomolny's formula for the transfer matrix [E. B. Bogomolny, Nonlinearity 5, 805 (1992); Chaos 2, 5 (1992)] and can be used to determine semiclassical energy levels. The method is adapted to levels associated with states of desired symmetries and applied to two two-dimensional quartic oscillator systems, one integrable and one mostly chaotic. For both cases, the technique is found to resolve all energy levels in the ranges investigated. The IVR method does not require a search for special trajectories obeying boundary conditions on the Poincaré surface of section and leads to more rapid convergence of Monte Carlo phase space integrations than a previously developed IVR technique. It is found that semiclassical energies can be extracted from the eigenvalues of transfer matrices of dimension close to the theoretical minimum determined by Bogomolny's theory. The results support the assertion that the present IVR theory provides a different semiclassical approximation to the transfer matrix than that of Bogomolny for ℏ≠0. For the chaotic system investigated the IVR energies are found to be generally more accurate than those predicted by Bogomolny's theory.

Electronically nonadiabatic dynamics in complex molecular systems: an efficient and accurate semiclassical solution

Chemical reaction dynamics is always a central theme in chemistry research. In many important chemical processes, reaction dynamics is electronically nonadiabatic, i.e., dynamics involves coupled multiple electronic states. We demonstrate in this paper that a semiclassical (SC) treatment based on an initial value representation methodology and a classical mapping formalism for the electronic degrees of freedom is now able to provide a rigorous and practical solution to electronically nonadiabatic dynamics in complex molecular systems. The key component of this treatment is to incorporate a correlated importance sampling protocol in nonadiabatic SC calculations, which results in a speedup factor of 100 or more in comparison with that using the standard sampling approach. This is illustrated by application to a two-state model coupled with up to 10 nuclear bath modes for a benchmark nonadiabatic excitation energy transfer problem. This work provides great opportunities for the effectively theoretical investigations on reaction mechanisms in complex molecular systems, in which electronically nonadiabatic dynamics plays an importance role.

Communication: importance sampling including path correlation in semiclassical initial value representation calculations for time correlation functions

Full semiclassical (SC) initial value representation (IVR) for time correlation functions involves a double phase space average over a set of two phase points, each of which evolves along a classical path. Conventionally, the two initial phase points are sampled independently for all degrees of freedom (DOF) in the Monte Carlo procedure. Here, we present an efficient importance sampling scheme by including the path correlation between the two initial phase points for the bath DOF, which greatly improves the performance of the SC-IVR calculations for large molecular systems. Satisfactory convergence in the study of quantum coherence in vibrational relaxation has been achieved for a benchmark system-bath model with up to 21 DOF.

Insights in quantum dynamical effects in the infrared spectroscopy of liquid water from a semiclassical study with an ab initio-based flexible and polarizable force field

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.

FORWARD-BACKWARD SEMICLASSICAL SIMULATION OF DYNAMICAL PROCESSES IN LIQUIDS

Forward-backward semiclassical dynamics (FBSD) provides a practical methodology for including quantum mechanical effects in classical trajectory simulations of polyatomic systems. FBSD expressions for time-dependent expectation values or correlation functions take the form of phase space integrals with respect to trajectory initial conditions, weighted by the coherent state transform of a corrected density operator. Quantization through a discretized path integral representation of the Boltzmann operator ensures a proper treatment of zero point energy effects and of imaginary components in finite-temperature correlation functions, and extension to systems obeying Bose statistics is possible. Accelerated convergence is achieved via Monte Carlo or molecular dynamics sampling techniques and through the construction of improved imaginary time propagators. The accuracy of the methodology is demonstrated on several model systems, including models of Bose and Fermi particles. Applications to liquid argon, neon and para-hydrogen are presented.

A NEW SEMICLASSICAL DYNAMICS FROM THE INTERACTION REPRESENTATION

This paper develops a new semiclassical mechanics from an exact quantum prescription. In this formulation, a zeroth order propagation, rather than Hamiltonian, is specified. The exact, full evolution operator is then given from a specific interaction representation of the evolved "perturbation" Hamiltonian. We then investigate a variety of approximate, semiclassical, and mixed Quantum/Classical methods, along with exact methodologies to evaluate this time dependent interaction Hamiltonian. The approximate full evolution operator can be described in a variety of ways including an iterated Lippman-Schwinger like equation, and an expansion of the perturbation propagator generated from the time evolved Hamiltonian.