pith. sign in

arxiv: 2201.08478 · v3 · submitted 2022-01-20 · ❄️ cond-mat.stat-mech · astro-ph.GA

Statistical properties of the gravitational force through ordering statistics

Constantin Payerne , Vincent Rossetto This is my paper

Pith reviewed 2026-05-24 12:40 UTC · model grok-4.3

classification ❄️ cond-mat.stat-mech astro-ph.GA
keywords gravitational forceorder statisticsHoltsmark distributionnearest neighborspoint massesPoisson processvariancerandom gas
0
0 comments X

The pith

The divergent variance of the Holtsmark gravitational force distribution arises entirely from the nearest neighbor.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper applies order statistics to distances between a test particle and successive nearest neighbors among randomly placed point masses. It derives the probability density functions for these distances in arbitrary dimensions. In three dimensions the analysis shows that the infinite variance of the total force comes only from the first neighbor, while every more distant neighbor contributes a finite variance. This separates local from distant contributions to the force statistics. A reader would care because the result identifies which particles control the singular behavior of the classic Holtsmark distribution.

Core claim

Using order statistics, the probability density functions of distances to the n-th nearest neighbors are derived in arbitrary spatial dimensions. In three dimensions, the formally divergent variance of the Holtsmark probability density function arises entirely from the first nearest neighbor, while contributions from more distant neighbors remain finite.

What carries the argument

Order statistics applied to the sequence of nearest-neighbor distances in a homogeneous Poisson point process of point masses.

If this is right

  • The total gravitational force variance decomposes into independent contributions ranked by neighbor order.
  • In three dimensions the variance contributed by the second and all further neighbors stays finite.
  • The same ordering approach yields finite moments for distant neighbors in any dimension.
  • Nearest-neighbor dominance explains the origin of the Holtsmark variance divergence.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Numerical gravity simulations may need only local high-resolution sampling to reproduce the full force variance.
  • Introducing spatial correlations among the masses would likely shift which neighbor order controls the divergence.
  • The same ordering decomposition could be applied to other inverse-square or long-range potentials.

Load-bearing premise

The point masses form an infinite, homogeneous, and uncorrelated random gas.

What would settle it

A Monte Carlo realization of random point masses in a large volume; the partial force variance computed from all particles except the nearest one should remain bounded as the volume tends to infinity.

Figures

Figures reproduced from arXiv: 2201.08478 by Constantin Payerne, Vincent Rossetto.

Figure 1
Figure 1. Figure 1: FIG. 1. Left - Relative error [PITH_FULL_IMAGE:figures/full_fig_p008_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Left - Probability density function [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗
read the original abstract

We investigate the statistical distribution of Newtonian gravitational forces acting on a test particle embedded in an infinite, homogeneous, and uncorrelated random gas of point masses. Using order statistics, we derive the probability density functions of distances to the $n$-th nearest neighbors in arbitrary spatial dimensions and analyze their contributions to the total gravitational force. We show that, in three dimensions, the formally divergent variance of the Holtsmark probability density function arises entirely from the first nearest neighbor, while contributions from more distant neighbors remain finite. Our results provide a clear decomposition of local versus distant contributions to the gravitational forces exerted on a test particle, and explore the dominant role of nearest neighbors in the divergence of the Holtsmark distribution's variance.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

0 major / 2 minor

Summary. The paper applies order statistics to the distances of point masses drawn from a homogeneous Poisson point process in d dimensions. It derives the PDFs of the n-th nearest-neighbor distances and computes their separate contributions to the Newtonian gravitational force (and its variance) on a test particle. The central claim is that, in three dimensions, the formally divergent variance of the Holtsmark distribution is produced entirely by the nearest neighbor (n=1), while the contributions from all n≥2 remain finite.

Significance. If the derivation holds, the work supplies a parameter-free decomposition that isolates the source of the Holtsmark variance divergence to the nearest-neighbor term alone. This clarifies the local versus distant contributions to gravitational force fluctuations under the standard Poisson assumption and is consistent with the small-r scaling of the ordered-distance PDFs (r_n ~ r^{3n-1} near zero, so that E[r_n^{-4}] diverges only for n=1). The result strengthens the analytic understanding of force statistics without additional fitted parameters or self-citations.

minor comments (2)
  1. [§2] §2 (or wherever the ordered-distance PDFs are stated): the normalization constants for the n-th neighbor PDF in general d could be written explicitly to facilitate direct comparison with the force integrals that follow.
  2. [Figure 1] Figure 1 (or equivalent): the plotted force-component distributions would benefit from an inset or caption note confirming that the n=1 curve alone produces the logarithmic divergence while n≥2 curves converge.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for their careful reading, positive assessment of the work, and recommendation to accept the manuscript. We are pleased that the central result on the nearest-neighbor origin of the Holtsmark variance divergence is viewed as a useful clarification.

Circularity Check

0 steps flagged

No significant circularity; derivation self-contained

full rationale

The paper starts from the standard Poisson point process in infinite homogeneous space and applies order statistics to derive the pdf of the n-th nearest-neighbor distance r_n ~ r^{3n-1} near zero. The variance integral E[r_n^{-4}] then diverges only for n=1 by direct power-counting (exponent 3n-5 = -2), while remaining finite for n≥2; vector cross terms are likewise finite. This reduction uses only the known void probability of the PPP and does not invoke fitted parameters, self-citations, or any ansatz that presupposes the target result. The central claim is therefore an independent consequence of the input measure.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the Poisson point process model and the application of order statistics to compute force contributions; no free parameters or new entities are introduced in the abstract.

axioms (1)
  • domain assumption The distribution of point masses is a homogeneous Poisson point process in infinite space.
    This is the standard setup for deriving the Holtsmark distribution and allows closed-form expressions for neighbor distance PDFs.

pith-pipeline@v0.9.0 · 5643 in / 1220 out tokens · 40763 ms · 2026-05-24T12:40:29.683888+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

9 extracted references · 9 canonical work pages

  1. [1]

    , " * write output.state after.block = add.period write newline

    ENTRY address archivePrefix author booktitle chapter doi edition editor eprint howpublished institution journal key month note number organization pages publisher school series title type volume year label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts #0 'before.all :=...

    work page

  2. [2]

    write newline

    " write newline "" before.all 'output.state := FUNCTION format.doi doi empty "" "doi:" doi * if FUNCTION format.eprint eprint empty "" archivePrefix empty "" archivePrefix ":" * if eprint field.or.null * if FUNCTION format.pid eprint empty format.doi format.eprint if FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = ...

    work page

  3. [3]

    "":K.Xn;wDxx8|M K. <x0222h ,˖ m СCtR<S 6m 9

    11em plus .33em minus .07em @technote 4000 4000 100 4000 4000 500 `\.=1000 = #1 #1 #1 0pt [0pt][0pt] #1 * \| ** #1 \@IEEEauthorblockNstyle \@IEEEauthorblockAstyle \@IEEEauthordefaulttextstyle \@IEEEauthorblockconfadjspace -0.25em \@IEEEauthorblockNtopspace 0.0ex \@IEEEauthorblockAtopspace 0.0ex \@IEEEauthorblockNinterlinespace 2.6ex \@IEEEauthorblockAinte...

    work page

  4. [4]

    1943, Rev

    Chandrasekhar, S. 1943, Rev. Mod. Phys., 15, 1

    work page 1943

  5. [5]

    2009, The European Physical Journal B volume, 70, 413–433

    Chavanis, P. 2009, The European Physical Journal B volume, 70, 413–433

    work page 2009

  6. [6]

    1992, The Journal of Chemical Physics, 97, 9276

    Mazur, S. 1992, The Journal of Chemical Physics, 97, 9276

    work page 1992

  7. [7]

    2002, Journal of Physics: Condensed Matter, 14, 2141–2152

    Pietronero, L., Bottaccio, M., Mohayaee, R., & Montuori, M. 2002, Journal of Physics: Condensed Matter, 14, 2141–2152

    work page 2002

  8. [8]

    2008, The European Physical Journal B, 64, 615–623

    Sylos Labini, F., & Pietronero, L. 2008, The European Physical Journal B, 64, 615–623

    work page 2008

  9. [9]

    write newline

    " write newline "" before.all 'output.state := FUNCTION fin.entry write newline FUNCTION new.block output.state before.all = 'skip after.block 'output.state := if FUNCTION new.sentence output.state after.block = 'skip output.state before.all = 'skip after.sentence 'output.state := if if FUNCTION not #0 #1 if FUNCTION and 'skip pop #0 if FUNCTION or pop #1...

    work page