Gross-Pitaevskii Equation Research Articles

Uncertainty Relations in Hydrodynamics

The qualitative behaviors of uncertainty relations in hydrodynamics are numerically studied for fluids with low Reynolds numbers in 1+1 dimensional system. We first give a review for the formulation of the generalized uncertainty relations in the stochastic variational method (SVM), following the work by two of the present authors [Phys. Lett. A 382, 1472 (2018)]. In this approach, the origin of the finite minimum value of uncertainty is attributed to the non-differentiable (virtual) trajectory of a quantum particle and then both of the Kennard and Robertson-Schrödinger inequalities in quantum mechanics are reproduced. The same non-differentiable trajectory is applied to the motion of fluid elements in the Navier-Stokes-Fourier equation or the Navier-Stokes-Korteweg equation. By introducing the standard deviations of position and momentum for fluid elements, the uncertainty relations in hydrodynamics are derived. These are applicable even to the Gross-Pitaevskii equation and then the field-theoretical uncertainty relation is reproduced. We further investigate numerically the derived relations and find that the behaviors of the uncertainty relations for liquid and gas are qualitatively different. This suggests that the uncertainty relations in hydrodynamics are used as a criterion to classify liquid and gas in fluid.

Ground states for a nonlocal mixed order cubic-quartic Gross-Pitaevskii equation

We prove existence and qualitative properties of ground state solutions to a generalized nonlocal 3rd-4th order Gross-Pitaevskii equation. Using a mountain pass argument on spheres and constructing appropriately localized Palais-Smale sequences we are able to prove existence of real positive ground states as saddle points. The analysis is deployed in the set of possible states, thus overcoming the problem that the energy is unbounded below. We also prove a corresponding nonlocal Pohozaev identity with no rest term, a crucial part of the analysis.

Deflation-based identification of nonlinear excitations of the three-dimensional Gross-Pitaevskii equation

We present previously unknown solutions to the 3D Gross--Pitaevskii equation describing atomic Bose-Einstein condensates. This model supports elaborate patterns, including excited states bearing vorticity. The discovered coherent structures exhibit striking topological features, involving combinations of vortex rings and multiple, possibly bent vortex lines. Although unstable, many of them persist for long times in dynamical simulations. These solutions were identified by a state-of-the-art numerical technique called deflation, which is expected to be applicable to many problems from other areas of physics.

Supersolid striped droplets in a Raman spin-orbit-coupled system

We analyze the role played by quantum fluctuations on a Raman Spin-Orbit Coupled system in the stripe phase. We show that beyond mean-field effects stabilize the collapse predicted by mean-field theory and induce the emergence of two phases: a gas and a liquid, which also show spatial periodicity along a privileged direction. We show that the energetically favored phase is determined by the Raman coupling and the spin-dependent scattering lengths. We obtain the ground-state solution of the finite system by solving the extended Gross-Pitaevskii equation and find self-bound, droplet-like solutions that feature internal structure through a striped pattern. We estimate the critical number for binding associated to these droplets and show that their value is experimentally accessible. We report an approximate energy functional in order to ease the evaluation of the Lee-Huang-Yang correction in practical terms.

Collective excitations of a spherical ultradilute quantum droplet

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

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.

Conserved energies for the one dimensional Gross-Pitaevskii equation

We prove the well-posedness results for the one dimensional Gross-Pitaevskii equation in the energy space, which is a complete metric space equipped with a newly introduced metric and with the energy norm describing the Hs regularities of the solutions. We establish a family of conserved energies for the one dimensional Gross-Pitaevskii equation, such that all the energy norms of the solutions are conserved globally in time. This family of energies is also conserved by the complex modified Korteweg-de Vries flow.

Quantum Gases of Dipoles, Quadrupoles and Octupoles in Gross–Pitaevskii Formalism with Form Factor

Here, classical and quantum field theory of dipolar, axisymmetric quadrupolar and octupolar Bose gases is considered within a general approach. Dipole, axisymmetric quadrupole and octupole interaction potentials in the momentum representation are calculated. These results clearly demonstrate attraction and repulsion areas in corresponding gases. Then the Gross–Pitaevskii (GP) equation, which plays a key role in the present paper, is derived from the corresponding functional. The zoology of the form factors appearing in the GP equation is studied in details. The proper classes for the description of spatially non-uniform condensates form factors are chosen. In the Thomas–Fermi approximation a general solution of the GP equation with a quasilocal form factor is obtained. This solution has an interesting form in terms of a double rapidly converging series that universally includes all the interactions considered. Plots of condensate density functions for the exponential-trigonometric form factor are given. For the sake of completeness, in this paper we consider the GP equation with an optical lattice potential in the limit of small condensate densities. This limit does not distinguish between dipolar, quadrupolar and octupolar gases. An important analysis of the condensate stability, in other words the study of condensate excitations, is also performed in this paper. In the Gaussian approximation (from the Gross–Pitaevskii functional), a functional describing the perturbations of the condensate is derived in detail. This problem is an analog of the Bogolubov transformation used in the study of quantum Bose gases in operator formalism. For a probe wave function in the form of a plane wave, a spectrum of (Bogoliubov) excitations was obtained, from which an equation describing the threshold momentum for the emergence of instability was derived. An important result of this paper is the dependence of the threshold on the momentum of a stationary condensate. For completeness of the presentation, the approximating expression in the form of a rapidly converging series is obtained for the corresponding dependence, and plots of the corresponding series for the exponential-trigonometric form factor are given. Finally, in the conclusion a quantum hydrodynamic theory for dipolar, axisymmetric quadrupolar and octupolar gases is briefly presented, giving a clue to the experimental determination of the form factors.

Kinetics and luminescence of the excitations of a nonequilibrium polariton condensate

We extend the description of the nonequilibrium Gross-Pitaevskii equation for microcavity polaritons by including not only the kinetics of the nonresonantly pumped exciton reservoir but also that of the collective excitations---called bogolons---above the condensate. Quasiequilibrium Hartree-Fock-Bogoliubov-Popov self-energies are used to describe the condensate and the bogolons. With nonequilibrium Green functions, the quantum Boltzmann equation of the bogolons in a frame moving with the condensate, together with the coupled Gross-Pitaevskii equation can be derived. The resulting scattering rates contain vertex corrections in form of symmetrized products of the $u$ and $v$ pairing functions of all involved scattering states. Corresponding results have been derived for ultracold atoms by A. Griffin, T. Nikuni, and E. Zaremba [Bose-Condensed Gases at Finite Temperatures, 1st ed. (Cambridge University Press, Cambridge, New York, 2009)] Furthermore, we show that the results of the quantum kinetics also determine the luminescence spectrum including the ghost branch. The kinetics is solved numerically, the resulting condensate, the bogolon distribution, and the luminescence are calculated temporally and spatially resolved after the switch-on of a $\mathit{cw}$ nonresonant pump beam with a Gaussian profile. It is studied in detail to which extent the bogolons approach a local thermal equilibrium.

Double-layer Bose-Einstein condensates: A quantum phase transition in the transverse direction, and reduction to two dimensions.

We revisit the problem of the reduction of the three-dimensional (3D) dynamics of Bose-Einstein condensates, under the action of strong confinement in one direction (z), to a 2D mean-field equation. We address this problem for the confining potential with a singular term, viz., V_{z}(z)=2z^{2}+ζ^{2}/z^{2}, with constant ζ. A quantum phase transition is induced by the latter term, between the ground state (GS) of the harmonic oscillator and the 3D condensate split in two parallel noninteracting layers, which is a manifestation of the "superselection" effect. A realization of the respective physical setting is proposed, making use of resonant coupling to an optical field, with the resonance detuning modulated along z. The reduction of the full 3D Gross-Pitaevskii equation (GPE) to the 2D nonpolynomial Schrödinger equation (NPSE) is based on the factorized ansatz, with the z -dependent multiplier represented by an exact GS solution of the 1D Schrödinger equation with potential V_{z}(z). For both repulsive and attractive signs of the nonlinearity, the 2D NPSE produces GS and vortex states, that are virtually indistinguishable from the respective numerical solutions provided by full 3D GPE. In the case of the self-attraction, the threshold for the onset of the collapse, predicted by the 2D NPSE, is also virtually identical to its counterpart obtained from the 3D equation. In the same case, stability and instability of vortices with topological charge S=1, 2, and 3 are considered in detail. Thus, the procedure of the spatial-dimension reduction, 3D → 2D, produces very accurate results, and it may be used in other settings.

Dynamics of mode-locked nanolasers based on Hermite-Gaussian modes

The different dynamical behaviors of the Hermite-Gaussian (HG) modes of mode-locked nanolasers based on a harmonic photonic cavity are investigated in detail using a model based on a modified Gross-Pitaevskii Equation. Such nanolasers are shown to exhibit mode-locking with a repetition rate independent of the cavity length, which is a strong asset for compactness.The differences with respect to conventional lasers are shown to originate from the peculiar gain competition between HG modes, which is investigated in details. In the presence of a saturable absorber, the different regimes, i. e. Q-switching, Q-switched mode-locking, and continuous-wave (cw) mode locking, are isolated in a phase diagram and separately described. Mode-locking is found to be robust against phase-intensity coupling and to be achievable in a scheme with spatially separated gain and absorber.

Bose-Einstein condensate dark matter model with three-particle interaction and two-phase structure

We explore the consequences of including the repulsive three-particle interaction in the model of Bose-Einstein condensate dark matter model or fuzzy dark matter. Such a model based on properly modified Gross-Pitaevskii equation is intended to describe the distribution of dark matter particles in the highly dense regions, which correspond to the galaxy core and/or to the overlap of colliding galaxies. Specifically, we deal with the $\phi^6$-model in terms of the macroscopic wave function of the condensate, where a locality of interaction is guaranteed by a large correlation length assumed to hold. After calculation of main thermodynamical characteristics, we find strong evidence of the existence of two distinct phases of dark matter, within its core, separated by the instability region lying between two differing special values of the pressure acting in the model. Some implications stemming from the existence of two phases and the related first-order phase transition are discussed.

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

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

Microscopic derivation of the extended Gross-Pitaevskii equation for quantum droplets in binary Bose mixtures

An ultradilute quantum droplet is a self-bound liquid-like state recently observed in ultracold Bose-Einstein condensates. In all previous theoretical studies, it is described by a phenomenological low-energy effective theory, termed as the extended Gross\textendash Pitaevskii equation. Here, we microscopically derive the Gross\textendash Pitaevskii equation for the condensate and also for a pairing field in an inhomogeneous quantum droplet realized by Bose-Bose mixtures with attractive inter-species interaction. We show that the inclusion of the pairing field is essential, in order to have a consistent description of the droplet state. We clarify that, the extended Gross\textendash Pitaevskii equation used earlier should be understood as the equation of motion for the pairing field, rather than the condensate. The fluctuations of the pairing field give rise to low-energy collective excitations of the droplet. We also present the Bogoliubov equations for gapless phonon modes and gapped modes due to pairing in real space, which characterizes single-particle-like excitations of the droplet. The equations of motion derived in this work for the condensate and the pairing field serve an ideal starting point to understand the structure and collective excitations of non-uniform ultradilute quantum droplets in on-going cold-atom experiments.

Effective self-similar expansion of a Bose-Einstein condensate: Free space versus confined geometries

We compare the exact evolution of an expanding three-dimensional Bose-Einstein condensate with that obtained from the effective scaling approach introduced in D. Gu\'ery-Odelin [Phys. Rev. A 66, 033613 (2002)]. This approach, which consists in looking for self-similar solutions to be satisfied on average, is tested here in different geometries and configurations. We find that, in case of almost isotropic traps, the effective scaling reproduces with high accuracy the exact evolution dictated by the Gross-Pitaevskii equation for arbitrary values of the interactions, in agreement with the proof-of-concept of M. Modugno, G. Pagnini, and M. A. Valle-Basagoiti [Phys. Rev. A 97, 043604 (2018)]. Conversely, it is shown that the hypothesis of universal self-similarity breaks down in case of strong anisotropies and trapped geometries.

Interaction-Enhanced Group Velocity of Bosons in the Flat Band of an Optical Kagome Lattice.

Geometric frustration of particle motion in a kagome lattice causes the single-particle band structure to have a flat s-orbital band. We probe this band structure by placing a Bose-Einstein condensate into excited Bloch states of an optical kagome lattice, and then measuring the group velocity through the atomic momentum distribution. We find that interactions renormalize the band structure, greatly increasing the dispersion of the third band, which is nearly non-dispersing the single-particle treatment. Calculations based on the lattice Gross-Pitaevskii equation indicate that band structure renormalization is caused by the distortion of the overall lattice potential away from the kagome geometry by interactions.

Commensurate-incommensurate supersolid ground state of a spin-orbit-coupled Bose-Einstein condensate in one-dimensional optical lattices

Supersolid, a fascinating quantum state of matter, features novel phenomena such as the non-classical rotational inertia and transport anomalies. It is a long standing issue of the coexistence of superfluidity and broken translational symmetry in condensed matter physics. By recent experimental advances to create tunable synthetic spin-orbit coupling in ultracold gases, such highly controllable atomic systems would provide new possibilities to access supersolidity with no counterpart in solids. Here we report that the combination of anharmonicity of trapping potential and spin-orbit coupling will provide a new paradigm to achieve supersolids. By means of imaginary time evolution of the Gross-Pitaevskii equation, we demonstrate that a supersolid state can be found when considering a trapped Rashba-type spin-orbit coupled bosonic atoms loaded in a one-dimensional optical lattice. Furthermore, a skyrmion-anti-skyrmion lattice is associated with the appearance of such supersoildity, indicating the topological nontrivial properties of our proposed supersolids.

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

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

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.

Numerical method for the projected Gross-Pitaevskii equation in an infinite rotating two-dimensional Bose gas.

We present a method for evolving the projected Gross-Pitaevskii equation in an infinite rotating Bose-Einstein condensate, the ground state of which is a vortex lattice. We use quasiperiodic boundary conditions to investigate the behavior of the bulk superfluid in this system, in the absence of boundaries and edge effects. We also give the Landau gauge expression for the phase of a BEC subjected to these boundary conditions. Our spectral representation uses the eigenfunctions of the one-body Hamiltonian as basis functions. Since there is no known exact quadrature rule for these basis functions we approximately implement the projection associated with the energy cutoff, but we show that by choosing a suitably fine spatial grid the resulting error can be made negligible. We show how the convergence of this model is affected by simulation parameters such as the size of the spatial grid and the number of Landau levels. Adding dissipation, we use our method to find the lattice ground state for N vortices. We can then perturb the ground-state, to investigate the melting of the lattice.