Time-dependent Gross-Pitaevskii Equation Research Articles

RKFD Scheme for Quantum Reflection Model of Bose-Einstein Condensates (BEC) from Silicon Surface

We applied the numerical combination of Runge-Kutta and Finite Difference (RKFD) scheme for a quantum reflection model of Bose-Einstein condensate (BEC) from a silicon surface. It is by the time-dependent Gross-Pitaevskii equation (GPE), a non-linear Schrödinger equation (NLSE) in the context of quantum mechanics. The role of cut-off potential δ and negative imaginary potential Vim is essential to estimating non-interacting BEC reflection models. Relying on these features, we performed a numerical simulation of the BEC quantum reflection model and calculated the effect of reflection probability R versus incident speed vx. The model is based on the three rapid potential variations: positive-step potential +Vstep, negative-step potential -Vstep, and Casimir-Polder potential VCP. As a result, the RKFP numerical scheme was successfully set up and applied to the quantum reflection model of BEC from the silicon surface. The numerical simulations results show that the reflection probability R decays exponentially to the incident speed vx.

Dynamic of time-independent and time-dependent asymmetric Gross-Pitaevskii equation around exceptional point

Systems operating at exceptional points (EPs) are highly sensitive to small perturbations, making it challenging to work near an EP. Eigenvalue analysis of the Gross-Pitaevskii equation has shown that asymmetric nonlinearity can compensate for detuning effects. However, an experimentally feasible system based on asymmetric nonlinear coupled resonators has not yet been explored. Additionally, some intriguing features of such a system are hidden in time domain analysis, which is rarely investigated. In this study, we demonstrate this feature using a full-wave simulation of an asymmetric nonlinear coupled resonator based on the finite-element method in Comsol. The time-dependent analysis reveals that detuning can shift the system from PT-symmetric to broken PT-symmetric (or vice versa), and nonlinearity can reverse this dynamic. This study provides an experimental framework for examining exceptional points (EPs) in nonlinear detuned coupled resonators and opens up new avenues for fundamental research into the influence of nonlinearity and detuning on the system’s state during EP encirclement.

Realization of photonic time crystals via spin bifurcations of exciton-polariton condensates

We have theoretically realized a photonic time crystal in a homogeneous spinor exciton-polariton condensate subject to a cavity strain induced energy splitting between the x- and y-polarized polaritons with distinct loss rates. Numerical modeling based on the time-dependent open dissipative Gross-Pitaevskii equation results in the phase transition from linearly-polarized spinor condensate at low pump powers into a circularly-polarized one and then a photonic time crystal at higher pump powers. This phenomenon occurs only under the condition that the lower-energy x-polarized condensate has a higher loss rate than the y-polaried condensate, which can exist intrinsically in the semiconductor microcavities.

Polygons of quantized vortices in Bose-Einstein condensates with a circular trap.

We consider a disk-shaped cold atom Bose-Einstein condensate with repulsive atom-atom interactions within a circular trap, described by a two-dimensional time-dependent Gross-Pitaevskii equationwith cubic nonlinearity and a circular box potential. In this setup, we discuss the existence of a type of stationary nonlinear waves with propagation-invariant density profiles, consisting of vortices located at the vertices of a regular polygon with or without an antivortex at its center. These polygons rotate around the center of the system and we provide approximate expressions for their angular velocity. For any size of the trap, we find a unique regular polygon solution that is static and is seemingly stable for long evolutions. It consists of a triangle of vortices with unit charge placed around a singly charged antivortex, with the size of the triangle fixed by the cancellation of competing effects on its rotation. There exist other geometries with discrete rotational symmetry that yield static solutions, even if they turn out to be unstable. By numerically integrating in real time the Gross-Pitaevskii equation, we compute the evolution of the vortex structures and discuss their stability and the fate of the instabilities that can unravel the regular polygon configurations. Such instabilities can be driven by the instability of the vortices themselves, by vortex-antivortex annihilation or by the eventual breaking of the symmetry due to the motion of the vortices.

Dynamic stability in spinor Bose gases in moiré lattices with square and hexagonal symmetries

A broad range of phenomena in correlated electrons traveling in moir\'e lattices has emerged in both scenarios, experiments and theory. In this paper we report the observation of a dynamic stability that arises in an analogous system to that of electrons, a weakly interacting spinor Bose gas of ultracold $^{23}\mathrm{Na}$ atoms lying in a single layer having a moir\'e pattern with square and hexagonal symmetries. Our paper is based on the dynamical description of two magnetic domains represented by two hyperfine spin components of a Bose condensate initially localized in the left and right halves of a moir\'e lattice defined by a specific angle $\ensuremath{\theta}$ plus a harmonic confinement. To demonstrate the persistence of such an initial condition under the competence of the moir\'e pattern and harmonic confinements we studied both single noninteracting and double-domain interacting cases. We solve the time dependent Gross-Pitaevskii equations, and track the time evolution of several observables on each half as a function of the twisting angle. In the case of square moir\'e lattices we found a dynamic stability for angles larger than a special one ${\ensuremath{\theta}}_{s}$, except for the Pythagorean angles. The value of such an angle depends on the existence of a harmonic trap when interactions are absent, while this dependence is negligible for the interacting case. Hexagonal moir\'e lattices exhibit the dynamic stability starting from a certain angle that also depends on the harmonic confinement for the noninteracting case.

OpenMP Fortran programs for solving the time-dependent dipolar Gross-Pitaevskii equation

OpenMP Fortran programs for solving the time-dependent dipolar Gross-Pitaevskii equation

Numerical solution to the time-dependent Gross-Pitaevskii equation

In this work we employ the split-step technique combined with a Legendre pseudospectral representation to solve various time-dependent GrossPitaevskii equations (GPE). Our findings based on the numerical accuracy of this approach applied for one-dimensional (1D) and two-dimensional (2D) problems show that it can provide accurate and stable solutions. Moreover, this approach has been applied to study the dynamics of the Bose-Einstein condensate which is modeled with the GPE. The breathing of condensate with a repulsive and attractive interactions trapped in 1D and 2D harmonic potentials has been simulated as well.

Poincaré index formula and analogy with the Kosterlitz-Thouless transition in a non-rotated cold atom Bose-Einstein condensate

A dilute gas of Bose-Einstein condensed atoms in a non-rotated and axially symmetric harmonic trap is modelled by the time dependent Gross-Pitaevskii equation. When the angular momentum carried by the condensate does not vanish, the minimum energy state describes vortices (or antivortices) that propagate around the trap center. The number of (anti)vortices increases with the angular momentum, and they repel each other to form Abrikosov lattices. Besides vortices and antivortices there are also stagnation points where the superflow vanishes; to our knowledge the stagnation points have not been analyzed previously, in the context of the Gross-Pitaevskii equation. The Poincaré index formula states that the difference in the number of vortices and stagnation points can never change. When the number of stagnation points is small, they tend to aggregate into degenerate propagating structures. But when the number becomes sufficiently large, the stagnation points tend to pair up with the vortex cores, to propagate around the trap center in regular lattice arrangements. There is an analogy with the geometry of the Kosterlitz-Thouless transition, with the angular momentum of the condensate as the external control parameter instead of the temperature.

Reduction-based strategy for optimal control of Bose-Einstein condensates.

Applications of Bose-Einstein condensates (BEC) often require that the condensate be prepared in a specific complex state. Optimal control is a reliable framework to prepare such a state while avoiding undesirable excitations, and, when applied to the time-dependent Gross-Pitaevskii equation(GPE) model of BEC in multiple space dimensions, results in a large computational problem. We propose a control method based on first reducing the problem, using a Galerkin expansion, from a partial differential equationto a low-dimensional Hamiltonian ordinary differential equationsystem. We then apply a two-stage hybrid control strategy. At the first stage, we approximate the control using a second Galerkin-like method known as the chopped random basis to derive a finite-dimensional nonlinear programing problem, which we solve with a differential evolution algorithm. This search method then yields a candidate local minimum which we further refine using a variant of gradient descent. This hybrid strategy allows us to greatly reduce excitations both in the reduced model and the full GPE system.

Moment of inertia and dynamical rotational response of a supersolid dipolar gas

We show that the knowledge of the time dependent response of a trapped gas, subject to a sudden rotation of a confining harmonic potential, allows for the determination of the moment of inertia of dipolar supersolid configurations. While in the presence of one-dimensional arrays of droplets the frequency of the resulting scissors oscillation provides accurate access to the value of the moment of inertia, two-dimensional like configurations are characterized by a multi-frequency structure in the rotating signal, reflecting the presence of significant rigid body components in the rotational motion. Using the formalism of response function theory and simulations based on the so-called extended time dependent Gross-Pitaevskii equation, we point out the crucial role played by the low frequency components in the determination of the moment of inertia and of its deviations from the irrotational value. We also propose a protocol based on the stationary rotation of the trap, followed by its sudden stop, which might provide a promising alternative to the experimental evaluation of the moment of inertia.

THE IMPACT OF INTERACTION ON BOSE-EINSTEIN CONDENSATION OF ULTRACOLD ATOMIC GASES IN TRAPPED POTENTIALS

In this paper, we discuss theoretically the behavior of ideal Bose-Einstein gases in a 3D box and in a harmonic trap. We make use of the mean filed theory based on the time dependent Gross-Pitaevskii equation which describes the main features of Bose-Einstein condensation. We investigate and discuss the density profiles and the ground state properties of a trapped condensate confined in harmonic trap. We found out that in an attractive interaction, the gas tends to increase its central density in order to lower the interaction energy.

Dynamics of Bose-Einstein condensates under anharmonic trap

The dynamics of weakly interacting three-dimensional Bose-Einstein condensates (BECs), trapped in external axially symmetric plus anharmonic distortion potential are studied. Within a variational approach and time-dependent Gross-Pitaevskii equation, the coupled condensate width equations are derived. By modulating anharmonic distortion of the trapping potential, nonlinear features are studied numerically and illustrated analytically, such as mode coupling of oscillation modes, and resonances. Furthermore, the stability of attractive interaction BEC in both repulsive and attractive anharmonic distortion is examined. We demonstrate that a small repulsive and attractive anharmonic distortion is effective in reducing (extending) the condensate stability region since it decreases (increases) the critical number of atoms in the trapping potential.

Interaction-modulated tunneling dynamics of a mixture of Bose-Fermi superfluid

In this paper, we study the interaction-modulated tunneling dynamics of a Bose-Fermi superfluid mixture, where a Bose-Einstein condensate (BEC) with weak repulsive interaction is confined in a symmetric deep double-well potential and an equally populated two-component Fermi gas in a harmonic potential symmetrically is positioned in the center of the double-well potential. The tunneling between the two wells is modulated by fermions trapped in a harmonic potential. When the temperature is adequately low and the bosonic particle number is adequately large, we can employ the mean-field theory to describe the evolution of the BEC in the double-well potential through the time-dependent Gross-Pitaevskii equation. For the Fermi gas in the harmonic potential trap, we consider the case where the inter-fermion interaction is tuned on the deep Bose-Einstein condensate of the inter-fermion Feshbach resonance, where two fermions of spin-up and spin-down form a two-body bound state. Within the regime, the Fermi gas is well described by a condensate of these fermionic dimers, and hence can be simulated as well by a Gross-Pitaevskii equation of dimers. The inter-species interactions couple the dynamics of the two species, which results in interesting features in the tunneling oscillations. The dynamic equations of the BEC in the double-well potential is described by a two-mode approximation. Coupling it with time-dependent Gross-Pitaevskii equation of the harmonically potential trapped molecular BEC, we numerically investigate the dynamical evolution of the Boson-Fermi hybrid system under different initial conditions. It is found that the interaction among fermions in a harmonic potential leads to strong non-linearity in the oscillations of the bosons in the double-well potential and enriches the tunneling dynamics of the bosons. Especially, it strengthens macroscopic quantum self-trapping. And the macroscopic quantum self-trapping can be expressed in three forms: the phase tends to be negative and monotonically decreases with time, the phase evolves with time, and the phase tends to be positive and increases monotonically with time. This means that it is possible the tunneling dynamics of the BEC in double-well potential is adjustable. Our results can be verified experimentally in a Bose-Fermi superfluid mixture by varying different interaction parameters via Feshbach resonance and confinement-induced resonance.

Stability of a Bose-condensed mixture on a bubble trap

Stability and dynamical behavior of binary Bose-Einstein condensed mixtures trapped on the surface of a rigid spherical shell are investigated in the mean-field level, exploring the miscibility with and without vortex charges, considering repulsive and attractive interactions. In order to compute the critical points for the stability, we follow the Bogoliubov-de Gennes method for the analysis of perturbed solutions, with the constraint that initially the stationary states are in a complete miscible configuration. For the perturbed equal density mixture, of a homogeneous uniform gas and when hidden vorticity is verified, with the species having opposite azimuthal circulation, we consider small perturbation analysis for each unstable mode, providing a complete diagram with the intra- and inter-species interaction role on the stability of the miscible system. Finally, beyond small perturbation analysis, we explore the dynamics of some repulsive and attractive inter-species states by full numerical solutions of the time-dependent Gross-Pitaevskii equation.

Characterization of the energy level-structure of a trapped dipolar Bose gas via mean-field parametric resonances

Analyzing the energy levels of a trapped Bose–Einstein condensate (BEC) can be difficult when dipole-dipole interactions (DDIs) are present. To address this issue, this study focuses on the parametric resonances (PRs) in the mean-field dynamics of a one-dimensional dipolar BEC (DBEC) over widely varying trapping geometries, with the primary objective of characterizing the energy levels of this system via analytical methods. This is achieved by matching the PR energies to the energy levels of the confining trap using perturbative methods. Further, this research reveals the role of the interplay between DDIs and the trapping geometry in defining the energies and amplitudes of the PRs. The PRs are induced by a negative Gaussian potential with a depth that oscillates with respect to time; DDIs also play a role in this induction. The dynamics of this system are modeled using the time-dependent Gross–Pitaevskii equation (TDGPE), which is numerically solved via the Crank–Nicolson method. The PRs are discussed based on analytical methods. First, we show that PRs similar to the ones obtained from the TDGPE can be reproduced via the Lagrangian variational method. Second, the energies at which the PRs occur are closely matched with the energy levels of the corresponding trap, calculated using the time-independent perturbation theory. Third, the most probable transitions between the trap energy levels yielding PRs are determined based on the time-dependent perturbation theory. The primary contribution of this research is that the energy levels of a DBEC within a complex trapping potential could be characterized.

Vortex precession and exchange in a Bose-Einstein condensate

Vortices in a Bose-Einstein condensate are modelled as spontaneously symmetry breaking minimum energy solutions of the time dependent Gross-Pitaevskii equation, using the method of constrained optimization. In a non-rotating axially symmetric trap, the core of a single vortex precesses around the trap center and, at the same time, the phase of its wave function shifts at a constant rate. The precession velocity, the speed of phase shift, and the distance between the vortex core and the trap center, depend continuously on the value of the conserved angular momentum that is carried by the entire condensate. In the case of a symmetric pair of identical vortices, the precession engages an emergent gauge field in their relative coordinate, with a flux that is equal to the ratio between the precession and shift velocities.

OpenMP solver for rotating spin-1 spin–orbit- and Rabi-coupled Bose–Einstein condensates

We present OpenMP version of a Fortran program for solving the Gross–Pitaevskii equation for a harmonically trapped three-component rotating spin-1 spinor Bose–Einstein condensate (BEC) in two spatial dimensions with or without spin–orbit (SO) and Rabi couplings. The program uses either Rashba or Dresselhaus SO coupling. We use the split-step Crank–Nicolson discretization scheme for imaginary- and real-time propagation to calculate stationary states and BEC dynamics, respectively. New version program summaryProgram title: BEC-GP-SPINOR-ROT-OMP, a program package containing programs spin-SO-rot-imre2d-omp.f90, with util.f90.CPC Library link to program files:https://doi.org/10.17632/j3wr4wn946.2Licensing provisions: Apache License 2.0Programming language: OpenMP Fortran 90/95. The program is tested with the GNU, Intel, PGI, and Oracle (former Sun) compilers.Supplementary material: File Supp.pdf gives additional details about the new program version and the underlying physical system.Journal Reference of previous version: Comput. Phys. Commun. 259 (2021) 107657.Does the new version supersede the previous version?: Only partially. The program spin-SO-rot-imre2d-omp.f90 supersedes spin-SO-imre2d-omp.f90, while the one-dimensional program is not part of this package.Nature of problem: The present Open Multi-Processing (OpenMP) Fortran program solves the time-dependent nonlinear partial differential Gross–Pitaevskii (GP) equation for a trapped rotating spinor Bose–Einstein condensate (BEC) in two spatial dimensions.Solution method: We employ the split-step Crank–Nicolson scheme to discretize the time-dependent GP equation in space and time. The discretized equation is then solved by imaginary- or real-time propagation, employing adequately small space and time steps, to yield the solution of stationary and non-stationary problems, respectively.Reason for new version: The BEC is a special form of matter called superfluid. A hallmark of superfluidity is the formation of quantized vortices in a rotating BEC. The present program can be used to study the generation of quantized vortices in a rotating spin-1 trapped BEC and hence should be of general interest to researchers from various fields.Summary of revisions: Previously we published Fortran [1] and C [2] programs for solving the mean-field GP equation for a BEC, which are now enjoying widespread use. Later we extended these programs to the more complex scenario of dipolar BECs [3], spin-1 spinor BECs [4], and of rotating BECs [5]. The OpenMP [6, 7] and CUDA/MPI [8, 9, 10] versions of these programs, designed to make these faster and more efficient in multi-core computers, are also available. In this paper we present Fortran 90/95 program for solving the GP equation of a two-dimensional (2D) rotating spin-1 spinor BEC with Rashba [11] and Dresselhaus [12] spin–orbit (SO) coupling and Rabi coupling, involving a modification over the same for a spin-1 spinor BEC [4]. A new input parameter OMEG, which represents the angular velocity of rotation Ω of the spin-1 spinor BEC, has been introduced in the program, following Ref. [5]. Besides this new parameter, the execution of the present program follows the same procedure as the 2D program of Ref. [4]. All other input parameters in the two programs are identical and the reader is advised to consult that reference for further details. For some values of input parameters the quantized vortices of a rotating BEC could be arranged in the form of a lattice with a certain spatial symmetry, e.g., triangular or square lattice [5]. In our numerical study, we established recently such a symmetric lattice structure for a Rashba SO-coupled rotating spin-1 BEC in the simplest case, without the Rabi coupling [13]. A Dresselhaus SO-coupled rotating spin-1 BEC should also lead to identical structure, provided the sign of the angular velocity of rotation is changed. For the sake of completeness, in the supplementary material related to this article that can be found online at URL we provide the corresponding GP equations for a rotating spin-1 BEC with some instructive numerical examples. The program package BEC-GP-SPINOR-ROT-OMP contains the programs spin-SO-rot-imre2d-omp.f90 and util.f90 in the directory src, as well as the files makefile and README.md. The makefile allows automated compilation of the program using different supported compilers (GNU, Intel, PGI, Oracle) by a simple make command, as in Ref. [4]. The file README.md contains instructions on how to compile and run the programs. The directory output contains examples of matching outputs of imaginary- and real-time propagation programs in sub-directories with a generic name rotxgamyferro or rotxgamyantiferro, where x denotes the value of the angular velocity of rotation Ω and y denotes the strength of the SO coupling γ for ferromagnetic (c0=482,c2=15) and antiferromagnetic (c0=669,c2=−3.1) cases. The results in imaginary-time sub-directories rot.3gam.5ferro and rot.3gam.5antiferro are calculated using the respective converged imaginary-time wave functions with zero angular velocity. The real-time sub-directories rot.3gam.5ferro and rot.3gam.5polar contain real-time results calculated using the respective converged imaginary-time wave functions as inputs. These sub-directories also contain gnuplot programs fig*.gnu which can be used to generate fig*.eps figure files of component densities.

Collective excitations of a spherical ultradilute quantum droplet

In three dimensions, exotic new state of matter of self-bound ultradilute quantum droplets can be realized in free space, when the mean-field attraction (i.e., with mean-field energy $E_{\textrm{MF}}\propto-n^{2}$ at the density $n$) is balanced by the repulsive beyond-mean-field quantum fluctuations (i.e., $E_{\textrm{BMF}}\propto n^{2+\gamma}$). The parameter $\gamma>0$ typically takes the value $1/2$ if we consider the Lee-Huang-Yang (LHY) energy functional, but it can vary when the beyond-LHY-effect becomes important or the three-body interaction becomes dominant. Here, we theoretically investigate how collective excitations of a three-dimensional quantum droplet are affected by the parameter $\gamma$ and a weak harmonic trapping potential, both of which could be tuned in experiments. We use both the approximate approach based on a Gaussian variational ansatz and the exact numerical solution of the Bogoliubov equations resulting from the linearized time-dependent extended Gross-Pitaevskii equation. We show that one of the key features of quantum droplets, i.e., the existence of the surface modes with dispersion relation $\omega_{s}\propto k^{3/2}$ is very robust with respect to the changes either in the parameter $\gamma$ or in the harmonic trapping potential. We predict the excitation spectrum of the droplet realized by binary $^{39}$K mixtures under the typical experimental conditions, which might be readily measured in current cold-atom laboratories.

Spin-1 spin–orbit- and Rabi-coupled Bose–Einstein condensate solver

We present OpenMP versions of FORTRAN programs for solving the Gross–Pitaevskii equation for a harmonically trapped three-component spin-1 spinor Bose–Einstein condensate (BEC) in one (1D) and two (2D) spatial dimensions with or without spin–orbit (SO) and Rabi couplings. Several different forms of SO coupling are included in the programs. We use the split-step Crank–Nicolson discretization for imaginary- and real-time propagation to calculate stationary states and BEC dynamics, respectively. The imaginary-time propagation programs calculate the lowest-energy stationary state. The real-time propagation programs can be used to study the dynamics. The simulation input parameters are provided at the beginning of each program. The programs propagate the condensate wave function and calculate several relevant physical quantities. Outputs of the programs include the wave function, energy, root-mean-square sizes, different density profiles (linear density for the 1D program, linear and surface densities for the 2D program). The imaginary- or real-time propagation can start with an analytic wave function or a pre-calculated numerical wave function. The imaginary-time propagation usually starts with an analytic wave function, while the real-time propagation is often initiated with the previously calculated converged imaginary-time wave function. Program summaryProgram title: BEC-GP-SPINOR, consisting of: BEC-GP-SPINOR-OMP package, containing programs spin-SO-imre1d-omp.f90 and spin-SO-imre2d-omp.f90, with util.f90.CPC Library link to program files:https://doi.org/10.17632/j3wr4wn946.1Licensing provisions: Apache License 2.0Programming language: OpenMP FORTRAN. The FORTRAN programs are tested with the GNU, Intel, PGI, and Oracle compiler.Nature of problem: The present Open Multi-Processing (OpenMP) FORTRAN programs solve the time-dependent nonlinear partial differential Gross–Pitaevskii (GP) equation for a trapped spinor Bose–Einstein condensate, with or without spin–orbit coupling, in one and two spatial dimensions.Solution method: We employ the split-step Crank–Nicolson rule to discretize the time-dependent GP equation in space and time. The discretized equation is then solved by imaginary- or real-time propagation, employing adequately small space and time steps, to yield the solution of stationary and non-stationary problems, respectively.

Enhanced transport of spin-orbit-coupled Bose gases in disordered potentials.

Anderson localization is a single-particle localization phenomena in disordered media that is accompanied by an absence of diffusion. Spin-orbit coupling (SOC) describes an interaction between a particle's spin and its momentum that directly affects its energy dispersion, for example, creating dispersion relations with gaps and multiple local minima. We show theoretically that combining one-dimensional spin-orbit coupling with a transverse Zeeman field suppresses the effects of disorder, thereby increasing the localization length and conductivity. This increase results from a suppression of backscattering between states in the gap of the SOC dispersion relation. Here, we focus specifically on the interplay of disorder from an optical speckle potential and SOC generated by two-photon Raman processes in quasi-one-dimensional Bose-Einstein condensates. We first describe backscattering by using a Fermi golden rule approach, and then numerically confirm this picture by solving the time-dependent one-dimensional Gross-Pitaevskii equation for a weakly interacting Bose-Einstein condensate with SOC and disorder. We find that on the tens of millisecond timescale of typical cold atom experiments moving in harmonic traps, initial states with momentum in the zero-momentum SOC gap evolve with negligible backscattering, while without SOC these same states rapidly localize.