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

International Journal of Modern Physics B

Generalization of the Gross–Pitaevskii equation for multi-particle cases

Nicolai Bogoliubov, M. Yu. Rasulova, Tohir Vohidovich Akramov

Published Nov 3, 2021• 2 citations
Publisher

Abstract

A new method is proposed to obtain Gross–Pitaevskii equation by the chain of Bogoliubov–Born–Green–Kirkwood–Yvon (BBGKY) quantum kinetic equations. In that sense, we investigate the dynamics of a quantum system including infinite number of identical particles which interact via a (special) pair potential on the form of Dirac delta-function.

DOI
10.1142/s0217979221503112

Related papers

Evolution of a quantum system of many particles interacting via the generalized Yukawa potential

Nikolai Nikolaevich Bogoliubov, M. Yu. Rasulova, Umarbek Avazov

We study the evolution of a system of N particles that have identical masses and charges and interact via the generalized Yukawa potential. The system is placed in a bounded region. The evolution of such a system is described by the Bogoliubov–Born–Green–Kirkwood–Yvon (BBGKY) chain of quantum kinetic equations. Using semigroup theory, we prove the existence of a unique solution of the BBGKY chain of quantum kinetic equations with the generalized Yukawa potential.

Theoretical and Mathematical Physics•Dec 1, 2016•2 citations

Particle-number-conserving Bogoliubov method which demonstrates the validity of the time-dependent Gross-Pitaevskii equation for a highly condensed Bose gas

Crispin W. Gardiner

The Bogoliubov method for the excitation spectrum of a Bose-condensed gas is generalized to apply to a gas with an exact large number $N$ of particles. This generalization yields a description of the Schr\"odinger picture field operators as the product of an annihilation operator $A$ for the total number of particles and the sum of a ``condensate wave function'' $\ensuremath{\xi}(\mathbf{x})$ and a phonon field operator $\ensuremath{\chi}(\mathbf{x})$ in the form $\ensuremath{\psi}(\mathbf{x})\ensuremath{\approx}A{\ensuremath{\xi}(\mathbf{x})+\ensuremath{\chi}(\mathbf{x})/\sqrt{N}}$ when the field operator acts on the $N$ particle subspace. It is then possible to expand the Hamiltonian in decreasing powers of $\sqrt{N}$, and thus obtain solutions for eigenvalues and eigenstates as an asymptotic expansion of the same kind. It is also possible to compute all matrix elements of field operators between states of different $N$. The excitation spectrum can be obtained by essentially the same method as Bogoliubov only if $\ensuremath{\xi}(\mathbf{x})$ is a solution of the time-independent Gross-Pitaevskii equation for $N$ particles and any chemical potential $\ensuremath{\mu}$, which yields a valid and stable solution of the Gross-Pitaevskii equation. The treatment within a subspace of fixed $N$ is identical in form to that usually used, but the interpretation of the operators is slightly different. A time-dependent generalization is then made, yielding an asymptotic expansion in decreasing powers of $\sqrt{N}$ for the equations of motion. In this expansion the condensate wave function has the time-dependent form $\ensuremath{\xi}(\mathbf{x},t)$, and the condition for the validity of the expansion is that $\ensuremath{\xi}(\mathbf{x},t)$ satisfies the time-dependent Gross-Pitaevskii equation $\ensuremath{\partial}\ensuremath{\xi}/\ensuremath{\partial}t=\ensuremath{-}({\ensuremath{\Elzxh}}^{2}/2m){\ensuremath{\nabla}}^{2}\ensuremath{\xi}+V\ensuremath{\xi}+Nu|\ensuremath{\xi}{|}^{2}\ensuremath{\xi}$. The physics is then described in a kind of interaction picture, called the condensate picture, in which the phonon operator can be expressed as $\ensuremath{\chi}(\mathbf{x},t)={\ensuremath{\sum}}_{k}{\ensuremath{\xi}}_{k}(\mathbf{x},t){\ensuremath{\alpha}}_{k},$ where the operators ${\ensuremath{\alpha}}_{k}$ are time-independent annihilation operators, and the state vector has a time evolution described by a Schr\"odinger equation in which the Hamiltonian is a time-dependent quadratic form in the phonon creation and annihilation operators, whose coefficients are explicitly determined in terms of the time-dependent condensate wave function $\ensuremath{\xi}(\mathbf{x},t)$.

arXiv:quant-ph/9703005•Physical Review A•Aug 1, 1997•256 citations

Correlation Structures, Many-Body Scattering Processes, and the Derivation of the Gross–Pitaevskii Hierarchy

Xuwen Chen, Justin Holmer

We consider the dynamics of |$N$| bosons in three dimensions. We assume that the pair interaction is given by |$N^{3\beta -1}V(N^{\beta }\cdot )$|⁠. By studying an associated many-body wave operator, we introduce a Bogoliubov–Born–Green–Kirkwood–Yvon (BBGKY) hierarchy which takes into account all of the interparticle singular correlation structures developed by the many-body evolution from the beginning. Assuming energy conditions on the |$N$|-body wave function, for |$\beta \in ( 0,1] $|⁠, we derive the Gross–Pitaevskii hierarchy with 2-body interaction. In particular, we establish that, in the |$N\rightarrow \infty $| limit, all |$k$|-body scattering processes vanish if |$k\geqslant 3$| and thus provide a direct answer to a question raised by Erdös et al. [30]. Moreover, this new BBGKY hierarchy shares the limit points with the ordinary BBGKY hierarchy strongly for |$\beta \in ( 0,1) $| and weakly for |$\beta =1$|⁠. Since this new BBGKY hierarchy converts the problem from a two-body estimate to a weaker three-body estimate for which we have the estimates to achieve |$\beta <1$|⁠, it then allows us to prove that all limit points of the ordinary BBGKY hierarchy satisfy the space–time bound conjectured by Klainerman and Machedon [48] for |$\beta \in (0,1)$|⁠.

arXiv:1409.1425•International Mathematics Research Notices•Aug 5, 2015•56 citations

Infinite Statistics and the Gross Pitaevskii Equation

Dewi Yustikasari, Mirza Satriawan

We clarify that an ideal gas obeying infinite statistics cannot undergo condensation. Then we derive the dynamic equation for an identical particle system obeying infinite statistics under external potential and inter-particle interaction. The derivation utilizes the Hamiltonian written in terms of the number operators and the transition number operators. At a very low temperature, where one can discard the dynamics of the excited occupation level, the dynamic of an infinite statistics system can be described by the Gross Pitaevskii equation, similar to the Bose-Einstein case.

arXiv:2210.09514•Oct 15, 2022•0 citations

Rigorous Derivation of the Gross-Pitaevskii Equation: the Case of a Strong Potential

László Erdős, Benjamin Schlein, Horng‐Tzer Yau

Consider a system of $N$ bosons in three dimensions interacting via a repulsive short range pair potential $N^2V(N(x_i-x_j))$, where $\bx=(x_1, >..., x_N)$ denotes the positions of the particles. Let $H_N$ denote the Hamiltonian of the system and let $\psi_{N,t}$ be the solution to the Schrodinger equation. Suppose that the initial data $\psi_{N,0}$ satisfies the energy condition \[ \leq C N >. \] and that the one-particle density matrix converges to a projection as $N \to \infty$. Then, we prove that the $k$-particle density matrices of $\psi_{N,t}$ factorize in the limit $N \to \infty$. Moreover, the one particle orbital wave function solves the time-dependent Gross-Pitaevskii equation, a cubic non-linear Schrodinger equation with the coupling constant proportional to the scattering length of the potential $V$. In \cite{ESY}, we proved the same statement under the condition that the interaction potential $V$ is sufficiently small; in the present work we extend this result to large potentials.

•Feb 26, 2008•2 citations

Derivation of the Gross-Pitaevskii equation for the dynamics of Bose-Einstein condensate

László Erdős, Benjamin Schlein, Horng‐Tzer Yau

Consider a system of $N$ bosons in three dimensions interacting via a repulsive short range pair potential $N^2V(N(x_i-x_j))$, where $\bx=(x_1, >..., x_N)$ denotes the positions of the particles. Let $H_N$ denote the Hamiltonian of the system and let $ψ_{N,t}$ be the solution to the Schrödinger equation. Suppose that the initial data $ψ_{N,0}$ satisfies the energy condition \[ < ψ_{N,0}, H_N^k ψ_{N,0} > \leq C^k N^k \] for $k=1,2,... $. We also assume that the $k$-particle density matrices of the initial state are asymptotically factorized as $N\to\infty$. We prove that the $k$-particle density matrices of $ψ_{N,t}$ are also asymptotically factorized and the one particle orbital wave function solves the Gross-Pitaevskii equation, a cubic non-linear Schrödinger equation with the coupling constant given by the scattering length of the potential $V$. We also prove the same conclusion if the energy condition holds only for $k=1$ but the factorization of $ψ_{N,0}$ is assumed in a stronger sense.

arXiv:math-ph/0606017•Annals of Mathematics•Jun 27, 2010•328 citations
A new method is proposed to obtain Gross–Pitaevskii equation by the chain of Bogoliubov–Born–Green–Kirkwood–Yvon (BBGKY) quantum kinetic equations. In that sense, we investigate the dynamics of a quantum system including infinite number of identical particles which interact via a (special) pair potential on the form of Dirac delta-function.
The Bogoliubov method for the excitation spectrum of a Bose-condensed gas is generalized to apply to a gas with an exact large number NNN of particles. This generalization yields a description of the Schr\"odinger picture field operators as the product of an annihilation operator AAA for the total number of particles and the sum of a ``condensate wave function'' $\ensuremath{\xi}(\mathbf{x})$ and a phonon field operator $\ensuremath{\chi}(\mathbf{x})$ in the form $\ensuremath{\psi}(\mathbf{x})\ensuremath{\approx}A{\ensuremath{\xi}(\mathbf{x})+\ensuremath{\chi}(\mathbf{x})/\sqrt{N}}$ when the field operator acts on the NNN particle subspace. It is then possible to expand the Hamiltonian in decreasing powers of N\sqrt{N}N​, and thus obtain solutions for eigenvalues and eigenstates as an asymptotic expansion of the same kind. It is also possible to compute all matrix elements of field operators between states of different NNN. The excitation spectrum can be obtained by essentially the same method as Bogoliubov only if $\ensuremath{\xi}(\mathbf{x})$ is a solution of the time-independent Gross-Pitaevskii equation for NNN particles and any chemical potential $\ensuremath{\mu}$, which yields a valid and stable solution of the Gross-Pitaevskii equation. The treatment within a subspace of fixed NNN is identical in form to that usually used, but the interpretation of the operators is slightly different. A time-dependent generalization is then made, yielding an asymptotic expansion in decreasing powers of N\sqrt{N}N​ for the equations of motion. In this expansion the condensate wave function has the time-dependent form $\ensuremath{\xi}(\mathbf{x},t)$, and the condition for the validity of the expansion is that $\ensuremath{\xi}(\mathbf{x},t)$ satisfies the time-dependent Gross-Pitaevskii equation $\ensuremath{\partial}\ensuremath{\xi}/\ensuremath{\partial}t=\ensuremath{-}({\ensuremath{\Elzxh}}^{2}/2m){\ensuremath{\nabla}}^{2}\ensuremath{\xi}+V\ensuremath{\xi}+Nu|\ensuremath{\xi}{|}^{2}\ensuremath{\xi}$. The physics is then described in a kind of interaction picture, called the condensate picture, in which the phonon operator can be expressed as $\ensuremath{\chi}(\mathbf{x},t)={\ensuremath{\sum}}_{k}{\ensuremath{\xi}}_{k}(\mathbf{x},t){\ensuremath{\alpha}}_{k},$ where the operators ${\ensuremath{\alpha}}_{k}$ are time-independent annihilation operators, and the state vector has a time evolution described by a Schr\"odinger equation in which the Hamiltonian is a time-dependent quadratic form in the phonon creation and annihilation operators, whose coefficients are explicitly determined in terms of the time-dependent condensate wave function $\ensuremath{\xi}(\mathbf{x},t)$.
We consider the dynamics of |NNN| bosons in three dimensions. We assume that the pair interaction is given by |N3β−1V(Nβ⋅)N^{3\beta -1}V(N^{\beta }\cdot )N3β−1V(Nβ⋅)|⁠. By studying an associated many-body wave operator, we introduce a Bogoliubov–Born–Green–Kirkwood–Yvon (BBGKY) hierarchy which takes into account all of the interparticle singular correlation structures developed by the many-body evolution from the beginning. Assuming energy conditions on the |NNN|-body wave function, for |$\beta \in ( 0,1] ∣⁠,wederivetheGross–Pitaevskiihierarchywith2−bodyinteraction.Inparticular,weestablishthat,inthe∣|⁠, we derive the Gross–Pitaevskii hierarchy with 2-body interaction. In particular, we establish that, in the |∣⁠,wederivetheGross–Pitaevskiihierarchywith2−bodyinteraction.Inparticular,weestablishthat,inthe∣N\rightarrow \infty ∣limit,all∣| limit, all |∣limit,all∣k∣−bodyscatteringprocessesvanishif∣|-body scattering processes vanish if |∣−bodyscatteringprocessesvanishif∣k\geqslant 3∣andthusprovideadirectanswertoaquestionraisedbyErdo¨setal.[30].Moreover,thisnewBBGKYhierarchysharesthelimitpointswiththeordinaryBBGKYhierarchystronglyfor∣| and thus provide a direct answer to a question raised by Erdös et al. [30]. Moreover, this new BBGKY hierarchy shares the limit points with the ordinary BBGKY hierarchy strongly for |∣andthusprovideadirectanswertoaquestionraisedbyErdo¨setal.[30].Moreover,thisnewBBGKYhierarchysharesthelimitpointswiththeordinaryBBGKYhierarchystronglyfor∣\beta \in ( 0,1) ∣andweaklyfor∣| and weakly for |∣andweaklyfor∣\beta =1∣⁠.SincethisnewBBGKYhierarchyconvertstheproblemfromatwo−bodyestimatetoaweakerthree−bodyestimateforwhichwehavetheestimatestoachieve∣|⁠. Since this new BBGKY hierarchy converts the problem from a two-body estimate to a weaker three-body estimate for which we have the estimates to achieve |∣⁠.SincethisnewBBGKYhierarchyconvertstheproblemfromatwo−bodyestimatetoaweakerthree−bodyestimateforwhichwehavetheestimatestoachieve∣\beta <1∣⁠,itthenallowsustoprovethatalllimitpointsoftheordinaryBBGKYhierarchysatisfythespace–timeboundconjecturedbyKlainermanandMachedon[48]for∣|⁠, it then allows us to prove that all limit points of the ordinary BBGKY hierarchy satisfy the space–time bound conjectured by Klainerman and Machedon [48] for |∣⁠,itthenallowsustoprovethatalllimitpointsoftheordinaryBBGKYhierarchysatisfythespace–timeboundconjecturedbyKlainermanandMachedon[48]for∣\beta \in (0,1)$|⁠.
We clarify that an ideal gas obeying infinite statistics cannot undergo condensation. Then we derive the dynamic equation for an identical particle system obeying infinite statistics under external potential and inter-particle interaction. The derivation utilizes the Hamiltonian written in terms of the number operators and the transition number operators. At a very low temperature, where one can discard the dynamics of the excited occupation level, the dynamic of an infinite statistics system can be described by the Gross Pitaevskii equation, similar to the Bose-Einstein case.
Consider a system of NNN bosons in three dimensions interacting via a repulsive short range pair potential N2V(N(xi−xj))N^2V(N(x_i-x_j))N2V(N(xi​−xj​)), where $\bx=(x_1, >..., x_N)$ denotes the positions of the particles. Let HNH_NHN​ denote the Hamiltonian of the system and let ψN,t\psi_{N,t}ψN,t​ be the solution to the Schrodinger equation. Suppose that the initial data ψN,0\psi_{N,0}ψN,0​ satisfies the energy condition ≤CN>. \leq C N >. ≤CN>. and that the one-particle density matrix converges to a projection as N→∞N \to \inftyN→∞. Then, we prove that the kkk-particle density matrices of ψN,t\psi_{N,t}ψN,t​ factorize in the limit N→∞N \to \inftyN→∞. Moreover, the one particle orbital wave function solves the time-dependent Gross-Pitaevskii equation, a cubic non-linear Schrodinger equation with the coupling constant proportional to the scattering length of the potential VVV. In \cite{ESY}, we proved the same statement under the condition that the interaction potential VVV is sufficiently small; in the present work we extend this result to large potentials.
Consider a system of NNN bosons in three dimensions interacting via a repulsive short range pair potential N2V(N(xi−xj))N^2V(N(x_i-x_j))N2V(N(xi​−xj​)), where $\bx=(x_1, >..., x_N)$ denotes the positions of the particles. Let HNH_NHN​ denote the Hamiltonian of the system and let ψN,tψ_{N,t}ψN,t​ be the solution to the Schrödinger equation. Suppose that the initial data ψN,0ψ_{N,0}ψN,0​ satisfies the energy condition <ψN,0,HNkψN,0>≤CkNk < ψ_{N,0}, H_N^k ψ_{N,0} > \leq C^k N^k <ψN,0​,HNk​ψN,0​>≤CkNk for $k=1,2,... $. We also assume that the kkk-particle density matrices of the initial state are asymptotically factorized as N→∞N\to\inftyN→∞. We prove that the kkk-particle density matrices of ψN,tψ_{N,t}ψN,t​ are also asymptotically factorized and the one particle orbital wave function solves the Gross-Pitaevskii equation, a cubic non-linear Schrödinger equation with the coupling constant given by the scattering length of the potential VVV. We also prove the same conclusion if the energy condition holds only for k=1k=1k=1 but the factorization of ψN,0ψ_{N,0}ψN,0​ is assumed in a stronger sense.