Two-stage continuation algorithms for Bloch waves of Bose-Einstein condensates in optical lattices

Bloch waves of Bose-Einstein condensates (BEC) ill optical lattices are extremum nonlinear eigenstates which satisfy the time-independent Gross-Pitaevskii equation (GPE). We describe an efficient Taylor predictor-Newton corrector continuation algorithm for tracing solution curves of parameter-dependent problems. Based on this algorithm, a novel two-stage continuation algorithm is developed for computing Bloch waves of ID and 2D Bose-Einstein condensates (BEC) in optical lattices. We split the complex wave function into the sum of its real and imaginary parts. The original GPE becomes a couple of two nonlinear eigenvalue problems defined in the real domain with periodic boundary conditions. At the first stage we use the chemical potential mu as the continuation parameter. The Bloch wavenumber k (k(x,) k(y)), and the coefficient of the cubic term are treated as the second and third continuation parameters, respectively. Then we compute the Bloch bands/surfaces for the 1D/2D problem with linear counterparts. At the second stage we use mu and k/k(x) or k(y) as the continuation parameters simultaneously with two constraint conditions. The states without linear counterparts in the GPE can be obtained via states with linear counterparts. Numerical results are reported for both 1D and 2D problems. (c) 2010 Elsevier B.V. All rights reserved.http://gateway.webofknowledge.com/gateway/Gateway.cgi?GWVersion=2&SrcApp=PARTNER_APP&SrcAuth=LinksAMR&KeyUT=WOS:000281945600009&DestLinkType=FullRecord&DestApp=ALL_WOS&UsrCustomerID=8e1609b174ce4e31116a60747a720701Computer Science, Interdisciplinary ApplicationsPhysics, MathematicalSCI(E)EI8ARTICLE101727-173718