Calibrating the Brody exponent as a quantitative measure of short-range exclusion in 2D spatial point processes

This paper establishes a calibrated measurement framework for quantifying short-range exclusion in 2D spatial point processes by correcting the Brody exponent's baseline from 0.96 (instead of the inappropriate 1D Poisson value) and validating a strong empirical correlation between the exponent and the effective hard-core radius across diverse datasets including manufactured surfaces, prime number embeddings, and fractal structures.

The Gist

Imagine you are looking at a crowded dance floor. You want to know: How much personal space does each dancer demand?

In the world of physics and mathematics, scientists study "point processes"—which is just a fancy way of describing dots scattered on a surface. These dots could be peaks on a rough metal surface, stars in a galaxy, or even the number 2, 3, 5, 7 (prime numbers) placed on a grid.

For a long time, scientists had a tool called the Brody distribution to measure how much these dots "hated" being close to each other. However, they were using an old, broken ruler. They were measuring a 2D dance floor using a rulebook designed for a 1D line.

This paper is about fixing the ruler and teaching us how to measure "personal space" (exclusion) correctly in two dimensions.

Here is the breakdown of what the paper does, using simple analogies:

1. The Broken Ruler (The Calibration Problem)

Imagine you are trying to measure how crowded a room is.

• The Old Way: The scientists previously assumed that if people were standing randomly (no one pushing anyone), the "crowdedness score" would be 0.
• The New Discovery: The author, Dawid Kucharski, realized that in a 2D room, even if people are standing completely randomly, the geometry of the room makes them look slightly more spaced out than on a 1D line.
• The Fix: The paper proves that the "random baseline" isn't 0. It's actually 0.96.
  • Analogy: Think of it like a thermometer. If the old thermometer said "0 degrees" for a freezing day, but it was actually 32 degrees, every reading was wrong. The author recalibrated the thermometer so that "Random" now reads 0.96. Anything above 0.96 means the dots are actively pushing each other away.

2. The "Personal Space" Score (The Brody Exponent $\beta$)

The paper introduces a single number, called $\beta$ (beta), to describe how much the dots avoid each other.

• $\beta \approx 0.96$: The dots are like strangers at a party who don't care who they stand next to. They are randomly scattered.
• $\beta \approx 2.0$: The dots are like polite guests who give each other a comfortable arm's length of space.
• $\beta > 5$: The dots are like soldiers in a perfect marching formation. They are extremely rigid and refuse to be close to anyone.

The paper creates a translation guide (a calibration curve) that turns this abstract math number ($\beta$) into a physical measurement: How big is their "no-entry zone"?

• If $\beta$ is high, the "no-entry zone" (the radius where a dot won't let another dot enter) is large.
• If $\beta$ is low, the "no-entry zone" is tiny or non-existent.

3. Testing the Ruler on Three Different "Dance Floors"

To prove the new ruler works, the author tested it on three very different types of "dots":

A. Rough Metal Surfaces (Manufacturing)

Imagine looking at a piece of metal under a microscope. It looks like a mountain range with peaks and valleys.

• The Test: The author measured the distance between the highest peaks.
• The Result: Different manufacturing processes created different "personal space" scores.
  • Burnished (polished) metal: The peaks were very organized and kept far apart (High $\beta$).
  • Ground or turned metal: The peaks were messy and closer together (Low $\beta$, close to random).
• Takeaway: The tool successfully told the difference between a smooth, organized surface and a rough, chaotic one.

B. Prime Numbers (Math)

This is the most surprising part. The author took prime numbers (2, 3, 5, 7, 11...) and mapped them onto a 2D grid.

• The Mystery: Do prime numbers have a hidden "personal space" rule?
• The Twist: The answer depends entirely on how you draw the map.
  • Map A (Row-by-Row): If you write numbers left-to-right, top-to-bottom, the primes show a strong "personal space" rule ($\beta = 2.15$). They seem to avoid each other.
  • Map B (Cantor Spiral): If you write them in a diagonal zig-zag pattern, the "personal space" disappears ($\beta = 1.40$). They look random.
• The Lesson: The "exclusion" isn't a magic property of the numbers themselves; it's created by the geometry of the map. If you change the map, you change the result. This proves the tool is sensitive enough to detect that the arrangement matters, not just the numbers.

C. Optical Measurements (Interferometry)

The author also tested a certified roundness standard (a perfectly round metal ring) using light waves.

• The Result: The peaks of the light waves showed a "personal space" score of 2.00. This confirmed the tool works on optical data, not just metal or math.

4. The "Density" Trap

A major discovery in the paper is that density matters.

• The Analogy: Imagine a crowded subway car vs. an empty park.
  • In the subway (high density), people are forced to be close. Even if they hate being close, they can't help it.
  • In the park (low density), people can spread out.
• The Finding: The author found that if you take the same group of dots and randomly remove half of them (thinning), the $\beta$ score changes.
• The Rule: You cannot compare a crowded surface to a sparse one directly. You must compare "apples to apples" (same density) or use the new "radius" translation guide to normalize the results.

Summary: What Did We Learn?

1. We fixed the baseline: Random 2D dots aren't "0"; they are 0.96.
2. We have a translator: We can now turn the abstract math number $\beta$ into a physical "exclusion radius" (how much space a dot demands).
3. Context is King:
   • For metal, it tells us how smooth or rough the process was.
   • For math, it tells us that the way we arrange numbers changes their statistical "personality."
   • For optics, it confirms the method works on light data.

In short: The paper gives scientists a new, calibrated ruler to measure how much "personal space" dots demand in a 2D world, correcting a long-standing mistake and showing that how you arrange those dots is just as important as what the dots are.

Technical Summary

Technical Summary: Calibrating the Brody Exponent for 2D Spatial Point Processes

Problem Statement
The statistical characterization of spatial point patterns often lacks a single, interpretable parameter to quantify short-range exclusion strength on a continuous scale. While the Brody distribution, originally developed as a phenomenological interpolation between Poisson and Wigner statistics in quantum chaos, has been widely applied to 1D spectra, its extension to two-dimensional (2D) spatial point processes faces two critical unresolved issues:

1. Incorrect Baseline: The standard 1D Poisson reference ($\beta = 0$) is inappropriate for 2D. Under Complete Spatial Randomness (CSR) in 2D, the nearest-neighbour spacing distribution follows a Rayleigh distribution, which is not a member of the Brody family. Consequently, the appropriate 2D CSR baseline had not been established.
2. Lack of Physical Calibration: There was no empirical framework to translate the dimensionless Brody parameter ($\beta$) into a physically interpretable measure of exclusion, such as an effective hard-core radius.

Methodology
The study develops a calibrated measurement framework comprising three methodological pillars:

1. Brody $\beta$ Estimation: The Brody distribution $P_\beta(s) = (\beta + 1)as^\beta \exp(-as^{\beta+1})$ is fitted to 2D nearest-neighbour spacing (NNS) distributions using Maximum Likelihood Estimation (MLE). The parameter $\beta$ is constrained to $\beta \geq 0$, where $\beta=0$ represents the Poisson limit (uncorrelated) and higher values indicate increasing short-range order.
2. Synthetic Calibration: A monotonic mapping between $\beta$ and the normalized exclusion radius ($r_{excl}/\langle NN \rangle$) is established using synthetic point processes. This includes sequential rejection sampling for hard-core processes and Strauss soft-core processes. The 2D Ginibre ensemble (determinantal repulsion) and 2D Poisson processes serve as reference benchmarks.
3. Control Protocols and Validation: To isolate genuine exclusion from confounding factors, the framework employs:
   • Density-thinning experiments: To decouple exclusion strength from point density.
   • Embedding-robustness tests: Using deterministic 2D mappings (Row-major, Ulam spiral, Cantor pairing) on arithmetic sequences to test if exclusion signals are intrinsic or geometry-dependent.
   • Null-model ensembles: Including Gaussian, Poisson, 1/f noise, Brownian, and matched-density Bernoulli processes.
   • Complementary Descriptors: For binary patterns where peak detection fails, the correlation dimension ($D_2$) and Banach-space invariants are used.

Key Results

• 2D CSR Baseline Correction: The study recalibrates the 2D CSR baseline to $\beta = 0.96 \pm 0.15$, correcting the inappropriate use of $\beta = 0$. This baseline is density-stable across the tested range ($\rho \in [0.01, 0.10]$). For binary fields at low fill fractions, a distinct baseline of $\beta_{CSR}^{binary} \approx 1.46$ is identified.
• Empirical $\beta$–$r_{excl}$ Calibration: A strong empirical correlation (Spearman $\rho = 0.988$) is validated between $\beta$ and the effective hard-core radius. The mapping is monotonic: $\beta \approx 0.96$ corresponds to no exclusion, $\beta \approx 2$ aligns with Ginibre-level repulsion ($r_{excl} \approx 0.55 \langle NN \rangle$), and $\beta \gtrsim 5$ indicates near-crystalline exclusion.
• Manufactured Surfaces: Analysis of 58 manufactured surfaces (10 materials, 10 processes) reveals a wide range of $\beta$ values ($0.00$ to $4.61$). Burnished and honed surfaces exhibit higher exclusion ($\beta \approx 2.2–4.3$), while turned and ground surfaces often approach the CSR baseline.
• Prime Number Embeddings:
  • Row-major embedding: Yields $\beta = 2.15$, indicating strong exclusion.
  • Cantor pairing: Yields $\beta = 1.40$, which is statistically indistinguishable from the Bernoulli null ($\beta \approx 1.47$).
  • Sparse-integer control: Random integers embedded in the same row-major geometry yield $\beta \approx 1.47$, confirming that the high $\beta$ for primes is a genuine arithmetic signal (due to gap constraints) rather than a sparsity artifact.
  • Conclusion: 2D exclusion in arithmetic embeddings is not intrinsic to the sequence but is created or destroyed by the embedding geometry.
• Density Dependence: Absolute $\beta$ values are density-dependent (e.g., prime $\beta$ drops from 2.15 at $\rho=0.096$ to 1.46 at $\rho=0.032$). However, the exclusion strength relative to the density-matched CSR baseline remains consistent. Cross-domain comparisons require density matching or the use of the density-invariant $r_{excl}/\langle NN \rangle$ calibration.
• Interferometric Profilometry: A proof-of-concept application on a certified roundness standard yields $\beta = 2.00$, consistent with the exclusion observed in other domains.

Significance and Claims
The paper claims to provide a calibrated measurement framework for the reproducible characterization of short-range exclusion in 2D spatial point processes. Its primary contributions are:

1. Correction of the Baseline: Establishing that the 2D CSR baseline is $\beta \approx 0.96$, not 0, which is essential for avoiding spurious exclusion claims.
2. Physical Interpretation: Providing an empirical $\beta$–$r_{excl}$ mapping that translates the abstract Brody parameter into a physically interpretable exclusion radius, analogous to the hard-sphere packing fraction in liquid-state theory.
3. Methodological Rigor: Demonstrating that $\beta$ captures exclusion strength rather than point density or embedding artifacts, provided that density matching and appropriate baselines are applied.
4. Phenomenological Scope: The authors explicitly state that $\beta$ is used as an empirical descriptor of spatial order, not as a claim of underlying random-matrix symmetry classes (e.g., $\beta > 2$ does not imply higher-order symmetry classes but simply stronger exclusion).

The framework is demonstrated to be applicable across physical (manufactured surfaces), mathematical (prime numbers), and optical (interferometry) domains, provided the necessary corrections for density, embedding geometry, and baseline selection are applied. The study does not claim to reveal a universal exclusion mechanism but rather to supply the tools required to quantify it quantitatively and reproducibly.