---
title: Pressure-Bounded Virial Equilibrium (PVE)
url: https://www.emergentmind.com/topics/pressure-bounded-virial-equilibrium-pve
type: topic
---

# Pressure-Bounded Virial Equilibrium (PVE)

Pressure-bounded virial equilibrium (PVE) is the form of virial equilibrium appropriate to finite systems whose boundaries are confined by an ambient medium. Instead of balancing only internal kinetic support against self-gravity, PVE adds a surface-pressure term, so the canonical equilibrium condition is \(2T + W + 3P_{\rm ext}V = 0\), or, for an approximately spherical cloud, \(3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0\). In this framework, a structure can be gravitationally unbound in the narrow sense of simple virial equilibrium and yet remain quasi-static because the surrounding medium supplies the missing confinement. The same formal logic has been applied to Galactic molecular clouds, Central Molecular Zone (CMZ) clumps, starless cores, and virialized regions of galaxy clusters treated as finite subsystems rather than closed systems [1106.3017; 2202.13987; 2210.13961].

## 1. Formal definition and dynamical content

In simple virial equilibrium, the equilibrium condition for a self-gravitating cloud is \(2T + W = 0\). For a roughly spherical object with mass \(M\), radius \(R\), one-dimensional velocity dispersion \(\sigma\), and density-profile factor \(\Gamma\), this is \(3M\sigma^2 - \Gamma GM^2/R = 0\). PVE generalizes that balance by adding the work term of the surrounding medium, \(4\pi P_{\rm e}R^3 \equiv 3P_{\rm e}V\), so that equilibrium becomes
\[
\frac{1}{2}\ddot I = 3M\sigma^2 - \Gamma \frac{GM^2}{R} + 4\pi P_{\rm e}R^3,
\]
with equilibrium defined by \(\ddot I=0\). The added term is compressive: it allows clouds with lower \(\Sigma\) or apparently too large \(\sigma\) for self-gravitating equilibrium to remain virialized if the boundary pressure is nonzero [1106.3017].

The same logic appears in finite-radius cluster analyses. For a subsystem bounded at \(r_{\rm max}\), the scalar virial theorem becomes \(2K + V + 3P_rV_r = 0\), where \(P_r\) is the radial pressure at the boundary and \(V_r = 4\pi r_{\rm max}^3/3\). Expressed in terms of the line-of-sight velocity dispersion, the result is
\[
\sigma_{v,r<r_{\rm max}}^2
=
-\frac{V(r_{\rm max})}{3M(r_{\rm max})}
+
\frac{P_rV_r}{M(r_{\rm max})},
\]
so the surface term is an additive correction to the usual closed-system virial estimate [2210.13961].

Several implementations extend the pressure term further. In CMZ clump models, magnetic support is included through a multiplicative factor \(C_M = 1-\mu^{-2}\), where \(\mu\) is the normalized mass-to-flux ratio. In Milky Way molecular-cloud models, the pressure balance is written in terms of cloud surface density, ambient gas surface density, nearby stellar surface density, and a magnetic mass-to-flux ratio \(\lambda\), with the magnetic correction small for supercritical clouds with \(\lambda \sim 2\)–\(3\) [2202.13987; 2508.05826].

## 2. Critical states, virial parameters, and scaling relations

A central consequence of PVE is the existence of a pressure-determined critical state. For a pressure-bounded isothermal sphere, the critical mass and radius scale as
\[
M_{\rm c} = 1.2\,\frac{\sigma^4}{G^{3/2}P_{\rm e}^{1/2}},
\qquad
R_{\rm c} = 0.49\,\frac{\sigma^2}{\sqrt{GP_{\rm e}}},
\]
with the same \(\sigma^4/P_{\rm e}^{1/2}\) dependence appearing in the uniform-density derivation. This critical mass is the maximum stable mass for fixed \(\sigma\) and \(P_{\rm e}\): below it, stable equilibria exist; at it, the configuration is marginally stable; above it, no hydrostatic equilibrium exists and collapse or fragmentation follows. The same framework implies that the size–linewidth normalization is not universal: \(\sigma^2/R_{\rm c} \propto P_{\rm e}^{1/2}\), so a Larson-like \(\sigma \propto R^{1/2}\) relation can retain its slope while acquiring an environment-dependent normalization [1106.3017].

The standard virial parameter remains
\[
\alpha_{\rm vir} = \frac{5\sigma^2R}{GM},
\]
but its interpretation changes in PVE. In the CMZ formalism, a critical column density
\[
N_0 = \left[\frac{20P_s}{\pi m^2GC_M}\right]^{1/2}
\]
defines \(\nu = N/N_0\), and the virial parameter becomes
\[
\alpha = C_M\left(1+\frac{1}{3\nu^2}\right).
\]
This yields a pressure-dominated branch at \(\nu \ll 1\), a critical point at \(\nu = 1\), and a gravity-dominated branch at \(\nu \gg 1\). In that formulation, the gravitationally bound range is \(1 \le \alpha/C_M \le 2\), while \(\alpha/C_M \to \infty\) on the low-column, strongly pressure-confined branch [2202.13987].

For molecular clouds with nearly constant velocity dispersion and external pressure, the pressure-bounded limit also predicts a distinctive mass–column relation. Eliminating \(R\) in favor of \(M\) and \(N\) gives
\[
N \propto M^{1/3},
\]
which is the characteristic signature of a pressure-confined ensemble with similar mean density. In the more phenomenological Milky Way disk treatment, the mean internal pressure is
\[
P_{\rm int}=\frac{3\pi G}{20}\,\alpha_{\rm vir}\,\Sigma_{\rm cld}^2,
\]
and the critical PVE state of an unmagnetized uniform sphere corresponds to \(\alpha_{\rm vir,p}=4/3\). That work then uses a practical critical range \(\alpha_{\rm vir,crit}\sim 1\)–\(2 \approx 2\), distinguishing unstable configurations from stable pressure-confined ones [2508.05826].

## 3. Molecular-cloud and clump manifestations

The observational literature uses PVE to explain several otherwise inconsistent cloud populations. The table summarizes three representative regimes.

| System | PVE indicator | Quantitative result |
|---|---|---|
| GRS molecular clouds [1106.3017] | Clouds lie above the simple-VE line in the \(\sigma^2/R\)–\(\Sigma\) plane | Different clouds are consistent with \(P_{\rm e}/k \sim 10^4\)–\(10^7\ {\rm cm^{-3}\,K}\); no single pressure fits the sample |
| CMZ clumps [2202.13987] | \(N\)–\(M\) slope near the pressure-bounded limit | 755 clumps in 22 clouds; \(N \propto M^{0.38\pm0.03}\); nine-cloud model gives \(\sigma = 1\)–\(2\ {\rm km\,s^{-1}}\), \(P_s/k = 0.5\)–\(4\times10^8\ {\rm cm^{-3}\,K}\), bound fraction \(0.06\), and typical \(\alpha = 4\)–\(15\) |
| Milky Way CO clouds [2508.05826] | Radial virial-parameter trend requires environmental pressure | \(\alpha_{\rm vir}\) increases by a factor \(\sim 2\) from \(R_{\rm gal}=4\) to \(15\ {\rm kpc}\), and the fit requires comparable cloud and nearby-star surface densities |

In the Galactic Ring Survey interpretation, the key point is not that all clouds share one ambient pressure, but that most can be placed on PVE curves if the external pressure varies from cloud to cloud. The data occupy the region between theoretical curves for \(P_{\rm e}/k = 10^4\), \(10^5\), and \(10^6\ {\rm K\,cm^{-3}}\), with the abstract extending that range to \(10^7\ {\rm cm^{-3}\,K}\). The same analysis argues that clouds cluster near the critical line rather than filling the entire PVE-allowed region, which is consistent with the idea that many clouds sit near their pressure-defined critical mass [1106.3017].

The CMZ result is more specific. Across 22 clouds, nearly all clumps follow \(N \sim M^s\) with \(s = 0.38 \pm 0.03\), close to the pressure-bounded prediction \(s=1/3\). The nine-cloud virial models further indicate a largely unbound population: 213 clumps have \(\sigma = 1\)–\(2\ {\rm km\,s^{-1}}\), mean external pressure \(0.5\)–\(4\times10^8\ {\rm cm^{-3}\,K}\), bound fraction \(0.06\), and typical \(\alpha = 4\)–\(15\). The interpretation is that most CMZ clumps are not freely dispersing but are pressure-confined by the surrounding turbulent medium. Sgr B2 forms a contrasting regime: 43 of 55 well-measured SMA clumps host ALMA sources or H II regions, the inferred bound fraction is \(f_{b,\mathrm{SgrB2}}\sim 0.7\pm0.1\), ten low-mass clumps follow a slope \(s=0.29\pm0.04\), and 73 more massive clumps follow \(s=0.53\pm0.05\), which the authors interpret as a sequence of critically bound clumps with increasing velocity dispersion [2202.13987].

The Milky Way disk application connects PVE explicitly to environment. In two CO surveys, \(\alpha_{\rm vir}\) increases by a factor \(\sim 2\) between 4 and 15 kpc. A fiducial fit gives \(\alpha_{\rm vir,0}(4\,{\rm kpc}) \approx 1.7\) and \(\alpha_{\rm vir,0}(15\,{\rm kpc}) \approx 3.6\), corresponding to \(P_{\rm int}/P_{\rm ext} \simeq 2.3\) and \(1.4\), respectively. The model reproduces the radial trend only when the external pressure includes the stellar term; a gas-only model is nearly flat with radius, and zero external pressure gives \(\alpha_{\rm vir}=1\) everywhere. In that sense, PVE is not only a confinement model but also an environmental diagnostic, with nearby stellar surface density acting as a major control variable [2508.05826].

## 4. Starless cores, external pressure, and magnetic support

In dense-core studies, PVE is often formulated in terms of energy-like virial components. For the B10 region of Taurus, the analysis uses
\[
2\Omega_K + \Omega_G + \Omega_P + \Omega_B = 0,
\]
with \(\Omega_K\) the internal plus bulk kinetic term, \(\Omega_G\) the gravitational binding term, \(\Omega_P = -4\pi P_{\rm out}R_{\rm out}^3\) the external-pressure term, and \(\Omega_B\) a magnetic contribution. The external pressure is estimated as \(P_{\rm out} = \mu m_H n_{\rm out}\sigma_{\rm tot,out}^2\) using the modeled density just outside the core and the NH\(_3\)-derived outer velocity dispersion [2303.09574].

The B10 sample contains 14 starless cores with central densities from \(5\times10^4\) to \(1\times10^6\ {\rm cm^{-3}}\), with mean \(2.6\times10^5\ {\rm cm^{-3}}\). Ignoring magnetic fields, none of the 14 are gravitationally bound by self-gravity alone. Once the surface-pressure term is included, however, 9 of the 14 cores (\(64\%\)) are either in virial equilibrium or bound by gravity and external pressure: \(36\%\) lie in the equilibrium region and \(28\%\) on the bound side. The outer densities used in the pressure estimate range from \(0.07\times10^3\) to \(10.21\times10^3\ {\rm cm^{-3}}\), with mean \(4.45\times10^3\ {\rm cm^{-3}}\), and the pressure term commonly exceeds the gravitational term in magnitude. The small cores f1 and f2 remain unbound even after pressure is included [2303.09574].

That study also shows why detailed structure matters for PVE diagnostics. Masses and gravitational energies are computed from 3D radiative-transfer density models rather than from projected column density alone. The authors report that line-of-sight mass estimates can differ from the 3D masses by factors of \(0.3\) to \(6\), with median \(1.3\), and that \(\Omega_{G,\rm obs}\) can differ from \(\Omega_{G,\rm 3D}\) by up to a factor of 14, with median factor 2. This makes the pressure-dominated nature of many B10 cores methodologically significant: a conventional gravity-versus-kinetics virial analysis would systematically understate the role of confinement. For the five over-bound cores, only a modest effective magnetic-field difference is needed to restore exact virial balance: \(\Delta B_{\rm eff} = 13\), 15, 16, 14, and \(14\ \mu{\rm G}\) for cores 6, 7-1, 9, 12, and 14, respectively, summarized in the paper as an effective magnetic field difference of only \(\sim 15\ \mu{\rm G}\) [2303.09574].

## 5. PVE in finite-radius galaxy clusters

In clusters of galaxies, PVE appears as a boundary correction rather than as classical interstellar confinement. The virial sphere at \(r_{200}\) or \(r_{500}\) is not a closed system: galaxies and gas cross the boundary, and the matter just outside the chosen radius exerts a confining radial pressure. The relevant virial relation is therefore
\[
2K(r_{\rm max}) + V(r_{\rm max}) + 3P_r(r_{\rm max})V_r = 0,
\]
with \(P_r(r_{\rm max}) = \sigma_{v,r=r_{\rm max}}^2\rho(r_{\rm max})\) for isotropic pressure, or
\[
P_r(r_{\rm max}) = \frac{3}{3-2\beta_a}\sigma_{v,r=r_{\rm max}}^2\rho(r_{\rm max})
\]
for anisotropic velocities. This produces a multiplicative pressure factor \(F_P\) in the mass–dispersion relation [2210.13961].

For Newtonian gravity with dark matter and an NFW halo, the analysis gives
\[
\sigma_{v,r_{200},{\rm Newton+DM}}
=
A_{\rm Newton,NFW}(C)
\left(\frac{M_{500}}{10^{14}M_\odot}\right)^{1/3},
\]
with weak concentration dependence over \(2.7<C<6.4\). For \(C=3\), the adopted values are \(A_{\rm Newton,NFW}=522\ {\rm km\,s^{-1}}\) and \(F_P=1.244\). In MOND, using baryons only with an isothermal \(\beta\)-model, the fitted form is
\[
\sigma_{v,r_{200},{\rm MOND}}
\approx
A(\beta,x_{500})
\left(\frac{M_{500}}{10^{14}M_\odot}\right)^{B(\beta,x_{500})},
\]
with \(B \approx 0.29\)–\(0.30\). For \(\beta=2/3\) and \(x_{500}=0.15\), the paper finds \(A \simeq 493\ {\rm km\,s^{-1}}\), \(B \simeq 0.295\), and \(F_P=1.357\); for \(\beta=0.65\) and \(x_{500}=0.3\), it finds \(A=553\ {\rm km\,s^{-1}}\), \(B=0.294\), and \(F_P=1.448\). The sample of 178 observed clusters has empirical best fit
\[
\sigma_{v,\rm best\ fit}
=
(613\pm22)
\left(\frac{M_{500}}{10^{14}M_\odot}\right)^{0.230\pm0.027}
{\rm km\,s^{-1}}.
\]
In that study, omitting pressure corrections yields MOND velocity dispersions \(15\)–\(25\%\) below Newton+DM for default parameters; because \(M\propto \sigma_v^3\), that corresponds to masses \(40\)–\(60\%\) lower, i.e. the factor-\(\sim 2\) discrepancy emphasized in earlier MOND work. The pressure term is therefore presented as essential to any virial analysis of a non-closed cluster subsystem [2210.13961].

## 6. Interpretation, misconceptions, and limitations

A persistent misconception is that a virial parameter above the isolated-cloud threshold automatically implies free expansion. PVE changes that interpretation. In the CMZ, the observed slope \(N \propto M^{0.38\pm0.03}\) is close to the pressure-bounded limit even though the characteristic virial parameters are \(\alpha \approx 8\)–11 under the authors’ slope-based interpretation and the modeled clump population has typical \(\alpha = 4\)–15. In the Milky Way disk analysis, clouds with \(\alpha_{\rm vir}\gtrsim 2\) are explicitly described as stable PVE configurations rather than dispersing objects [2202.13987; 2508.05826].

A second misconception is that high ambient pressure should automatically imply efficient star formation. The CMZ case argues the opposite. The same turbulent environment that provides confinement also raises \(\sigma\) and keeps most clumps gravitationally unbound. The CMZ paper combines its virial analysis with a stopped-accretion model and finds accretion and dispersal times of \(\sim 0.2\ {\rm Myr}\), concluding that most clumps are unbound and cannot grow significantly before turbulence disrupts them. The inferred bound fraction is only \(f_b \approx 0.06\), which is presented as a mechanism for star-formation suppression. By contrast, the Milky Way cloud study argues that many clouds with \(\alpha_{\rm vir}>2\) still form stars through local collapse in filaments and dense cores, not through simple global contraction; in that framework, a protostellar core mass fraction of order \(10^{-3}\) is sufficient to match the Milky Way star-formation rate [2202.13987; 2508.05826].

A third misconception is that one external pressure should characterize an entire observational sample. The GRS cloud analysis explicitly rejects that: a single \(P_{\rm e}\) cannot explain the data, whereas a distribution \(P_{\rm e}/k \sim 10^4\)–\(10^7\ {\rm cm^{-3}\,K}\) can. The implication is that PVE is not a one-parameter universal law but a framework in which the confining term is itself environment dependent [1106.3017].

The main limitations recur across applications. Molecular-cloud and clump models often assume spherical or uniform-density structure, even though real objects are filamentary and hierarchical. Magnetic fields are frequently parameterized rather than measured directly. Some CMZ conclusions rely on \(N\)–\(M\) slopes because direct \(\sigma\) measurements are unavailable for most clumps. Core analyses depend sensitively on dust opacity, core-boundary definitions, and 3D structural reconstruction. Cluster applications assume spherical symmetry, a fixed anisotropy parameter, slowly varying \(\sigma_v(r)\), and specific baryon or dark-matter profiles. These caveats do not nullify the framework, but they define the regime in which PVE is an effective description rather than an exact dynamical solution [2303.09574; 2210.13961].

Across these literatures, PVE functions as a unifying statement: virial balance in real astrophysical systems is often a balance among internal motions, self-gravity, and a boundary term supplied by the surrounding medium. Its main technical value is that it distinguishes gravitational binding from confinement, and therefore separates questions of equilibrium from questions of collapse, fragmentation, and star formation.

Source: https://www.emergentmind.com/topics/pressure-bounded-virial-equilibrium-pve