By employing the numerically accurate multiple Davydov Ansatz (mDA) formalism in combination with the thermo-field dynamics (TFD) representation of quantum mechanics, we systematically explore the influence of three parameters-temperature, photonic-mode detuning, and qubit-phonon coupling-on population dynamics and absorption spectra of the Holstein-Tavis-Cummings (HTC) model. It is found that elevated qubit-phonon couplings and/or temperatures have a similar impact on all dynamic observables: they suppress the amplitudes of Rabi oscillations in photonic populations as well as broaden the peaks and decrease their intensities in the absorption spectra. Our results unequivocally demonstrate that the HTC dynamics is very sensitive to the concerted variation of the three aforementioned parameters, and this finding can be used for fine-tuning polaritonic transport. The developed mDA-TFD methodology can be efficiently applied for modeling, predicting, optimizing, and comprehensively understanding dynamic and spectroscopic responses of actual molecular systems in microcavities.
By employing the numerically accurate multiple Davydov Ansatz (mDA) formalism in combination with the thermo-field dynamics (TFD) representation of quantum mechanics, we systematically explore the influence of three parameters-temperature, photonic-mode detuning, and qubit-phonon coupling-on population dynamics and absorption spectra of the Holstein-Tavis-Cummings (HTC) model. It is found that elevated qubit-phonon couplings and/or temperatures have a similar impact on all dynamic observables: they suppress the amplitudes of Rabi oscillations in photonic populations as well as broaden the peaks and decrease their intensities in the absorption spectra. Our results unequivocally demonstrate that the HTC dynamics is very sensitive to the concerted variation of the three aforementioned parameters, and this finding can be used for fine-tuning polaritonic transport. The developed mDA-TFD methodology can be efficiently applied for modeling, predicting, optimizing, and comprehensively understanding dynamic and spectroscopic responses of actual molecular systems in microcavities.
Davydov trial states, in combination with the Green's function method, we study the dynamics of the Tavis-Cummings model and the Holstein-Tavis-Cummings model in the presence of diagonal disorder and cavity-qubit coupling disorder. For the Tavis-Cummings model, time evolution of the photon population, the optical absorption spectra, and the hetero-entanglement between the qubits and the cavity mode are calculated by using the Green's function method to corroborate numerically exact results of Davydov's Ansätze. For the Holstein-Tavis-Cummings model, only the latter is utilized to simulate effects of disorder on the photon population dynamics and the absorption spectra. We have demonstrated that the complementary employment of analytical and numerical methods permits uncovering a fairly comprehensive picture of a variety of complex behaviors in disordered multidimensional polaritonic cavity quantum electrodynamics systems.
We examine the applicability of the numerically accurate method of time dependent variation with multiple Davydov Ansätze (mDA) to non-Hermitian systems. As illustrative examples, three systems of interest have been studied, a non-Hermitian system of dissipative Landau-Zener transitions, a non-Hermitian multimode Jaynes-Cummings model, and a dissipative Holstein-Tavis-Cummings model, all of which are shown to be effectively described by the mDA method. Our findings highlight the versatility of the mDA as a powerful numerical tool for investigating complex many-body non-Hermitian systems, which can be extended to explore diverse phenomena such as skin effects, excited-state dynamics, and spectral topology in the non-Hermitian field.
We perform benchmark simulations using the time-dependent variational approach with the multiple Davydov Ansatz (mDA) to study real-time nonequilibrium dynamics in a single qubit model coupled to two thermal baths with distinct temperatures. A broad region of the parameter space has been investigated, accompanied by a detailed analysis of the convergence behavior of the mDA method. In addition, we have compared our mD2 results to those from two widely adopted, numerically "exact" techniques: the methods of hierarchical equations of motion (HEOM) and the quasi-adiabatic path integral (QUAPI). It is found that the mDA approach in combination with thermal field dynamics yields numerically accurate, convergent results in nearly all regions of the parameter space examined, including those that pose serious challenges for QUAPI and HEOM. Our results reveal that mDA offers a highly adaptable framework capable of capturing long-time dynamics, even in challenging regimes where other methods face limitations. These findings underscore the potential of mDA as a versatile tool for exploring quantum thermodynamics, energy transfer processes, and non-equilibrium quantum systems.
This Perspective presents an overview of the development of the hierarchy of Davydov's Ansätze and a few of their applications in many-body problems in computational chemical physics. Davydov's solitons originated in the investigation of vibrational energy transport in proteins in the 1970s. Momentum-space projection of these solitary waves turned up to be accurate variational ground-state wave functions for the extended Holstein molecular crystal model, lending unambiguous evidence to the absence of formal quantum phase transitions in Holstein systems. The multiple Davydov Ansätze have been proposed, with increasing Ansatz multiplicity, as incremental improvements of their single-Ansatz parents. For a given Hamiltonian, the time-dependent variational formalism is utilized to extract accurate dynamic and spectroscopic properties using Davydov's Ansätze as its trial states. A quantity proven to disappear for large multiplicities, the Ansatz relative deviation is introduced to quantify how closely the Schrödinger equation is obeyed. Three finite-temperature extensions to the time-dependent variation scheme are elaborated, i.e., the Monte Carlo importance sampling, the method of thermofield dynamics, and the method of displaced number states. To demonstrate the versatility of the methodology, this Perspective provides applications of Davydov's Ansätze to the generalized Holstein Hamiltonian, variants of the spin-boson model, and systems of cavity-assisted singlet fission, where accurate dynamic and spectroscopic properties of the many-body systems are given by the Davydov trial states.
The field of Polaritonic Chemistry studies molecules strongly coupled to the light field of optical cavities. The resulting modification of the potential energy landscape has been demonstrated to enable control over the chemical and physical properties of materials and molecular ensembles. This thesis extends the Tavis–Cummings model to describe molecules under strong light–matter coupling, incorporating static dipole moments and dipole self-energy contributions derived from the Pauli–Fierz Hamiltonian. These additional terms are shown to be essential for reproducing cavity-modified molecular dynamics and spectra. Simulations on the nuclear wavefunction propagation of MgH+ molecules demonstrate the impact of these terms, while an effective two-level system model accurately captures the behavior of large ensembles. Beyond single excitations, the thesis also explores the higher excitation manifolds of the Tavis-Cummings model, which reveals that allowing for more than a single excitation makes the reaction of the involved polaritons entropically more favorable. Open-system simulations highlight how the coherence of the system's initial state governs decay pathways and enables long-lived energy storage. Finally, a multilevel vibrational model coupled to a two-mode cavity is introduced, where a Raman scheme enables selective mixing of vibrational states. The simulation results show that continuous pumping can drive molecules into specific vibrational states while suppressing electronic excitation, both for single molecules and small ensembles. Together, these results provide new theoretical insights into the mechanisms by which optical cavities reshape molecular energy landscapes. The results suggest approaches for controlling photonic, electronic, and vibrational dynamics in strongly coupled systems, contributing to the theoretical foundation of polariton chemistry.
Addressing needs of contemporary nonlinear femtosecond optical spectroscopy, we have developed a fully quantum, numerically accurate wave function-based approach for the calculation of third-order spectroscopic signals of polyatomic molecules and molecular aggregates at finite temperature. The systems are described by multimode nonadiabatic vibronic-coupling Hamiltonians, in which diagonal terms are treated in harmonic approximation, while off-diagonal interstate couplings are assumed to be coordinate independent. The approach is based on the Thermo Field Dynamics (TFD) representation of quantum mechanics and tensor-train (TT) machinery for efficient numerical simulation of quantum evolution of systems with many degrees of freedom. The developed TFD-TT approach is applied to the calculation of time- and frequency-resolved fluorescence spectra of the Fenna-Matthews-Olson (FMO) antenna complex at room temperature taking into account finite time-frequency resolution in fluorescence detection, orientational averaging, and static disorder.