GPELab, a Matlab toolbox to solve Gross-Pitaevskii equations. I: Computation of stationary solutions. (English) Zbl 1348.35003

Summary: This paper presents GPELab (Gross-Pitaevskii Equation Laboratory), an advanced easy-to-use and flexible Matlab toolbox for numerically simulating many complex physics situations related to Bose-Einstein condensation. The model equation that GPELab solves is the Gross-Pitaevskii equation. The aim of this first part is to present the physical problems and the robust and accurate numerical schemes that are implemented for computing stationary solutions, to show a few computational examples and to explain how the basic GPELab functions work. Problems that can be solved include: 1d, 2d and 3d situations, general potentials, large classes of local and nonlocal nonlinearities, multi-components problems, and fast rotating gases. The toolbox is developed in such a way that other physics applications that require the numerical solution of general Schrödinger-type equations can be considered.


35-04 Software, source code, etc. for problems pertaining to partial differential equations
82-04 Software, source code, etc. for problems pertaining to statistical mechanics
35Q55 NLS equations (nonlinear Schrödinger equations)
82-08 Computational methods (statistical mechanics) (MSC2010)
Full Text: DOI


[1] Gross, E. P., Structure of a quantized vortex in boson systems, Il Nuovo Cimento 10, 20, 3, 454-477, (1961) · Zbl 0100.42403
[2] M. Lewin, P.T. Nam, N. Rougerie, Derivation of Hartree’s theory for generic mean-field Bose systems, 2013, arXiv preprint arXiv:1303.0981. · Zbl 1316.81095
[3] Lieb, E. H.; Seiringer, R., Derivation of the Gross-Pitaevskii equation for rotating Bose gases, Comm. Math. Phys., 264, 2, 505-537, (2006) · Zbl 1233.82004
[4] Pethick, C. J.; Smith, H., Bose-Einstein condensation in dilute gases, (2002), Cambridge University Press
[5] Pitaevskii, L. P., Vortex lines in an imperfect Bose gas, Soviet Phys. JETP-USSR, 13, 2, (1961)
[6] 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, 5221, 198-201, (1995)
[7] Davis, K. B.; Mewes, M.-O.; van 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, 22, 3969-3973, (1995)
[8] Fried, D. G.; Killian, T. C.; Willmann, L.; Landhuis, D.; Moss, S. C.; Kleppner, D.; Greytak, T. J., Bose-Einstein condensation of atomic hydrogen, Phys. Rev. Lett., 81, 3811-3814, (1998)
[9] Jackson, B.; McCann, J. F.; Adams, C. S., Vortex formation in dilute inhomogeneous Bose-Einstein condensates, Phys. Rev. Lett., 80, 3903-3906, (1998)
[10] Kasamatsu, K.; Tsubota, M.; Ueda, M., Giant hole and circular superflow in a fast rotating Bose-Einstein condensate, Phys. Rev. A, 66, 053606, (2002)
[11] Kuga, T.; Torii, Y.; Shiokawa, N.; Hirano, T.; Shimizu, Y.; Sasada, H., Novel optical trap of atoms with a doughnut beam, Phys. Rev. Lett., 78, 4713-4716, (1997)
[12] Wen, L.; Xiong, H.; Wu, B., Hidden vortices in a Bose-Einstein condensate in a rotating double-well potential, Phys. Rev. A, 82, 5, 053627, (2010)
[13] Bao, W.; Cai, Y.; Wang, H., Efficient numerical methods for computing ground states and dynamics of dipolar Bose-Einstein condensates, J. Comput. Phys., 229, 20, 7874-7892, (2010) · Zbl 1198.82036
[14] Giovanazzi, S.; Görlitz, A.; Pfau, T., Tuning the dipolar interaction in quantum gases, Phys. Rev. Lett., 89, 130401, (2002)
[15] Góral, K.; Rza̧źewski, K.; Pfau, T., Bose-Einstein condensation with magnetic dipole-dipole forces, Phys. Rev. A, 61, 051601, (2000)
[16] Griesmaier, A.; Werner, J.; Hensler, S.; Stuhler, J.; Pfau, T., Bose-Einstein condensation of chromium, Phys. Rev. Lett., 94, 160401, (2005)
[17] Lahaye, T.; Menotti, C.; Santos, L.; Lewenstein, M.; Pfau, T., The physics of dipolar bosonic quantum gases, Rep. Progr. Phys., 72, 12, 126401, (2009)
[18] Pedri, P.; Santos, L., Two-dimensional bright solitons in dipolar Bose-Einstein condensates, Phys. Rev. Lett., 95, 200404, (2005)
[19] Abo-Shaeer, J. R.; Raman, C.; Vogels, J. M.; Ketterle, W., Observation of vortex lattices in Bose-Einstein condensates, Science, 292, 5516, 476-479, (2001)
[20] Madison, K. W.; Chevy, F.; Bretin, V.; Dalibard, J., Stationary states of a rotating Bose-Einstein condensate: routes to vortex nucleation, Phys. Rev. Lett., 86, 20, 4443-4446, (2001)
[21] Madison, K. W.; Chevy, F.; Wohlleben, W.; Dalibard, J., Vortex formation in a stirred Bose-Einstein condensate, Phys. Rev. Lett., 84, 5, 806-809, (2000)
[22] Madison, K. W.; Chevy, F.; Wohlleben, W.; Dalibard, J., Vortices in a stirred Bose-Einstein condensate, Journal of Modern Optics, 47, 14-15, 2715-2723, (2000)
[23] Raman, C.; Abo-Shaeer, J. R.; Vogels, J. M.; Xu, K.; Ketterle, W., Vortex nucleation in a stirred Bose-Einstein condensate, Phys. Rev. Lett., 87, 21, 210402, (2001)
[24] Abdullaev, F. Kh.; Baizakov, B. B.; Konotop, V. V., Dynamics of a Bose-Einstein condensate in optical trap, (Nonlinearity and Disorder: Theory and Applications, NATO Science Series, vol. 45, (2001), Springer Netherlands), 69-78 · Zbl 1032.82032
[25] Abdullaev, F. Kh.; Bronski, J. C.; Galimzyanov, R. M., Dynamics of a trapped 2d Bose-Einstein condensate with periodically and randomly varying atomic scattering length, Physica D, 184, 1-4, 319-332, (2003) · Zbl 1030.82006
[26] Abdullaev, F. Kh.; Bronski, J. C; Papanicolaou, G., Soliton perturbations and the random Kepler problem, Physica D, 135, 3-4, 369-386, (2000) · Zbl 0936.35171
[27] Garnier, J.; Abdullaev, F. Kh.; Baizakov, B. B., Collapse of a Bose-Einstein condensate induced by fluctuations of the laser intensity, Phys. Rev. A, 69, 053607, (2004)
[28] Gehm, M. E.; O’Hara, K. M.; Savard, T. A.; Thomas, J. E., Dynamics of noise-induced heating in atom traps, Phys. Rev. A, 58, 3914-3921, (1998)
[29] Savard, T. A.; O’Hara, K. M.; Thomas, J. E., Laser-noise-induced heating in far-off resonance optical traps, Phys. Rev. A, 56, R1095-R1098, (1997)
[30] Aftalion, A.; Mason, P., Phase diagrams and Thomas-Fermi estimates for spin-orbit-coupled Bose-Einstein condensates under rotation, Phys. Rev. A, 88, 023610, (2013)
[31] Bao, W., Ground states and dynamics of multicomponent Bose-Einstein condensates, Multiscale Model. Simul., 2, 2, 210-236, (2004) · Zbl 1062.82034
[32] Bao, W.; Lim, F. Y., Computing ground states of spin-1 Bose-Einstein condensates by the normalized gradient flow, SIAM J. Sci. Comput., 30, 4, 1925-1948, (2008) · Zbl 1173.81028
[33] Bao, W.; Wang, H., A mass and magnetization conservative and energy-diminishing numerical method for computing ground state of spin-1 Bose-Einstein condensates, SIAM J. Numer. Anal., 45, 5, 2177-2200, (2007) · Zbl 1149.82018
[34] Kawaguchi, Y.; Ueda, M., Spinor Bose-Einstein condensates, Phys. Rep., (2012)
[35] Law, C. K.; Pu, H.; Bigelow, N. P., Quantum spins mixing in spinor Bose-Einstein condensates, Phys. Rev. Lett., 81, 5257-5261, (1998)
[36] Lim, F. Y.; Bao, W., Numerical methods for computing the ground state of spin-1 Bose-Einstein condensates in a uniform magnetic field, Phys. Rev. E, 78, 6, 066704, (2008)
[37] Miesner, H.-J.; Stamper-Kurn, D. M.; Stenger, J.; Inouye, S.; Chikkatur, A. P.; Ketterle, W., Observation of metastable states in spinor Bose-Einstein condensates, Phys. Rev. Lett., 82, 2228-2231, (1999)
[38] Myatt, C. J.; Burt, E. A.; Ghrist, R. W.; Cornell, E. A.; Wieman, C. E., Production of two overlapping Bose-Einstein condensates by sympathetic cooling, Phys. Rev. Lett., 78, 586-589, (1997)
[39] Sadler, L. E.; Higbie, J. M.; Leslie, S. R.; Vengalattore, M.; Stamper-Kurn, D. M., Spontaneous symmetry breaking in a quenched ferromagnetic spinor Bose-Einstein condensate, Nature, 443, 7109, 312-315, (2006)
[40] Bao, W.; Chern, I.-L.; Zhang, Y., Efficient numerical methods for computing ground states of spin-1 Bose-Einstein condensates based on their characterizations, J. Comput. Phys., 253, 0, 189-208, (2013) · Zbl 1349.82069
[41] Xu, X.-Q.; Han, J. H., Spin-orbit coupled Bose-Einstein condensate under rotation, Phys. Rev. Lett., 107, 200401, (2011)
[42] Antoine, X.; Duboscq, R., Modeling and computation of Bose-Einstein condensates: stationary states, nucleation, dynamics, stochasticity, (Lecture Notes in Mathematics, (2014), Springer), (in press)
[43] Bao, W.; Chern, I.-L.; Lim, F. Y., Efficient and spectrally accurate numerical methods for computing ground and first excited states in Bose-Einstein condensates, J. Comput. Phys., 219, 2, 836-854, (2006) · Zbl 1330.82031
[44] Zeng, R.; Zhang, Y., Efficiently computing vortex lattices in rapid rotating Bose-Einstein condensates, Computer Physics Communications, 180, 6, 854-860, (2009) · Zbl 1198.82007
[45] X. Antoine, R. Duboscq, GPELab, A Matlab Toolbox to solve Gross-Pitaevskii Equations II: dynamics and stochastics. In preparation. · Zbl 1344.82004
[46] Antoine, X.; Bao, W.; Besse, C., Computational methods for the dynamics of the nonlinear Schrödinger/Gross-Pitaevskii equations, Comput. Phys. Comm., 184, 12, 2621-2633, (2013) · Zbl 1344.35130
[47] Agrawal, G., Nonlinear fiber optics, (2000), Springer · Zbl 1024.78514
[48] Agrawal, G., Applications of nonlinear fiber optics. optics and photonics, (2001), Elsevier Science
[49] Fouque, J.-P., Wave propagation and time reversal in randomly layered media, vol. 56, (2007), Springer
[50] Marty, R., On a splitting scheme for the nonlinear schroedinger equation in a random medium, Commun. Math. Sci., 4, 4, 363-376, (2006)
[51] Tiwari, R. P.; Shukla, A., A basis-set based Fortran program to solve the Gross-Pitaevskii equation for dilute Bose gases in harmonic and anharmonic traps, Computer Physics Communications, 174, 12, 966-982, (2006) · Zbl 1196.81039
[52] Muruganandam, P.; Adhikari, S. K., Fortran programs for the time-dependent Gross-Pitaevskii equation in a fully anisotropic trap, Computer Physics Communications, 180, 10, 1888-1912, (2009) · Zbl 1353.35002
[53] Vudragović, D.; Vidanović, I.; Balaž, A.; Muruganandam, P.; Adhikari, S. K., C programs for solving the time-dependent Gross-Pitaevskii equation in a fully anisotropic trap, Computer Physics Communications, 183, 9, 2021-2025, (2012) · Zbl 1353.35003
[54] Caliari, M.; Rainer, S., GSGPEs: a Matlab code for computing the ground state of systems of Gross-Pitaevskii equations, Computer Physics Communications, 184, 3, 812-823, (2013) · Zbl 1302.35004
[55] Hohenester, U., OCTBEC: a Matlab toolbox for optimal quantum control of Bose-Einstein condensates, Computer Physics Communications, (2013)
[56] Madarassy, E. J.M.; Toth, V. T., Numerical simulation code for self-gravitating Bose-Einstein condensates, Computer Physics Communications, 184, 4, 1339-1343, (2013)
[57] Adhikari, S. K., Numerical solution of the two-dimensional Gross-Pitaevskii equation for trapped interacting atoms, Phys. Lett. A, 265, 1, 91-96, (2000)
[58] Amara, P.; Hsu, D.; Straub, J. E., Global energy minimum searches using an approximate solution of the imaginary time Schrödinger equation, J. Phys. Chem., 97, 25, 6715-6721, (1993)
[59] Bao, W.; Du, Q., Computing the ground state solution of Bose-Einstein condensates by a normalized gradient flow, SIAM J. Sci. Comput., 25, 5, 1674-1697, (2004) · Zbl 1061.82025
[60] Baye, D.; Sparenberg, J.-M., Resolution of the Gross-Pitaevskii equation with the imaginary-time method on a Lagrange mesh, Phys. Rev. E, 82, 5, 056701, (2010)
[61] Cerimele, M. M.; Chiofalo, M. L.; Pistella, F.; Succi, S.; Tosi, M. P., Numerical solution of the Gross-Pitaevskii equation using an explicit finite-difference scheme: an application to trapped Bose-Einstein condensates, Phys. Rev. E, 62, 1, 1382, (2000)
[62] 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, 5, 7438, (2000)
[63] Edwards, M.; Burnett, K., Numerical solution of the nonlinear Schrödinger equation for small samples of trapped neutral atoms, Phys. Rev. A, 51, 1382-1386, (1995)
[64] Gammal, A.; Frederico, T.; Tomio, L., Improved numerical approach for the time-independent Gross-Pitaevskii nonlinear Schrödinger equation, Phys. Rev. E, 60, 2, 2421, (1999)
[65] Saad, Y., Iterative methods for sparse linear systems, (2003), SIAM · Zbl 1002.65042
[66] Saad, Y.; Schultz, M. H., GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems, SIAM J. Sci. Stat. Comput., 7, 3, 856-869, (1986) · Zbl 0599.65018
[67] Van der Vorst, H. A., Bi-CGSTAB: a fast and smoothly converging variant of bi-CG for the solution of nonsymmetric linear systems, SIAM J. Sci. Stat. Comput., 13, 2, 631-644, (1992) · Zbl 0761.65023
[68] 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., 258C, 509-523, (2014) · Zbl 1349.82027
[69] Pitaevskii, L. P.; Stringari, S., Bose-Einstein condensation, vol. 116, (2003), Clarendon Press
[70] Bao, W.; Tang, W., Ground-state solution of Bose-Einstein condensate by directly minimizing the energy functional, J. Comput. Phys., 187, 1, 230-254, (2003) · Zbl 1028.82500
[71] Dalfovo, F.; Stringari, S., Bosons in anisotropic traps: ground state and vortices, Phys. Rev. A, 53, 2477-2485, (1996)
[72] Holland, M. J.; Jin, D. S.; Chiofalo, M. L.; Cooper, J., Emergence of interaction effects in Bose-Einstein condensation, Phys. Rev. Lett., 78, 3801-3805, (1997)
[73] Bao, W.; Du, Q.; Zhang, Y., Dynamics of rotating Bose-Einstein condensates and its efficient and accurate numerical computation, SIAM J. Appl. Math., 66, 3, 758-786, (2006) · Zbl 1141.35052
[74] Jackson, A. D.; Kavoulakis, G. M.; Pethick, C. J., Solitary waves in clouds of Bose-Einstein condensed atoms, Phys. Rev. A, 58, 2417-2422, (1998)
[75] Leboeuf, P.; Pavloff, N., Bose-Einstein beams: coherent propagation through a guide, Phys. Rev. A, 64, 033602, (2001)
[76] X. Antoine, R. Duboscq, GPELab User Guide.
[77] Bao, W.; Cai, Y., Mathematical theory and numerical methods for Bose-Einstein condensation, Kinet. Relat. Models, 6, 1, 1-135, (2013) · Zbl 1266.82009
[78] Bao, W.; Wang, H., An efficient and spectrally accurate numerical method for computing dynamics of rotating Bose-Einstein condensates, J. Comput. Phys., 217, 2, 612-626, (2006) · Zbl 1160.82343
[79] Bao, W.; Wang, H.; Markowich, P. A., Ground, symmetric and central vortex states in rotating Bose-Einstein condensates, Commun. Math. Sci., 3, 1, 57-88, (2005) · Zbl 1073.82004
[80] Dunjko, V.; Lorent, V.; Olshanii, M., Bosons in cigar-shaped traps: Thomas-Fermi regime, Tonks-Girardeau regime, and in between, Phys. Rev. Lett., 86, 5413-5416, (2001)
[81] Choi, D.-I.; Niu, Q., Bose-Einstein condensates in an optical lattice, Phys. Rev. Lett., 82, 2022-2025, (1999)
[82] Bao, W.; Cai, Y., Ground states of two-component Bose-Einstein condensates with an internal atomic Josephson junction, East Asian J. Appl. Math., 1, 49-81, (2011) · Zbl 1290.35236
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. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.