Two-dimensional Bose-Einstein Condensate Research Articles

Stable-fixed-point description of square-pattern formation in driven two-dimensional Bose-Einstein condensates

Stable-fixed-point description of square-pattern formation in driven two-dimensional Bose-Einstein condensates

Superfluid fraction tensor of a two-dimensional supersolid

We investigate the superfluid fraction of crystalline stationary states within the framework of mean-field Gross–Pitaevskii theory. Our primary focus is on a two-dimensional Bose–Einstein condensate with a non-local soft-core interaction, where the superfluid fraction is described by a rank-2 tensor. We then calculate the superfluid fraction tensor for crystalline states exhibiting triangular, square, and stripe geometries across a broad range of interaction parameters. Factors leading to an anisotropic superfluid fraction tensor are also considered. We also refine the Leggett bounds for the superfluid fraction of the 2D system. We systematically compare these bounds to our full numerical results, and other results in the literature. This work is of direct relevance to other supersolid systems of current interest, such as supersolids produced using dipolar Bose–Einstein condensates.

Dissipation-driven superradiant phase transition of a two-dimensional BEC in a double-cavity

Abstract We study superradiant phase transitions in a hybrid system of a two-dimensional Bose-Einstein condensate of atoms and two cavities arranged with a tilt angle. By adjusting the loss rate of cavities, we map out the phase diagram of steady states within a mean field framework. We find that when the loss rates of the two cavities are different, superradiant transitions may not occur at the same time in the two cavities. A first-order phase transition is observed between the states with only one cavity in superradiance and both in superradiance. In the case of both cavities are superradiant, a net photon current is observed flowing from the cavity with small decay rate to the one with large decay rate. The photon current shows a non-monotonic dependence on the loss rate difference, owing to the competition of photon number difference and cavity field phase difference. Our findings can be realized and detected in experiments.

Linear limit continuation: Theory and an application to two-dimensional Bose–Einstein condensates

We present a coherent and effective theoretical framework to systematically construct numerically exact nonlinear solitary waves from their respective linear limits. First, all possible linear degenerate sets are classified for a harmonic potential using lattice planes. For a generic linear degenerate set, distinct wave patterns are identified in the near-linear regime using a random searching algorithm by suitably mixing the linear degenerate states, followed by a numerical continuation in the chemical potential extending the waves into the Thomas–Fermi regime. The method is applied to the two-dimensional, one-component Bose–Einstein condensates, yielding a spectacular set of waveforms. Our method opens a remarkably large program, and many more solitary waves are expected. Finally, the method can be readily generalized to three dimensions, and also multi-component condensates, providing a highly powerful technique for investigating solitary waves in future works.

U-net based vortex detection in Bose–Einstein condensates with automatic correction for manually mislabeled data

Quantum vortices in Bose–Einstein condensates (BECs) are essential phenomena in condensed matter physics, and precisely locating their positions, especially the vortex core, is a precondition for studying their properties. With the rise of machine learning, there is a possibility to expedite the localization process and provide accurate predictions. However, traditional machine learning requires particular considerable amount of manual data annotation, leading to uncontrollable accuracy. In this paper, we utilize the U-Net method to detect vortex positions accurately at the pixel level and propose an Automatic Correction Labeling (ACL) approach to optimize the acquisition of data sets for vortex localization in BECs. This approach addresses inaccuracies in the labeled vortex positions and improves the accuracy of vortex localization, especially the vortex core positions, while enhancing the tolerance for human mislabeling. The main process involves Rough Labeling →\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\rightarrow $$\\end{document} Machine Learning →\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\rightarrow $$\\end{document} Probability Region Search →\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\rightarrow $$\\end{document} Data Relabeling →\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\rightarrow $$\\end{document} Machine Learning again. The objective of ACL is to secure more accurate labeled data for model retraining. Through vortex localization experiments conducted in a two-dimensional Bose-Einstein condensate, our results establish the following: 1. Even under conditions of biased and missing manual annotations, U-Net can still accurately locate vortex positions; 2. Vortices exhibit certain regularities, and training U-Net with a small number of samples yields excellent predictive consequences; 3. The machine learning vortex locator based on the ACL method effectively corrects errors in manually annotated data, significantly improving the model’s performance metrics, thus enhancing the precision and metrics of vortex localization. This substantial advancement in the application of machine learning in vortex localization provides an effective way for vortex dynamics localization. Furthermore, this method of obtaining more accurate positions of approximate human labels through machine learning offers new insights for machine learning in other types of image recognition problems.

Dynamical nonlinear excitations induced by interaction quench in a two-dimensional box-trapped Bose–Einstein condensate

Manipulating nonlinear excitations, including solitons and vortices, is an essential topic in quantum many-body physics. A new progress in this direction is a protocol proposed in [Phys. Rev. Res. 2 043256 (2020)] to produce dark solitons in a one-dimensional atomic Bose–Einstein condensate (BEC) by quenching inter-atomic interaction. Motivated by this work, we generalize the protocol to a two-dimensional BEC and investigate the generic scenario of its post-quench dynamics. For an isotropic disk trap with a hard-wall boundary, we find that successive inward-moving ring dark solitons (RDSs) can be induced from the edge, and the number of RDSs can be controlled by tuning the ratio of the after- and before-quench interaction strength across different critical values. The role of the quench played on the profiles of the density, phase, and sound velocity is also investigated. Due to the snake instability, the RDSs then become vortex–antivortex pairs with peculiar dynamics managed by the initial density and the after-quench interaction. By tuning the geometry of the box traps, demonstrated as polygonal ones, more subtle dynamics of solitons and vortices are enabled. Our proposed protocol and the discovered rich dynamical effects on nonlinear excitations can be realized in near future cold-atom experiments.

Sonic horizon formation for two-dimensional Bose–Einstein condensates with higher-order nonlinear interaction

Based on the Gross–Pitaevskii equation model that incorporates higher-order nonlinear interaction, we studied sonic black hole and sonic horizon formation involved in nonlinear dynamical evolution for Bose–Einstein condensates with higher-order nonlinear interaction. On the basis of the modified variational method and the scenario where the system starts dynamic evolution from ground state, we derived the typical system distribution width function, which is analytically formulated as periodic oscillation solution and monotonically damped variational oscillation solution under different parametric settings. We also calculated the criteria formula for sonic black hole horizon formation with regard to the two evolution modes: oscillation mode and monotonically decay mode, pictorially demonstrating the time interval of sonic horizon appearance. The theoretical results obtained here can be used to guide relevant experimental studies of sonic horizon and sonic black hole formation for Bose–Einstein condensates incorporating the higher-order nonlinear interaction effects.

Precursory patterns to vortex nucleation in stirred Bose-Einstein condensates

Within the Gross-Pitaevskii theory, we study precursory pattern formations to the nucleation of vortex dipoles in a two-dimensional Bose-Einstein condensate stirred by a Gaussian potential. We introduce a smooth superfluid vorticity field and its conservative current, which capture very well the gradual process of vortex nucleation as a mechanism of topological singularities acquiring smooth cores. This is characterized by the localization of the superfluid vorticity into core regions which harbor the phase slips. For more impenetrable obstacles, we find that there are additional phase slips that do not acquire cores, thus remaining pinned as ghost vortices to the potential. We show that the vortex kinematics is slaved to the superfluid vorticity current, which determines not only the onset of nucleation but also the shedding dynamics.

Interplay of Kelvin-Helmholtz and superradiant instabilities of an array of quantized vortices in a two-dimensional Bose-Einstein condensate

We investigate the various physical mechanisms that underlie the dynamical instability of a quantized vortex array at the interface between two counter-propagating superflows in a two-dimensional Bose–Einstein condensate. Instabilities of markedly different nature are found to dominate in different flow velocity regimes. For moderate velocities where the two flows are subsonic, the vortex lattice displays a quantized version of the hydrodynamic Kelvin–Helmholtz instability (KHI), with the vortices rolling up and co-rotating. For supersonic flow velocities, the oscillation involved in the KHI can resonantly couple to acoustic excitations propagating away in the bulk fluid on both sides. This makes the KHI rate to be effectively suppressed and other mechanisms to dominate: For finite and relatively small systems along the transverse direction, the instability involves a repeated superradiant scattering of sound waves off the vortex lattice; for transversally unbound systems, a radiative instability dominates, leading to the simultaneous growth of a localized wave along the vortex lattice and of acoustic excitations propagating away in the bulk. Finally, for slow velocities, where the KHI rate is intrinsically slow, another instability associated to the rigid lateral displacement of the vortex lattice due to the vicinity of the system’s boundary is found to dominate.

Minimum critical velocity of a Gaussian obstacle in a Bose-Einstein condensate

When a superfluid flows past an obstacle, quantized vortices can be created in the wake above a certain critical velocity. In the experiment by Kwon et al. [Phys. Rev. A 91, 053615 (2015)], the critical velocity ${v}_{c}$ was measured for atomic Bose-Einstein condensates (BECs) using a moving repulsive Gaussian potential, and ${v}_{c}$ was minimized when the potential height ${V}_{0}$ of the obstacle was close to the condensate chemical potential $\ensuremath{\mu}$. Here we numerically investigate the evolution of the critical vortex shedding in a two-dimensional BEC with increasing ${V}_{0}$ and show that the minimum ${v}_{c}$ at the critical strength ${V}_{0c}\ensuremath{\approx}\ensuremath{\mu}$ results from the local density reduction and vortex-pinning effect of the repulsive obstacle. The spatial distribution of the superflow around the moving obstacle just below ${v}_{c}$ is examined. The particle density at the tip of the obstacle decreases as ${V}_{0}$ increases to ${V}_{c0}$, and at the critical strength, a vortex dipole is suddenly formed and dragged by the moving obstacle, indicating the onset of vortex pinning. The minimum ${v}_{c}$ exhibits power-law scaling with the obstacle size $\ensuremath{\sigma}$ as ${v}_{c}\ensuremath{\sim}{\ensuremath{\sigma}}^{\ensuremath{-}\ensuremath{\gamma}}$ with $\ensuremath{\gamma}\ensuremath{\approx}1/2$.

Quantum field simulator for dynamics in curved spacetime.

In most cosmological models, rapid expansion of space marks the first moments of the Universe and leads to the amplification of quantum fluctuations1. The description of subsequent dynamics and related questions in cosmology requires an understanding of the quantum fields of the standard model and dark matter in curved spacetime. Even the reduced problem of a scalar quantum field in an explicitly time-dependent spacetime metric is a theoretical challenge2-5, and thus a quantum field simulator can lead to insights. Here we demonstrate such a quantum field simulator in a two-dimensional Bose-Einstein condensate with a configurable trap6,7 and adjustable interaction strength to implement this model system. We explicitly show the realization of spacetimes with positive and negative spatial curvature by wave-packet propagation and observe particle-pair production in controlled power-law expansion of space, using Sakharov oscillations to extract amplitude and phase information of the produced state. We find quantitative agreement with analytical predictions for different curvatures in time and space. This benchmarks and thereby establishes a quantum field simulatorof a new class. In the future, straightforward upgrades offer the possibility to enter unexplored regimes that give further insight into relativistic quantum field dynamics.

A modified point-vortex model for the vortex dynamics in a two-dimensional superfluid Bose gas

In this paper, we investigate the dynamics of the vortices in the homogeneous two-dimensional Bose–Einstein condensate by considering the condensate density variation. It is shown that the condensate density variation plays an important role when two vortices are close together. By the time-dependent Lagrangian functional method and the numerical simulation, we study how the interaction strength of the vortices varies with the distance between two vortices with the same or different topological charges, and generalize the results for the vortex pairs to the multi-vortex case. The numerical simulations demonstrate the validity of our theory.

Asymmetric vortex dynamics in two-dimensional Bose–Einstein condensate with harmonic trap potential

Based on the two-dimensional Gross–Pitaevskii equation model, we investigated the asymmetric vortex evolution of two-dimensional Bose–Einstein condensates in a harmonic potential trap with polar direction perturbation. This corresponds to the generation of an asymmetric vortex evolution mode under initial vortex light manipulation together with polar direction perturbation. Unlike most of the prior work on this topic that uses a pure numerical method for the system under study, we use an exact analytical method rather than numerical simulation to investigate the key features of the system evolution dynamics in this study. Based on the variational method, and for different system parameter settings including the strength of the harmonic oscillator potential and the nonlinear interaction, we derived two evolution modes, namely, the periodic evolution oscillation mode and the monotonic decay mode, and pictorially demonstrated the evolution patterns of the system. In addition, we investigated the scenario when the system is in the rotating state, which corresponds to the action of the quadratic centrifugal potential. Notably, we also identified the damping effects for the two modes with increasing angular velocity, whose upper limit value corresponds to the constant vortex pattern of the system. Our theoretical results can be used to guide the experimental investigation of asymmetric vortex evolution in two-dimensional Bose–Einstein condensates.

Modulational instability and quantum droplets in a two-dimensional Bose-Einstein condensate

Modulational instability of a uniform two-dimensional binary Bose-Einstein condensate (BEC) in the presence of quantum fluctuations is studied. The analysis is based on the coupled Gross-Pitaevskii equations. It is shown that quantum fluctuations can induce instability when the BEC density is below a threshold. The dependence of the growth rate of modulations on the BEC parameters is found. It is observed that an asymmetry of the interaction parameters and/or initial densities of the components typically decreases the growth rate. Further development of the instability results in the BEC breaking up into a set of quantum droplets. These droplets merge dynamically with each other, so the total number of droplets decreases rapidly. The rate of this decrease is evaluated numerically for different initial parameters.

Modulation instability of two-dimensional Bose-Einstein condensates with helicoidal and a mixture of Rashba-Dresselhaus spin-orbit couplings

The modulational instability (MI) process is exclusively studied in a two-component Bose-Einstein condensate (BEC) which includes Rashba-Dresselhaus (RD) spin-orbit (SO) and helicoidal SO couplings. A generalized set of two-dimensional (2D) Gross-Pitaevskii (GP) equations is derived. The tunability of the helicoidal gauge potential is exploited to address BECs dynamics in a square lattice. The MI growth rate is derived, and parametric analyses of MI show the dependence of the instability on interatomic interaction strengths, the RD SO coupling, and helicoidal SO coupling, which combines the gauge amplitude and the helicoidal gauge potential. Direct numerical simulations are carried out to confirm the analytical predictions. Trains of solitons are obtained, and their behaviors are debated when the RD SO parameters are varied under different combinations between the gauge amplitude and the helicoidal gauge potential. The latter gives a potential way to manipulate the trapping capacities of the proposed BEC model.

Detection of roton and phonon excitations in a spin-orbit-coupled Bose-Einstein condensate with a moving barrier

We propose to detect phonon and roton excitations in a two-dimensional Bose-Einstein condensate with Raman-induced spin-orbit coupling by perturbing the atomic cloud with a weak barrier. The two excitation modes can be observed by moving the barrier along different directions in appropriate parameter regimes. Phonon excitations are identified by the appearance of solitary waves, while roton excitations lead to distinctive spatial density modulations. We show that this method can also be used to determine the anisotropic critical velocities of superfluid.

Stable striped state in a rotating two-dimensional spin–orbit coupled spin-1/2 Bose–Einstein condensate

We consider an effective two-dimensional Bose–Einstein condensate with some spin–orbit coupling (SOC) and a rotation term in an external harmonic potential. We find the striped state, and analyze the effects of SOC, the external potential, and the rotation frequency/direction on the profile and the stability of the striped state. Without the rotation term, the two spinor components exhibit striped pattern, and the numbers of stripes in the two components are always an odd–even or an even–odd. With the increase of the SOC strength, the number of stripes in both components increases, while the difference of the striped numbers is always one. After adding the rotation term, the profiles of the spinor components change qualitatively, and the change regulation of the striped numbers differs, while the difference of the striped numbers is still one. In addition, we find that the rotation direction only makes the striped state of the two spinor components exchange each other, though the clockwise and counterclockwise rotation directions are inequivalent with the presence of SOC. Such regulation is different from the previous study. And the rotation frequency gives rise to the transition from the striped state to a mixture of the striped state and vortex state. Furthermore, we prove the stability of these states by the evolution and linear stability analysis.

Quench Dynamics of Bose–Einstein Condensates in Boxlike Traps

By quenching the interatomic interactions, we investigate the nonequilibrium dynamics of two-dimensional Bose–Einstein condensates in boxlike traps with power-law potential boundaries. We show that ring dark solitons can be excited during the quench dynamics for both concave and convex potentials. The quench’s modulation strength and the steepness of the boundary are two major factors influencing the system’s evolution. In terms of the number of ring dark solitons excited in the condensate, five dynamic regimes have been identified. The condensate undergoes damped radius oscillation in the absence of ring dark soliton excitations. When it comes to the appearance of ring dark solitons, their decay produces interesting structures. The excitation patterns for the concave potential show a nested structure of vortex-antivortex pairs. The dynamic excitation patterns for the convex potential, on the other hand, show richer structures with multiple transport behaviors.

Sound-ring radiation of expanding vortex clusters

We investigate wave-vortex interaction emerging from an expanding compact vortex cluster in a two-dimensional Bose-Einstein condensate. We adapt techniques developed for compact gravitational objects to derive the characteristic modes of the wave-vortex interaction perturbatively around an effective vortex flow field. We demonstrate the existence of orbits or sound-rings, in analogy to gravitational light-rings, and compute the characteristic spectrum for the out-of-equilibrium vortex cluster. The spectrum obtained from numerical simulations of a stochastic Gross-Pitaevskii equation exhibiting an expanding vortex cluster is in excellent agreement with analytical predictions. Our findings are relevant for 2d-quantum turbulence, the semi-classical limit around fluid flows, and rotating compact objects exhibiting discrete circulation.

Vortex chains in rotating two-dimensional Bose-Einstein condensate in a harmonic plus optical lattices potential

<sec>Bose-Einstein condensate (BEC) is essentially a macroscopic quantum effect with quantum volatility, macroscopic quantum coherence and artificial controllability. Owing to its unique controllability, it becomes a new ideal platform for quantum simulations and studies of interacting quantum systems.</sec><sec>In this paper, the generation of vortices and the formation of vortex chains, as well as characteristics of vortex chains in rotating two-dimensional BEC in a potential composed of harmonic potential and optical lattice are studied numerically. Firstly, the generation of vortices, the formation and distribution of vortex chains and the effects of different physical parameters on the vortex chains in two-dimensional BEC are investigated by using the multigrid preconditioned conjugate gradient method. Secondly, the evolution of the vortex chains with time is studied by using the time-splitting spectral method. The results show that the generation of vortices in BEC trapped in the compound potential corresponds to the minimum value of the potential. When the depth of the optical lattice increases to a certain value, vortex chains are formed in the BEC. With the further increase of the depth of the optical lattice, the vortex depth in the vortex chain in the BEC decreases continuously, and finally the vortex chain disappears completely. When the interaction strength between atoms increases, the distribution range of the condensate expands, and the number of vortices and the number of vortex chains in the condensate also increase. When the interaction strength between atoms increases to a certain value, the symmetry of the vortex chains is broken. As the rotation frequency of the condensate increases, the distribution range of the condensate expands, and the number of vortices and the number of vortex chains in the condensate also increase. When the rotation frequency is close to the external trapping potential frequency, the linear alignment of the vortex chains is disrupted. It is also found that there are three stages in the evolution of the vortex chains in the BEC: in the first stage, vortex chains rotate together with the condensate, and the original chain distribution keeps unchanged; in the second stage, the phenomenon of vortex space extrusion appears, and the vortex chain is destroyed; in the third stage, the phenomenon of vortex space expansion occurs, and finally the vortex chains disappear. The results above show that the depth of the optical lattice, the interaction strength between atoms, and the rotation frequency of the condensate have important effects on the vortices and vortex chains in the condensate. By adjusting these physical quantities, the number of vortices and the shape of vortex chains in the BEC can be effectively manipulated. This may provide some theoretical reference and guidance for future experiments and applications.</sec>