Quantum Monte Carlo Simulation of the Davydov Model

Matthias Troyer, Philipp Werner, Adolfo Avella +1 more

Lu Wang, Yuta Fujihashi, Lipeng Chen +1 more

The Dirac-Frenkel time-dependent variational approach with Davydov Ansätze is a sophisticated, yet efficient technique to obtain an accurate solution to many-body Schrödinger equations for energy and charge transfer dynamics in molecular aggregates and light-harvesting complexes. We extend this variational approach to finite temperature dynamics of the spin-boson model by adopting a Monte Carlo importance sampling method. In order to demonstrate the applicability of this approach, we compare calculated real-time quantum dynamics of the spin-boson model with that from numerically exact iterative quasiadiabatic propagator path integral (QUAPI) technique. The comparison shows that our variational approach with the single Davydov Ansätze is in excellent agreement with the QUAPI method at high temperatures, while the two differ at low temperatures. Accuracy in dynamics calculations employing a multitude of Davydov trial states is found to improve substantially over the single Davydov Ansatz, especially at low temperatures. At a moderate computational cost, our variational approach with the multiple Davydov Ansatz is shown to provide accurate spin-boson dynamics over a wide range of temperatures and bath spectral densities.

Hyeonseong Jung, Shiho Kim

Lu Wang, Lipeng Chen, Nengji Zhou +1 more

Dynamics of the sub-Ohmic spin-boson model is investigated by employing a multitude of the Davydov D1 trial states, also known as the multi-D1 Ansatz. Accuracy in dynamics simulations is improved significantly over the single D1 Ansatz, especially in the weak system-bath coupling regime. The reliability of the multi-D1 Ansatz for various coupling strengths and initial conditions is also systematically examined, with results compared closely with those of the hierarchy equations of motion and the path integral Monte Carlo approaches. In addition, a coherent-incoherent phase crossover in the nonequilibrium dynamics is studied through the multi-D1 Ansatz. The phase diagram is obtained with a critical point sc = 0.4. For sc < s < 1, the coherent-to-incoherent crossover occurs at a certain coupling strength, while the coherent state recurs at a much larger coupling strength. For s < sc, only the coherent phase exists.

Yang Zhao

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.

Matthew Clifford Wilson