Benchmarking multiconfigurational Hartree by the exact wavefunction of two harmonically trapped bosons with contact interaction

Yeongjin Gwak, Oleksandr V. Marchukov, Uwe R. Fischer

J. Dukelsky, G.G. Dussel, R.P.J. Perazzo +1 more

Yeongjin Gwak, Oleksandr V. Marchukov, Uwe R. Fischer

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.

Jacky J. Chong, Zehua Zhao

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.

Luca Salasnich