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.
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.
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.
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 linear-response theory of the multiconfigurational time-dependent Hartree for bosons method (LR-MCTDHB) for computing many-body excitations of trapped Bose-Einstein condensates [Phys. Rev. A 88 , 023606 (2013); J. Chem. Phys. 140 , 034108 (2014)] is implemented, for the first time, for systems with general interparticle interaction. This allows us to investigate the many-body excitation spectrum of interacting bosons with, for instance, a long-range interaction. Illustrative numerical examples for repulsive and attractive bosons are provided. The LR-MCTDHB theory is capable of identifying all excitations, including the excitations which are not unraveled within Bogoliubov-de Gennes equations. The theory is herewith benchmarked against the exactly-solvable one-dimensional harmonic-interaction model. As a complementary result, using a complex transformation, we represent LR-MCTDHB in a compact block-diagonal form, opening up thereby an avenue for treating larger many-body systems. We expect the LR-MCTDHB theory and its implementation for general interparticle interaction to provide a proved probe into the many-body excitations involved in the out-ofequilibrium dynamics of trapped interacting bosons.
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 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.
Following a ``bottom-up approach'' in understanding many-particle effects and dynamics we provide a systematic ab initio study of the dependence of the breathing dynamics of ultracold bosons in a one-dimensional (1D) harmonic trap on the number of bosons ranging from few to many. To this end, we employ the multilayer multiconfiguration time-dependent Hartree method for bosons (ML-MCTDHB) which has been developed very recently [Kr\"onke, Cao, Vendrell, and Schmelcher, New J. Phys. 15, 063018 (2013)]. The beating behavior for two bosons is found numerically and consequently explained by an analytical approach. Drawing on this, we show how to compute the complete breathing mode spectrum in this case. We examine how the two-mode breathing behavior of two bosons evolves to the single-frequency behavior of the many-particle limit when adding more particles. In the limit of many particles, we numerically study the dependence of the breathing mode frequency on both the interaction strength as well as on the particle number. We provide an estimate for the parameter region where the mean-field description provides a valid approximation.
We employ the multiconfigurational time-dependent Hartree for bosons (MCTDHB) method to study excited states of interacting Bose-Einstein condensates confined by harmonic and double-well trap potentials. Two approaches to access excitations, one static and the other dynamic, are investigated and contrasted. In static simulations the low-lying excitations are computed by utilizing a linear-response theory constructed on top of a static MCTDHB solution (LR-MCTDHB). Complimentarily, we propose two dynamic protocols that address excitations by propagating the MCTDHB wave function. In particular, we investigate dipolelike oscillations induced by shifting the origin of the confining potential and breathinglike excitations by quenching the frequency of a parabolic part of the trap. To contrast static predictions and dynamic results we compute the time evolution and regard the respective Fourier transform of several local and nonlocal observables. Namely, we study the expectation value of the position operator $\ensuremath{\langle}x(t)\ensuremath{\rangle}$, its variance Var[x(t)], and a local density computed at selected positions. We find that the variance is the most sensitive and informative quantity: Along with excitations it contains information about deexcitations even in a linear regime of the induced dynamics. The dynamic protocols are found to access the many-body excitations predicted by the static LR-MCTDHB approach.