Watching dark solitons decay into vortex rings in a Bose-Einstein condensate

Semivortex solitons and their excited states in spin-orbit-coupled binary bosonic condensates

It is known that two-dimensional two-component fundamental solitons of the semivortex (SV) type, with vorticities (s_{+},s_{-})=(0,1) in their components, are stable ground states (GSs) in the spin-orbit-coupled (SOC) binary Bose-Einstein condensate with the contact self-attraction acting in both components, in spite of the possibility of the critical collapse in the system. However, excited states (ESs) of the SV solitons, with the vorticity set (s_{+},s_{-})=(S_{+},S_{+}+1) and S_{+}=1,2,3,..., are unstable in the same system. We construct ESs of SV solitons in the SOC system with opposite signs of the self-interaction in the two components. The main finding is stability of the ES-SV solitons, with the extra vorticity (at least) up to S_{+}=6. The threshold value of the norm for the onset of the critical collapse, N_{thr}, in these excited states is higher than the commonly known critical value, N_{c}≈5.85, associated with the single-component Townes solitons, N_{thr} increasing with the growth of S_{+}. A velocity interval for stable motion of the GS-SV solitons is found too. The results suggest a solution for the challenging problem of the creation of stable vortex solitons with high topological charges.

Self-trapping of vortex crystals via competing nonlinearities

We investigate the existence of self-trapped nonlinear waves with multiple phase singularities. Working with the cubic-quintic nonlinear Schrödinger equation, we focus on configurations with an antivortex surrounded by a triangular arrangement of vortices within a hosting soliton. We find stationary patterns that can be interpreted as stable self-trapped vortex crystals, constituting the first example of a configuration of this sort with space-independent potentials. Their stability is linked to their norm, transitioning from unstable to stable as their size increases, with an intermediate region where the structure is marginally unstable, undergoing a remarkable and puzzling self-reconstruction during its evolution.

Negative radiation pressure in Bose-Einstein condensates

In two-component nonlinear Schrödinger equations, the force exerted by incident monochromatic plane waves on an embedded dark soliton and on dark-bright-type solitons is investigated, both perturbatively and by numerical simulations. When the incoming wave is nonvanishing only in the orthogonal component to that of the embedded dark soliton, its acceleration is in the opposite direction to that of the incoming wave. This somewhat surprising phenomenon can be attributed to the well-known negative effective mass of the dark soliton. When a dark-bright soliton, whose effective mass is also negative, is hit by an incoming wave nonvanishing in the component corresponding to the dark soliton, the direction of its acceleration coincides with that of the incoming wave. This implies that the net force acting on it is in the opposite direction to that of the incoming wave. This rather counterintuitive effect is a yet another manifestation of negative radiation pressure exerted by the incident wave, observed in other systems. When a dark-bright soliton interacts with an incoming wave in the component of the bright soliton, it accelerates in the opposite direction; hence the force is pushing it now. We expect that these remarkable effects, in particular the negative radiation pressure, can be experimentally verified in Bose-Einstein condensates.

Vortex rings in paraxial laser beams

Interference of a fundamental vortex-free Gaussian beam with a co-propagating plane wave leads to nucleation of a series of vortex rings in the planes transverse to the optical axis; the number of rings grows with vanishing amplitude of the plane wave. In contrast, such interference with a beam carrying on-axis vortex with winding number l results in the formation of |l| rings elongated and gently twisted in propagation direction. The twist handedness of the vortex lines is determined by the interplay between dynamic and geometric phases of the Gaussian beam and the twist angle grows with vanishing amplitude of the plane wave. In the counter-propagating geometry the vortex rings nucleate and twist with half-wavelength period dominated by the interference grating in propagation direction.

Weak-measurement-induced heating in Bose-Einstein condensates

Ultracold atoms are an ideal platform for understanding system-reservoir dynamics of many-body systems. Here, we study quantum back-action in atomic Bose-Einstein condensates, weakly interacting with a far-from resonant, i.e., dispersively interacting, probe laser beam. The light scattered by the atoms can be considered as a part of quantum measurement process, whereby the change in the system state derives from measurement back-action. We experimentally quantify the resulting back-action in terms of the deposited energy. We model the interaction of the system and environment with a generalized measurement process, leading to a Markovian reservoir. Further, we identify two systematic sources of heating and loss: a stray optical lattice and probe-induced light-assisted collisions (an intrinsic atomic process). The observed heating and loss rates are larger for blue detuning than for red detuning, where they are oscillatory functions of detuning with increased loss at molecular resonances and reduced loss between molecular resonances.

Dynamical instability of 3D stationary and traveling planar dark solitons

Here we revisit the topic of stationary and propagating solitonic excitations in self-repulsive three-dimensional (3D) Bose-Einstein condensates by quantitatively comparing theoretical analysis and associated numerical computations with our experimental results. Motivated by numerous experimental efforts, including our own herein, we use fully 3D numerical simulations to explore the existence, stability, and evolution dynamics of planar dark solitons. This also allows us to examine their instability-induced decay products including solitonic vortices and vortex rings. In the trapped case and with no adjustable parameters, our numerical findings are in correspondence with experimentally observed coherent structures. Without a longitudinal trap, we identify numerically exact traveling solutions and quantify how their transverse destabilization threshold changes as a function of the solitary wave speed.

Observation of vortices and vortex stripes in a dipolar condensate

Quantized vortices are a prototypical feature of superfluidity that have been observed in multiple quantum gas experiments. But the occurrence of vortices in dipolar quantum gases-a class of ultracold gases characterized by long-range anisotropic interactions-has not been reported yet. Here we exploit the anisotropic nature of the dipole-dipole interaction of a dysprosium Bose-Einstein condensate to induce angular symmetry breaking in an otherwise cylindrically symmetric pancake-shaped trap. Tilting the magnetic field towards the radial plane deforms the cloud into an ellipsoid, which is then set into rotation. At stirring frequencies approaching the radial trap frequency, we observe the generation of dynamically unstable surface excitations, which cause angular momentum to be pumped into the system through vortices. Under continuous rotation, the vortices arrange into a stripe configuration along the field, in close agreement with numerical simulations.

Emergence and Ordering of Polygonal Breathers in Polariton Condensates

We show that the simultaneous driving of a polariton condensate with both nonresonant and nth order resonant pump frequencies allows for a generic mechanism of breather formation. From this we construct for the second order resonance a family of exotic breathers with nontrivial discrete order of rotational symmetry. Finally, we demonstrate the spontaneous emergence of both crystalline and glassy orderings of lattices of polygonal breathers, depending on the degree of polygonal excitations at the lattice sites.

Supersolid gap soliton in a Bose-Einstein condensate and optical ring cavity coupling system

The system of a transversely pumped Bose-Einstein condensate (BEC) coupled to a lossy ring cavity can favor a supersolid steady state. Here we find the existence of supersolid gap soliton in such a driven-dissipative system. By numerically solving the mean-field atom-cavity field coupling equations, gap solitons of a few different families have been identified. Their dynamical properties, including stability, propagation, and soliton collision, are also studied. Due to the feedback atom-intracavity field interaction, these supersolid gap solitons show numerous new features compared with the usual BEC gap solitons in static optical lattices.

Contour dynamics of two-dimensional dark solitons

Equations for contour dynamics of trough-shaped dark solitons are obtained for the general form of the nonlinearity function. Their self-similar solution which describes the nonlinear stage of the bending instability of dark solitons is studied in detail.

Transverse instability and dynamics of nonlocal bright solitons

We study the transverse instability and dynamics of bright soliton stripes in two-dimensional nonlocal nonlinear media. Using a multiscale perturbation method, we derive analytically the first-order correction to the soliton shape, which features an exponential growth in time-a signature of the transverse instability. The soliton's characteristic timescale associated with its exponential growth is found to depend on the square root of the nonlocality parameter. This, in turn, highlights the nonlocality-induced suppression of the transverse instability. Our analytical predictions are corroborated by direct numerical simulations, with the analytical results being in good agreement with the numerical ones.

Localization and spin dynamics of spin-orbit-coupled Bose-Einstein condensates in deep optical lattices

We analytically and numerically discuss the dynamics of two pseudospin components Bose-Einstein condensates (BECs) with spin-orbit coupling (SOC) in deep optical lattices. Rich localized phenomena, such as breathers, solitons, self-trapping, and diffusion, are revealed and strongly depend on the strength of the atomic interaction, SOC, Raman detuning, and the spin polarization (i.e., the initial population difference of atoms between the two pseudospin components of BECs). The critical conditions for the transition of localized states are derived analytically. Based on the critical conditions, the detailed dynamical phase diagram describing the different dynamical regimes is derived. When the Raman detuning satisfies a critical condition, localized states with a fixed initial spin polarization can be observed. When the critical condition is not satisfied, we use two quenching methods, i.e., suddenly and linearly quenching Raman detuning from the soliton or breather state, to discuss the spin dynamics, phase transition, and wave packet dynamics by numerical simulation. The sudden quenching results in a damped oscillation of spin polarization and transforms the system to a new polarized state. Interestingly, the linear quenching of Raman detuning induces a controllable phase transition from an unpolarized phase to an expected polarized phase, while the soliton or breather dynamics is maintained.

Self-Bayesian aberration removal via constraints for ultracold atom microscopy

High-resolution imaging of ultracold atoms typically requires custom high numerical aperture (NA) optics, as is the case for quantum gas microscopy. These high NA objectives involve many optical elements, each of which contributes to loss and light scattering, making them unsuitable for quantum backaction limited "weak" measurements. We employ a low-cost high NA aspheric lens as an objective for a practical and economical-although aberrated-high-resolution microscope to image ⁸⁷Rb Bose-Einstein condensates. Here, we present a methodology for digitally eliminating the resulting aberrations that is applicable to a wide range of imaging strategies and requires no additional hardware. We recover nearly the full NA of our objective, thereby demonstrating a simple and powerful digital aberration correction method for achieving optimal microscopy of quantum objects. This reconstruction relies on a high-quality measure of our imaging system's even-order aberrations from density-density correlations measured with differing degrees of defocus. We demonstrate our aberration compensation technique using phase-contrast imaging, a dispersive imaging technique directly applicable to quantum backaction limited measurements. Furthermore, we show that our digital correction technique reduces the contribution of photon shot noise to density-density correlation measurements which would otherwise contaminate the desired quantum projection noise signal in weak measurements.

Deep-learning-based quantum vortex detection in atomic Bose–Einstein condensates

Quantum vortices naturally emerge in rotating Bose–Einstein condensates (BECs) and, similarly to their classical counterparts, allow the study of a range of interesting out-of-equilibrium phenomena, such as turbulence and chaos. However, the study of such phenomena requires the determination of the precise location of each vortex within a BEC, which becomes challenging when either only the density of the condensate is available or sources of noise are present, as is typically the case in experimental settings. Here, we introduce a machine-learning-based vortex detector motivated by state-of-the-art object detection methods that can accurately locate vortices in simulated BEC density images. Our model allows for robust and real-time detection in noisy and non-equilibrium configurations. Furthermore, the network can distinguish between vortices and anti-vortices if the phase profile of the condensate is also available. We anticipate that our vortex detector will be advantageous for both experimental and theoretical studies of the static and dynamic properties of vortex configurations in BECs.

Stability analysis on dark solitons in quasi-1D Bose-Einstein condensate with three-body interactions

The stability properties of dark solitons in quasi-one-dimensional Bose–Einstein condensate (BEC) loaded in a Jacobian elliptic sine potential with three-body interactions are investigated theoretically. The solitons are obtained by the Newton-Conjugate Gradient method. A stationary cubic-quintic nonlinear Schrödinger equation is derived to describe the profiles of solitons via the multi-scale technique. It is found that the three-body interaction has distinct effect on the stability properties of solitons. Especially, such a nonlinear system supports the so-called dark solitons (kink or bubble), which can be excited not only in the gap, but also in the band. The bubbles are always linearly and dynamically unstable, and they cannot be excited if the three-body interaction is absent. Both stable and unstable kinks, depending on the physical parameters, can be excited in the BEC system.

Berezinskii-Kosterlitz-Thouless phase induced by dissipating quasisolitons

We theoretically study the sound propagation in a two-dimensional weakly interacting uniform Bose gas. Using the classical fields approximation we analyze in detail the properties of density waves generated both in a weak and strong perturbation regimes. While in the former case density excitations can be described in terms of hydrodynamic or collisionless sound, the strong disturbance of the system results in a qualitatively different response. We identify observed structures as quasisolitons and uncover their internal complexity for strong perturbation case. For this regime quasisolitons break into vortex pairs as time progresses, eventually reaching an equilibrium state. We find this state, characterized by only fluctuating in time averaged number of pairs of opposite charge vortices and by appearance of a quasi-long-range order, as the Berezinskii-Kosterlitz-Thouless (BKT) phase.

Observation of Magnetic Solitons in Two-Component Bose-Einstein Condensates

We experimentally investigate the dynamics of spin solitary waves (magnetic solitons) in a harmonically trapped, binary superfluid mixture. We measure the in situ density of each pseudospin component and their relative local phase via an interferometric technique we developed and as such, fully characterize the magnetic solitons while they undergo oscillatory motion in the trap. Magnetic solitons exhibit nondispersive, dissipationless longtime dynamics. By imprinting multiple magnetic solitons in our ultracold gas sample, we engineer binary collisions between solitons of either the same or opposite magnetization and map out their trajectories.

Magnetic Solitons in a Spin-1 Bose-Einstein Condensate

Vector solitons are a type of solitary or nonspreading wave packet occurring in a nonlinear medium composed of multiple components. As such, a variety of synthetic systems can be constructed to explore their properties, from nonlinear optics to ultracold atoms, and even in metamaterials. Bose-Einstein condensates have a rich panoply of internal hyperfine levels, or spin components, which make them a unique platform for exploring these solitary waves. However, existing experimental work has focused largely on binary systems confined to the Manakov limit of the nonlinear equations governing the soliton behavior, where quantum magnetism plays no role. Here we observe, using a "magnetic shadowing" technique, a new type of soliton in a spinor Bose-Einstein condensate, one that exists only when the underlying interactions are antiferromagnetic and which is deeply embedded within a full spin-1 quantum system. Our approach opens up a vista for future studies of "solitonic matter" whereby multiple solitons interact with one another at deterministic locations.

Creating solitons with controllable and near-zero velocity in Bose-Einstein condensates

Established techniques for deterministically creating dark solitons in repulsively interacting atomic Bose-Einstein condensates (BECs) can only access a narrow range of soliton velocities. Because velocity affects the stability of individual solitons and the properties of soliton-soliton interactions, this technical limitation has hindered experimental progress. Here we create dark solitons in highly anisotropic cigar-shaped BECs with arbitrary position and velocity by simultaneously engineering the amplitude and phase of the condensate wave function, improving upon previous techniques which explicitly manipulated only the condensate phase. The single dark soliton solution present in true one-dimensional (1D) systems corresponds to the kink soliton in anisotropic three-dimensional systems and is joined by a host of additional dark solitons, including vortex ring and solitonic vortex solutions. We readily create dark solitons with speeds from zero to half the sound speed. The observed soliton oscillation frequency suggests that we imprinted solitonic vortices, which for our cigar-shaped system are the only stable solitons expected for these velocities. Our numerical simulations of 1D BECs show this technique to be equally effective for creating kink solitons when they are stable. We demonstrate the utility of this technique by deterministically colliding dark solitons with domain walls in two-component spinor BECs.

Solitons in spin-orbit-coupled spin-2 spinor Bose-Einstein condensates

We investigate the different types of matter-wave solitons in spin-orbit-coupled spin-2 spinor Bose-Einstein condensates. Using mean-field theory and adopting the multiscale perturbation method, the original five-component Gross-Pitaevskii spin-orbit-coupled spin-2 spinor Bose-Einstein condensate model can be reduced to a single effective nonlinear Schrödinger equation, which allows us to find analytical soliton solutions of this system. In this way, for different regimes of the spin-orbit coupling, Raman coupling, and interatomic interactions, we find approximate bright and dark soliton solutions. Particularly, the type of solitons depends on the dispersion properties of the system. When the lowest-energy band has a single-well structure, we find there only exist positive mass bright or dark solitons due to the dispersion coefficient of effective nonlinear Shrödinger equation always positive. However, when the lowest-energy band has a double-well structure, there will appear positive (negative) mass bright or dark solitons because the sign of the dispersion coefficient can be positive (negative) under different momentum. We employ direct numerical simulation of the original five-component Gross-Pitaevskii equations to confirm the analytical results.

Flemish Strings of Magnetic Solitons and a Nonthermal Fixed Point in a One-Dimensional Antiferromagnetic Spin-1 Bose Gas

Thermalization in a quenched one-dimensional antiferromagnetic spin-1 Bose gas is shown to proceed via a nonthermal fixed point through annihilation of Flemish-string bound states of magnetic solitons. A possible experimental situation is discussed.

Impact of higher-order nonlinearity on modulational instability in two-component Bose-Einstein condensates

We investigate the effect of higher-order interactions induced by shape-dependent confinement in the modulational instability (MI) of a binary mixture of Bose-Einstein condensates. For this, we present and compute both analytically and numerically a system of coupled Gross-Pitaevskii equations with residual nonlinearity that rule the dynamics of the mixture. Using the linear stability approach, we obtain the instability criteria of the mixture and find that the MI can be excited in miscible condensates and altered in immiscible condensates due to the effect of residual nonlinearity. Direct numerical calculations are performed to support the analytical predictions, and a good agreement is found. The space-time evolution of the condensate density is displayed in both cases when the mixture is miscible and immiscible, showing the generation of bright solitons for modes predicted to be unstable.

Observation of Wall-Vortex Composite Defects in a Spinor Bose-Einstein Condensate

We report the observation of spin domain walls bounded by half-quantum vortices (HQVs) in a spin-1 Bose-Einstein condensate with antiferromagnetic interactions. A spinor condensate is initially prepared in the easy-plane polar phase, and then, suddenly quenched into the easy-axis polar phase. Domain walls are created via the spontaneous Z_{2} symmetry breaking in the phase transition and the walls dynamically split into composite defects due to snake instability. The end points of the defects are identified as HQVs for the polar order parameter and the mass supercurrent in their proximity is demonstrated using Bragg scattering. In a strong quench regime, we observe that singly charged quantum vortices are formed with the relaxation of free wall-vortex composite defects. Our results demonstrate a nucleation mechanism for composite defects via phase transition dynamics.

Measurement-induced dynamics and stabilization of spinor-condensate domain walls

Weakly measuring many-body systems and allowing for feedback in real time can simultaneously create and measure new phenomena in quantum systems. We theoretically study the dynamics of a continuously measured two-component Bose-Einstein condensate (BEC) potentially containing a domain wall and focus on the tradeoff between usable information obtained from measurement and quantum backaction. Each weakly measured system yields a measurement record from which we extract real-time dynamics of the domain wall. We show that quantum backaction due to measurement causes two primary effects: domain-wall diffusion and overall heating. The system dynamics and signal-to-noise ratio depend on the choice of measurement observable. We propose a feedback protocol to dynamically create a stable domain wall in the regime where domain walls are unstable, giving a prototype example of Hamiltonian engineering using measurement and feedback.

Dissipative shock waves generated by a quantum-mechanical piston

The piston shock problem is a prototypical example of strongly nonlinear fluid flow that enables the experimental exploration of fluid dynamics in extreme regimes. Here we investigate this problem for a nominally dissipationless, superfluid Bose-Einstein condensate and observe rich dynamics including the formation of a plateau region, a non-expanding shock front, and rarefaction waves. Many aspects of the observed dynamics follow predictions of classical dissipative-rather than superfluid dispersive-shock theory. The emergence of dissipative-like dynamics is attributed to the decay of large amplitude excitations at the shock front into turbulent vortex excitations, which allow us to invoke an eddy viscosity hypothesis. Our experimental observations are accompanied by numerical simulations of the mean-field, Gross-Pitaevskii equation that exhibit quantitative agreement with no fitting parameters. This work provides an avenue for the investigation of quantum shock waves and turbulence in channel geometries, which are currently the focus of intense research efforts.

Management of matter-wave solitons in Bose-Einstein condensates with time-dependent atomic scattering length in a time-dependent parabolic complex potential

In this paper, we consider a Gross-Pitaevskii (GP) equation with a time-dependent nonlinearity and a spatiotemporal complex linear term which describes the dynamics of matter-wave solitons in Bose-Einstein condensates (BECs) with time-dependent interatomic interactions in a parabolic potential in the presence of feeding or loss of atoms. We establish the integrability conditions under which analytical solutions describing the modulational instability and the propagation of both bright and dark solitary waves on a continuous wave background are obtained. The obtained integrability conditions also appear as the conditions under which the solitary waves of the BECs can be managed by controlling the functional gain or loss parameter. For specific BECs, the dynamics of bright and dark solitons are investigated analytically through the found exact solutions of the GP equation. Our results show that under the integrability conditions, the gain or loss parameter of the GP equation can be used to manage the motion of both bright and dark solitons. We show that for BECs with loss (gain) of atoms, the bright and dark solitons during their propagation have a compression (broadening) in their width. Furthermore, under a safe range of parameters and under the integrability conditions, it is possible to squeeze a bright soliton of BECs with loss of atoms into the assumed peak matter density, which can provide an experimental tool for investigating the range of validity of the 1D GP equation. Our results also reveal that under the conditions of the solitary wave management, neither the injection or the ejection of atoms from the condensate affects the soliton peak during its propagation.

Suppressed Solitonic Cascade in Spin-Imbalanced Superfluid Fermi Gas

Cold atoms experiments offer invaluable information on superfluid dynamics, including decay cascades of topological defects. While the cascade properties are well established for Bose systems, our understanding of their behavior in Fermi counterparts is very limited, in particular in spin-imbalanced systems, where superfluid (paired) and normal (unpaired) particles naturally coexist giving rise to complex spatial structure of the atomic cloud. Here we show, based on a newly developed microscopic approach, that the decay cascades of topological defects are dramatically modified by the spin polarization. We demonstrate that decay cascades end up at different stages: "dark soliton," "vortex ring," or "vortex line," depending on the polarization. We reveal that it is caused by sucking of unpaired particles into the soliton's internal structure. As a consequence vortex reconnections are hindered and we anticipate that quantum turbulence phenomenon can be significantly affected, indicating new physics induced by polarization effects.

Three-Component Soliton States in Spinor F=1 Bose-Einstein Condensates

Dilute-gas Bose-Einstein condensates are an exceptionally versatile test bed for the investigation of novel solitonic structures. While matter-wave solitons in one- and two-component systems have been the focus of intense research efforts, an extension to three components has never been attempted in experiments. Here, we experimentally demonstrate the existence of robust dark-bright-bright (DBB) and dark-dark-bright solitons in a multicomponent F=1 condensate. We observe lifetimes on the order of hundreds of milliseconds for these structures. Our theoretical analysis, based on a multiscale expansion method, shows that small-amplitude solitons of these types obey universal long-short wave resonant interaction models, namely, Yajima-Oikawa systems. Our experimental and analytical findings are corroborated by direct numerical simulations highlighting the persistence of, e.g., the DBB soliton states, as well as their robust oscillations in the trap.

First Observation of Bright Solitons in Bulk Superfluid ^{4}He

The existence of bright solitons in bulk superfluid ^{4}He is demonstrated by time-resolved shadowgraph imaging experiments and density functional theory (DFT) calculations. The initial liquid compression that leads to the creation of nonlinear waves is produced by rapidly expanding plasma from laser ablation. After the leading dissipative period, these waves transform into bright solitons, which exhibit three characteristic features: dispersionless propagation, negligible interaction in a two-wave collision, and direct dependence between soliton amplitude and the propagation velocity. The experimental observations are supported by DFT calculations, which show rapid evolution of the initially compressed liquid into bright solitons. At high amplitudes, solitons become unstable and break down into dispersive shock waves.

Quantum phase transition in an effective three-mode model of interacting bosons

In this work we study an effective three-mode model describing interacting bosons. These bosons can be considered as exciton-polaritons in a semiconductor microcavity at the magic angle. This model exhibits quantum phase transition (QPT) when the parameters of the corresponding Hamiltonian are continuously varied. The properties of the Hamiltonian spectrum (e.g., the distance between two adjacent energy levels) and the phase space structure of the thermodynamic limit of the model are used to indicate QPT. The relation between spectral properties of the Hamiltonian and the corresponding classical frame of the thermodynamic limit of the model is established as indicative of QPT. The average number of bosons in a specific mode and the entanglement properties of the ground state as functions of the parameters are used to characterize the order of the transition and also to construct a phase diagram. Finally, we verify our results for experimental data obtained for a setting of exciton-polaritons in a semiconductor microcavity.

Critical Spin Superflow in a Spinor Bose-Einstein Condensate

We investigate the critical dynamics of spin superflow in an easy-plane antiferromagnetic spinor Bose-Einstein condensate. Spin-dipole oscillations are induced in a trapped condensate by applying a linear magnetic field gradient and we observe that the damping rate increases rapidly as the field gradient increases above a certain critical value. The onset of dissipation is found to be associated with the generation of dark-bright solitons due to the modulation instability of the counterflow of two spin components. Spin turbulence emerges as the solitons decay because of their snake instability. We identify another critical point for spin superflow, in which transverse magnon excitations are dynamically generated via spin-exchanging collisions, which leads to the transient formation of axial polar spin domains.

Observation of Two-Dimensional Localized Jones-Roberts Solitons in Bose-Einstein Condensates

Jones-Roberts solitons are the only known class of stable dark solitonic solutions of the nonlinear Schrödinger equation in two and three dimensions. They feature a distinctive elongated elliptical shape that allows them to travel without change of form. By imprinting a triangular phase pattern, we experimentally generate two-dimensional Jones-Roberts solitons in a three-dimensional atomic Bose-Einstein condensate. We monitor their dynamics, observing that this kind of soliton is indeed not affected by dynamic (snaking) or thermodynamic instabilities, that instead make other classes of dark solitons unstable in dimensions higher than one. Our results confirm the prediction that Jones-Roberts solitons are stable solutions of the nonlinear Schrödinger equation and promote them for applications beyond matter wave physics, like energy and information transport in noisy and inhomogeneous environments.

Collective Modes of a Soliton Train in a Fermi Superfluid

We characterize the collective modes of a soliton train in a quasi-one-dimensional Fermi superfluid, using a mean-field formalism. In addition to the expected Goldstone and Higgs modes, we find novel long-lived gapped modes associated with oscillations of the soliton cores. The soliton train has an instability that depends strongly on the interaction strength and the spacing of solitons. It can be stabilized by filling each soliton with an unpaired fermion, thus forming a commensurate Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) phase. We find that such a state is always dynamically stable, which paves the way for realizing long-lived FFLO states in experiments via phase imprinting.

Adiabatic Invariant Approach to Transverse Instability: Landau Dynamics of Soliton Filaments

Consider a lower-dimensional solitonic structure embedded in a higher-dimensional space, e.g., a 1D dark soliton embedded in 2D space, a ring dark soliton in 2D space, a spherical shell soliton in 3D space, etc. By extending the Landau dynamics approach [Phys. Rev. Lett. 93, 240403 (2004)PRLTAO0031-900710.1103/PhysRevLett.93.240403], we show that it is possible to capture the transverse dynamical modes (the "Kelvin modes") of the undulation of this "soliton filament" within the higher-dimensional space. These are the transverse stability or instability modes and are the ones potentially responsible for the breakup of the soliton into structures such as vortices, vortex rings, etc. We present the theory and case examples in 2D and 3D, corroborating the results by numerical stability and dynamical computations.

Vortices and vortex lattices in quantum ferrofluids

The experimental realization of quantum-degenerate Bose gases made of atoms with sizeable magnetic dipole moments has created a new type of fluid, known as a quantum ferrofluid, which combines the extraordinary properties of superfluidity and ferrofluidity. A hallmark of superfluids is that they are constrained to rotate through vortices with quantized circulation. In quantum ferrofluids the long-range dipolar interactions add new ingredients by inducing magnetostriction and instabilities, and also affect the structural properties of vortices and vortex lattices. Here we give a review of the theory of vortices in dipolar Bose-Einstein condensates, exploring the interplay of magnetism with vorticity and contrasting this with the established behaviour in non-dipolar condensates. We cover single vortex solutions, including structure, energy and stability, vortex pairs, including interactions and dynamics, and also vortex lattices. Our discussion is founded on the mean-field theory provided by the dipolar Gross-Pitaevskii equation, ranging from analytic treatments based on the Thomas-Fermi (hydrodynamic) and variational approaches to full numerical simulations. Routes for generating vortices in dipolar condensates are discussed, with particular attention paid to rotating condensates, where surface instabilities drive the nucleation of vortices, and lead to the emergence of rich and varied vortex lattice structures. We also present an outlook, including potential extensions to degenerate Fermi gases, quantum Hall physics, toroidal systems and the Berezinskii-Kosterlitz-Thouless transition.

Brownian motion of solitons in a Bose-Einstein condensate

We observed and controlled the Brownian motion of solitons. We launched solitonic excitations in highly elongated [Formula: see text] Bose-Einstein condensates (BECs) and showed that a dilute background of impurity atoms in a different internal state dramatically affects the soliton. With no impurities and in one dimension (1D), these solitons would have an infinite lifetime, a consequence of integrability. In our experiment, the added impurities scatter off the much larger soliton, contributing to its Brownian motion and decreasing its lifetime. We describe the soliton's diffusive behavior using a quasi-1D scattering theory of impurity atoms interacting with a soliton, giving diffusion coefficients consistent with experiment.

Phantom vortices: hidden angular momentum in ultracold dilute Bose-Einstein condensates

Vortices are essential to angular momentum in quantum systems such as ultracold atomic gases. The existence of quantized vorticity in bosonic systems stimulated the development of the Gross-Pitaevskii mean-field approximation. However, the true dynamics of angular momentum in finite, interacting many-body systems like trapped Bose-Einstein condensates is enriched by the emergence of quantum correlations whose description demands more elaborate methods. Herein we theoretically investigate the full many-body dynamics of the acquisition of angular momentum by a gas of ultracold bosons in two dimensions using a standard rotation procedure. We demonstrate the existence of a novel mode of quantized vorticity, which we term the phantom vortex. Contrary to the conventional mean-field vortex, can be detected as a topological defect of spatial coherence, but not of the density. We describe previously unknown many-body mechanisms of vortex nucleation and show that angular momentum is hidden in phantom vortices modes which so far seem to have evaded experimental detection. This phenomenon is likely important in the formation of the Abrikosov lattice and the onset of turbulence in superfluids.

Sonic horizon formation for oscillating Bose-Einstein condensates in isotropic harmonic potential

We study the sonic horizon phenomena of the oscillating Bose-Einstein condensates in isotropic harmonic potential. Based on the Gross-Pitaevskii equation model and variational method, we derive the original analytical formula for the criteria and lifetime of the formation of the sonic horizon, demonstrating pictorially the interaction parameter dependence for the occur- rence of the sonic horizon and damping effect of the system distribution width. Our analytical results corroborate quantitatively the particular features of the sonic horizon reported in previous numerical study.

Traveling Majorana Solitons in a Low-Dimensional Spin-Orbit-Coupled Fermi Superfluid

We investigate traveling solitons of a one- or two-dimensional spin-orbit-coupled Fermi superfluid in both topologically trivial and nontrivial regimes by solving the static and time-dependent Bogoliubov-de Gennes equations. We find a critical velocity v_{h} for traveling solitons that is much smaller than the value predicted using the Landau criterion due to spin-orbit coupling. Above v_{h}, our time-dependent simulations in harmonic traps indicate that traveling solitons decay by radiating sound waves. In the topological phase, we predict the existence of peculiar Majorana solitons, which host two Majorana fermions and feature a phase jump of π across the soliton, irrespective of the velocity of travel. These unusual properties of Majorana solitons may open an alternative way to manipulate Majorana fermions for fault-tolerant topological quantum computations.

Dark soliton pair of ultracold Fermi gases for a generalized Gross-Pitaevskii equation model

We present the theoretical investigation of dark soliton pair solutions for one-dimensional as well as three-dimensional generalized Gross-Pitaevskii equation (GGPE) which models the ultracold Fermi gas during Bardeen-Cooper-Schrieffer-Bose-Einstein condensates crossover. Without introducing any integrability constraint and via the self-similar approach, the three-dimensional solution of GGPE is derived based on the one-dimensional dark soliton pair solution, which is obtained through a modified F-expansion method combined with a coupled modulus-phase transformation technique. We discovered the oscillatory behavior of the dark soliton pair from the theoretical results obtained for the three-dimensional case. The calculated period agrees very well with the corresponding reported experimental result [Weller et al., Phys. Rev. Lett. 101, 130401 (2008)PRLTAO0031-900710.1103/PhysRevLett.101.130401], demonstrating the applicability of the theoretical treatment presented in this work.

Dynamics of dipoles and vortices in nonlinearly coupled three-dimensional field oscillators

The dynamics of a pair of harmonic oscillators represented by three-dimensional fields coupled with a repulsive cubic nonlinearity is investigated through direct simulations of the respective field equations and with the help of the finite-mode Galerkin approximation (GA), which represents the two interacting fields by a superposition of 3+3 harmonic-oscillator p-wave eigenfunctions with orbital and magnetic quantum numbers l=1 and m=1, 0, -1. The system can be implemented in binary Bose-Einstein condensates, demonstrating the potential of the atomic condensates to emulate various complex modes predicted by classical field theories. First, the GA very accurately predicts a broadly degenerate set of the system's ground states in the p-wave manifold, in the form of complexes built of a dipole coaxial with another dipole or vortex, as well as complexes built of mutually orthogonal dipoles. Next, pairs of noncoaxial vortices and/or dipoles, including pairs of mutually perpendicular vortices, develop remarkably stable dynamical regimes, which feature periodic exchange of the angular momentum and periodic switching between dipoles and vortices. For a moderately strong nonlinearity, simulations of the coupled-field equations agree very well with results produced by the GA, demonstrating that the dynamics is accurately spanned by the set of six modes limited to l=1.

Solitons riding on solitons and the quantum Newton's cradle

The reduced dynamics for dark and bright soliton chains in the one-dimensional nonlinear Schrödinger equation is used to study the behavior of collective compression waves corresponding to Toda lattice solitons. We coin the term hypersoliton to describe such solitary waves riding on a chain of solitons. It is observed that in the case of dark soliton chains, the formulated reduction dynamics provides an accurate an robust evolution of traveling hypersolitons. As an application to Bose-Einstein condensates trapped in a standard harmonic potential, we study the case of a finite dark soliton chain confined at the center of the trap. When the central chain is hit by a dark soliton, the energy is transferred through the chain as a hypersoliton that, in turn, ejects a dark soliton on the other end of the chain that, as it returns from its excursion up the trap, hits the central chain repeating the process. This periodic evolution is an analog of the classical Newton's cradle.

Precession and nutation dynamics of nonlinearly coupled non-coaxial three-dimensional matter wave vortices

Nonlinearity is the driving force for numerous important effects in nature typically showing transitions between different regimes, regular, chaotic or catastrophic behavior. Localized nonlinear modes have been the focus of intense research in areas such as fluid and gas dynamics, photonics, atomic and solid state physics etc. Due to the richness of the behavior of nonlinear systems and due to the severe numerical demands of accurate three-dimensional (3D) numerical simulations presently only little knowledge is available on the dynamics of complex nonlinear modes in 3D. Here, we investigate the dynamics of 3D non-coaxial matter wave vortices that are trapped in a parabolic potential and interact via a repulsive nonlinearity. Our numerical simulations demonstrate the existence of an unexpected and fascinating nonlinear regime that starts immediately when the nonlinearity is switched-on and is characterized by a smooth dynamics representing torque-free precession with nutations. The reported motion is proven to be robust regarding various effects such as the number of particles, dissipation and trap deformations and thus should be observable in suitably designed experiments. Since our theoretical approach, i.e., coupled nonlinear Schrödinger equations, is quite generic, we expect that the obtained novel dynamical behavior should also exist in other nonlinear systems.

Cascade of Solitonic Excitations in a Superfluid Fermi gas: From Planar Solitons to Vortex Rings and Lines

We follow the time evolution of a superfluid Fermi gas of resonantly interacting ^{6}Li atoms after a phase imprint. Via tomographic imaging, we observe the formation of a planar dark soliton, its subsequent snaking, and its decay into a vortex ring, which, in turn, breaks to finally leave behind a single solitonic vortex. In intermediate stages, we find evidence for an exotic structure resembling the Φ soliton, a combination of a vortex ring and a vortex line. Direct imaging of the nodal surface reveals its undulation dynamics and its decay via the puncture of the initial soliton plane. The observed evolution of the nodal surface represents dynamics beyond superfluid hydrodynamics, calling for a microscopic description of unitary fermionic superfluids out of equilibrium.

Transitions from order to disorder in multiple dark and multiple dark-bright soliton atomic clouds

We have performed a systematic study quantifying the variation of solitary wave behavior from that of an ordered cloud resembling a "crystalline" configuration to that of a disordered state that can be characterized as a soliton "gas." As our illustrative examples, we use both one-component, as well as two-component, one-dimensional atomic gases very close to zero temperature, where in the presence of repulsive interatomic interactions and of a parabolic trap, a cloud of dark (dark-bright) solitons can form in the one- (two-) component system. We corroborate our findings through three distinct types of approaches, namely a Gross-Pitaevskii type of partial differential equation, particle-based ordinary differential equations describing the soliton dynamical system, and Monte Carlo simulations for the particle system. We define an "empirical" order parameter to characterize the order of the soliton lattices and study how this changes as a function of the strength of the "thermally" (i.e., kinetically) induced perturbations. As may be anticipated by the one-dimensional nature of our system, the transition from order to disorder is gradual without, apparently, a genuine phase transition ensuing in the intermediate regime.

Chladni solitons and the onset of the snaking instability for dark solitons in confined superfluids

Complex solitary waves composed of intersecting vortex lines are predicted in a channeled superfluid. Their shapes in a cylindrical trap include a cross, spoke wheels, and Greek Φ, and trace the nodal lines of unstable vibration modes of a planar dark soliton in analogy to Chladni's figures of membrane vibrations. The stationary solitary waves extend a family of solutions that include the previously known solitonic vortex and vortex rings. Their bifurcation points from the dark soliton indicating the onset of new unstable modes of the snaking instability are predicted from scale separation for Bose-Einstein condensates (BECs) and superfluid Fermi gases across the BEC-BCS crossover, and confirmed by full numerical calculations. Chladni solitons could be observed in ultracold gas experiments by seeded decay of dark solitons.

Dark solitons with Majorana fermions in spin-orbit-coupled Fermi gases

We show that a single dark soliton can exist in a spin-orbit-coupled Fermi gas with a high spin imbalance, where spin-orbit coupling favors uniform superfluids over nonuniform Fulde-Ferrell-Larkin-Ovchinnikov states, leading to dark soliton excitations in highly imbalanced gases. Above a critical spin imbalance, two topological Majorana fermions without interactions can coexist inside a dark soliton, paving a way for manipulating Majorana fermions through controlling solitons. At the topological transition point, the atom density contrast across the soliton suddenly vanishes, suggesting a signature for identifying topological solitons.