The dynamics of a coupled electron–boson system is investigated by employing a multitude of the Davydov D 1 trial states, also known as the multi-D 1 Ansatz, and a second trial state based on a superposition of the time-dependent generalized coherent state (GCS Ansatz ). The two Ansätze are applied to study population dynamics in the spin-boson model and the Holstein molecular crystal model, and a detailed comparison with numerically exact results obtained by the (multilayer) multiconfiguration time-dependent Hartree method and the hierarchy equations of motion approach is drawn. It is found that the two methodologies proposed here have significantly improved over that with the single D 1 Ansatz, yielding quantitatively accurate results even in the critical cases of large energy biases and large transfer integrals. The two methodologies provide new effective tools for accurate, efficient simulation of many-body quantum dynamics thanks to a relatively small number of parameters which characterize the electron–nuclear wave functions. The wave-function-based approaches are capable of tracking explicitly detailed bosonic dynamics, which is absent by construct in approaches based on the reduced density matrix. The efficiency and flexibility of our methods are also advantages as compared with numerically exact approaches such as QUAPI and HEOM, especially at low temperatures and in the strong coupling regime.
The dynamics of a coupled electron–boson system is investigated by employing a multitude of the Davydov D 1 trial states, also known as the multi-D 1 Ansatz, and a second trial state based on a superposition of the time-dependent generalized coherent state (GCS Ansatz ). The two Ansätze are applied to study population dynamics in the spin-boson model and the Holstein molecular crystal model, and a detailed comparison with numerically exact results obtained by the (multilayer) multiconfiguration time-dependent Hartree method and the hierarchy equations of motion approach is drawn. It is found that the two methodologies proposed here have significantly improved over that with the single D 1 Ansatz, yielding quantitatively accurate results even in the critical cases of large energy biases and large transfer integrals. The two methodologies provide new effective tools for accurate, efficient simulation of many-body quantum dynamics thanks to a relatively small number of parameters which characterize the electron–nuclear wave functions. The wave-function-based approaches are capable of tracking explicitly detailed bosonic dynamics, which is absent by construct in approaches based on the reduced density matrix. The efficiency and flexibility of our methods are also advantages as compared with numerically exact approaches such as QUAPI and HEOM, especially at low temperatures and in the strong coupling regime.
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.
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.
Dynamics of the two-spin spin-boson model in the presence of Ohmic and sub-Ohmic baths is investigated by employing a multitude of the Davydov D1 trial states, also known as the multi-D1 Ansatz. Its accuracy in dynamics simulations of the two-spin SBM is improved significantly over the single D1 Ansatz, especially in the weak to moderately strong coupling regime. Validity of the multi-D1 Ansatz for various coupling strengths is also systematically examined by making use of the deviation vector which quantifies how faithfully the trial state obeys the Schrödinger equation. The time evolution of population difference and entanglement has been studied for various initial conditions and coupling strengths. Careful comparisons are carried out between our approach and three other methods, i.e., the time-dependent numerical renormalization group (TD-NRG) approach, the Bloch-Redfield theory, and a method based on a variational master equation. For strong coupling, the multi-D1 trial state yields consistent results as the TD-NRG approach in the Ohmic regime while the two disagree in the sub-Ohmic regime, where the multi-D1 trial state is shown to be more accurate. For weak coupling, the multi-D1 trial state agrees with the two master-equation methods in the presence of both Ohmic and sub-Ohmic baths, but shows considerable differences with the TD-NRG approach in the presence of a sub-Ohmic bath, calling into question the validity of the TD-NRG results at long times in the literature.
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.
Dynamics of the sub-Ohmic spin-boson model is examined using three numerical approaches, namely the Dirac-Frenkel time-dependent variation with the Davydov D(1) ansatz, the adaptive time-dependent density matrix renormalization group method within the representation of orthogonal polynomials, and a perturbative approach based on a unitary transformation. In order to probe the validity regimes of the three approaches, we study the dynamics of a qubit coupled to a bosonic bath with and without a local field. Comparison of the up-state population evolution shows that the three approaches are in agreement in the weak-coupling regime but exhibit marked differences when the coupling strength is large. The Davydov D(1) ansatz and the time-dependent density matrix renormalization group can both be reliably employed in the weak-coupling regime, while the former is also valid in the strong-coupling regime as judged by how faithfully the trial state follows the Schrödinger equation. We further explore the bipartite entanglement dynamics between two qubits coupled with individual bosonic baths which reveals entanglement sudden death and revival.
Abstract We propose a new version of the Davydov Ansatz for the approximate solution of the time-dependent Schrödinger equation (TDSE) and apply it to a spin 1/2 system coupled to a single bosonic mode. In complexity the new Ansatz is between the well-studied D1 and D2 Ansätze and is therefore termed D1.5 Ansatz. It is tailored to specific but still rather general initial conditions and is free of numerical inconsistencies that plague the D1 Ansatz. Finally, a generalization to nontrivial multiplicity in the oscillator state space is given that leads, in principle, to an exact solution of the TDSE for arbitrary system parameters. By employing Bargmann coherent states, the equations of motion could be made explicit, which is a big advantage, computationally.