A Simple and Efficient Numerical Method for Computing the Dynamics of Rotating Bose--Einstein Condensates via Rotating Lagrangian Coordinates
KAUST Grant NumberKUK-I1-007-43
MetadataShow full item record
AbstractWe propose a simple, efficient, and accurate numerical method for simulating the dynamics of rotating Bose-Einstein condensates (BECs) in a rotational frame with or without longrange dipole-dipole interaction (DDI). We begin with the three-dimensional (3D) Gross-Pitaevskii equation (GPE) with an angular momentum rotation term and/or long-range DDI, state the twodimensional (2D) GPE obtained from the 3D GPE via dimension reduction under anisotropic external potential, and review some dynamical laws related to the 2D and 3D GPEs. By introducing a rotating Lagrangian coordinate system, the original GPEs are reformulated to GPEs without the angular momentum rotation, which is replaced by a time-dependent potential in the new coordinate system. We then cast the conserved quantities and dynamical laws in the new rotating Lagrangian coordinates. Based on the new formulation of the GPE for rotating BECs in the rotating Lagrangian coordinates, a time-splitting spectral method is presented for computing the dynamics of rotating BECs. The new numerical method is explicit, simple to implement, unconditionally stable, and very efficient in computation. It is spectral-order accurate in space and second-order accurate in time and conserves the mass on the discrete level. We compare our method with some representative methods in the literature to demonstrate its efficiency and accuracy. In addition, the numerical method is applied to test the dynamical laws of rotating BECs such as the dynamics of condensate width, angular momentum expectation, and center of mass, and to investigate numerically the dynamics and interaction of quantized vortex lattices in rotating BECs without or with the long-range DDI.Copyright © by SIAM.
CitationBao W, Marahrens D, Tang Q, Zhang Y (2013) A Simple and Efficient Numerical Method for Computing the Dynamics of Rotating Bose--Einstein Condensates via Rotating Lagrangian Coordinates. SIAM Journal on Scientific Computing 35: A2671–A2695. Available: http://dx.doi.org/10.1137/130911111.
SponsorsThis research was supported by the Singapore A*STAR SERC "Complex Systems" Research Programme grant 1224504056.This author's work was supported by research award KUK-I1-007-43 from King Abdullah University of Science and Technology (KAUST).This author's work was supported by the University of Missouri Research Board, Simons Foundation award 210138, and NSF grant DMS-1217000.
CollectionsPublications Acknowledging KAUST Support
Showing items related by title, author, creator and subject.
Enhancing Rotational Diffusion Using Oscillatory ShearLeahy, Brian D.; Cheng, Xiang; Ong, Desmond C.; Liddell-Watson, Chekesha; Cohen, Itai (American Physical Society (APS), 2013-05-29)Taylor dispersion - shear-induced enhancement of translational diffusion - is an important phenomenon with applications ranging from pharmacology to geology. Through experiments and simulations, we show that rotational diffusion is also enhanced for anisotropic particles in oscillatory shear. This enhancement arises from variations in the particle's rotation (Jeffery orbit) and depends on the strain amplitude, rate, and particle aspect ratio in a manner that is distinct from the translational diffusion. This separate tunability of translational and rotational diffusion opens the door to new techniques for controlling positions and orientations of suspended anisotropic colloids. © 2013 American Physical Society.
High‐order rotated staggered finite difference modeling of 3D elastic wave propagation in general anisotropic mediaChu, Chunlei; Stoffa, Paul L.; Seif, Roustam (Society of Exploration Geophysicists, 2009-01)We analyze the dispersion properties and stability conditions of the high‐order convolutional finite difference operators and compare them with the conventional finite difference schemes. We observe that the convolutional finite difference method has better dispersion properties and becomes more efficient than the conventional finite difference method with the increasing order of accuracy. This makes the high‐order convolutional operator a good choice for anisotropic elastic wave simulations on rotated staggered grids since its enhanced dispersion properties can help to suppress the numerical dispersion error that is inherent in the rotated staggered grid structure and its efficiency can help us tackle 3D problems cost‐effectively.
Influence of strain gradients on lattice rotation in nano-indentation experiments: A numerical studyDemiral, Murat; Roy, Anish; El Sayed, Tamer S.; Silberschmidt, Vadim V. (Elsevier BV, 2014-07)In this paper the texture evolution in nano-indentation experiments was investigated numerically. To achieve this, a three-dimensional implicit finite-element model incorporating a strain-gradient crystal-plasticity theory was developed to represent accurately the deformation of a body-centred cubic metallic material. A hardening model was implemented to account for strain hardening of the involved slip systems. The surface topography around indents in different crystallographic orientations was compared to corresponding lattice rotations. The influence of strain gradients on the prediction of lattice rotations in nano-indentation was critically assessed. © 2014 Elsevier B.V..