Reduced equations of motion for semiclassical dynamics in phase space

A Gaussian Wave Packet Propagation Approach to Vibrationally Resolved Optical Spectra at Non-Zero Temperatures

An effective time dependent approach based on a method that is similar to the Gaussian wave packet propagation (GWP) technique of Heller is developed for the computation of vibrationally resolved electronic spectra at finite temperatures in the harmonic, Franck-Condon/Hertzberg-Teller approximations. Since the vibrational thermal density matrix of the ground electronic surface and the time evolution operator on that surface commute, it is possible to write the spectrum generating correlation function as a trace of the time evolved doorway state. In the stated approximations, the doorway state is a superposition of the harmonic oscillator zero and one quantum eigenfunctions and thus can be propagated by the GWP. The algorithm has an O(N(3)) dependence on the number of vibrational modes. An application to pyrene absorption spectrum at two temperatures is presented as a proof of the concept.

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.

Anharmonic electron-phonon coupling in condensed media: 1. Formalism

Three different schemes for calculating anharmonic line shape functions are reported and discussed for the first time in this article using eigenstate representation. First, the linear dipole-moment time correlation function (DMTCF), homogeneous (single-site) absorption line shape function, and the respective Franck-Condon factors (FCF) are derived and explored as a molecule makes a transition from a harmonic to an anharmonic (Morse potential) electronic state. Second, the linear DMTCF, homogeneous absorption line shape function, and FCFs are also derived as a molecule makes a transition from one anharmonic to another linearly displaced anharmonic state; FCFs in this case are reported in an exact closed-form expressed in terms of Appell's hypergeometric function. Third, same as the latter set of results are reported but with both linearly displaced and distorted shape of the upper Morse potential. FCFs of the zero-phonon line in all three cases are reported. The first case is rather mathematically complex as a result of taking the overlap integral of the Morse oscillator eigenfunctions, whose spatial decay is a simple exponential, with those of harmonic oscillator, whose decay is a Gaussian. This form of a functional disparity gives rise to some challenges. Model calculations are presented and discussed.

Second-order quantized Hamilton dynamics coupled to classical heat bath

Starting with a quantum Langevin equation describing in the Heisenberg representation a quantum system coupled to a quantum bath, the Markov approximation and, further, the closure approximation are applied to derive a semiclassical Langevin equation for the second-order quantized Hamilton dynamics (QHD) coupled to a classical bath. The expectation values of the system operators are decomposed into products of the first and second moments of the position and momentum operators that incorporate zero-point energy and moderate tunneling effects. The random force and friction as well as the system-bath coupling are decomposed to the lowest classical level. The resulting Langevin equation describing QHD-2 coupled to classical bath is analyzed and applied to free particle, harmonic oscillator, and the Morse potential representing the OH stretch of the SPC-flexible water model.

A canonical averaging in the second-order quantized Hamilton dynamics

Quantized Hamilton dynamics (QHD) is a simple and elegant extension of classical Hamilton dynamics that accurately includes zero-point energy, tunneling, dephasing, and other quantum effects. Formulated as a hierarchy of approximations to exact quantum dynamics in the Heisenberg formulation, QHD has been used to study evolution of observables subject to a single initial condition. In present, we develop a practical solution for generating canonical ensembles in the second-order QHD for position and momentum operators, which can be mapped onto classical phase space in doubled dimensionality and which in certain limits is equivalent to thawed Gaussian. We define a thermal distribution in the space of the QHD-2 variables and show that the standard beta=1/kT relationship becomes beta'=2/kT in the high temperature limit due to an overcounting of states in the extended phase space, and a more complicated function at low temperatures. The QHD thermal distribution is used to compute total energy, kinetic energy, heat capacity, and other canonical averages for a series of quartic potentials, showing good agreement with the quantum results.