Neutron Star Properties

• Neutron star properties are defined by extreme densities, compact sizes, and strong relativistic effects, as described by the TOV equations.
• Their microphysics involve layered compositions with degenerate matter, exotic states, and transitions governed by dense-matter equations of state.
• Observations across radio, X-ray, and gravitational waves tightly constrain mass, radius, and tidal deformability, advancing dense-matter research.

A neutron star is a gravitationally bound remnant of massive stellar evolution, distinguished by extreme densities ($\gtrsim 10^{14}$–$10^{15}\, \mathrm{g/cm}^3$), compact radii ($\sim 10$–$15\ \mathrm{km}$), and the emergence of quantum, nuclear, and relativistic effects at macroscopic scale. Its astrophysical and microphysical properties are interrelated through general relativity, dense-matter equations of state (EoS), and multi-messenger observations, making neutron stars unique laboratories for supranuclear matter and strong-field gravity.

1. Equations of State and Macroscopic Structure

The equilibrium structure of a nonrotating, spherically symmetric neutron star is governed by the general relativistic Tolman–Oppenheimer–Volkoff (TOV) equations:

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

$$\frac{dP}{dr} = -\frac{G\,\bigl[\rho(r)+P(r)/c^2\bigr]\,\bigl[m(r)+4\pi r^3 P(r)/c^2\bigr]}{r\bigl[r-2G m(r)/c^2\bigr]}$$

where $r$ is the radial coordinate, $m(r)$ is the enclosed gravitational mass, $P(r)$ is the pressure, $\rho(r)$ is the energy density, $10^{15}\, \mathrm{g/cm}^3$0 is the gravitational constant, and $10^{15}\, \mathrm{g/cm}^3$1 is the speed of light. A physically meaningful stellar sequence is obtained by supplying a cold, β-equilibrated EoS $10^{15}\, \mathrm{g/cm}^3$2 and integrating outward from a central density $10^{15}\, \mathrm{g/cm}^3$3 until $10^{15}\, \mathrm{g/cm}^3$4 defines the stellar radius.

The solution space produces a one-parameter family of mass-radius $10^{15}\, \mathrm{g/cm}^3$5 configurations, typically with $10^{15}\, \mathrm{g/cm}^3$6 increasing with $10^{15}\, \mathrm{g/cm}^3$7 up to a critical maximum ($10^{15}\, \mathrm{g/cm}^3$8–$10^{15}\, \mathrm{g/cm}^3$9 for realistic nuclear EoS), beyond which $\sim 10$0 and stable equilibrium breaks down (Pizzochero, 2010). General relativistic corrections (e.g., pressure as a source of gravity) reduce $\sim 10$1 compared to Newtonian and Chandrasekhar–Stoner–Landau limits, enforcing an upper bound set by causality, $\sim 10$2 (Ekşi, 2015, Reisenegger et al., 2015).

Principal macroscopic properties:

Parameter	Typical Value	Notes
Mass, $\sim 10$3	$\sim 10$4–$\sim 10$5	Observed up to $\sim 10$6 (e.g., PSR J0740+6620)
Radius, $\sim 10$7	$\sim 10$8–$\sim 10$9 km	Canonical star: $15\ \mathrm{km}$0–$15\ \mathrm{km}$1 km
Central density, $15\ \mathrm{km}$2	$15\ \mathrm{km}$3–$15\ \mathrm{km}$4	$15\ \mathrm{km}$5 g cm$15\ \mathrm{km}$6
Surface gravitational redshift, $15\ \mathrm{km}$7	$15\ \mathrm{km}$8–$15\ \mathrm{km}$9	$\frac{dm}{dr} = 4\pi r^2 \rho(r)$0 (Pizzochero, 2010, Zhao, 2017)

2. Microphysical Composition and Equation of State

At the core of a neutron star, matter is arranged in distinct zones defined by density:

• Outer crust ($\frac{dm}{dr} = 4\pi r^2 \rho(r)$1 g cm$\frac{dm}{dr} = 4\pi r^2 \rho(r)$2): Fully ionized nuclei in a Coulomb lattice, degenerate electrons.
• Inner crust ($\frac{dm}{dr} = 4\pi r^2 \rho(r)$3 g cm$\frac{dm}{dr} = 4\pi r^2 \rho(r)$4): Neutron-rich nuclei with a free neutron superfluid and electrons.
• Outer core ($\frac{dm}{dr} = 4\pi r^2 \rho(r)$5): Homogeneous nuclear matter, mainly neutrons with a fraction of protons, electrons, and muons in beta equilibrium.
• Inner core (few$\frac{dm}{dr} = 4\pi r^2 \rho(r)$6): Potential presence of hyperons (e.g., $\frac{dm}{dr} = 4\pi r^2 \rho(r)$7, $\frac{dm}{dr} = 4\pi r^2 \rho(r)$8), meson condensates, deconfined quarks, depending on the EoS (Miyatsu et al., 2012, Zhao, 2017, Zhu et al., 27 Oct 2025).

Degenerate Fermi pressure is central:

$\frac{dm}{dr} = 4\pi r^2 \rho(r)$9

for ultra-relativistic particles at $$\frac{dP}{dr} = -\frac{G\,\bigl[\rho(r)+P(r)/c^2\bigr]\,\bigl[m(r)+4\pi r^3 P(r)/c^2\bigr]}{r\bigl[r-2G m(r)/c^2\bigr]}$$0.

Beta equilibrium among neutrons, protons, and electrons ($$\frac{dP}{dr} = -\frac{G\,\bigl[\rho(r)+P(r)/c^2\bigr]\,\bigl[m(r)+4\pi r^3 P(r)/c^2\bigr]}{r\bigl[r-2G m(r)/c^2\bigr]}$$1) and charge neutrality ($$\frac{dP}{dr} = -\frac{G\,\bigl[\rho(r)+P(r)/c^2\bigr]\,\bigl[m(r)+4\pi r^3 P(r)/c^2\bigr]}{r\bigl[r-2G m(r)/c^2\bigr]}$$2) shift the composition to increasingly neutron-rich matter as density increases (Pizzochero, 2010). Above the neutronization threshold ($$\frac{dP}{dr} = -\frac{G\,\bigl[\rho(r)+P(r)/c^2\bigr]\,\bigl[m(r)+4\pi r^3 P(r)/c^2\bigr]}{r\bigl[r-2G m(r)/c^2\bigr]}$$3 MeV, $$\frac{dP}{dr} = -\frac{G\,\bigl[\rho(r)+P(r)/c^2\bigr]\,\bigl[m(r)+4\pi r^3 P(r)/c^2\bigr]}{r\bigl[r-2G m(r)/c^2\bigr]}$$4 g/cm$$\frac{dP}{dr} = -\frac{G\,\bigl[\rho(r)+P(r)/c^2\bigr]\,\bigl[m(r)+4\pi r^3 P(r)/c^2\bigr]}{r\bigl[r-2G m(r)/c^2\bigr]}$$5) electron capture becomes energetically favorable.

At supranuclear densities, the EoS is determined by nucleonic interactions (as modeled by RBHF, chiral EFT, RMF, etc.), three-body forces, and possible exotic components (hyperons, quarks) (Wang et al., 2020, Logoteta et al., 2018, Sen et al., 2018, Miyatsu et al., 2012). Direct incorporation of hyperons or phase transitions softens the EoS and typically reduces the maximum supported mass—a manifestation of the "hyperon puzzle."

3. Superfluidity, Magnetism, and Transport Properties

Several superfluid and superconducting phenomena are realized in neutron stars:

• Superfluid neutrons in the inner crust and core (dominant $$\frac{dP}{dr} = -\frac{G\,\bigl[\rho(r)+P(r)/c^2\bigr]\,\bigl[m(r)+4\pi r^3 P(r)/c^2\bigr]}{r\bigl[r-2G m(r)/c^2\bigr]}$$6 and $$\frac{dP}{dr} = -\frac{G\,\bigl[\rho(r)+P(r)/c^2\bigr]\,\bigl[m(r)+4\pi r^3 P(r)/c^2\bigr]}{r\bigl[r-2G m(r)/c^2\bigr]}$$7 pairing), characterized by energy gaps $$\frac{dP}{dr} = -\frac{G\,\bigl[\rho(r)+P(r)/c^2\bigr]\,\bigl[m(r)+4\pi r^3 P(r)/c^2\bigr]}{r\bigl[r-2G m(r)/c^2\bigr]}$$8–$$\frac{dP}{dr} = -\frac{G\,\bigl[\rho(r)+P(r)/c^2\bigr]\,\bigl[m(r)+4\pi r^3 P(r)/c^2\bigr]}{r\bigl[r-2G m(r)/c^2\bigr]}$$9 MeV. Proton superconductivity ($r$0) emerges at lower densities in the core.
• Observable consequences include suppression of the specific heat $r$1 and of neutrino emissivity $r$2 below the critical temperature, directly influencing the cooling trajectory and the shape of surface temperature curves.
• Rotational glitches arise from the differential rotation of superfluid neutron vortices pinned to the crustal lattice, facilitating sudden angular momentum transfer during vortex unpinning (Pizzochero, 2010).

The external dipolar magnetic field spans $r$3–$r$4 G, with stronger fields ($r$5 G) in magnetars. The magnetic field structure modifies the charged-particle energy spectra (Landau quantization) and introduces pressure anisotropy (Malik et al., 2018, Menezes, 2021), affecting mass-radius curves by a few percent at $r$6 G.

Transport and emission properties:

• Thermal evolution is governed by neutrino emission (modified and direct Urca processes) and photon cooling, with early neutrino-dominated stage ($r$7 yr) rapidly reducing core temperatures to $r$8 K.
• High-energy emission arises from surface thermal emission (blackbody, $r$9–$m(r)$0 keV), non-thermal magnetospheric processes (radio, X-ray, gamma-ray), and field decay in magnetars.

4. Observational Constraints and Global Properties

Neutron star macroscopic properties are tightly constrained by multi-messenger data:

• Masses via radio timing and Shapiro delay in binary pulsars (e.g., PSR J0348+0432, $m(r)$1).
• Radii from X-ray pulse profiling (NICER) and thermal emission fits ($m(r)$2–$m(r)$3 km).
• Tidal deformabilities ($m(r)$4) from gravitational-wave measurement of binary mergers (GW170817): at $m(r)$5, current bounds $m(r)$6 (90% C.I.) (Kumar et al., 2019, Golomb et al., 2024), restricting the EoS stiffness at $m(r)$7 (Collazos, 2023, Wang et al., 2020, Zhang et al., 2019).
• Moment of inertia ($m(r)$8), gravitational redshift ($m(r)$9), and compactness ($P(r)$0), which are inferred via pulsar timing and modeling of high-mass neutron stars.
• Maximum mass ($P(r)$1) inferred jointly from heavy pulsar observations and GW mergers: $P(r)$2 (Golomb et al., 2024), consistent with the $P(r)$3 lower bound.

Canonical neutron star properties (as constrained by current theory and observation):

Property	Range/Value	Reference Cases
$P(r)$4	$P(r)$5–$P(r)$6	PSR J0740+6620, GW170817
$P(r)$7	$P(r)$8–$P(r)$9 km	NICER, GW170817
$\rho(r)$0	$\rho(r)$1	GW170817

Precise correlations between nuclear matter parameters and neutron-star observables have been quantified:

$\rho(r)$2

with $\rho(r)$3, $\rho(r)$4, $\rho(r)$5. $\rho(r)$6 and $\rho(r)$7 thus tightly constrain the EoS parameter space (Providência, 2019).

5. Universal Relations and Population Systematics

Dimensionless relations—universal across a wide variety of microphysical EoS—link the moment of inertia ($\rho(r)$8), quadrupole moment ($\rho(r)$9), and tidal deformability ($10^{15}\, \mathrm{g/cm}^3$00):

$10^{15}\, \mathrm{g/cm}^3$01

for $10^{15}\, \mathrm{g/cm}^3$02 being $10^{15}\, \mathrm{g/cm}^3$03, $10^{15}\, \mathrm{g/cm}^3$04, or $10^{15}\, \mathrm{g/cm}^3$05, with fit coefficients universal to within $10^{15}\, \mathrm{g/cm}^3$06 (Kunz, 2022, Kumar et al., 2019).

These I-Love-Q and binary-Love relations enable the inference of properties (e.g., $10^{15}\, \mathrm{g/cm}^3$07, $10^{15}\, \mathrm{g/cm}^3$08, $10^{15}\, \mathrm{g/cm}^3$09) from measured $10^{15}\, \mathrm{g/cm}^3$10, and facilitate joint GW/EM constraints on the EoS without detailed microphysical model dependence (Kumar et al., 2019).

Population synthesis and hierarchical Bayesian analyses now distinguish between Galactic-EM (radio/X-ray) and GW-inferred NS populations, finding consistent maximum masses ($10^{15}\, \mathrm{g/cm}^3$11) but statistically distinct mass distributions. The radius inference ($10^{15}\, \mathrm{g/cm}^3$12 km, 90% HPD) is still primarily determined by microphysics, while details of the mass distribution reflect formation and astrophysical selection effects (Golomb et al., 2024).

6. Bulk, Surface, and Symmetry Properties

Bulk and surface nuclear properties of neutron star matter, such as incompressibility ($10^{15}\, \mathrm{g/cm}^3$13), symmetry energy ($10^{15}\, \mathrm{g/cm}^3$14), slope parameter ($10^{15}\, \mathrm{g/cm}^3$15), and curvature ($10^{15}\, \mathrm{g/cm}^3$16), vary systematically with EoS stiffness and stellar mass (Kumar et al., 2021). Surface-averaged values are calculated via local density approximation and the coherent density fluctuation model:

$10^{15}\, \mathrm{g/cm}^3$17

For soft EoS (G3), $10^{15}\, \mathrm{g/cm}^3$18 km, $10^{15}\, \mathrm{g/cm}^3$19 MeV, $10^{15}\, \mathrm{g/cm}^3$20 MeV for $10^{15}\, \mathrm{g/cm}^3$21, rising with stiffer EoS (NL3) and higher mass. Observed trends:

Model	$10^{15}\, \mathrm{g/cm}^3$22 (MeV)	$10^{15}\, \mathrm{g/cm}^3$23 (MeV)	$10^{15}\, \mathrm{g/cm}^3$24 (MeV)	$10^{15}\, \mathrm{g/cm}^3$25 (MeV)
NL3, $10^{15}\, \mathrm{g/cm}^3$26	$10^{15}\, \mathrm{g/cm}^3$27	$10^{15}\, \mathrm{g/cm}^3$28	$10^{15}\, \mathrm{g/cm}^3$29	$10^{15}\, \mathrm{g/cm}^3$30
G3, $10^{15}\, \mathrm{g/cm}^3$31	$10^{15}\, \mathrm{g/cm}^3$32	$10^{15}\, \mathrm{g/cm}^3$33	$10^{15}\, \mathrm{g/cm}^3$34	$10^{15}\, \mathrm{g/cm}^3$35
IU-FSU, $10^{15}\, \mathrm{g/cm}^3$36	$10^{15}\, \mathrm{g/cm}^3$37	$10^{15}\, \mathrm{g/cm}^3$38	$10^{15}\, \mathrm{g/cm}^3$39	$10^{15}\, \mathrm{g/cm}^3$40

High $10^{15}\, \mathrm{g/cm}^3$41 correlates with larger radii and higher central proton fractions, influencing cooling and Urca threshold conditions (Providência, 2019, Kumar et al., 2021).

7. Open Problems and Future Prospects

The maximum mass and radii of neutron stars continue to serve as critical discriminants among candidate EoSs, especially regarding the transition to exotic phases (hyperonic matter, quark deconfinement). Recent advances employ Gaussian process and spectral approaches to reconstruct the EoS with minimal model bias (Golomb et al., 2024). Model-agnostic multi-messenger inference is expected to constrain $10^{15}\, \mathrm{g/cm}^3$42 to within $10^{15}\, \mathrm{g/cm}^3$43 km and $10^{15}\, \mathrm{g/cm}^3$44 within $10^{15}\, \mathrm{g/cm}^3$45 as GW detector catalogs grow.

Key uncertainties persist in the microphysics at $10^{15}\, \mathrm{g/cm}^3$46–$10^{15}\, \mathrm{g/cm}^3$47—notably, the high-density symmetry energy, hyperonic and quark matter interactions, and the character of the hadron–quark transition (Zhu et al., 27 Oct 2025). Global properties (mass, radius, tidal response) remain largely insensitive to low-energy nuclear data once EoS models fit symmetry energy and its slope at saturation, highlighting the need for direct constraints at supranuclear densities (Carlson et al., 2022, Golomb et al., 2024).

Astrophysical observations—precision pulse profile modeling, moment of inertia measurements in double pulsars, and next-generation GW events—along with theoretical advances in dense-matter modeling will continue to refine our understanding of the most extreme nuclear-physical environments accessible to experiment.