Chapter 25 — Nuclear Physics of Neutron Stars: The Densest Matter in the Universe

"A neutron star is like a giant atomic nucleus — except that it is held together by gravity rather than the nuclear force, and it contains more particles than all the nuclei in a thousand suns." — Gordon Baym, Nuclear Physics A (1979)

In 1932, two years after James Chadwick discovered the neutron, Walter Baade and Fritz Zwicky proposed an audacious idea: when a massive star exhausts its nuclear fuel and collapses, the result might be a star composed almost entirely of neutrons — a "neutron star" — with the density of an atomic nucleus and a radius of only about 10 km. For over three decades, the idea seemed barely more than a theoretical curiosity. Then, in 1967, Jocelyn Bell Burnell and Antony Hewish detected a pulsing radio source with an astonishingly precise period of 1.337 seconds. Within a year, Thomas Gold identified the source as a rapidly rotating, magnetized neutron star — a pulsar. Neutron stars were real.

Today, neutron stars are recognized as the densest objects in the universe that are not black holes. They pack roughly 1.4 solar masses ($M_\odot = 1.989 \times 10^{30}$ kg) into a sphere about 12 km across. Their central densities reach several times the nuclear saturation density $\rho_0 = 2.7 \times 10^{14}$ g/cm$^3$ — the density of ordinary nuclear matter, the same density we measured in Chapter 2 using electron scattering. No laboratory on Earth can compress matter to such densities and hold it there. Neutron stars are the only places in the universe where we can study the nuclear equation of state at densities of $2\rho_0$, $5\rho_0$, perhaps $10\rho_0$.

This chapter is about what nuclear physics tells us about neutron stars, and what neutron stars tell us about nuclear physics. The two questions are inseparable. The structure of a neutron star — its mass, its radius, how it deforms under tidal forces — is determined entirely by a single function: the equation of state (EOS) of dense matter, which specifies the pressure $P$ as a function of the energy density $\varepsilon$ (or equivalently, the baryon density $n_B$). If we knew the EOS from nuclear theory, we could predict every observable neutron star property. Conversely, if we could measure neutron star properties precisely enough, we could infer the EOS — and thereby learn about nuclear matter at densities that cannot be reached any other way.

This is the central problem of neutron star physics, and it is a nuclear physics problem.

📊 Spaced Review (Chapter 3): The nuclear force saturates — each nucleon interacts primarily with its nearest neighbors, producing a nearly constant nuclear density $\rho_0 \approx 0.16$ fm$^{-3}$. The repulsive core at $r \lesssim 0.5$ fm prevents nucleons from collapsing onto one another. At the center of a neutron star, matter is compressed to several times $\rho_0$. The repulsive core of the nuclear force is what prevents the entire star from collapsing into a black hole.

📊 Spaced Review (Chapter 23): In our study of explosive nucleosynthesis, we encountered neutron star mergers as the confirmed site of r-process nucleosynthesis. The gravitational wave event GW170817 — detected on August 17, 2017, by LIGO and Virgo — was our anchor example. That same event provides some of the tightest constraints on the neutron star equation of state, through the measurement of tidal deformability. GW170817 returns in this chapter as a probe of dense matter.

📊 Spaced Review (Chapter 4): The semi-empirical mass formula treated nuclear matter as an approximately incompressible liquid with saturation density $\rho_0$. The volume, surface, Coulomb, asymmetry, and pairing terms all reappear in the physics of the neutron star crust, where nuclei exist in a sea of free neutrons and electrons. The competition between Coulomb and surface energy — familiar from the SEMF — produces the exotic "nuclear pasta" phases at the base of the crust.

25.1 Neutron Stars: The Basic Picture

Formation

A neutron star is born in the gravitational collapse of the iron core of a massive star ($M \gtrsim 8 M_\odot$), as described in Chapter 23. When the core exceeds the Chandrasekhar mass ($\sim 1.4 M_\odot$ for an iron core), electron degeneracy pressure can no longer support it. The core collapses in less than a second, reaching nuclear density and bouncing. The resulting shock wave, revived by neutrino heating, expels the outer layers of the star as a core-collapse supernova. What remains is a hot, neutrino-emitting proto-neutron star that cools over the next $\sim$10--30 seconds into a cold neutron star.

Bulk Properties

The observed properties of neutron stars are extraordinary:

A useful way to appreciate these numbers: if you could take a sugar cube of neutron star material and place it on Earth, it would weigh roughly 400 million tonnes — more than the mass of all humans alive today, compressed into a centimeter-sized cube.

A Sense of Scale: Numerical Estimates

Let us build some physical intuition with order-of-magnitude estimates.

Average density. For $M = 1.4\,M_\odot = 2.78 \times 10^{30}$ kg and $R = 12$ km:

$$\bar{\rho} = \frac{3M}{4\pi R^3} = \frac{3 \times 2.78 \times 10^{30}}{4\pi (1.2 \times 10^4)^3} \approx 4.1 \times 10^{17} \text{ kg/m}^3 = 4.1 \times 10^{14} \text{ g/cm}^3 \approx 1.5\,\rho_0$$

The average density is already above nuclear saturation density. Since the density increases toward the center, the central density must be significantly higher — typically $\rho_c \sim (4\text{--}8)\rho_0$ depending on the mass and EOS.

Surface gravity. The Newtonian surface gravity is:

$$g = \frac{GM}{R^2} = \frac{6.674 \times 10^{-11} \times 2.78 \times 10^{30}}{(1.2 \times 10^{4})^2} \approx 1.3 \times 10^{12} \text{ m/s}^2$$

This is $\sim 1.3 \times 10^{11}$ times Earth's surface gravity. An object dropped from 1 meter above the surface would hit at $\sim 1.6 \times 10^6$ m/s — about 0.5% of the speed of light.

Compactness. The dimensionless compactness parameter $\xi = GM/(Rc^2)$ measures how relativistic the star is:

$$\xi = \frac{GM}{Rc^2} = \frac{6.674 \times 10^{-11} \times 2.78 \times 10^{30}}{1.2 \times 10^4 \times (3 \times 10^8)^2} \approx 0.17$$

For a black hole, $\xi = 0.5$. The neutron star is about a third of the way to becoming a black hole. This compactness is large enough that general relativity is essential — the GR corrections to the structure equations are of order $\xi \sim 20\%$, not negligible.

Gravitational redshift. A photon emitted from the surface of the neutron star is redshifted by:

$$1 + z = (1 - 2\xi)^{-1/2} \approx (1 - 0.34)^{-1/2} \approx 1.23$$

An emission line at 6.4 keV (iron K$\alpha$) from the surface would be observed at $6.4/1.23 \approx 5.2$ keV. This redshift is routinely observed in X-ray spectra of accreting neutron stars.

Number of baryons. The total number of baryons is approximately:

$$N_B \approx \frac{M}{m_N} = \frac{2.78 \times 10^{30}}{1.67 \times 10^{-27}} \approx 1.7 \times 10^{57}$$

For comparison, the most massive known nucleus ($^{294}$Og) has 294 nucleons. A neutron star contains about $6 \times 10^{54}$ times more nucleons than the largest nucleus.

The Giant Nucleus Analogy (and Its Limits)

It is common to describe a neutron star as a "giant atomic nucleus," and the analogy has some truth: both are held at approximately nuclear density, both are composed primarily of nucleons, and both are governed by the strong interaction at the relevant scales. However, the analogy has important limits:

1. Gravity vs. the nuclear force. Ordinary nuclei are held together by the nuclear force; gravity is utterly negligible (the gravitational energy of $^{208}$Pb is $\sim 10^{-37}$ MeV). Neutron stars are held together by gravity. The nuclear pressure and the repulsive core of the nuclear force resist gravity, preventing collapse.

2. Composition. Ordinary nuclei contain roughly equal numbers of protons and neutrons (or slightly more neutrons for heavy nuclei). Neutron star matter is overwhelmingly neutron-rich: roughly 95% neutrons, 5% protons, with enough electrons and muons to maintain charge neutrality.

3. Beta equilibrium. In a nucleus, the proton fraction is set by the competition between Coulomb and asymmetry energy at fixed $A$. In a neutron star, the proton fraction is set by beta equilibrium: the condition that the beta decay $n \to p + e^- + \bar\nu_e$ and the inverse electron capture $p + e^- \to n + \nu_e$ proceed at equal rates, so there is no net change in composition. The equilibrium condition is:

$$\mu_n = \mu_p + \mu_e$$

where $\mu_i$ are the chemical potentials. At nuclear saturation density, this gives a proton fraction of only about 4--5%, far below the $Z/A \approx 0.4$ of heavy nuclei.

4. Scale. A nucleus contains at most $\sim 300$ nucleons. A neutron star contains $\sim 10^{57}$ baryons. The surface effects that dominate light nuclei (the surface term in the SEMF) are negligible for a neutron star.

🔄 Check Your Understanding: Why is neutron star matter so neutron-rich? Starting from the beta equilibrium condition $\mu_n = \mu_p + \mu_e$, argue qualitatively why the proton fraction must be small when the density is near $\rho_0$.

25.2 The Nuclear Equation of State

What Is the EOS?

The equation of state (EOS) of dense matter is the relationship between the pressure $P$ and the energy density $\varepsilon$ (or equivalently, the baryon number density $n_B$). For a zero-temperature neutron star (a good approximation for all but the first $\sim$30 seconds after formation), the EOS is a single curve $P(\varepsilon)$ that encodes all the thermodynamic information about the matter.

The EOS is the central unknown of neutron star physics. If we knew $P(\varepsilon)$ at all densities from the surface to the center, we could compute the structure of any neutron star — its mass, radius, moment of inertia, tidal deformability, and oscillation modes. Every observable property of a neutron star is a functional of the EOS.

The EOS at Nuclear Saturation Density

At densities near $\rho_0$, the EOS is reasonably well constrained by nuclear physics experiments. The key quantities are:

Saturation density and binding energy. Nuclear matter saturates at $n_0 \approx 0.16$ fm$^{-3}$ (corresponding to $\rho_0 = m_N n_0 \approx 2.7 \times 10^{14}$ g/cm$^3$) with a binding energy per nucleon of $E/A \approx -16$ MeV for symmetric nuclear matter ($n_p = n_n$). These are the same numbers that underlie the volume term of the SEMF (Chapter 4).

Nuclear incompressibility $K$. The incompressibility modulus characterizes how the energy changes with density near saturation:

$$K = 9 n_0^2 \frac{d^2 (E/A)}{d n_B^2}\bigg|_{n_0}$$

Experimentally, $K \approx 230 \pm 20$ MeV, determined from the energies of giant monopole resonances (isoscalar breathing modes) in medium and heavy nuclei. This is the "stiffness" of nuclear matter at $\rho_0$.

Symmetry energy $S$ and its slope $L$. The energy per nucleon of asymmetric matter can be expanded as:

$$\frac{E}{A}(n_B, \delta) \approx \frac{E}{A}(n_B, 0) + S(n_B) \, \delta^2 + \mathcal{O}(\delta^4)$$

where $\delta = (n_n - n_p)/n_B$ is the isospin asymmetry parameter. For symmetric matter, $\delta = 0$; for pure neutron matter, $\delta = 1$. The symmetry energy at saturation is $S(n_0) = J \approx 32 \pm 2$ MeV — this is the same quantity as the asymmetry coefficient $a_\text{sym}$ in the SEMF.

The density dependence of the symmetry energy is characterized by the slope parameter:

$$L = 3 n_0 \frac{dS}{dn_B}\bigg|_{n_0} \approx 50\text{--}70 \text{ MeV}$$

The parameter $L$ is critically important for neutron star physics because it controls the pressure of neutron-rich matter near $\rho_0$, which in turn determines the radius and crust thickness of a neutron star. A larger $L$ means a stiffer symmetry energy at high density, producing larger neutron stars.

The EOS at High Density: Terra Incognita

Above about $2\rho_0$, the EOS is genuinely uncertain. This is the regime where nuclear theory loses its firm experimental grounding, and where the most interesting physics may lurk:

1. Many-body forces. The three-nucleon forces that were essential in Chapter 3 for reproducing light nuclear spectra and the oxygen drip line become increasingly important at high density. Their effects at $3$--$5\rho_0$ are poorly constrained.

2. Hyperons. At sufficiently high density, the neutron Fermi energy rises high enough that it becomes energetically favorable to convert neutrons into hyperons — strange baryons such as $\Lambda$ ($uds$), $\Sigma^-$ ($dds$), and $\Xi^-$ ($dss$). Hyperon thresholds are typically predicted to occur around $2$--$3\rho_0$. The appearance of hyperons softens the EOS (adds new degrees of freedom that reduce the pressure), which tends to lower the maximum neutron star mass. This creates the hyperon puzzle: many EOS models that include hyperons predict maximum masses below $2 M_\odot$, in conflict with observations (Section 25.4). Resolving the hyperon puzzle requires either strong hyperon-nucleon repulsion at high density (possibly from three-body forces involving hyperons) or some mechanism that prevents hyperons from appearing.

3. Meson condensates. Bose-Einstein condensation of pions ($\pi^-$) or kaons ($K^-$) has been predicted by several groups but remains controversial. A condensate would also soften the EOS.

4. Quark matter. At sufficiently high density, nucleons may dissolve into their constituent quarks, forming a quark-gluon plasma or, at low temperature, a color superconductor. The transition from hadronic matter to quark matter could be a sharp first-order phase transition (producing a discontinuity in the density), a smooth crossover, or something in between. If a first-order transition occurs, neutron stars could have a pure quark matter core, a mixed phase region, or even be entirely composed of quark matter ("strange stars"). Lattice QCD, which works beautifully at zero baryon density, cannot currently be applied at the high baryon density relevant for neutron stars because of the fermion sign problem.

💡 Physical Insight: The equation of state at densities above $2\rho_0$ is one of the great open questions in nuclear physics. We have theory (with large uncertainties), we have observation (improving rapidly), and we have essentially no direct experimental access. Neutron star observations are, for now, the only way to probe the EOS in this regime. This is why neutron star physics is nuclear physics.

The Symmetry Energy and Neutron Star Radii: A Direct Connection

One of the most important results in neutron star physics is the connection between the nuclear symmetry energy and neutron star radii, established by James Lattimer and colleagues. The pressure of pure neutron matter near saturation density is:

$$P_{\text{NM}}(n_0) \approx \frac{n_0 L}{3}$$

where $L$ is the slope of the symmetry energy defined above. Since the neutron star crust and outer core are at densities near $n_0$, and the pressure at these densities controls the size of the star (through the TOV equation), there is an approximately linear correlation between $L$ and $R_{1.4}$:

$$R_{1.4} \approx (9.5 + 0.045 \, L) \text{ km}$$

where $L$ is in MeV. This remarkable relation means that a nuclear physics measurement of $L$ — from nuclear masses, neutron skin thicknesses, giant resonances, or heavy-ion collisions — directly predicts the radius of a neutron star 3,000 light-years away.

The PREX-II experiment. The Lead Radius Experiment (PREX-II) at Jefferson Lab measured the neutron skin thickness of $^{208}$Pb using parity-violating electron scattering:

$$\Delta r_{np} = r_n - r_p = 0.283 \pm 0.071 \text{ fm}$$

A thicker neutron skin means a larger symmetry energy slope $L$, which in turn implies larger neutron star radii. The PREX-II central value suggests $L \approx 106$ MeV, which would imply $R_{1.4} \approx 14.3$ km — somewhat in tension with the GW170817 and NICER constraints. This tension between laboratory measurements and astrophysical observations is an active area of investigation. The CREX experiment (on $^{48}$Ca) provided a complementary but somewhat discrepant result, deepening the puzzle.

📊 Spaced Review (Chapter 2): The neutron skin thickness $\Delta r_{np}$ is the difference between the RMS neutron and proton radii. We first encountered nuclear radii in Chapter 2, measured by electron scattering. The proton radius is well determined; the neutron radius is harder to measure because neutrons are electrically neutral. Parity-violating electron scattering exploits the fact that the $Z^0$ boson couples preferentially to neutrons, providing a clean probe of the neutron distribution.

25.3 The Tolman-Oppenheimer-Volkoff Equation

Hydrostatic Equilibrium in General Relativity

The structure of a non-rotating neutron star is determined by the condition of hydrostatic equilibrium: at every point, the inward pull of gravity must be exactly balanced by the outward pressure gradient. In Newtonian gravity, this gives:

$$\frac{dP}{dr} = -\frac{G \, m(r) \, \rho(r)}{r^2}$$

where $m(r) = 4\pi \int_0^r \rho(r') r'^2 dr'$ is the mass enclosed within radius $r$. But neutron stars are strongly relativistic — the escape velocity at the surface is $\sim 0.6c$, and general relativity modifies the equilibrium condition significantly.

In 1939, Richard Tolman and, independently, J. Robert Oppenheimer and George Volkoff derived the general-relativistic generalization. The Tolman-Oppenheimer-Volkoff (TOV) equation is:

$$\boxed{\frac{dP}{dr} = -\frac{G}{r^2}\left[\varepsilon(r) + \frac{P(r)}{c^2}\right]\left[m(r) + \frac{4\pi r^3 P(r)}{c^2}\right]\left[1 - \frac{2 G m(r)}{r c^2}\right]^{-1}}$$

supplemented by the mass equation:

$$\frac{dm}{dr} = 4\pi r^2 \frac{\varepsilon(r)}{c^2}$$

Here $\varepsilon$ is the total energy density (rest mass energy plus internal energy), $P$ is the pressure, and $m(r)$ is the gravitational mass enclosed within radius $r$. The equation is written in Schwarzschild coordinates.

Understanding the GR Corrections

Compare the TOV equation to its Newtonian counterpart. Three general-relativistic correction factors appear, each making gravity stronger:

1. $\varepsilon + P/c^2$ replaces $\rho$. In GR, pressure contributes to the gravitational source. This is purely relativistic — in Newtonian gravity, pressure supports against gravity but does not itself gravitate. In a neutron star, $P$ can be a significant fraction of $\varepsilon$ (at the center of a $2 M_\odot$ star, $P_c/\varepsilon_c \sim 0.1$--$0.3$), so this correction is non-negligible.

2. $m + 4\pi r^3 P/c^2$ replaces $m$. The pressure within the volume enclosed by radius $r$ also contributes to the effective gravitational mass.

3. $(1 - 2Gm/rc^2)^{-1}$ is the metric factor from the Schwarzschild geometry. As $2Gm/rc^2 \to 1$ (the condition for a black hole horizon), this factor diverges — gravity becomes irresistibly strong.

All three corrections act in the same direction: they make gravity stronger than the Newtonian prediction. This means that for a given EOS, the maximum mass of a neutron star is lower in GR than it would be in Newtonian gravity. GR makes it harder for pressure to support a star against collapse.

Solving the TOV Equation

The TOV equation is an initial-value problem. Given an EOS $P(\varepsilon)$, one specifies the central energy density $\varepsilon_c$ (or equivalently the central pressure $P_c$) and integrates outward:

1. Boundary conditions at the center ($r = 0$): $m(0) = 0$, $P(0) = P_c$, $\varepsilon(0) = \varepsilon_c$.

2. Integration: Step outward in $r$, using the TOV equation for $dP/dr$ and the mass equation for $dm/dr$. At each step, the EOS provides the relationship $\varepsilon(P)$.

3. Surface: The surface of the star is at $r = R$ where $P(R) = 0$. The total gravitational mass is $M = m(R)$.

Each choice of $\varepsilon_c$ gives a different star with a specific mass $M$ and radius $R$. Varying $\varepsilon_c$ traces out the mass-radius curve $M(R)$ for that EOS.

The Mass-Radius Relation

The mass-radius curve is the key observable prediction of any EOS model. Its qualitative features are universal:

• Low central density ($\varepsilon_c \lesssim 2\varepsilon_0$): The star has low mass and relatively large radius. For very low masses, the star approaches the behavior of a white dwarf.

• Intermediate central density ($\varepsilon_c \sim 2$--$5\varepsilon_0$): The mass increases as $\varepsilon_c$ increases, while the radius typically decreases. The curve bends toward higher $M$ and smaller $R$.

• Maximum mass: At some critical central density $\varepsilon_c^{\text{max}}$, the mass reaches a maximum $M_{\text{max}}$. This is the Oppenheimer-Volkoff limit for that EOS. Beyond this density, adding more matter actually decreases the mass because the GR corrections overwhelm the pressure support.

• Unstable branch: Stars with $\varepsilon_c > \varepsilon_c^{\text{max}}$ are dynamically unstable — a small perturbation causes them to either collapse to a black hole or expand back to the stable branch. The stability criterion is $dM/d\varepsilon_c > 0$.

The maximum mass is the single most important number predicted by an EOS. If an observed neutron star has mass $M_{\text{obs}}$, then any EOS with $M_{\text{max}} < M_{\text{obs}}$ is ruled out. This is the most powerful observational constraint on the EOS.

📜 Historical Context: Oppenheimer and Volkoff (1939)

Oppenheimer and Volkoff's 1939 paper, "On Massive Neutron Cores," used the simplest possible EOS — a free Fermi gas of neutrons — and found $M_{\text{max}} = 0.71 \, M_\odot$. This was below the Chandrasekhar mass for white dwarfs ($\sim 1.4 \, M_\odot$), which seemed paradoxical: how could a neutron star be lighter than a white dwarf? The resolution is that the free Fermi gas is unrealistically soft — it ignores the repulsive nuclear interaction at short range. Realistic EOS models, which include nuclear forces, give $M_{\text{max}} \approx 2.0$--$2.5 \, M_\odot$. The stiffening effect of the nuclear force is what allows neutron stars heavier than $\sim 1 M_\odot$ to exist.

The Polytropic EOS

For pedagogical purposes and quick estimates, a useful parametric EOS is the polytrope:

$$P = K \varepsilon^\Gamma$$

where $K$ is a constant and $\Gamma$ is the adiabatic index (a measure of stiffness). A stiffer EOS has larger $\Gamma$ (more pressure at a given density), producing larger radii and higher maximum masses. A softer EOS has smaller $\Gamma$, producing more compact stars with lower maximum masses.

For $\Gamma = 2$ (a reasonable approximation for nuclear matter), the TOV equation can be integrated numerically to give $M_{\text{max}} \approx 2\,M_\odot$ with appropriate choices of $K$. The polytropic EOS is used in the chapter project (tov_solver.py).

A Worked Example: Scaling Estimates

Before turning to numerical solutions, we can extract useful scaling relations analytically. For a polytropic EOS $P = K\varepsilon^\Gamma$ with $\Gamma = 2$, dimensional analysis of the TOV equation yields:

$$M_{\text{max}} \propto K^{1/2} \left(\frac{c^4}{G^{3/2}}\right)$$

$$R_{\text{max}} \propto K^{1/2} \left(\frac{c^2}{G^{1/2}}\right)$$

The dependence on $K$ is the same for both: doubling the stiffness parameter $K$ increases both the maximum mass and the corresponding radius by $\sqrt{2}$.

For the free neutron Fermi gas (the Oppenheimer-Volkoff original calculation), the maximum mass can be derived analytically:

$$M_{\text{max}}^{\text{OV}} = 0.71 \, M_\odot$$

This is the absolute minimum prediction: any repulsive interaction at high density pushes $M_\text{max}$ upward. The fact that observed neutron stars have masses up to $\sim 2.1\,M_\odot$ means the nuclear force increases $M_\text{max}$ by a factor of $\sim 3$ above the free gas prediction. This is the nuclear repulsive core at work.

Numerical example: Central pressure. For a $1.4\,M_\odot$ neutron star with $R = 12$ km, we can estimate the central pressure by balancing the weight of the overlying material against the pressure gradient. The Newtonian estimate gives:

$$P_c \sim \frac{GM\bar{\rho}}{R} \sim \frac{6.674 \times 10^{-11} \times 2.78 \times 10^{30} \times 4.1 \times 10^{17}}{1.2 \times 10^4} \approx 6.3 \times 10^{33} \text{ Pa}$$

Converting: $P_c \approx 6.3 \times 10^{34}$ dyn/cm$^2 \approx 35$ MeV/fm$^3$. For comparison, the pressure at nuclear saturation density in symmetric matter is approximately zero (by definition — saturation is the equilibrium density). The central pressure of the neutron star is enormous on the scale of nuclear physics.

🔄 Check Your Understanding: Why does increasing the pressure at a given density (a stiffer EOS) lead to a larger maximum mass? Would you not expect stronger pressure to make the star more compact? (Hint: the answer involves the GR correction — pressure gravitates.)

25.4 Observational Constraints on the EOS

Mass Measurements: The Heaviest Neutron Stars

The most robust observational constraint on the EOS comes from mass measurements. Neutron star masses are measured most precisely in binary pulsar systems, using the Shapiro delay and other post-Keplerian parameters from pulsar timing.

Key mass measurements:

The existence of neutron stars with masses near or above $2\,M_\odot$ is a hard constraint: any viable EOS must support at least $\sim 2\,M_\odot$. This immediately rules out many "soft" EOS models, including most models that predict a large population of hyperons without compensating stiffening mechanisms.

The measurement of PSR J0740+6620 at $2.08 \pm 0.07 \, M_\odot$ is currently the most precise measurement of a very massive neutron star. The even more massive candidate PSR J0952-0607, at $2.35 \pm 0.17 \, M_\odot$, would be extraordinarily constraining if confirmed at higher precision — it would rule out essentially all EOS models that predict $M_{\text{max}} < 2.2 \, M_\odot$.

💡 Physical Insight: A single mass measurement above $2 M_\odot$ tells us something profound about nuclear physics: the nuclear force must be sufficiently repulsive at high density to support this mass against gravitational collapse. In particular, it constrains the three-nucleon force at densities of $3$--$5\rho_0$, where no terrestrial experiment has probed. This is nuclear physics done with a radio telescope.

Gravitational Waves: Tidal Deformability

When two neutron stars spiral toward each other in a binary system, the tidal gravitational field of each star deforms its companion. The degree of deformation depends on the internal structure — and hence the EOS — of the star. A stiffer EOS produces a larger, more easily deformed star; a softer EOS produces a more compact, harder-to-deform star.

The tidal deformability $\Lambda$ quantifies this: it is the dimensionless ratio of the induced quadrupole moment to the applied tidal field, scaled by appropriate powers of the stellar radius:

$$\Lambda = \frac{2}{3} k_2 \left(\frac{R c^2}{G M}\right)^5$$

where $k_2$ is the tidal Love number (a dimensionless number between 0 and $\sim$0.15 that depends on the density profile of the star).

GW170817 — the neutron star merger event we studied in Chapter 23 as a probe of r-process nucleosynthesis — also provided the first measurement of tidal deformability from gravitational wave observations. The LIGO/Virgo analysis of the late inspiral waveform yielded a constraint on the combined tidal deformability parameter:

$$\tilde{\Lambda} = \frac{16}{13} \frac{(m_1 + 12 m_2) m_1^4 \Lambda_1 + (m_2 + 12 m_1) m_2^4 \Lambda_2}{(m_1 + m_2)^5}$$

The result, $\tilde{\Lambda} = 300^{+420}_{-230}$ (at 90% confidence), ruled out the stiffest EOS models — those predicting very large neutron star radii ($R \gtrsim 14$ km for a $1.4\,M_\odot$ star). Combined with the $2\,M_\odot$ mass constraint (which rules out the softest models), GW170817 squeezed the allowed EOS into a relatively narrow band.

📊 Numerical estimate: For a $1.4\,M_\odot$ neutron star, the GW170817 constraint implies $\Lambda_{1.4} \lesssim 800$. Using the approximate relation $\Lambda \propto R^{5\text{--}6}$, this translates to $R_{1.4} \lesssim 13.5$ km.

NICER: Simultaneous Mass and Radius Measurements

NASA's Neutron Star Interior Composition Explorer (NICER), mounted on the International Space Station since 2017, measures the X-ray emission from hot spots on the surfaces of millisecond pulsars. As the star rotates, the observed X-ray pulse profile is shaped by the star's mass and radius through gravitational light bending: photons emitted from the far side of the star can curve around and reach the observer, and the degree of bending depends on the compactness $M/R$.

By modeling the pulse profiles with sophisticated ray-tracing codes, the NICER team has obtained simultaneous mass-radius measurements for two pulsars:

The remarkable feature is that the radius is approximately the same ($\sim 12$--$13$ km) for a $1.3\,M_\odot$ star and a $2.1\,M_\odot$ star. This places strong constraints on the density dependence of the EOS: the pressure must increase steeply enough with density that even a much heavier star is only slightly more compact.

Combining Constraints: The Allowed EOS Band

No single observation pins down the EOS. But the combination of mass measurements ($M_\text{max} \geq 2\,M_\odot$), tidal deformability from GW170817 ($\tilde\Lambda \lesssim 700$--$800$), and NICER mass-radius measurements significantly narrows the range of viable EOS models. The current best constraints give, for a canonical $1.4\,M_\odot$ neutron star:

$$R_{1.4} \approx 11.5\text{--}13.5 \text{ km}$$

This corresponds to a "moderately stiff" EOS — not so soft that it cannot support $2\,M_\odot$, not so stiff that the tidal deformability exceeds the GW170817 bound.

🔄 Check Your Understanding: A very stiff EOS predicts large radii and high maximum masses. GW170817 constrains the tidal deformability (favoring smaller radii), while the $2 M_\odot$ mass measurement requires a sufficiently high maximum mass. Explain why these two constraints act in opposite directions on the stiffness of the EOS, and why their combination is so powerful.

25.5 The Neutron Star Crust

Below the thin atmosphere (a few centimeters of hydrogen or helium plasma), a neutron star has a clearly differentiated internal structure. The outermost region is the crust, which extends from the surface to a depth of about 1--2 km. Though it contains only $\sim 1$--$2\%$ of the star's mass, the crust is where the most diverse nuclear physics occurs.

The Outer Crust: Nuclei in an Electron Gas

At the surface, the matter is iron-group nuclei (primarily $^{56}$Fe) in a lattice, surrounded by a degenerate electron gas — essentially a metallic solid, but at enormous density. As the density increases below the surface, the electron Fermi energy $\mu_e$ rises. When $\mu_e$ exceeds the threshold for electron capture,

$$e^- + (Z, A) \to (Z-1, A) + \nu_e$$

the nuclei become progressively more neutron-rich. The sequence of equilibrium nuclei is determined by the condition of beta equilibrium at each density layer:

At $\rho \approx 10^6$ g/cm$^3$: $^{56}$Fe At $\rho \approx 10^9$ g/cm$^3$: $^{62}$Ni, $^{64}$Ni, ... At $\rho \approx 10^{10}$ g/cm$^3$: $^{80}$Zn, $^{82}$Ge, ... At $\rho \approx 10^{11}$ g/cm$^3$: $^{118}$Kr, $^{120}$Sr, ...

The nuclei become increasingly exotic — far more neutron-rich than anything that exists stably on Earth, and in many cases more neutron-rich than anything yet produced at radioactive beam facilities like FRIB. The outer crust is a lattice of these exotic nuclei, arranged in a body-centered cubic (BCC) structure by the Coulomb interaction, immersed in a relativistic degenerate electron gas.

The Neutron Drip Transition

At a density of $\rho_{\text{drip}} \approx 4.3 \times 10^{11}$ g/cm$^3$ — about $1.5 \times 10^{-3}\,\rho_0$ — the nuclei become so neutron-rich that the neutron separation energy $S_n$ drops to zero. Additional neutrons cannot be bound inside the nucleus and "drip" out into the surrounding space. This is the neutron drip point, and it marks the boundary between the outer crust and the inner crust.

Below the neutron drip density, matter consists of neutron-rich nuclear clusters (sometimes called "nuclei," though they are far larger and more neutron-rich than any terrestrial nucleus) coexisting with a gas of free (unbound) neutrons, all immersed in the electron gas. The free neutron gas is superfluid below a critical temperature of $\sim 10^9$--$10^{10}$ K, as a result of the attractive $^1S_0$ pairing interaction between neutrons — the same interaction that produces the pairing term in the SEMF (Chapter 4).

The Inner Crust: Exotic Nuclear Clusters

As the density increases further through the inner crust, the nuclear clusters grow in size and the neutron gas fills an increasing fraction of the volume. The proton fraction decreases; the clusters become more diffuse; the distinction between "inside" and "outside" the cluster blurs.

The physics of the inner crust connects directly to the nuclear physics of exotic neutron-rich nuclei studied at radioactive beam facilities (Chapter 10). The neutron-rich clusters in the inner crust, with $Z/A \approx 0.1$--$0.3$, resemble the very neutron-rich isotopes near the neutron drip line. Indeed, the neutron drip line on the chart of nuclides corresponds physically to the neutron drip density in the neutron star crust: both mark the point where the neutron separation energy vanishes.

Superfluid neutrons in the inner crust. The free neutrons in the inner crust pair via the attractive $^1S_0$ interaction and form a superfluid below a critical temperature of $T_c \sim 5 \times 10^9$ K. Since even young neutron stars cool below this temperature within $\sim 100$ years, the inner crust neutrons are almost always superfluid. The energy gap is:

$$\Delta(k_F) \approx 1\text{--}3 \text{ MeV}$$

peaking at a neutron Fermi momentum $k_F \approx 0.8$ fm$^{-1}$ (corresponding to $n_n \approx 0.02$ fm$^{-3}$, about $0.12\,n_0$). The superfluid neutron component plays a crucial role in pulsar glitch dynamics (Section 25.8).

At the base of the inner crust, around $\rho \approx 0.5\,\rho_0$ ($\sim 1.3 \times 10^{14}$ g/cm$^3$), the competition between Coulomb energy and nuclear surface energy drives the formation of exotic non-spherical shapes — the nuclear pasta phases.

25.6 Nuclear Pasta: When Nuclear Matter Makes Exotic Shapes

The term "nuclear pasta" was coined by Ravenhall, Pethick, and Wilson in 1983 (and independently by Hashimoto, Seki, and Yamada) because the shapes of nuclear matter at the base of the neutron star crust resemble Italian food. This is not a joke — the names are used in the research literature.

The Physics of Pasta Formation

The key to understanding pasta is the competition between two energies:

1. Coulomb energy — favors dispersing the positively charged protons as widely as possible (minimizing the electrostatic energy). This favors nuclear matter being spread out into thin sheets, rods, or a uniform distribution.

2. Nuclear surface energy — favors minimizing the surface area of the nuclear matter interface (just as surface tension makes liquid drops spherical). This favors compact, spherical shapes.

At low density (the outer crust), Coulomb energy is small and surface energy wins: nuclei are approximately spherical. At high density (the uniform core), there is no surface at all. In between, the competition produces a sequence of shapes as the volume fraction of nuclear matter $u$ increases from $\sim 0.1$ to $\sim 0.5$ and beyond:

The Pasta Sequence

The transition from spherical nuclei to uniform matter is not abrupt — it passes through these intermediate geometries. The sequence has a beautiful mathematical structure: it mirrors the classic Wigner-Seitz cell construction, and the shapes correspond to minimal-surface geometries that minimize the total (Coulomb + surface) energy at each volume fraction.

Energy Calculation: Why Spaghetti Beats Spheres

A simplified energy estimate illustrates why non-spherical shapes become favorable. Consider nuclear matter at volume fraction $u$ inside a Wigner-Seitz cell. The total energy per unit volume (relative to uniform matter) has the form:

$$e_{\text{non-uniform}} = e_{\text{surface}} + e_{\text{Coulomb}}$$

For spherical nuclei of radius $r$ in a cell of radius $R_c$ (so $u = (r/R_c)^3$):

$$e_{\text{surface}}^{\text{sphere}} \propto \frac{\sigma}{r}, \qquad e_{\text{Coulomb}}^{\text{sphere}} \propto (n_p e)^2 r^2 f_d(u)$$

where $\sigma$ is the surface tension, $n_p$ is the proton density, and $f_d(u)$ is a dimensionality-dependent filling function. Minimizing the sum with respect to $r$ gives an optimal size. The key result is that for $d$-dimensional nuclear structures (spheres: $d=3$, cylinders: $d=2$, slabs: $d=1$), the filling function changes, and above a critical volume fraction, the cylindrical ($d=2$) configuration has lower total energy than the spherical ($d=3$) configuration.

The critical volume fraction for the sphere-to-cylinder transition is approximately $u_{\text{crit}} \approx 0.2$, consistent with detailed Hartree-Fock and molecular dynamics calculations.

Observational Consequences

Nuclear pasta, if it exists, would affect several observable properties of neutron stars:

1. Crust breaking strain. The pasta layer is predicted to be the strongest material in the universe, with a shear modulus roughly $10^{30}$ erg/cm$^3$ and a breaking strain of $\sim 0.1$ — about 10 billion times stronger than steel. The breaking of the pasta layer could power starquakes and may be related to magnetar flares.

2. Neutrino opacity. During the first $\sim$30 seconds after a neutron star is born in a supernova, the crust is still hot and transparent to neutrinos. The pasta phases, with their long-range order, could scatter neutrinos coherently, affecting the neutrino signal from the next galactic supernova.

3. Gravitational wave emission. Deformations of the neutron star crust (including the pasta layer) could produce a nonzero ellipticity, making the star a source of continuous gravitational waves detectable by LIGO/Virgo.

🔄 Check Your Understanding: The nuclear pasta sequence can be understood as the competition between Coulomb and surface energy, the same competition that appears in the SEMF for ordinary nuclei. In a heavy nucleus, Coulomb energy favors fission (stretching the nucleus into elongated shapes), while surface energy favors spherical shapes. How is the pasta physics analogous to the fission barrier problem of Chapter 20?

25.7 The Neutron Star Core

Below the crust, at densities above $\sim 0.5\rho_0$, the pasta phases give way to uniform nuclear matter — the core of the neutron star. The core accounts for $\sim 99\%$ of the mass and extends from $\sim 1$ km depth to the center.

The Outer Core: npe$\mu$ Matter

The well-established composition of the outer core (at densities from $\sim\rho_0$ to $\sim 2\rho_0$) is neutron-proton-electron-muon (npe$\mu$) matter in beta equilibrium. The composition is determined by the conditions:

1. Beta equilibrium: $\mu_n = \mu_p + \mu_e$ (as derived in Section 25.1)
2. Charge neutrality: $n_p = n_e + n_\mu$
3. Muon threshold: Muons ($\mu^-$) appear when $\mu_e \geq m_\mu c^2 = 105.7$ MeV

At $\rho_0$, the typical composition is approximately 95.5% neutrons, 4% protons, and 0.5% electrons, with a negligible muon fraction. As density increases, the proton fraction gradually rises — reaching perhaps 10--15% at $3\rho_0$ — because the symmetry energy increases with density.

The protons in the outer core are predicted to form a superconductor (Type II) through the $^1S_0$ pairing channel, with a critical temperature $T_c \sim 10^9$--$10^{10}$ K. The neutrons pair in the $^3P_2$-$^3F_2$ channel (the $^1S_0$ channel becomes repulsive at these densities). Both superfluidity and superconductivity have observable consequences: they affect the cooling rate of the neutron star (through the suppression of neutrino emission processes) and the rotational dynamics (through the formation of quantized vortices in the superfluid neutron component).

The Inner Core: Unknown Territory

At densities above $\sim 2$--$3\rho_0$, we enter genuinely unknown territory. The central density of a heavy ($\sim 2\,M_\odot$) neutron star may reach $5$--$8\rho_0$, and the composition at these densities is one of the biggest open questions in nuclear physics.

Scenario 1: Nucleons only. The simplest possibility is that matter remains composed of neutrons and protons at all densities in the core, with the EOS determined by the nuclear force at high density. Modern chiral EFT calculations and quantum Monte Carlo simulations suggest that this scenario can produce maximum masses consistent with the $2\,M_\odot$ constraint, provided the three-nucleon force provides sufficient repulsion at high density.

Scenario 2: Hyperons. As discussed in Section 25.2, the appearance of $\Lambda$, $\Sigma$, and $\Xi$ hyperons would soften the EOS. The hyperon puzzle remains unresolved, though hyperon three-body forces and quark-level repulsion are active areas of research.

Scenario 3: Quark matter. At sufficiently high density, a phase transition to deconfined quark matter may occur. In the simplest picture, the quarks form a Fermi liquid of up ($u$), down ($d$), and strange ($s$) quarks — strange quark matter — that is charge-neutral and in weak equilibrium. At low temperatures, quarks can form Cooper pairs, producing a color superconductor. Several phases have been proposed:

• Color-flavor locked (CFL) phase: all three quark flavors pair symmetrically. This is the theoretically best-understood phase, expected at asymptotically high densities.
• 2SC (two-flavor superconductor): only $u$ and $d$ quarks pair. Possible at intermediate densities where the strange quark mass is significant.
• Crystalline color superconductor: analogous to the LOFF phase in condensed matter physics, with spatially varying order parameter.

The presence of quark matter in neutron star cores would manifest observationally through modifications to the mass-radius curve. A strong first-order phase transition from hadronic to quark matter would produce a "twin star" phenomenon: two stable branches of neutron stars with similar masses but different radii. Detecting twin stars would be strong evidence for a phase transition.

Scenario 4: Something else. The history of physics suggests humility. New phases of matter — not yet imagined — may exist at the extreme densities found in neutron star cores.

⚠️ Key point: We genuinely do not know what exists at the centers of neutron stars. This is not a gap in our knowledge that will be easily filled — it requires either a breakthrough in lattice QCD at finite baryon density, or sufficiently precise astronomical observations to tightly constrain the EOS from the outside.

25.8 Magnetars and Pulsars: Neutron Stars as Laboratories

Pulsars: Rotating, Magnetized Neutron Stars

Pulsars are neutron stars that emit beams of electromagnetic radiation from their magnetic poles. As the star rotates, these beams sweep across the sky like a lighthouse, producing the observed pulsed signal.

The rotation periods of pulsars range from $\sim 1.4$ ms (the fastest known: PSR J1748-2446ad) to $\sim 23$ s (the slowest known radio pulsars). The extraordinary rotational stability of millisecond pulsars — some with period derivatives as small as $\dot{P} \sim 10^{-21}$ s/s — makes them among the most precise clocks in the universe, rivaling atomic clocks.

Pulsars spin down over time as they lose rotational energy through magnetic dipole radiation and relativistic particle winds. The characteristic age is estimated as:

$$\tau_c = \frac{P}{2\dot{P}}$$

and the spin-down luminosity is:

$$\dot{E} = 4\pi^2 I \frac{\dot{P}}{P^3}$$

where $I \approx 10^{45}$ g cm$^2$ is the moment of inertia. For the Crab pulsar ($P = 33$ ms, $\dot{P} = 4.2 \times 10^{-13}$), $\dot{E} \approx 5 \times 10^{38}$ erg/s $\approx 1.3 \times 10^5 L_\odot$ — the rotational energy powers the entire Crab Nebula.

The moment of inertia $I$ depends on the EOS (through the density profile of the star). Future measurements of $I$ — potentially through spin-orbit coupling effects in binary pulsars — would provide an additional constraint on the EOS, complementary to mass and radius.

Magnetars: The Strongest Magnetic Fields

Magnetars are neutron stars with surface magnetic fields of $B \sim 10^{14}$--$10^{15}$ T — the strongest magnetic fields in the known universe. For comparison:

At these field strengths, matter behaves in qualitatively different ways:

1. Landau quantization. Electrons are confined to the lowest Landau level when $eB \gg m_e^2 c^3/\hbar$, which occurs for $B \gg 4.4 \times 10^9$ T (the quantum critical field for electrons). In magnetar fields, atoms are elongated into needle-like shapes aligned with the magnetic field.

2. Proton cyclotron frequency. The proton cyclotron energy $\hbar\omega_c = \hbar eB/(m_p c)$ reaches $\sim 1$ keV at $B \sim 10^{14}$ T, producing observable spectral features in magnetar X-ray spectra.

3. Nuclear physics effects. At $B \sim 10^{18}$ T (the QCD scale), the magnetic field would modify the nuclear force itself. Magnetar fields are still several orders of magnitude below this threshold, but some theorists have explored the effects of strong fields on the neutron star EOS.

Magnetars are observed as Soft Gamma Repeaters (SGRs) and Anomalous X-ray Pulsars (AXPs). Their energy source is the decay of the magnetic field (magnetic energy $E_B \sim B^2 R^3/6 \sim 10^{46}$--$10^{47}$ erg for $B \sim 10^{15}$ T), rather than rotation as in ordinary pulsars. Giant flares from magnetars — such as the extraordinary event from SGR 1806-20 on December 27, 2004 — can release $\sim 10^{46}$ erg in a fraction of a second, briefly outshining the entire Galaxy in gamma rays.

Neutron Star Cooling: Nuclear Physics at Finite Temperature

The cooling of a neutron star after its birth in a supernova is governed by nuclear physics processes. The star is born hot ($T \sim 10^{11}$ K) and cools by neutrino emission from the interior and photon emission from the surface. The cooling curve — surface temperature $T_s$ versus age $t$ — is sensitive to the composition and superfluid properties of the core.

Neutrino cooling processes. During the first $\sim 10^5$--$10^6$ years, neutrino emission from the core dominates. The main processes are:

1. Modified Urca process (slow): $n + n \to n + p + e^- + \bar\nu_e$ and the reverse. This requires a spectator neutron to conserve momentum. The rate scales as $T^8$.

2. Direct Urca process (fast): $n \to p + e^- + \bar\nu_e$ and $p + e^- \to n + \nu_e$. This is allowed only if momentum conservation can be satisfied: $k_F^n \leq k_F^p + k_F^e$, which requires a proton fraction $x_p > x_{\text{DU}} \approx 11\%$--$15\%$ (the exact threshold depends on whether muons are present). The direct Urca rate scales as $T^6$ and is roughly $10^5$--$10^6$ times faster than the modified Urca rate at the same temperature.

The direct Urca threshold is a sharp transition: if the proton fraction anywhere in the core exceeds $x_{\text{DU}}$, the star cools dramatically faster. The proton fraction is determined by the symmetry energy $S(n_B)$ at high density, so the cooling curve is another probe of the EOS.

Observational test. The Cassiopeia A neutron star (in the 330-year-old supernova remnant) was observed by Chandra to have a surface temperature of $T_s \approx 2 \times 10^6$ K in 2000. A decade of monitoring suggested a $\sim 4\%$ decline in surface temperature over 10 years (though this is debated). If confirmed, the rapid cooling is consistent with the onset of neutron $^3P_2$ superfluidity in the core — the superfluid pairing gap suppresses neutrino emission at $T < T_c$ but enhances it through pair-breaking-formation (PBF) processes near $T_c$.

💡 Physical Insight: A single temperature measurement of a neutron star of known age constrains the dominant neutrino emission process in the core, which in turn constrains the proton fraction, which depends on the symmetry energy at high density. Once again, an astrophysical observation becomes a nuclear physics measurement.

Glitches: Probing Superfluidity

Pulsar glitches — sudden spin-ups of typically $\Delta\Omega/\Omega \sim 10^{-9}$--$10^{-6}$ — are observed in many pulsars, most frequently in the Vela pulsar ($\sim$once every 3 years). The leading model explains glitches as the sudden transfer of angular momentum from the superfluid neutron component in the inner crust to the rigid crust.

As the star spins down due to electromagnetic braking, the crust decelerates. The superfluid component, threaded by quantized vortices, can lag behind the crust because vortices can be pinned to nuclear lattice sites in the crust. When the angular velocity lag exceeds a critical value, the vortices unpin catastrophically and transfer angular momentum to the crust, producing the observed glitch.

The magnitude of glitches constrains the fraction of the stellar moment of inertia that resides in the superfluid — and hence the extent of the inner crust, which depends on the EOS. The Vela glitch data require that at least $\sim 1.6\%$ of the total moment of inertia be in the superfluid crust, a constraint known as the glitch activity constraint.

💡 Physical Insight: Neutron stars are not merely extreme objects to marvel at — they are precision laboratories. Pulsar timing provides mass measurements with precision limited only by systematic effects. Gravitational waves measure tidal deformability. X-ray pulse profiles give mass-radius constraints. Glitches probe superfluidity. Each observation is a window into nuclear physics at densities that no accelerator can reach.

25.9 The Frontier: Current and Future Constraints

The study of neutron stars is entering a golden age, driven by several concurrent advances:

1. Third-generation gravitational wave detectors. Einstein Telescope (Europe) and Cosmic Explorer (USA), planned for the 2030s, will detect neutron star mergers at much higher signal-to-noise ratio, improving tidal deformability measurements by an order of magnitude.

2. Next-generation X-ray telescopes. The STROBE-X concept would dramatically improve the precision of mass-radius measurements from pulse profile modeling.

3. Heavy-ion collision experiments. Facilities such as FAIR (Germany), NICA (Russia), and FRIB (USA) probe the nuclear EOS at high density, though under different conditions (hot, symmetric matter) than neutron star interiors (cold, asymmetric). The CBM experiment at FAIR is specifically designed to explore the QCD phase diagram at high baryon density.

4. Chiral EFT and nuclear theory. Systematic calculations of the nuclear EOS using chiral effective field theory, with controlled uncertainty estimates, are now available up to about $2\rho_0$. These provide rigorous lower bounds on the pressure and can be combined with observational constraints at higher densities.

5. Multi-messenger astronomy. The combination of gravitational waves, electromagnetic observations (X-ray, optical, radio), and potentially neutrinos from a galactic supernova or nearby neutron star merger will provide unprecedented constraints.

The ultimate goal is to determine $P(\varepsilon)$ to a precision of $\sim 10\%$ or better across the entire density range relevant for neutron stars ($\rho_0$ to $\sim 8\rho_0$). This would resolve the question of whether exotic matter (hyperons, quarks, condensates) exists in neutron star cores, and it would provide a definitive test of nuclear many-body theory at extreme conditions.

Chapter Summary

• Neutron stars are the remnants of core-collapse supernovae, packing $\sim 1.4 M_\odot$ into a radius of $\sim 12$ km at densities up to several times nuclear saturation density $\rho_0 = 2.7 \times 10^{14}$ g/cm$^3$.

• The nuclear equation of state (EOS) — the relation $P(\varepsilon)$ — is the central unknown. Given the EOS, the Tolman-Oppenheimer-Volkoff equation determines the mass-radius relationship and the maximum neutron star mass.

• Three GR corrections in the TOV equation (pressure gravitates, volume pressure contributes to mass, metric curvature factor) all strengthen gravity relative to the Newtonian prediction.

• The maximum mass $M_\text{max}$ is the single most powerful observable: any EOS with $M_\text{max}$ below the mass of the heaviest observed neutron star ($\sim 2.1\,M_\odot$) is ruled out.

• GW170817 constrained the tidal deformability, ruling out the stiffest EOS models. NICER provided simultaneous mass-radius measurements. The combination constrains the $1.4\,M_\odot$ radius to $R_{1.4} \approx 11.5$--$13.5$ km.

• The neutron star crust consists of the outer crust (neutron-rich nuclei in a BCC lattice + electron gas), the neutron drip transition at $\rho \approx 4 \times 10^{11}$ g/cm$^3$, and the inner crust (nuclear clusters + free neutron superfluid + electrons).

• Nuclear pasta phases (gnocchi, spaghetti, lasagna, anti-spaghetti, anti-gnocchi) form at the base of the crust through the competition between Coulomb and surface energy.

• The neutron star core is npe$\mu$ matter at moderate densities ($\rho_0$--$2\rho_0$) and unknown matter at higher densities. Hyperons, meson condensates, and quark matter are all possibilities, but none are confirmed.

• Magnetars ($B \sim 10^{14}$--$10^{15}$ T) and pulsars provide observational probes of the EOS, superfluidity, and extreme magnetic field physics. Pulsar glitches constrain the superfluid fraction and the extent of the crust.

• The field is entering a golden age: gravitational wave detectors, X-ray telescopes, and nuclear theory are converging toward a precision determination of the EOS.

Next: Chapter 26 — Nuclear Energy: Reactors, Fuel Cycles, and the Future of Fission