**Principal Investigator:**

Axel U.J. Lode

**Affiliation:**

Department of Physics, University of Basel (Switzerland)

**Local Project ID:**

MCTDHB

**HPC Platform used:**

Hermit and Hornet of HLRS

**Date published:**

It has been a long-standing goal to understand the quantum physics of interacting quantum many-body systems. From the experimental side, the realization of Bose-Einstein condensates in 1995 provided scientists with such a system that is under almost perfect control: the confinement, the strength of interactions, and even the dimensionality of the Bose-Einstein condensate – a vapor of ultracold bosonic atoms held and cooled with lasers and magnetic fields – can be altered almost at will.

Recently, an in principle exact description of such systems of interacting ultracold bosonic atoms has become available with MCTDHB (multiconfigurational time-dependent Hartree method for bosons). MCTDHB basically is a clever way of translating the time-dependent many-boson Schrödinger equation which is governing the dynamics of interacting many-body systems to a form in which it is tractable with computers. The application and further development of MCTDHB is the goal of this project. The emergence of correlations seems to be a generic feature of interacting many-body systems – the interactions result in a collective behavior of the ultracold atoms that goes in sync with this emergence of correlations which requires many-body methods like MCTDHB.

The interacting quantum many-body problem is complicated to solve because of its high dimensionality: the dimension of the state of the system increases by the number of spatial dimensions for every particle added. For four particles in three spatial dimensions, a twelve-dimensional space is required – this is already at the border of todays supercomputers' capabilities. A direct solution of the problem is hence not possible for more than a few particles. With MCTDHB, a clever, time-dependent representation for the high dimensional many-body state has been found to describe up to a few thousand interacting quantum particles numerically exactly. The required computational resources are still substantial and the availability of supercomputing facilities such as the Cray XE6 „Hermit“ and Cray XC40 „Hornet“ at the HLRS Stuttgart has been crucial to obtain results with a larger number of atoms, to a higher degree of accuracy, for systems in two and three spatial dimensions and ultralong grids in one spatial dimension.

Recently, MCTDHB simulations enabled the discovery of a new way how angular momentum can be stored in interacting two-dimensional Bose-Einstein condensates: the phantom vortex. Strikingly, the transferred angular momentum manifests in the collective behavior of the atoms but is obscured from their density (See figure and videos). The obtained results on phantom vortices and the way they are created show that a theory that incorporates multiple superfluid order parameters instead of just one is likely to be necessary for a more concise description of the physics of superconducting materials in a magnetic field or the dynamics of superfluid Helium.