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

arXiv:math-ph/0606017 · Annals of Mathematics

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

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

Published Jun 27, 2010• 328 citations
arXivHTMLPDF

Abstract

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.

DOI
10.4007/annals.2010.172.291

Related papers

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 hierarchy for the dynamics of Bose‐Einstein condensate

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

Abstract Consider a system of N bosons on the three‐dimensional unit torus interacting via a pair potential N2V(N(xi − xj)) where x = (x1, …, xN) denotes the positions of the particles. Suppose that the initial data ψN, 0 satisfies the condition where HN is the Hamiltonian of the Bose system. This condition is satisfied if ψN, 0 = WNϕN, 0 where WN is an approximate ground state to HN and ϕN, 0 is regular. Let ψN, t denote the solution to the Schrödinger equation with Hamiltonian HN. Gross and Pitaevskii proposed to model the dynamics of such a system by a nonlinear Schrödinger equation, the Gross‐Pitaevskii (GP) equation. The GP hierarchy is an infinite BBGKY hierarchy of equations so that if ut solves the GP equation, then the family of k‐particle density matrices ⊗k|ut〉 〈ut| solves the GP hierarchy. We prove that as N → ∞ the limit points of the k‐particle density matrices of ψN, t are solutions of the GP hierarchy. Our analysis requires that the N‐boson dynamics be described by a modified Hamiltonian that cuts off the pair interactions whenever at least three particles come into a region with diameter much smaller than the typical interparticle distance. Our proof can be extended to a modified Hamiltonian that only forbids at least n particles from coming close together for any fixed n. © 2006 Wiley Periodicals, Inc.

arXiv:math-ph/0410005•Communications on Pure and Applied Mathematics•Feb 2, 2006•171 citations

Derivation of the 2d Gross–Pitaevskii Equation for Strongly Confined 3d Bosons

Lea Boßmann

Abstract We study the dynamics of a system of N interacting bosons in a disc-shaped trap, which is realised by an external potential that confines the bosons in one spatial dimension to an interval of length of order $$\varepsilon $$ ε . The interaction is non-negative and scaled in such a way that its scattering length is of order $$\varepsilon /N$$ ε / N , while its range is proportional to $$(\varepsilon /N)^{\beta }$$ ( ε / N ) β with scaling parameter $$\beta \in (0,1]$$ β ∈ ( 0 , 1 ] . We consider the simultaneous limit $$(N,\varepsilon )\rightarrow (\infty ,0)$$ ( N , ε ) → ( ∞ , 0 ) and assume that the system initially exhibits Bose–Einstein condensation. We prove that condensation is preserved by the N-body dynamics, where the time-evolved condensate wave function is the solution of a two-dimensional non-linear equation. The strength of the non-linearity depends on the scaling parameter $$\beta $$ β . For $$\beta \in (0,1)$$ β ∈ ( 0 , 1 ) , we obtain a cubic defocusing non-linear Schrödinger equation, while the choice $$\beta =1$$ β = 1 yields a Gross–Pitaevskii equation featuring the scattering length of the interaction. In both cases, the coupling parameter depends on the confining potential.

arXiv:1907.04547•Archive for Rational Mechanics and Analysis•Jul 18, 2020•10 citations

Derivation of the Gross-Pitaevskii Equation from Quantum Dynamics of Many-Body Systems

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

We prove rigorously that the one-particle density matrix of interacting Bose systems with a short-scale repulsive pair interaction converges to the solution of the Gross-Pitaevskii equation in suitable scaling limits. The result is extended to k-particle density matrices for all positive integer k.

•Aug 2, 2005•1 citation

Gross-Pitaevskii non-linear dynamics for pseudo-spinor condensates

Alessandro Michelangeli, Alessandro Olgiati

We derive the equations for the non-linear effective dynamics of a so called pseudo-spinor Bose-Einstein condensate, which emerges from the linear many-body Schrödinger equation at the leading order in the number of particles. The considered system is a three-dimensional diluted gas of identical bosons with spin, possibly confined in space, and coupled with an external time-dependent magnetic field; particles also interact among themselves through a short-scale repulsive interaction. The limit of infinitely many particles is monitored in the physically relevant Gross-Pitaevskii scaling. In our main theorem, if at time zero the system is in a phase of complete condensation (at the level of the reduced one-body marginal) and with energy per particle fixed by the Gross-Pitaevskii functional, then such conditions persist also at later times, with the one-body orbital of the condensate evolving according to a system of non-linear cubic Schrödinger equations coupled among themselves through linear (Rabi) terms. The proof relies on an adaptation to the spinor setting of Pickl's projection counting method developed for the scalar case. Quantitative rates of convergence are available, but not made explicit because evidently non-optimal. In order to substantiate the formalism and the assumptions made in the main theorem, in an introductory section we review the mathematical formalisation of modern typical experiments with pseudo-spinor condensates.

arXiv:1704.00150•Apr 1, 2017•5 citations

Gross-Pitaevskii non-linear dynamics for pseudo-spinor condensates

Alessandro Michelangeli, Alessandro Olgiati

We derive the equations for the non-linear effective dynamics of a so called pseudo-spinor Bose-Einstein condensate, which emerges from the linear many-body Schrödinger equation at the leading order in the number of particles.The considered system is a three-dimensional diluted gas of identical bosons with spin, possibly confined in space, and coupled with an external time-dependent magnetic field; particles also interact among themselves through a short-scale repulsive interaction.The limit of infinitely many particles is monitored in the physically relevant Gross-Pitaevskii scaling.In our main theorem, if at time zero the system is in a phase of complete condensation (at the level of the reduced one-body marginal) and with energy per particle fixed by the Gross-Pitaevskii functional, then such conditions persist also at later times, with the one-body orbital of the condensate evolving according to a system of non-linear cubic Schrödinger equations coupled among themselves through linear (Rabi) terms.The proof relies on an adaptation to the spinor setting of Pickl's projection counting method developed for the scalar case.Quantitative rates of convergence are available, but not made explicit because evidently non-optimal.In order to substantiate the formalism and the assumptions made in the main theorem, in an introductory section we review the mathematical formalisation of modern typical experiments with pseudo-spinor condensates.

Journal of Nonlinear Mathematical Physics•Jun 28, 2017•18 citations
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=(x1,>...,xN)\bx=(x_1, >..., x_N)\bx=(x1​,>...,xN​) 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}, H_N^k ψ_{N,0} > \leq C^k N^k ] for k=1,2,...k=1,2,... 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.
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.
Abstract Consider a system of N bosons on the three‐dimensional unit torus interacting via a pair potential N2V(N(xi − xj)) where x = (x1, …, xN) denotes the positions of the particles. Suppose that the initial data ψN, 0 satisfies the condition where HN is the Hamiltonian of the Bose system. This condition is satisfied if ψN, 0 = WNϕN, 0 where WN is an approximate ground state to HN and ϕN, 0 is regular. Let ψN, t denote the solution to the Schrödinger equation with Hamiltonian HN. Gross and Pitaevskii proposed to model the dynamics of such a system by a nonlinear Schrödinger equation, the Gross‐Pitaevskii (GP) equation. The GP hierarchy is an infinite BBGKY hierarchy of equations so that if ut solves the GP equation, then the family of k‐particle density matrices ⊗k|ut〉 〈ut| solves the GP hierarchy. We prove that as N → ∞ the limit points of the k‐particle density matrices of ψN, t are solutions of the GP hierarchy. Our analysis requires that the N‐boson dynamics be described by a modified Hamiltonian that cuts off the pair interactions whenever at least three particles come into a region with diameter much smaller than the typical interparticle distance. Our proof can be extended to a modified Hamiltonian that only forbids at least n particles from coming close together for any fixed n. © 2006 Wiley Periodicals, Inc.
Abstract We study the dynamics of a system of N interacting bosons in a disc-shaped trap, which is realised by an external potential that confines the bosons in one spatial dimension to an interval of length of order ε\varepsilon ε ε . The interaction is non-negative and scaled in such a way that its scattering length is of order ε/N\varepsilon /Nε/N ε / N , while its range is proportional to (ε/N)β(\varepsilon /N)^{\beta }(ε/N)β ( ε / N ) β with scaling parameter β∈(0,1]\beta \in (0,1]β∈(0,1] β ∈ ( 0 , 1 ] . We consider the simultaneous limit (N,ε)→(∞,0)(N,\varepsilon )\rightarrow (\infty ,0)(N,ε)→(∞,0) ( N , ε ) → ( ∞ , 0 ) and assume that the system initially exhibits Bose–Einstein condensation. We prove that condensation is preserved by the N-body dynamics, where the time-evolved condensate wave function is the solution of a two-dimensional non-linear equation. The strength of the non-linearity depends on the scaling parameter β\beta β β . For β∈(0,1)\beta \in (0,1)β∈(0,1) β ∈ ( 0 , 1 ) , we obtain a cubic defocusing non-linear Schrödinger equation, while the choice β=1\beta =1β=1 β = 1 yields a Gross–Pitaevskii equation featuring the scattering length of the interaction. In both cases, the coupling parameter depends on the confining potential.
We derive the equations for the non-linear effective dynamics of a so called pseudo-spinor Bose-Einstein condensate, which emerges from the linear many-body Schrödinger equation at the leading order in the number of particles. The considered system is a three-dimensional diluted gas of identical bosons with spin, possibly confined in space, and coupled with an external time-dependent magnetic field; particles also interact among themselves through a short-scale repulsive interaction. The limit of infinitely many particles is monitored in the physically relevant Gross-Pitaevskii scaling. In our main theorem, if at time zero the system is in a phase of complete condensation (at the level of the reduced one-body marginal) and with energy per particle fixed by the Gross-Pitaevskii functional, then such conditions persist also at later times, with the one-body orbital of the condensate evolving according to a system of non-linear cubic Schrödinger equations coupled among themselves through linear (Rabi) terms. The proof relies on an adaptation to the spinor setting of Pickl's projection counting method developed for the scalar case. Quantitative rates of convergence are available, but not made explicit because evidently non-optimal. In order to substantiate the formalism and the assumptions made in the main theorem, in an introductory section we review the mathematical formalisation of modern typical experiments with pseudo-spinor condensates.
We derive the equations for the non-linear effective dynamics of a so called pseudo-spinor Bose-Einstein condensate, which emerges from the linear many-body Schrödinger equation at the leading order in the number of particles.The considered system is a three-dimensional diluted gas of identical bosons with spin, possibly confined in space, and coupled with an external time-dependent magnetic field; particles also interact among themselves through a short-scale repulsive interaction.The limit of infinitely many particles is monitored in the physically relevant Gross-Pitaevskii scaling.In our main theorem, if at time zero the system is in a phase of complete condensation (at the level of the reduced one-body marginal) and with energy per particle fixed by the Gross-Pitaevskii functional, then such conditions persist also at later times, with the one-body orbital of the condensate evolving according to a system of non-linear cubic Schrödinger equations coupled among themselves through linear (Rabi) terms.The proof relies on an adaptation to the spinor setting of Pickl's projection counting method developed for the scalar case.Quantitative rates of convergence are available, but not made explicit because evidently non-optimal.In order to substantiate the formalism and the assumptions made in the main theorem, in an introductory section we review the mathematical formalisation of modern typical experiments with pseudo-spinor condensates.