Skip to main content

Chapter 10
Chandrasekhar Limit and Star Collapse

Cosmological and Global Coherence
Author: Mark Anthony Higgins
© Mark Anthony Higgins 11 October 2025
Revised Date: 23 February 2026
Audio Overview
Flow Algebra
flow "ChandrasekharLimit_StarCollapse"

# ---------- STATES ----------
state:
  # bulk properties
  star.mass.M: float                 # total mass (M_sun)
  core.mass.Mc: float                # core mass (M_sun)
  core.radius.Rc: float              # core radius (km)
  core.density.rho_c: float          # central density (g/cm^3)
  core.temperature.Tc: float         # central temperature (K)
  core.entropy.S: float              # specific entropy (k_B/baryon)
  core.Ye: float                     # electron fraction (e-/baryon)
  composition.mu_e: float            # mean molecular weight per e- (≈2 for C/O, ≈2 for He)

  # pressures (dyn/cm^2)
  pressure.degen_e: float            # electron degeneracy pressure
  pressure.therm: float              # thermal gas + radiation pressure
  pressure.total: float

  # gravity / curvature
  gravity.potential: float           # |Φ| at core
  gravity.compactness.C: float       # GM/(Rc c^2)
  boundary.integrity: float          # 0..1 global/internal coherence
  coherence.tension: float           # curvature tension density

  # neutrino channel
  neutrino.lum: float                # neutrino luminosity (erg/s)
  neutrino.cooling_rate: float

  # thresholds / references
  limit.Mch: float                   # Chandrasekhar mass (M_sun)
  limit.Mns_max: float               # max cold NS mass (M_sun) ~ 2.0–2.3
  limit.compactness_bh: float        # C_crit ~ 0.5 (order) for BH trapping

  # time
  time.t: float

# ---------- PARAMETERS ----------
params:
  G: 6.674e-8                        # cgs
  c: 2.998e10                        # cm/s
  hbar: 1.055e-27                    # erg·s
  m_e: 9.109e-28                     # g
  m_p: 1.673e-24                     # g
  a_rad: 7.5657e-15                  # radiation const (erg cm^-3 K^-4)
  M_sun: 1.989e33                    # g
  K_NR: float                        # NR e- deg. EOS coefficient
  K_REL: float                       # UR e- deg. EOS coefficient
  Mch_coeff: 5.83                    # M_ch ≈ 5.83 / μ_e^2  (in M_sun)
  cooling.coeff: float               # neutrino cooling scale
  capture.rate0: float               # e- capture baseline (Ye drain)
  dt: float

# ---------- INVARIANTS ----------
invariant:
  - limit.Mch = params.Mch_coeff / (composition.mu_e^2)
  - pressure.degen_e = P_deg(core.density.rho_c, core.Ye, regime="NR/UR")
  - pressure.therm = P_gas_rad(core.density.rho_c, core.temperature.Tc, params.a_rad)
  - pressure.total = pressure.degen_e + pressure.therm
  - gravity.compactness.C = (params.G * core.mass.Mc * params.M_sun) / (max(core.radius.Rc,1e-5)*1e5 * params.c**2)
  - quasi-hydrostatic stability requires dM/d(rho_c) > 0 and Mc < limit.Mch for WD branch
  - energy balanced (cooling/heating):
    dE_core/dt = - neutrino.lum + compressional_heating(core.density.rho_c, dV/dt)
  - boundary.integrity ∈ [0,1]; collapse triggers as integrity falls with rising compactness

# ---------- TRANSFORMS ----------
transform:

  initialise_core:
    do:
      composition.mu_e = 2.0          # e.g., C/O core
      core.mass.Mc = min(star.mass.M, 1.2)
      core.radius.Rc = 8000           # km (WD-scale start)
      core.density.rho_c = 1e6        # g/cm^3
      core.temperature.Tc = 1e8       # K (post-burning core)
      core.Ye = 0.5
      boundary.integrity = 0.995
      neutrino.lum = 0.0
      limit.Mns_max = 2.1
      limit.compactness_bh = 0.5
      time.t = 0

  update_chandrasekhar_limit:
    do:
      limit.Mch = params.Mch_coeff / (composition.mu_e**2)

  update_eos_and_pressures:
    do:
      # choose NR vs UR degeneracy regime internally
      pressure.degen_e = eos_degeneracy(core.density.rho_c, core.Ye, params.K_NR, params.K_REL)
      pressure.therm   = eos_gas_rad(core.density.rho_c, core.temperature.Tc, params.a_rad)
      pressure.total   = pressure.degen_e + pressure.therm

  quasi_hydro_adjust:
    do:
      # toy hydrostatic adjustment: R_c responds to P_total vs self-gravity
      grav_req = grav_pressure_needed(core.mass.Mc, core.radius.Rc, params.G)
      δR = response_radius(pressure.total, grav_req, core.radius.Rc)
      core.radius.Rc = max(core.radius.Rc + δR, 50)  # floor to avoid zero
      core.density.rho_c = mass_to_central_density(core.mass.Mc, core.radius.Rc)

  weak_interactions_update:
    do:
      # electron captures at high rho_c reduce Ye, undermining degeneracy support
      capture = params.capture.rate0 * f_capture(core.density.rho_c, core.temperature.Tc)
      core.Ye = max(core.Ye - capture*params.dt, 0.03)

  neutrino_cooling:
    do:
      neutrino.lum = params.cooling.coeff * f_nu(core.temperature.Tc, core.density.rho_c)
      neutrino.cooling_rate = neutrino.lum
      core.temperature.Tc = max(core.temperature.Tc - cooling_delta(neutrino.lum, core.mass.Mc)*params.dt, 1e7)

  compute_compactness:
    do:
      gravity.compactness.C = (params.G * core.mass.Mc * params.M_sun) / (core.radius.Rc*1e5 * params.c**2)
      gravity.potential = gravity_potential(core.mass.Mc, core.radius.Rc, params.G)

  integrity_drift:
    do:
      # rising compactness + Ye drain → integrity loss (onset of collapse)
      boundary.integrity = clamp(boundary.integrity - g_integrity(gravity.compactness.C, core.Ye), 0, 1)

  decide_fate:
    do:
      if core.mass.Mc < limit.Mch and boundary.integrity > 0.98:
        mark("white_dwarf_branch")
        settle_white_dwarf()
      elif core.mass.Mc >= limit.Mch:
        onset_core_collapse()

  settle_white_dwarf:
    do:
      # stable WD: set mass–radius relation (R ∝ M^{-1/3} NR; shrinks in UR)
      core.radius.Rc = wd_mass_radius(core.mass.Mc, composition.mu_e)
      core.density.rho_c = mass_to_central_density(core.mass.Mc, core.radius.Rc)
      pressure.degen_e = eos_degeneracy(core.density.rho_c, core.Ye, params.K_NR, params.K_REL)

  onset_core_collapse:
    do:
      # rapid contraction; degeneracy pressure can’t balance gravity (UR limit)
      collapse_rate = collapse_speed(gravity.compactness.C, core.density.rho_c)
      core.radius.Rc = max(core.radius.Rc * (1 - collapse_rate*params.dt), 10)
      core.density.rho_c = mass_to_central_density(core.mass.Mc, core.radius.Rc)
      boundary.integrity = max(boundary.integrity - 5e-3, 0.0)
      if core.density.rho_c > 2e14: bounce_or_continue()

  bounce_or_continue:
    do:
      # if below NS maximum mass → neutron star; else → black hole
      if core.mass.Mc <= limit.Mns_max:
        form_neutron_star()
      else:
        form_black_hole()

  form_neutron_star:
    do:
      mark("neutron_star")
      core.radius.Rc = 12               # km (canonical)
      gravity.compactness.C = (params.G * core.mass.Mc * params.M_sun) / (core.radius.Rc*1e5 * params.c**2)
      boundary.integrity = 0.9          # re-stabilised at nuclear density (new phase)
      neutrino.lum = neutrino_burst(core.mass.Mc)  # deleptonisation burst

  form_black_hole:
    do:
      mark("event_horizon")
      gravity.compactness.C = 0.5
      boundary.integrity = 0.0
      core.radius.Rc = schwarzschild_radius(core.mass.Mc, params.G, params.c)

  tick:
    do:
      time.t += params.dt

# ---------- MONITORS ----------
monitor:
  stability_track:
    observe: [core.mass.Mc, limit.Mch, core.density.rho_c]
    report: "Compares core mass to Chandrasekhar limit and follows central density."

  support_channels:
    observe: [pressure.degen_e, pressure.therm, gravity.compactness.C]
    report: "Balance of electron degeneracy + thermal support against self-gravity."

  composition_drift:
    observe: [core.Ye]
    report: "Electron capture reduces Ye, lowering degeneracy support and pushing toward collapse."

  fate_indicator:
    observe: [boundary.integrity]
    report: "High integrity → stable WD; sharp declines → collapse → NS/BH branch."

# ---------- NOTES (coherence view) ----------
note:
  - The Chandrasekhar limit emerges when **electron degeneracy support** saturates in the ultra-relativistic regime:
      M_ch ≈ (5.83 / μ_e^2) M_⊙.
  - In HST terms, degeneracy pressure is **phase-exclusion coherence**; UR electrons cannot raise pressure faster than gravity’s curvature demand.
  - As Ye falls via electron capture, the available coherence (Pauli support) decreases → **integrity loss** → collapse.
  - Outcomes:
      • **White Dwarf (WD):** Mc < M_ch, boundary re-closes at low compactness; mass–radius anti-correlation.
      • **Neutron Star (NS):** Mc ≲ M_NS,max, boundary re-stabilises at nuclear density (new coherence phase, neutrino burst).
      • **Black Hole (BH):** Mc > M_NS,max, boundary coherence ruptures; horizon forms (integrity → 0).
  - Neutrino emission is the **coherence release channel** during collapse/bounce, rapidly cooling and signalling phase transition.

What this flow captures

ChandrasekharLimit_StarCollapse describes how stars reach and cross the critical threshold where degeneracy support can no longer maintain coherence against self-gravity.
In classical astrophysics, this threshold — the Chandrasekhar limit — marks the maximum mass (~1.4 M⊙ for μₑ = 2) that a white dwarf can sustain before collapse.
In the HST framework, this limit is not simply an equation of state boundary but a transition of coherence: a shift in how the star’s internal boundary structure stores and resists curvature tension.

This flow models that process as a continuous feedback between degeneracy pressure (coherence retention), gravity (curvature compression), and electron capture (phase loss).
Collapse is triggered when boundary integrity can no longer maintain equilibrium within the curvature field.


Core intuition (HST view)

Degeneracy pressure is not an arbitrary quantum force — it is the expression of coherence exclusion: the geometry of electron wavefunctions packed into the minimal available curvature volume.
As stellar mass increases, curvature deepens; coherence modes tighten until electrons become relativistic.
At that point, their ability to sustain phase coherence (pressure) grows more slowly than the curvature demand.
This geometric bottleneck is the true origin of the Chandrasekhar limit.

When the internal boundary (core) reaches this critical curvature, integrity begins to fall — the phase structure of the star’s electrons can no longer adapt smoothly.
Electron capture reduces the electron fraction (Yₑ), further degrading coherence.
The core collapses, releasing stored curvature energy as neutrinos — the system’s way of shedding incoherent tension and seeking a new boundary equilibrium.


Mathematical structure
The classical limit:

\[ M_{\mathrm{Ch}} \approx \frac{5.83}{\mu_{e}^{2}}\,M_{\odot}. \]

emerges when relativistic degeneracy and gravitational pressure balance exactly.

Degeneracy pressure PeP_e follows the curvature-scaling relation:

\[ P_{e} \propto \rho_{c}^{5/3} \quad (\text{non-relativistic}), \qquad P_{e} \propto \rho_{c}^{4/3} \quad (\text{ultra-relativistic}). \]

\[ \text{The transition from } \tfrac{5}{3} \text{ to } \tfrac{4}{3} \text{ in the exponent reduces the stabilising feedback,} \quad \frac{dM}{d\rho} \;\rightarrow\; 0, \quad \text{marking the onset of collapse.} \]

In HST form, the same condition expresses as loss of coherence stiffness:

\[ \frac{\partial P_{\mathrm{coh}}}{\partial \rho} \;\rightarrow\; 0 \quad \Rightarrow \quad \text{boundary integrity} \;\downarrow\; \text{collapse.} \]


Interpretation in coherence physics
StageConventional viewHST coherence view
White dwarf support
Electron degeneracy pressure
Phase-exclusion
coherence maintaining
boundary integrity
Chandrasekhar limit
Maximum mass before collapse
Critical curvature
beyond which coherence
can’t adjust
Electron capture
Ye loss reduces pressure
Phase annihilation —
loss of coherence carriers
Core collapse
Gravitational instability
Curvature–coherence
imbalance; boundary
integrity rupture
Neutron star formation
Nuclear density equilibrium
New boundary phase
— re-closure of coherence
at nuclear scale
Black hole formation
Horizon formation
Global boundary rupture —
coherence collapse to
zero integrity

Boundary logic and stellar fate

White Dwarf (WD)

\[ M_{c} < M_{\mathrm{Ch}} : \]

\[ \text{Coherence integrity remains high } (\approx 0.98\text{–}1.0), \quad \text{degeneracy pressure balances curvature tension, forming a stable boundary closure,} \] \[ R \propto M^{-1/3} \quad \text{showing how the curvature field compresses coherence as mass increases.} \]

Neutron Star (NS)

For Mc≈MCh–MNS,maxM_c ≈ M_{\mathrm{Ch}}–M_{\mathrm{NS,max}}:

Collapse proceeds until nuclear densities are reached.
The boundary re-stabilises through a new coherence phase — neutrons replacing electrons as primary coherence carriers.
A neutrino burst signals the system’s rapid phase reorganisation (deleptonisation).
Boundary integrity settles near 0.9 — not perfect, but dynamically stable.

Black Hole (BH)

\[ M_{c} \approx M_{\mathrm{Ch}} - M_{\mathrm{NS,max}} : \]

Curvature exceeds all internal coherence capacities.
No phase structure can sustain the boundary; integrity → 0.
The boundary collapses into a self-contained curvature resonance — the event horizon.
In HST terms, the black hole is not a singularity but a coherence vacuum: all internal curvature locked beneath a closed boundary of zero phase communication.


Energy and coherence feedback

During collapse, curvature potential energy is redistributed as:

  • Neutrino emission: coherence discharge.
  • Gravitational waves: boundary oscillation damping.
  • Thermal energy: residual curvature smoothing.

The process obeys global coherence conservation:

Curvature energy=stored coherence tension+released neutrinos+boundary work.

Each pathway is a form of curvature feedback, redistributing tension until a stable coherence state re-emerges.


Philosophical coherence

The Chandrasekhar limit represents not just a number, but a law of structural coherence:
every boundary in nature — atomic, stellar, cosmic — has a limit beyond which phase structure cannot maintain equilibrium with curvature.
Beyond that limit, the system must transform.
Star collapse is thus not destruction but re-organisation of coherence — the boundary re-encoding itself at a deeper harmonic level.


Significance

This flow completes the stellar arc within HST cosmology:

  • Formation (curvature capture),
  • Equilibrium (degeneracy coherence),
  • Collapse (coherence phase transition),
  • Re-closure (neutron stars),
  • Re-closure (neutron stars),

It unites the quantum mechanical, relativistic, and thermodynamic descriptions into a single geometric narrative:
mass, pressure, and collapse are all expressions of how coherence resists or yields to curvature.


In essence

The Chandrasekhar limit is the curvature threshold of coherence itself.
Stars do not die; they change phase — their boundaries rewriting the memory of gravity in a new geometry.


Kireth Research
38 Thorndale Croft, Wetwang, Driffield
East Riding of Yorkshire YO25 9XZ
Designed By