We encode the many-body wave function of a Bose-Einstein condensate (BEC) in the $N$-particle sector of an extended catalytic state. This catalytic state is a coherent state for the condensate mode and an arbitrary state for the modes orthogonal to the condensate mode. Going to a time-dependent interaction picture where the state of the condensate mode is displaced to the vacuum, we can organize the effective Hamiltonian by powers of ${N}^{\ensuremath{-}1/2}$. Requiring the terms of order ${N}^{1/2}$ to vanish gives the Gross-Pitaevskii equation. Going to the next order, ${N}^{0}$, we derive equations for the number-conserving Bogoliubov approximation, first given by Castin and Dum [Phys. Rev. A 57, 3008 (1998)]. In contrast to other approaches, ours is well suited to calculating the state evolution in the Schr\"odinger picture; moreover, it is straightforward to generalize our method to multicomponent BECs and to higher-order corrections.
We encode the many-body wave function of a Bose-Einstein condensate (BEC) in the N-particle sector of an extended catalytic state. This catalytic state is a coherent state for the condensate mode and an arbitrary state for the modes orthogonal to the condensate mode. Going to a time-dependent interaction picture where the state of the condensate mode is displaced to the vacuum, we can organize the effective Hamiltonian by powers of N\ensuremath−1/2. Requiring the terms of order N1/2 to vanish gives the Gross-Pitaevskii equation. Going to the next order, N0, we derive equations for the number-conserving Bogoliubov approximation, first given by Castin and Dum [Phys. Rev. A 57, 3008 (1998)]. In contrast to other approaches, ours is well suited to calculating the state evolution in the Schr"odinger picture; moreover, it is straightforward to generalize our method to multicomponent BECs and to higher-order corrections.
The impact of interparticle correlations on the behavior of Bose-Einstein Condensates (BECs) is discussed using two approaches. In the first approach, the wavefunction of a BEC is encoded in the N-particle sector of an extended "catalytic state". Going to a time-dependent interaction picture, we can organize the effective Hamiltonian by powers of N−1/2. Requiring the terms of order N1/2 to vanish, we get the Gross-Pitaevskii Equation. Going to the next order, N0, we obtain the number-conserving Bogoliubov approximation. Our approach allows one to stay in the Schrödinger picture and to apply many techniques from quantum optics. Moreover, it is easier to track different orders in the Hamiltonian and to generalize to the multi-component case. In the second approach, I consider a state of N=l×n bosons that is derived by symmetrizing the n-fold tensor product of an arbitrary l-boson state. Particularly, we are interested in the pure state case for l=2, which we call the Pair-Correlated State (PCS). I show that PCS reproduces the number-conserving Bogoliubov approximation; moreover, it also works in the strong interaction regime where the Bogoliubov approximation fails. For the two-site Bose-Hubbard model, I find numerically that the error (measured by the trace distance of the two-particle reduced density matrices) of PCS is less than two percent over the entire parameter space, thus making PCS a bridge between the superfluid and Mott insulating phases. Amazingly, the error of PCS does not increase, in the time-dependent case, as the system evolves for longer times. I derive both time-dependent and -independent equations for the ground state and the time evolution of the PCS ansatz. The time complexity of simulating PCS does not depend on N and is linear in the number of orbitals in use.
The Bogoliubov method for the excitation spectrum of a Bose-condensed gas is generalized to apply to a gas with an exact large number N of particles. This generalization yields a description of the Schr\"odinger picture field operators as the product of an annihilation operator A for the total number of particles and the sum of a ``condensate wave function'' $\ensuremath{\xi}(\mathbf{x})$ and a phonon field operator $\ensuremath{\chi}(\mathbf{x})$ in the form $\ensuremath{\psi}(\mathbf{x})\ensuremath{\approx}A{\ensuremath{\xi}(\mathbf{x})+\ensuremath{\chi}(\mathbf{x})/\sqrt{N}}$ when the field operator acts on the N particle subspace. It is then possible to expand the Hamiltonian in decreasing powers of N, and thus obtain solutions for eigenvalues and eigenstates as an asymptotic expansion of the same kind. It is also possible to compute all matrix elements of field operators between states of different N. The excitation spectrum can be obtained by essentially the same method as Bogoliubov only if $\ensuremath{\xi}(\mathbf{x})$ is a solution of the time-independent Gross-Pitaevskii equation for N particles and any chemical potential $\ensuremath{\mu}$, which yields a valid and stable solution of the Gross-Pitaevskii equation. The treatment within a subspace of fixed N is identical in form to that usually used, but the interpretation of the operators is slightly different. A time-dependent generalization is then made, yielding an asymptotic expansion in decreasing powers of N for the equations of motion. In this expansion the condensate wave function has the time-dependent form $\ensuremath{\xi}(\mathbf{x},t)$, and the condition for the validity of the expansion is that $\ensuremath{\xi}(\mathbf{x},t)$ satisfies the time-dependent Gross-Pitaevskii equation $\ensuremath{\partial}\ensuremath{\xi}/\ensuremath{\partial}t=\ensuremath{-}({\ensuremath{\Elzxh}}^{2}/2m){\ensuremath{\nabla}}^{2}\ensuremath{\xi}+V\ensuremath{\xi}+Nu|\ensuremath{\xi}{|}^{2}\ensuremath{\xi}$. The physics is then described in a kind of interaction picture, called the condensate picture, in which the phonon operator can be expressed as $\ensuremath{\chi}(\mathbf{x},t)={\ensuremath{\sum}}_{k}{\ensuremath{\xi}}_{k}(\mathbf{x},t){\ensuremath{\alpha}}_{k},$ where the operators ${\ensuremath{\alpha}}_{k}$ are time-independent annihilation operators, and the state vector has a time evolution described by a Schr\"odinger equation in which the Hamiltonian is a time-dependent quadratic form in the phonon creation and annihilation operators, whose coefficients are explicitly determined in terms of the time-dependent condensate wave function $\ensuremath{\xi}(\mathbf{x},t)$.
We develop the number-conserving approach, which has previously been used in a single-component Bose-Einstein condensed dilute atomic gas, to describe consistent coupled condensate and noncondensate number dynamics, to an n-component condensate. The resulting system of equations is comprised, for each component, of a generalized Gross-Pitaevskii equation coupled to modified Bogoliubov--de Gennes equations. Lower order approximations yield general formulations for multicomponent Gross-Pitaevskii equations, and systems of multicomponent Gross-Pitaevskii equations coupled to multicomponent modified number-conserving Bogoliubov--de Gennes equations. The analysis is left general, such that, in the n-component condensate, there may or may not be mutually coherent components. An expansion in powers of the ratio of noncondensate-to-condensate particle numbers for each coherent set is used to derive the self-consistent, second-order, dynamical equations of motion. The advantage of the analysis developed in this article is in its applications to dynamical instabilities that appear when two (or more) components are in conflict and where a significant noncondensed fraction of atoms is expected to appear.
Starting from first-principle many-body quantum dynamics, we show that the dynamics of Bose-Einstein condensates can be approximated by the time-dependent nonlinear Gross-Pitaevskii equation, giving a bound on the rate of the convergence. Initial data are constructed on the bosonic Fock space applying an appropriate Bogoliubov transformation on a coherent state with expected number of particles N. The Bogoliubov transformation plays a crucial role; it produces the correct microscopic correlations among the particles. Our analysis shows that, on the level of the one-particle reduced density, the form of the initial data is preserved by the many-body evolution, up to a small error that vanishes as N−1/2 in the limit of large N.
A recent paper of Gardiner [Phys. Rev. A 56, 1414 (1997)] introduces a particle-number-conserving Bogoliubov method for the excitation spectrum of a Bose-condensed gas, for use in theories of recently experimentally produced trapped atomic Bose condensates. Gardiner's approach is compared and contrasted with the 1959 Girardeau-Arnowitt theory [Phys. Rev. 113, 755 (1959)], to which it is closely related and which is also fully number conserving. The number-conserving Bogoliubov quasiparticle operators of the Girardeau-Arnowitt theory satisfy Bose commutation relations exactly so long as states with the condensate totally depleted are neglected, whereas those of Gardiner satisfy Bose commutation relations only in an approximation that deteriorates progressively as the condensate is depleted.
We describe a number-conserving approach to the dynamics of Bose-Einstein condensed dilute atomic gases. This builds upon the works of Gardiner [Phys. Rev. A 56, 1414 (1997)] and Castin and Dum [Phys. Rev. A 57, 3008 (1998)]. We consider what is effectively an expansion in powers of the ratio of noncondensate to condensate particle numbers, rather than inverse powers of the total number of particles. This requires the number of condensate particles to be a majority, but not necessarily almost equal to the total number of particles in the system. We argue that a second-order treatment of the relevant dynamical equations of motion is the minimum order necessary to provide consistent coupled condensate and noncondensate number dynamics for a finite total number of particles, and show that such a second-order treatment is provided by a suitably generalized Gross-Pitaevskii equation, coupled to the Castin-Dum number-conserving formulation of the Bogoliubov--de Gennes equations. The necessary equations of motion can be generated from an approximate third-order Hamiltonian, which effectively reduces to second order in the steady state. Such a treatment as described here is suitable for dynamics occurring at finite temperature, where there is a significant noncondensate fraction from the outset, or dynamics leading to dynamical instabilities, where depletion of the condensate can also lead to a significant noncondensate fraction, even if the noncondensate fraction is initially negligible.