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

Applied Mathematics

A Spectral Integral Equation Solution of the Gross-Pitaevskii Equation

George H. Rawitscher

Published Jan 1, 2013• 1 citation
Publisher

Abstract

The Gross-Pitaevskii equation (GPE), that describes the wave function of a number of coherent Bose particles contained in a trap, contains the cube of the normalized wave function, times a factor proportional to the number of coherent atoms. The square of the wave function, times the above mentioned factor, is defined as the Hartree potential. A method implemented here for the numerical solution of the GPE consists in obtaining the Hartree potential iteratively, starting with the Thomas Fermi approximation to this potential. The energy eigenvalues and the corresponding wave functions for each successive potential are obtained by a spectral method described previously. After approximately 35 iterations a stability of eight significant figures for the energy eigenvalues is obtained. This method has the advantage of being physically intuitive, and could be extended to the calculation of a shell-model potential in nuclear physics, once the Pauli exclusion principle is allowed for.

DOI
10.4236/am.2013.410a3009

Related papers

Solution of the Gross-Pitaevskii equation in terms of the associated non-linear Hartree potential

George H. Rawitscher

The Gross-Pitaevskii equation (GP), that describes the wave function of a number of coherent Bose particles contained in a trap, contains the cube of the normalized wave function, times a factor proportional to the number of coherent atoms. The square of the wave function, times the above mentioned factor, is defined as the Hartree potential. A method implemented here for the numerical solution of the GP equation consists in obtaining the Hartree potential iteratively, starting with the Thomas Fermi approximation to this potential. The energy eigenvalues and the corresponding wave functions for each successive potential are obtained by a method described previously. After approximately 35 iterations a stability of eight significant figures for the energy eigenvalues is obtained.This method has the advantage of being physically intuitive, and could be extended to the calculation of a shell-model potential in nuclear physics, once the Pauli exclusion principle is allowed for.

arXiv:1307.4031•Jul 15, 2013•0 citations

Solution of the ground state wave function of Bose-condensed gas in a harmonic trap based on the Gross-Pitaevskii function

Xu Zhi-Jun, Shi Jian-Qing, Li Zhen +2 more

Based on the Gross-Pitaevskii(G-P) energy functional and variational method, we present a new method to find the ground state wave function for Bose-condensed gas in a harmonic trap at zero temperature. With this method we are able to find the analytic expression of ground state wave function and explore the relevant quantities, such as energy, chemical potential and so on. These results agree well with previous ground state numerical solutions of the G-P equation given by Dalfovo et al. Under the large number of particles N limit for repulsive interaction, our results and Thomas-Fermi approximation model approach to reaching unanimity. This new method is simple compared to other methods used to solve numerically the G-P equation, and one can obtained the analytic and reliable results.

Acta Physica Sinica•Jan 1, 2006•1 citation

Simple model for interactions and corrections to the Gross-Pitaevskii equation

Hagar Veksler, Shmuel Fishman, Wolfgang Ketterle

One of the assumptions leading to the Gross-Pitaevskii equation (GPE) is that the interaction between atom pairs can be written effectively as a $\ensuremath{\delta}$ function so that the interaction range of the particles is assumed to vanish. A simple model that takes into account the extension of the interparticle potential is introduced. The correction to the GPE predictions for the energy of a condensate confined by a harmonic trap in the Thomas-Fermi regime is estimated. Although it is found to be small, we believe that in some situations it can be measured using its dependance on the frequency of the confining trap. Due to the simplicity of the model, it may have a wide range of applications.

arXiv:1404.2586•Physical Review A•Aug 13, 2014•36 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

Hartree-Fock calculation for trapped, two-dimensional interacting Bose gases of finite size

Jae Gil Kim, Kab Seok Kang, Bong Soo Kim +1 more

The mean-field properties of trapped, two-dimensional interacting Bose gases of finite size are surveyed numerically. The coupled nonlinear equations of the Gross-Pitaevskii (GP) equation and the Hartree-Fock equation for the condensate and the thermal component, respectively, are calculated by the hybrid method, in which the GP equation and the low-lying excitations less than cut-off energy are calculated quantum mechanically using the basis representation, while the high-lying states are obtained within the local density approximation. Consistent solutions are obtained for the whole range of temperatures including the critical temperature by determining the thermodynamic chemical potential µ from the condition that the total atomic number must be the desired value, instead of identifying µ with the eigenvalue of the GP equation. Comparisons with the results of the semiclassical two-fluid model are presented.

Journal of Physics B Atomic Molecular and Optical Physics•Jun 22, 2000•7 citations

Green function method and the Gross-Pitaevskii equation

Carlos A Trallero-Herrero, Victor Lopez‐Richard, Julio C. Drake-Pérez +1 more

"The time-independent one-dimensional Gross-Pitaevskii (GP) equation for a Bose-Einstein condensate is solved by using Green functionformalism atT=0 K. From the integral GP equation, the analytical solution for the order parameter,©(x), and the chemical potential,π, has been formally obtained as functions of a single dimensionless parameter. The obtained solutions are valid in the whole range ofnegative and positive values of the non-linear interaction term that is used to characterize the condensates of alkali atoms. A function hasbeen introduced in order to describe the accumulated error that allows a comparative analysis with respect to known methods of solutionsuch as the Thomas-Fermi approximation and the “bright” soliton. This study sets universal ranges of parameters where analytical solutionscan be implemented."

Revista Mexicana de Física•Jan 1, 2007•0 citations
The Gross-Pitaevskii equation (GPE), that describes the wave function of a number of coherent Bose particles contained in a trap, contains the cube of the normalized wave function, times a factor proportional to the number of coherent atoms. The square of the wave function, times the above mentioned factor, is defined as the Hartree potential. A method implemented here for the numerical solution of the GPE consists in obtaining the Hartree potential iteratively, starting with the Thomas Fermi approximation to this potential. The energy eigenvalues and the corresponding wave functions for each successive potential are obtained by a spectral method described previously. After approximately 35 iterations a stability of eight significant figures for the energy eigenvalues is obtained. This method has the advantage of being physically intuitive, and could be extended to the calculation of a shell-model potential in nuclear physics, once the Pauli exclusion principle is allowed for.
The Gross-Pitaevskii equation (GP), that describes the wave function of a number of coherent Bose particles contained in a trap, contains the cube of the normalized wave function, times a factor proportional to the number of coherent atoms. The square of the wave function, times the above mentioned factor, is defined as the Hartree potential. A method implemented here for the numerical solution of the GP equation consists in obtaining the Hartree potential iteratively, starting with the Thomas Fermi approximation to this potential. The energy eigenvalues and the corresponding wave functions for each successive potential are obtained by a method described previously. After approximately 35 iterations a stability of eight significant figures for the energy eigenvalues is obtained.This method has the advantage of being physically intuitive, and could be extended to the calculation of a shell-model potential in nuclear physics, once the Pauli exclusion principle is allowed for.
Based on the Gross-Pitaevskii(G-P) energy functional and variational method, we present a new method to find the ground state wave function for Bose-condensed gas in a harmonic trap at zero temperature. With this method we are able to find the analytic expression of ground state wave function and explore the relevant quantities, such as energy, chemical potential and so on. These results agree well with previous ground state numerical solutions of the G-P equation given by Dalfovo et al. Under the large number of particles N limit for repulsive interaction, our results and Thomas-Fermi approximation model approach to reaching unanimity. This new method is simple compared to other methods used to solve numerically the G-P equation, and one can obtained the analytic and reliable results.
One of the assumptions leading to the Gross-Pitaevskii equation (GPE) is that the interaction between atom pairs can be written effectively as a $\ensuremath{\delta}$ function so that the interaction range of the particles is assumed to vanish. A simple model that takes into account the extension of the interparticle potential is introduced. The correction to the GPE predictions for the energy of a condensate confined by a harmonic trap in the Thomas-Fermi regime is estimated. Although it is found to be small, we believe that in some situations it can be measured using its dependance on the frequency of the confining trap. Due to the simplicity of the model, it may have a wide range of applications.
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.
The mean-field properties of trapped, two-dimensional interacting Bose gases of finite size are surveyed numerically. The coupled nonlinear equations of the Gross-Pitaevskii (GP) equation and the Hartree-Fock equation for the condensate and the thermal component, respectively, are calculated by the hybrid method, in which the GP equation and the low-lying excitations less than cut-off energy are calculated quantum mechanically using the basis representation, while the high-lying states are obtained within the local density approximation. Consistent solutions are obtained for the whole range of temperatures including the critical temperature by determining the thermodynamic chemical potential µ from the condition that the total atomic number must be the desired value, instead of identifying µ with the eigenvalue of the GP equation. Comparisons with the results of the semiclassical two-fluid model are presented.
"The time-independent one-dimensional Gross-Pitaevskii (GP) equation for a Bose-Einstein condensate is solved by using Green functionformalism atT=0 K. From the integral GP equation, the analytical solution for the order parameter,©(x), and the chemical potential,π, has been formally obtained as functions of a single dimensionless parameter. The obtained solutions are valid in the whole range ofnegative and positive values of the non-linear interaction term that is used to characterize the condensates of alkali atoms. A function hasbeen introduced in order to describe the accumulated error that allows a comparative analysis with respect to known methods of solutionsuch as the Thomas-Fermi approximation and the “bright” soliton. This study sets universal ranges of parameters where analytical solutionscan be implemented."