Vibronic spectra at nonzero temperatures from Herman-Kluk coherence thermofield dynamics

Abstract: We combine the semiclassical Herman--Kluk approximation with the coherence thermofield dynamics in order to evaluate vibrationally resolved electronic spectra at nonzero temperatures. In coherence thermofield dynamics, the dipole time correlation function is rewritten exactly as a wavepacket autocorrelation function, and the corresponding wavepacket is a solution to a zero-temperature time-dependent Schr\"odinger equation on an augmented configuration space of doubled dimension. We derive the Herman--Kluk representation for the thermofield wavepacket autocorrelation function and demonstrate how it can be computed from individual trajectories. To analyze this method, we compare spectra of Morse potentials of increasing anharmonicity evaluated at various temperatures with a numerically exact approach, with the Herman--Kluk coherence thermofield dynamics, and with the single-trajectory thawed Gaussian coherence thermofield dynamics. At low anharmonicity, both approximate methods yield accurate spectra. However, in a Morse potential with higher anharmonicity, the thawed Gaussian thermofield dynamics, based on the local harmonic approximation, fails to capture emerging hot bands, whereas the Herman--Kluk thermofield approach successfully reproduces them.