Communications in Computational Physics

Pseudo-Arclength Continuation Algorithms for Binary Rydberg-Dressed Bose-Einstein Condensates

Sirilak Sriburadet, Y.-S. Wang, C.‐S. Chien, Yin-Tzer Shih

Published Jan 1, 2016• 3 citations

Abstract

We study pseudo-arclength continuation methods for both Rydberg-dressed Bose-Einstein condensates (BEC), and binary Rydberg-dressed BEC which are governed by the Gross-Pitaevskii equations (GPEs). A divide-and-conquer technique is proposed for rescaling the range/ranges of nonlocal nonlinear term/terms, which gives enough information for choosing a proper stepsize. This guarantees that the solution curve we wish to trace can be precisely approximated. In addition, the ground state solution would successfully evolve from one peak to vortices when the affect of the rotating term is imposed. Moreover, parameter variables with different number of components are exploited in curve-tracing. The proposed methods have the advantage of tracing the ground state solution curve once to compute the contours for various values of the coefficients of the nonlocal nonlinear term/terms. Our numerical results are consistent with those published in the literatures.

DOI: 10.4208/cicp.151214.021015a

Related papers

Multi-parameter continuation and collocation methods for rotating multi-component Bose–Einstein condensates

S.-Y. Chen, Y.-S. Wang, B.‐W. Jeng +1 more

We describe multi-parameter continuation methods combined with spectral collocation methods for computing numerical solutions of rotating two-component Bose–Einstein condensates (BECs), which are governed by the Gross–Pitaevskii equations (GPEs). Various types of orthogonal polynomials are used as the basis functions for the trial function space. A novel multi-parameter/multiscale continuation algorithm is proposed for computing the solutions of the governing GPEs, where the chemical potential of each component and angular velocity are treated as the continuation parameters simultaneously. The proposed algorithm can effectively compute numerical solutions with abundant physical phenomena. Numerical results on rotating two-component BECs are reported.

International Journal of Computer MathematicsApr 25, 20142 citations

A Two-Parameter Continuation Method for Rotating Two-Component Bose-Einstein Condensates in Optical Lattices

Y.-S. Wang, B.‐W. Jeng, C.‐S. Chien

We study efficient spectral-collocation and continuation methods (SCCM) for rotating two-component Bose-Einstein condensates (BECs) and rotating two-component BECs in optical lattices, where the second kind Chebyshev polynomials are used as the basis functions for the trial function space. A novel two-parameter continuation algorithm is proposed for computing the ground state and first excited state solutions of the governing Gross-Pitaevskii equations (GPEs), where the classical tangent vector is split into two constraint conditions for the bordered linear systems. Numerical results on rotating two-component BECs and rotating two-component BECs in optical lattices are reported. The results on the former are consistent with the published numerical results.

Communications in Computational PhysicsJun 27, 201220 citations

Efficient Continuation Methods for Computing Ground States of Quasi-2D Rotating Dipolar Bose-Einstein Condensates

B.‐W. Jeng, P.-H. Su, M.-I. Char

We study efficient continuation methods for computing the ground state solution of quasi-2D rotating dipolar Bose-Einstein condensates (BECs). First, the highly accurate spectral collocation method is used to discretize the governing Gross-Pitaevskii equation (GPE). Then, we modify the two-level continuation scheme for 3D dipolar BECs described in Jeng et al. (2014) to develop a single-parameter continuation method for quasi-2D rotating dipolar BECs, where the chemical potential is treated as the continuation parameter. Further, by adding the ratio of dipolar interaction strength to contact interaction strength as the second continuation parameter, we propose an efficient two-parameter continuation method which can effectively show the change of the ground-state vortex structures as the dipolar interaction strength gradually increases. Moreover, we also study linear stability analysis for the GPE. Sample numerical results on quasi-2D rotating dipolar BECs are reported.

Advances in Mathematical PhysicsAug 4, 20222 citations

A numerical scheme for the ground state of rotating spin-1 Bose–Einstein condensates

Sirilak Sriburadet, Yin-Tzer Shih, B.‐W. Jeng +2 more

We study the existence of nontrivial solution branches of three-coupled Gross-Pitaevskii equations (CGPEs), which are used as the mathematical model for rotating spin-1 Bose-Einstein condensates (BEC). The Lyapunov-Schmidt reduction is exploited to test the branching of nontrivial solution curves from the trivial one in some neighborhoods of bifurcation points. A multilevel continuation method is proposed for computing the ground state solution of rotating spin-1 BEC. By properly choosing the constraint conditions associated with the components of the parameter variable, the proposed algorithm can effectively compute the ground states of spin-1 [Formula: see text] and [Formula: see text] under rapid rotation. Extensive numerical results demonstrate the efficiency of the proposed algorithm. In particular, the affect of the magnetization on the CGPEs is investigated.

Scientific ReportsNov 23, 20211 citation

Efficient spectral collocation and continuation methods for symmetry-breaking solutions of rotating Bose–Einstein condensates

S.-Y. Chen, Zhenyang Li, Y.-S. Wang +1 more

We describe an efficient spectral collocation method (SCM) for symmetry-breaking solutions of rotating Bose–Einstein condensates (BECs) which is governed by the Gross–Pitaevskii equation (GPE). The Lagrange interpolants using the Legendre–Gauss–Lobatto points are used as the basis functions for the trial function space. Some formulas for the derivatives of the basis functions are given so that the GPE can be efficiently computed. The SCMs are incorporated in the context of a predictor–corrector continuation algorithm for tracing primary and secondary solution branches of the GPE. Symmetry-breaking solutions are numerically presented for both rotating BECs, BECs in optical lattices, and two-component BECs in optical lattices. Our numerical results show that the numerical algorithm we propose in this paper outperforms the classical orthogonal Legendre polynomials.

International Journal of Computer MathematicsOct 9, 20121 citation

Mathematical Models and Numerical Methods for Bose-Einstein Condensation

Weizhu Bao

The achievement of Bose-Einstein condensation (BEC) in ultracold vapors of alkali atoms has given enormous impulse to the theoretical and experimental study of dilute atomic gases in condensed quantum states inside magnetic traps and optical lattices. This article offers a short survey on mathematical models and theories as well as numerical methods for BEC based on the mean field theory. We start with the Gross-Pitaevskii equation (GPE) in three dimensions (3D) for modeling one-component BEC of the weakly interacting bosons, scale it to obtain a three-parameter model and show how to reduce it to two dimensions (2D) and one dimension (1D) GPEs in certain limiting regimes. Mathematical theories and numerical methods for ground states and dynamics of BEC are provided. Extensions to GPE with an angular momentum rotation term for a rotating BEC, to GPE with long-range anisotropic dipole-dipole interaction for a dipolar BEC and to coupled GPEs for spin-orbit coupled BECs are discussed. Finally, some conclusions are drawn and future research perspectives are discussed.

arXiv:1403.3884Mar 16, 201410 citations