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:1602.05791 · Physical Review A

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

Published Jun 1, 2016• 98 citations
arXivHTMLPDF

Abstract

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.''

DOI
10.1103/physreva.93.063601

Related papers

Vorticity, Variance, and the Vigor of Many-Body Phenomena in Ultracold Quantum Systems: MCTDHB and MCTDH-X

Ofir E. Alon, Raphael Beinke, L. S. Cederbaum +9 more

During the past year of the MCTDHB project at the HLRS, we continued to strive and conquest further applications, developments, and expansion of the MultiConfigurational Time-Dependent Hartree for Bosons (MCTDHB) method in the context of ultracold atomic systems. We also announce the MCTDH-X package, the Multiconfigurational Time-Dependent Hartree for Indistinguishable Particles X package, which is able to treat identical bosons and fermions, with or without spin/internal degrees of freedom, alike. Here we report on a plethora of results and versatile applications which include: (i) single-shot imaging of fluctuating vortices in a fragmented Bose-Einstein condensate (BEC); (ii) the many-body tunneling and fragmetnation of vortices in 2D trapped BECs; (iii) the transition from vortices to solitonic vortices in 2D trapped BECs; (iv) the variance of a many-particle system being very sensitive to correlations even in the infinite-particle limit; (v) the consequences of the latter on the out-of-equilibrium uncertainty product of an evolving BEC; (vi) the mechanism of tunneling to open space of a few interacting polarized fermions; and (vii) composite fragmentation of multi-components BECs (i.e., with internal degrees of freedom). These are all exciting results made throughout the allocation of computer time by the HLRS to the MCTDHB project. Finally, further perspectives and future research plans are briefly discussed.

Springer eBooks•Jan 1, 2016•4 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

Accurate multi-boson long-time dynamics in triple-well periodic traps

Alexej I. Streltsov, Kaspar Sakmann, Ofir E. Alon +1 more

In order to solve the many-boson Schr\"odinger equation we utilize the multiconfigurational time-dependent Hartree method for bosons (MCTDHB). To be able to attack larger systems and/or to propagate the solution for longer times, we implement a parallel version of the MCTDHB method, thereby realizing the recently proposed [Streltsov et al., Phys. Rev. A 81, 022124 (2010)] idea on how to construct efficiently the result of the action of the Hamiltonian on a bosonic state vector. As an illustrative example of its own interest, we study the real-space dynamics of repulsive bosonic systems made of $N=12$, 51, and 3003 bosons in triple-well periodic potentials. The ground state of this system is threefold fragmented. By suddenly strongly distorting the trap potential, the system performs complex many-body quantum dynamics. At long times it reveals a tendency to an oscillatory behavior around a threefold fragmented state. These oscillations are strongly suppressed and damped by quantum depletions. In spite of the richness of the observed dynamics, the three time-adaptive orbitals of MCTDHB($M=3$) are capable of describing the many-boson quantum dynamics of the system for short and intermediate times. For longer times, however, more self-consistent time-adaptive orbitals are needed to correctly describe the nonequilibrium many-body physics. The convergence of the MCTDHB($M$) method with the number $M$ of self-consistent time-dependent orbitals used is demonstrated.

arXiv:0910.5916•Physical Review A•Apr 6, 2011•50 citations

Multiconfigurational time-dependent Hartree method for fermions: Implementation, exactness, and few-fermion tunneling to open space

Elke Faßhauer, Axel U. J. Lode

We report on an implementation of the multiconfigurational time-dependent Hartree method (MCTDH) for spin-polarized fermions (MCTDHF). Our approach is based on a mapping for operators in Fock space that allows a compact and efficient application of the Hamiltonian and solution of the MCTDHF equations of motion. Our implementation extends, builds on, and exploits the recursive implementation of MCTDH for bosons (r-mctdhb) package. Together with r-mctdhb, the present implementation of MCTDHF forms the mctdh-x package. We benchmark the accuracy of the algorithm with the harmonic interaction model and a time-dependent generalization thereof. These models consider parabolically trapped particles that interact through a harmonic interaction potential. We demonstrate that MCTDHF is capable of solving the time-dependent many-fermion Schr\"odinger equation to an arbitrary degree of precision and can hence yield numerically exact results even in the case of Hamiltonians with time-dependent one-body and two-body potentials. We study the problem of two initially parabolically confined and charged fermions tunneling through a barrier to open space. We demonstrate the validity of a model proposed previously for the many-body tunneling to open space of bosonic particles with contact interactions [Proc. Natl. Acad. Sci. USA 109, 13521 (2012)]. The many-fermion tunneling can be built up from sequentially happening single-fermion tunneling processes. The characteristic momenta of these processes are determined by the chemical potentials of trapped subsystems of smaller particle numbers: The escaped fermions convert the different chemical potentials into kinetic energy. Using the two-body correlation function, we present a detailed picture of the sequentiality of the process and are able to tell tunneling from over-the-barrier escape.

arXiv:1510.02984•Physical Review A•Mar 21, 2016•102 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

Dynamics of Hubbard Hamiltonians with the multiconfigurational time-dependent Hartree method for indistinguishable particles

Axel U. J. Lode, Christoph Bruder

We apply the multiconfigurational time-dependent Hartree method for indistinguishable particles (MCTDH-X) to systems of bosons or fermions in lattices described by Hubbard-type Hamiltonians with long-range or short-range interparticle interactions. The wave function is expanded in a variationally optimized time-dependent many-body basis generated by a set of effective creation operators that are related to the original particle creation operators by a time-dependent unitary transform. We use the time-dependent variational principle for the coefficients of this transform and the expansion coefficients of the wave function expressed in the time-dependent many-body basis as variational parameters to derive equations of motion. The convergence of MCTDH-X is shown by comparing its results to the exact diagonalization of one-, two-, and three-dimensional lattices filled with bosons with contact interactions. We use MCTDH-X to study the buildup of correlations in the long-time splitting dynamics of a Bose-Einstein condensate loaded into a large two-dimensional lattice subject to a barrier that is ramped up in the center. We find that the system is split into two parts with emergent time-dependent correlations that depend on the ramping time---for most barrier-raising times the system becomes twofold fragmented, but for some of the very fast ramps, the system shows revivals of coherence.

arXiv:1604.08809•Physical Review A•Jul 22, 2016•38 citations
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.''
During the past year of the MCTDHB project at the HLRS, we continued to strive and conquest further applications, developments, and expansion of the MultiConfigurational Time-Dependent Hartree for Bosons (MCTDHB) method in the context of ultracold atomic systems. We also announce the MCTDH-X package, the Multiconfigurational Time-Dependent Hartree for Indistinguishable Particles X package, which is able to treat identical bosons and fermions, with or without spin/internal degrees of freedom, alike. Here we report on a plethora of results and versatile applications which include: (i) single-shot imaging of fluctuating vortices in a fragmented Bose-Einstein condensate (BEC); (ii) the many-body tunneling and fragmetnation of vortices in 2D trapped BECs; (iii) the transition from vortices to solitonic vortices in 2D trapped BECs; (iv) the variance of a many-particle system being very sensitive to correlations even in the infinite-particle limit; (v) the consequences of the latter on the out-of-equilibrium uncertainty product of an evolving BEC; (vi) the mechanism of tunneling to open space of a few interacting polarized fermions; and (vii) composite fragmentation of multi-components BECs (i.e., with internal degrees of freedom). These are all exciting results made throughout the allocation of computer time by the HLRS to the MCTDHB project. Finally, further perspectives and future research plans are briefly discussed.
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.
In order to solve the many-boson Schr\"odinger equation we utilize the multiconfigurational time-dependent Hartree method for bosons (MCTDHB). To be able to attack larger systems and/or to propagate the solution for longer times, we implement a parallel version of the MCTDHB method, thereby realizing the recently proposed [Streltsov et al., Phys. Rev. A 81, 022124 (2010)] idea on how to construct efficiently the result of the action of the Hamiltonian on a bosonic state vector. As an illustrative example of its own interest, we study the real-space dynamics of repulsive bosonic systems made of N=12N=12N=12, 51, and 3003 bosons in triple-well periodic potentials. The ground state of this system is threefold fragmented. By suddenly strongly distorting the trap potential, the system performs complex many-body quantum dynamics. At long times it reveals a tendency to an oscillatory behavior around a threefold fragmented state. These oscillations are strongly suppressed and damped by quantum depletions. In spite of the richness of the observed dynamics, the three time-adaptive orbitals of MCTDHB(M=3M=3M=3) are capable of describing the many-boson quantum dynamics of the system for short and intermediate times. For longer times, however, more self-consistent time-adaptive orbitals are needed to correctly describe the nonequilibrium many-body physics. The convergence of the MCTDHB(MMM) method with the number MMM of self-consistent time-dependent orbitals used is demonstrated.
We report on an implementation of the multiconfigurational time-dependent Hartree method (MCTDH) for spin-polarized fermions (MCTDHF). Our approach is based on a mapping for operators in Fock space that allows a compact and efficient application of the Hamiltonian and solution of the MCTDHF equations of motion. Our implementation extends, builds on, and exploits the recursive implementation of MCTDH for bosons (r-mctdhb) package. Together with r-mctdhb, the present implementation of MCTDHF forms the mctdh-x package. We benchmark the accuracy of the algorithm with the harmonic interaction model and a time-dependent generalization thereof. These models consider parabolically trapped particles that interact through a harmonic interaction potential. We demonstrate that MCTDHF is capable of solving the time-dependent many-fermion Schr\"odinger equation to an arbitrary degree of precision and can hence yield numerically exact results even in the case of Hamiltonians with time-dependent one-body and two-body potentials. We study the problem of two initially parabolically confined and charged fermions tunneling through a barrier to open space. We demonstrate the validity of a model proposed previously for the many-body tunneling to open space of bosonic particles with contact interactions [Proc. Natl. Acad. Sci. USA 109, 13521 (2012)]. The many-fermion tunneling can be built up from sequentially happening single-fermion tunneling processes. The characteristic momenta of these processes are determined by the chemical potentials of trapped subsystems of smaller particle numbers: The escaped fermions convert the different chemical potentials into kinetic energy. Using the two-body correlation function, we present a detailed picture of the sequentiality of the process and are able to tell tunneling from over-the-barrier escape.
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.
We apply the multiconfigurational time-dependent Hartree method for indistinguishable particles (MCTDH-X) to systems of bosons or fermions in lattices described by Hubbard-type Hamiltonians with long-range or short-range interparticle interactions. The wave function is expanded in a variationally optimized time-dependent many-body basis generated by a set of effective creation operators that are related to the original particle creation operators by a time-dependent unitary transform. We use the time-dependent variational principle for the coefficients of this transform and the expansion coefficients of the wave function expressed in the time-dependent many-body basis as variational parameters to derive equations of motion. The convergence of MCTDH-X is shown by comparing its results to the exact diagonalization of one-, two-, and three-dimensional lattices filled with bosons with contact interactions. We use MCTDH-X to study the buildup of correlations in the long-time splitting dynamics of a Bose-Einstein condensate loaded into a large two-dimensional lattice subject to a barrier that is ramped up in the center. We find that the system is split into two parts with emergent time-dependent correlations that depend on the ramping time---for most barrier-raising times the system becomes twofold fragmented, but for some of the very fast ramps, the system shows revivals of coherence.