Stochastic Gross-Pitaevskii Equation Research Articles

Comparative analysis of numerical with optical soliton solutions of stochastic Gross–Pitaevskii equation in dispersive media

This article deals with the stochastic Gross–Pitaevskii equation (SGPE) perturbed with multiplicative time noise. The numerical solutions of the governing model are carried out with the proposed stochastic non-standard finite difference (SNSFD) scheme. The stability of the scheme is proved by using the Von-Neumann criteria and the consistency is shown in the mean square sense. To seek exact solutions, we applied the Sardar subequation (SSE) and modified exponential rational functional (MERF) techniques. The exact solutions are constructed in the form of exponential, hyperbolic, and trigonometric forms. Finally, the comparison of the exact solutions with numerical solutions is drawn in the 3D and line plots for the different values of parameters.

Formation of exciton-polaritonic BEC in the non-Markovian regime

Dynamics of exciton-polaritonic Bose-Einstein condensate is considered. We construct the stochastic Gross-Pitaevskii equation, where condensate is coupled to the excitonic reservoir in the non-Markovian regime. This equation is utilized for study of condensate formation. It is found that increasing of temperature results in transition from spatially uniform to fragmented pattern with onset of numerous vortices. The transition temperature corresponds to minimum of condensate density. Below the transition temperature, density decreases with increasing of temperature due to noise-induced loss of phase coherence. Increasing of density above the transition temperature is caused by suppression of memory that facilitates nearly-exponential density growth consistent with the Markovian regime.

Phase transition dimensionality crossover from two to three dimensions in a trapped ultracold atomic Bose gas

The equilibrium properties of a weakly interacting atomic Bose gas across the Berezinskii-Kosterlitz-Thouless (BKT) and Bose-Einstein condensation (BEC) phase transitions are numerically investigated through a dimensionality crossover from two to three dimensions. The crossover is realized by confining the gas in an experimentally feasible hybridized trap which provides homogeneity along the planar $xy$ directions through a box potential in tandem with a harmonic transverse potential along the transverse $z$ direction. The dimensionality is modified by varying the frequency of the harmonic trap from tight to loose transverse trapping. Our findings, based on a stochastic (projected) Gross-Pitaevskii equation, showcase a continuous shift in the character of the phase transition from BKT to BEC, and a monotonic increase of the identified critical temperature as a function of dimensionality, with the strongest variation exhibited for small chemical potential values up to approximately twice the transverse confining potential.

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.

Finite-temperature spin dynamics of a two-dimensional Bose-Bose atomic mixture

We examine the role of thermal fluctuations in uniform two-dimensional binary Bose mixtures of dilute ultracold atomic gases. We use a mean-field Hartree-Fock theory to derive analytical predictions for the miscible-immiscible transition. A nontrivial result of this theory is that a fully miscible phase at $T=0$ may become unstable at $T\neq0$, as a consequence of a divergent behaviour in the spin susceptibility. We test this prediction by performing numerical simulations with the Stochastic (Projected) Gross-Pitaevskii equation, which includes beyond mean-field effects. We calculate the equilibrium configurations at different temperatures and interaction strengths and we simulate spin oscillations produced by a weak external perturbation. Despite some qualitative agreement, the comparison between the two theories shows that the mean-field approximation is not able to properly describe the behavior of the two-dimensional mixture near the miscible-immiscible transition, as thermal fluctuations smoothen all sharp features both in the phase diagram and in spin dynamics, except for temperature well below the critical temperature for superfluidity.

Optical orientation, polarization pinning, and depolarization dynamics in optically confined polariton condensates

We investigate the optical orientation, polarization pinning, and depolarization of optically confined semiconductor exciton-polariton condensates. We perform a complete mapping of the condensate polarization as a function of incident nonresonant excitation polarization and power. We utilize a ring-shaped excitation pattern to generate an exciton-induced potential that spatially confines polariton condensates into a single mode. We observe that formation of circular polarization in the condensate persists even for a weakly cocircularly polarized pump. By varying the excitation ring diameter we realize a transition from the condensate polarization being pinned along the coordinate-dependent cavity-strain axes, to a regime of zero degree of condensate polarization. Analysis through the driven-dissipative stochastic Gross-Pitaevskii equation reveals that this depolarization stems from a competition between sample induced in-plane polarization splitting and the condensate-reservoir overlap. An increase in the role of the latter results in weakening of the condensate fixed-point phase space attractors, and enhanced random phase space walk and appearance of limit cycle trajectories, reducing the degree of time-integrated polarization.

Kibble-Zurek dynamics in a trapped ultracold Bose gas

The dynamical evolution of an inhomogeneous ultracold atomic gas quenched at different controllable rates through the Bose-Einstein condensation phase transition is studied numerically in the premise of a recent experiment in an anisotropic harmonic trap. Our findings based on the stochastic (projected) Gross-Pitaevskii equation are shown to be consistent at early times with the predictions of the homogeneous Kibble-Zurek mechanism. This is demonstrated by collapsing the early dynamical evolution of densities, spectral functions and correlation lengths for different quench rates, based on an appropriate characterization of the distance to criticality felt by the quenched system. The subsequent long-time evolution, beyond the identified dynamical critical region, is also investigated by looking at the behaviour of the density wavefront evolution and the corresponding phase ordering dynamics.

Collisionless Sound in a Uniform Two-Dimensional Bose Gas.

Using linear response theory within the random phase approximation, we investigate the propagation of sound in a uniform two dimensional (2D) Bose gas in the collisionless regime. We show that the sudden removal of a static density perturbation produces a damped oscillatory behavior revealing that sound can propagate also in the absence of collisions, due to mean-field interaction effects. We provide explicit results for the sound velocity and damping as a function of temperature, pointing out the crucial role played by Landau damping. We support our predictions by performing numerical simulations with the stochastic (projected) Gross-Pitaevskii equation. The results are consistent with the recent experimental observation of sound in a weakly interacting 2D Bose gas both below and above the superfluid Berezinskii-Kosterlitz-Thouless transition.

Continuum of classical-field ensembles in Bose gases from canonical to grand canonical and the onset of their equivalence

The canonical and grand-canonical ensembles are two usual marginal cases for ultracold Bose gases, but real collections of experimental runs commonly have intermediate properties. Here we study the continuum of intermediate cases and look into the appearance of ensemble equivalence as interaction rises for mesoscopic $1d$ systems. We demonstrate how at sufficient interaction strength the distributions of condensate and excited atoms become practically identical, regardless of the ensemble used. Importantly, we find that features that are fragile in the ideal gas and appear only in a strict canonical ensemble can become robust in all ensembles when interactions become strong. As evidence, the steep cliff in the distribution of the number of excited atoms is preserved. To make this study, a straightforward approach for generating canonical and intermediate classical field ensembles using a modified stochastic Gross-Pitaevskii equation is developed.

Thermal quenches in the stochastic Gross-Pitaevskii equation: Morphology of the vortex network

We study the evolution of 3d weakly interacting bosons at finite chemical potential with the stochastic Gross-Pitaevskii equation. We fully characterise the vortex network in an out of equilibrium. At high temperature the filament statistics are the ones of fully-packed loop models. The vortex tangle undergoes a geometric percolation transition within the thermodynamically ordered phase. After infinitely fast quenches across the thermodynamic critical point deep into the ordered phase, we identify a first approach towards a state that is numerically indistinguishable from the one of the critical threshold, a later coarsening process that does not alter the fractal properties of the long vortex loops, and a final approach to equilibrium. Our results are also relevant to the statistics of linear defects in type-II superconductors, magnetic materials and cosmological models.

Quantum fluctuation effects on the quench dynamics of thermal quasicondensates

We study the influence of quantum fluctuations on the phase, density, and pair correlations in a trapped quasicondensate after a quench of the interaction strength. To do so, we derive a description similar to the stochastic Gross–Pitaevskii equation (SGPE) but keeping a fully quantum description of the low-energy fields using the positive-P representation. This allows us to treat both the quantum and thermal fluctuations together in an integrated way. A plain SGPE only allows for thermal fluctuations. The approach is applicable to such situations as finite temperature quantum quenches, but not equilibrium calculations due to the time limitations inherent in positive-P descriptions of interacting gases. One sees the appearance of antibunching, the generation of counter-propagating atom pairs, and increased phase fluctuations. We show that the behavior can be estimated by adding the T = 0 quantum fluctuation contribution to the thermal fluctuations described by the plain SGPE.

Quantum Langevin model for nonequilibrium condensation

We develop a quantum model for nonequilibrium Bose-Einstein condensation of photons and polaritons in planar microcavity devices. The model builds on laser theory and includes the spatial dynamics of the cavity field, a saturation mechanism, and some frequency dependence of the gain: quantum Langevin equations are written for a cavity field coupled to a continuous distribution of externally pumped two-level emitters with a well-defined frequency. As an example of application, the method is used to study the linearized quantum fluctuations around a steady-state condensed state. In the good-cavity regime, an effective equation for the cavity field only is proposed in terms of a stochastic Gross-Pitaevskii equation. Perspectives in view of a full quantum simulation of the nonequilibrium condensation process are finally sketched.

Stochastic projected Gross-Pitaevskii equation for spinor and multicomponent condensates

A stochastic Gross-Pitaevskii equation is derived for partially condensed Bose gas systems subject to binary contact interactions. The theory we present provides a classical-field theory suitable for describing dissipative dynamics and phase transitions of spinor and multi-component Bose gas systems comprised of an arbitrary number of distinct interacting Bose fields. A new class of dissipative processes involving distinguishable particle interchange between coherent and incoherent regions of phase-space is identified. The formalism and its implications are illustrated for two-component mixtures and spin-1 Bose-Einstein condensates. For systems comprised of atoms of equal mass, with thermal reservoirs that are close to equilibrium, the dissipation rates of the theory are reduced to analytical expressions that may be readily evaluated. The unified treatment of binary contact interactions presented here provides a theory with broad relevance for quasi-equilibrium and far-from-equilibrium Bose-Einstein condensates.

Finite-temperature dynamics of vortices in Bose-Einstein condensates

We study the dynamics of a single and a pair of vortices in quasi two-dimensional Bose-Einstein condensates at finite temperatures. We use the stochastic Gross-Pitaevskii equation, which is the Langevin equation for the Bose-Einstein condensate, to this end. For a pair of vortices, we study the dynamics of both the vortex-vortex and vortex-antivortex pairs, which are generated by rotating the trap and moving the Gaussian obstacle potential, respectively. Due to thermal fluctuations, the constituent vortices are not symmetrically generated with respect to each other at finite temperatures. This initial asymmetry coupled with the presence of random thermal fluctuations in the system can lead to different decay rates for the component vortices of the pair, especially in the case of two corotating vortices.

Ab initiomethods for finite-temperature two-dimensional Bose gases

The stochastic Gross-Pitaevskii equation and modified Popov theory are shown to provide an ab initio description of finite temperature, weakly-interacting two-dimensional Bose gas experiments. Using modified Popov theory, a systematic approach is developed in which the momentum cut-off inherent to classical field methods is removed as a free parameter. This is shown to yield excellent agreement with the recent experiment of Hung et al. [Nature, 470, 236 (2011)], verifying that the stochastic Gross-Pitaevskii equation captures the observed universality and scale-invariance.

Phase coherence in quasicondensate experiments: Anab initioanalysis via the stochastic Gross-Pitaevskii equation

We perform an ab initio analysis of the temperature dependence of the phase coherence length of finite temperature, quasi-one-dimensional Bose gases measured in the experiments of Richard et al. (Phys. Rev. Lett. 91, 010405 (2003)) and Hugbart et al. (Eur. Phys. J. D 35, 155-163 (2005)), finding very good agreement across the entire observed temperature range ($0.8<T/T_{\phi}<28$). Our analysis is based on the one-dimensional stochastic Gross-Pitaevskii equation, modified to self-consistently account for transverse, quasi-one-dimensional effects, thus making it a valid model in the regime $\mu ~ few \hbar \omega_\perp$. We also numerically implement an alternative identification of $T_{\phi}$, based on direct analysis of the distribution of phases in a stochastic treatment.

Winding up superfluid in a torus via Bose Einstein condensation.

Phase transitions are usually treated as equilibrium phenomena, which yields telltale universality classes with scaling behavior of relaxation time and healing length. However, in second-order phase transitions relaxation time diverges near the critical point (“critical slowing down”). Therefore, every such transition traversed at a finite rate is a non-equilibrium process. Kibble-Zurek mechanism (KZM) captures this basic physics, predicting sizes of domains – fragments of broken symmetry – and the density of topological defects, long-lived relics of symmetry breaking that can survive long after the transition. To test KZM we simulate Bose-Einstein condensation in a ring using stochastic Gross-Pitaevskii equation and show that BEC formation can spontaneously generate quantized circulation of the newborn condensate. The magnitude of the resulting winding numbers and the time-lag of BEC density growth – both experimentally measurable – follow scalings predicted by KZM. Our results may also facilitate measuring the dynamical critical exponent for the BEC transition.

Fluctuating and dissipative dynamics of dark solitons in quasicondensates

The fluctuating and dissipative dynamics of matter-wave dark solitons within harmonically trapped, partially condensed Bose gases is studied both numerically and analytically. A study of the stochastic Gross-Pitaevskii equation, which correctly accounts for density and phase fluctuations at finite temperatures, reveals dark-soliton decay times to be lognormally distributed at each temperature, thereby characterizing the previously predicted long-lived soliton trajectories within each ensemble of numerical realizations [S. P. Cockburn et al., Phys. Rev. Lett. 104, 174101 (2010)]. Expectation values for the average soliton lifetimes extracted from these distributions are found to agree well with both numerical and analytic predictions based upon the dissipative Gross-Pitaevskii model (with the same ab initio damping). Probing the regime for which 0.8 ${k}_{B}T&lt;\ensuremath{\mu}&lt;1.6$ ${k}_{B}T$, we find average soliton lifetimes to scale with temperature as $\ensuremath{\tau}\ensuremath{\sim}{T}^{\ensuremath{-}4}$, in agreement with predictions previously made for the low-temperature regime ${k}_{B}T\ensuremath{\ll}\ensuremath{\mu}$. The model is also shown to capture the experimentally relevant decrease in the visibility of an oscillating soliton due to the presence of background fluctuations.

Quantitative study of quasi-one-dimensional Bose gas experiments via the stochastic Gross-Pitaevskii equation

The stochastic Gross-Pitaevskii equation is shown to be an excellent model for quasi-one-dimensional Bose gas experiments, accurately reproducing the in situ density profiles recently obtained in the experiments of Trebbia et al. [Phys. Rev. Lett. 97, 250403 (2006)] and van Amerongen et al. [Phys. Rev. Lett. 100, 090402 (2008)], and the density fluctuation data reported by Armijo et al. [Phys. Rev. Lett. 105, 230402 (2010)]. To facilitate such agreement, we propose and implement a quasi-one-dimensional stochastic equation for the low-energy, axial modes, while atoms in excited transverse modes are treated as independent ideal Bose gases.

Comparison between microscopic methods for finite-temperature Bose gases

We analyze the equilibrium properties of a weakly interacting, trapped quasi-one-dimensional Bose gas at finite temperatures and compare different theoretical approaches. We focus in particular on two stochastic theories: a number-conserving Bogoliubov (ncB) approach and a stochastic Gross-Pitaevskii equation (sGPe) that have been extensively used in numerical simulations. Equilibrium properties like density profiles, correlation functions, and the condensate statistics are compared to predictions based upon a number of alternative theories. We find that due to thermal phase fluctuations, and the corresponding condensate depletion, the ncB approach loses its validity at relatively low temperatures. This can be attributed to the change in the Bogoliubov spectrum, as the condensate gets thermally depleted, and to large fluctuations beyond perturbation theory. Although the two stochastic theories are built on different thermodynamic ensembles (ncB: canonical, sGPe: grand-canonical), they yield the correct condensate statistics in a large BEC (strong enough particle interactions). For smaller systems, the sGPe results are prone to anomalously large number fluctuations, well-known for the grand-canonical, ideal Bose gas. Based on the comparison of the above theories to the modified Popov approach, we propose a simple procedure for approximately extracting the Penrose-Onsager condensate from first- and second-order correlation functions that is computationally convenient. This also clarifies the link between condensate and quasi-condensate in the Popov theory of low-dimensional systems.