×

A normalized gradient flow method with attractive-repulsive splitting for computing ground states of Bose-Einstein condensates with higher-order interaction. (English) Zbl 1415.65198

Summary: In this paper, we generalize the normalized gradient flow method to compute the ground states of Bose-Einstein condensates (BEC) with higher order interactions (HOI), which is modeled via the modified Gross-Pitaevskii equation (MGPE). Schemes constructed in naive ways suffer from severe stability problems due to the high restrictions on time steps. To build an efficient and stable scheme, we split the HOI term into two parts with each part treated separately. The part corresponding to a repulsive/positive energy is treated semi-implicitly while the one corresponding to an attractive/negative energy is treated fully explicitly. Based on the splitting, we construct the BEFD-splitting and BESP-splitting schemes. A variety of numerical experiments show that the splitting will improve the stability of the schemes significantly. Besides, we will show that the methods can be applied to multidimensional problems and to the computation of the first excited state as well.

MSC:

65M06 Finite difference methods for initial value and initial-boundary value problems involving PDEs
65M70 Spectral, collocation and related methods for initial value and initial-boundary value problems involving PDEs
81V70 Many-body theory; quantum Hall effect
35Q40 PDEs in connection with quantum mechanics

Software:

Wesseling; GPELab
PDFBibTeX XMLCite
Full Text: DOI arXiv

References:

[1] Anderson, M. H.; Ensher, J. R.; Matthews, M. R.; Wieman, C. E.; Cornell, E. A., Observation of Bose-Einstein condensation in a dilute atomic vapor, Science, 269, 198-201 (1995)
[2] Antoine, X.; Besse, C.; Duboscq, R.; Rispoli, V., Acceleration of the imaginary time method for spectrally computing the stationary states of Gross-Pitaevskii equations, Comput. Phys. Commun., 219, 70-78 (2017) · Zbl 1411.81026
[3] Antoine, X.; Duboscq, R., GPELab, a Matlab toolbox to solve Gross-Pitaevskii equations I: computation of stationary solutions, Comput. Phys. Commun., 185, 2969-2991 (2014) · Zbl 1348.35003
[4] Antoine, X.; Duboscq, R., Robust and efficient preconditioned Krylov spectral solvers for computing the ground states of fast rotating and strongly interacting Bose-Einstein condensates, J. Comput. Phys., 258, 509-523 (2014) · Zbl 1349.82027
[5] Antoine, X.; Levitt, A.; Tang, Q., Efficient spectral computation of the stationary states of rotating Bose-Einstein condensates by the preconditioned nonlinear conjugate gradient method, J. Comput. Phys., 343, 92-109 (2017) · Zbl 1380.81496
[6] Bao, W., Ground states and dynamics of multicomponent Bose-Einstein condensates, Multiscale Model. Simul., 2, 210-236 (2004) · Zbl 1062.82034
[7] Bao, W., Mathematical models and numerical methods for Bose-Einstein condensation, (Proceeding of International Congress of Mathematicians, IV (2014)), 971-996 · Zbl 1373.35283
[8] Bao, W.; Cai, Y., Mathematical theory and numerical methods for Bose-Einstein condensation, Kinet. Relat. Models, 6, 1-135 (2012) · Zbl 1266.82009
[9] Bao, W.; Cai, Y., Ground states and dynamics of spin-orbit-coupled Bose-Einstein condensates, SIAM J. Appl. Math., 75, 492-517 (2015) · Zbl 1335.35227
[10] Bao, W.; Cai, Y.; Ruan, X., Ground states and dynamics of Bose-Einstein condensates with higher order interactions
[11] Bao, W.; Chern, I-L.; Zhang, Y., Efficient methods for computing ground states of spin-1 Bose-Einstein condensates based on their characterizations, J. Comput. Phys., 253, 189-208 (2013) · Zbl 1349.82069
[12] Bao, W.; Du, Q., Computing the ground state solution of Bose-Einstein condensates by a normalized gradient flow, SIAM J. Sci. Comput., 25, 1674-1697 (2004) · Zbl 1061.82025
[13] Bao, W.; Jaksch, D.; Markowich, P. A., Numerical solution of the Gross-Pitaevskii equation for Bose-Einstein condensation, J. Comput. Phys., 187, 318-342 (2003) · Zbl 1028.82501
[14] Bao, W.; Lim, F. Y.; Zhang, Y., Energy and chemical potential asymptotics for the ground state of Bose-Einstein condensates in the semiclassical regime, Bull. Inst. Math. Acad. Sin., 2, 495-532 (2007) · Zbl 1132.82003
[15] Bao, W.; Ruan, X., Fundamental gaps of the Gross-Pitaevskii equation with repulsive interaction · Zbl 1417.35175
[16] Bao, W.; Tang, W., Ground state solution of Bose-Einstein condensate by directly minimizing the energy functional, J. Comput. Phys., 187, 230-254 (2003) · Zbl 1028.82500
[17] Bao, W.; Wang, H.; Markowich, P. A., Ground, symmetric and central vortex states in rotating Bose-Einstein condensates, Commun. Math. Sci., 3, 57-88 (2005) · Zbl 1073.82004
[18] Bradley, C. C.; Sackett, C. A.; Tollett, J. J.; Hulet, R. G., Evidence of Bose-Einstein condensation in an atomic gas with attractive interaction, Phys. Rev. Lett., 75, 1687-1690 (1995)
[19] Briggs, W. L.; Henson, V. E.; McCormick, S. F., A Multigrid Tutorial (2000), SIAM: SIAM Philadelphia · Zbl 0958.65128
[20] Cancès, E.; Chakir, R.; He, L.; Maday, Y., Two-grid methods for a class of nonlinear elliptic eigenvalue problems, IMA J. Numer. Anal., hal-01088947 (2017)
[21] Chang, S. M.; Lin, W. W.; Shieh, S. F., Gauss-Seidel-type methods for energy states of a multi-component Bose-Einstein condensate, J. Comput. Phys., 202, 367-390 (2005) · Zbl 1056.81088
[22] Chiofalo, M. L.; Succi, S.; Tosi, M. P., Ground state of trapped interacting Bose-Einstein condensates by an explicit imaginary-time algorithm, Phys. Rev. E, 62, 7438-7444 (2000)
[23] Collin, A.; Massignan, P.; Pethick, C. J., Energy-dependent effective interactions for dilute many-body systems, Phys. Rev. A, 75, Article 013615 pp. (2007)
[24] Danaila, I.; Hecht, F., A finite element method with mesh adaptivity for computing vortex states in fast-rotating Bose-Einstein condensates, J. Comput. Phys., 229, 6946-6960 (2010) · Zbl 1198.82035
[25] Danaila, I.; Protas, B., Computation of Ground States of the Gross-Pitaevskii Functional via Riemannian Optimization, SIAM J. Sci. Comput., 39, B1102-B1129 (2017) · Zbl 1378.81174
[26] Davis, K. B.; Mewes, M. O.; Andrews, M. R.; van Druten, N. J.; Durfee, D. S.; Kurn, D. M.; Ketterle, W., Bose-Einstein condensation in a gas of sodium atoms, Phys. Rev. Lett., 75, 3969-3973 (1995)
[27] Esry, B. D.; Greene, C. H., Validity of the shape-independent approximation for Bose-Einstein condensates, Phys. Rev. A, 60, 1451-1462 (1999)
[28] Fu, H.; Wang, Y.; Gao, B., Beyond the Fermi pseudopotential: a modified Gross-Pitaevskii equation, Phys. Rev. A, 67, Article 053612 pp. (2003)
[29] Garcia-Ripoll, J. J.; Konotop, V. V.; Malomed, B. A.; Perez-Garcia, V. M., A quasi-local Gross-Pitaevskii equation for attractive Bose-Einstein condensates, Math. Comput. Simul., 62, 21-30 (2003) · Zbl 1034.35126
[30] Gross, E. P., Structure of a quantized vortex in boson systems, Nuovo Cimento, 20, 454-457 (1961) · Zbl 0100.42403
[31] Henning, P.; Målqvist, A.; Peterseim, D., Two-level discretization techniques for ground state computations Of Bose-Einstein condensates, SIAM J. Numer. Anal., 52, 1525-1550 (2014) · Zbl 1308.35272
[32] Lieb, E. H.; Seiringer, R.; Solovej, J. P.; Yngvason, J., The Mathematics of the Bose Gas and Its Condensation, Oberwolfach Seminars, vol. 34 (2005), Birkhäuser: Birkhäuser Basel · Zbl 1104.82012
[33] Lieb, E. H.; Seiringer, R.; Yngvason, J., Bosons in a trap: a rigorous derivation of the Gross-Pitaevskii energy functional, Phys. Rev. A, 61, Article 043602 pp. (2000)
[34] Muruganandam, P.; Adhikari, S. K., Bose-Einstein condensation dynamics in three dimensions by the pseudospectral and finite-difference methods, J. Phys. B, 36, 2501-2513 (2003)
[35] Pitaevskii, L. P., Vortex lines in an imperfect Bose gas, Sov. Phys. JETP, 13, 451-454 (1961)
[36] Pitaevskii, L. P.; Stringari, S., Bose-Einstein Condensation (2003), Clarendon Press: Clarendon Press Oxford · Zbl 1110.82002
[37] Ruan, X., Mathematical Theory and Numerical Methods for Bose-Einstein Condensation with Higher Order Interactions (2017), National University of Singapore, PhD Thesis
[38] Ruan, X.; Bao, W.; Cai, Y., Mean-field regime and Thomas-Fermi approximations of trapped Bose-Einstein condensates with higher order interactions in one and two dimensions, J. Phys. B, At. Mol. Opt. Phys., 49, Article 125304 pp. (2016)
[39] Shen, J.; Wang, C.; Wang, X.; Wise, S. M., Second-order convex splitting schemes for gradient flows with Ehrlich-Schwoebel type energy: application to thin film epitaxy, SIAM J. Numer. Anal., 50, 105-125 (2012) · Zbl 1247.65088
[40] Thøgersen, M.; Zinner, N. T.; Jensen, A. S., Thomas-Fermi approximation for a condensate with higher-order interactions, Phys. Rev. A, 80, Article 043625 pp. (2009)
[41] Veksler, H.; Fishman, S.; Ketterle, W., Simple model for interactions and corrections to the Gross-Pitaevskii equation, Phys. Rev. A, 90, Article 023620 pp. (2014)
[42] Wesseling, P., An Introduction to Multigrid Methods (1992), Wiley: Wiley Chichester · Zbl 0760.65092
[43] Wu, X.; Wen, Z.; Bao, W., A regularized Newton method for computing ground states of Bose-Einstein condensates, J. Sci. Comput., 73, 303-329 (2017) · Zbl 1433.82025
[44] Wu, X.; van Zwieten, G. J.; van der Zee, K. G., Stabilized second-order convex splitting schemes for Cahn-Hilliard models with application to diffuse-interface tumor-growth models, Int. J. Numer. Methods Biomed. Eng., 30, 2, 180-203 (2014)
[45] Xu, J.; Zhou, A., Local and parallel finite element algorithms based on two-grid discretizations, Math. Comput., 69, Article 881909 pp. (2000) · Zbl 0948.65122
[46] Yuille, A.; Rangarajan, A., The concave convex procedure, Neural Comput., 15, 915-936 (2003) · Zbl 1022.68112
[47] Zinner, N. T.; Thøgersen, M., Stability of a Bose-Einstein condensate with higher-order interactions near a Feshbach resonance, Phys. Rev. A, 80, Article 023607 pp. (2009)
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.