Commun. Comput. Phys.,
A Uniformly Convergent Numerical Method for Singularly Perturbed Nonlinear Eigenvalue Problems
Weizhu Bao 1*, Ming-Huang Chai 21 Department of Mathematics and Center for Computational Science and Engineering, National University of Singapore, 117543, Singapore.
2 Department of Mathematics, National University of Singapore, 117543, Singapore; and NUS High School, National University of Singapore, 117543, Singapore.
Received 22 August 2007; Accepted (in revised version) 15 October 2007
Available online 27 February 2008
In this paper we propose a uniformly convergent numerical method for discretizing singularly perturbed nonlinear eigenvalue problems under constraints with applications in Bose-Einstein condensation and quantum chemistry. We begin with the time-independent Gross-Pitaevskii equation and show how to reformulate it into a singularly perturbed nonlinear eigenvalue problem under a constraint. Matched asymptotic approximations for the problem are presented to locate the positions and characterize the widths of boundary layers and/or interior layers in the solution. A uniformly convergent numerical method is proposed by using the normalized gradient flow and piecewise uniform mesh techniques based on the asymptotic approximations for the problem. Extensive numerical results are reported to demonstrate the effectiveness of our numerical method for the problems. Finally, the method is applied to compute ground and excited states of Bose-Einstein condensation in the semiclassical regime and some conclusive findings are reported.AMS subject classifications: 35P30, 35Q55, 65N25, 65N06, 81Q05, 81V45
Notice: Undefined variable: pac in /var/www/html/issue/abstract/readabs.php on line 164
Key words: Nonlinear eigenvalue problem, Bose-Einstein condensation, ground state, excited state, energy, chemical potential, piecewise uniform mesh.
Email: email@example.com (W. Z. Bao), firstname.lastname@example.org (M. H. Chai)