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.
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.
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.''
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.
Abstract We introduce and describe the multiconfigurational time-depenent Hartree for indistinguishable particles (MCTDH-X) software, which is hosted, documented, and distributed at http://ultracold.org . This powerful tool allows the investigation of ground state properties and dynamics of interacting quantum many-body systems in different spatial dimensions. The MCTDH-X software is a set of programs and scripts to compute, analyze, and visualize solutions for the time-dependent and time-independent many-body Schrödinger equation for indistinguishable quantum particles. As the MCTDH-X software represents a general solver for the Schrödinger equation, it is applicable to a wide range of problems in the fields of atomic, optical, molecular physics, and condensed matter systems. In particular, it can be used to study light–matter interactions, correlated dynamics of electrons in the solid state as well as some aspects related to quantum information and computing. The MCTDH-X software solves a set of nonlinear coupled working equations based on the application of the time-dependent variational principle to the Schrödinger equation. These equations are obtained by using an ansatz for the many-body wavefunction that is a expansion in a set of time-dependent, fully symmetrized bosonic (X = B) or fully anti-symmetrized fermionic (X = F) many-body basis states. It is the time-dependence of the basis set that enables MCTDH-X to deal with quantum dynamics at a superior accuracy as compared to, for instance, exact diagonalization approaches with a static basis, where the number of basis states necessary to capture the dynamics of the wavefunction typically grows rapidly with time. Herein, we give an introduction to the MCTDH-X software via an easy-to-follow tutorial with a focus on accessibility. The illustrated exemplary problems are hosted at http://ultracold.org/tutorial and consider the physics of a few interacting bosons or fermions in a double-well potential. We explore computationally the position-space and momentum-space density, the one-body reduced density matrix, Glauber correlation functions, phases, (dynamical) phase transitions, and the imaging of the quantum systems in single-shot images. Although a few particles in a double well potential represent a minimal model system, we are able to demonstrate a rich variety of phenomena with it. We use the double well to illustrate the fermionization of bosonic particles, the crystallization of fermionic particles, characteristics of the superfluid and Mott-insulator quantum phases in Hubbard models, and even dynamical phase transitions. We provide a complete set of input files and scripts to redo all computations in this paper at http://ultracold.org/data/tutorial_input_files.zip , accompanied by tutorial videos at https://tinyurl.com/tjx35sq . Our tutorial should guide the potential users to apply the MCTDH-X software also to more complex systems.
We consider two bosons in a one-dimensional harmonic trap, interacting by a contact potential, and compare the exact solution of this problem to a self-consistent numerical solution by using the multiconfigurational time-dependent Hartree (MCTDH) method. We thereby benchmark the predictions of the MCTDH method, and the statement that it is numerically exact, with a few-body problem that has an analytical solution for the most commonly experimentally realized interaction potential in ultracold quantum gases. It is found that exact ground state energy and first order correlations are accurately reproduced by MCTDH up to the intermediate dimensionless coupling strengths corresponding to typical background scattering lengths of magnetically trapped ultracold dilute Bose gases. For larger couplings, established for example by (a combination of) Feshbach resonances and optical trapping, the MCTDH approach overestimates the depth of the trap-induced correlation dip of first order correlations in position space, and underestimates the fragmentation, defined as the average relative occupation of orbitals other than the energetically lowest one. We anticipate that qualitatively similar features in the correlation function may arise for larger particle numbers, paving the way for a quantitative assessment of the accuracy of MCTDH by experiments with ultracold atoms.
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 show that the successful and formally exact multiconfigurational time-dependent Hartree method (MCTDH) takes on a unified and compact form when specified for systems of identical particles (MCTDHF for fermions MCTDHB for bosons). In particular the equations of motion for the orbitals depend explicitly and solely on the reduced one- and two-body density matrices of the system's many-particle wave function. We point out that this appealing representation of the equations of motion opens up further possibilities for approximate propagation schemes.