Simulating vibronic spectra via Matsubara-like dynamics: Coping with the sign problem

Abstract

Measuring the vibronic spectrum probes dynamical processes in molecular systems.When interpreted via suitable theoretical tools, the experimental data provides comprehensive information about the system in question. For complex many-body problems, such an approach usually requires the formulation of proper classical-like approximations, which is particularly challenging if multiple electronic states are involved. In this manuscript, we express the imaginary-time shifted time correlation function (TCF) and, thus, the vibronic spectrum in terms of the so-called Matsubara dynamics, which combines quantum statistics and classical-like dynamics. In contrast to the existing literature, we invoke a local harmonic approximation to the potential allowing an analytical evaluation of integrals. By subsequently applying the Matsubara approximation, we derive a generalization of the existing Matsubara method to multiple potential energy surfaces (PESs), which, however, suffers from the sign problem as its single-PES counterpart does. The mathematical analysis for two shifted harmonic oscillators suggests a new modified method to simulate the standard correlation function via classical-like dynamics. Importantly, this modified method samples the thermal Wigner function without suffering from the sign problem and it yields an accurate approximation to the vibronic absorption spectrum, not only for the harmonic system, but also for an anharmonic one.

Polarizable embedding QM/MM: the future gold standard for complex (bio)systems?

“…The approach we review here is based on expressing the vibronic coupling in terms of a "spectral density" function (SD), which describes the frequency-dependent linear coupling between the electronic excitation and the nuclear degrees of freedom. [153][154][155][156] The electronic lineshapes can be expressed in a timedomain approach within the second-order cumulant expansion formalism. 153 From a Born-Oppenheimer ground-state trajectory, the SD can be obtained from the Fourier transform of the ACF of the excitation energy fluctuations, i.e., C cl (t) = U(t)U(0) , where U(t) =…”

On the “Matsubara heating” of overtone intensities and Fermi splittings

“…The HD approximation is a modification of an approximation introduced in ref. 69. It is based on the observation that, in a harmonic system, the terms…”

Theoretical X‐ray spectroscopy of transition metal compounds

“…45 For instance for XAS, these are two-point and for RIXS four-point dipole correlation functions. [62][63][64][65] Beyond this perturbational regime, that is, for strong-field excitation, or in cases of (shaped) ultrashort laser pulses, one has to resort to explicit solution of the time-dependent Schrödinger equation including the field-matter interaction. Compared to the QC methods, real-time (RT) propagation approaches have received less attention so far.…”