Second order averaging for the nonlinear Schroedinger equation with strongly anisotropic potential

International audienceWe consider the three dimensional Gross-Pitaevskii equation (GPE) describing a Bose-Einstein Condensate (BEC) which is highly confi ned in vertical z direction. The highly confi ned potential induces high oscillations in time. If the confi nement in the z direction is a harmonic trap (which is widely used in physical experiments), the very special structure of the spectrum of the confi nement operator will imply that the oscillations are periodic in time. Based on this observation, it can be proved that the GPE can be averaged out with an error of order of epsilon, which is the typical period of the oscillations. In this article, we construct a more accurate averaged model, which approximates the GPE up to errors of order epsilon squared. Then, expansions of this model over the eigenfunctions (modes) of the vertical Hamiltonian Hz are given in convenience of numerical application. Effi cient numerical methods are constructed for solving the GPE with cylindrical symmetry in 3D and the approximation model with radial symmetry in 2D, and numerical results are presented for various kinds of initial data