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

Chaos An Interdisciplinary Journal of Nonlinear Science

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, E. G. Charalampidis, Avadh Saxena

Published Apr 1, 2024• 2 citations
Publisher

Abstract

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.

DOI
10.1063/5.0183167

Related papers

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 +2 more

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.

Physical review. E•Jun 5, 2023•5 citations

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

Soliton dynamics for trapped Bose-Einstein condensate with higher-order interaction

Ying Wang, Shuyu Zhou

We model the three-dimensional (3D) trapped Bose-Einstein condensate with the Gross-Pitaevskii equation (GPE) incorporating higher-order nonlinear interaction effects. Based on the F-expansion method and self-similar approach, we analytically derive the exact soliton solutions of the 3D GPE, identifying the oscillatory mode under certain experimental settings. The oscillatory kinetic parameter values we obtained coincide very well with those reported in previous work, which is based on the variational approach, showing the applicability of the theoretical treatment presented in this work.

AIP Advances•Aug 1, 2017•11 citations

Gap solitons in elongated geometries: The one-dimensional Gross-Pitaevskii equation and beyond

Antonio Muñoz Mateo, V. Delgado, Boris A. Malomed

We report results of a systematic analysis of matter-wave gap solitons (GSs) in three-dimensional self-repulsive Bose-Einstein condensates (BECs) loaded into a combination of a cigar-shaped trap and axial optical-lattice (OL) potential. Basic cases of the strong, intermediate, and weak radial (transverse) confinement are considered, as well as settings with shallow and deep OL potentials. Only in the case of the shallow lattice combined with tight radial confinement, which actually has little relevance to realistic experimental conditions, does the usual one-dimensional (1D) cubic Gross-Pitaevskii equation (GPE) furnish a sufficiently accurate description of GSs. However, the effective 1D equation with the nonpolynomial nonlinearity, derived in Ref. [Phys. Rev. A 77, 013617 (2008)], provides for quite an accurate approximation for the GSs in all cases, including the situation with weak transverse confinement, when the soliton's shape includes a considerable contribution from higher-order transverse modes, in addition to the usual ground-state wave function of the respective harmonic oscillator. Both fundamental GSs and their multipeak bound states are considered. The stability is analyzed by means of systematic simulations. It is concluded that almost all the fundamental GSs are stable, while their bound states may be stable if the underlying OL potential is deep enough.

arXiv:1104.0962•Physical Review A•May 9, 2011•18 citations

Exact exact solutions of the Gross-Pitaevskii equation for stable vortex modes

Lei Wu, Lu Li, Jie-Fang Zhang +3 more

We construct exact solutions of the Gross-Pitaevskii equation for solitary vortices, and approximate ones for fundamental solitons, in 2D models of Bose-Einstein condensates with a spatially modulated nonlinearity of either sign and a harmonic trapping potential. The number of vortex-soliton (VS) modes is determined by the discrete energy spectrum of a related linear Schrödinger equation. The VS families in the system with the attractive and repulsive nonlinearity are mutually complementary. \emph{% Stable} VSs with vorticity $S\geq 2$ and those corresponding to higher-order radial states are reported for the first time, in the case of the attraction and repulsion, respectively.

arXiv:1006.3695•Jun 18, 2010•1 citation

THREE-BODY INTERACTIONS BEYOND THE GROSS–PITAEVSKII EQUATION AND MODULATIONAL INSTABILITY OF BOSE–EINSTEIN CONDENSATES

D. Belobo Belobo, G. H. Ben-Bolie, Timoléon Crépin Kofané

Beyond the mean-field theory, a new model of the Gross–Pitaevskii equation (GPE) that describes the dynamics of Bose–Einstein condensates (BECs) is derived using an appropriate phase-imprint on the old wavefunction. This modified version of the GPE in addition to the two-body interactions term, also takes into account effects of the three-body interactions. The three-body interactions consist of a quintic term and the delayed nonlinear response of the condensate system term. Then, the modulational instability (MI) of the new GPE confined in an attractive harmonic potential is investigated. The analytical study shows that the three-body interactions destabilize more the condensate system while the external potential alleviates the instability. Numerical results confirm the theoretical predictions. Further numerical investigations of the behavior of solitons reveal that the three-body interactions enhance the appearance of solitons, increase the number of solitons generated and deeply change the lifetime of solitons. Moreover, the external potential delays the appearance of solitons. Besides, a new initial condition is introduced which enables to increase the number of solitons created and deeply affects the trail of chains of solitons generated. Moreover, the MI of a condensate without the external potential, and in a repulsive potential is also investigated.

International Journal of Modern Physics B•Nov 8, 2012•4 citations
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 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.
We model the three-dimensional (3D) trapped Bose-Einstein condensate with the Gross-Pitaevskii equation (GPE) incorporating higher-order nonlinear interaction effects. Based on the F-expansion method and self-similar approach, we analytically derive the exact soliton solutions of the 3D GPE, identifying the oscillatory mode under certain experimental settings. The oscillatory kinetic parameter values we obtained coincide very well with those reported in previous work, which is based on the variational approach, showing the applicability of the theoretical treatment presented in this work.
We report results of a systematic analysis of matter-wave gap solitons (GSs) in three-dimensional self-repulsive Bose-Einstein condensates (BECs) loaded into a combination of a cigar-shaped trap and axial optical-lattice (OL) potential. Basic cases of the strong, intermediate, and weak radial (transverse) confinement are considered, as well as settings with shallow and deep OL potentials. Only in the case of the shallow lattice combined with tight radial confinement, which actually has little relevance to realistic experimental conditions, does the usual one-dimensional (1D) cubic Gross-Pitaevskii equation (GPE) furnish a sufficiently accurate description of GSs. However, the effective 1D equation with the nonpolynomial nonlinearity, derived in Ref. [Phys. Rev. A 77, 013617 (2008)], provides for quite an accurate approximation for the GSs in all cases, including the situation with weak transverse confinement, when the soliton's shape includes a considerable contribution from higher-order transverse modes, in addition to the usual ground-state wave function of the respective harmonic oscillator. Both fundamental GSs and their multipeak bound states are considered. The stability is analyzed by means of systematic simulations. It is concluded that almost all the fundamental GSs are stable, while their bound states may be stable if the underlying OL potential is deep enough.
We construct exact solutions of the Gross-Pitaevskii equation for solitary vortices, and approximate ones for fundamental solitons, in 2D models of Bose-Einstein condensates with a spatially modulated nonlinearity of either sign and a harmonic trapping potential. The number of vortex-soliton (VS) modes is determined by the discrete energy spectrum of a related linear Schrödinger equation. The VS families in the system with the attractive and repulsive nonlinearity are mutually complementary. \emph{% Stable} VSs with vorticity S≥2S\geq 2S≥2 and those corresponding to higher-order radial states are reported for the first time, in the case of the attraction and repulsion, respectively.
Beyond the mean-field theory, a new model of the Gross–Pitaevskii equation (GPE) that describes the dynamics of Bose–Einstein condensates (BECs) is derived using an appropriate phase-imprint on the old wavefunction. This modified version of the GPE in addition to the two-body interactions term, also takes into account effects of the three-body interactions. The three-body interactions consist of a quintic term and the delayed nonlinear response of the condensate system term. Then, the modulational instability (MI) of the new GPE confined in an attractive harmonic potential is investigated. The analytical study shows that the three-body interactions destabilize more the condensate system while the external potential alleviates the instability. Numerical results confirm the theoretical predictions. Further numerical investigations of the behavior of solitons reveal that the three-body interactions enhance the appearance of solitons, increase the number of solitons generated and deeply change the lifetime of solitons. Moreover, the external potential delays the appearance of solitons. Besides, a new initial condition is introduced which enables to increase the number of solitons created and deeply affects the trail of chains of solitons generated. Moreover, the MI of a condensate without the external potential, and in a repulsive potential is also investigated.