Gross–Pitaevskii Hydrodynamics

• Gross–Pitaevskii hydrodynamics is the study of the fluid-like behavior of dilute Bose–Einstein condensates, interpreted via the Madelung transformation of the GP equation into conservation laws.
• It reveals nonlinear phenomena such as quantized vortices and KPZ scaling in low-dimensional systems, offering insights into quantum turbulence.
• Extensions beyond mean field and multi-component systems connect GP hydrodynamics with classical fluid models, enhancing both numerical and experimental investigations.

Gross–Pitaevskii hydrodynamics is the study of the collective fluid-like behavior of dilute Bose–Einstein condensates (BECs) and other nonlinear quantum systems, as described by the Gross–Pitaevskii equation (GPE) or nonlinear Schrödinger equation (NLS). By transforming the GPE via density and phase variables (the Madelung representation), one obtains macroscopic conservation laws and effective equations analogous to classical hydrodynamics, while encoding quantum effects such as quantized vortices and dispersive (quantum pressure) regularization. Modern research situates Gross–Pitaevskii hydrodynamics at the crossroads of nonlinear PDEs, statistical physics, and the theory of universal scaling limits in low-dimensional systems.

1. Hydrodynamic Formulation via the Madelung Transform

The Gross–Pitaevskii equation for a complex order parameter $\psi(t,{\bf x})$ (typically for $T\approx0$ BECs) is

$$i\hbar \partial_t \psi = -\frac{\hbar^2}{2m} \Delta\psi + V({\bf x})\psi + g|\psi|^2 \psi,$$

where $m$ is particle mass, $V({\bf x})$ an external potential, and $g$ the contact interaction strength. Introducing the Madelung transform $\psi = \sqrt{\rho}\,e^{iS/\hbar}$ yields density $\rho$ and velocity ${\bf v} = \nabla S/m$ fields. The equations for $\rho$ and ${\bf v}$ separate into \begin{align*} \text{Continuity:} & \quad \partial_t \rho + \nabla\cdot(\rho {\bf v}) = 0,\ \text{Euler–like:} & \quad m\big[\partial_t {\bf v} + ({\bf v}\cdot\nabla){\bf v}\big] = -\nabla V - g\nabla\rho + \nabla Q, \end{align*} where $$Q(\rho) = -(\hbar^2/2m)\Delta\sqrt{\rho}/\sqrt{\rho}$$ is the “quantum pressure.” This identifies GP hydrodynamics as a compressible Euler-Korteweg system (with quantum capillarity term) (Carles et al., 2011, Andreev, 2012).

2. Nonlinear and Fluctuating Hydrodynamics in Low Dimensions

In one spatial dimension, linearizing Gross–Pitaevskii hydrodynamics yields chiral sound modes (Luttinger liquid limit) with undamped phonons. At finite temperature, nonlinear effects such as mode–mode interactions and stochasticity cause superdiffusive broadening of dynamical correlation functions that are invisible to linear theory. Systematic addition of diffusive and stochastic noise contributions, constrained by conservation laws and fluctuation–dissipation, leads to a nonlinear fluctuating hydrodynamics. Linearizing and diagonalizing yields two normal modes which, due to retained cubic nonlinearities, are governed by the 1D stochastic Burgers—and thus the KPZ—equation. This mapping quantitatively predicts that the dynamical structure factor $S(k,\omega)$ exhibits KPZ scaling, i.e., sound-peak broadening with dynamical exponent $z=3/2$, and universal scaling lineshapes computed numerically (Kulkarni et al., 2015).

Summary of the mode-coupling KPZ mapping:

System	Hydrodynamic limit	Nonlinear universality class
1D GP, low-$T$	Nonlinear fluctuating HD	KPZ, $z=3/2$
Integrable NLS	Ballistic with $\delta$-peaks	Not KPZ

The mapping is robust and applies to multicomponent GPEs (e.g., spinor BECs or birefringent fibers), where richer universality classes of coupled KPZ-type phenomena are expected (Kulkarni et al., 2015).

3. Vortex Dynamics and Hydrodynamic Limits

In dimensions $d\geq2$, the hydrodynamic limit of the GPE reveals quantized vortices as topological singularities of the phase, around which the circulation of ${\bf v}$ is quantized. Vortex dynamics are governed by singular hydrodynamic equations for the vorticity field, which, in the limit of vanishing healing length (core size $\varepsilon\to0$), can be rigorously connected to classical point-vortex ODEs and incompressible Euler flow:

• In 2D, with many vortices, the system under the suitable scaling limit converges to the classical incompressible Euler equations, as the discrete vorticity “melts” into a continuous distribution (Jerrard et al., 2013). The instantaneous position of the vortex cores follows the Kirchhoff–Onsager law.
• In 3D, the dynamics of (nearly) parallel vortex filaments emerging from the GPE can be rigorously reduced (under controlled rescaling) to Hamiltonian equations for curves in space, combining linearized self-induction (binormal flow) and logarithmic mutual interaction (Jerrard et al., 2020). The limiting dynamics recover and extend classical filament models, with rigorous error bounds.

Table: GP Vortex Limit Results

Dimension	Limit	Vortex Dynamics
2D	$\varepsilon\to0$, $n\gg 1$	Kirchhoff–Onsager ODE, Euler PDE
3D	$\varepsilon\to0$, parallel	Coupled Hamiltonian filaments

Unlike classical fluid models, the GPE framework provides a physically intrinsic core size and allows rigorous asymptotics.

4. Quantum Hydrodynamics Beyond Mean Field

The canonical GP hydrodynamics emerges in the first order of the many-body interaction radius expansion. Systematic corrections (third-order in interaction radius) introduce additional terms in the Euler equation:

$$m n(\partial_t + {\bf v}\cdot\nabla) v^\alpha + \partial^\alpha p_q + \partial^\alpha (g n^2) + \partial^\alpha [g_2 \nabla^2(n^2)] + \frac{1}{30}\partial_\beta [g_3 I^{\alpha\beta\gamma\delta} \partial_\gamma\partial_\delta (n^2)] = 0$$

Isotropic and anisotropic corrections (in $g_2,g_3$) arise from nonlocality and potential anisotropy, respectively (Andreev, 2020, Andreev, 2012).

Going beyond mean-field, a four-equation hydrodynamic hierarchy accounts for quantum fluctuations: density, velocity, pressure tensor, and third-rank pressure flux. Quantum fluctuation sources appear explicitly in pressure-tensor evolution, microscopically connected to Lee–Huang–Yang corrections (Andreev, 2020, Stringari, 2018).

5. Mathematical Well-Posedness and Functional Frameworks

Global-in-time well-posedness for the coupled hydrodynamic system (density, velocity) directly derived from GPE has been rigorously analyzed under no-vacuum assumptions (Wegner, 2023). By establishing local bilipschitz equivalence of suitable Sobolev-type function spaces under the Madelung transform, one transfers well-posedness of the GPE to the hydrodynamic variables. In 1D, for energy strictly below the critical threshold for vacuum formation, the Cauchy problem in $(1+H^s)\times H^{s-1}$ is globally well posed, and the mappings between wavefunction and hydrodynamic variables are quantitatively controlled.

6. Gross–Pitaevskii Hydrodynamics in Multi-Component and Coupled Systems

Gross–Pitaevskii hydrodynamics naturally extends to multi-component BECs, spinor condensates, and systems with dipolar or nonlocal interactions, leading to multi-field hydrodynamic formulations. Coupled systems, such as the two-fluid model for superfluid helium (Navier–Stokes + GP), require careful regularization and consistent mutual-friction terms, derived via Madelung hydrodynamics and coarse-graining. Fully consistent coupled GP–Navier–Stokes simulation frameworks, formulated with regularized vorticity and covariant gradients, enable the study of quantum turbulence, vortex reconnections, and friction across the quantum–classical interface (Brachet et al., 2022).

7. Experimental and Numerical Realizations

Experimentally, the hydrodynamics of the GPE is accessible in a wide range of settings, from ultracold atom BECs to controlled optical and fluid systems. Mechanical analogs (e.g., surface water waves with tunable currents) enable direct emulation of GP dynamics, including soliton formation, trapping, and modulation instability, by mapping physical control parameters to hydrodynamic coefficients (Rousseaux et al., 2018). Numerical solutions leverage both direct PDE solvers (finite-difference, spectral, and pseudo-spectral) and mesh-free Lagrangian methods—e.g., smoothed particle hydrodynamics (SPH)—capable of capturing vortex nucleation and lattice formation in rotating BECs, with benchmarked agreement to static and dynamic hydrodynamic behaviors (Tsuzuki, 2023).

---

References:

• (Carles et al., 2011) Madelung, Gross-Pitaevskii and Korteweg
• (Andreev, 2012) The quantum hydrodynamic description of quantum gases with different interactions
• (Jerrard et al., 2013) Hydrodynamic limit of the Gross-Pitaevskii equation
• (Kulkarni et al., 2015) Fluctuating hydrodynamics for a discrete Gross-Pitaevskii equation: mapping to Kardar-Parisi-Zhang universality class
• (Rousseaux et al., 2018) Modelling of Bose - Einstein condensation in a water tank
• (Stringari, 2018) Quantum fluctuations and Gross-Pitaevskii theory
• (Jerrard et al., 2020) Dynamics of nearly parallel vortex filaments for the Gross-Pitaevskii equation
• (Andreev, 2020) Hydrodynamics of the atomic Bose-Einstein condensate beyond the mean-field approximation: a mini-review
• (Brachet et al., 2022) Coupling Navier-Stokes and Gross-Pitaevskii equations for the numerical simulation of two-fluid quantum flows
• (Tsuzuki, 2023) Numerical model of the Gross-Pitaevskii equation for rotating Bose-Einstein condensates using smoothed-particle hydrodynamics
• (Wegner, 2023) Global-in-time Well-posedness of the One-dimensional Hydrodynamic Gross-Pitaevskii Equations without Vacuum