Quantum effects in intermediate-temperature dipole-dipole correlation-functions in the presence of an environment

Finite-temperature vibronic spectra from the split-operator coherence thermofield dynamics

We present a numerically exact approach for evaluating vibrationally resolved electronic spectra at finite temperatures using the coherence thermofield dynamics. In this method, which avoids implementing an algorithm for solving the von Neumann equation for coherence, the thermal vibrational ensemble is first mapped to a pure-state wavepacket in an augmented space, and this wavepacket is then propagated by solving the standard, zero-temperature Schrödinger equation with the split-operator Fourier method. We show that the finite-temperature spectra obtained with the coherence thermofield dynamics in a Morse potential agree exactly with those computed by Boltzmann-averaging the spectra of individual vibrational levels. Because the split-operator thermofield dynamics on a full tensor-product grid is restricted to low-dimensional systems, we briefly discuss how the accessible dimensionality can be increased by various techniques developed for the zero-temperature split-operator Fourier method.

Numerically "exact" approach to open quantum dynamics: The hierarchical equations of motion (HEOM)

An open quantum system refers to a system that is further coupled to a bath system consisting of surrounding radiation fields, atoms, molecules, or proteins. The bath system is typically modeled by an infinite number of harmonic oscillators. This system-bath model can describe the time-irreversible dynamics through which the system evolves toward a thermal equilibrium state at finite temperature. In nuclear magnetic resonance and atomic spectroscopy, dynamics can be studied easily by using simple quantum master equations under the assumption that the system-bath interaction is weak (perturbative approximation) and the bath fluctuations are very fast (Markovian approximation). However, such approximations cannot be applied in chemical physics and biochemical physics problems, where environmental materials are complex and strongly coupled with environments. The hierarchical equations of motion (HEOM) can describe the numerically "exact" dynamics of a reduced system under nonperturbative and non-Markovian system-bath interactions, which has been verified on the basis of exact analytical solutions (non-Markovian tests) with any desired numerical accuracy. The HEOM theory has been used to treat systems of practical interest, in particular, to account for various linear and nonlinear spectra in molecular and solid state materials, to evaluate charge and exciton transfer rates in biological systems, to simulate resonant tunneling and quantum ratchet processes in nanodevices, and to explore quantum entanglement states in quantum information theories. This article presents an overview of the HEOM theory, focusing on its theoretical background and applications, to help further the development of the study of open quantum dynamics.

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.

Thermal weights for semiclassical vibrational response functions

Semiclassical approximations to response functions can allow the calculation of linear and nonlinear spectroscopic observables from classical dynamics. Evaluating a canonical response function requires the related tasks of determining thermal weights for initial states and computing the dynamics of these states. A class of approximations for vibrational response functions employs classical trajectories at quantized values of action variables and represents the effects of the radiation-matter interaction by discontinuous transitions. Here, we evaluate choices for a thermal weight function which are consistent with this dynamical approximation. Weight functions associated with different semiclassical approximations are compared, and two forms are constructed which yield the correct linear response function for a harmonic potential at any temperature and are also correct for anharmonic potentials in the classical mechanical limit of high temperature. Approximations to the vibrational linear response function with quantized classical trajectories and proposed thermal weight functions are assessed for ensembles of one-dimensional anharmonic oscillators. This approach is shown to perform well for an anharmonic potential that is not locally harmonic over a temperature range encompassing the quantum limit of a two-level system and the limit of classical dynamics.