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

Measuring vibronic spectra 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 work, we express the imaginary-time shifted time correlation function and, thus, the vibronic spectrum in terms of the so-called Matsubara dynamics, which combines quantum statistics and classical-like dynamics. By applying the Matsubara approximation in the adiabatic limit, we derive a formal generalization of the existing Matsubara dynamics formalism to multiple potential energy surfaces (PESs), which, however, does not feature all the defining properties of its single-PES counterpart though suffering equally from the sign problem. The mathematical analysis for two shifted harmonic oscillators suggests a new modified method to practically simulate the standard correlation function via Matsubara-like dynamics. Importantly, this modified method samples the thermal Wigner function without suffering from the sign problem and yields an accurate approximation to the vibronic absorption spectrum, not only for the harmonic system but also for the anharmonic one.

Quasi-classical approaches to vibronic spectra revisited.

The framework to approach quasi-classical dynamics in the electronic ground state is well established and is based on the Kubo-transformed time correlation function (TCF), being the most classical-like quantum TCF. Here we discuss whether the choice of the Kubo-transformed TCF as a starting point for simulating vibronic spectra is as unambiguous as it is for vibrational ones. Employing imaginary-time path integral techniques in combination with the interaction representation allowed us to formulate a method for simulating vibronic spectra in the adiabatic regime that takes nuclear quantum effects and dynamics on multiple potential energy surfaces into account. Further, a generalized quantum TCF is proposed that contains many well-established TCFs, including the Kubo one, as particular cases. Importantly, it also provides a framework to construct new quantum TCFs. Applying the developed methodology to the generalized TCF leads to a plethora of simulation protocols, which are based on the well-known TCFs as well as on new ones. Their performance is investigated on 1D anharmonic model systems at finite temperatures. It is shown that the protocols based on the new TCFs may lead to superior results with respect to those based on the common ones. The strategies to find the optimal approach are discussed.

A time-correlation function approach to nuclear dynamical effects in X-ray spectroscopy.

Modern X-ray spectroscopy has proven itself as a robust tool for probing the electronic structure of atoms in complex environments. Despite working on energy scales that are much larger than those corresponding to nuclear motions, taking nuclear dynamics and the associated nuclear correlations into account may be of importance for X-ray spectroscopy. Recently, we have developed an efficient protocol to account for nuclear dynamics in X-ray absorption and resonant inelastic X-ray scattering spectra [Karsten et al., J. Phys. Chem. Lett. 8, 992 (2017)], based on ground state molecular dynamics accompanied with state-of-the-art calculations of electronic excitation energies and transition dipoles. Here, we present an alternative derivation of the formalism and elaborate on the developed simulation protocol using gas phase and bulk water as examples. The specific spectroscopic features stemming from the nuclear motions are analyzed and traced down to the dynamics of electronic energy gaps and transition dipole correlation functions. The observed tendencies are explained on the basis of a simple harmonic model, and the involved approximations are discussed. The method represents a step forward over the conventional approaches that treat the system in full complexity and provides a reasonable starting point for further improvements.

Nuclear Dynamical Correlation Effects in X-ray Spectroscopy from a Theoretical Time-Domain Perspective.

To date X-ray spectroscopy has become a routine tool that can reveal highly local and element-specific information on the electronic structure of atoms in complex environments. Here, we focus on nuclear dynamical correlation effects in X-ray spectra and develop a rigorous time-correlation function method employing ground state classical molecular dynamics simulations. The importance of nuclear correlation phenomena is demonstrated by comparison against the results from the conventional sampling approach performed on the same data set for gas phase water. In contrast to the first-order absorption, second-order resonant inelastic scattering spectra exhibit pronounced fingerprints of nuclear motions. The developed methodology is not biased to a particular electronic structure method and, owing to its generality, can be applied to, e.g., X-ray photoelectron and Auger spectroscopies.

Parametrizing linear generalized Langevin dynamics from explicit molecular dynamics simulations.

Fundamental understanding of complex dynamics in many-particle systems on the atomistic level is of utmost importance. Often the systems of interest are of macroscopic size but can be partitioned into a few important degrees of freedom which are treated most accurately and others which constitute a thermal bath. Particular attention in this respect attracts the linear generalized Langevin equation, which can be rigorously derived by means of a linear projection technique. Within this framework, a complicated interaction with the bath can be reduced to a single memory kernel. This memory kernel in turn is parametrized for a particular system studied, usually by means of time-domain methods based on explicit molecular dynamics data. Here, we discuss that this task is more naturally achieved in frequency domain and develop a Fourier-based parametrization method that outperforms its time-domain analogues. Very surprisingly, the widely used rigid bond method turns out to be inappropriate in general. Importantly, we show that the rigid bond approach leads to a systematic overestimation of relaxation times, unless the system under study consists of a harmonic bath bi-linearly coupled to the relevant degrees of freedom.