---
title: Isothermal Euler–Poisson System
url: https://www.emergentmind.com/topics/isothermal-euler-poisson-system
type: topic
---

# Isothermal Euler–Poisson System

The isothermal Euler–Poisson system is a class of compressible fluid models in which Euler dynamics is coupled to a Poisson field and the pressure law is isothermal, \(p(\rho)=c_s^2\rho\). In plasma formulations, the Poisson equation describes electrostatic self-consistency and the isothermal law produces the logarithmic force \(c_s^2\nabla\log \rho\) in nonconservative variables; in astrophysical formulations, the same Euler–Poisson structure appears with the opposite, attractive Poisson sign and models self-gravitating gases. The literature therefore does not use a single canonical equation, but rather several closely related systems: single-fluid electron models with an immobile ion background, ion-acoustic models with Boltzmann electrons, two-fluid electrostatic plasma models, damped periodic models, and self-gravitating isothermal flows [1109.5980][1701.06370].

## 1. Model classes, closures, and sign conventions

The basic isothermal closure is
\[
p(n)=c_s^2\,n,\qquad \frac{\nabla p(n)}{n}=c_s^2\,\nabla\log n,
\]
and this closure is explicitly identified in the two-dimensional Euler–Poisson analysis of Li and Wu, which also notes that the same global theory applies to \(\gamma=1\) after normalization [1109.5980]. In the electron-fluid plasma model with a uniform immobile ion background of density \(n_0\), the governing equations before normalization are
\[
\partial_t n+\nabla\!\cdot(nu)=0,\qquad
m_e n(\partial_t u+u\!\cdot\nabla u)+\nabla p(n)=e n\nabla\phi,\qquad
\Delta\phi=4\pi e (n-n_0),
\]
together with the neutrality condition \(\int_{\mathbb R^2}(n-n_0)\,dx=0\) [1109.5980].

A second standard formulation is the ion-acoustic model with Boltzmann electrons. In the one-dimensional nondimensional isothermal Euler–Poisson system studied for solitary-wave stability,
\[
\partial_t n + \partial_s\big((1+n)u\big)=0,\qquad
\partial_t u + u\,\partial_s u + K\,\partial_s\log(1+n) = -\,\partial_s\phi,\qquad
-\,\partial_s^2\phi=(1+n)-e^\phi,
\]
where \(K=T_i/T_e>0\), \(1+n\) is the ion density, and the far-field equilibrium is \((n,u,\phi)\to(0,0,0)\) [2012.07687]. Closely related one-dimensional ion-acoustic models also appear in the small-amplitude solitary-wave and KdV-limit analyses, with the nondimensional equations
\[
n_t+(nu)_x=0,\qquad
u_t+u\,u_x+\sigma (n_x/n)=-\phi_x,\qquad
\phi_{xx}=e^\phi-n,
\]
where \(\sigma=T_i/T_e\) and the equilibrium is \(n\to1\), \(u\to0\), \(\phi\to0\) [1805.00197].

A third formulation is the electrostatic two-fluid model. In normalized variables around a neutral constant background, the isothermal two-fluid Euler–Poisson system in three dimensions is
\[
\partial_t n+\nabla\!\cdot[(1+n)v]=0,\qquad
\partial_t v+v\!\cdot\nabla v+T\nabla n-\nabla\phi=0,
\]
\[
\partial_t p+\nabla\!\cdot[(1+p)u]=0,\qquad
\partial_t u+u\!\cdot\nabla u+\nabla p+\nabla\phi=0,\qquad
-\Delta\phi=p-n,
\]
with \(T=T_e/T_i\) in the isothermal interpretation and irrotational initial data \(\nabla\times v^0=\nabla\times u^0=0\) [1303.1060].

Boundary-value and quasineutral studies usually employ the ion formulation with massless Boltzmann electrons. In the half-space \(\mathbb R^3_+\), Gérard-Varet, Han-Kwan, and Rousset write
\[
\partial_t n+\nabla\!\cdot(nu)=0,\qquad
\partial_t u+(u\!\cdot\nabla)u+T^i\nabla(\ln n)=-\nabla\phi,\qquad
-\varepsilon^2\Delta\phi=n-e^{-\phi},
\]
with impermeability \(u_3|_{x_3=0}=0\) and Dirichlet data \(\phi|_{x_3=0}=\phi_b\) [1112.0164].

In astrophysical usage the sign of Poisson coupling is reversed. Makino’s coordinate derivations use
\[
\partial_t\rho+\nabla\!\cdot(\rho u)=0,\qquad
\partial_t(\rho u)+\nabla\!\cdot(\rho u\otimes u)+\nabla P=-\rho\nabla\Phi,\qquad
\Delta\Phi=4\pi G\rho,
\]
with the isothermal equation of state \(P=c_s^2\rho\) [1701.06370]. In the spherically symmetric exterior-domain problem of Li, Li, and Zhu, the normalized self-gravitating system is
\[
\partial_t\rho+\nabla\!\cdot m=0,\qquad
\partial_t m+\nabla\!\cdot(m\otimes m/\rho)+\nabla p(\rho)=-\rho\nabla\Phi,\qquad
-\Delta\Phi=\rho,
\]
with \(p(\rho)=\rho\) [2308.12567]. A common misconception is that all Euler–Poisson papers study the same sign convention; the literature instead separates repulsive plasma models from attractive self-gravitating ones, and their stability mechanisms differ accordingly [1109.5980].

## 2. Linearization, dispersion, and normal forms

For the repulsive plasma problem, linearization around a homogeneous equilibrium often produces a Klein–Gordon or Klein–Gordon-type structure. In the two-dimensional electron-fluid model, after normalization and irrotational reduction \(v=\nabla\varpi\), Li and Wu introduce
\[
h(t,x)=u(t,x)-i\,|\nabla|\,\varpi(t,x),
\]
so that the linear flow is
\[
\widehat h_{\mathrm{lin}}(t,\xi)=e^{it\langle\xi\rangle}\widehat h_0(\xi),\qquad
\langle\xi\rangle=\sqrt{1+|\xi|^2}.
\]
The Klein–Gordon mass is normalized to \(1\), the linear decay in two dimensions is \(\sim t^{-1}\), and the quadratic phase
\[
\Phi_0(\xi,\eta)=\langle\xi\rangle-\langle\xi-\eta\rangle-\langle\eta\rangle
\]
has no exact time resonance, which makes a Shatah-type normal form possible [1109.5980]. Jang, Li, and Zhang formulate the same two-dimensional problem as a quasi-linear Klein–Gordon system
\[
(\Box+1)u=G_1(u,v),\qquad (\Box+1)v=G_2(u,v),
\]
with nonlocal quadratic terms coming from Poisson coupling and with the same nonintegrable \(t^{-1}\) decay as the principal obstruction [1109.3882].

Not all isothermal Euler–Poisson dispersions are Klein–Gordon in the strict sense. In the three-dimensional ion model with Boltzmann electrons, Guo and Pausader derive the linearized scalar equation
\[
\partial_t^2 p-\Delta p-\Delta(-\Delta+1)^{-1}p=0,
\]
whose Fourier frequency is
\[
\omega(\xi)^2=|\xi|^2\frac{2+|\xi|^2}{1+|\xi|^2}.
\]
This dispersion is wave-like at high frequency, nondegenerate at low frequency, and has an inflection point in \(p''(r)\); the corresponding linear group satisfies \(L^\infty\) decay \((1+|t|)^{-1}\) and \(L^{10}\) decay \((1+|t|)^{-1/8}\) in the norms used in the paper [1003.3653].

For the three-dimensional two-fluid isothermal system, Ionescu and Lie diagonalize the irrotational linearized equations into two dispersive modes \(U_e\) and \(U_i\) with dispersion relations \(\Lambda_e(\xi)\) and \(\Lambda_i(\xi)\). Both are Klein–Gordon-type in the sense that \(\Lambda_s(\xi)\approx 1+O(|\xi|^2)\) as \(|\xi|\to0\), and the linear decay is \(t^{-3/2}\) in \(L^\infty\) in three dimensions. Their nonlinear analysis further isolates null-type cancellations in the ion mode and classifies bilinear resonances into several cases, including a strong low-frequency degeneracy for the ion-like dispersion that is compensated by derivative structure [1303.1060].

A recurring analytical point is therefore that the isothermal law does not determine the dispersive class by itself. In the electron-fluid model with immobile ions, Poisson coupling produces an exact Klein–Gordon mass term after normalization; in ion-acoustic and two-fluid models, the resulting dispersions are more varied, but still sufficiently structured to permit normal forms, resonance analysis, or asymptotic reductions [1109.5980][1003.3653].

## 3. Cauchy theory: global existence, scattering, and damping

The most complete small-data result in two dimensions is Li and Wu’s solution of the full small irrotational Cauchy problem for the repulsive Euler–Poisson system. For \(N\ge 8\), \(N'=N-\tfrac12\), and \(0<\delta<\delta_*\), they define a control norm \(X_T\) combining high Sobolev regularity, weighted dispersive bounds, and an \(L_t^1L_x^{2+\delta}\) component. If
\[
\|e^{it\langle\nabla\rangle}h_0\|_{X_0}\le \varepsilon,
\]
then there is a unique global smooth solution
\[
h\in C([0,\infty);H^N)\cap C([0,\infty);H^{N'})
\]
with \(\|h\|_{X_\infty}\lesssim\varepsilon\), and there exists \(h_\infty\in H^{N'}\) such that \(e^{-it\langle\nabla\rangle}h(t)\to h_\infty\) in \(H^{N'}\). The method combines diagonalization, normal forms converting quadratic interactions into cubic ones, hidden-derivative identities, phase-derivative transformations, multiscale frequency decompositions, and Coifman–Meyer control of the nonlocal Riesz multipliers. The paper states explicitly that the analysis extends to arbitrary \(\gamma\ge1\), including the isothermal case \(\gamma=1\) [1109.5980].

An earlier two-dimensional result by Jang, Li, and Zhang constructed global smooth solutions in a final-data formulation. Their analysis subtracts a resonant correction built from the linear self-interaction, localizes the quasi-linear part in time, and then uses a modified energy scheme with an infinite-time Gronwall argument. The paper emphasizes that the same quasi-linear Klein–Gordon mechanism persists under the isothermal specialization \(p(n)=c_s^2 n\) [1109.3882].

In three dimensions, the isothermal two-fluid problem admits a global small-data theory for localized irrotational perturbations of a neutral background. Theorem D.1 of Ionescu and Lie fixes \(T\in[1/10,100]\), \(N_0=104\), and assumes
\[
\|(n^0,v^0,p^0,u^0)\|_{H^{N_0}}+\|(n^0,v^0,p^0,u^0)\|_Z\le \varepsilon_0,
\]
together with \(\nabla\times v^0=\nabla\times u^0=0\). The resulting solution is global in \(H^{N_0}\), remains irrotational for all time, and satisfies the decay estimate
\[
\sup_{|\alpha|\le 4}(1+t)^{1+\beta/2}\|D^\alpha(n(t),v(t),p(t),u(t))\|_{L^\infty}\le \varepsilon_0,\qquad \beta=\frac1{100},
\]
together with scattering in the diagonalized variables [1303.1060].

For the three-dimensional ion-acoustic system with Boltzmann electrons, Guo and Pausader prove global smooth irrotational small-amplitude solutions in Sobolev spaces \(H^{2k+1}\times H^{2k+2}\), \(k\ge5\). Their proof relies on a normal form with singular phase denominators, refined bilinear multiplier estimates, and an energy functional augmented by an \(H^{-1}\)-component to absorb the most singular normal-form contribution [1003.3653].

These dispersive small-data theories do not exhaust the isothermal regime. On the periodic torus \(\mathbb T^d\), \(2\le d\le4\), the damped isothermal Euler–Poisson system
\[
\partial_t \rho+\nabla\!\cdot(\rho u)=0,\qquad
\partial_t(\rho u)+\nabla\!\cdot(\rho u\otimes u)+\nabla\rho+\rho\nabla\phi=-\nu\rho u,\qquad
-\Delta\phi=\rho-1
\]
admits unique global smooth solutions for arbitrarily large \(H^{m+1}\) data with density bounded away from vacuum, provided \(\nu\) is sufficiently large. The solutions relax exponentially to \((\rho,u,\phi)=(1,0,0)\) at rate \(\theta=c\min\{\nu,\nu^{-1}\}\), and in the large-damping regime \(\theta=c/\nu\) [2606.04693]. This is a distinct, large-data, non-dispersive mechanism based on damping, hidden parabolicity in the density equation, and comparison with a drift–diffusion–Poisson system.

## 4. Solitary waves, embedded neutral modes, and Evans-function stability

One-dimensional isothermal Euler–Poisson equations support solitary waves in the super–ion–sonic regime. For the ion/Boltzmann-electron model with ion sound speed \(\mathsf V=\sqrt{1+K}\), the traveling coordinate \(x=s-ct\), and \(\varepsilon=c-\mathsf V>0\), Kodama, Miki, and Nakanishi show that for each \(\varepsilon\in(0,\varepsilon_K)\) there exists a unique, up to translation, even solitary wave \((n_c,u_c,\phi_c)\) with \(\partial_x n_c>0\), \(\partial_x u_c>0\), and \(\partial_x\phi_c>0\) on \((-\infty,0)\), and with exponential decay as \(|x|\to\infty\). The profiles satisfy
\[
(1+n)(c-u)=c,\qquad
\phi=H(n,c):=\frac{c^2}{2}\Big(1-\frac{1}{(1+n)^2}\Big)-K\ln(1+n),
\]
and admit a Hamiltonian first integral through the function \(g(n,c)\) [2012.07687].

Linearization about such a solitary wave produces a nonlocal operator \(\mathcal L\) whose spectral structure is delicate because \(\lambda=0\) is embedded in the essential spectrum in unweighted \(L^2\). Translation invariance and the speed parameter generate the Jordan chain
\[
\mathcal L\,\partial_x(n_c,u_c)^T=0,\qquad
\mathcal L\,\partial_c(n_c,u_c)^T=-\,\partial_x(n_c,u_c)^T,
\]
so the geometric multiplicity is one and the algebraic multiplicity is at least two [2012.07687]. In unweighted \(L^2\), the essential spectrum fills the imaginary axis; in exponentially weighted spaces \(L^2_\beta\), conjugation by the weight moves the essential spectrum strictly into the left half-plane when \(c>\sqrt{1+K}\). On the weighted spectral domain \(\Omega^\beta\), the Evans function is analytic and vanishes only at \(\lambda=0\), with order two, for sufficiently small amplitude. The same paper proves spectral stability in \(L^2\), and in weighted space it obtains asymptotic linear stability modulo the two neutral modes:
\[
\|e^{\mathcal Lt}(\dot n_0,\dot u_0)^T\|_{(L^2_\beta)^2}\le C_1e^{-C_2 t}\|(\dot n_0,\dot u_0)^T\|_{(L^2_\beta)^2}
\]
for data orthogonal to the spectral projection onto \(\lambda=0\) [2012.07687].

The small-amplitude structure of solitary waves is tied rigorously to KdV. Chae and Hwang study the stretched moving frame
\[
\xi=\varepsilon^{1/2}(x-(V+\gamma)t),\qquad V=\sqrt{1+\sigma},
\]
and prove existence of smooth solitary waves for the one-dimensional Euler–Poisson system that converge to the KdV soliton
\[
n_{\mathrm{KdV}}(\xi)=\frac{3\gamma}{V}\,\mathrm{sech}^2\!\Big(\sqrt{\frac{V\gamma}{2}}\,\xi\Big).
\]
For every integer \(k\ge0\), the remainder satisfies an exponentially weighted \(C^k\) bound of order \(O(\varepsilon^2)\) [1805.00197]. In the Evans-function language of the stability paper, the long-wave/small-amplitude scaling also leads to uniform convergence of the Euler–Poisson Evans function to the explicit KdV Evans function, which vanishes only at the origin with order two [2012.07687]. The isothermal solitary-wave problem is thus one of the clearest settings in which profile construction, spectral stability, and asymptotic reduction to KdV are simultaneously rigorous.

## 5. Singular limits, asymptotic reductions, and boundary layers

A major branch of the theory concerns singular limits in which the isothermal Euler–Poisson system approaches a reduced dispersive, diffusive, or quasineutral model. In the quasineutral limit on \(\mathbb R^3_+\), letting \(\varepsilon\to0\) in
\[
\partial_t n+\nabla\!\cdot(nu)=0,\qquad
\partial_t u+(u\!\cdot\nabla)u+T^i\nabla(\ln n)=-\nabla\phi,\qquad
-\varepsilon^2\Delta\phi=n-e^{-\phi},
\]
formally yields neutrality \(n=e^{-\phi}\) and the limiting isothermal Euler system
\[
\partial_t n+\nabla\!\cdot(nu)=0,\qquad
\partial_t u+(u\!\cdot\nabla)u+(T^i+1)\nabla(\ln n)=0,
\]
with only the impermeability condition \(u_3|_{x_3=0}=0\) at the boundary [1112.0164]. Because the Dirichlet datum \(\phi|_{x_3=0}=\phi_b\) is generally incompatible with neutrality, Gérard-Varet, Han-Kwan, and Rousset construct a Debye sheath of thickness \(O(\varepsilon)\) through a high-order boundary-layer expansion. The leading profile \(\Phi^0\) solves
\[
\partial_{zz}\Phi^0+\Gamma n^0\,S(\Phi^0)=0,\qquad
S(\Phi)=e^{-\Phi}-e^{\Phi/T^i},
\]
with decay as \(z\to+\infty\). They prove \(L^2\)-convergence to the limiting Euler flow with rate \(O(\sqrt\varepsilon)\), and sharper convergence to the full boundary-layer corrected approximation [1112.0164].

The supersonic outflow case requires a different stabilization mechanism. In the half-space formulation with
\[
\varepsilon^2\Delta\phi^\varepsilon+e^{-\phi^\varepsilon}=n^\varepsilon,
\]
Gérard-Varet, Han-Kwan, and Rousset analyze the Bohm regime \(u_3^0(y,0)<-\sqrt{T_i+1}\), where the limiting Euler characteristics are all outgoing and the boundary layer includes an \(O(1)\) normal-velocity correction. The Debye sheath is again of thickness \(O(\varepsilon)\), but its stability is obtained through Goodman-type weights adapted to outward convection, yielding convergence to the Euler solution plus explicit layer correctors [1403.1189].

In one space dimension, the long-wave/small-amplitude scaling
\[
x'=\varepsilon^{1/2}(x-Vt),\qquad t'=\varepsilon^{3/2}t
\]
reduces the isothermal ion-acoustic Euler–Poisson system to KdV. Wang proves that the leading amplitude \(n^{(1)}\) satisfies
\[
\partial_t n^{(1)}+V\,n^{(1)}\partial_x n^{(1)}+c_2\,\partial_x^3 n^{(1)}=0,
\]
with \(V^2=(T_e+T_i)/M\) and \(c_2=\lambda_D^2/(2V)\), and establishes global-in-time convergence of the Euler–Poisson solution to the KdV asymptotics in \(H^2\), with an \(O(\varepsilon^3)\) remainder in the slow variables [1202.1830]. Under a uniform magnetic field, Colin and Lannes derive the Zakharov–Kuznetsov equation from the magnetized Euler–Poisson system, first for cold plasma and then with the isothermal pressure term, while also proving local well-posedness for the Cauchy problem in dimensions two and three [1205.5080].

The large-friction limit produces a different reduced model. On \(\mathbb T^d\), \(2\le d\le4\), the damped isothermal Euler–Poisson equations admit the slow-time rescaling
\[
s=t/\nu,\qquad
\rho_\nu(s,x)=\rho^\nu(\nu s,x),\qquad
J_\nu(s,x)=\nu\rho^\nu u^\nu(\nu s,x),
\]
under which the limit \(\nu\to\infty\) is the drift–diffusion–Poisson system
\[
\partial_s\bar\rho-\nabla\!\cdot(\bar\rho\nabla\bar\Phi)=\Delta\bar\rho,\qquad
-\Delta\bar\Phi=\bar\rho-1.
\]
The density converges with rate \(O(\nu^{-2})\), and after subtracting a damped-heat initial layer the rescaled flux converges with the same rate to the drift–diffusion flux \(\bar J=-\nabla\bar\rho-\bar\rho\nabla\bar\Phi\) [2606.04693].

## 6. Singularity formation and self-gravitating isothermal flows

The isothermal closure does not imply global smoothness in every geometry or sign regime. For the one-dimensional ion/Boltzmann-electron model
\[
\partial_t\rho+\partial_x(\rho u)=0,\qquad
\partial_t u+u\,\partial_x u+K\,\partial_x(\log\rho)=-\partial_x\phi,\qquad
-\partial_{xx}\phi=\rho-e^\phi,
\]
Hou, Huang, and Li give a constructive proof of finite-time cusp formation from smooth data near a Burgers self-similar blow-up profile. In Riemann variables \(w=u+\sqrt K\log\rho\), \(z=u-\sqrt K\log\rho\), they build a renormalized dynamics around the Burgers profile \(W\) and show that the solution can develop a \(C^1\) singularity in finite time while remaining \(C^{1/3}\) at the blow-up time. More precisely, for any \(\beta\le1/3\),
\[
\sup_{t<T^*}\|(\rho,u)(\cdot,t)\|_{C^\beta}<\infty,
\]
whereas for any \(\beta>1/3\),
\[
\|\rho(\cdot,t)\|_{C^\beta}\sim (T^*-t)^{-(3\beta-1)/2},\qquad
\|u(\cdot,t)\|_{C^\beta}\sim (T^*-t)^{-(3\beta-1)/2}
\]
as \(t\uparrow T^*\) [2405.02557]. This result is specific to a one-dimensional ion-acoustic regime, but it shows that electrostatic coupling does not universally suppress finite-time singularity formation.

For attractive self-gravity, the mathematical emphasis shifts from dispersive stabilization to geometric structure, entropy solutions, and equilibrium configurations. Makino gives coordinate-system derivations of the self-gravitating isothermal Euler–Poisson equations in Cartesian, cylindrical, and spherical coordinates, together with the enthalpy \(u_{\text{enth}}(\rho)=c_s^2\ln\rho\) and the hydrostatic reduction leading to the isothermal Lane–Emden equation
\[
\frac{1}{\xi^2}\frac{d}{d\xi}\!\left(\xi^2\frac{d\psi}{d\xi}\right)=e^{-\psi}
\]
under the scaling \(\rho(r)=\rho_c e^{-\psi(r)}\) [1701.06370]. In the exterior of a ball, Li, Li, and Zhu prove global existence of spherically symmetric entropy solutions for the isothermal self-gravitating problem by a fractional-step Lax–Friedrichs scheme and compensated compactness. In weighted radial variables \(Q=r^{N-1}\rho\), \(w=r^{N-1}m\), they obtain global bounds
\[
0\le Q(r,t)\le C(T),\qquad
|w(r,t)|\le Q(r,t)\big(C(T)+|\ln Q(r,t)|\big)
\]
for almost every \((r,t)\), allowing vacuum and shocks [2308.12567].

Rotation can also regularize the attractive isothermal flow. Kwong and Yuen construct exact periodic solutions of the two-dimensional self-gravitating isothermal Euler–Poisson equations with
\[
\rho(t,r)=a(t)^{-2}\exp(f(r/a(t))),\qquad
u(t,x)=\frac{\dot a(t)}{a(t)}x+\frac{\xi}{a(t)^2}x^\perp,
\]
where \(a(t)\) solves the Emden-type ODE
\[
a''(t)=-\,\Lambda/a(t)+\xi^2/a(t)^3.
\]
For \(\Lambda>0\) and \(\xi\neq0\), the solutions are nontrivially periodic unless \(a(0)=|\xi|/\sqrt\Lambda\) and \(\dot a(0)=0\), in which case they are steady. The paper interprets the \(\xi^2/a^3\) term as the rotational mechanism preventing the blowup found earlier in non-rotating two-dimensional isothermal solutions [1408.0716].

Taken together, these results show that the phrase “isothermal Euler–Poisson system” names a structurally coherent but mathematically diverse family. In repulsive plasma settings it supports Klein–Gordon dispersion, scattering, Evans-function stability, KdV and Zakharov–Kuznetsov asymptotics, and overdamped drift–diffusion limits; in attractive settings it supports boundary layers, hydrostatic Lane–Emden equilibria, entropy solutions with vacuum, periodic rotating states, and, in some one-dimensional plasma regimes, constructive cusp singularities [1109.5980][2012.07687][1112.0164][2405.02557].

Source: https://www.emergentmind.com/topics/isothermal-euler-poisson-system