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

Technische Samenvatting: Kalibratie van de Brody-exponent voor 2D ruimtelijke puntprocessen

Probleemstelling
De statistische karakterisering van ruimtelijke puntenpatronen mist vaak een enkele, interpreteerbare parameter om kortstondige uitsluitingssterkte op een continue schaal te kwantificeren. De Brody-distributie, oorspronkelijk ontwikkeld als een fenomenologische interpolatie tussen Poisson- en Wigner-statistieken in de kwantumchaos, is breed toegepast op 1D-spectra, maar de uitbreiding naar twee-dimensionale (2D) ruimtelijke puntprocessen staat voor twee kritieke onopgeloste problemen:

1. Incorrecte Baseline: De standaard 1D Poisson-referentie ($\beta = 0$) is ongepast voor 2D. Onder Complete Spatiale Willekeur (CSR) in 2D volgt de naburige afstand (nearest-neighbour spacing) een Rayleigh-distributie, die geen lid is van de Brody-familie. Bijgevolg was de geschikte 2D CSR-baseline niet eerder vastgesteld.
2. Gebrek aan Fysieke Kalibratie: Er bestond geen empirisch kader om de dimensieloze Brody-parameter ($\beta$) te vertalen naar een fysiek interpreteerbare maat voor uitsluiting, zoals een effectieve harde kernstraal (hard-core radius).

Methodologie
De studie ontwikkelt een gekalibreerd meetkader bestaande uit drie methodologische pijlers:

1. Brody $\beta$-schatting: De Brody-distributie $P_\beta(s) = (\beta + 1)as^\beta \exp(-as^{\beta+1})$ wordt gefit aan 2D naburige afstand-verdelingen (NNS) met behulp van Maximum Likelihood Estimation (MLE). De parameter $\beta$ is begrensd op $\beta \geq 0$, waarbij $\beta=0$ de Poisson-limiet (ongecorreleerd) vertegenwoordigt en hogere waarden duiden op toenemende kortstondige orde.
2. Synthetische Kalibratie: Een monotone mapping tussen $\beta$ en de genormaliseerde uitsluitingsstraal ($r_{excl}/\langle NN \rangle$) wordt vastgesteld met behulp van synthetische puntprocessen. Dit omvat sequentiële reject-sampling voor hard-core processen en Strauss soft-core processen. De 2D Ginibre-ensemble (determinantiale repulsie) en 2D Poisson-processen dienen als referentiemodellen.
3. Controleprotocollen en Validatie: Om echte uitsluiting te isoleren van verwarrende factoren, gebruikt het kader:
   • Dichtheid-verdunningsexperimenten: Om de uitsluitingssterkte te ontkoppelen van de puntendichtheid.
   • Inbeddings-robuustheidstests: Gebruikmakend van deterministische 2D-mappings (Row-major, Ulam-spiraal, Cantor-paring) op rekenkundige reeksen om te testen of uitsluitingssignalen intrinsiek zijn of geometrie-afhankelijk.
   • Null-model ensembles: Inclusief Gaussische, Poisson, 1/f ruis, Brownse en gematchte-densiteit Bernoulli-processen.
   • Complementaire Descriptoren: Voor binaire patronen waar detectiepieken falen, worden de correlatiedimensie ($D_2$) en Banach-ruimte invarianten gebruikt.

Belangrijkste Resultaten

• Correctie van de 2D CSR-baseline: De studie herkalibreert de 2D CSR-baseline naar $\beta = 0.96 \pm 0.15$, waarmee het ongepaste gebruik van $\beta = 0$ wordt gecorrigeerd. Deze baseline is dichtheidsstabiel over het geteste bereik ($\rho \in [0.01, 0.10]$). Voor binaire velden bij lage vulfracties is een afzonderlijke baseline van $\beta_{CSR}^{binary} \approx 1.46$ geïdentificeerd.
• Empirische $\beta$–$r_{excl}$ Kalibratie: Een sterke empirische correlatie (Spearman $\rho = 0.988$) wordt gevalideerd tussen $\beta$ en de effectieve harde kernstraal. De mapping is monotoon: $\beta \approx 0.96$ komt overeen met geen uitsluiting, $\beta \approx 2$ komt overeen met Ginibre-niveau repulsie ($r_{excl} \approx 0.55 \langle NN \rangle$), en $\beta \gtrsim 5$ duidt op bijna kristallijne uitsluiting.
• Gemanufactureerde Oppervlakken: Analyse van 58 gemanufactureerde oppervlakken (10 materialen, 10 processen) onthult een breed scala aan $\beta$-waarden ($0.00$ tot $4.61$). Gepolijste (burnished) en geslepen (honed) oppervlakken vertonen hogere uitsluiting ($\beta \approx 2.2–4.3$), terwijl gedraaide (turned) en geslepen (ground) oppervlakken vaak de CSR-baseline benaderen.
• Priemgetal Inbeddingen:
  • Row-major inbedding: Levert $\beta = 2.15$ op, wat wijst op sterke uitsluiting.
  • Cantor-paring: Levert $\beta = 1.40$ op, wat statistisch ononderscheidbaar is van de Bernoulli-null ($\beta \approx 1.47$).
  • Sparse-integer controle: Willekeurige gehele getallen ingebed in dezelfde row-major geometrie leveren $\beta \approx 1.47$ op, wat bevestigt dat de hoge $\beta$ voor priemgetallen een echt rekenkundig signaal is (door gap-restricties) in plaats van een artefact van de schaarste.
  • Conclusie: 2D uitsluiting in rekenkundige inbeddingen is niet intrinsiek aan de reeks, maar wordt gecreëerd of vernietigd door de inbeddingsgeometrie.
• Dichtheidsafhankelijkheid: Absolute $\beta$-waarden zijn dichtheidsafhankelijk (bijv. priem $\beta$ daalt van 2.15 bij $\rho=0.096$ naar 1.46 bij $\rho=0.032$). Echter, de uitsluitingssterkte relatief aan de dichtheids-gematchte CSR-baseline blijft consistent. Vergelijkingen tussen domeinen vereisen dichtheidsmatching of het gebruik van de dichtheids-invariante $r_{excl}/\langle NN \rangle$ kalibratie.
• Interferometrische Profilometrie: Een proof-of-concept toepassing op een gecertificeerde rondheidstandaard levert $\beta = 2.00$ op, consistent met de uitsluiting waargenomen in andere domeinen.

Betekenis en Claims
Het artikel beweert een gekalibreerd meetkader te bieden voor de reproduceerbare karakterisering van kortstondige uitsluiting in 2D ruimtelijke puntprocessen. De primaire bijdragen zijn:

1. Correctie van de Baseline: Het vaststellen dat de 2D CSR-baseline $\beta \approx 0.96$ is, en niet 0, wat essentieel is om spookmatige uitsluitingsclaims te voorkomen.
2. Fysieke Interpretatie: Het bieden van een empirische $\beta$–$r_{excl}$ mapping die de abstracte Brody-parameter vertaalt naar een fysiek interpreteerbare uitsluitingsstraal, analoog aan de hard-sphere packing fraction in de vloeistof-theorie.
3. Methodologische Rigor: Het aantonen dat $\beta$ uitsluitingssterkte vangt in plaats van puntendichtheid of inbeddingsartefacten, mits dichtheidsmatching en geschikte baselines worden toegepast.
4. Fenomenologische Reikwijdte: De auteurs stellen expliciet dat $\beta$ wordt gebruikt als een empirische descriptor van ruimtelijke orde, en niet als een claim van onderliggende random-matrix symmetrieklassen (bijv. $\beta > 2$ impliceert geen hogere-orde symmetrieklassen, maar simpelweg sterkere uitsluiting).

Het kader is aangetoond toepasbaar te zijn over fysieke (gefabriceerde oppervlakken), wiskundige (priemgetallen) en optische (interferometrie) domeinen, mits de noodzakelijke correcties voor dichtheid, inbeddingsgeometrie en baseline-selectie worden toegepast. De studie claimt geen universeel uitsluitingsmechanisme te onthullen, maar levert de instrumenten om dit kwantitatief en reproduceerbaar te kwantificeren.