Void-Shoring: Boundary-Dominated Dynamics

• Void-Shoring is the study of how void boundaries and environmental interactions stabilize and regulate empty regions across fields like cosmology and mechanics.
• It analyzes the interface where contact forces, pressure balances, or reinforcement compensate for mass or rigidity deficits, impacting structural dynamics.
• Applications range from cosmic void evolution and elastic cloaking in plates to geological fold voiding and energetic stabilization in materials.

Searching arXiv for relevant papers on "Void-Shoring" and closely related void studies to ground the article in published work. Across the cited literature, “void-shoring” is used, or invoked by close analogy, for mechanisms by which a void is bounded, supported, compensated, or dynamically regulated by its surroundings. In cosmology, this concerns the coupling between underdense regions and the surrounding cosmic web, the geometry of void boundaries, and the use of void observables for parameter inference. In continuum and structural mechanics, it refers more literally to support or compensation around an opening, as in geological fold voiding or reinforcement-based neutralization of holes in elastic plates. In materials and astrophysical gas dynamics, it denotes the energetic or pressure balance that stabilizes, grows, or contains cavities. This suggests a cross-domain concept centered on boundary conditions, environmental backreaction, and the nontrivial mechanics of “empty” regions rather than emptiness alone.

1. Conceptual range and recurring mechanics

A recurring feature of the literature is that the void itself is not the complete object of analysis. The decisive physics lies at the interface between the void and its environment. In cosmic-void studies, the relevant interface is the wall-filament-halo ridge delimiting an underdensity; in fold mechanics it is the free boundary between a detached elastic layer and a rigid obstacle; in elastic cloaking it is a reinforcing frame that restores the mass, moment of inertia, and stiffness removed by a hole; in graphene rupture it is the void perimeter whose line tension competes with stress relief; and in self-similar gas dynamics it is a pressure-balanced contact discontinuity separating a low-mass cavity from a massive envelope (Cautun et al., 2015).

This family resemblance is clearest when the void is treated as an interface problem rather than a volume deficit. The geological model of an elastic layer over a V-shaped obstacle formulates voiding as an obstacle problem with a free boundary and yields a unique symmetric non-contact interval (Dodwell et al., 2011). The structural-dynamics model for flexural cloaking treats a void as a deficit in inertia and rigidity that must be compensated by stiffening and added masses (Misseroni et al., 2020). The astrophysical cavity model treats the void boundary as a moving contact discontinuity whose velocity matches the flow and whose pressure is continuous across the interface (Lou et al., 2011). A plausible implication is that “void-shoring” names a class of boundary-dominated problems whose governing variables differ by field, but whose logic is homologous.

2. Cosmic-void boundaries, environments, and shell dynamics

The most explicit cosmological formulation of void-shoring is the claim that void evolution is dynamically coupled to the surrounding large-scale structure rather than determined by the interior alone. In “Clues on void evolution III,” voids are classified into R-type and S-type systems using the integrated density contrast at $r_{\max}\approx 3R_{\rm void}$: $\Delta(r_{\max})<0$ for R-type and $\Delta(r_{\max})>0$ for S-type. Spherical shells from $R_{\rm void}$ to $3R_{\rm void}$, pixelized with HEALPix, show that R-type voids exhibit mostly coherent outflow, whereas S-type voids show a patchy angular mixture of infall and outflow. The linear comparison uses

$$V_{\rm lin}(r) = -\frac{H}{3}\, r\, \Delta_{\rm m}(r)\,\Omega_{\rm m}^{0.6},$$

and the measured departures depend on tracer density, shell morphology, and environment (Ruiz et al., 2015).

Several quantitative conclusions follow. Velocity smoothness increases with void size, indicating that laminar flow dominates for voids larger than about $18\,h^{-1}\,\mathrm{Mpc}$. High-density shell structures show the largest deviations from linear theory, and elongated structures depart more strongly than compact ones. The paper identifies a transition regime around $R_{\rm void}\approx 18\,h^{-1}\,\mathrm{Mpc}$, where linear theory changes from overpredicting to underpredicting the shell velocity depending on environment and tracer choice. The stated interpretation is that nonlinearities are not produced only by inner underdensity, but by scale coupling between the void interior and the surrounding structures (Ruiz et al., 2015).

A complementary reformulation appears in “The view from the boundary,” which argues that voids should be stacked relative to their boundary rather than a chosen center because cosmic voids are highly non-spherical and store most of their displaced mass at the edge. The signed boundary distance is defined by the minimum Euclidean distance to the boundary, negative inside the void and positive outside. With this coordinate, the density profile separates into a low-density interior, a sharp overdense ridge at $\mathrm{dist}=0$, and a gradual relaxation outside. The boundary profile is nearly self-similar when rescaled by the ridge thickness $t_{\rm in}$, and the weak-lensing shear and convergence increase by a factor of two relative to conventional spherical stacking (Cautun et al., 2015).

The dynamical content of the boundary formulation is equally important. The boundary velocity profile shows outflow in the inner void, a rapid turnover to inflow into the filaments and walls, and a size-dependent amplitude of the interior outflow. Small voids are, on average, contracting, large voids are, on average, expanding, and the transition occurs around $\Delta(r_{\max})<0$0. Size alone is not sufficient for classifying an individual object, but boundary velocities provide a direct criterion: positive boundary velocity indicates expansion, negative boundary velocity indicates contraction (Cautun et al., 2015).

3. Void-in-cloud, geometric systematics, and methodological controversies

A central controversy in cosmological void-shoring concerns the strength of the void-in-cloud mechanism. In the Sheth–van de Weygaert two-barrier excursion-set model, voids can either merge within larger underdensities or be crushed when embedded in overdense environments. The standard spherical thresholds are $\Delta(r_{\max})<0$1 for void shell crossing and $\Delta(r_{\max})<0$2 for collapse. “Void Formation: Does the Void-in-Cloud Process Matter?” tests this picture with five high-resolution $\Delta(r_{\max})<0$3CDM $\Delta(r_{\max})<0$4-body simulations spanning box sizes from $\Delta(r_{\max})<0$5 to $\Delta(r_{\max})<0$6 and finds that the measured void abundance below $\Delta(r_{\max})<0$7 exceeds the standard prediction by up to 2 orders of magnitude. A one-barrier Press–Schechter-like model using $\Delta(r_{\max})<0$8, or a naive power law for aspherical voids, gives a much better description of the small-void abundance (Chan et al., 2019).

The same paper argues that small voids are not generally “crushed out of existence.” Density and velocity profiles indicate that many are partially collapsing underdensities rather than fully eliminated structures. Environmental classification by the tidal tensor shows that many small voids live in filaments and even clusters, while the data are interpreted as being dominated by the void-in-void effect, with void-in-cloud only weakly influencing the abundance, even in dense environments. This directly challenges the strong version of the standard shoring-up picture and supports the Eulerian partial-collapse interpretation associated with Paranjape et al. (Chan et al., 2019).

A separate but related methodological issue is the treatment of void geometry. “Significance of void shape: Neutrino mass from Voronoi void halos?” revisits the proposed void-halo mass function (VHMF), introduced by Zhang et al. as a “void-shoring” analogue of the ordinary halo mass function. The earlier construction selected halos within a sphere of radius $\Delta(r_{\max})<0$9, where $\Delta(r_{\max})>0$0 is the effective radius of a VIDE void. The 2024 reanalysis compares this spherical-cut VHMF with a full-shape VHMF defined by the actual watershed void boundary and finds that a single anomalously non-spherical void in the $\Delta(r_{\max})>0$1 eV MassiveNuS realization drives the apparent extra neutrino signal. That void has radius about $\Delta(r_{\max})>0$2, ellipticity $\Delta(r_{\max})>0$3, contributes $\Delta(r_{\max})>0$4 of the total void volume, and dominates the spherical sample while barely affecting the full-shape sample. After removing it, there is no evidence that the VHMF contains information beyond the global halo mass function for $\Delta(r_{\max})>0$5 (Bayer et al., 2024).

The same study then proposes the VorHMF, which bins halos by the Voronoi cell volume $\Delta(r_{\max})>0$6 of their local environment rather than by membership in a grouped watershed void. Larger $\Delta(r_{\max})>0$7 corresponds to a more isolated, emptier environment, and the ratio between the $\Delta(r_{\max})>0$8 eV and $\Delta(r_{\max})>0$9 eV cases departs more strongly from unity as $R_{\rm void}$0 increases. This is presented as a more local, fine-grained summary statistic and as a more robust alternative to grouped-void definitions, precisely because it reduces sensitivity to pathological void geometries (Bayer et al., 2024).

4. Precision cosmology with voids: abundances, degeneracy breaking, and lensing

Void-shoring in cosmological inference is most explicit in the use of voids as complementary probes rather than stand-alone counts. “Cluster-Void Degeneracy Breaking” combines the largest claimed galaxy cluster, ACT-CL J0102-4915 (“El Gordo”), with the Cold-Spot void (“CS Void”) and shows that cluster and void abundances generate mutually orthogonal degeneracies in the $R_{\rm void}$1–$R_{\rm void}$2 plane. The abundance model is written in terms of the expected number of objects above a threshold observable, and the extreme-value likelihood for the maximum object in a Poisson sample is

$R_{\rm void}$3

Using the Watson et al. halo mass function for clusters and the Sheth–van de Weygaert excursion-set void mass function for voids, together with log-normal measurement models for Eddington bias, the joint flat-$R_{\rm void}$4CDM inference gives

$R_{\rm void}$5

The paper states that the joint analysis detects dark energy at $R_{\rm void}$6 using only these two extreme objects and emphasizes the complementarity of clusters and voids in scale, density, and non-linearity (Sahlén et al., 2015).

The physical reason for the degeneracy breaking is that cluster scale is inferred from mass and collapse physics, whereas void scale is inferred from observed angular and radial extent and depends geometrically on $R_{\rm void}$7. Clusters are more sensitive to growth at higher redshift, while a void such as the CS Void is sensitive to low-redshift growth. The inferred self-calibrated effective dark-matter void radius is $R_{\rm void}$8, and the authors argue that future surveys such as eROSITA, DES, DESI, Euclid, and LSST could improve constraints by factors of $R_{\rm void}$9–$3R_{\rm void}$0 (Sahlén et al., 2015).

A second observational strand is Void $3R_{\rm void}$1 CMB lensing. Using validated Roman mock catalogs, “Towards precision cosmology with Void x CMB correlations (II)” finds that the lensing signal is less sensitive to mock-catalog construction than galaxy and void statistics. The highest signal-to-noise is achieved for 2D voids with rescaled profiles. The forecast values are $3R_{\rm void}$2 for 2D voids and $3R_{\rm void}$3 for 3D voids when Roman is combined with Planck, increasing to $3R_{\rm void}$4 and $3R_{\rm void}$5 for SO, and $3R_{\rm void}$6 and $3R_{\rm void}$7 for CMB-S4-like surveys. The paper is explicit that the cosmological dependence of the observable remains to be quantified, but it positions Void $3R_{\rm void}$8 CMB lensing as a route toward direct cosmological constraints (Sar et al., 21 May 2026).

These two strands are methodologically connected. Boundary-based stacking sharpens void lensing by aligning the physically relevant ridge (Cautun et al., 2015), while cluster–void abundance complementarity exploits the fact that overdense and underdense extremes respond differently to the same background cosmology (Sahlén et al., 2015). A plausible implication is that cosmological void-shoring is most powerful when voids are treated as boundary- and environment-sensitive objects rather than as spherical deficits defined only by size.

5. The Local Void as a dynamical example

The nearby universe provides a concrete case in which void-shoring is observable as both topology and kinematics. “Cosmicflows-3: Cosmography of the Local Void” reconstructs the local density and velocity field from 18,000 CF3 galaxy distances using Wiener filtering with constrained realizations and a Bayesian/MCMC reconstruction. The inferred density contrast is related to the velocity divergence by

$3R_{\rm void}$9

The Local Void begins only about 1 Mpc from the Local Group, subtends about 40% of the sky, and is bounded by major structures including the Perseus-Pisces filament, the Norma–Pavo–Indus/Great Attractor complex, the Local Supercluster/Virgo region, and the Southern Supercluster/Fornax–Eridanus region. At density level $$V_{\rm lin}(r) = -\frac{H}{3}\, r\, \Delta_{\rm m}(r)\,\Omega_{\rm m}^{0.6},$$0, its approximate dimensions are $$V_{\rm lin}(r) = -\frac{H}{3}\, r\, \Delta_{\rm m}(r)\,\Omega_{\rm m}^{0.6},$$1, $$V_{\rm lin}(r) = -\frac{H}{3}\, r\, \Delta_{\rm m}(r)\,\Omega_{\rm m}^{0.6},$$2, and $$V_{\rm lin}(r) = -\frac{H}{3}\, r\, \Delta_{\rm m}(r)\,\Omega_{\rm m}^{0.6},$$3 Mpc, with volume roughly $$V_{\rm lin}(r) = -\frac{H}{3}\, r\, \Delta_{\rm m}(r)\,\Omega_{\rm m}^{0.6},$$4 Mpc$$V_{\rm lin}(r) = -\frac{H}{3}\, r\, \Delta_{\rm m}(r)\,\Omega_{\rm m}^{0.6},$$5 (Tully et al., 2019).

The same reconstruction emphasizes that the Local Void is not an isolated cavity. It has multiple minima, is connected through low-density passages to the Hercules and Sculptor voids, and contains thin filaments crossing its interior, including a sparse Virgo–Perseus chain and the Pegasus Cloud/Pegasus Spur. The V-web analysis confirms that voids are regions of expansion along three axes and that sheets, filaments, and knots delimit their topology. The Local Void contributes a deviant motion of roughly $$V_{\rm lin}(r) = -\frac{H}{3}\, r\, \Delta_{\rm m}(r)\,\Omega_{\rm m}^{0.6},$$6–$$V_{\rm lin}(r) = -\frac{H}{3}\, r\, \Delta_{\rm m}(r)\,\Omega_{\rm m}^{0.6},$$7 to the Local Group, and together with Virgo’s pull accounts for about 50% of the Local Group’s motion in the CMB rest frame (Tully et al., 2019).

Direct kinematic confirmation of Local Void drainage comes from “Draining the Local Void.” The dwarf galaxies ALFAZOAJ1952+1428 and KK246 have secure HST TRGB distances of $$V_{\rm lin}(r) = -\frac{H}{3}\, r\, \Delta_{\rm m}(r)\,\Omega_{\rm m}^{0.6},$$8 Mpc and $$V_{\rm lin}(r) = -\frac{H}{3}\, r\, \Delta_{\rm m}(r)\,\Omega_{\rm m}^{0.6},$$9 Mpc, respectively. With $18\,h^{-1}\,\mathrm{Mpc}$0, their inferred line-of-sight peculiar velocities are $18\,h^{-1}\,\mathrm{Mpc}$1 and $18\,h^{-1}\,\mathrm{Mpc}$2, respectively, in the Local Sheet frame. Because the Milky Way itself moves away from the Local Void by about $18\,h^{-1}\,\mathrm{Mpc}$3 to $18\,h^{-1}\,\mathrm{Mpc}$4, these negative peculiar velocities correspond to substantially larger outward motions from the void center: roughly $18\,h^{-1}\,\mathrm{Mpc}$5 for ALFAZOAJ1952+1428 and $18\,h^{-1}\,\mathrm{Mpc}$6 for KK246 (Rizzi et al., 2016).

The numerical action model used in that work includes 1385 entities within 38 Mpc and is embedded in a Wiener Filter density field. It reproduces the general downward motion in supergalactic $18\,h^{-1}\,\mathrm{Mpc}$7 and supports the interpretation that both galaxies participate in the general drainage of the Local Void, although it slightly underpredicts their observed void-directed motions (Rizzi et al., 2016). In this local setting, void-shoring is not static support but the dynamic relation between an underdense basin, the filamentary structures that delimit it, and the peculiar velocities generated by evacuation.

6. Geological, structural, materials, and astrophysical formulations

In geological folding, void-shoring is literal. “Self-similar voiding solutions of a single layered model of folding rocks” studies an elastic layer with vertical displacement $18\,h^{-1}\,\mathrm{Mpc}$8 constrained above a rigid V-shaped obstacle $18\,h^{-1}\,\mathrm{Mpc}$9. The total potential energy is the sum of bending energy and work against overburden pressure,

$R_{\rm void}\approx 18\,h^{-1}\,\mathrm{Mpc}$0

The minimizer is convex, symmetric, and has one detached interval $R_{\rm void}\approx 18\,h^{-1}\,\mathrm{Mpc}$1; the free boundary $R_{\rm void}\approx 18\,h^{-1}\,\mathrm{Mpc}$2 is determined by the load and stiffness rather than prescribed. The central scaling law is

$R_{\rm void}\approx 18\,h^{-1}\,\mathrm{Mpc}$3

so larger overburden pressure suppresses voiding while larger bending stiffness permits larger voids (Dodwell et al., 2011).

In structural dynamics, “Omnidirectional flexural invisibility of multiple interacting voids in vibrating elastic plates” gives a different but closely related meaning. A void removes mass, moment of inertia, bending stiffness, and torsional stiffness. The proposed cloak restores these missing contributions by a reinforcing frame and added masses so that the perforated region behaves dynamically like the intact plate. The method is not transformation-based cloaking but direct compensation by stiffness matching and inertia matching. For multiple voids, the reported scattering-reduction coefficient reaches near-perfect values corresponding to reductions of $R_{\rm void}\approx 18\,h^{-1}\,\mathrm{Mpc}$4, $R_{\rm void}\approx 18\,h^{-1}\,\mathrm{Mpc}$5, and $R_{\rm void}\approx 18\,h^{-1}\,\mathrm{Mpc}$6 in selected examples, and the plate with cloaked voids shares the first four vibration modes of the intact plate (Misseroni et al., 2020).

In energetic materials, void-shoring appears as the dependence of hotspot formation on void geometry and neighborhood effects. “Modeling meso-scale energy localization in shocked HMX, Part II” introduces modifier functions for ignition and growth rates that depend on void aspect ratio $R_{\rm void}\approx 18\,h^{-1}\,\mathrm{Mpc}$7, orientation $R_{\rm void}\approx 18\,h^{-1}\,\mathrm{Mpc}$8, and void fraction $R_{\rm void}\approx 18\,h^{-1}\,\mathrm{Mpc}$9. Elongated voids with $\mathrm{dist}=0$0 and acute orientations $\mathrm{dist}=0$1 strongly amplify ignition, while increasing $\mathrm{dist}=0$2 from $\mathrm{dist}=0$3 to $\mathrm{dist}=0$4 changes the interaction modifiers from approximately $\mathrm{dist}=0$5 and $\mathrm{dist}=0$6 at low porosity to $\mathrm{dist}=0$7 and $\mathrm{dist}=0$8 at high porosity. The paper argues that the strongest interaction effects are fundamentally hydrodynamic–chemical rather than purely mechanical (Roy et al., 2019).

In two-dimensional amorphous materials, the issue is the energetic stabilization of a cavity boundary. “Rupture of amorphous graphene via void formation” writes the free energy of a void of radius $\mathrm{dist}=0$9 under pressure $t_{\rm in}$0 as

$t_{\rm in}$1

with finite-size-corrected line tension

$t_{\rm in}$2

The critical radius is $t_{\rm in}$3. At $t_{\rm in}$4, the reported values are $t_{\rm in}$5 for flat amorphous graphene and $t_{\rm in}$6 for buckled graphene. The analytical polygon model links the asymptotic lower bound of the line tension to the shear modulus, so small voids are stabilized by boundary cost until the stress-relief term dominates (Jain et al., 2018).

A further astrophysical formulation appears in “Dynamic Voids Surrounded by Shocked Conventional Polytropic Gas Envelopes.” Here the void is a low-mass central cavity in a self-gravitating polytropic sphere, bounded by a pressure-balanced contact discontinuity on or near the zero-mass line $t_{\rm in}$7. With the similarity transformation $t_{\rm in}$8 and the conventional-polytrope condition $t_{\rm in}$9, the model admits global void solutions that cross the sonic critical surface either smoothly or by shocks and propagate into static, inflow, outflow, breeze, or contraction envelopes. In the supernova interpretation, a cavity is first carved by an extremely powerful neutrinosphere and later driven by electromagnetic radiation and/or electron-positron pair plasma; at sufficiently late times the boundary pressure decays and the envelope continues to expand by inertia (Lou et al., 2011).

Taken together, these non-cosmological literatures sharpen the general meaning of void-shoring. The void is maintained, neutralized, enlarged, or eliminated not by its emptiness but by the surrounding medium’s ability to supply contact force, pressure balance, inertia compensation, chemical amplification, or perimeter energy. A plausible synthesis is that void-shoring is the study of how boundaries make voids mechanically and dynamically consequential.