---
title: Extended Gross–Pitaevskii Framework
url: https://www.emergentmind.com/topics/extended-gross-pitaevskii-framework
type: topic
---

# Extended Gross–Pitaevskii Framework

The expression *extended Gross–Pitaevskii framework* denotes a family of non-equivalent generalizations of the Gross–Pitaevskii equation constructed to relax one or more of the standard GP assumptions: a dilute Bose gas at \(T\approx 0\), dominance of \(s\)-wave scattering, zero-range contact interactions modeled by a \(\delta\)-pseudopotential, and a mean-field Hartree description in which all atoms occupy the same orbital \(\psi\). Across the literature, the same label is used for derivative corrections that encode finite interaction range, Lee–Huang–Yang-augmented equations for droplets, rigorous condensate-plus-Bogoliubov descriptions in the Gross–Pitaevskii regime, matrix-valued quantum equations on the continuous matrix product state manifold, nonlocal self-interaction-subtracted theories for excitonic solids, and Gross–Pitaevskii–Poisson systems for self-gravitating superfluids [1404.2586], [1501.06575], [2308.11687], [2507.20236], [2512.16936], [2604.04460].

## 1. Baseline Gross–Pitaevskii structure and the scope of extension

In its canonical form, the time-dependent Gross–Pitaevskii equation is
\[
i\hbar\,\partial_t\psi(\mathbf{r},t)
=
\left[
-\frac{\hbar^2}{2m}\nabla^2 + V(\mathbf{r}) + g\,|\psi(\mathbf{r},t)|^2
\right]\psi(\mathbf{r},t),
\qquad
g=\frac{4\pi\hbar^2 a_s}{m},
\]
with energy functional
\[
E[\psi]
=
\int d^3r\,
\left[
\frac{\hbar^2}{2m}|\nabla\psi|^2
+
V(\mathbf{r})|\psi|^2
+
\frac{g}{2}|\psi|^4
\right].
\]
The order parameter may be normalized either to the particle number \(N\) or to \(1\), with the factor \(N\) absorbed into the nonlinear coefficient. This baseline is the reference point from which all “extended” GP constructions depart [1404.2586].

The extensions in the literature are not a single hierarchy. Some retain a scalar order parameter but modify the interaction term; some keep the mean-field PDE but change its equation of state; some replace the scalar wavefunction by matrix-valued fields; and some treat the GP equation only as the condensate sector of a larger condensate-plus-fluctuation dynamics. A compact classification is therefore useful.

| Extension class | Defining modification | Representative sources |
|---|---|---|
| Finite-range or effective-range GP | Adds \(g_2\nabla^2(|\psi|^2)\psi\) or equivalent short-range nonlocal correction | [1404.2586] |
| LHY-augmented eGPE | Adds a local beyond-mean-field term proportional to \(n^{3/2}\) or \(|\psi|^3\psi\) | [1805.11325], [2604.04460] |
| Rigorous GP-plus-Bogoliubov dynamics | Keeps GP evolution for the condensate and adds microscopic correlation dressing and quadratic fluctuation dynamics | [2308.11687], [2310.11347] |
| Multi-component hard-core GP derivation | Derives coupled GP equations with uniform-in-time convergence bounds | [2501.18787] |
| Pairing-field droplet theory | Interprets the eGPE as the equation of motion for a pairing field rather than the condensate | [2006.05769] |
| Quantum/cMPS GP | Replaces the scalar field by matrix-valued cMPS variables and derives a non-commutative QGPE | [1501.06575] |
| Nonlocal correlated and self-interaction-subtracted GP | Uses a nonlocal kernel and explicitly excludes intra-cell self-interactions | [2507.20236] |
| Self-gravitating GPPE | Couples GP dynamics to Poisson gravity, and in neutron-star applications to Maxwell fields and multi-component order parameters | [2512.16936] |

This plurality is the central feature of the subject. A statement about an “extended GP equation” is therefore only meaningful together with the mechanism of extension.

## 2. Finite-range, nonlocal, and long-range interaction extensions

One of the most direct generalizations of the GP equation replaces the zero-range contact approximation by an effective-range correction. A simple model uses a repulsive core plus an attractive thin shell, which in the zero-thickness limit generates a sum of \(\delta\)-functions and captures the essential structure of van der Waals interactions while remaining analytically tractable. In the resulting stationary equation, the leading finite-range correction is a derivative term independent of the microscopic details of the interaction potential,
\[
\mu\,\psi(\mathbf{r})
=
-\frac{\hbar^2}{2m}\nabla^2\psi
+
U(\mathbf{r})\psi
+
N g_{3D}|\psi|^2\psi
+
N g_{3D} g_2 \nabla^2\!\left(|\psi|^2\right)\psi,
\]
with
\[
g_2 = a\left(\frac{a}{3}-\frac{r_e}{2}\right).
\]
In an elongated trap, the same physics can be recast in a simple 1D three-\(\delta\) model,
\[
V(x)=2g\,\delta(x)-\frac{g}{2}\left[\delta(x+l)+\delta(x-l)\right],
\qquad
\frac{1}{2}l^2=-g_2,
\]
so that effective-range information is encoded in a single length scale \(l\) [1404.2586].

In the Thomas–Fermi regime of a 1D harmonic trap, this extension yields closed-form leading-order corrections. The density acquires a constant offset, the chemical potential shifts by
\[
\Delta\mu=\frac{1}{4}m\omega^2 l^2,
\]
and the per-particle energy gains
\[
\Delta E(l)=\frac{1}{4}m l^2\omega^2.
\]
The relative correction is
\[
\frac{\Delta E(l)}{E(l=0)}=\frac{5}{6}\frac{l^2}{R^2}=-\frac{5}{3}\frac{g_2}{R^2}.
\]
A distinctive feature is that \(\Delta E\propto \omega^2\), whereas the standard GP Thomas–Fermi energy scales as \(\omega^{2/3}\) in the 1D setting analyzed there. This scaling separation is the proposed experimental handle for finite-range physics [1404.2586].

The same derivative viewpoint fits into a broader nonlocal perspective. For short-range nonlocal kernels, the leading effect can be written as \(g_2\nabla^2(|\psi|^2)\psi\), whereas for genuinely long-range kernels one keeps an explicit convolution. A different long-range extension appears in self-gravitating superfluids, where the Gross–Pitaevskii–Poisson equation takes the form
\[
i\hbar\,\partial_t \psi
=
-\frac{\hbar^2}{2m}\nabla^2\psi
+
g|\psi|^2\psi
+
m\Phi\,\psi,
\qquad
\nabla^2\Phi
=
4\pi G\,(m|\psi|^2-\rho_{\rm bg}).
\]
This framework is used for boson and axion stars, dark-matter haloes, and minimal neutron-star models; for axions it is extended further to a cubic–quintic GPPE, and in neutron-star applications it is coupled to a proton superconductor, the Maxwell equations for the vector potential \(\mathbf{A}\), and Poisson gravity [2512.16936].

These interaction-based extensions remain closest to the original GP philosophy: a macroscopic order parameter is preserved, but the contact Hartree term is replaced by a more structured interaction functional.

## 3. Quantum-fluctuation corrections and droplet eGPEs

A second major branch of the extended GP literature incorporates beyond-mean-field quantum fluctuations through a local Lee–Huang–Yang correction. In uniform matter, linearized time-dependent GP theory reproduces the Bogoliubov dispersion,
\[
E_k=\sqrt{\varepsilon_k(\varepsilon_k+2gn)},
\qquad
\varepsilon_k=\frac{\hbar^2k^2}{2m},
\]
the static structure factor \(S(q)=\varepsilon_q/E_q\), the momentum distribution
\[
n_k=\frac{\varepsilon_k+gn}{2E_k}-\frac{1}{2},
\]
and the quantum depletion fraction
\[
\frac{N_{\rm ex}}{N}
=
\frac{8}{3\sqrt{\pi}}\sqrt{\rho a^3}+o\!\big(\sqrt{\rho a^3}\big).
\]
Within the same framework, the chemical potential acquires the LHY form
\[
\mu(n)=gn\left[1+\frac{32}{3\sqrt{\pi}}(na^3)^{1/2}\right],
\]
which corresponds to an additional local term proportional to \(|\psi|^3\psi\) in the extended GP equation [1805.11325].

In dimensionless notation, a widely used \(d\)-dimensional eGPE with LHY correction is
\[
i\partial_t\psi(\mathbf{x},t)
=
-\frac{1}{2}\nabla^2\psi
+
V(\mathbf{x})\psi
+
\beta |\psi|^2\psi
+
\lambda |\psi|^3\psi,
\]
with energy
\[
E(\psi)
=
\int_{\mathbb R^d}
\left(
\frac{1}{2}|\nabla\psi|^2
+
V(\mathbf{x})|\psi|^2
+
\frac{\beta}{2}|\psi|^4
+
\frac{2\lambda}{5}|\psi|^5
\right)d\mathbf{x}.
\]
For \(\lambda>0\) and free space, the ground-state problem exhibits dimension-dependent existence properties: if \(\beta\ge 0\), \(\gamma(c)=0\) and there is no minimizer; if \(\beta<0\), minimizers exist for all \(c>0\) in \(d=1\), while in \(d=2,3\) there is a threshold \(c_d^*\) below which no minimizer exists and above which minimizers exist. Under a confining potential \(V\ge 0\) with \(V(\mathbf{x})\to\infty\) as \(|\mathbf{x}|\to\infty\), a minimizer exists for every \(c>0\) and \(\beta\in\mathbb R\) [2604.04460].

This equation supports distinct free-space regimes that the paper organizes as no-ground-state, soliton-like, and droplet-like. In the droplet regime, a simple flat-top approximation neglects kinetic energy and yields
\[
a=-\frac{5\beta}{6\lambda},
\qquad
E_{\mathrm{app}}^{\min}
=
\frac{\beta c}{6}\left(\frac{5\beta}{6\lambda}\right)^2,
\]
which captures the high-density core of self-bound droplets [2604.04460].

Dipolar quantum droplets provide a stringent test of this eGPE logic. The eGPE used for dysprosium droplets adds dipolar mean-field interactions and a local LHY chemical potential
\[
\mu_{\mathrm{LHY}}(n)=\gamma_{\mathrm{QF}}\,n^{3/2},
\qquad
\gamma_{\mathrm{QF}}
=
\frac{32\,g\,a_s^{3/2}}{3\sqrt{\pi}}
\left(1+\frac{3}{2}\varepsilon_{\mathrm{dd}}^2\right),
\]
together with a three-body loss term in dynamical simulations. Although the central density of a saturated Dy droplet corresponds to \(n a_s^3\approx 2\times 10^{-4}\), the measured critical atom numbers are systematically lower than the eGPE prediction, while quantum Monte Carlo calculations agree well with experiment. This directly questions the validity of the current theoretical state-of-the-art description of dilute dipolar droplets within the eGPE framework and indicates that correlations are significant [1904.10349].

A plausible implication is that “extended GP” in the LHY sense is quantitatively reliable only within a limited correlation window. The formal addition of a local \(n^{5/2}\) energy density is not automatically sufficient when nonlocality, finite-range effects, or strong dipolar correlations become important.

## 4. Microscopic derivations, rigorous many-body theory, and coupled GP limits

A different meaning of extended GP emerges when the GP equation is derived as only one sector of a more complete many-body description. In the three-dimensional Gross–Pitaevskii regime, one rigorous framework lets the condensate evolve according to a time-dependent GP equation, factors out the microscopic two-body correlation structure through the zero-energy scattering solution, and describes the orthogonal excitations by a quadratic Bogoliubov dynamics. In this setting the many-body dynamics is approximated in norm with rate \(N^{-1/8}\), the Bogoliubov dynamics converges to an \(N\)-independent limit with rate \(N^{-1/4}\), and a central limit theorem for bounded observables is obtained with rate \(N^{-1/16}\) [2308.11687].

A complementary rigorous development diagonalizes dilute Bose gases algebraically in the GP regime and beyond. In the GP case it recovers the leading-order energy \(4\pi a\,N(N-1)\), yields the discrete Lee–Huang–Yang correction, identifies the excitation dispersion
\[
E(p)=\sqrt{p^4+16\pi a\,\rho\,p^2},
\]
and provides an upper bound with a logarithmic remainder of the conjectured next-order scale. This same analysis motivates an extended GP functional with an LHY term and an associated Euler–Lagrange equation valid in the dilute regime \(\rho a^3\ll 1\) [2310.11347].

For mixtures, the extension can be made fully quantitative. In the hard-core GP limit of a two-component Bose gas, the effective equations are
\[
\begin{aligned}
i\partial_t u &= -\Delta u + 8\pi a_{\mathrm{sc},1}\,\mathfrak{n}_1|u|^2u + 8\pi a_{\mathrm{sc},12}\,\mathfrak{n}_2|v|^2u,\\
i\partial_t v &= -\Delta v + 8\pi a_{\mathrm{sc},2}\,\mathfrak{n}_2|v|^2v + 8\pi a_{\mathrm{sc},12}\,\mathfrak{n}_1|u|^2v,
\end{aligned}
\]
with a modified \(N\)-dependent system incorporating Neumann scattering correctors used in the derivation. The main theorem proves a uniform-in-time trace-norm bound
\[
\operatorname{Tr}\left|
\Gamma_{N,t}^{(1)}-|\boldsymbol{\phi}_t\rangle\!\langle \boldsymbol{\phi}_t|
\right|
\le
\frac{C}{\sqrt N},
\]
and the same strategy extends to the one-component case and to any finite number of components \(p\) [2501.18787].

Binary quantum droplets introduce a further microscopic refinement. In a homonuclear Bose–Bose mixture with attractive inter-species interaction, the condensate wavefunction obeys
\[
\left[
-\frac{\hbar^2 \nabla^2}{2m}
+
V_T(x)
-
\mu
+
g |\phi_c(x)|^2
\right]\phi_c(x)
-
\tilde A(x)\phi_c^*(x)=0,
\]
where \(\tilde A(x)\) is a bosonic pairing field generated by a Hubbard–Stratonovich decoupling of the inter-species channel. The paper’s central claim is that the “extended Gross–Pitaevskii equation” used earlier for droplets should be understood as the equation of motion for the pairing field, rather than the condensate. In uniform matter this coupled theory yields a gapless branch
\[
E_-(k)=\sqrt{\varepsilon_k(\varepsilon_k+2C+2A)}
\]
and a gapped branch
\[
E_+(k)=\sqrt{(\varepsilon_k+2C)(\varepsilon_k+2A)},
\]
so that low-energy collective excitations are fluctuations of the pairing field, while the atomic Bogoliubov equations describe single-particle-like excitations [2006.05769].

This many-body branch of the subject therefore extends GP in a stricter sense than a modified nonlinear coefficient. It retains the condensate PDE but embeds it within a fluctuation theory, a scattering renormalization, or a coupled condensate–pairing system.

## 5. Strong-correlation and non-mean-field reformulations

In one dimension, the failure of the ordinary GPE is not merely quantitative. A non-commutative *quantum Gross–Pitaevskii equation* is obtained by applying the time-dependent variational principle to the variational manifold of continuous matrix product states. The dynamical variables are matrix-valued fields \(R(x,t)\), \(Q(x,t)\), and \(P(x,t)\), together with left and right density matrices \(\rho_L\) and \(\rho_R\). The gauge-covariant QGPE reads
\[
\mathscr{D}_t R
=
\big(-\mathscr{D}_x^2+v(x)\big)R
+
g(\rho_L^{-1}R^\dagger \rho_L)R^2
+
gR^2(\rho_R R^\dagger \rho_R^{-1})
-
(\rho_L^{-1}R^\dagger \rho_L)[R,\mathscr{D}_xR]
-
[R,\mathscr{D}_xR](\rho_R R^\dagger \rho_R^{-1}),
\]
supplemented by an equation for \(Q\) and by gauge-density equations. When the bond dimension is \(D=1\), all commutators vanish and the QGPE reduces to the ordinary GPE; for \(D>1\), it retains entanglement and correlations. Linearization yields quantum Bogoliubov–de Gennes equations whose steady-state response shows a pronounced feature near \(k=2k_F\), consistent with the Tonks–Girardeau limit and type II excitations [1501.06575].

A still more radical extension appears in the theory of an excitonic incompressible Bose solid in double-layer semiconductor heterostructures. There the effective interaction kernel is nonlocal and correlated,
\[
t_{\mathrm{XX}}(r)=
\begin{cases}
V_{\mathrm{XX}}(r), & r>R_{\mathrm{HC}},\\
0, & r\le R_{\mathrm{HC}},
\end{cases}
\]
with the hard-core radius \(R_{\mathrm{HC}}\) extracted from QMC structure factors of 2D dipoles. The generalized stationary GP equation is
\[
-\frac{\hbar^2}{2M_{\mathrm X}}\nabla^2\Psi(\mathbf r)
+
U_{\mathrm{XX}}(\mathbf r)\Psi(\mathbf r)
=
\mu\,\Psi(\mathbf r),
\]
where
\[
U_{\mathrm{XX}}(\mathbf r)
=
\int d^2\mathbf r'\, t_{\mathrm{XX}}(|\mathbf r'-\mathbf r|)\,|\Psi(\mathbf r')|^2
-
\int_{\mathbf r'\in i_{\mathbf r}} d^2\mathbf r'\, t_{\mathrm{XX}}(|\mathbf r'-\mathbf r|)\,|\Psi(\mathbf r')|^2.
\]
The second term removes the particle’s own-cell contribution and explicitly excludes self-interactions. This prescription stabilizes an incompressible supersolid with exactly one boson per lattice site and a finite homogeneous background \(\rho_B\). Numerically, at fixed low density \(r_0=30\), the solutions are uniform superfluid-like at \(d=1\), become supersolid at \(d=7\), and the energy eigenvalues cross around \(d\simeq 6\); the stationary problem is solved by imaginary-time propagation on a \(256^2\) grid with convergence at relative energy change \(<10^{-8}\) [2507.20236].

This framework also clarifies a common misconception. Cluster supersolids in dipolar gases are generally compressible and arise from a competition between attractive and repulsive forces, whereas the exciton supersolid uses purely repulsive interactions together with explicit self-interaction exclusion to enforce single occupancy and incompressibility [2507.20236].

## 6. Solvable constructions, perturbative schemes, computation, and conceptual ambiguities

Some extended GP frameworks are analytic rather than physical extensions. One construction defines a family of exactly solvable stationary GP-type equations with variable coefficients and arbitrary nonlinear powers by the transform
\[
x=|\xi|^\sigma,
\qquad
\psi(x)=|\xi|^{\frac{\sigma-1}{2}}\phi(\xi),
\]
together with explicit coefficient maps
\[
V_{\mathrm{eff}}(\xi)
=
\sigma^2\left[
\frac{1-\sigma^{-2}}{4\xi^2}
+
|\xi|^{2(\sigma-1)}U_{\mathrm{eff}}(|\xi|^\sigma)
\right],
\qquad
G_l(\xi)
=
\sigma^2 |\xi|^{(\sigma-1)\frac{l+4}{2}} g_l(|\xi|^\sigma).
\]
The resulting “family” transfers exact solvability from one member to all others and covers the nonlinear Schrödinger equation, quintic GP, cubic–quintic GP, generalized polynomial nonlinearities, and a 3D spherically symmetric radial GP-type equation [2004.09972].

Another analytic route is perturbative. For the time-independent GP equation in \(\mathbb R^D\) with general external potential, a Green’s-function/eigenbasis expansion yields first-order corrections to the wavefunction, density, chemical potential, internal energy, and mean-square radius. The true small parameter is \(\tilde g_D N\), and the validity criterion is
\[
(\tilde g_D N)^2 S^{(1)}\ll 1.
\]
Applied to harmonic traps in \(D=1,2,3\), the method gives explicit first-order formulas and predicts breakdown once the interaction-driven admixture of excited states is no longer small; the paper notes, for example, breakdown around \(N\gtrsim 10^2\) in the 2D and 3D examples discussed there [2303.02755].

The numerical literature extends the framework in a different direction. For steady symmetric and central vortex states of nonlocal GP equations coupled to microwaves, two accelerated solvers are proposed: the accelerated stabilized normalized gradient flow (ASGF) method and the perturbed, projected nonlinear conjugate gradient (PPNCG) method. The governing equations are discretized by a Legendre–Galerkin spectral method in polar coordinates, and the PPNCG method performs the best in the reported benchmarks, consistently reducing iteration counts and CPU times relative to ASGF and earlier normalized-gradient methods [2201.12343].

Two recurring conceptual cautions follow from this body of work. First, the label *extended Gross–Pitaevskii equation* is not unique: it may denote a gradient-corrected mean-field PDE, an LHY-augmented droplet equation, a pairing-field equation, a condensate-plus-Bogoliubov dynamical scheme, or a matrix-valued quantum evolution. Second, several papers explicitly challenge over-identification of the extension with a single scalar condensate equation. In binary Bose-mixture droplets, the earlier eGPE is reinterpreted as the equation of motion for the pairing field [2006.05769]. In dilute dipolar droplets, disagreement between experiment, eGPE, and quantum Monte Carlo shows that the current theoretical state-of-the-art description within the eGPE framework can miss significant correlations [1904.10349].

Taken together, these developments suggest that the extended Gross–Pitaevskii framework is best understood not as one equation but as a research program: a set of controlled modifications of the GP paradigm, each designed to encode a specific physical ingredient—finite range, quantum fluctuations, hard-core correlations, pairing, multicomponent coupling, nonlocal kernels, self-gravity, or exact solvability—while preserving as much as possible of the variational and dynamical structure that made the original GP equation useful.

Source: https://www.emergentmind.com/topics/extended-gross-pitaevskii-framework