Fundamental gaps of the Gross-Pitaevskii equation with repulsive interaction

We study asymptotically and numerically the fundamental gaps (i.e. the difference between the first excited state and the ground state) in energy and chemical potential of the Gross-Pitaevskii equation (GPE) -- nonlinear Schrodinger equation with cubic nonlinearity -- with repulsive interaction under different trapping potentials including box potential and harmonic potential. Based on our asymptotic and numerical results, we formulate a gap conjecture on the fundamental gaps in energy and chemical potential of the GPE on bounded domains with the homogeneous Dirichlet boundary condition, and in the whole space with a convex trapping potential growing at least quadratically in the far field. We then extend these results to the GPE on bounded domains with either the homogeneous Neumann boundary condition or periodic boundary condition.Comment: 26 pages, 16 figure

A Regularized Newton Method for Computing Ground States of Bose-Einstein condensates

In this paper, we propose a regularized Newton method for computing ground states of Bose-Einstein condensates (BECs), which can be formulated as an energy minimization problem with a spherical constraint. The energy functional and constraint are discretized by either the finite difference, or sine or Fourier pseudospectral discretization schemes and thus the original infinite dimensional nonconvex minimization problem is approximated by a finite dimensional constrained nonconvex minimization problem. Then an initial solution is first constructed by using a feasible gradient type method, which is an explicit scheme and maintains the spherical constraint automatically. To accelerate the convergence of the gradient type method, we approximate the energy functional by its second-order Taylor expansion with a regularized term at each Newton iteration and adopt a cascadic multigrid technique for selecting initial data. It leads to a standard trust-region subproblem and we solve it again by the feasible gradient type method. The convergence of the regularized Newton method is established by adjusting the regularization parameter as the standard trust-region strategy. Extensive numerical experiments on challenging examples, including a BEC in three dimensions with an optical lattice potential and rotating BECs in two dimensions with rapid rotation and strongly repulsive interaction, show that our method is efficient, accurate and robust.Comment: 25 pages, 6 figure