Generalized Gross-Pitaevskii Equation Research Articles

Propagation properties and stability of dark solitons in weakly interacting Bose–Bose droplets

We investigate dark solitons in two-component Bose systems with competing interactions in one dimension. Such a system hosts a liquid phase stabilized by the beyond-mean field corrections. Using the generalized Gross–Pitaevskii equation, we reveal the presence of two families of solitonic solutions. The solitons in both of them can be engineered to be arbitrarily wide. One family of solutions, however, has an anomalous dispersion relation, and our analyses show one of its branches is unstable. We find a critical velocity that demarcates the stable from unstable solutions. Nonetheless, gray anomalous solitons can exist inside quantum droplets and can be treated as solitonic excitations thereof.

Quantum droplet molecules in Bose-Bose mixtures

We investigate the interaction between two self-bound droplets separated by a finite distance in one-dimensional symmetric Bose-Bose mixtures. The main objective is to identify the regime where such two identical droplets support a stable bound-state known as droplet molecule. Our analysis is based on the generalized Gross-Pitaevskii equation, combining cubic and quadratic nonlinearities, which is solved by both variational method and numerical simulation. We calculate in particular the droplet-droplet interaction potential, the binding energy, the density profiles and the width of the droplet. It is shown that the developed variational approximation properly predicts the stability and the static and dynamic properties of such droplet molecules. Our results reveal that under appropriate adjustment of the system parameters such as the droplet separation and the number of atoms, a stable bound-state is formed in the regime of Gaussian-like shape. Nevertheless, in the flat top regime, the two droplets merge into one, leading to an unstable molecule.

Quantum droplets in one-dimensional Bose-Bose mixtures: beyond the Lee-Huang-Yang description

The static and dynamic properties of self-bound quantum droplets in a one-dimensional Bose-Bose mixture are discussed in the spirit of the Hartree–Fock-Bogoliubov theory. This latter enables us to provide beyond the Lee-Huang-Yang (LHY) quantum corrections to the equation of state at both zero and finite temperatures. In the uniform case our results for the ground-state energy and the critical temperature are confirmed through comparison with Quantum Monte-Carlo simulation and with available theoretical results. The density profiles are supported by numerical simulations of the generalized Gross-Pitaevskii equation which selfconsistently includes higher-order terms originating from the normal and anomalous fluctuations under the local density approximation. We show that the density exhibits a dip near its center in the flat-top plateau region for large interspecies interactions. We exemplify the impact of the beyond LHY corrections on the spatiotemporal evolution of the self-bound droplet in the presence of excitation induced by periodic density modulation. It is found that higher-order corrections may lead to the formation of a train of small droplets. We then extend our study for the case of inhomogeneous droplets in quasi one-dimensional Bose mixtures.

Quantum vortex formation in the "rotating bucket" experiment with polariton condensates.

The appearance of quantized vortices in the classical "rotating bucket" experiments of liquid helium and ultracold dilute gases provides the means for fundamental and comparative studies of different superfluids. Here, we realize the rotating bucket experiment for optically trapped quantum fluid of light based on exciton-polariton Bose-Einstein condensate in semiconductor microcavity. We use the beating note of two frequency-stabilized single-mode lasers to generate an asymmetric time-periodic rotating, nonresonant excitation profile that both injects and stirs the condensate through its interaction with a background exciton reservoir. The pump-induced external rotation of the condensate results in the appearance of a corotating quantized vortex. We investigate the rotation frequency dependence and reveal the range of stirring frequencies (from 1 to 4 GHz) that favors quantized vortex formation. We describe the phenomenology using the generalized Gross-Pitaevskii equation. Our results enable the study of polariton superfluidity on a par with other superfluids, as well as deterministic, all-optical control over structured nonlinear light.

Kardar-Parisi-Zhang universality in discrete two-dimensional driven-dissipative exciton polariton condensates

The statistics of the fluctuations of quantum many-body systems are highly revealing of their nature. In driven-dissipative systems displaying macroscopic quantum coherence, as exciton polariton condensates under incoherent pumping, the phase dynamics can be mapped to the stochastic Kardar-Parisi-Zhang (KPZ) equation. However, in two dimensions (2D), it was theoretically argued that the KPZ regime may be hindered by the presence of vortices, and a non-equilibrium BKT behavior was reported close to condensation threshold. We demonstrate here that, when a discretized 2D polariton system is considered, universal KPZ properties can emerge. We support our analysis by extensive numerical simulations of the discrete stochastic generalized Gross-Pitaevskii equation. We show that the first-order correlation function of the condensate exhibits stretched exponential behaviors in space and time with critical exponents characteristic of the 2D KPZ universality class, and find that the related scaling function accurately matches the KPZ theoretical one, stemming from functional Renormalization Group. We also obtain the distribution of the phase fluctuations and find that it is non-Gaussian, as expected for a KPZ stochastic process.

Global Behavior of Solutions to Generalized Gross-Pitaevskii Equation

This paper is concerned with time global behavior of solutions to nonlinear Schrödinger equation with a non-vanishing condition at the spatial infinity. Under a non-vanishing condition, it would be expected that the behavior is determined by the shape of the nonlinear term around the non-vanishing state. To observe this phenomenon, we introduce a generalized version of the Gross-Pitaevskii equation, which is a typical equation involving a non-vanishing condition, by modifying the shape of nonlinearity around the non-vanishing state. It turns out that, if the nonlinearity decays fast as a solution approaches to the non-vanishing state, then the equation admits a global solution which scatters to the non-vanishing element for both time directions.

Formation and fragmentation of quantum droplets in a quasi-one-dimensional dipolar Bose gas

We theoretically investigate the droplets formation in a tightly trapped one-dimensional dipolar gas of bosonic atoms. When the strength of the dipolar interaction becomes sufficiently attractive compared to the contact one, we show how a solitoniclike density profile evolves into a liquidlike droplet on increasing the number of particles in the trap. The incipient gas-liquid transition is also signaled by a steep increase of the breathing mode and a change in sign of the chemical potential. Upon a sudden release of the trap, varying the number of trapped atoms and the scattering length, the numerical solution of a time-dependent generalized Gross-Pitaevskii equation shows either an evaporation of the cloud, the formation of a single self-bound droplet, or a fragmentation in multiple droplets. These results can be probed with lanthanide atoms and help in characterizing the effect of the dipolar interaction in a quasi-one-dimensional geometry.

Soliton solutions of a generalized Gross-Pitaevskii equation via Darboux transform

The nonlinear Schrödinger equation (NLS) is widely used in different branches of physics. The GP equation is a kind of NLS equation with potential function term. In this paper, the soliton solutions of a generalized GP equation are explored. Firstly, we construct the Lax pair and Darboux transform(DT) of the equation. Then, we solve the single and double soliton solutions of the equation. Lastly, we draw the images of the single and double soliton solutions, and investigate the properties of the solitons. It is found that the amplitude of a single solitary wave does not change, and the shape and amplitude of the double solitary wave remain unchanged after the collision.

Vortex properties in the extended supersolid phase of dipolar Bose-Einstein condensates

We study the properties of singly-quantized linear vortices in the supersolid phase of a dipolar Bose-Einstein condensate at zero temperature modeling $^{164}$Dy atoms. The system is extended in the $x-y$ plane and confined by a harmonic trap in the the polarization direction $z$. Our study is based on a generalized Gross-Pitaevskii equation. We characterize the ground state of the system in terms of spatial order and superfluid fraction and compare the properties of a single vortex and of a vortex dipole in the superfluid phase (SFP) and in the supersolid phase (SSP). At variance with a vortex in the SFP, which is free to move in the superfluid, a vortex in the SSP is localized at the interstitial sites and does not move freely. We have computed the energy barrier for motion from an equilibrium site to another. The fact that the vortex is submitted to a periodic potential has a dramatic effect on the dynamics of a vortex dipole made of two counter rotating parallel vortices; instead of rigidly translating as in the SFP, the vortex and anti-vortex approach each other by a series of jumps from one site to another until they annihilate in a very short time and their energy is transferred to bulk excitations.

Spatiotemporal scaling of two-dimensional nonequilibrium exciton-polariton systems with weak interactions

We perform a theoretical study on the spatiotemporal correlation of two-dimensional (2D) nonequilibrium exciton-polariton (EP) systems driven by incoherent pumping. Based on the stochastic generalized Gross-Pitaevskii equation, we calculate the density fluctuation, coherence function, and scaling function. It is found that the correlations of short range agree well with the Bogoliubov linear theory, justifying the existence of an algebraic order. At large distances, both dynamic and static correlations are characterized by the nonlinear scaling behaviors of the Kardar-Parisi-Zhang universality class, especially when the interaction is weak. The coexistence of different types of correlations with their respective length scales is a distinctive feature of 2D EP systems. In addition, the interplay between vortices is modified in the regime of weak interaction and strong pumping, leading to complex nonequilibrium vortex patterns.

Kolmogorov and Kelvin wave cascades in a generalized model for quantum turbulence

We performed numerical simulations of decaying quantum turbulence by using a generalized Gross-Pitaevskii equation, that includes a beyond mean field correction and a nonlocal interaction potential. The nonlocal potential is chosen in order to mimic He II by introducing a roton minimum in the excitation spectrum. We observe that at large scales the statistical behavior of the flow is independent of the interaction potential, but at scales smaller than the intervortex distance a Kelvin wave cascade is enhanced in the generalized model. In this range, the incompressible kinetic energy spectrum obeys the weak wave turbulence prediction for Kelvin wave cascade not only for the scaling with wave numbers but also for the energy fluxes and the intervortex distance.

Temperature dependence of the density and excitations of dipolar droplets

Droplet states of ultracold gases which are stabilized by fluctuations have recently been observed for dipolar and two component Bose gases. These systems present a novel form of equilibrium where an instability at the mean field level is arrested by higher order correlations making the droplet states sensitive probes of fluctuations. In a recent paper, we argued that thermal fluctuations can play an important role for droplets even at low temperatures where the non-condensed density is much smaller than the condensate density. We used the Hartree-Fock-Bogoliubov theory together with local density approximation for fluctuations to obtain a generalized Gross Pitaevskii (GP) equation and solved it with a Gaussian variational ansatz to show that the transition between the low density and droplet states can be significantly modified by the temperature. In this paper, we first solve the same GP equation numerically with a time splitting spectral method to check the validity of the Gaussian variational ansatz. Our numerical results are in good agreement with the Gaussian ansatz for a large parameter regime and show that the density of the gas is most strongly modified by temperature near the abrupt transition between a pancake shaped cloud and the droplet. For cigar shaped condensates, as in the recent Er experiments, the dependence of the density on temperature remains quite small throughout the smooth transition. We then consider the effect of temperature on the collective oscillation frequencies of the droplet using both a time dependent Gaussian variational ansatz and real time numerical evolution. We find that the oscillation frequencies depend significantly on the temperature close to the transition for the experimentally relevant temperature regime ($\simeq 100$nK).

Quantized vortices in dipolar supersolid Bose-Einstein-condensed gases

We investigate the properties of quantized vortices in a dipolar Bose-Einstein condensed gas by means of a generalized Gross-Pitaevskii equation. The size of the vortex core hugely increases by increasing the weight of the dipolar interaction, approaching the transition to the supersolid phase. The critical angular velocity for the existence of an energetically stable vortex decreases in the supersolid, due to the reduced value of the density in the interdroplet region. The angular momentum per particle associated with the vortex line is shown to be smaller than $\ensuremath{\hbar}$, reflecting the reduction of the global superfluidity. The real-time vortex nucleation in a rotating trap is shown to be triggered, as for a standard condensate, by the softening of the quadrupole mode. For large angular velocities, when the distance between vortices becomes comparable to the interdroplet distance, the vortices are arranged into a honeycomb structure, which coexists with the triangular geometry of the supersolid lattice and persists during the free expansion of the atomic cloud.

Quantum droplets in one-dimensional Bose mixtures: A quantum Monte Carlo study

We use exact Quantum Monte Carlo techniques to study the properties of quantum droplets in two-component bosonic mixtures with contact interactions in one spatial dimension. We systematically study the surface tension, the density profile and the breathing mode as a function of the number of particles in the droplet and of the ratio of coupling strengths between intra-species repulsion and inter-species attraction. We find that deviations from the predictions of the generalized Gross-Pitaevskii equation are small in most cases of interest.

Balancing long-range interactions and quantum pressure: Solitons in the Hamiltonian mean-field model.

The Hamiltonian mean-field (HMF) model describes particles on a ring interacting via a cosine interaction, or equivalently, rotors coupled by infinite-range XY interactions. Conceived as a generic statistical mechanical model for long-range interactions such as gravity (of which the cosine is the first Fourier component), it has recently been used to account for self-organization in experiments on cold atoms with long-range optically mediated interactions. The significance of the HMF model lies in its ability to capture the universal effects of long-range interactions and yet be exactly solvable in the canonical ensemble. In this work we consider the quantum version of the HMF model in one dimension and provide a classification of all possible stationary solutions of its generalized Gross-Pitaevskii equation (GGPE), which is both nonlinear and nonlocal. The exact solutions are Mathieu functions that obey a nonlinear relation between the wave function and the depth of the mean-field potential, and we identify them as bright solitons. Using a Galilean transformation these solutions can be boosted to finite velocity and are increasingly localized as the mean-field potential becomes deeper. In contrast to the usual local GPE, the HMF case features a tower of solitons, each with a different number of nodes. Our results suggest that long-range interactions support solitary waves in a novel manner relative to the short-range case.

Violent relaxation in quantum fluids with long-range interactions.

Violent relaxation is a process that occurs in systems with long-range interactions. It has the peculiar feature of dramatically amplifying small perturbations, and rather than driving the system to equilibrium, it instead leads to slowly evolving configurations known as quasistationary states that fall outside the standard paradigm of statistical mechanics. Violent relaxation was originally identified in gravity-driven stellar dynamics; here, we extend the theory into the quantum regime by developing a quantum version of the Hamiltonian mean field (HMF) model which exemplifies many of the generic properties of long-range interacting systems. The HMF model can either be viewed as describing particles interacting via a cosine potential, or equivalently as the kinetic XY model with infinite-range interactions, and its quantum fluid dynamics can be obtained from a generalized Gross-Pitaevskii equation. We show that singular caustics that form during violent relaxation are regulated by interference effects in a universal way described by Thom's catastrophe theory applied to waves and this leads to emergent length scales and timescales not present in the classical problem. In the deep quantum regime we find that violent relaxation is suppressed altogether by quantum zero-point motion. Our results are relevant to laboratory studies of self-organization in cold atomic gases with long-range interactions.

Kardar-Parisi-Zhang universality in the phase distributions of one-dimensional exciton-polaritons

Exciton-polaritons under driven-dissipative conditions exhibit a condensation transition which belongs to a different universality class than equilibrium Bose-Einstein condensates. By numerically solving the generalized Gross-Pitaevskii equation with realistic experimental parameters, we show that one-dimensional exciton-polaritons display fine features of Kardar-Parisi-Zhang (KPZ) dynamics. Beyond the scaling exponents, we show that their phase distribution follows the Tracy-Widom form predicted for KPZ growing interfaces. We moreover evidence a crossover to the stationary Baik-Rains statistics. We finally show that these features are unaffected on a certain timescale by the presence of a smooth disorder often present in experimental setups.

Tunable band-gap structure and gap solitons in the generalized Gross-Pitaevskii equation with a periodic potential

The tunable band-gap structure is fundamentally important in the dynamics of both linear and nonlinear modes trapped in a lattice because Bloch modes can only exist in the bands of the periodic system and nonlinear modes associating with them are usually confined to the gaps. We reveal that when a momentum operator is introduced into the Gross-Pitaevskii equation (GPE), the bandgap spectra of the periodic system can be shifted upward parabolically by the growth of the constant momentum coefficient. During this process, the band edges become asymmetric, in sharp contrast to the standard GPE with an external periodic potential. Extended complex Bloch modes with asymmetric profiles can be derived by applying a phase transformation to the symmetric profiles. We find that the inherent parity-time symmetry of the complex system is never broken with increasing momentum coefficient. Under repulsive interactions, solitons with different numbers of peaks bifurcating from the band edges are found in finite gaps. We also address the existence of embedded solitons in the generalized two-dimensional GPE. Linear stability analysis corroborated by direct evolution simulations demonstrates that multi-peaked solitons are almost completely stable in their entire existence domains.

Stable switching among high-order modes in polariton condensates

We report multistate optical switching among high-order bouncing-ball modes ("ripples") and whispering-gallerying modes ("petals") of exciton-polariton condensates in a laser-generated annular trap. By tailoring the diameter and power of the annular trap, the polariton condensate can be switched among different trapped modes, accompanied by redistribution of spatial densities and superlinear increase in the emission intensities, implying that polariton condensates in this geometry could be exploited for a multistate switch. A model based on non-Hermitian modes of the generalized Gross-Pitaevskii equation reveals that this mode switching arises from competition between pump-induced gain and in-plane polariton loss. The parameters for reproducible switching among trapped modes have been measured experimentally, giving us a phase diagram for mode switching. Taken together, the experimental result and theoretical modeling advances our fundamental understanding of the spontaneous emergence of coherence and move us toward its practical exploitation.

Control of propagation of spatially localized polariton wave packets in a Bragg mirror with embedded quantum wells

We considered the nonlinear dynamics of Bragg polaritons in a specially designed stratified semiconductor structure with embedded quantum wells, which possesses a convex dispersion. The model for the ensemble of single periodically arranged quantum wells coupled with the Bragg photon fields has been developed. In particular, the generalized Gross-Pitaevskii equation with the non-parabolic dispersion has been obtained for the Bragg polariton wave function. We revealed a number of dynamical regimes for polariton wave packets resulting from competition of the convex dispersion and the repulsive nonlinearity effects. Among the regimes are spreading, breathing and soliton propagation. When the control parameters including the exciton-photon detuning, the matter-field coupling and the nonlinearity are manipulated, the dynamical regimes switch between themselves.