Particle-number-conserving Bogoliubov method which demonstrates the validity of the time-dependent Gross-Pitaevskii equation for a highly condensed Bose gas

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 Schrodinger 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'' xi(x) and a phonon field operator chi(x) in the form psi(x) approximate to A{xi(x) + chi(x)/root N} when the field operator acts on the N particle subspace. It is then possible to expand the Hamiltonian in decreasing powers of root 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 xi(x) is a solution of the time-independent Gross-Pitaevskii equation for N particles and any chemical potential 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 root N for the equations of motion. In this expansion the condensate wave function has the time-dependent form xi(x,t), and the condition for the validity of the expansion is that xi(x,t) satisfies the time-dependent Gross-Pitaevskii equation partial derivative xi/partial derivative t = -(h(2)/2m)del(2) xi + V xi + Nu\xi\(2) 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 chi(x,t) = Sigma(k) xi(k)(x,t)alpha(k), where the operators alpha(k) are time-independent annihilation operators, and the state vector has a time evolution described by a Schrodinger 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 xi(x,t).