Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Physical review. E

Uniform-density Bose-Einstein condensates of the Gross-Pitaevskii equation found by solving the inverse problem for the confining potential

Fred Cooper, Avinash Khare, John F. Dawson, E. G. Charalampidis, Avadh Saxena

Published Jun 5, 2023• 5 citations
Publisher

Abstract

In this work, we study the existence and stability of constant density (flat-top) solutions to the Gross-Pitaevskii equation (GPE) in confining potentials. These are constructed by using the "inverse problem" approach which corresponds to the identification of confining potentials that make flat-top waveforms exact solutions to the GPE. In the one-dimensional case, the exact solution is the sum of stationary kink and antikink solutions, and in the overlapping region, the density is constant. In higher spatial dimensions, the exact solutions are generalizations of this wave function. In the absence of self-interactions, the confining potential is similar to a smoothed-out finite square well with minima also at the edges. When self-interactions are added, terms proportional to ±gψ^{*}ψ and ±gM with M representing the mass or number of particles in Bose-Einstein condensates get added to the confining potential and total energy, respectively. In the realm of stability analysis, we find (linearly) stable solutions in the case with repulsive self-interactions which also are stable to self-similar deformations. For attractive interactions, however, the minima at the edges of the potential get deeper and a barrier in the center forms as we increase the norm. This leads to instabilities at a critical value of M. Comparing the stability criteria from Derrick's theorem with Bogoliubov-de Gennes (BdG) analysis stability results, we find that both predict stability for repulsive self-interactions and instability at a critical mass M for attractive interactions. However, the numerical analysis gives a much lower critical mass. This is due to the emergence of symmetry-breaking instabilities that were detected by the BdG analysis and violate the symmetry x→-x assumed by Derrick's theorem.

DOI
10.1103/physreve.107.064202

Related papers

Uniform Bose-Einstein Condensates as Kovaton solutions of the Gross-Pitaevskii Equation through a Reverse-Engineered Potential

Fred Cooper, Avinash Khare, Dawson, John F. +2 more

In this work, we consider a ``reverse-engineering'' approach to construct confining potentials that support exact, constant density kovaton solutions to the classical Gross-Pitaevskii equation (GPE) also known as the nonlinear Schrödinger equation (NLSE). In the one-dimensional case, the exact solution is the sum of stationary kink and anti-kink solutions, i.e. a kovaton, and in the overlapping region, the density is constant. In higher dimensions, the exact solutions are generalizations of this wave function. In the absence of self-interactions, the confining potential is similar to a smoothed out finite square well with minima also at the edges. When self-interactions are added, a term proportional to $\pm g ψ^{\ast}ψ$ gets added to the confining potential and $\pm g M$, where $M$ is the norm, gets added to the total energy. In the realm of stability analysis, we find (linearly) stable solutions in the case with repulsive self-interactions which also are stable to self-similar deformations. For attractive interactions, however, the minima at the edges of the potential get deeper and a barrier in the center forms as we increase the norm. This leads to instabilities at a critical value of $M$ (related to the number of particles in the BEC). Comparing the stability criteria from Derrick's theorem and Bogoliubov-de Gennes analysis stability results, we find that both predict stability for repulsive self-interactions and instability at a critical mass $M$ for attractive interactions. However, the numerical analysis gives a much lower critical mass. The numerical analysis shows further that the initial instabilities violate the symmetry $x\rightarrow-x$ assumed by Derrick's theorem.

arXiv:2303.02275•Mar 4, 2023•0 citations

The application of the “inverse problem” method for constructing confining potentials that make N-soliton waveforms exact solutions in the Gross–Pitaevskii equation

Fred Cooper, Avinash Khare, John F. Dawson +2 more

In this work, we discuss an application of the "inverse problem" method to find the external trapping potential, which has particular N trapped soliton-like solutions of the Gross-Pitaevskii equation (GPE) also known as the cubic nonlinear Schrödinger equation (NLSE). This inverse method assumes particular forms for the trapped soliton wave function, which then determines the (unique) external (confining) potential. The latter renders these assumed waveforms exact solutions of the GPE (NLSE) for both attractive (g<0) and repulsive (g>0) self-interactions. For both signs of g, we discuss the stability with respect to self-similar deformations and translations. For g<0, a critical mass Mc or equivalently the number of particles for instabilities to arise can often be found analytically. On the other hand, for the case with g>0 corresponding to repulsive self-interactions which is often discussed in the atomic physics realm of Bose-Einstein condensates, the bound solutions are found to be always stable. For g<0, we also determine the critical mass numerically by using linear stability or Bogoliubov-de Gennes analysis, and compare these results with our analytic estimates. Various analytic forms for the trapped N-soliton solutions in one, two, and three spatial dimensions are discussed, including sums of Gaussians or higher-order eigenfunctions of the harmonic oscillator Hamiltonian.

Chaos An Interdisciplinary Journal of Nonlinear Science•Apr 1, 2024•2 citations

An analytical solution for the critical number of particles for stable bose-einstein condensation under the influence of an anisotropic potential

Alvaro S. de Sousa, Vanderlei Salvador Bagnato, Albérico B. F. da Silva

We have considered a Bose gas in an anisotropic potential. Applying the the Gross-Pitaevskii Equation (GPE) for a confined dilute atomic gas, we have used the methods of optimized perturbation theory and self-similar root approximants, to obtain an analytical formula for the critical number of particles as a function of the anisotropy parameter for the potential. The spectrum of the GPE is also discussed.

Brazilian Journal of Physics•Mar 1, 2008•6 citations

Theoretical Studies of dilute Bose-Einstein condensates in a double-well potential

E. Sakellari

In this Thesis we apply the Gross-Pitaevskii equation (GPE) to describe properties of a dilute, near zero temperature Bose gas for various confining geometries. We start by reviewing some basic information about the density, the chemical potential and elementary excitations of a dilute atomic condensate confined in a single harmonic trap for a Bose condensate with repulsive and attractive interactions and we also discuss the stability in the case of attractive interactions. We extend our study to a one and three dimensional double-well trap. We investigate the eigenenergy levels and show that the nonlinearity leads to triangular structures which appear either in the ground or excited states for the case of a Bose condensate with attractive or repulsive interactions respectively. We apply the eigenenergy level picture to analyse Josephson effects induced when the barrier IS moved at a constant velocity across the trapping potential or by the application of a time-dependent potential gradient. The GPE simulations are compared to the predictions of a nonlinear two state model. Above a critical velocity there is a transition to a superposition of ground and excited states which leads to sudden changes in the population difference. The direction of Josephson flow depends critically on the initial state of the system and we discuss the feasibility of experimental control of the atomic flow using phase-imprinting. The stability of a low temperature Bose-Einstein condensate with attract interactions in one and three dimensional double-well potentiate is discussed. The condensate is shown to collapse at a critical potential gradient which corresponds to a critical number of atoms in one of the two wells. Finally we investigate the stability and tunnelling effects in a multi-well system.

•Jan 1, 2004•0 citations

Similarity solutions and collapse in the attractive Gross-Pitaevskii equation

Andrei V. Rybin, G. G. Varzugin, Markus Lindberg +2 more

We analyze a generalized Gross-Pitaevskii (GP) equation involving a paraboloidal trap potential in D space dimensions and generalized to a nonlinearity of order 2n+1. For attractive coupling constants collapse of the particle density occurs for Dn>/=2 and typically to a delta function centered at the origin of the trap. By introducing a special variable for the spherically symmetric solutions, we show that all such solutions are self-similar close to the center of the trap. Exact self-similar solutions occur if, and only if, Dn=2, and for this case of Dn=2 we exhibit an exact but rather special D=1 analytical self-similar solution collapsing to a delta function which, however, recovers and collapses periodically, while the ordinary GP equation in two space dimensions also has a special solution with periodic delta function collapses and revivals of the density. The relevance of these various results to attractive Bose-Einstein condensation in spherically symmetric traps is discussed.

arXiv:cond-mat/0001059•Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Nov 1, 2000•32 citations

Reentrant stability of Bose-Einstein-condensate standing-wave patterns

Ryan M. Kalas, Dmitry Solenov, Eddy Timmermans

We describe standing-wave patterns induced by an attractive finite-ranged external potential inside a large Bose-Einstein Condensate (BEC). As the potential depth increases, the time-independent Gross-Pitaevskii equation develops pairs of solutions that have nodes in their wave function. We elucidate the nature of these states and study their dynamical stability. Although we study the problem in a two-dimensional BEC subject to a cylindrically symmetric square well potential of a radius that is comparable to the coherence length of the BEC, our analysis reveals general trends that are valid in two and three dimensions, independent of the symmetry of the localized potential well, and suggestive of the behavior in general short-range and large-range potentials. One set of nodal BEC wave functions resembles the single particle $n$-node bound-state wave function of the potential well, the other wave functions resemble the $n\ensuremath{-}1$-node bound-state wave function with a kink state pinned by the potential. The second state, though corresponding to the lower free energy value of the pair of $n$-node BEC states, is always unstable, whereas the first can be dynamically stable in intervals of the potential well depth, implying that the standing-wave BEC can evolve from being dynamically unstable to stable and back to unstable as the potential well is adiabatically deepened---a phenomenon that we refer to as ``reentrant dynamical stability.''

arXiv:0910.2711•Physical Review A•May 20, 2010•1 citation
In this work, we study the existence and stability of constant density (flat-top) solutions to the Gross-Pitaevskii equation (GPE) in confining potentials. These are constructed by using the "inverse problem" approach which corresponds to the identification of confining potentials that make flat-top waveforms exact solutions to the GPE. In the one-dimensional case, the exact solution is the sum of stationary kink and antikink solutions, and in the overlapping region, the density is constant. In higher spatial dimensions, the exact solutions are generalizations of this wave function. In the absence of self-interactions, the confining potential is similar to a smoothed-out finite square well with minima also at the edges. When self-interactions are added, terms proportional to ±gψ^{*}ψ and ±gM with M representing the mass or number of particles in Bose-Einstein condensates get added to the confining potential and total energy, respectively. In the realm of stability analysis, we find (linearly) stable solutions in the case with repulsive self-interactions which also are stable to self-similar deformations. For attractive interactions, however, the minima at the edges of the potential get deeper and a barrier in the center forms as we increase the norm. This leads to instabilities at a critical value of M. Comparing the stability criteria from Derrick's theorem with Bogoliubov-de Gennes (BdG) analysis stability results, we find that both predict stability for repulsive self-interactions and instability at a critical mass M for attractive interactions. However, the numerical analysis gives a much lower critical mass. This is due to the emergence of symmetry-breaking instabilities that were detected by the BdG analysis and violate the symmetry x→-x assumed by Derrick's theorem.
In this work, we consider a ``reverse-engineering'' approach to construct confining potentials that support exact, constant density kovaton solutions to the classical Gross-Pitaevskii equation (GPE) also known as the nonlinear Schrödinger equation (NLSE). In the one-dimensional case, the exact solution is the sum of stationary kink and anti-kink solutions, i.e. a kovaton, and in the overlapping region, the density is constant. In higher dimensions, the exact solutions are generalizations of this wave function. In the absence of self-interactions, the confining potential is similar to a smoothed out finite square well with minima also at the edges. When self-interactions are added, a term proportional to ±gψ∗ψ\pm g ψ^{\ast}ψ±gψ∗ψ gets added to the confining potential and ±gM\pm g M±gM, where MMM is the norm, gets added to the total energy. In the realm of stability analysis, we find (linearly) stable solutions in the case with repulsive self-interactions which also are stable to self-similar deformations. For attractive interactions, however, the minima at the edges of the potential get deeper and a barrier in the center forms as we increase the norm. This leads to instabilities at a critical value of MMM (related to the number of particles in the BEC). Comparing the stability criteria from Derrick's theorem and Bogoliubov-de Gennes analysis stability results, we find that both predict stability for repulsive self-interactions and instability at a critical mass MMM for attractive interactions. However, the numerical analysis gives a much lower critical mass. The numerical analysis shows further that the initial instabilities violate the symmetry x→−xx\rightarrow-xx→−x assumed by Derrick's theorem.
In this work, we discuss an application of the "inverse problem" method to find the external trapping potential, which has particular N trapped soliton-like solutions of the Gross-Pitaevskii equation (GPE) also known as the cubic nonlinear Schrödinger equation (NLSE). This inverse method assumes particular forms for the trapped soliton wave function, which then determines the (unique) external (confining) potential. The latter renders these assumed waveforms exact solutions of the GPE (NLSE) for both attractive (g<0) and repulsive (g>0) self-interactions. For both signs of g, we discuss the stability with respect to self-similar deformations and translations. For g<0, a critical mass Mc or equivalently the number of particles for instabilities to arise can often be found analytically. On the other hand, for the case with g>0 corresponding to repulsive self-interactions which is often discussed in the atomic physics realm of Bose-Einstein condensates, the bound solutions are found to be always stable. For g<0, we also determine the critical mass numerically by using linear stability or Bogoliubov-de Gennes analysis, and compare these results with our analytic estimates. Various analytic forms for the trapped N-soliton solutions in one, two, and three spatial dimensions are discussed, including sums of Gaussians or higher-order eigenfunctions of the harmonic oscillator Hamiltonian.
In this Thesis we apply the Gross-Pitaevskii equation (GPE) to describe properties of a dilute, near zero temperature Bose gas for various confining geometries. We start by reviewing some basic information about the density, the chemical potential and elementary excitations of a dilute atomic condensate confined in a single harmonic trap for a Bose condensate with repulsive and attractive interactions and we also discuss the stability in the case of attractive interactions. We extend our study to a one and three dimensional double-well trap. We investigate the eigenenergy levels and show that the nonlinearity leads to triangular structures which appear either in the ground or excited states for the case of a Bose condensate with attractive or repulsive interactions respectively. We apply the eigenenergy level picture to analyse Josephson effects induced when the barrier IS moved at a constant velocity across the trapping potential or by the application of a time-dependent potential gradient. The GPE simulations are compared to the predictions of a nonlinear two state model. Above a critical velocity there is a transition to a superposition of ground and excited states which leads to sudden changes in the population difference. The direction of Josephson flow depends critically on the initial state of the system and we discuss the feasibility of experimental control of the atomic flow using phase-imprinting. The stability of a low temperature Bose-Einstein condensate with attract interactions in one and three dimensional double-well potentiate is discussed. The condensate is shown to collapse at a critical potential gradient which corresponds to a critical number of atoms in one of the two wells. Finally we investigate the stability and tunnelling effects in a multi-well system.
We analyze a generalized Gross-Pitaevskii (GP) equation involving a paraboloidal trap potential in D space dimensions and generalized to a nonlinearity of order 2n+1. For attractive coupling constants collapse of the particle density occurs for Dn>/=2 and typically to a delta function centered at the origin of the trap. By introducing a special variable for the spherically symmetric solutions, we show that all such solutions are self-similar close to the center of the trap. Exact self-similar solutions occur if, and only if, Dn=2, and for this case of Dn=2 we exhibit an exact but rather special D=1 analytical self-similar solution collapsing to a delta function which, however, recovers and collapses periodically, while the ordinary GP equation in two space dimensions also has a special solution with periodic delta function collapses and revivals of the density. The relevance of these various results to attractive Bose-Einstein condensation in spherically symmetric traps is discussed.
We describe standing-wave patterns induced by an attractive finite-ranged external potential inside a large Bose-Einstein Condensate (BEC). As the potential depth increases, the time-independent Gross-Pitaevskii equation develops pairs of solutions that have nodes in their wave function. We elucidate the nature of these states and study their dynamical stability. Although we study the problem in a two-dimensional BEC subject to a cylindrically symmetric square well potential of a radius that is comparable to the coherence length of the BEC, our analysis reveals general trends that are valid in two and three dimensions, independent of the symmetry of the localized potential well, and suggestive of the behavior in general short-range and large-range potentials. One set of nodal BEC wave functions resembles the single particle nnn-node bound-state wave function of the potential well, the other wave functions resemble the $n\ensuremath{-}1$-node bound-state wave function with a kink state pinned by the potential. The second state, though corresponding to the lower free energy value of the pair of nnn-node BEC states, is always unstable, whereas the first can be dynamically stable in intervals of the potential well depth, implying that the standing-wave BEC can evolve from being dynamically unstable to stable and back to unstable as the potential well is adiabatically deepened---a phenomenon that we refer to as ``reentrant dynamical stability.''