A time-dependent multiconfigurational self-consistent field theory is presented to describe the many-body dynamics of a gas of identical bosonic atoms confined to an external trapping potential at zero temperature from first principles. A set of generalized evolution equations are developed, through the time-dependent variational principle, which account for the complete and self-consistent coupling between the expansion coefficients of each configuration and the underlying one-body wave functions within a restricted two state Fock space basis that includes the full effects of the condensate's mean field as well as atomic correlation. The resulting dynamical equations are a classical Hamiltonian system and, by construction, form a well-defined initial value problem. They are implemented in an efficient numerical algorithm. An example is presented, highlighting the generality of the theory, in which the ballistic expansion of a fragmented condensate ground state is compared to that of a macroscopic quantum superposition state, taken here to be a highly entangled number state, upon releasing the external trapping potential. Strikingly different many-body matter-wave dynamics emerge in each case, accentuating the role of both atomic correlation and mean-field effects in the two condensate states.
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.