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:1305.3862 · The Journal of Chemical Physics

The multi-layer multi-configuration time-dependent Hartree method for bosons: Theory, implementation, and applications

Lushuai Cao, Sven Krönke, Oriol Vendrell, Peter Schmelcher

Published Oct 1, 2013• 139 citations
arXivHTMLPDF

Abstract

We develop the multi-layer multi-configuration time-dependent Hartree method for bosons (ML-MCTDHB), a variational numerically exact ab initio method for studying the quantum dynamics and stationary properties of general bosonic systems. ML-MCTDHB takes advantage of the permutation symmetry of identical bosons, which allows for investigations of the quantum dynamics from few to many-body systems. Moreover, the multi-layer feature enables ML-MCTDHB to describe mixed bosonic systems consisting of arbitrary many species. Multi-dimensional as well as mixed-dimensional systems can be accurately and efficiently simulated via the multi-layer expansion scheme. We provide a detailed account of the underlying theory and the corresponding implementation. We also demonstrate the superior performance by applying the method to the tunneling dynamics of bosonic ensembles in a one-dimensional double well potential, where a single-species bosonic ensemble of various correlation strengths and a weakly interacting two-species bosonic ensemble are considered.

DOI
10.1063/1.4821350

Related papers

Non-equilibrium quantum dynamics of ultra-cold atomic mixtures: the multi-layer multi-configuration time-dependent Hartree method for bosons

Sven Krönke, Lushuai Cao, Oriol Vendrell +1 more

We develop and apply the multi-layer multi-configuration time-dependent Hartree method for bosons, which represents an ab initio method for investigating the non-equilibrium quantum dynamics of multi-species bosonic systems. Its multi-layer feature allows for tailoring the wave function ansatz to describe intra- and inter-species correlations accurately and efficiently. To demonstrate the beneficial scaling and efficiency of the method, we explored the correlated tunneling dynamics of two species with repulsive intra- and inter-species interactions, to which a third species with vanishing intra-species interaction was weakly coupled. The population imbalances of the first two species can feature a temporal equilibration and their time evolution significantly depends on the coupling to the third species. Bosons of the first and second species exhibit a bunching tendency, whose strength can be influenced by their coupling to the third species.

arXiv:1210.8063•New Journal of Physics•Jun 14, 2013•112 citations

A unified ab initio approach to the correlated quantum dynamics of ultracold fermionic and bosonic mixtures

L. Cao, V. Bolsinger, S. I. Mistakidis +4 more

We extent the recently developed Multi-Layer Multi-Configuration Time-Dependent Hartree method for Bosons for simulating the correlated quantum dynamics of bosonic mixtures to the fermionic sector and establish a unifying approach for the investigation of the correlated quantum dynamics of a mixture of indistinguishable particles, be it fermions or bosons. Relying on a multi-layer wave-function expansion, the resulting Multi-Layer Multi-Configuration Time-Dependent Hartree method for Mixtures (ML-MCTDHX) can be adapted to efficiently resolve system-specific intra- and inter-species correlations. The versatility and efficiency of ML-MCTDHX are demonstrated by applying it to the problem of colliding few-atom mixtures of both Bose-Fermi and Fermi-Fermi types. Thereby, we elucidate the role of correlations in the transmission and reflection properties of the collisional events. In particular, we present examples where the reflection (transmission) at the other atomic species is a correlation-dominated effect, i.e., it is suppressed in the mean-field approximation.

arXiv:1704.03810•The Journal of Chemical Physics•Jul 26, 2017•134 citations

Multiconfigurational time-dependent Hartree method for bosons: Many-body dynamics of bosonic systems

Ofir E. Alon, Alexej I. Streltsov, L. S. Cederbaum

The evolution of Bose-Einstein condensates is amply described by the time-dependent Gross-Pitaevskii mean-field theory which assumes all bosons to reside in a single time-dependent one-particle state throughout the propagation process. In this work, we go beyond mean field and develop an essentially exact many-body theory for the propagation of the time-dependent Schr\"odinger equation of $N$ interacting identical bosons. In our theory, the time-dependent many-boson wave function is written as a sum of permanents assembled from orthogonal one-particle functions, or orbitals, where both the expansion coefficients and the permanents (orbitals) themselves are time-dependent and fully determined according to a standard time-dependent variational principle. By employing either the usual Lagrangian formulation or the Dirac-Frenkel variational principle we arrive at two sets of coupled equations of motion, one for the orbitals and one for the expansion coefficients. The first set comprises of first-order differential equations in time and nonlinear integrodifferential equations in position space, whereas the second set consists of first-order differential equations with time-dependent coefficients. We call our theory multiconfigurational time-dependent Hartree for bosons, or $\text{MCTDHB}(M)$, where $M$ specifies the number of time-dependent orbitals used to construct the permanents. Numerical implementation of the theory is reported and illustrative numerical examples of many-body dynamics of trapped Bose-Einstein condensates are provided and discussed. The convergence of the method with a growing number $M$ of orbitals is demonstrated in a specific example of four interacting bosons in a double well.

arXiv:cond-mat/0703237•Physical Review A•Mar 14, 2008•355 citations

Beyond mean-field dynamics of ultra-cold bosonic atoms in higher dimensions: facing the challenges with a multi-configurational approach

Valentin Bolsinger, Sven Krönke, Peter Schmelcher

Abstract Exploring the impact of dimensionality on the quantum dynamics of interacting bosons in traps including particle correlations is an interesting but challenging task. Due to the different participating length scales, the modelling of the short-range interactions in three dimensions plays a special role. We review different approaches for the latter and elaborate that for multi-configurational computational strategies, finite-range potentials are adequate resulting in the need for large grids to resolve the relevant length scales. This results in computational challenges, which include the exponential scaling of complexity with the number of atoms. We show that the recently developed ab initio multi-layer multi-configurational time-dependent Hartee method for bosons (ML–MCTDHB) (2013 J. Chem. Phys. 139 134103 ) can face both numerical challenges and present an efficient numerical implementation of ML–MCTDHB in three spatial dimensions, particularly suited to describe the quantum dynamics for elongated traps. The beneficial scaling of our approach is demonstrated by studying the tunnelling dynamics of bosonic ensembles in a double well. Comparing three-dimensional with quasi-one dimensional simulations, we find dimensionality-induced effects in the density. Furthermore, we study the crossover from weak transversal confinement, where a mean-field description of the system is sufficient, towards tight transversal confinement, where particle correlations and beyond mean-field effects are pronounced.

arXiv:1608.04710•Journal of Physics B Atomic Molecular and Optical Physics•Jan 17, 2017•17 citations

Multiconfigurational time-dependent Hartree method for bosons with internal degrees of freedom: Theory and composite fragmentation of multicomponent Bose-Einstein condensates

Axel U. J. Lode

In this paper, the multiconfigurational time-dependent Hartree for bosons (MCTDHB) method is derived for the case of $N$ identical bosons with internal degrees of freedom. The theory for bosons with internal degrees of freedom constitutes a generalization of the MCTDHB method that substantially enriches the many-body physics that can be described. We demonstrate that the numerically exact solution of the time-dependent many-body Schr\"odinger equation for interacting bosonic particles with internal degrees of freedom is now feasible. We report on the MCTDHB equations of motion for bosons with internal degrees of freedom and their implementation for a general many-body Hamiltonian with one-body and two-body terms, both of which may depend on the internal states of the considered particles and time. To demonstrate the capabilities of the theory and its software implementation integrated in the mctdh-x software, we apply MCTDHB to the emergence of fragmentation of parabolically trapped bosons with two internal states: we study the ground state of $N=100$ bosons as a function of the separation between the state-dependent minima of the two parabolic potentials. To quantify the coherence of the system, we compute its normalized first-order correlation function. We find that the coherence within each internal state of the atoms is maintained, while it is lost between the different internal states. This is a hallmark of a kind of fragmentation absent in bosons without internal structure. We term the emergent phenomenon ``composite fragmentation.''

arXiv:1602.05791•Physical Review A•Jun 1, 2016•98 citations

Numerically exact quantum dynamics of bosons with time-dependent interactions of harmonic type

Axel U. J. Lode, Kaspar Sakmann, Ofir E. Alon +2 more

The exactly solvable quantum many-particle model with harmonic one- and two-particle interaction terms is extended to include time dependency. We show that when the external trap potential and interparticle interaction have a time dependency, the numerically exact solutions of the corresponding time-dependent many-boson Schr\"odinger equation are still available. We use these exact solutions to benchmark the recently developed multiconfigurational time-dependent Hartree method for bosons (MCTDHB) [Phys. Rev. Lett. 99, 030402 (2007); Phys. Rev. A 77, 033613 (2008)]. In particular, we benchmark the MCTDHB method for (i) the ground state; (ii) the breathing many-body dynamics activated by a quench scenario where the interparticle interaction strength is suddenly turned on to a finite value; (iii) the nonequilibrium dynamic for driven scenarios where both the trap- and interparticle-interaction potentials are time-dependent. Excellent convergence of the ground state and dynamics is demonstrated. The great relevance of the self-consistency and time adaptivity, which are the intrinsic features of the MCTDHB method, is demonstrated by contrasting the MCTDHB predictions and those obtained within the standard full configuration interaction method spanning a Fock space of the same size, but utilizing as one-particle basis set the fixed-shape eigenstates of the one-particle potential. Connections of the model's results to ultracold Bose-Einstein condensed systems are addressed.

arXiv:1207.5128•Physical Review A•Dec 6, 2012•117 citations
We develop the multi-layer multi-configuration time-dependent Hartree method for bosons (ML-MCTDHB), a variational numerically exact ab initio method for studying the quantum dynamics and stationary properties of general bosonic systems. ML-MCTDHB takes advantage of the permutation symmetry of identical bosons, which allows for investigations of the quantum dynamics from few to many-body systems. Moreover, the multi-layer feature enables ML-MCTDHB to describe mixed bosonic systems consisting of arbitrary many species. Multi-dimensional as well as mixed-dimensional systems can be accurately and efficiently simulated via the multi-layer expansion scheme. We provide a detailed account of the underlying theory and the corresponding implementation. We also demonstrate the superior performance by applying the method to the tunneling dynamics of bosonic ensembles in a one-dimensional double well potential, where a single-species bosonic ensemble of various correlation strengths and a weakly interacting two-species bosonic ensemble are considered.
We develop and apply the multi-layer multi-configuration time-dependent Hartree method for bosons, which represents an ab initio method for investigating the non-equilibrium quantum dynamics of multi-species bosonic systems. Its multi-layer feature allows for tailoring the wave function ansatz to describe intra- and inter-species correlations accurately and efficiently. To demonstrate the beneficial scaling and efficiency of the method, we explored the correlated tunneling dynamics of two species with repulsive intra- and inter-species interactions, to which a third species with vanishing intra-species interaction was weakly coupled. The population imbalances of the first two species can feature a temporal equilibration and their time evolution significantly depends on the coupling to the third species. Bosons of the first and second species exhibit a bunching tendency, whose strength can be influenced by their coupling to the third species.
We extent the recently developed Multi-Layer Multi-Configuration Time-Dependent Hartree method for Bosons for simulating the correlated quantum dynamics of bosonic mixtures to the fermionic sector and establish a unifying approach for the investigation of the correlated quantum dynamics of a mixture of indistinguishable particles, be it fermions or bosons. Relying on a multi-layer wave-function expansion, the resulting Multi-Layer Multi-Configuration Time-Dependent Hartree method for Mixtures (ML-MCTDHX) can be adapted to efficiently resolve system-specific intra- and inter-species correlations. The versatility and efficiency of ML-MCTDHX are demonstrated by applying it to the problem of colliding few-atom mixtures of both Bose-Fermi and Fermi-Fermi types. Thereby, we elucidate the role of correlations in the transmission and reflection properties of the collisional events. In particular, we present examples where the reflection (transmission) at the other atomic species is a correlation-dominated effect, i.e., it is suppressed in the mean-field approximation.
The evolution of Bose-Einstein condensates is amply described by the time-dependent Gross-Pitaevskii mean-field theory which assumes all bosons to reside in a single time-dependent one-particle state throughout the propagation process. In this work, we go beyond mean field and develop an essentially exact many-body theory for the propagation of the time-dependent Schr\"odinger equation of NNN interacting identical bosons. In our theory, the time-dependent many-boson wave function is written as a sum of permanents assembled from orthogonal one-particle functions, or orbitals, where both the expansion coefficients and the permanents (orbitals) themselves are time-dependent and fully determined according to a standard time-dependent variational principle. By employing either the usual Lagrangian formulation or the Dirac-Frenkel variational principle we arrive at two sets of coupled equations of motion, one for the orbitals and one for the expansion coefficients. The first set comprises of first-order differential equations in time and nonlinear integrodifferential equations in position space, whereas the second set consists of first-order differential equations with time-dependent coefficients. We call our theory multiconfigurational time-dependent Hartree for bosons, or MCTDHB(M)\text{MCTDHB}(M)MCTDHB(M), where MMM specifies the number of time-dependent orbitals used to construct the permanents. Numerical implementation of the theory is reported and illustrative numerical examples of many-body dynamics of trapped Bose-Einstein condensates are provided and discussed. The convergence of the method with a growing number MMM of orbitals is demonstrated in a specific example of four interacting bosons in a double well.
Abstract Exploring the impact of dimensionality on the quantum dynamics of interacting bosons in traps including particle correlations is an interesting but challenging task. Due to the different participating length scales, the modelling of the short-range interactions in three dimensions plays a special role. We review different approaches for the latter and elaborate that for multi-configurational computational strategies, finite-range potentials are adequate resulting in the need for large grids to resolve the relevant length scales. This results in computational challenges, which include the exponential scaling of complexity with the number of atoms. We show that the recently developed ab initio multi-layer multi-configurational time-dependent Hartee method for bosons (ML–MCTDHB) (2013 J. Chem. Phys. 139 134103 ) can face both numerical challenges and present an efficient numerical implementation of ML–MCTDHB in three spatial dimensions, particularly suited to describe the quantum dynamics for elongated traps. The beneficial scaling of our approach is demonstrated by studying the tunnelling dynamics of bosonic ensembles in a double well. Comparing three-dimensional with quasi-one dimensional simulations, we find dimensionality-induced effects in the density. Furthermore, we study the crossover from weak transversal confinement, where a mean-field description of the system is sufficient, towards tight transversal confinement, where particle correlations and beyond mean-field effects are pronounced.
In this paper, the multiconfigurational time-dependent Hartree for bosons (MCTDHB) method is derived for the case of NNN identical bosons with internal degrees of freedom. The theory for bosons with internal degrees of freedom constitutes a generalization of the MCTDHB method that substantially enriches the many-body physics that can be described. We demonstrate that the numerically exact solution of the time-dependent many-body Schr\"odinger equation for interacting bosonic particles with internal degrees of freedom is now feasible. We report on the MCTDHB equations of motion for bosons with internal degrees of freedom and their implementation for a general many-body Hamiltonian with one-body and two-body terms, both of which may depend on the internal states of the considered particles and time. To demonstrate the capabilities of the theory and its software implementation integrated in the mctdh-x software, we apply MCTDHB to the emergence of fragmentation of parabolically trapped bosons with two internal states: we study the ground state of N=100N=100N=100 bosons as a function of the separation between the state-dependent minima of the two parabolic potentials. To quantify the coherence of the system, we compute its normalized first-order correlation function. We find that the coherence within each internal state of the atoms is maintained, while it is lost between the different internal states. This is a hallmark of a kind of fragmentation absent in bosons without internal structure. We term the emergent phenomenon ``composite fragmentation.''
The exactly solvable quantum many-particle model with harmonic one- and two-particle interaction terms is extended to include time dependency. We show that when the external trap potential and interparticle interaction have a time dependency, the numerically exact solutions of the corresponding time-dependent many-boson Schr\"odinger equation are still available. We use these exact solutions to benchmark the recently developed multiconfigurational time-dependent Hartree method for bosons (MCTDHB) [Phys. Rev. Lett. 99, 030402 (2007); Phys. Rev. A 77, 033613 (2008)]. In particular, we benchmark the MCTDHB method for (i) the ground state; (ii) the breathing many-body dynamics activated by a quench scenario where the interparticle interaction strength is suddenly turned on to a finite value; (iii) the nonequilibrium dynamic for driven scenarios where both the trap- and interparticle-interaction potentials are time-dependent. Excellent convergence of the ground state and dynamics is demonstrated. The great relevance of the self-consistency and time adaptivity, which are the intrinsic features of the MCTDHB method, is demonstrated by contrasting the MCTDHB predictions and those obtained within the standard full configuration interaction method spanning a Fock space of the same size, but utilizing as one-particle basis set the fixed-shape eigenstates of the one-particle potential. Connections of the model's results to ultracold Bose-Einstein condensed systems are addressed.