Self-Gravitating Sheet Model

• Self-Gravitating Sheet Model is a reduced-dimensional framework using planar symmetry to simplify gravitational interactions while retaining key features like long-range self-gravity and collective oscillations.
• It unifies distinct formulations—classical Hamiltonian, shearing sheet, and hydrostatic slab models—to analyze dynamics, equilibrium, and non-equilibrium relaxation processes.
• Applications include numerical studies of violent relaxation, core-halo structure formation, and disc dynamics, providing insights into thermalization, ergodicity, and fragmentation phenomena.

The self-gravitating sheet model denotes a family of reduced-dimensional gravitational models built from planar symmetry. In its classical statistical-mechanics usage, it consists of identical infinite parallel sheets moving only along one coordinate and interacting through a pair potential proportional to $|x_i-x_j|$; in disc dynamics, the term often refers to the local two-dimensional self-gravitating shearing sheet; and in hydrostatic work it denotes plane-parallel isothermal or polytropic slabs. Across these variants, the attraction of the model is that it retains long-range self-gravity, collective oscillations, quasi-stationary structure, and tractable field equations while eliminating much of the geometric complexity of full three-dimensional gravity (Joyce et al., 2010, Teles et al., 2011, Paardekooper, 2012).

1. Terminology and model classes

The phrase “self-gravitating sheet model” is not unique in the literature. The classical long-range-interaction model uses infinite parallel mass sheets in three dimensions, constrained so that only their perpendicular displacement matters. A distinct astrophysical usage denotes the local rotating patch of a gravitating disc, namely the self-gravitating shearing sheet. A third usage appears in hydrostatic studies of plane-parallel slabs, especially isothermal and polytropic equilibria. These usages share planar symmetry but differ in dynamics, thermodynamic closure, and intended applications (Souza et al., 2020, Paardekooper, 2012, Ito, 2023).

Usage	Geometry and degrees of freedom	Characteristic equations
Classical Hamiltonian sheet model	$N$ infinite parallel sheets moving along $x$	$H=\sum_i \frac{p_i^2}{2m}+g m^2\sum_{i<j}|x_i-x_j|$
Self-gravitating shearing sheet	Local Cartesian patch $(x,y)$ of a rotating disc	Shearing-sheet hydrodynamics with cooling and thin-disc self-gravity
Hydrostatic slab model	Plane-parallel collisional gas depending on $z$	$$\frac{d^2\Phi}{dz^2}=4\pi G\rho,\ \frac{dp}{dz}=-\rho\,\frac{d\Phi}{dz}$$

In the classical one-dimensional model, the sheets are infinite in the transverse directions, so the interaction becomes effectively one-dimensional. In the shearing-sheet model, the planar symmetry is local rather than global: one follows a corotating patch of a differentially rotating disc. In hydrostatic slab models, planar symmetry is exact, but the system is treated as a collisional fluid rather than a collisionless many-body Hamiltonian system. This suggests that “self-gravitating sheet model” is best understood as a family of reductions rather than a single equation set.

2. Classical one-dimensional Hamiltonian formulation

In the standard equal-mass sheet model, each particle represents an infinite parallel sheet of mass. One common convention writes the pair force as

$$F_{ij}=- gm^{2}\frac{x_i-x_j}{|x_i-x_j|}\equiv -gm^2\,\mathrm{sgn}(x_i-x_j),$$

so the $N$-body Hamiltonian is

$$H=\sum_{i=1}^N \frac{p_i^2}{2m} + g m^2 \sum_{i<j}|x_i-x_j|.$$

A dimensionless convention used for the one-dimensional self-gravitating sheet model (ODSGM) instead writes the potential of a unit point mass at the origin as $N$0, with Hamiltonian

$N$1

Different papers therefore differ by normalization, but not by the underlying ordering-dependent $N$2 interaction (Teles et al., 2011, Joyce et al., 2010).

A distinctive property of the model is that the force on a given sheet depends only on how many sheets lie on each side. If $N$3 and $N$4 are the numbers of sheets to the right and left of $N$5, then

$N$6

or, in the $N$7-normalized convention,

$N$8

Between crossings the accelerations are constant, so the dynamics is piecewise ballistic with uniform acceleration. This makes the model numerically attractive: event-driven algorithms integrate trajectories exactly between crossings, and the longest runs reported in one study conserve total energy to about $N$9 or $x$0, depending on implementation details (Joyce et al., 2010, Joyce et al., 2010).

The continuum limit is the Vlasov or mean-field limit

$x$1

for which the one-particle distribution $x$2 obeys

$x$3

with density

$x$4

and field determined by the one-dimensional Poisson relation

$x$5

or, in the alternative normalization,

$x$6

The collisionless incompressibility of Vlasov flow implies that a one-level waterbag cannot exceed its initial phase-space density at later times (Teles et al., 2011, Joyce et al., 2010).

3. Violent relaxation, Lynden-Bell theory, and core-halo structure

Most of the modern statistical-mechanical work on the classical sheet model concerns the non-equilibrium evolution from waterbag initial data. A symmetric waterbag is written as

$x$7

and, after setting $x$8, its initial energy is

$x$9

The single control parameter left by scaling is the virial number

$H=\sum_i \frac{p_i^2}{2m}+g m^2\sum_{i<j}|x_i-x_j|$0

If $H=\sum_i \frac{p_i^2}{2m}+g m^2\sum_{i<j}|x_i-x_j|$1, macroscopic oscillations are suppressed; if $H=\sum_i \frac{p_i^2}{2m}+g m^2\sum_{i<j}|x_i-x_j|$2, the mismatch between kinetic and potential energies excites collective oscillations (Teles et al., 2011).

Lynden-Bell’s violent-relaxation theory gives, for a one-level waterbag, a Fermi-Dirac-like coarse-grained distribution

$H=\sum_i \frac{p_i^2}{2m}+g m^2\sum_{i<j}|x_i-x_j|$3

obtained by maximizing the mixing entropy

$H=\sum_i \frac{p_i^2}{2m}+g m^2\sum_{i<j}|x_i-x_j|$4

with $H=\sum_i \frac{p_i^2}{2m}+g m^2\sum_{i<j}|x_i-x_j|$5. In the sheet model this prediction depends, up to rescalings, on a single dimensionless parameter

$H=\sum_i \frac{p_i^2}{2m}+g m^2\sum_{i<j}|x_i-x_j|$6

where $H=\sum_i \frac{p_i^2}{2m}+g m^2\sum_{i<j}|x_i-x_j|$7 is the minimum possible energy for given mass and phase-space density. Numerical tests show that Lynden-Bell theory is reasonably good in the low-energy region and for sufficiently gentle relaxation, but it is generally not even qualitatively correct at higher energies. At $H=\sum_i \frac{p_i^2}{2m}+g m^2\sum_{i<j}|x_i-x_j|$8 and $H=\sum_i \frac{p_i^2}{2m}+g m^2\sum_{i<j}|x_i-x_j|$9, different one-level waterbags cluster near the Lynden-Bell prediction in the $(x,y)$0 plane; at $(x,y)$1, different initial shapes produce quasi-stationary states that scatter widely, even with different signs of the order parameters (Joyce et al., 2010).

A complementary theory explains this failure in terms of resonance-driven core-halo formation. The rms envelope

$(x,y)$2

obeys, for the initial waterbag and under an affine short-time approximation,

$(x,y)$3

The oscillating envelope then defines a time-dependent mean field for test particles,

$(x,y)$4

and this dynamics exhibits resonance islands. Resonant particles absorb energy from the bulk oscillation, move to large amplitudes, and form a dilute halo; the remaining low-energy region contracts into a maximally occupied core. The resulting parameter-free stationary ansatz is

$(x,y)$5

with

$(x,y)$6

The remaining parameters $(x,y)$7 and $(x,y)$8 are fixed self-consistently by mass and energy conservation. The corresponding Poisson closure is piecewise: $(x,y)$9 Direct $z$0-body tests with $z$1 sheets run to $z$2 show very good agreement for the final density and velocity distributions, including examples with $z$3 and $z$4, and the paper emphasizes that this agreement is achieved without fitting parameters (Teles et al., 2011).

4. Thermalization, ergodicity, and exact steady states

The collisionless quasi-stationary state is not the end of the story at finite $z$5. For the equal-mass sheet model, the exact thermal equilibrium in the mean-field limit is known: $z$6 with

$z$7

This equilibrium is separable in $z$8 and $z$9, and that fact motivates the use of mixed-moment order parameters

$$\frac{d^2\Phi}{dz^2}=4\pi G\rho,\ \frac{dp}{dz}=-\rho\,\frac{d\Phi}{dz}$$0

They vanish in thermal equilibrium but remain nonzero in generic quasi-stationary states. For rectangular waterbags, the relaxation time to equilibrium scales approximately linearly in $$\frac{d^2\Phi}{dz^2}=4\pi G\rho,\ \frac{dp}{dz}=-\rho\,\frac{d\Phi}{dz}$$1, depends strongly on the initial virial ratio $$\frac{d^2\Phi}{dz^2}=4\pi G\rho,\ \frac{dp}{dz}=-\rho\,\frac{d\Phi}{dz}$$2, and is shorter for colder initial conditions. Reported values include $$\frac{d^2\Phi}{dz^2}=4\pi G\rho,\ \frac{dp}{dz}=-\rho\,\frac{d\Phi}{dz}$$3 for $$\frac{d^2\Phi}{dz^2}=4\pi G\rho,\ \frac{dp}{dz}=-\rho\,\frac{d\Phi}{dz}$$4, $$\frac{d^2\Phi}{dz^2}=4\pi G\rho,\ \frac{dp}{dz}=-\rho\,\frac{d\Phi}{dz}$$5 for $$\frac{d^2\Phi}{dz^2}=4\pi G\rho,\ \frac{dp}{dz}=-\rho\,\frac{d\Phi}{dz}$$6, and $$\frac{d^2\Phi}{dz^2}=4\pi G\rho,\ \frac{dp}{dz}=-\rho\,\frac{d\Phi}{dz}$$7 for $$\frac{d^2\Phi}{dz^2}=4\pi G\rho,\ \frac{dp}{dz}=-\rho\,\frac{d\Phi}{dz}$$8. In the best-sampled case, the decay of the order parameter is better fitted by a stretched exponential,

$$\frac{d^2\Phi}{dz^2}=4\pi G\rho,\ \frac{dp}{dz}=-\rho\,\frac{d\Phi}{dz}$$9

than by a hyperbolic tangent form (Joyce et al., 2010).

Ergodicity is likewise state-dependent. One study defines the time-averaged single-particle observables

$$F_{ij}=- gm^{2}\frac{x_i-x_j}{|x_i-x_j|}\equiv -gm^2\,\mathrm{sgn}(x_i-x_j),$$0

and corresponding dispersions

$$F_{ij}=- gm^{2}\frac{x_i-x_j}{|x_i-x_j|}\equiv -gm^2\,\mathrm{sgn}(x_i-x_j),$$1

For non-homogeneous states, these quantities decay on a timescale comparable to the relaxation time to equilibrium; a representative open-boundary simulation with $$F_{ij}=- gm^{2}\frac{x_i-x_j}{|x_i-x_j|}\equiv -gm^2\,\mathrm{sgn}(x_i-x_j),$$2, $$F_{ij}=- gm^{2}\frac{x_i-x_j}{|x_i-x_j|}\equiv -gm^2\,\mathrm{sgn}(x_i-x_j),$$3, and $$F_{ij}=- gm^{2}\frac{x_i-x_j}{|x_i-x_j|}\equiv -gm^2\,\mathrm{sgn}(x_i-x_j),$$4 gives $$F_{ij}=- gm^{2}\frac{x_i-x_j}{|x_i-x_j|}\equiv -gm^2\,\mathrm{sgn}(x_i-x_j),$$5. For homogeneous states with periodic boundary conditions treated through an Ewald construction and then the limit

$$F_{ij}=- gm^{2}\frac{x_i-x_j}{|x_i-x_j|}\equiv -gm^2\,\mathrm{sgn}(x_i-x_j),$$6

the gravitational contribution to the force vanishes and the one-particle distribution has zero collision term. The paper therefore concludes that homogeneous states are effectively non-ergodic, whereas non-homogeneous states are ergodic only on a time window of the order of the relaxation time to equilibrium (Souza et al., 2020).

The model also admits a special exact Vlasov equilibrium. In a monoenergetic steady state with

$$F_{ij}=- gm^{2}\frac{x_i-x_j}{|x_i-x_j|}\equiv -gm^2\,\mathrm{sgn}(x_i-x_j),$$7

and units

$$F_{ij}=- gm^{2}\frac{x_i-x_j}{|x_i-x_j|}\equiv -gm^2\,\mathrm{sgn}(x_i-x_j),$$8

the self-consistent half-orbit is

$$F_{ij}=- gm^{2}\frac{x_i-x_j}{|x_i-x_j|}\equiv -gm^2\,\mathrm{sgn}(x_i-x_j),$$9

The turning points are at $N$0, the maximal excursion is

$N$1

and the motion is periodic with period

$N$2

The density is

$N$3

so near the turning point it has the fold-caustic divergence

$N$4

This is not a generic equilibrium of the sheet model, but an exactly solvable cold two-stream stationary solution in the Vlasov limit, useful as a benchmark and pedagogical example (Nityananda, 4 Jan 2026).

5. The self-gravitating shearing sheet in disc dynamics

In accretion-disc theory, the self-gravitating shearing sheet is a local rotating patch centered on

$N$5

with local coordinates

$N$6

The fluid is treated as a razor-thin two-dimensional hydrodynamic sheet with surface density $N$7, pressure $N$8, internal energy density $N$9, and gravitational potential $$H=\sum_{i=1}^N \frac{p_i^2}{2m} + g m^2 \sum_{i<j}|x_i-x_j|.$$0. A standard formulation is

$$H=\sum_{i=1}^N \frac{p_i^2}{2m} + g m^2 \sum_{i<j}|x_i-x_j|.$$1

$$H=\sum_{i=1}^N \frac{p_i^2}{2m} + g m^2 \sum_{i<j}|x_i-x_j|.$$2

$$H=\sum_{i=1}^N \frac{p_i^2}{2m} + g m^2 \sum_{i<j}|x_i-x_j|.$$3

The Toomre parameter is

$$H=\sum_{i=1}^N \frac{p_i^2}{2m} + g m^2 \sum_{i<j}|x_i-x_j|.$$4

with $$H=\sum_{i=1}^N \frac{p_i^2}{2m} + g m^2 \sum_{i<j}|x_i-x_j|.$$5 for Keplerian rotation, and local thermal balance implies

$$H=\sum_{i=1}^N \frac{p_i^2}{2m} + g m^2 \sum_{i<j}|x_i-x_j|.$$6

For a periodic Fourier mode, the thin-sheet potential is

$$H=\sum_{i=1}^N \frac{p_i^2}{2m} + g m^2 \sum_{i<j}|x_i-x_j|.$$7

This local model is physically distinct from the classical one-dimensional Hamiltonian sheet model, but it inherits the same emphasis on reduced geometry and self-consistent long-range gravity (Paardekooper, 2012).

A major result in this setting is that fragmentation in the two-dimensional local $$H=\sum_{i=1}^N \frac{p_i^2}{2m} + g m^2 \sum_{i<j}|x_i-x_j|.$$8-cooling model is stochastic rather than sharply thresholded. At $$H=\sum_{i=1}^N \frac{p_i^2}{2m} + g m^2 \sum_{i<j}|x_i-x_j|.$$9 resolution and short runtimes $N$00, fragmentation is recovered for $N$01, reproducing the traditional result. At the same resolution but with long runtimes $N$02, fragmentation is observed up to $N$03. At $N$04, $N$05 becomes stochastic, with two of four runs fragmenting at $N$06 and $N$07. Transient clumps that can survive shear are seen for all $N$08 considered, and the highest $N$09 for which actual fragmentation is observed is $N$10, requiring $N$11 resolution. The paper therefore argues that in the simple two-dimensional local $N$12-cooling sheet there is no sharp fragmentation boundary; gravito-turbulence is only metastable, continuously producing transient self-gravitating clumps whose survival is a rare-event problem (Paardekooper, 2012).

The same local sheet framework has been extended to dust-gas systems. In one set of simulations, passive drag-coupled solids with stopping times

$N$13

are evolved in a gravito-turbulent gas sheet. Spiral density waves concentrate particles most strongly for $N$14, producing local surface-density enhancements of order $N$15 to $N$16 over the mean and a narrow velocity dispersion, implying more ordered motion and lower relative collision speeds than absolute particle speeds alone would suggest (Gibbons et al., 2012).

When particle self-gravity and back-reaction are added, overdense particle filaments formed inside gas spiral crests contract into highly dense, gravitationally bound particle clumps. The total mass contained in bound structures appears nearly independent of the cooling time, and the paper concludes that dust trapping by self-gravitating density waves can trigger collapse over a larger range of radii than trapping efficiency alone would suggest (Gibbons et al., 2014). A further extension shows that transient anticyclonic vortices in the same local self-gravitating sheet can trap small and intermediate particles even more efficiently than relatively large-scale density waves. Before particle self-gravity is switched on, concentrations rise to nearly $N$17 the mean; after self-gravity is enabled, many bound clumps form, with typical masses

$N$18

and the majority of clumps form within vortices rather than wave crests (Gibbons et al., 2015).

6. Hydrostatic slabs, exact sheet equilibria, and gravitational polarization

A different branch of the literature studies self-gravitating sheets as collisional hydrostatic slabs. For an isothermal sheet, one writes

$N$19

With dimensionless variables

$N$20

the plane-parallel isothermal Lane-Emden equation becomes

$N$21

with exact solution

$N$22

If a test mass sheet is inserted at the midplane and the gas is assumed to relax to a new thermal equilibrium at the same temperature, the perturbed potential satisfies

$N$23

Writing

$N$24

the regularized response equation is

$N$25

The amplification factor is defined as

$N$26

and in the isothermal sheet the maximal amplification is $N$27, with

$N$28

The interpretation is gravitational amplification rather than screening: the medium strengthens the bare field of the inserted sheet (Ito, 2023).

For polytropic slabs with

$N$29

the unperturbed planar Lane-Emden equation is

$N$30

and for $N$31 the sheet has finite height $N$32. Adding a central test mass sheet gives

$N$33

which linearizes to

$N$34

Defining

$N$35

one obtains the regularized equation

$N$36

The response is summarized by

$N$37

For $N$38, there is no amplification. For $N$39,

$N$40

$N$41

$N$42

so

$N$43

For general $N$44, the paper reports that $N$45 increases with $N$46, asymptotically approaching about $N$47, $N$48 is maximized near $N$49, and the relative shortening obeys

$N$50

at large $N$51. These studies are explicitly static and collisional; they do not address sheet crossings, violent relaxation, or finite-$N$52 kinetic evolution (Ito, 2023).

The hydrostatic and Hamiltonian branches therefore meet mainly at the level of geometry. The former emphasizes equilibrium structure and collective response of collisional slabs; the latter emphasizes Vlasov dynamics, quasi-stationary states, ergodicity, and very slow relaxation. Taken together, they show that the self-gravitating sheet model is not a single theory but a compact framework in which planar symmetry supports exact solutions, controlled numerical tests, and sharply formulated questions about long-range gravity.