| Abstract: |
| We present a robust and efficient numerical method to compute the dynamics of the rotating two-component
dipolar Bose-Einstein condensates (BEC).
Using the rotating Lagrangian coordinates transform, we reformulate the original coupled Gross-Pitaevskii equations (CGPE) into new
equations where the rotating term vanishes and the potential becomes time-dependent. A time-splitting Fourier pseudospectral method is proposed to numerically solve the new equations
where the nonlocal Dipole-Dipole Interactions (DDI) are computed by a newly-developed Gaussian-sum (GauSum) solver which helps
achieve spectral accuracy in space within $O(N\log N)$ operations ($N$ is the total number of grid points).
The new method is spectrally accurate in space and second order accurate in time - these accuracies are confirmed numerically.
Dynamical properties of some physical quantities, including the total mass, energy, center of mass and angular momentum expectation, are presented and confirmed numerically. Interesting dynamical phenomena that are peculiar to the rotating two-component dipolar BECs, such as dynamics of center of mass,
quantized vortex lattices dynamics and the collapse dynamics in 3D, are presented. |
|