Finite-temperature vibronic spectra from the split-operator coherence thermofield dynamics
arXiv:2311.10004 · doi:10.1063/5.0187823
Abstract
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.
arXiv comment5 pages, 4 figures
References in corpus (8)
- On-the-fly ab initio semiclassical evaluation of absorption spectra of polyatomic molecules beyond the Condon approximation
- On-the-fly ab initio semiclassical evaluation of vibronic spectra at finite temperature
- Semiclassical Approach to Photophysics Beyond Kasha's Rule and Vibronic Spectroscopy Beyond the Condon Approximation. The Case of Azulene
- Finite-temperature, anharmonicity, and Duschinsky effects on the two-dimensional electronic spectra from ab initio thermo-field Gaussian wavepacket dynamics
- Efficient geometric integrators for nonadiabatic quantum dynamics. II. The diabatic representation
- Efficient geometric integrators for nonadiabatic quantum dynamics. I. The adiabatic representation
- Internal Conversion Rates from the Extended Thawed Gaussian Approximation: Theory and Validation
- A time-reversible integrator for the time-dependent Schrödinger equation on an adaptive grid