Relativistic Matter and Compact Objects

Relativistic matter turns local constitutive physics into observable compact objects. A stress tensor and conservation laws are only the beginning: one must decide which currents exist, how pressure and temperature depend on state, how radiation or magnetic fields exchange energy with matter, and whether the resulting equations are stable, causal, and valid on the scales of interest. Those choices then determine stellar equilibrium, modes, collapse, merger remnants, and the inferences drawn from timing, X-ray, gravitational-wave, and laboratory data.

This volume develops that chain without hiding its weak links. It begins with continuum and kinetic matter, builds the closures used for dense stars, follows them through equilibrium and dynamics, and ends by asking which material properties observations actually constrain. Established derivations are separated from computational approximations and from conclusions that remain model-dependent or disputed.

Helpful background. Hilbert stress–energy from metric variation supplies the matter tensor and conservation law; the Einstein equation supplies the geometric response to that tensor. Individual chapters link the additional preparation they use, including special relativity, differential forms, ordinary differential equations, spacetime symmetries, and tetrads.

From a local material model to a compact object

A material model is not just a relation $p=p(\varepsilon)$. In the simplest perfect-fluid description,

$$T^{\mu\nu}=(\varepsilon+p)u^\mu u^\nu+p g^{\mu\nu}, \qquad u^\mu u_\mu=-1,$$

but closing the equations requires thermodynamics, conserved currents, composition, and a declaration of which variables equilibrate on the timescale under study. Magnetic fields add Maxwell stresses and characteristic modes. Radiation and neutrinos add distribution functions or moment hierarchies whose truncation needs a closure. Elasticity, viscosity, heat flow, superfluidity, and reactions introduce still more state variables and relaxation scales. Andersson and Comer 2021, §§ 5–6 and 8–10 develops the covariant single-fluid and multifluid variables, while Rezzolla and Zanotti 2013, chs. 3–4, pp. 133–257 gives a complementary perfect-fluid and characteristic treatment.

The central modeling question is therefore: which degrees of freedom remain resolved in this regime? A perfect fluid may be excellent in an opaque, collision-dominated core and wrong in a decoupling layer. Ideal magnetohydrodynamics may describe flux-frozen bulk motion but not a resistive current sheet. A cold beta-equilibrated equation of state may define an old neutron-star sequence but not a hot, lepton-rich proto-neutron star. The pages repeatedly derive the reduced equations from their ordering assumptions so that a closure is not used beyond its domain by accident.

Gravity converts this local input into a global boundary-value problem. For a static spherical star, the same energy density and pressure enter

$$\frac{dm}{dr}=4\pi r^2\varepsilon, \qquad \frac{dp}{dr} =-\frac{G(\varepsilon+p)(m+4\pi r^3p)} {r^2(1-2Gm/r)}.$$

The first equation counts the energy enclosed within radius $r$; the second tells us how rapidly pressure must fall outward to hold that matter in equilibrium. Compared with Newtonian hydrostatic balance, pressure contributes both to inertia through $\varepsilon+p$ and to gravity through $m+4\pi r^3p$. The factor $1-2Gm/r$ adds a geometric correction. When $p\ll\varepsilon$ and $Gm/r\ll1$, these corrections disappear and the familiar Newtonian equation returns. Pressure is therefore both part of the support and part of the gravitational problem.

Regular center data and a surface-matching condition select a stellar model. Varying the central state produces a mass–radius sequence; perturbing one model produces radial, fluid, spacetime, magnetic, or interface modes; placing two models in a binary produces a material tidal response. At every step, a local constitutive choice leaves a global signature, but that signature is filtered through the equilibrium family, boundary conditions, and observational forward model.

The thirteen-chapter path

The chapters are ordered so that each new physical system inherits a controlled material model rather than an unexplained formula.

1. Relativistic Fluids and Kinetic Theory derives perfect-fluid stress–energy, Euler flow, relativistic thermodynamics and sound cones, Boltzmann transport, Einstein–Vlasov matter, and entropy-admissible shocks.
2. Relativistic Magnetohydrodynamics and Plasmas adds electromagnetic stress, flux freezing, characteristic structure, force-free and resistive limits, reconnection, turbulence, and kinetic scale separation.
3. Radiation, Neutrinos, Elasticity, and Transport treats covariant transfer, moment closures, opacity, radiation hydrodynamics, neutrino transport, relativistic elasticity, and causal dissipative theories.
4. Dense Matter and Equations of State builds cold and finite-temperature thermodynamics, composition and beta equilibrium, nuclear and quark phases, tabulated-EOS validation, and uncertainty domains.
5. Relativistic Stellar Structure derives the Tolman–Oppenheimer–Volkoff equations, mass–radius families, baryon number and binding energy, radial stability, anisotropic or multifluid extensions, and exact-interior admissibility tests.
6. White Dwarfs and Neutron Stars compares degeneracy-supported families, neutron-star composition, crusts and superfluidity, maximum mass, thermal evolution, glitches, and magnetic evolution.
7. Rotating and Magnetized Stars develops stationary equilibria, mass-shedding and secular limits, numerical construction, magnetic deformation, relativistic multipoles, and stability boundaries.
8. Pulsars and Magnetospheres follows rotating fields through vacuum and force-free magnetospheres, pair-creation gaps, spin-down, timing tests, braking indices, magnetars, and timing irregularities.
9. Stellar Oscillations and Instabilities classifies radial and nonradial modes, spacetime and interface modes, the Chandrasekhar–Friedman–Schutz mechanism, r-mode damping and saturation, and compact-star asteroseismology.
10. Core Collapse and Supernovae connects progenitor structure, electron capture, trapping, bounce, shock formation, neutrino heating, multidimensional instability, proto-neutron-star evolution, and multimessenger signals.
11. Compact-Star Binaries and Tidal Dynamics develops Love numbers, dynamical tides, binary initial data, merger regimes, ejecta and kilonova inputs, postmerger spectra, and remnant thresholds.
12. Boson Stars and Exotic Compact Objects tests bosonic fields, self-interactions, formation and stability, thin-shell models, energy conditions, light rings, echoes, and proposed observational distinctions.
13. Dense-Matter Multimessenger Inference combines mass, X-ray pulse-profile, tidal, nuclear, and population information while retaining calibration, selection, prior, parameterization, and systematic uncertainty.

The sequence is not a claim that nature is modular. It is a way to expose where one approximation hands information to the next. A merger simulation, for example, simultaneously uses an EOS table, neutrino treatment, magnetic closure, relativistic initial data, and numerical discretization. Reading the contributing chapters separately makes their assumptions inspectable before they are coupled.

For sustained reference treatments, Oertel et al. 2017, §§ II–V surveys EOS construction across compact-star and supernova regimes; Paschalidis and Stergioulas 2017, §§ 2–5 develops stationary rotating-star structure and limits. The transport problems used in collapse and merger calculations are compared in Mezzacappa et al. 2020, §§ 4.7, 5.2–5.3 and 6 and Foucart 2023, §§ 3–4. These reviews organize the source landscape; the chapter pages give the derivations, model boundaries, and primary citations needed for each result.

Conventions that keep the chain coherent

The site uses metric signature $(-,+,+,+)$, sets $c=\hbar=k_{\mathrm B}=1$, and keeps $G$ explicit. Thus $GM/R$ is dimensionless compactness, whereas a physical mass $M$ alone is not a length. Restoring ordinary units replaces $GM/R$ by $GM/(Rc^2)$. The future-directed material velocity satisfies $u^\mu u_\mu=-1$, and $\varepsilon$ denotes total rest-frame energy density unless a page explicitly declares another quantity.

For spherical stars, $m(r)$ is the physical enclosed gravitational mass defined by the radial metric factor $1-2Gm(r)/r$. Baryonic mass, proper material energy, and gravitational mass are kept distinct. The circumferential stellar radius is $R$. For the static electric quadrupolar response of a nonrotating star in four-dimensional general relativity,

$$C=\frac{GM}{R}, \qquad \Lambda=\frac{2}{3}k_2 C^{-5},$$

after the Love number $k_2$, perturbing field, exterior solution, and matching convention have been fixed. Rotating-star pages reserve $\Omega$ for a background angular velocity and $\omega$ for a mode or Fourier frequency. With perturbations proportional to $e^{-i\omega t}$, $\operatorname{Im}\omega<0$ means damping and $\operatorname{Im}\omega>0$ means growth.

These declarations do more than standardize notation. They supply fast independent checks. A Newtonian limit must recover the usual hydrostatic, sound-wave, and degenerate-gas results. A zero-field limit must return hydrodynamics. Optically thick and thin limits must agree with diffusion and free streaming in their respective regimes. A static tide must emerge from the low-frequency dynamical response. A nonrotating limit must remove rotational multipoles and splittings. When an imported result uses different signs, units, frames, or Love-number normalization, the page translates it before comparison.

Recurring physical checks

Five tests recur across the volume.

Thermodynamic consistency. Pressure, entropy, chemical potentials, and energy density must derive compatibly from a thermodynamic potential. Maxwell relations, convexity, composition constraints, and phase constructions are checked. Interpolation of a table is part of the model because inaccurate derivatives can create spurious sound speeds or entropy production.

Conservation and balance. Stress–energy, baryon number, charge, lepton number, magnetic flux, angular momentum, and energy are tracked whenever the chosen model conserves them. Across a shock, pointwise differentiability is replaced by jump conditions and an entropy inequality. In transport, exchange terms must cancel between coupled sectors. In a transient, an energy budget must include binding, neutrino, electromagnetic, kinetic, thermal, and ejecta channels appropriate to the problem.

Characteristics and causality. The derivative $c_s^2=(\partial p/\partial\varepsilon)_{s,Y_i}$ at fixed entropy per baryon $s$ and composition $Y_i$ is a necessary sound-speed check for an ideal fluid. A multi-component or dissipative theory has additional characteristic modes and relaxation variables. Hyperbolicity, characteristic cones, degeneracies, and the domain of the constitutive approximation are examined together. An energy condition is not a substitute for this analysis.

Boundary and limiting behavior. Center regularity, stellar-surface matching, atmosphere treatment, kinetic mass-shell measure, transport boundary conditions, and binary asymptotics are part of the solution. Constant-density stars, polytropes, shock tubes, diffusion and free-streaming tests, slow rotation, and known exterior solutions provide benchmarks. A passed local residual cannot compensate for a wrong boundary condition.

Uncertainty and identifiability. A measured mass, radius-like posterior, tidal phase, spectral peak, or nuclear observable is not identical to an EOS parameter. The inference must carry its likelihood, interim priors, calibration, selection, population model, and nuisance parameters through the forward map. Results are reported over the density and model domain actually informed by the data. Pages concerning active observations or disputed mechanisms state an evidence date and a concrete reason to revisit the conclusion.

Reader routes

A first graduate route is Chapters 1, 4, 5, 6, 9, and 13: it moves from a relativistic fluid through dense thermodynamics and stellar equilibrium to modes and inference. Readers chiefly interested in closure should take Chapters 1–4 in order, then choose a physical application. The relation between continuum assumptions and characteristic speeds is especially important before using a numerical evolution.

For compact-star structure, start with Chapters 4 and 5, continue through white dwarfs and neutron stars in Chapter 6, and then add rotation and magnetism in Chapter 7. Chapter 8 turns those equilibria into pulsar observables. For dynamical transients, Chapters 3, 9, 10, and 11 connect transport and modes to collapse and merger matter. For horizonless alternatives, Chapters 5 and 9 provide the admissibility and stability tests used critically in Chapter 12.

Observers and nuclear theorists may enter through Chapter 13, but the backward links matter. A radius likelihood depends on the stellar sequence that maps an EOS to $(M,R)$; a tidal likelihood depends on the response convention and mass distribution; a nuclear prior occupies a stated density and composition regime. The inference pages therefore point back to the relevant forward problem rather than treating published one-dimensional intervals as independent measurements of the same scalar.

Boundaries with neighboring volumes

This volume treats material closure, compact-star equilibrium, stellar response, and the physical interpretation of dense matter. It uses the Einstein equation but does not rederive the foundations of general relativity. It develops continuum initial data and reproducible teaching benchmarks, while production evolution systems, discretizations, and simulation workflows belong to Formulations of General Relativity and Computational Relativity.

In collapse, the material physics of an iron core, bounce, neutrino transport, and proto-neutron-star evolution belongs here; finite-amplitude spacetime collapse, trapped regions, and horizon formation continue in Gravitational Collapse and Black-Hole Dynamics. For binaries, this volume supplies stellar Love response, modes, matter-rich merger regimes, and remnant properties. The insertion of response into orbital dynamics belongs to Relativistic Two-Body Dynamics, while radiative variables and waveform construction belong to Gravitational Radiation and Waveforms.

Detector calibration, search pipelines, and event-level likelihood production are not silently reproduced here. The multimessenger chapter develops only the likelihood structure needed to understand dense-matter inference and identifies the versioned products a reproduction would require. Likewise, black-hole accretion environments and population interpretation are treated elsewhere; Chapter 12 concentrates on the equilibrium, formation, stability, and observables of proposed horizonless material models.

Distinguishing these stages prevents double counting of both physical effects and observations. A continuum benchmark, waveform calculation, or likelihood analysis answers a different question; each treatment explains which additional results are needed to draw a broader conclusion.

Check the model-to-observable chain

A cold spherical model has $M=1.4\,M_\odot$, circumferential radius $R=12\,\mathrm{km}$, and quadrupolar Love number $k_2=0.080$. Using $GM_\odot=1.4766\,\mathrm{km}$ in units with $c=1$, compute $C$ and $\Lambda$. Then name three checks needed before interpreting this number as an observational EOS constraint.

Solution

The compactness is

$$C=\frac{1.4(1.4766\,\mathrm{km})}{12\,\mathrm{km}} \simeq0.1723.$$

Therefore

$$\Lambda=\frac{2}{3}(0.080)(0.1723)^{-5}\simeq3.5\times10^2.$$

The calculation first requires a thermodynamically stable, causal EOS in the densities sampled by this star and a regular TOV solution with a verified surface match. The perturbation calculation must use a declared $k_2$ normalization and independently reproduce the regular-center and exterior-matching limits. Finally, comparison with data requires a calibrated waveform and likelihood with its mass, spin, prior, frequency-band, and selection assumptions. Agreement of the scalar $\Lambda$ with a quoted interval alone would not validate that full chain.

References

• Andersson, Nils, and Gregory L. Comer. “Relativistic Fluid Dynamics: Physics for Many Different Scales.” Living Reviews in Relativity 24 (2021): 3. DOI.
• Foucart, Francois. “Neutrino Transport in General Relativistic Neutron Star Merger Simulations.” Living Reviews in Computational Astrophysics 9 (2023): 1. DOI.
• Mezzacappa, Anthony, et al. “Physical, Numerical, and Computational Challenges of Modeling Neutrino Transport in Core-Collapse Supernovae.” Living Reviews in Computational Astrophysics 6 (2020): 4. DOI.
• Oertel, Micaela, et al. “Equations of State for Supernovae and Compact Stars.” Reviews of Modern Physics 89 (2017): 015007. DOI.
• Paschalidis, Vasileios, and Nikolaos Stergioulas. “Rotating Stars in Relativity.” Living Reviews in Relativity 20 (2017): 7. DOI.
• Rezzolla, Luciano, and Olindo Zanotti. Relativistic Hydrodynamics. Oxford: Oxford University Press, 2013. DOI.