Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:1512.05161 · The Journal of Chemical Physics

Variational dynamics of the sub-Ohmic spin-boson model on the basis of multiple Davydov D1 states

Lu Wang, Lipeng Chen, Nengji Zhou, Yang Zhao

Published Jan 8, 2016• 75 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1063/1.4939144

Related papers

Dynamics of the spin-boson model: The effect of bath initial conditions

Lipeng Chen, Yiying Yan, Maxim F. Gelin +1 more

The dynamics of the (sub-)Ohmic spin-boson model under various bath initial conditions is investigated by employing the Dirac–Frenkel time-dependent variational approach with the multiple Davydov D1 Ansatz in the interaction picture. The validity of our approach is carefully checked by comparing the results with those of the hierarchy equations of motion method. By analyzing the features of nonequilibrium dynamics, we identify the phase diagrams for different bath initial conditions. We find that for the spectral exponent s < sc, there exists a transition from coherent to quasicoherent dynamics with increasing coupling strengths. For sc < s ≤ 1, the coherent to incoherent crossover occurs at a certain coupling strength and the quasicoherent dynamics emerges at much larger couplings. The initial preparation of the bath has a considerable influence on the dynamics.

arXiv:2212.05723•The Journal of Chemical Physics•Feb 21, 2023•16 citations

Dynamics of the two-spin spin-boson model with a common bath

Tianrui Deng, Yiying Yan, Lipeng Chen +1 more

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 Journal of Chemical Physics•Apr 8, 2016•42 citations

Dynamics of the sub-Ohmic spin-boson model: A time-dependent variational study

Ning Wu, Liwei Duan, Xin Li +1 more

The Dirac-Frenkel time-dependent variation is employed to probe the dynamics of the zero temperature sub-Ohmic spin-boson model with strong friction utilizing the Davydov D1 ansatz. It is shown that initial conditions of the phonon bath have considerable influence on the dynamics. Counterintuitively, even in the very strong coupling regime, quantum coherence features still manage to survive under the polarized bath initial condition, while such features are absent under the factorized bath initial condition. In addition, a coherent-incoherent transition is found at a critical coupling strength α ≈ 0.1 for s = 0.25 under the factorized bath initial condition. We quantify how faithfully our ansatz follows the Schrödinger equation, finding that the time-dependent variational approach is robust for strong dissipation and deep sub-Ohmic baths (s ≪ 1).

arXiv:1302.1682•The Journal of Chemical Physics•Feb 28, 2013•42 citations

Dynamics of Coupled Electron–Boson Systems with the Multiple Davydov D1 Ansatz and the Generalized Coherent State

Lipeng Chen, Raffaele Borrelli, Yang Zhao

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 Journal of Physical Chemistry A•Oct 24, 2017•36 citations

Dynamics of the sub-Ohmic spin-boson model: A comparison of three numerical approaches

Yao Yao, Liwei Duan, Zhiguo Lü +2 more

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.

Physical Review E•Aug 16, 2013•44 citations

Quantum phase transition of sub-Ohmic spin-boson models: An approach by the multiple Davydov D2 Ansatz

Junzhe Tan, Nengji Zhou, Yang Zhao

The ground state properties and quantum phase transitions of sub-Ohmic spin-boson models are investigated using the multiple Davydov D2Ansatz in conjunction with the variational principle. Three variants of the model are studied: (i) a single bath with diagonal coupling, (ii) two independent baths with diagonal and off-diagonal couplings, and (iii) a single bath with simultaneous diagonal and off-diagonal couplings. For the purely diagonal model, the multiple Davydov D2Ansatz yields critical coupling strengths that are consistent with other methodologies, validating its accuracy and efficiency. In the two-bath model, the competition between diagonal and off-diagonal couplings drives a first-order transition for both symmetric and asymmetric spectral exponents, with the von Neumann exhibiting a continuous peak only under exact symmetry. Finally, for a single bath with simultaneous diagonal and off-diagonal couplings, we demonstrate that a rotational transformation maps the system to an equivalent purely diagonal model, enabling simpler and intuitive physical interpretation and reduced computational complexity.

arXiv:2510.11130•The Journal of Chemical Physics•Nov 4, 2025•2 citations
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.
The dynamics of the (sub-)Ohmic spin-boson model under various bath initial conditions is investigated by employing the Dirac–Frenkel time-dependent variational approach with the multiple Davydov D1 Ansatz in the interaction picture. The validity of our approach is carefully checked by comparing the results with those of the hierarchy equations of motion method. By analyzing the features of nonequilibrium dynamics, we identify the phase diagrams for different bath initial conditions. We find that for the spectral exponent s < sc, there exists a transition from coherent to quasicoherent dynamics with increasing coupling strengths. For sc < s ≤ 1, the coherent to incoherent crossover occurs at a certain coupling strength and the quasicoherent dynamics emerges at much larger couplings. The initial preparation of the bath has a considerable influence on the dynamics.
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 variation is employed to probe the dynamics of the zero temperature sub-Ohmic spin-boson model with strong friction utilizing the Davydov D1 ansatz. It is shown that initial conditions of the phonon bath have considerable influence on the dynamics. Counterintuitively, even in the very strong coupling regime, quantum coherence features still manage to survive under the polarized bath initial condition, while such features are absent under the factorized bath initial condition. In addition, a coherent-incoherent transition is found at a critical coupling strength α ≈ 0.1 for s = 0.25 under the factorized bath initial condition. We quantify how faithfully our ansatz follows the Schrödinger equation, finding that the time-dependent variational approach is robust for strong dissipation and deep sub-Ohmic baths (s ≪ 1).
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 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.
The ground state properties and quantum phase transitions of sub-Ohmic spin-boson models are investigated using the multiple Davydov D2Ansatz in conjunction with the variational principle. Three variants of the model are studied: (i) a single bath with diagonal coupling, (ii) two independent baths with diagonal and off-diagonal couplings, and (iii) a single bath with simultaneous diagonal and off-diagonal couplings. For the purely diagonal model, the multiple Davydov D2Ansatz yields critical coupling strengths that are consistent with other methodologies, validating its accuracy and efficiency. In the two-bath model, the competition between diagonal and off-diagonal couplings drives a first-order transition for both symmetric and asymmetric spectral exponents, with the von Neumann exhibiting a continuous peak only under exact symmetry. Finally, for a single bath with simultaneous diagonal and off-diagonal couplings, we demonstrate that a rotational transformation maps the system to an equivalent purely diagonal model, enabling simpler and intuitive physical interpretation and reduced computational complexity.