TY - JOUR

T1 - Multiconfigurational time-dependent Hartree method for bosons

T2 - Many-body dynamics of bosonic systems

AU - Alon, Ofir E.

AU - Streltsov, Alexej I.

AU - Cederbaum, Lorenz S.

PY - 2008/3/14

Y1 - 2008/3/14

N2 - 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ödinger 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.

AB - 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ödinger 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.

UR - http://www.scopus.com/inward/record.url?scp=40949158372&partnerID=8YFLogxK

U2 - 10.1103/PhysRevA.77.033613

DO - 10.1103/PhysRevA.77.033613

M3 - Article

AN - SCOPUS:40949158372

SN - 1050-2947

VL - 77

JO - Physical Review A - Atomic, Molecular, and Optical Physics

JF - Physical Review A - Atomic, Molecular, and Optical Physics

IS - 3

M1 - 033613

ER -