Pressure-Bounded Virial Equilibrium (PVE)

• Pressure-bounded virial equilibrium (PVE) is a framework that extends classical virial equilibrium by incorporating an ambient surface-pressure term to achieve quasi-static stability in finite systems.
• It is applied to various astrophysical structures such as molecular clouds, starless cores, and galaxy clusters, using metrics like velocity dispersion, mass, and external pressure to determine stability.
• PVE distinguishes gravitational binding from pressure confinement, revealing that high virial parameters may arise from external pressure support rather than indicating free expansion.

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 [(Field et al., 2011); (Myers et al., 2022); (Lopez-Corredoira et al., 2022)].

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 $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$0. The added term is compressive: it allows clouds with lower $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$1 or apparently too large $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$2 for self-gravitating equilibrium to remain virialized if the boundary pressure is nonzero (Field et al., 2011).

The same logic appears in finite-radius cluster analyses. For a subsystem bounded at $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$3, the scalar virial theorem becomes $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$4, where $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$5 is the radial pressure at the boundary and $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$6. Expressed in terms of the line-of-sight velocity dispersion, the result is

$3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$7

so the surface term is an additive correction to the usual closed-system virial estimate (Lopez-Corredoira et al., 2022).

Several implementations extend the pressure term further. In CMZ clump models, magnetic support is included through a multiplicative factor $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$8, where $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$9 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 $2T + W = 0$0, with the magnetic correction small for supercritical clouds with $2T + W = 0$1–$2T + W = 0$2 (Myers et al., 2022, Myers et al., 7 Aug 2025).

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

$2T + W = 0$3

with the same $2T + W = 0$4 dependence appearing in the uniform-density derivation. This critical mass is the maximum stable mass for fixed $2T + W = 0$5 and $2T + W = 0$6: 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: $2T + W = 0$7, so a Larson-like $2T + W = 0$8 relation can retain its slope while acquiring an environment-dependent normalization (Field et al., 2011).

The standard virial parameter remains

$2T + W = 0$9

but its interpretation changes in PVE. In the CMZ formalism, a critical column density

$M$0

defines $M$1, and the virial parameter becomes

$M$2

This yields a pressure-dominated branch at $M$3, a critical point at $M$4, and a gravity-dominated branch at $M$5. In that formulation, the gravitationally bound range is $M$6, while $M$7 on the low-column, strongly pressure-confined branch (Myers et al., 2022).

For molecular clouds with nearly constant velocity dispersion and external pressure, the pressure-bounded limit also predicts a distinctive mass–column relation. Eliminating $M$8 in favor of $M$9 and $R$0 gives

$R$1

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

$R$2

and the critical PVE state of an unmagnetized uniform sphere corresponds to $R$3. That work then uses a practical critical range $R$4–$R$5, distinguishing unstable configurations from stable pressure-confined ones (Myers et al., 7 Aug 2025).

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 (Field et al., 2011)	Clouds lie above the simple-VE line in the $R$6–$R$7 plane	Different clouds are consistent with $R$8–$R$9; no single pressure fits the sample
CMZ clumps (Myers et al., 2022)	$\sigma$0–$\sigma$1 slope near the pressure-bounded limit	755 clumps in 22 clouds; $\sigma$2; nine-cloud model gives $\sigma$3–$\sigma$4, $\sigma$5–$\sigma$6, bound fraction $\sigma$7, and typical $\sigma$8–$\sigma$9
Milky Way CO clouds (Myers et al., 7 Aug 2025)	Radial virial-parameter trend requires environmental pressure	$\Gamma$0 increases by a factor $\Gamma$1 from $\Gamma$2 to $\Gamma$3, 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 $\Gamma$4, $\Gamma$5, and $\Gamma$6, with the abstract extending that range to $\Gamma$7. 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 (Field et al., 2011).

The CMZ result is more specific. Across 22 clouds, nearly all clumps follow $\Gamma$8 with $\Gamma$9, close to the pressure-bounded prediction $3M\sigma^2 - \Gamma GM^2/R = 0$0. The nine-cloud virial models further indicate a largely unbound population: 213 clumps have $3M\sigma^2 - \Gamma GM^2/R = 0$1–$3M\sigma^2 - \Gamma GM^2/R = 0$2, mean external pressure $3M\sigma^2 - \Gamma GM^2/R = 0$3–$3M\sigma^2 - \Gamma GM^2/R = 0$4, bound fraction $3M\sigma^2 - \Gamma GM^2/R = 0$5, and typical $3M\sigma^2 - \Gamma GM^2/R = 0$6–$3M\sigma^2 - \Gamma GM^2/R = 0$7. 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 $3M\sigma^2 - \Gamma GM^2/R = 0$8, ten low-mass clumps follow a slope $3M\sigma^2 - \Gamma GM^2/R = 0$9, and 73 more massive clumps follow $4\pi P_{\rm e}R^3 \equiv 3P_{\rm e}V$0, which the authors interpret as a sequence of critically bound clumps with increasing velocity dispersion (Myers et al., 2022).

The Milky Way disk application connects PVE explicitly to environment. In two CO surveys, $4\pi P_{\rm e}R^3 \equiv 3P_{\rm e}V$1 increases by a factor $4\pi P_{\rm e}R^3 \equiv 3P_{\rm e}V$2 between 4 and 15 kpc. A fiducial fit gives $4\pi P_{\rm e}R^3 \equiv 3P_{\rm e}V$3 and $4\pi P_{\rm e}R^3 \equiv 3P_{\rm e}V$4, corresponding to $4\pi P_{\rm e}R^3 \equiv 3P_{\rm e}V$5 and $4\pi P_{\rm e}R^3 \equiv 3P_{\rm e}V$6, 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 $4\pi P_{\rm e}R^3 \equiv 3P_{\rm e}V$7 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 (Myers et al., 7 Aug 2025).

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

$4\pi P_{\rm e}R^3 \equiv 3P_{\rm e}V$8

with $4\pi P_{\rm e}R^3 \equiv 3P_{\rm e}V$9 the internal plus bulk kinetic term, $$\frac{1}{2}\ddot I = 3M\sigma^2 - \Gamma \frac{GM^2}{R} + 4\pi P_{\rm e}R^3,$$0 the gravitational binding term, $$\frac{1}{2}\ddot I = 3M\sigma^2 - \Gamma \frac{GM^2}{R} + 4\pi P_{\rm e}R^3,$$1 the external-pressure term, and $$\frac{1}{2}\ddot I = 3M\sigma^2 - \Gamma \frac{GM^2}{R} + 4\pi P_{\rm e}R^3,$$2 a magnetic contribution. The external pressure is estimated as $$\frac{1}{2}\ddot I = 3M\sigma^2 - \Gamma \frac{GM^2}{R} + 4\pi P_{\rm e}R^3,$$3 using the modeled density just outside the core and the NH$$\frac{1}{2}\ddot I = 3M\sigma^2 - \Gamma \frac{GM^2}{R} + 4\pi P_{\rm e}R^3,$$4-derived outer velocity dispersion (Scibelli et al., 2023).

The B10 sample contains 14 starless cores with central densities from $$\frac{1}{2}\ddot I = 3M\sigma^2 - \Gamma \frac{GM^2}{R} + 4\pi P_{\rm e}R^3,$$5 to $$\frac{1}{2}\ddot I = 3M\sigma^2 - \Gamma \frac{GM^2}{R} + 4\pi P_{\rm e}R^3,$$6, with mean $$\frac{1}{2}\ddot I = 3M\sigma^2 - \Gamma \frac{GM^2}{R} + 4\pi P_{\rm e}R^3,$$7. 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 ($$\frac{1}{2}\ddot I = 3M\sigma^2 - \Gamma \frac{GM^2}{R} + 4\pi P_{\rm e}R^3,$$8) are either in virial equilibrium or bound by gravity and external pressure: $$\frac{1}{2}\ddot I = 3M\sigma^2 - \Gamma \frac{GM^2}{R} + 4\pi P_{\rm e}R^3,$$9 lie in the equilibrium region and $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$00 on the bound side. The outer densities used in the pressure estimate range from $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$01 to $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$02, with mean $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$03, and the pressure term commonly exceeds the gravitational term in magnitude. The small cores f1 and f2 remain unbound even after pressure is included (Scibelli et al., 2023).

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 $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$04 to $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$05, with median $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$06, and that $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$07 can differ from $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$08 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: $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$09, 15, 16, 14, and $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$10 for cores 6, 7-1, 9, 12, and 14, respectively, summarized in the paper as an effective magnetic field difference of only $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$11 (Scibelli et al., 2023).

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 $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$12 or $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$13 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

$3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$14

with $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$15 for isotropic pressure, or

$3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$16

for anisotropic velocities. This produces a multiplicative pressure factor $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$17 in the mass–dispersion relation (Lopez-Corredoira et al., 2022).

For Newtonian gravity with dark matter and an NFW halo, the analysis gives

$3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$18

with weak concentration dependence over $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$19. For $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$20, the adopted values are $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$21 and $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$22. In MOND, using baryons only with an isothermal $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$23-model, the fitted form is

$3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$24

with $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$25–$3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$26. For $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$27 and $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$28, the paper finds $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$29, $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$30, and $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$31; for $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$32 and $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$33, it finds $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$34, $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$35, and $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$36. The sample of 178 observed clusters has empirical best fit

$3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$37

In that study, omitting pressure corrections yields MOND velocity dispersions $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$38–$3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$39 below Newton+DM for default parameters; because $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$40, that corresponds to masses $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$41–$3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$42 lower, i.e. the factor-$3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$43 discrepancy emphasized in earlier MOND work. The pressure term is therefore presented as essential to any virial analysis of a non-closed cluster subsystem (Lopez-Corredoira et al., 2022).

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 $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$44 is close to the pressure-bounded limit even though the characteristic virial parameters are $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$45–11 under the authors’ slope-based interpretation and the modeled clump population has typical $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$46–15. In the Milky Way disk analysis, clouds with $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$47 are explicitly described as stable PVE configurations rather than dispersing objects (Myers et al., 2022, Myers et al., 7 Aug 2025).

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 $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$48 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 $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$49, concluding that most clumps are unbound and cannot grow significantly before turbulence disrupts them. The inferred bound fraction is only $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$50, which is presented as a mechanism for star-formation suppression. By contrast, the Milky Way cloud study argues that many clouds with $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$51 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 $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$52 is sufficient to match the Milky Way star-formation rate (Myers et al., 2022, Myers et al., 7 Aug 2025).

A third misconception is that one external pressure should characterize an entire observational sample. The GRS cloud analysis explicitly rejects that: a single $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$53 cannot explain the data, whereas a distribution $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$54–$3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$55 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 (Field et al., 2011).

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 $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$56–$3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$57 slopes because direct $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$58 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 $3M\sigma^2 - \Gamma GM^2/R + 4\pi P_{\rm e}R^3 = 0$59, 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 (Scibelli et al., 2023, Lopez-Corredoira et al., 2022).

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.