Multiconfigurational time-dependent Hartree method for bosons with internal degrees of freedom: Theory and composite fragmentation of multicomponent Bose-Einstein condensates

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