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 Chemical Physics

Symmetries in the multi-configurational time-dependent Hartree wavefunction representation and propagation

Thomas Weike, Uwe Manthe

Published May 20, 2021• 12 citations
Publisher

Abstract

In multi-configurational time-dependent Hartree (MCTDH) approaches, different multi-layered wavefunction representations can be used to represent the same physical wavefunction. Transformations between different equivalent representations of a physical wavefunction that alter the tree structure used in the multi-layer MCTDH wavefunction representation interchange the role of single-particle functions (SPFs) and single-hole functions (SHFs) in the MCTDH formalism. While the physical wavefunction is invariant under these transformations, this invariance does not hold for the standard multi-layer MCTDH equations of motion. Introducing transformed SPFs, which obey normalization conditions typically associated with SHFs, revised equations of motion are derived. These equations do not show the singularities resulting from the inverse single-particle density matrix and are invariant under tree transformations. Based on the revised equations of motion, a new integration scheme is introduced. The scheme combines the advantages of the constant mean-field approach of Beck and Meyer [Z. Phys. D 42, 113 (1997)] and the singularity-free integrator suggested by Lubich [Appl. Math. Res. Express 2015, 311]. Numerical calculations studying the spin boson model in high dimensionality confirm the favorable properties of the new integration scheme.

DOI
10.1063/5.0054105

Related papers

On the multi-layer multi-configurational time-dependent Hartree approach for bosons and fermions

Uwe Manthe, Thomas Weike

A multi-layer multi-configurational time-dependent Hartree (MCTDH) approach using a second quantization representation (SQR) based on optimized time-dependent orbitals is introduced. The approach combines elements of the multi-layer MCTDH-SQR approach of Wang and Thoss, which employs a preselected time-independent orbital basis, and the MCTDH for bosons and multi-configuration time-dependent Hartree-Fock approaches, which do not use multi-layering but employ time-dependent orbital bases. In contrast to existing MCTDH-type approaches, the results of the present approach for a given number of configurations are not invariant with respect to unitary transformations of the time-dependent orbital basis. Thus a natural orbital representation is chosen to achieve fast convergence with respect to the number of configurations employed. Equations of motion for the present ansatz, called (multi-layer) MCTDH in optimized second quantization representation, are derived. Furthermore, a scheme for the calculation of optimized unoccupied single-particle functions is given which can be used to avoid singularities in the equations of motion.

The Journal of Chemical Physics•Feb 14, 2017•39 citations

The multi-configurational time-dependent Hartree approach in optimized second quantization: Imaginary time propagation and particle number conservation

Thomas Weike, Uwe Manthe

The multilayer multiconfigurational time-dependent Hartree (MCTDH) in optimized second quantization representation (oSQR) approach combines the tensor contraction scheme of the multilayer MCTDH approach with the use of an optimized time-dependent orbital basis. Extending the original work on the subject [U. Manthe and T. Weike, J. Chem. Phys. 146, 064117 (2017)], here MCTDH-oSQR propagation in imaginary time and properties related to particle number conservation are studied. Differences between the orbital equation of motion in real and imaginary time are highlighted and a new gauge operator, which facilitates efficient imaginary time propagation, is introduced. Studying Bose-Hubbard models, particle number conservation in MCTDH-oSQR calculations is investigated in detail. Interesting properties of the single-particle functions used in the multilayer MCTDH representation are identified. Based on these results, a tensor contraction scheme, which explicitly utilizes particle number conservation, is suggested.

arXiv:1912.00954•The Journal of Chemical Physics•Jan 16, 2020•15 citations

A multi-layer multi-configurational time-dependent Hartree approach to lattice models beyond one dimension

Tristan Niermann, Hannes Hoppe, Uwe Manthe

The multi-layer multi-configurational time-dependent Hartree (MCTDH) approach is an efficient method to study quantum dynamics in real and imaginary time. The present work explores its potential to describe quantum fluids. The multi-layer MCTDH approach in second quantization representation is used to study lattice models beyond one dimension at finite temperatures. A scheme to map the lattice sites onto the MCTDH tree representation for multi-dimensional lattice models is proposed. A statistical sampling scheme previously used in MCTDH calculations is adapted to facilitate an efficient description of the thermal ensemble. As example, a two-dimensional hard-core Bose-Hubbard model is studied considering up to 64 × 64 lattice sites. The single particle function basis set size required to obtain converged results is found to not increase with the lattice size. The numerical results properly simulate the finite temperature Berezinskii-Kosterlitz-Thouless phase transition.

The Journal of Chemical Physics•Oct 3, 2024•3 citations

Time evolution of ML-MCTDH wavefunctions. II. Application of the projector splitting integrator

Lachlan P. Lindoy, Benedikt Kloss, David R. Reichman

The multi-layer multi-configuration time-dependent Hartree (ML-MCTDH) approach can suffer from numerical instabilities whenever the wavefunction is weakly entangled. These instabilities arise from singularities in the equations of motion (EOMs) and necessitate the use of regularization of the EOMs. The Projector Splitting Integrator (PSI) has previously been presented as an approach for evolving ML-MCTDH wavefunctions that is free of singularities. Here, we will discuss the implementation of the multi-layer PSI with a particular focus on how the steps required relate to those required to implement standard ML-MCTDH. We demonstrate the efficiency and stability of the PSI for large ML-MCTDH wavefunctions containing up to hundreds of thousands of nodes by considering a series of spin-boson models with up to 106 bath modes and find that for these problems, the PSI requires roughly 3–4 orders of magnitude fewer Hamiltonian evaluations and 2–3 orders of magnitude fewer Hamiltonian applications than standard ML-MCTDH and 2–3/1–2 orders of magnitude fewer evaluations/applications than approaches that use improved regularization schemes. Finally, we consider a series of significantly more challenging multi-spin-boson models that require much larger numbers of single-particle functions with wavefunctions containing up to ∼1.3×109 parameters to obtain accurate dynamics.

arXiv:2109.03134•The Journal of Chemical Physics•Oct 11, 2021•20 citations

Implementation of a novel projector-splitting integrator for the multi-configurational time-dependent Hartree approach

Benedikt Kloss, Irène Burghardt, Christian Lubich

The variational equations of motion of the multi-configuration time-dependent Hartree (MCTDH) approach contain the inverse of reduced density matrices which are typically ill-conditioned and therefore lead to small stepsizes for numerical time integration. This problem is usually dealt with via regularization of the density matrices, which works well in many cases but still calls for systematic improvement schemes. Recently this problem, its implications and possible solutions have become the subject of increased interest. Notably, a projector splitting integrator for the MCTDH approach that does not require the inversion of reduced density matrices has been proposed [C. Lubich, Appl. Math. Res. Express 2015, 311]. Here, we present the first implementation of this integration scheme. Results for low-dimensional benchmark systems are presented, and the case of initially unoccupied single-particle functions is discussed.

The Journal of Chemical Physics•May 2, 2017•42 citations

Colloquium: Multiconfigurational time-dependent Hartree approaches for indistinguishable particles

Axel U. J. Lode, Camille Lévêque, Lars Bojer Madsen +2 more

In this Colloquium, the wavefunction-based Multiconfigurational Time-Dependent Hartree approaches to the dynamics of indistinguishable particles (MCTDH-F for Fermions and MCTDH-B for Bosons) are reviewed. MCTDH-B and MCTDH-F or, together, MCTDH-X are methods for describing correlated quantum systems of identical particles by solving the time-dependent Schrödinger equation from first principles. MCTDH-X is used to accurately model the dynamics of real-world quantum many-body systems in atomic, molecular, and optical physics. The key feature of these approaches is the time-dependence and optimization of the single-particle states employed for the construction of a many-body basis set, which yields nonlinear working equations. We briefly describe the historical developments that have lead to the formulation of the MCTDH-X methods and motivate the necessity for wavefunction-based approaches. We sketch the derivation of the unified MCTDH-F and MCTDH-B equations of motion for complete and also specific restricted configuration spaces. The strengths and limitations of the MCTDH-X approach are assessed via benchmarks against an exactly solvable model and via convergence checks. We highlight some applications to instructive and experimentally-realized quantum many-body systems: the dynamics of atoms in Bose-Einstein condensates in magneto-optical and optical traps and of electrons in atoms and molecules. We discuss the current development and frontiers in the field of MCTDH-X: theories and numerical methods for indistinguishable particles, for mixtures of multiple species of indistinguishable particles, the inclusion of nuclear motion for the nonadia-batic dynamics of atomic and molecular systems, the so-called multilayer generalizations to the MCTDH-F and MCTDH-B methods, and the time-dependent orbital-adaptive coupled cluster theory are discussed.

•Nov 19, 2019•3 citations
In multi-configurational time-dependent Hartree (MCTDH) approaches, different multi-layered wavefunction representations can be used to represent the same physical wavefunction. Transformations between different equivalent representations of a physical wavefunction that alter the tree structure used in the multi-layer MCTDH wavefunction representation interchange the role of single-particle functions (SPFs) and single-hole functions (SHFs) in the MCTDH formalism. While the physical wavefunction is invariant under these transformations, this invariance does not hold for the standard multi-layer MCTDH equations of motion. Introducing transformed SPFs, which obey normalization conditions typically associated with SHFs, revised equations of motion are derived. These equations do not show the singularities resulting from the inverse single-particle density matrix and are invariant under tree transformations. Based on the revised equations of motion, a new integration scheme is introduced. The scheme combines the advantages of the constant mean-field approach of Beck and Meyer [Z. Phys. D 42, 113 (1997)] and the singularity-free integrator suggested by Lubich [Appl. Math. Res. Express 2015, 311]. Numerical calculations studying the spin boson model in high dimensionality confirm the favorable properties of the new integration scheme.
A multi-layer multi-configurational time-dependent Hartree (MCTDH) approach using a second quantization representation (SQR) based on optimized time-dependent orbitals is introduced. The approach combines elements of the multi-layer MCTDH-SQR approach of Wang and Thoss, which employs a preselected time-independent orbital basis, and the MCTDH for bosons and multi-configuration time-dependent Hartree-Fock approaches, which do not use multi-layering but employ time-dependent orbital bases. In contrast to existing MCTDH-type approaches, the results of the present approach for a given number of configurations are not invariant with respect to unitary transformations of the time-dependent orbital basis. Thus a natural orbital representation is chosen to achieve fast convergence with respect to the number of configurations employed. Equations of motion for the present ansatz, called (multi-layer) MCTDH in optimized second quantization representation, are derived. Furthermore, a scheme for the calculation of optimized unoccupied single-particle functions is given which can be used to avoid singularities in the equations of motion.
The multilayer multiconfigurational time-dependent Hartree (MCTDH) in optimized second quantization representation (oSQR) approach combines the tensor contraction scheme of the multilayer MCTDH approach with the use of an optimized time-dependent orbital basis. Extending the original work on the subject [U. Manthe and T. Weike, J. Chem. Phys. 146, 064117 (2017)], here MCTDH-oSQR propagation in imaginary time and properties related to particle number conservation are studied. Differences between the orbital equation of motion in real and imaginary time are highlighted and a new gauge operator, which facilitates efficient imaginary time propagation, is introduced. Studying Bose-Hubbard models, particle number conservation in MCTDH-oSQR calculations is investigated in detail. Interesting properties of the single-particle functions used in the multilayer MCTDH representation are identified. Based on these results, a tensor contraction scheme, which explicitly utilizes particle number conservation, is suggested.
The multi-layer multi-configurational time-dependent Hartree (MCTDH) approach is an efficient method to study quantum dynamics in real and imaginary time. The present work explores its potential to describe quantum fluids. The multi-layer MCTDH approach in second quantization representation is used to study lattice models beyond one dimension at finite temperatures. A scheme to map the lattice sites onto the MCTDH tree representation for multi-dimensional lattice models is proposed. A statistical sampling scheme previously used in MCTDH calculations is adapted to facilitate an efficient description of the thermal ensemble. As example, a two-dimensional hard-core Bose-Hubbard model is studied considering up to 64 × 64 lattice sites. The single particle function basis set size required to obtain converged results is found to not increase with the lattice size. The numerical results properly simulate the finite temperature Berezinskii-Kosterlitz-Thouless phase transition.
The multi-layer multi-configuration time-dependent Hartree (ML-MCTDH) approach can suffer from numerical instabilities whenever the wavefunction is weakly entangled. These instabilities arise from singularities in the equations of motion (EOMs) and necessitate the use of regularization of the EOMs. The Projector Splitting Integrator (PSI) has previously been presented as an approach for evolving ML-MCTDH wavefunctions that is free of singularities. Here, we will discuss the implementation of the multi-layer PSI with a particular focus on how the steps required relate to those required to implement standard ML-MCTDH. We demonstrate the efficiency and stability of the PSI for large ML-MCTDH wavefunctions containing up to hundreds of thousands of nodes by considering a series of spin-boson models with up to 106 bath modes and find that for these problems, the PSI requires roughly 3–4 orders of magnitude fewer Hamiltonian evaluations and 2–3 orders of magnitude fewer Hamiltonian applications than standard ML-MCTDH and 2–3/1–2 orders of magnitude fewer evaluations/applications than approaches that use improved regularization schemes. Finally, we consider a series of significantly more challenging multi-spin-boson models that require much larger numbers of single-particle functions with wavefunctions containing up to ∼1.3×109 parameters to obtain accurate dynamics.
The variational equations of motion of the multi-configuration time-dependent Hartree (MCTDH) approach contain the inverse of reduced density matrices which are typically ill-conditioned and therefore lead to small stepsizes for numerical time integration. This problem is usually dealt with via regularization of the density matrices, which works well in many cases but still calls for systematic improvement schemes. Recently this problem, its implications and possible solutions have become the subject of increased interest. Notably, a projector splitting integrator for the MCTDH approach that does not require the inversion of reduced density matrices has been proposed [C. Lubich, Appl. Math. Res. Express 2015, 311]. Here, we present the first implementation of this integration scheme. Results for low-dimensional benchmark systems are presented, and the case of initially unoccupied single-particle functions is discussed.
In this Colloquium, the wavefunction-based Multiconfigurational Time-Dependent Hartree approaches to the dynamics of indistinguishable particles (MCTDH-F for Fermions and MCTDH-B for Bosons) are reviewed. MCTDH-B and MCTDH-F or, together, MCTDH-X are methods for describing correlated quantum systems of identical particles by solving the time-dependent Schrödinger equation from first principles. MCTDH-X is used to accurately model the dynamics of real-world quantum many-body systems in atomic, molecular, and optical physics. The key feature of these approaches is the time-dependence and optimization of the single-particle states employed for the construction of a many-body basis set, which yields nonlinear working equations. We briefly describe the historical developments that have lead to the formulation of the MCTDH-X methods and motivate the necessity for wavefunction-based approaches. We sketch the derivation of the unified MCTDH-F and MCTDH-B equations of motion for complete and also specific restricted configuration spaces. The strengths and limitations of the MCTDH-X approach are assessed via benchmarks against an exactly solvable model and via convergence checks. We highlight some applications to instructive and experimentally-realized quantum many-body systems: the dynamics of atoms in Bose-Einstein condensates in magneto-optical and optical traps and of electrons in atoms and molecules. We discuss the current development and frontiers in the field of MCTDH-X: theories and numerical methods for indistinguishable particles, for mixtures of multiple species of indistinguishable particles, the inclusion of nuclear motion for the nonadia-batic dynamics of atomic and molecular systems, the so-called multilayer generalizations to the MCTDH-F and MCTDH-B methods, and the time-dependent orbital-adaptive coupled cluster theory are discussed.