Neutron-Star Composition

A neutron star is not a uniform ball of neutrons. It is a sequence of phases selected by charge neutrality, weak-interaction equilibrium, and the minimum free energy at each pressure. Nuclei and electrons dominate the outer crust; neutron-rich clusters, a neutron superfluid, and electrons occupy the inner crust; and a strongly interacting liquid with at least neutrons, protons, electrons, and muons forms the core. Beyond roughly nuclear saturation density, the particle content is a research question rather than an observed inventory. This page separates robust equilibrium requirements from model-dependent claims about hyperons, condensates, and deconfined quarks.

Required background. The Einstein Equation supplies the gravitational field equation used to turn a candidate EOS into a stellar density profile.

Helpful background. Flat-Background Hilbert Stress-Energy from Metric Variation supplies the energy-density and pressure conventions; Relativistic Matter and Compact Objects supplies the thermodynamic setting; Computational Relativity supplies numerical sequence and interpolation checks.

The layer map is a thermodynamic construction

For cold catalyzed matter, the preferred state at a specified pressure minimizes the Gibbs free energy per baryon subject to conserved baryon number and vanishing macroscopic charge. In the zero-temperature, weak-field idealization, the surface is condensed matter; a real star may instead have a gaseous atmosphere, depending on temperature and magnetic field. As pressure rises, atoms are ionized and increasingly neutron-rich nuclei form a lattice in a degenerate electron gas. The exact nuclide sequence depends on nuclear masses, but the organizing calculation is well defined; the classic Baym–Pethick–Sutherland construction gives the equilibrium condition and an explicit sequence in Baym, Pethick, and Sutherland 1971, pp. 300–306 and Tables 1–2.

Near a mass density of a few $10^{11}\,\mathrm{g\,cm^{-3}}$, the neutron chemical potential reaches the rest-energy threshold and neutrons drip out of nuclei. The inner crust then contains nuclear clusters, electrons, and a neutron fluid; close to the crust–core transition, nonspherical clusters may be competitive. The locations of neutron drip and the transition are not universal constants: nuclear mass models, symmetry energy, pairing, and finite-size terms move them. The phase-by-phase derivation and its model dependence are reviewed in Chamel and Haensel 2008, §§3–4 and Tables 1–3.

Weak equilibrium and the minimal core

On times long compared with weak reactions and short compared with stellar evolution, neutrino-transparent matter obeys beta equilibrium. For the reaction $n\leftrightarrow p+e^-+\bar\nu_e$, escaping neutrinos have zero chemical potential, so

$$\mu_n=\mu_p+\mu_e, \qquad n_p=n_e+n_\mu .$$

The first equation is chemical equilibrium and the second is local charge neutrality. Muons appear when the electron chemical potential reaches $m_\mu$; thereafter equilibrium also requires $\mu_\mu=\mu_e$. These relations determine particle fractions only after an interaction model fixes the chemical potentials as derivatives of the energy density. They therefore constrain composition without choosing it.

The minimal outer-core model is a liquid of neutrons, a smaller fraction of protons, and leptons. Both nucleon species may pair, changing heat capacities, neutrino rates, vortex dynamics, and transport without changing the conserved charges. At higher density, additional particles can lower the free energy. Hyperons, pion or kaon condensates, and quark matter are possible examples, but each introduces interactions that are incompletely constrained. The thermodynamic and microphysical alternatives are compared in Oertel et al. 2017, §V.A.2.

What evidence can and cannot identify

Evidence cutoff: 2026-08-18. Laboratory nuclear data and microscopic many-body calculations constrain matter near and somewhat above saturation density. Stellar masses, radii, tidal responses, and cooling histories probe integrated consequences of the EOS. None directly labels the quasiparticles at a particular radius. A stiff nucleonic interaction and a softer phase transition followed by restiffening can therefore produce similar bulk observables.

The principal evidence axes are: the nuclear Hamiltonian and regulator; the many-body truncation and density range; the assumed particle set; the phase construction and surface tension; enforcement of causality and thermodynamic stability; and the stellar likelihood used for comparison. Evidence for new degrees of freedom comes from energetic calculations in which they lower the free energy. Competing evidence against a particular realization comes when it softens the EOS so much that stable stars cannot reach securely measured large masses. That tension is model-specific: it does not prove that every hyperonic or quark phase is absent.

The unresolved problem is to connect controlled strong-interaction theory to the several-times-saturation densities sampled by massive cores while quantifying truncation and phase uncertainties. Re-review should be triggered by a new controlled many-body calculation with quantified errors, a new laboratory constraint on isospin-rich matter, an official multimessenger data release that materially changes the EOS likelihood, or a systematic reanalysis that changes the inferred density range.

Reading a composition profile correctly

A plotted particle fraction is conditional on the chosen EOS and equilibrium assumptions. A discontinuity may represent a first-order transition; a smooth mixed-phase interval requires an explicit treatment of charge screening and surface energy. Similarly, a radial layer boundary follows only after the EOS is inserted into a stellar structure solution. Density alone does not specify radius, and central composition is not representative of the entire star.

In cold neutral neutron–proton–electron–muon matter, write $\mu_i=B_i\mu_B+q_i\mu_Q$, where $B_i$ and $q_i$ are baryon number and electric charge in proton-charge units. Weak equilibrium permits only the conserved-charge chemical potentials; charge neutrality then cancels $\mu_Q\sum_iq_i n_i$. Consequently

$$\varepsilon+p=\sum_i\mu_i n_i=\mu_B n_B.$$

This explains why a chemically complicated equilibrium can supply a one-parameter cold stellar EOS. It does not erase the composition degrees of freedom during a fast disturbance.

One reliable check is thermodynamic consistency: along an equilibrium branch, $p=\sum_i\mu_i n_i-\varepsilon$, and the sound speed derived from the same interpolation must remain real and causal. Another is limiting behavior: removing interactions should recover the appropriate degenerate-gas chemical potentials, while suppressing muons must force charge neutrality to $n_p=n_e$. These tests reject inconsistent tables; they do not select the correct high-density phase.

Check your understanding

Cold matter contains neutrons, protons, electrons, and no trapped neutrinos. Suppose the electron chemical potential is $90\,\mathrm{MeV}$. Should a zero-temperature equilibrium model contain a muon Fermi sea?

Solution

No. Beta equilibrium would require $\mu_\mu=\mu_e$, but creating the first muon costs its rest energy, $m_\mu\simeq105.7\,\mathrm{MeV}$. At $\mu_e=90\,\mathrm{MeV}$ the zero-density muon state is favored. Once $\mu_e\geq m_\mu$, muons can appear and share the negative charge, changing the relation between proton and electron fractions.

References

• Baym, Gordon, Christopher Pethick, and Peter Sutherland. “The Ground State of Matter at High Densities: Equation of State and Stellar Models.” The Astrophysical Journal 170 (1971): 299–317. DOI.
• Chamel, Nicolas, and Pawel Haensel. “Physics of Neutron Star Crusts.” Living Reviews in Relativity 11 (2008): 10. DOI.
• Oertel, Micaela, Matthias Hempel, Thomas Klähn, and Stefan Typel. “Equations of State for Supernovae and Compact Stars.” Reviews of Modern Physics 89 (2017): 015007. DOI.