Direct evaluation of boson dynamics via finite-temperature time-dependent variation with multiple Davydov states

Non-Markovian Vibrational Relaxation Dynamics at Surfaces

Vibrational dynamics of adsorbates near surfaces plays both an important role for applied surface science and as model lab for studying fundamental problems of open quantum systems. We employ a previously developed model for the relaxation of a D-Si-Si bending mode at a D:Si(100)-(2$\times$1) surface, induced by a ``bath' of more than $2000$ phonon modes [U. Lorenz, P. Saalfrank, Chem. Phys. {\bf 482}, 69 (2017)], to extend previous work along various directions. First, we use a Hierarchical Effective Mode (HEM) model [E.W. Fischer, F. Bouakline, M. Werther, P. Saalfrank, J. Chem. Phys. {\bf 153}, 064704 (2020)] to study relaxation of higher excited vibrational states than hitherto done, by solving a high-dimensional system-bath time-dependent Schr\"odinger equation (TDSE). In the HEM approach, (many) real bath modes are replaced by (much less) effective bath modes. Accordingly, we are able to examine scaling laws for vibrational relaxation lifetimes for a realistic surface science problem. Second, we compare the performance of the multilayer multiconfigurational time-dependent Hartree (ML-MCTDH) approach with the recently developed coherent-state based multi-Davydov D2 {\it ansatz} [N. Zhou, Z. Huang, J. Zhu, V. Chernyak, Y. Zhao, {J. Chem. Phys.} {\bf 143}, 014113 (2015)]. Both approaches work well, with some computational advantages for the latter in the presented context. Third, we apply open-system density matrix theory in comparison with basically ``exact' solutions of the multi-mode TDSEs. Specifically, we use an open-system Liouville-von Neumann (LvN) equation treating vibration-phonon coupling as Markovian dissipation in Lindblad form to quantify effects beyond the Born-Markov approximation.

Lamb Shift and the Vacuum Rabi Splitting in a Strongly Dissipative Environment

We study the vacuum Rabi splitting of a qubit ultrastrongly coupled to a high-Q cavity mode and a radiation reservoir. Three methods are employed: a numerically exact variational approach with a multiple Davydov ansatz, the rotating-wave approximation (RWA), and the transformed RWA. Agreement between the variational results and the transformed RWA results is found in the regime of validity of the latter, where the RWA breaks down completely. We illustrate that the Lamb shift plays an essential role in modifying the vacuum Rabi splitting in the ultrastrong coupling regime, leading to off-resonant qubit-cavity coupling even though the cavity frequency equals the bare transition frequency of the qubit. Specifically, the emission spectrum exhibits one broad low-frequency peak and one narrow high-frequency peak in the presence of relatively weak cavity-qubit coupling. As the cavity-qubit coupling increases, the low-frequency peak narrows while the high-frequency peak broadens until they have similar widths.

Schrödinger-Cat States in Landau-Zener-Stückelberg-Majorana Interferometry: A Multiple Davydov Ansatz Approach

Employing the time-dependent variational principle combined with the multiple Davydov D₂ Ansatz, we investigate Landau-Zener (LZ) transitions in a qubit coupled to a photon mode with various initial photon states at zero temperature. Thanks to the multiple Davydov trial states, exact photonic dynamics taking place in the course of the LZ transition is also studied efficiently. With the qubit driven by a linear external field and the photon mode initialized with Schrödinger-cat states, asymptotic behavior of the transition probability beyond the rotating-wave approximation is uncovered for a variety of initial states. Using a sinusoidal external driving field, we also explore the photon-assisted dynamics of Landau-Zener-Stückelberg-Majorana interferometry. Transition pathways involving multiple energy levels are unveiled by analyzing the photon dynamics.

Efficient simulation of time- and frequency-resolved four-wave-mixing signals with a multiconfigurational Ehrenfest approach

We have extended the multiconfigurational Ehrenfest approach to the simulation of four-wave-mixing signals of systems involving multiple electronic and vibrational degrees of freedom. As an illustration, we calculate signals of three widely used spectroscopic techniques, time- and frequency-resolved fluorescence spectroscopy, transient absorption spectroscopy, and two-dimensional (2D) electronic spectroscopy, for a two-electronic-state, twenty-four vibrational-mode conical intersection model. It has been shown that all these three spectroscopic signals characterize fast population transfer from the higher excited electronic state to the lower excited electronic state. While the time- and frequency-resolved spectrum maps the wave packet propagation exclusively on the electronically excited states, the transient absorption and 2D electronic spectra reflect the wave packet dynamics on both electronically excited states and the electronic ground state. Combining trajectory-guided Gaussian basis functions and the nonlinear response function formalism, the present approach provides a promising general technique for the applications of various Gaussian basis methods to the calculations of four-wave-mixing spectra of polyatomic molecules.