index

Binding Energy, Critical Radii, & Information Maintenance Tax

· 17min

Every atom in a bound system pays continuous energy to maintain its structure against entropy. Nuclear binding energy, chemical bonds, gravitational self-energy—these represent ongoing thermodynamic costs, an information maintenance tax that all organized systems pay to resist thermal randomization.

White dwarfs accreting toward the Chandrasekhar limit reveal the mechanism. Geometric gravitational compression increases by a factor of 2.2 while organizational complexity explodes by a factor of 2200—a 1000-fold disparity suggesting information bankruptcy drives collapse. Independent analysis of two white dwarf catalogs identifies r/rcrit103r/r_{\mathrm{crit}} \approx 10^3 as a discrete phase transition boundary: objects crossing this threshold appear systematically older than mass-and-temperature matched references, with the Montreal sample showing +103 Myr excess aging at 3.59σ significance and Gaia confirming at 34.6σ. The effect concentrates below r/rcrit1000r/r_{\mathrm{crit}} \approx 1000 and vanishes above, consistent with sector saturation forcing accelerated entropy production. The correspondence between binding energies and thermodynamic maintenance costs spans 57 orders of magnitude, from elementary particles to black holes.

The Thermodynamic Tax

Landauer’s principle establishes the minimum energy cost for information maintenance against thermal noise—kBTln2k_B T \ln 2 per bit, where kBk_B is Boltzmann’s constant and TT is temperature1. Binding energy across all scales represents this ongoing thermodynamic maintenance against entropy. It quantifies the continuous work required to keep structure coherent. Nuclear binding maintains quark confinement, electromagnetic binding maintains atomic configurations, gravitational binding maintains bulk coherence—each through sustained energy expenditure against thermal randomization.

Any organized system requires bits to specify its configuration. For a bound system with NN particles of mass mm confined to radius RR with characteristic momentum pp, the bit count equals,

Nbits=log2(Ωh3N),N_{\text{bits}} = \log_2\left(\frac{\Omega}{h^{3N}}\right),

where Ω\Omega is the accessible phase space volume and hh is Planck’s constant. Each particle contributes (Rp/)3\sim (Rp/\hbar)^3 states, so the total scales as,

Ω(Rp)3N,\Omega \sim \left(\frac{Rp}{\hbar}\right)^{3N},

where pGMm/Rp \sim \sqrt{GMm/R} for gravitationally bound systems at virial equilibrium.

Thermal and quantum fluctuations constantly randomize configurations, requiring continuous energy expenditure,

Em=Nbits×kBTln2.E_m = N_{\text{bits}} \times k_B T \ln 2.

This maintenance energy equals binding energy for the dominant force at each scale. Dividing by rest-mass energy defines the universal maintenance fraction

ξEmMc2=rcritr,\xi \equiv \frac{E_m}{Mc^2} = \frac{r_{\mathrm{crit}}}{r},

where rcritr_{\mathrm{crit}} is the bankruptcy radius set by the dominant force—QCD confinement, the classical electron radius, the Chandrasekhar scale, or the Schwarzschild radius, depending on regime. The complete energy budget partition follows,

Em+Ea=Mc2,E_m + E_a = Mc^2,

where Ea=Mc2(1ξ)E_a = Mc^2(1 - \xi) represents energy available for work beyond maintenance. For gravitational horizons specifically, equating EmE_m with binding energy GM2/RGM^2/R and incorporating the Landauer-Bekenstein-Hawking factor of 2 from dimensional reduction at horizons identifies the Schwarzschild radius RS=2GM/c2R_S = 2GM/c^2 as rcritr_{\mathrm{crit}} for that rung.

The Geometric Origin of Maintenance

The thermodynamic tax has a geometric origin. Any organized system enforces feasibility constraints on its dynamics: electron shells confine electrons to quantized orbitals, nuclear binding confines quarks within hadrons, gravitational virial relations confine matter within characteristic radii. These constraints project the system’s dynamics onto admissible configurations.

When constraints are state-independent—the same restriction everywhere in phase space—the projected dynamics derive from a scalar potential. The system descends toward equilibrium along a well-defined gradient. When constraints vary with state—when the admissible directions depend on where the system currently sits—the projection generically introduces curl into the effective dynamics.

Curl measures irreducible circulation: work that must be continuously supplied because no global potential exists. The curl-maintenance functional,

Mcurl=12dα2dV,\mathcal{M}_{\mathrm{curl}} = \frac{1}{2} \int |d\alpha|^2 \, dV,

where α=F\alpha = F^\flat is the 1-form dual to the correction field FF, quantifies this cost. On compact manifolds with trivial first cohomology, a spectral lower bound ensures that non-integrable projections carry an irreducible maintenance floor proportional to the projection defect magnitude.

Binding energies are curl-maintenance costs—the thermodynamic price of enforcing non-integrable constraints against entropy. Nuclear binding maintains quark confinement against QCD fluctuations. Electromagnetic binding maintains electron configurations against thermal noise. Gravitational binding maintains bulk coherence against dispersal. Each represents continuous energy expenditure to sustain state-dependent feasibility projections that lack global potentials.

Force-Specific Bankruptcy Radii

Each fundamental force imposes a characteristic bankruptcy—the scale where maintenance costs equal total energy. For the strong force at the QCD confinement scale,

rQCD=ΛQCD1 fm,r_{\text{QCD}} = \frac{\hbar}{\Lambda_{\text{QCD}}} \approx 1 \text{ fm},

where ΛQCD200\Lambda_{\text{QCD}} \approx 200 MeV sets quark confinement energy. Electromagnetic binding reaches bankruptcy at the classical electron radius,

rEM=ke2mec22.8×1015 m,r_{\text{EM}} = \frac{ke^2}{m_ec^2} \approx 2.8 \times 10^{-15} \text{ m},

where kk is Coulomb’s constant, ee is elementary charge, mem_e is electron mass, and cc is light speed. Gravity carries two bankruptcy scales: the Chandrasekhar scale for the electron-degeneracy rung, and the Schwarzschild radius

RS=2GMc2R_S = \frac{2GM}{c^2}

for the horizon rung.

The maintenance fraction ξ\xi quantifies the fraction of a system’s energy budget devoted to curvature maintenance (§5 of Triadic Tension, Decade Symmetry, & Dissipation Flow). Observed values cluster near decade-spaced walking plateaus: elementary particles at ξ106\xi \sim 10^{-6}, atoms at 103\sim 10^{-3}, molecules at 102\sim 10^{-2}, and biological systems at 101\sim 10^{-1}. The decade spacing reads directly off the KK semigroup 6,10,15\langle 6, 10, 15\rangle. The generator 10 is lcm(2,5)\mathrm{lcm}(2, 5) from the icosahedral axis orders, equivalently the C10C_{10} factor partitioning the ground-state curvature I=4πφ2I = 4\pi\varphi^2 into the per-shell coupling u=I/103.29u^* = I/10 \approx 3.29. Each decade in ξ\xi is one complete C10C_{10} shell cycle.

The bottom rung’s placement at ξ106\xi \sim 10^{-6} is a boundary condition. The formal UV fixed point (ξ,D)=(0,3)(\xi, \mathcal{D}) = (0, 3) is spectrally perfect and physically inaccessible. QED vacuum coupling—self-energy, vacuum polarization, and Fermi-golden-rule coupling to photon and phonon modes—lifts the flow off (0,3)(0, 3) at the lowest accessible decade, anchoring the ladder at ξα2me/M\xi \sim \alpha^2\sqrt{m_e/M}. The fine-structure constant α\alpha enters as an input, setting the anchor scale. Below ξ=0.5\xi = 0.5 the full KK spectrum is admissible; above ξ=0.5\xi = 0.5 the gap structure of 6,10,15\langle 6, 10, 15\rangle carves the compact-remnant-to-horizon rungs, culminating at the Coxeter boundary lmax=30l_{\max} = 30.

Forcercritr_{\text{crit}}Typical r/rcritr/r_{\text{crit}}Maintenance Cost
QCD1 fm~100~1% rest mass
EM2.8×10152.8 \times 10^{-15} m~10510^5~10510^{-5} rest mass
Gravity (electron-degeneracy rung)Chandrasekhar scalevaries with M/MChM/M_{\mathrm{Ch}}ξ×\xi \times rest mass
Gravity (horizon rung)RS=2GM/c2R_S = 2GM/c^2variesξ×\xi \times rest mass

QCD and electromagnetic forces operate far from bankruptcy, enabling stable structures with minimal overhead. These bankruptcy radii mark where accumulated curvature has consumed all available capacity—the system has exhausted integrable degrees of freedom, and maintenance cost equals total energy budget. At this threshold, no further organizational complexity can be sustained without catastrophic reorganization.

Information Bankruptcy Trajectory

White dwarfs accreting mass toward the Chandrasekhar limit MCh=1.36MM_{\text{Ch}} = 1.36 M_{\odot}2 provide a natural laboratory for studying information bankruptcy.

The complexity multiplier quantifying the overhead beyond baseline requirements follows,

C(ξ,D)=φ2D2×(1ξ)u,\mathcal{C}(\xi,\mathcal{D}) = \varphi^{2^{\mathcal{D}-2}} \times (1-\xi)^{-u^*},

where φ=(1+5)/2\varphi = (1+\sqrt{5})/2 is the golden ratio and u=4πφ2/103.29u^* = 4\pi\varphi^2/10 \approx 3.29 is the coupling constant derived in constraint geometry. This contains two competing terms. The dimensional factor φ2D2\varphi^{2^{\mathcal{D}-2}} decreases mildly as effective dimension D\mathcal{D} drops from 3 toward 2, representing reduced interference in lower dimensions. The bankruptcy factor (1ξ)u(1-\xi)^{-u^*} diverges catastrophically as maintenance fraction ξ\xi approaches unity.

Numerical integration from stable white dwarfs through collapse yields the trajectory (using constant radius R5000R \approx 5000 km from electron degeneracy pressure):

MM (M)(M_{\odot})rcrit/rr_{\mathrm{crit}}/r (Schwarzschild)ξ\xiD\mathcal{D}φ2D2\varphi^{2^{\mathcal{D}-2}}(1ξ)u(1-\xi)^{-u^*}C(ξ,D)\mathcal{C}(\xi,\mathcal{D})ξ×C\xi \times \mathcal{C}Status
0.603.6×1043.6 \times 10^{-4}0.0662.972.611.243.20.21Stable
1.006.0×1046.0 \times 10^{-4}0.272.872.522.907.32.0Normal
1.177.0×1047.0 \times 10^{-4}0.462.782.425.6613.76.3Anomaly
1.307.8×1047.8 \times 10^{-4}0.632.702.3512.429.118.3Critical
1.358.0×1048.0 \times 10^{-4}0.972.532.15229492477Collapse

The trajectory reveals the mechanism. Geometric compression RS/RR_S/R rises by a factor of 2.2 from M=0.60M = 0.60 to 1.35M1.35 M_{\odot}—mild gravitational strengthening. Organizational complexity ξ×C\xi \times \mathcal{C} rises by a factor of 2200 over the same range. The 1000-fold disparity identifies the organizational channel as the dominant driver of instability, with gravitational compression contributing a small share of the rising cost. The dimensional factor φ2D2\varphi^{2^{\mathcal{D}-2}} drops modestly from 2.61 to 2.15 as D\mathcal{D} flows from 2.97 to 2.53—barely 20% variation. The bankruptcy factor (1ξ)u(1-\xi)^{-u^*} generates the explosion: from 1.24 at stable masses to 229 near collapse, a 185-fold increase driven by the same coupling constant u3.29u^* \approx 3.29 that governs the β\beta-function’s RG flow.

Basin of Attraction Entrance

The observational anomaly at r/rcrit=103r/r_{\mathrm{crit}} = 10^3—with rcrit=RSr_{\mathrm{crit}} = R_S for gravitational compactness at the Schwarzschild-horizon rung—corresponds to M1.17MM \approx 1.17 M_{\odot} where ξ=0.46\xi = 0.46 and (1ξ)u=5.66(1-\xi)^{-u^*} = 5.66. This marks the boundary where thermodynamic bankruptcy becomes inevitable rather than merely possible—the entrance to the basin of attraction toward organizational collapse. The basin language is precise: as formalized in Finite Residence, Feasibility Projections, and Quartic Transport, it suffices that the system’s dynamics land within the basin of attraction of the feasibility operator Πx\Pi_x, with exact feasibility enforced by the projector rather than the proposal dynamics3. Below r/rcrit103r/r_{\mathrm{crit}} \approx 10^3, proposals that would maintain the current organizational state no longer converge under the feasibility map—the system has exited the projector’s attraction basin, and no correction schedule can restore convergence without structural reorganization. The residual-gradient correction method from constrained diffusion for ACOPF3 provides the operational analog: once the infeasibility residual R(x)\mathcal{R}(x) exceeds the basin boundary, gradient corrections x=xλxR(x)x' = x - \lambda\,\nabla_x \mathcal{R}(x) no longer converge to a feasible state, and the system must jump to a qualitatively different organizational mode.

Before this threshold, complexity overhead grows slowly—a factor of 3.6 from M=0.6M = 0.6 to 1.17M1.17 M_{\odot}. After crossing r/rcrit=103r/r_{\mathrm{crit}} = 10^3, overhead explodes—a factor of 36 from M=1.17M = 1.17 to 1.35M1.35 M_{\odot}. The (1ξ)u(1-\xi)^{-u^*} value of 5.66 at the anomaly threshold represents the onset of nonlinear divergence. Systems maintain ξ<0.5\xi < 0.5 through moderate overhead factors. Beyond ξ0.5\xi \sim 0.5, divergence accelerates.

Analysis of 7,515 Gaia DR3 white dwarfs4 reveals that objects in the anomaly zone (r/rcrit[700,1000)r/r_{\mathrm{crit}} \in [700, 1000)) appear systematically older than mass- and temperature-matched references, with a mean excess of +76+76 Myr (14.5σ14.5\sigma, Wilcoxon p=4.27×1084p = 4.27 \times 10^{-84}, 95% bootstrap CI: [+67,+85][+67, +85] Myr). Binned analysis reveals a monotonic compactness gradient: +345+345 Myr at r/rcrit=500r/r_{\mathrm{crit}} = 500700700 (using relativistic ONe-core tracks for the ultra-massive subpopulation), +108+108 Myr at 700700900900, +6+6 Myr at 90090011001100, and negligible anomaly above r/rcrit=1100r/r_{\mathrm{crit}} = 1100. The full analysis, including dual methodology for different mass regimes and systematic checks on core composition, is developed in Repair as Local Optimization in Constraint Geometry.

The discrete, threshold-like character of the effect—rather than a continuous gradient with mass—suggests a phase transition in constraint geometry. Once a white dwarf crosses into the r/rcrit<1000r/r_{\mathrm{crit}} < 1000 regime, it enters a qualitatively different organizational basin characterized by accelerated entropy production and faster evolution through state space. This is consistent with sector saturation in the constraint eigenvalue decomposition—at r/rcrit103r/r_{\mathrm{crit}} \approx 10^3, one organizational sector reaches capacity, forcing reconfiguration whose thermodynamic cost manifests as ~100 Myr of accelerated aging.

Independent observational work has recently identified anomalous behavior in ultra-massive white dwarfs concentrated near the same compactness scale. A 2026 analysis of Gaia-selected white dwarfs along the Q branch5 reports extended cooling delays localized to this regime, inferred from kinematic clustering and photometric evolution. The conventional stellar-evolution interpretation preserves the discreteness and localization of the effect, mirroring the threshold behavior identified here. The convergence of independent observables—cooling dwell time on the Q branch and apparent age excess relative to matched controls—identifies a structural transition in white dwarf evolution, reinforcing r/rcrit103r/r_{\mathrm{crit}} \sim 10^3 as a physically meaningful regime boundary across independent observational parameterizations.

The Discontinuous Jump to Neutron Degeneracy

White dwarfs reach the (ξ=1,D=2)(\xi=1, \mathcal{D}=2) black hole fixed point through a discontinuous organizational jump forced by information bankruptcy. At MMChM \approx M_{\text{Ch}}, the system reaches ξ0.97\xi \approx 0.97, D2.5\mathcal{D} \approx 2.5 with maintenance cost ξ×C477\xi \times \mathcal{C} \approx 477. This overhead exceeds sustainable levels, triggering catastrophic reorganization—the white dwarf jumps to neutron degeneracy at ξ0.3\xi \sim 0.3, D2.5\mathcal{D} \sim 2.5 with complexity ξ×C2.3\xi \times \mathcal{C} \approx 2.3.

The organizational complexity drops by a factor of 207, requiring massive information restructuring. The energy cost of this reorganization follows from Landauer’s principle: counting the bits required to reorganize phase space information from electron to neutron degeneracy (ΔNbits4.5×1058\Delta N_{\text{bits}} \approx 4.5 \times 10^{58}) at the shock temperature T109T \sim 10^9 K gives a transition energy of Etrans=4.3×1044E_{\text{trans}} = 4.3 \times 10^{44} J—matching observed Type Ia supernova energies6 to within measurement uncertainty. The full derivation, which requires only four observational inputs (Chandrasekhar mass, white dwarf radius, neutron star radius, and shock temperature) and no parameters from the constraint geometry, is developed in the companion paper. The white dwarf collapses because maintaining organizational complexity at ξ0.97\xi \approx 0.97 requires infinite energy through the (1ξ)u(1-\xi)^{-u^*} divergence.

The (ξ=1,D=2)(\xi=1, \mathcal{D}=2) black hole state requires additional compression beyond neutron star density, achievable only by exceeding the Oppenheimer-Volkoff limit7. White dwarf collapse represents information bankruptcy at electron degeneracy—a different failure mode with lower organizational overhead than neutron star collapse to black holes. This mechanism operates across all scales: ionization at atomic bankruptcy (~100 eV per particle), molecular dissociation at chemical limits (10 eV per bond), stellar collapse at gravitational thresholds (104410^{44} J for solar masses).

The white-dwarf-to-neutron-star jump is the first of a two-rung gravitational cascade. The second rung is the neutron-star-to-black-hole transition at the E8E_8 Coxeter threshold. Each rung is driven by a distinct failure mode. The first is information bankruptcy, governed by the (1ξ)u(1-\xi)^{-u^*} divergence and priced by Landauer reorganization at shock temperature (the Type Ia 4.3×10444.3 \times 10^{44} J above). The second is geometric bankruptcy, governed by the geometric slack mapping lmax=h(r/rcrit1)=h(1/ξ1)l_{\max} = \lfloor h(r/r_{\mathrm{crit}} - 1)\rfloor = \lfloor h(1/\xi - 1)\rfloor—when ξ\xi reaches 0.50.5, the accessible KK spin cutoff drops to the Coxeter number h=30h = 30. PSR J0740+6620 at ξ=0.493\xi = 0.493 lands on this boundary exactly. Beyond it, the system requests lmax<30l_{\max} < 30, descending into the 15 forbidden levels of the 6,10,15\langle 6, 10, 15\rangle gap set; the internal geometry cannot fulfill the request and executes a topological phase transition to a two-dimensional horizon. Both rungs are instances of a single principle—an infinite information request against a finite substrate, structurally parallel to the 1900 electromagnetic ultraviolet catastrophe, with dimensional cutoff (gravity, 3D \to 2D) replacing spectral cutoff (electromagnetism, hνh\nu). The Bekenstein–Hawking area law is the geometric circuit breaker that closes the cascade.

Information Processing Regimes

Black holes occupy the IR fixed point of the β\beta-function at (ξ=1,D=2)(\xi = 1, \mathcal{D} = 2), where all available energy is dedicated to maintaining horizon structure—maximum entropy as pure curvature maintenance with zero excess capacity8. The dimensional reduction to D=2\mathcal{D} = 2 at this fixed point is consistent with the holographic scaling of Bekenstein-Hawking entropy with boundary area rather than volume. All other organized systems operate at ξ<1\xi < 1, with residual capacity (1ξ)(1-\xi) available for organizational dynamics beyond pure pattern preservation.

Recent observations of the tidal disruption event AT2020afhd provide an empirical window into this regime9. The system exhibits synchronized X-ray and radio oscillations persisting over many dynamical times, consistent with Lense–Thirring precession driven by black hole spin. Angular momentum information remains coherently encoded in boundary geometry across extended timescales—behavior consistent with the information maintenance saturation expected as ξ1\xi \to 1. That such coherence emerges naturally in the strong-field limit suggests the dimensional reduction and boundary processing architecture are generic features of gravitational information storage.

Critical mass phenomena mark bankruptcy thresholds. White dwarfs approaching the Chandrasekhar limit see the maintenance fraction diverge through the saturating form,

ξ(M)=ξ0(1M/MCh)u1+ξ0(1M/MCh)u,\xi(M) = \frac{\xi_0(1 - M/M_{\text{Ch}})^{-u^*}}{1 + \xi_0(1 - M/M_{\text{Ch}})^{-u^*}},

with baseline ξ0=102\xi_0 = 10^{-2}—the molecular-plateau rung of the C10C_{10} ladder. The Coulomb-crystallized white-dwarf core carries condensed-matter organizational overhead, pinning ξ0\xi_0 to this rung. This prevents unphysical ξ>1\xi > 1 while capturing divergence as MMChM \to M_{\text{Ch}}. Combined with the complexity multiplier, this generates the observed 2200-fold increase in organizational overhead over the mass range from stable white dwarfs to collapse. Conversely, black holes face inverted criticality through Hawking evaporation8,

TH=c38πGMkB1M,T_H = \frac{\hbar c^3}{8\pi G M k_B} \propto \frac{1}{M},

where THT_H is Hawking temperature. Below approximately 101210^{12} kg, rising temperature creates thermal noise overwhelming error correction, triggering runaway evaporation.

Binding Energy as Universal Maintenance Cost

Binding energy is the thermodynamic signature of non-integrable feasibility projection. The universal maintenance fraction ξrcrit/r\xi \equiv r_{\mathrm{crit}}/r quantifies it operationally across all forces, and the coupling constant u=4πφ2/10u^* = 4\pi\varphi^2/10 controls bankruptcy divergence through (1ξ)u(1-\xi)^{-u^*}. The identity holds over 57 orders of magnitude — QED vacuum dressing at ξ106\xi \sim 10^{-6}, biological metabolism at ξ101\xi \sim 10^{-1}, gravitational horizons at ξ=1\xi = 1 — a single thermodynamic tax collected by each fundamental force at its characteristic rcritr_{\mathrm{crit}}.

Under this reading, binding energy is the ongoing cost of structural coherence. The Landauer bit-counting at Chandrasekhar mass recovers Type Ia supernova energies to within observational scatter, and the discrete cooling-age threshold at r/rcrit=103r/r_{\mathrm{crit}} = 10^3 in Gaia DR3 white dwarfs marks the boundary where the bankruptcy factor enters nonlinear divergence. Both point to the β\beta-function’s IR fixed point at ξ=1\xi = 1 as the thermodynamic endpoint of bound structure, with every system below it paying maintenance set by the same coupling.

Limitations and Falsifiability

The binding-as-maintenance interpretation rests on two empirical claims and one theoretical identification. If the white dwarf cooling anomaly at r/rcrit103r/r_{\mathrm{crit}} \approx 10^3 is explained by conventional physics—enhanced neutrino cooling, modified equations of state near the Chandrasekhar limit, or systematic biases in progenitor populations—then the sector saturation interpretation loses its primary evidence, though the Landauer energy calculation would remain independent. If the Type Ia supernova energy match (4.3×10444.3 \times 10^{44} J) is coincidental rather than causal—if the agreement between Landauer bit-counting and observed energies reflects fortuitous cancellation of errors in the four input parameters—then the theory’s strongest quantitative prediction dissolves. The theoretical identification of binding energy with curl-maintenance cost depends on the constraint geometry’s claim that state-dependent projections introduce irreducible circulation (Theorem 4 of the self-correction paper); if that theorem’s hypotheses fail on the relevant physical manifolds, the geometric foundation weakens. The methodology appendix below details the empirical analysis and its limitations.

Appendix: White Dwarf Anomaly Methodology

The empirical analysis supporting the r/rcrit103r/r_{\mathrm{crit}} \approx 10^3 threshold uses 7,515 white dwarfs from the Gaia DR3 catalog4, pre-filtered for r/rcrit[500,1500]r/r_{\mathrm{crit}} \in [500, 1500] and M1.1MM \geq 1.1 \, M_\odot (hydrogen-atmosphere, thick-envelope fits). Cooling ages are derived from Montréal DA thick-hydrogen-envelope tracks (0.200.201.30M1.30 \, M_\odot), supplemented by La Plata relativistic ONe-core and CO-core GR tracks extending coverage to 1.369M1.369 \, M_\odot and 1.382M1.382 \, M_\odot respectively.

Two statistical methods address different mass regimes. For r/rcrit[700,1000)r/r_{\mathrm{crit}} \in [700, 1000) versus the reference zone [1000,1500][1000, 1500], kk-nearest-neighbor matching (k=10k = 10) in standardized (M,logTeff)(M, \log T_{\mathrm{eff}}) space constructs cooling-age residuals Δti=tit~ref\Delta t_i = t_i - \tilde{t}_{\mathrm{ref}}, where Δt>0\Delta t > 0 indicates the anomaly-zone object appears older. For r/rcrit[500,700)r/r_{\mathrm{crit}} \in [500, 700), where all objects have M>1.32MM > 1.32 \, M_\odot and no reference population exists at comparable mass, a direct Mann-Whitney comparison against the 700–900 bin replaces kNN matching.

The anomaly shows a monotonic compactness gradient (bins in r/rcritr/r_{\mathrm{crit}}): +345+345 Myr median excess at 500–700 (p=1.5×106p = 1.5 \times 10^{-6}), +108+108 Myr at 700–900 (p=4.6×1029p = 4.6 \times 10^{-29}), +6+6 Myr at 900–1100, and null above 1300. Globally across 700–1000, the mean residual is +76+76 Myr (14.5σ14.5\sigma, Wilcoxon p=4.27×1084p = 4.27 \times 10^{-84}). The Spearman correlation ρ(M2,Δt)=+0.350\rho(M^2, \Delta t) = +0.350 (p=4.93×1030p = 4.93 \times 10^{-30}) confirms that more compact objects show larger positive residuals. The result is robust to core-composition assumptions: CO-core tracks yield +697+697 Myr for the 500–700 bin (p=5.9×1010p = 5.9 \times 10^{-10}). Full methodological details, including track construction, interpolation procedures, and systematic checks, are given in the repair-as-local-optimization companion post.

High-mass white dwarfs may have systematic differences in progenitor populations or formation channels, and the causal mechanism cannot yet be distinguished from alternative explanations such as enhanced neutrino cooling, crystallization energy release, or GR corrections to the cooling rate near the Chandrasekhar limit.

Footnotes

  1. Landauer, R. (1961). Irreversibility and heat generation in the computing process. IBM Journal of Research and Development, 5(3), 183-191.

  2. Chandrasekhar, S. (1931). The maximum mass of ideal white dwarfs. Astrophysical Journal, 74, 81-82.

  3. Shekhar, S., Karn, A., Keshav, K., Bansal, S., & Pareek, P. (2026). “Fast Diffusion with Physics-Correction for ACOPF.” arXiv

    .03020. https://arxiv.org/abs/2602.03020 2

  4. Dufour, P., et al. (2017). The Montreal White Dwarf Database. Proceedings of the 20th European White Dwarf Workshop, 509, 3. See also Gentile Fusillo, N. P., et al. (2021). A catalogue of white dwarfs in Gaia EDR3. Monthly Notices of the Royal Astronomical Society, 508(3), 3877-3896. 2

  5. Ould Rouis, L. B., Hermes, J. J., Guidry, J. A., et al. (2026). White Dwarf Merger Remnants with Cooling Delays on the Q Branch Lack Strong Magnetism. Astrophysical Journal, submitted. arXiv

    .02670.

  6. Hillebrandt, W., & Niemeyer, J. C. (2000). Type Ia supernova explosion models. Annual Review of Astronomy and Astrophysics, 38(1), 191-230.

  7. Oppenheimer, J. R., & Volkoff, G. M. (1939). On massive neutron cores. Physical Review, 55(4), 374.

  8. Hawking, S. W. (1975). Particle creation by black holes. Communications in Mathematical Physics, 43(3), 199-220. 2

  9. Pasham, D. R., et al. (2025). Synchronized X-ray and radio variability from the tidal disruption event AT2020afhd consistent with Lense–Thirring precession. Science Advances, 11(49), eady9068.