High-order numerical methods with mass and energy conservation for spin–orbit-coupled Bose–Einstein condensates. Issue 4 (3rd April 2021)

• Record Type: Journal Article
• Authors: Li, Xiang-Gui; Li, Shu-Cun
• Is Part Of: International journal of computer mathematics. Volume 98:Issue 4(2021)
• Page Start: 738
• Page End: 757
• Publication Date: 2021-04-03
• Subjects: Spin–orbit-coupled Bose–Einstein condensates -- coupled Gross–Pitaevskii equations -- operator-compensation method -- high-order conservative method -- dynamics; 65M06 -- 35L05 -- 81-08; Computers -- Periodicals; Numerical analysis -- Periodicals; Automation -- Periodicals; 004.0151
• DOI: 10.1080/00207160.2020.1777407
• Languages: English
• ISSNs: 0020-7160

Abstract

ABSTRACT: In this work, a high-order numerical method with mass and energy conservation obtained by the operator-compensation (OC) technique is proposed for the coupled Gross–Pitaevskii equations (CGPEs) modelling the spin–orbit-coupled Bose–Einstein condensates. The OC technique is utilized to seek high-order approximations to the first- and the second-order spatial derivatives in the CGPEs, after which a time-dependent semi-discrete system is obtained. Then we obtain a fully discrete scheme by applying the Crank–Nicolson method to discretize this system in time. We theoretically show that the semi-discrete system and the fully discrete scheme conserve both mass and energy on the discrete level. In addition, the stability analysis is discussed for the proposed method. Numerical experiments are given to test accuracy and verify conservative properties. The dynamics of the two-dimensional CGPEs and the effects of the parameters in the CGPEs are simulated.