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

The Journal of Physical Chemistry A

Iterative Calculation of Energy Eigenstates Employing the Multilayer Multiconfiguration Time-Dependent Hartree Theory

Haobin Wang

Published May 15, 2014• 32 citations
Publisher

Abstract

The improved relaxation approach of Meyer and co-workers is extended by employing the multilayer multiconfiguration time-dependent Hartree (ML-MCTDH) theory. Within this approach the trial wave function is expressed by a recursive, multilayered expansion that is the same as in a normal ML-MCTDH simulation. Calculation of an energy eigenstate then proceeds by iteratively diagonalizing the Boltzmann operator in the top layer representation using a Lanczos/Arnoldi method while relaxing the single particle functions of all layers using the ML-MCTDH imaginary time propagation. Numerical iteration is stopped until the relative energy and/or density of the target state meet the convergence criteria. The application of this multilayer improved relaxation method is illustrated by computing the energy splitting for the lowest pair of eigenstates as well as the spin observables and entanglement entropy of the spin-boson model. A large number of degrees of freedom, up to 10,000 bath modes, can be treated, and thus the appropriate phase transition at a critical system-bath coupling strength may be predicted via appropriate extrapolations.

DOI
10.1021/jp503351t

Related papers

Multilayer Multiconfiguration Time-Dependent Hartree Theory

Haobin Wang

Multilayer multiconfiguration time-dependent Hartree (ML-MCTDH) theory is a rigorous and powerful method to simulate quantum dynamics in complex many-body systems. This approach extends the original MCTDH theory of Meyer, Manthe, and Cederbaum to include dynamically contracted layers in a recursive way, within which the equations of motion are determined from the Dirac-Frenkel variational principle. This paper presents the general derivation of the theory and analyzes the important features that make the ML-MCTDH method numerically efficient. Furthermore, we discuss the generalization of the theory to treat many-body identical particles (fermions or bosons) as well as calculating energy eigenstates via the improved relaxation method.

The Journal of Physical Chemistry A•May 28, 2015•253 citations

Quantum Phase Transition in the Spin-Boson Model: A Multilayer Multiconfiguration Time-Dependent Hartree Study

Haobin Wang, Jiushu Shao

The multilayer improved relaxation is applied to study the delocalization-localization transition in the spin-boson model at zero temperature-a well-known example of quantum phase transition. Calculations of energy eigenstates are obtained by iteratively diagonalizing the matrix of the Boltzmann operator in the top layer representation, using a Lanczos/Arnoldi method while relaxing the single particle functions of all layers using the multilayer multiconfiguration time-dependent Hartree imaginary time propagation. Two properties are used to examine the quantum phase transition: the energy splitting for the lowest pair of eigenstates and the magnetic susceptibility. Consistent findings are obtained with appropriate scaling parameters.

The Journal of Physical Chemistry A•Feb 8, 2019•18 citations

Multiconfiguration time-dependent Hartree impurity solver for nonequilibrium dynamical mean-field theory

Karsten Balzer, Zheng Li, Oriol Vendrell +1 more

Nonequilibrium dynamical mean-field theory (DMFT) solves correlated lattice models by obtaining their local correlation functions from an effective model consisting of a single impurity in a self-consistently determined bath. The recently developed mapping of this impurity problem from the Keldysh time contour onto a time-dependent single-impurity Anderson model (SIAM) [C. Gramsch et al., Phys. Rev. B 88, 235106 (2013)] allows one to use wave-function-based methods in the context of nonequilibrium DMFT. Within this mapping, long times in the DMFT simulation become accessible by an increasing number of bath orbitals, which requires efficient representations of the time-dependent SIAM wave function. These can be achieved by the multiconfiguration time-dependent Hartree (MCTDH) method and its multilayer extensions. We find that MCTDH outperforms exact diagonalization for large baths in which the latter approach is still within reach and allows for the calculation of SIAMs beyond the system size accessible by exact diagonalization. Moreover, we illustrate the computation of the self-consistent two-time impurity Green's function within the MCTDH second quantization representation.

arXiv:1407.6578•Physical Review B•Jan 26, 2015•37 citations

Importance of Appropriately Regularizing the ML-MCTDH Equations of Motion

Haobin Wang, Hans‐Dieter Meyer

The multiconfiguration time-dependent Hartree (MCTDH) method and its generalization, the multilayer MCTDH (ML-MCTDH), result in equations of motion (EOMs) that are singular when there are virtual orbitals-the unoccupied single-particle functions-in the wave function expansion. For decades this singularity had been numerically removed by regularizing the reduced density matrix. In this Perspective we discuss our recent proposal to regularize the coefficient tensor instead, which has significant impact on both the efficiency and correctness of the EOMs in MCTDH and ML-MCTDH for challenging problems. We further demonstrate that when the system becomes large such that it is necessary to use ML-MCTDH with many layers, it is much more important to employ this new regularization scheme. We illustrate this point by studying a spin-boson model with a large bath that contains up to 100 000 modes. We show that even in the weak coupling regime the new regularization scheme is required to quickly rotate the virtual orbitals into the correct directions in Hilbert space. We argue that this situation can be common for applying a time-dependent tensor network approach to any large enough system.

The Journal of Physical Chemistry A•Feb 23, 2021•16 citations

Exploring disordered quantum spin models with a multilayer multiconfigurational approach

Fabian Köhler, Rick Mukherjee, Peter Schmelcher

Numerical simulations of quantum spin models are crucial for a profound understanding of many-body phenomena in a variety of research areas in physics. An outstanding problem is the availability of methods to tackle systems that violate area laws of entanglement entropy. Such scenarios cover a wide range of compelling physical situations including disordered quantum spin systems among others. In this paper, we employ a numerical technique referred to as multilayer multiconfiguration time-dependent Hartree (ML-MCTDH) to evaluate the ground state of several disordered spin models. ML-MCTDH has previously been used to study problems of high-dimensional quantum dynamics in molecular and ultracold physics but is here applied to study spin systems. We exploit the inherent flexibility of the method to present results in one and two spatial dimensions and treat challenging setups that incorporate long-range interactions as well as disorder. Our results suggest that the hierarchical multilayering inherent to ML-MCTDH allows to tackle a wide range of quantum many-body problems such as spin dynamics of varying dimensionality.

arXiv:2212.02254•Physical Review Research•May 30, 2023•4 citations

Accurate calculation of equilibrium reduced density matrix for the system-bath model: A multilayer multiconfiguration time-dependent Hartree approach and its comparison to a multi-electronic-state path integral molecular dynamics approach

Haobin Wang, Xinzijian Liu, Jian Liu

An efficient and accurate method for computing the equilibrium reduced density matrix is presented for treating open quantum systems characterized by the system-bath model. The method employs the multilayer multiconfiguration time-dependent Hartree theory for imaginary time propagation and an importance sampling procedure for calculating the quantum mechanical trace. The method is applied to the spin-boson Hamiltonian, which leads to accurate results in agreement with those produced by the multi-electronic-state path integral molecular dynamics method.

Chinese Journal of Chemical Physics•Aug 1, 2018•11 citations
The improved relaxation approach of Meyer and co-workers is extended by employing the multilayer multiconfiguration time-dependent Hartree (ML-MCTDH) theory. Within this approach the trial wave function is expressed by a recursive, multilayered expansion that is the same as in a normal ML-MCTDH simulation. Calculation of an energy eigenstate then proceeds by iteratively diagonalizing the Boltzmann operator in the top layer representation using a Lanczos/Arnoldi method while relaxing the single particle functions of all layers using the ML-MCTDH imaginary time propagation. Numerical iteration is stopped until the relative energy and/or density of the target state meet the convergence criteria. The application of this multilayer improved relaxation method is illustrated by computing the energy splitting for the lowest pair of eigenstates as well as the spin observables and entanglement entropy of the spin-boson model. A large number of degrees of freedom, up to 10,000 bath modes, can be treated, and thus the appropriate phase transition at a critical system-bath coupling strength may be predicted via appropriate extrapolations.
Multilayer multiconfiguration time-dependent Hartree (ML-MCTDH) theory is a rigorous and powerful method to simulate quantum dynamics in complex many-body systems. This approach extends the original MCTDH theory of Meyer, Manthe, and Cederbaum to include dynamically contracted layers in a recursive way, within which the equations of motion are determined from the Dirac-Frenkel variational principle. This paper presents the general derivation of the theory and analyzes the important features that make the ML-MCTDH method numerically efficient. Furthermore, we discuss the generalization of the theory to treat many-body identical particles (fermions or bosons) as well as calculating energy eigenstates via the improved relaxation method.
The multilayer improved relaxation is applied to study the delocalization-localization transition in the spin-boson model at zero temperature-a well-known example of quantum phase transition. Calculations of energy eigenstates are obtained by iteratively diagonalizing the matrix of the Boltzmann operator in the top layer representation, using a Lanczos/Arnoldi method while relaxing the single particle functions of all layers using the multilayer multiconfiguration time-dependent Hartree imaginary time propagation. Two properties are used to examine the quantum phase transition: the energy splitting for the lowest pair of eigenstates and the magnetic susceptibility. Consistent findings are obtained with appropriate scaling parameters.
Nonequilibrium dynamical mean-field theory (DMFT) solves correlated lattice models by obtaining their local correlation functions from an effective model consisting of a single impurity in a self-consistently determined bath. The recently developed mapping of this impurity problem from the Keldysh time contour onto a time-dependent single-impurity Anderson model (SIAM) [C. Gramsch et al., Phys. Rev. B 88, 235106 (2013)] allows one to use wave-function-based methods in the context of nonequilibrium DMFT. Within this mapping, long times in the DMFT simulation become accessible by an increasing number of bath orbitals, which requires efficient representations of the time-dependent SIAM wave function. These can be achieved by the multiconfiguration time-dependent Hartree (MCTDH) method and its multilayer extensions. We find that MCTDH outperforms exact diagonalization for large baths in which the latter approach is still within reach and allows for the calculation of SIAMs beyond the system size accessible by exact diagonalization. Moreover, we illustrate the computation of the self-consistent two-time impurity Green's function within the MCTDH second quantization representation.
The multiconfiguration time-dependent Hartree (MCTDH) method and its generalization, the multilayer MCTDH (ML-MCTDH), result in equations of motion (EOMs) that are singular when there are virtual orbitals-the unoccupied single-particle functions-in the wave function expansion. For decades this singularity had been numerically removed by regularizing the reduced density matrix. In this Perspective we discuss our recent proposal to regularize the coefficient tensor instead, which has significant impact on both the efficiency and correctness of the EOMs in MCTDH and ML-MCTDH for challenging problems. We further demonstrate that when the system becomes large such that it is necessary to use ML-MCTDH with many layers, it is much more important to employ this new regularization scheme. We illustrate this point by studying a spin-boson model with a large bath that contains up to 100 000 modes. We show that even in the weak coupling regime the new regularization scheme is required to quickly rotate the virtual orbitals into the correct directions in Hilbert space. We argue that this situation can be common for applying a time-dependent tensor network approach to any large enough system.
Numerical simulations of quantum spin models are crucial for a profound understanding of many-body phenomena in a variety of research areas in physics. An outstanding problem is the availability of methods to tackle systems that violate area laws of entanglement entropy. Such scenarios cover a wide range of compelling physical situations including disordered quantum spin systems among others. In this paper, we employ a numerical technique referred to as multilayer multiconfiguration time-dependent Hartree (ML-MCTDH) to evaluate the ground state of several disordered spin models. ML-MCTDH has previously been used to study problems of high-dimensional quantum dynamics in molecular and ultracold physics but is here applied to study spin systems. We exploit the inherent flexibility of the method to present results in one and two spatial dimensions and treat challenging setups that incorporate long-range interactions as well as disorder. Our results suggest that the hierarchical multilayering inherent to ML-MCTDH allows to tackle a wide range of quantum many-body problems such as spin dynamics of varying dimensionality.
An efficient and accurate method for computing the equilibrium reduced density matrix is presented for treating open quantum systems characterized by the system-bath model. The method employs the multilayer multiconfiguration time-dependent Hartree theory for imaginary time propagation and an importance sampling procedure for calculating the quantum mechanical trace. The method is applied to the spin-boson Hamiltonian, which leads to accurate results in agreement with those produced by the multi-electronic-state path integral molecular dynamics method.