The hierarchy of Davydov's Ansätze: From guesswork to numerically "exact" many-body wave functions

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.

Inspecting molecular aggregate quadratic vibronic coupling effects using squeezed coherent states

We present a systematic comparison of three quantum mechanical approaches describing excitation dynamics in molecular complexes using the time-dependent variational principle (TDVP) with increasing sophistication trial wavefunctions (ansatze): Davydov D₂, squeezed D₂ (sqD₂) and a numerically exact multiple D₂ (mD₂) ansatz in order to characterize validity of the sqD₂ ansatz. Numerical simulations of molecular aggregate absorption and fluorescence spectra with intra- and intermolecular vibrational modes, including quadratic electronic-vibrational (vibronic) coupling term, which is due to vibrational frequency shift upon pigment excitation are presented. Simulated absorption and fluorescence spectra of a J type molecular dimer with high frequency intramolecular vibrational modes obtained with D₂ and sqD₂ ansatze match the spectra of mD₂ ansatz only in the single pigment model without quadratic vibronic coupling. In general, the use of mD₂ ansatz is required to model an accurate dimer and larger aggregate's spectra. For a J dimer aggregate coupled to a low frequency intermolecular phonon bath, absorption and fluorescence spectra are qualitatively similar using all three ansatze. The quadratic vibronic coupling term in both absorption and fluorescence spectra manifests itself as a lineshape peak amplitude redistribution, static frequency shift and an additional shift, which is temperature dependent. Overall the squeezed D₂ model does not result in a considerable improvement of the simulation results compared to the simplest Davydov D₂ approach.

Simulation of Ab Initio Optical Absorption Spectrum of β-Carotene with Fully Resolved S0 and S2 Vibrational Normal Modes

The electronic absorption spectrum of β-carotene (β-Car) is studied using quantum chemistry and quantum dynamics simulations. Vibrational normal modes were computed in optimized geometries of the electronic ground state S₀ and the optically bright excited S₂ state using the time-dependent density functional theory. By expressing the S₂-state normal modes in terms of the ground-state modes, we find that no one-to-one correspondence between the ground- and excited-state vibrational modes exists. Using the ab initio results, we simulated the β-Car absorption spectrum with all 282 vibrational modes in a model solvent at 300 K using the time-dependent Dirac-Frenkel variational principle and are able to qualitatively reproduce the full absorption line shape. By comparing the 282-mode model with the prominent 2-mode model, widely used to interpret carotenoid experiments, we find that the full 282-mode model better describes the high-frequency progression of carotenoid absorption spectra; hence, vibrational modes become highly mixed during the S₀ → S₂ optical excitation. The obtained results suggest that electronic energy dissipation is mediated by numerous vibrational modes.

Vibration-mediated energy transport in bacterial reaction center: Simulation study

Exciton energy relaxation in a bacterial Reaction Center (bRC) pigment-protein aggregate presumably involves emission of high energy vibrational quanta to cover wide energy gaps between excitons. Here, we assess this hypothesis utilizing vibronic two-particle theory in modeling of the excitation relaxation process in bRC. Specific high frequency molecular vibrational modes are included explicitly one at a time in order to check which high frequency vibrations are involved in the excitation relaxation process. The low frequency bath modes are treated perturbatively within Redfield relaxation theory. The analysis of the population relaxation rate data indicates energy flow pathways in bRC and suggests that specific vibrations may be responsible for the excitation relaxation process.