Fact-checked by Grok 2 weeks ago

Neutron star

A neutron star is a stellar remnant formed from the core collapse of a massive star—typically one with an initial mass between 8 and 20 times that of the Sun—following a supernova explosion, resulting in an extremely dense object composed primarily of neutrons.^[1] These stars pack roughly 1.4 times the Sun's mass into a sphere only about 12 kilometers (7.5 miles) in diameter, making them among the densest objects observable in the universe.^[2] Their surface gravity is immense, approximately 10¹¹ times stronger than Earth's, and a teaspoon of neutron star material would weigh billions of tons on our planet.^[3] Neutron stars form when a progenitor star exhausts its nuclear fuel, causing the core to exceed the Chandrasekhar limit of about 1.44 solar masses and collapse under gravity until neutron degeneracy pressure halts the implosion.^[3] This process ejects the star's outer layers in a type II supernova, leaving behind the ultra-compact core where protons and electrons fuse into neutrons, dominating the composition.^[2] The resulting object has a typical mass range of 1.0 to 2.0 solar masses, with radii between 10 and 14 kilometers, though theoretical models suggest maximum masses up to 2.5 solar masses before further collapse into a black hole.^[3] Formation occurs in the final stages of stellar evolution for stars too massive for white dwarfs but insufficient for direct black hole formation without significant mass loss.^[1] The internal structure of a neutron star features a thin outer crust of neutron-rich nuclei, a denser inner crust with a superfluid neutron component, and a core of degenerate neutron matter possibly including exotic particles like hyperons or quarks at densities exceeding 10¹⁵ grams per cubic centimeter—several times nuclear saturation density.^[3] Many neutron stars rotate rapidly, with periods from milliseconds to seconds, and possess magnetic fields ranging from 10⁸ to 10¹⁵ gauss, which can accelerate particles and produce observable emissions.^[2] They cool from initial temperatures of 10¹⁰ to 10¹¹ Kelvin primarily through neutrino emission in the first million years, followed by photon radiation from the surface.^[3] Prominent subtypes include pulsars, rapidly spinning neutron stars that emit beamed radiation detectable as periodic pulses when aligned with Earth, first discovered in 1967 and numbering over 3,700 known examples as of 2025.^[4] Magnetars represent an extreme variant with magnetic fields up to 10¹⁵ gauss, capable of powering giant flares that release more energy in seconds than the Sun emits over millennia.^[1] Neutron stars serve as unique laboratories for studying extreme physics, including general relativity, quantum chromodynamics, and the equation of state of supranuclear matter, with observations from radio telescopes, X-ray satellites, and gravitational wave detectors providing constraints on their properties.^[3]

Formation

Core-collapse supernovae

Neutron stars primarily form through the core-collapse of massive stars with initial masses ranging from approximately 8 to 20 solar masses. These stars exhaust their nuclear fuel through successive stages of hydrogen, helium, carbon, oxygen, and silicon burning, culminating in the formation of an iron-nickel core with a mass of 1.2 to 2 solar masses. Iron cannot fuse exothermically, so the core, initially supported by electron degeneracy pressure, becomes unstable when electron captures on iron-group nuclei reduce the electron number and pressure, dropping the effective Chandrasekhar limit below the core mass. This triggers a rapid implosion, with the core collapsing in milliseconds to nuclear densities of about 10^{14} to 10^{15} g/cm³, where repulsive nuclear forces halt the infall, causing a bounce that forms a proto-neutron star with an initial radius of ~20-40 km, which contracts to 10-15 km during subsequent evolution. During the collapse, electron capture drives neutronization of the core, converting protons and electrons into neutrons and emitting electron neutrinos, which escape and carry away significant entropy and energy, accelerating the infall to relativistic speeds. Neutrino emission plays a pivotal role in the explosion energetics: post-bounce, a shock wave forms but stalls, and subsequent neutrino diffusion from the hot proto-neutron star deposits energy via absorption in the overlying material, heating it in a "gain" region and potentially reviving the shock through convection-enhanced transport. This neutrino-driven mechanism releases a total of about 3 × 10^{53} erg in neutrinos across all flavors, with only 1-2% absorbed to impart the observed supernova kinetic energy of 10^{51} erg, while the rest escapes, enabling the core's neutronization to proceed efficiently.^[5] Immediately after bounce, the proto-neutron star is hot (temperatures ~10-50 MeV), lepton-rich (electron fraction Ye ~0.3-0.4), and opaque to neutrinos due to high densities. Over the subsequent deleptonization phase, lasting 10 to 20 seconds, trapped neutrinos gradually diffuse outward, cooling the interior and reducing the lepton abundance as the star contracts and becomes more neutron-rich. This phase sets the initial conditions for the neutron star's evolution, with the proto-neutron star contracting to a stable configuration while ejecting the star's envelope in the supernova blast.^[6] Observationally, neutron stars are tightly linked to Type II, Ib, and Ic supernovae, which mark core-collapse events in massive stars retaining or stripping their hydrogen/helium envelopes. Prominent examples include the Crab pulsar (PSR B0531+21) at the heart of the Crab Nebula remnant from a Type II supernova in 1054 CE, and the young neutron star candidate in Cassiopeia A, the remnant of a Type IIb event circa 1680 CE. These associations confirm that core-collapse supernovae produce the vast majority of Galactic neutron stars, with remnants often showing X-ray and radio signatures of the compact object.^[7] Recent three-dimensional general relativistic simulations of core-collapse from low-mass progenitors (9-10 solar masses) have elucidated the formation of the lightest neutron stars through fallback accretion. In these models, weak explosions allow some ejecta to reverse course and accrete back onto the proto-neutron star, reducing its final mass to as low as 1.19 solar masses. Recent 2025 studies, including supernova simulations, confirm that the 1.174 solar mass companion in PSR J0453+1559 is a neutron star formed via core-collapse with fallback accretion. Such simulations, incorporating multi-group neutrino transport, highlight how progenitor structure and explosion asymmetry influence outcomes, challenging earlier limits from electron-capture supernovae and extending the lower mass boundary for core-collapse remnants.

Alternative formation channels

While the majority of neutron stars form through the core collapse of massive stars, alternative pathways involve the collapse of white dwarfs or the electron capture in degenerate cores of lower-mass progenitors. These channels typically produce neutron stars with distinct properties, such as lower masses or reduced natal kick velocities, and are thought to contribute a small fraction to the overall neutron star population.^[8] Accretion-induced collapse occurs when a low-mass oxygen-neon-magnesium white dwarf in a binary system accretes sufficient material from its companion to exceed the Chandrasekhar limit, triggering gravitational collapse into a neutron star rather than a Type Ia supernova. This process is particularly relevant for white dwarfs with initial masses around 1.2–1.4 solar masses, potentially yielding neutron stars with masses as low as 0.7–1.0 solar masses. Recent simulations indicate that such collapses can occur within planetary nebulae, preserving evidence of the progenitor system's evolution.^[9]^[10] Merger-induced collapse of double white dwarf binaries represents another pathway, where the coalescence of two oxygen-neon white dwarfs—each with masses exceeding 1.3 solar masses—leads directly to neutron star formation without an accompanying bright supernova. This mechanism often results in rapidly rotating neutron stars due to the angular momentum from the merger, and it may explain isolated neutron stars or those in unusual binary configurations. Population synthesis models suggest these events produce single neutron stars that evade detection as luminous transients.^[11]^[12] Electron capture supernovae arise from the collapse of oxygen-neon-magnesium cores in intermediate-mass stars with initial masses of approximately 7–10 solar masses, where degenerate electron capture on neon and magnesium nuclei destabilizes the core, leading to an explosion and neutron star remnant. These events are fainter than standard core-collapse supernovae, with energies around 10^50 ergs, and produce neutron stars with masses typically below 1.3 solar masses. Unlike iron-core collapses, the more symmetric explosion in electron capture events imparts lower natal kick velocities, often below 100 km/s, facilitating the retention of wide binaries.^[13]^[14] Theoretical estimates from binary population synthesis indicate that alternative formation channels account for less than 1% of all neutron stars, primarily influencing the low-mass end of the mass distribution and contributing to older population age profiles in globular clusters or the galactic halo. This rarity stems from the specific evolutionary requirements, such as close binary interactions or precise core compositions.^[15] Recent hydrodynamic simulations from 2024 have linked these channels to underluminous supernovae, with accretion-induced and merger-induced collapses producing dim, neutrino-driven outflows and peculiar kick velocities that differ from canonical values, potentially observable in gamma-ray bursts or fast radio bursts from low-mass remnants. These models highlight how such events could explain outliers in neutron star demographics without invoking standard core-collapse mechanisms.^[16]^[10]

Physical Properties

Mass and radius

Neutron stars possess masses generally spanning 1.1 to 2.0 solar masses (M_\odot), with the majority clustered around 1.4 M_\odot, reflecting the typical outcome of core-collapse supernovae from progenitors of 8–20 M_\odot. ^[17] The highest precisely measured masses approach 2.08 M_\odot, as observed in the pulsar PSR J0740+6620, while theoretical models informed by nuclear physics and observations suggest a maximum stable mass between 2.1 and 2.5 M_\odot before collapse to a black hole. Radii for these objects are remarkably compact, typically 10–14 km, yielding average densities exceeding $10^{17} kg/m³ and highlighting their extreme compactness. ^[18] Masses are primarily determined through pulsar timing observations in binary systems, where relativistic effects like the Shapiro delay—caused by the pulsar's signal passing through the companion's gravitational field—allow precise inference of both components' masses; for instance, the double neutron star binary PSR J0737−3039A/B yields masses of approximately 1.338 M_\odot and 1.250 M_\odot. ^[17] Radii measurements rely on X-ray observations, including pulse profile modeling with the Neutron Star Interior Composition Explorer (NICER) telescope, which analyzes hotspots on the stellar surface to constrain geometry, and cooling tails of thermonuclear X-ray bursts from accreting neutron stars, where the burst luminosity traces the emitting area. ^[19] NICER's analyses of PSR J0030+0451 and PSR J0740+6620 provide some of the tightest constraints, with radii of 12.71 ± 1.14 km and 12.39 +1.30 -0.98 km, respectively, assuming masses near 1.4 M_\odot and 1.4–2.0 M_\odot. ^[19] The mass-radius (M-R) relation for neutron stars arises from integrating the equation of state (EOS) of dense matter under general relativity, producing a characteristic curve where radius decreases with increasing mass up to a maximum, beyond which no stable configurations exist; nuclear physics models predict a nearly universal low-mass branch, independent of specific EOS details at subnuclear densities. ^[20] Recent NICER data from 2023–2025, combined with Bayesian inference across multiple sources, have narrowed the radius for a 1.4 M_\odot neutron star to approximately 12.3 km (with 90% confidence bounds of 11.0–13.7 km), excluding softer EOS models that would yield radii below 11 km. ^[21] This relation's shape depends sensitively on the high-density EOS, with stiffer equations supporting larger maximum masses and radii. Stability against gravitational collapse imposes a Chandrasekhar-like limit via solutions to the Tolman–Oppenheimer–Volkoff (TOV) equation, capping stable neutron star masses at around 2–3 M_\odot depending on the EOS stiffness, beyond which pressure fails to counter self-gravity. ^[22] Multimessenger observations, including the gravitational wave event GW170817—a binary neutron star merger—have further refined these bounds; analyses of the inspiral tidal deformability and post-merger remnants suggest radii greater than 11 km for typical 1.4 M_\odot stars, ruling out overly compact models and supporting a moderately stiff EOS at nuclear densities. ^[23]

Object	Mass (M_\odot)	Radius (km)	Method	Reference
PSR J0030+0451	~1.4	12.71 ± 1.14	NICER pulse profiles	^[19]
PSR J0740+6620	1.4–2.0	12.39 +1.30 -0.98	NICER pulse profiles	^[19]
PSR J0437−4715	~1.4	13.0 ± 1.0	NICER & X-ray bursts	^[24]
GW170817 components	~1.17–1.60 each	>11 (for 1.4 M_\odot)	Gravitational waves	^[23]

Density and composition

Neutron stars exhibit extraordinary densities, with average values of approximately $10^{14} g/cm³ and central densities up to $10^{15} g/cm³, far exceeding that of atomic nuclei. At the core, densities reach nuclear saturation density, around $2.8 \times 10^{14} g/cm³, where matter is compressed to extremes that challenge our understanding of quantum chromodynamics. The outer crust, extending from the surface to densities of about $10^{11} g/cm³, consists of a lattice of neutron-rich iron-like nuclei immersed in a degenerate electron gas. Transitioning inward, the inner crust spans densities from roughly $4 \times 10^{11} g/cm³ to $10^{14} g/cm³, featuring a lattice of increasingly neutron-drip nuclei accompanied by a sea of free superfluid neutrons and degenerate electrons. This superfluid component arises from neutron pairing, enabling zero-viscosity flow at low temperatures. In the core, which occupies the central region at densities exceeding $10^{14} g/cm³, the dominant constituent is a superfluid of degenerate neutrons, with contributions from protons and electrons to maintain beta equilibrium. At higher densities, beyond about twice nuclear saturation (\sim 5.6 \times 10^{14} g/cm³), exotic phases may emerge, including hyperons such as lambda or sigma particles, which soften the matter due to their appearance via weak interactions. Even denser regimes, potentially up to several times nuclear saturation, could involve deconfined quark matter or condensates of pions and kaons, where hadronic degrees of freedom dissolve into quark-gluon plasma or bosonic fields.^[25] First-order phase transitions to these exotic states can cause abrupt softening of the equation of state, reducing pressure at fixed density and potentially leading to "twin stars"—compact objects of identical mass but differing radii and internal compositions, one purely hadronic and the other hybrid.^[25] Such transitions are modeled in hybrid equations of state, often occurring at critical densities of 2–4 times nuclear saturation.^[25] Multimessenger observations, including gravitational waves from binary mergers like GW170817 and X-ray measurements from NICER of pulsar radii, impose stringent limits that disfavor strongly softening phase transitions to exotic matter in many models.^[26] For instance, the tidal deformability and radius constraints for 1.4 M_\odot stars (11–14 km) rule out equations of state with early onset of hyperonic or quark phases that would predict softer interiors, though some hybrid models remain viable if the quark phase is sufficiently stiff.^[25]^[26] These data favor neutron-dominated cores up to at least twice nuclear density, with exotic contributions possible only in the innermost regions.^[26]

Magnetic fields

Neutron stars exhibit magnetic fields spanning a wide range of strengths, typically from $10^8 to $10^{15} Gauss (G), which play a crucial role in their observability and emission properties. Ordinary pulsars possess fields around $10^{11} to $10^{13} G, while recycled millisecond pulsars have weaker fields of $10^8 to $10^{10} G; magnetars, a subclass with exceptionally strong fields of $10^{14} to $10^{15} G, represent the upper extreme.^[27] The origins of these magnetic fields are attributed to several mechanisms rooted in the star's formation process. Fossil fields from the progenitor massive star can be amplified by factors up to $5 \times 10^{11} during the core-collapse supernova, transforming progenitor fields of ~10 kG into neutron star fields reaching $5 \times 10^{15} G. Additionally, dynamo action driven by convection in the proto-neutron star stage can generate poloidal fields up to $5 \times 10^{16} G, accompanied by toroidal components 100–300 times stronger. Compression of the progenitor's field during gravitational collapse provides another pathway, enhancing the field through flux conservation.^[27] Over time, neutron star magnetic fields evolve and decay primarily through ohmic dissipation, where magnetic energy is dissipated via finite electrical conductivity in the stellar crust, and Hall drift, a nonlinear process that cascades magnetic energy to smaller scales. Ohmic decay operates on timescales of $10^3 to $10^6 years, influenced by crustal conductivity from phonons and impurities, while Hall drift proceeds on shorter scales of $10^3 to $10^5 years, potentially leading to instabilities. These processes result in overall field decay over $10^4 to $10^6 years, shaping the long-term spin-down of the star.^[27]^[28] Observational constraints on magnetic fields come from X-ray spectroscopy, where cyclotron resonance lines—arising from electron or proton cyclotron absorption—reveal field strengths, such as ~$10^{12} G in accreting pulsars or $10^{14}–$10^{15} G in magnetars. Pulsar spin-down measurements yield braking indices n = 3 - \frac{\dot{P} P}{\ddot{P}}, where deviations from the canonical n=3 (for pure dipole radiation) indicate field evolution, as seen in sources like PSR J1734-3333 with n \approx 0.9–3. In 2025, XRISM observations detected unexpected cosmic winds from an accreting neutron star in the binary system GX 13+1, such as sluggish outflows at ~1 million km/h, challenging models of wind launching and providing new insights into field-wind interactions.^[27]^[29]

Temperature and cooling

Newly formed neutron stars, known as proto-neutron stars, emerge from core-collapse supernovae with extremely high internal temperatures on the order of $10^{11} K.^[30] This intense heat arises from the gravitational binding energy released during collapse, with the star initially trapped in a dense, lepton-rich state.^[30] Over the first few seconds to minutes, the proto-neutron star undergoes deleptonization, a rapid phase where trapped neutrinos diffuse outward, reducing the lepton fraction and releasing significant thermal energy.^[30] This process lasts approximately the first minute, marking the transition from the birth state to a more transparent configuration for further cooling.^[30] Following deleptonization, the neutron star enters the neutrino cooling phase, which dominates thermal evolution for roughly the first $10^3 to $10^5 years.^[30] During this period, neutrino emission from the dense core provides the primary cooling mechanism, far exceeding photon losses from the surface.^[30] The modified Urca process, involving neutron-proton scattering and weak interactions, becomes the dominant neutrino emission channel after about one year in standard non-superfluid models without direct Urca processes.^[30] This two-body reaction, enhanced by the involvement of a bystander nucleon, efficiently carries away energy as neutrino-antineutrino pairs, leading to a core temperature drop to around $10^9 K within the first decade.^[30] By approximately $10^3 years, the effective surface temperature has cooled to about $10^6 K through these neutrino-dominated losses.^[30] After the neutrino cooling phase, the star transitions to the photon-dominated era, where surface emission of photons becomes the leading heat loss mechanism, typically beyond $10^5 years.^[30] At this stage, the core temperature stabilizes at around $10^8 K or lower, and the surface radiates primarily as a blackbody in soft X-rays, with luminosity scaling as T_s^4 where T_s is the surface temperature.^[30] Observational constraints on cooling come from fitting X-ray spectra of isolated neutron stars to atmospheric or blackbody models, revealing surface temperatures that correlate with spin-down ages; for instance, young pulsars like the Crab show T_s \approx 2 \times 10^6 K at age \sim 10^3 years, aligning with modified Urca predictions.^[31] These fits, often using data from Chandra and XMM-Newton, confirm the expected T_s \propto t^{-1/2} decline in the photon era for standard cooling curves.^[31] Recent observations with the James Webb Space Telescope (JWST) have provided evidence for a compact object, likely a neutron star, at the heart of the Supernova 1987A remnant through detection of emission lines from ionized gas, based on 2022 data analyzed as of 2024 for the \sim 38-year-old object.^[32] This detection refines theoretical cooling curves by constraining early neutrino emissivities and envelope properties in massive progenitors, offering insights into the transition from neutrino to photon dominance.^[32]

Internal Structure

Layered composition

Neutron stars exhibit a radially stratified internal structure, determined primarily by the equation of state of dense matter and gravitational equilibrium, extending from a tenuous outer atmosphere to a dense central core.^[33] The outermost layer is the atmosphere, a thin envelope composed predominantly of hydrogen and helium accreted from the interstellar medium or companion stars, with a typical thickness of approximately 10 cm. This layer has extremely low density, on the order of 10^5 g/cm³ at the surface, and serves as the site where thermal radiation escapes, strongly influenced by the star's intense magnetic fields. Beneath the atmosphere lies the crust, divided into outer and inner regions, with a total thickness of about 1 km for a typical neutron star of radius 10–12 km. The outer crust, extending from densities of roughly 10^4 g/cm³ up to the neutron drip point at approximately 10^11 g/cm³, consists of a lattice of neutron-poor nuclei such as iron-56, immersed in a sea of degenerate electrons. At the neutron drip density, free neutrons begin to appear as the binding energy of nuclei allows neutrons to drip out, marking the transition to the inner crust. The inner crust spans densities from about 10^11 g/cm³ to around 10^14 g/cm³ and features a lattice of neutron-rich nuclear clusters surrounded by a gas of free neutrons that become superfluid, along with degenerate electrons and a smaller proton component. At the higher densities within this layer, the nuclear matter adopts exotic non-spherical configurations known as nuclear pasta phases, including rod-like (spaghetti) and slab-like (lasagna) structures, which influence the crust's mechanical properties. The crust-core boundary occurs at a density of roughly 10^14 g/cm³, where the nuclear clusters dissolve. The core constitutes the bulk of the neutron star, with a radial extent of approximately 10 km, and is subdivided into outer and inner regions based on density and composition.^[33] The outer core, from about 10^14 g/cm³ up to several times the nuclear saturation density (∼2.8 × 10^14 g/cm³), is a uniform fluid of mostly neutrons with a minority of protons and electrons (npe matter), where baryons are no longer confined to nuclei.^[33] Deeper in the inner core, at densities exceeding twice the nuclear saturation density (above ∼5.6 × 10^14 g/cm³), phase transitions to exotic forms of matter—such as hyperonic matter, pion or kaon condensates, or deconfined quark matter—may occur, though their exact nature remains uncertain and depends on the underlying equation of state.^[33] These layers collectively support the enormous pressure required to counteract gravitational collapse in the compact object.^[33]

Equation of state

The equation of state (EOS) for neutron star matter describes the thermodynamic relation between pressure P and baryon density \rho, denoted P(\rho), which determines the structure and stability of these compact objects. This relation arises from the strong nuclear interactions at densities ranging from nuclear saturation (\rho_0 \approx 0.16 fm^{-3}) to several times \rho_0 in the core, ultimately rooted in quantum chromodynamics (QCD) for the asymptotic high-density regime where quark degrees of freedom dominate, while at lower densities it is modeled using effective nuclear interactions constrained by atomic nuclei data.^[34] Theoretical models for the EOS vary in their treatment of matter composition and interactions. The relativistic mean field (RMF) model approximates nucleons as point particles exchanging isoscalar and isovector mesons, providing a semi-microscopic description that reproduces nuclear saturation properties and supports neutron star masses up to about 2 solar masses (M_\odot).^[35] The quark-meson coupling (QMC) model extends this by treating nucleons as composite quark bags coupled to meson fields, effectively incorporating quark structure while maintaining hadronic phenomenology at moderate densities.^[36] Hybrid EOS models combine a hadronic phase described by RMF or similar at lower densities with a quark-matter core at higher densities, often featuring a first-order phase transition where deconfinement occurs, leading to potentially softer EOS transitions that affect star compactness.^[25] Uncertainties in the EOS are pronounced at supranuclear densities, particularly the stiffness—characterized by the adiabatic index \Gamma = (d \ln P / d \ln \rho) + 1—which governs how rapidly pressure rises with density and directly impacts the maximum stable neutron star mass, with stiffer EOS allowing masses exceeding $2 M_\odot while softer ones predict lower limits.^[25] These ambiguities stem from incomplete knowledge of hyperon emergence, pion condensation, and quark deconfinement, complicating predictions for core compositions. Laboratory constraints help mitigate this: heavy-ion collisions at facilities like RHIC and LHC recreate hot, dense conditions to probe the speed of sound c_s^2 = dP/d\epsilon (where \epsilon is energy density), revealing that c_s approaches but does not exceed the causal limit c at extreme densities, while neutron skin thickness measurements in heavy nuclei (e.g., ^{208}Pb via the PREX-II experiment) quantify the nuclear symmetry energy's slope L, linking it to the isospin asymmetry in neutron-rich matter and favoring stiffer low-density EOS.^[37] The I-Love-Q relations offer EOS-independent insights, establishing near-universal correlations among the dimensionless tidal deformability \Lambda, moment of inertia \bar{I}, and quadrupole tidal moment \bar{Q} for slowly rotating neutron stars, with fitting functions like \bar{C}_\Lambda = a_1 (\bar{I})^{a_2} (where a_1, a_2 are fitted constants) holding to within 1% accuracy across diverse EOS, enabling indirect probes of interior physics via gravitational wave measurements without relying on specific P(\rho).^[38] In 2025, advanced numerical simulations incorporating gravitational wave data from binary neutron star mergers, such as GW170817, have assessed hybrid EOS viability, finding that phase-transition signatures in post-merger remnants—manifesting as distinct oscillation modes or ejecta properties—remain consistent with observations but require next-generation detectors like Einstein Telescope to confirm quark-core presence in massive stars exceeding $2 M_\odot.^[39]

Tolman–Oppenheimer–Volkoff equation

The Tolman–Oppenheimer–Volkoff (TOV) equation provides the relativistic generalization of hydrostatic equilibrium for spherically symmetric, static configurations of matter under strong gravity, such as neutron stars. Derived within the framework of general relativity, it describes the balance between gravitational collapse and internal pressure support, incorporating corrections from spacetime curvature that become crucial at nuclear densities. This equation, along with the continuity equation for mass, forms the core mathematical framework for modeling neutron star interiors. The TOV equation originates from the Oppenheimer–Volkoff analysis in 1939, which extended the Newtonian hydrostatic equilibrium equation \frac{dP}{dr} = -\frac{G M(r) \rho(r)}{r^2} by incorporating general relativistic effects through the Schwarzschild metric for interior solutions. In their seminal work, Oppenheimer and Volkoff solved the Einstein field equations for a static, isotropic fluid, yielding the pressure gradient as:

\frac{dP}{dr} = -\frac{G M(r) \rho(r)}{r^2} \left(1 + \frac{P}{\rho c^2}\right) \left(1 + \frac{4\pi r^3 P}{M(r) c^2}\right) \left(1 - \frac{2 G M(r)}{r c^2}\right)^{-1},

where P(r) is the pressure, \rho(r) is the total energy density (including rest mass), M(r) is the enclosed mass within radius r, G is the gravitational constant, and c is the speed of light. This is coupled to the mass continuity equation:

\frac{dM}{dr} = 4\pi r^2 \rho(r),

with boundary conditions P(0) finite at the center, M(0) = 0, and P(R) = 0 at the stellar radius R. The three parenthetical factors in the TOV equation represent corrections for the gravitational redshift of energy density, the contribution of internal pressure to the gravitational field, and the enhancement of gravitational attraction due to spacetime curvature, respectively.^[40] Numerical solutions to the TOV equation are obtained by integrating outward from the stellar center using a specified equation of state (EOS) that relates pressure to density, P = P(\rho). Typically, this involves discretizing the coupled differential equations via methods like the Runge–Kutta integrator, starting with a trial central density and pressure, and iterating until the surface condition is met. Stability criteria emerge from analyzing sequences of solutions parameterized by central density: stable configurations correspond to increasing mass with density up to a maximum, beyond which the mass-radius curve turns over, signaling dynamical instability against collapse, as dictated by the relativistic structure.^[41]^[42] The TOV equation implies a maximum mass for neutron stars, beyond which no stable equilibrium exists, arising from the relativistic corrections that amplify gravitational binding. In the non-relativistic limit (weak gravity, soft EOS), this yields a Chandrasekhar-like limit of approximately 1.4 solar masses (M_\odot), analogous to white dwarf stability. For realistic neutron star EOS, which are stiffer at high densities due to nuclear repulsion, the limit increases to 2–3 M_\odot, depending on the microphysics; softer EOS produce lower maxima, while stiffer ones allow heavier stars before instability sets in.^[42]^[43] In recent applications as of 2025, solutions to the TOV equation have incorporated gravitational wave data from binary neutron star mergers, such as GW170817, to constrain EOS parameters through tidal deformability measurements, which quantify how internal structure responds to external tidal fields and refine the maximum mass limit to around 2.1–2.5 M_\odot with improved precision. These efforts use Bayesian inference on TOV-integrated models to jointly fit multi-messenger observations, including pulsar timing masses and kilonova ejecta, tightening bounds on neutron star compactness.^[44]^[45]

Rotation

Spin evolution

Neutron stars are born with initial spin periods typically around tens of milliseconds, with the distribution peaking at approximately 40–50 ms, as a result of angular momentum conservation from the collapsing progenitor core during a core-collapse supernova. This inheritance from the pre-collapse stellar core, which rotates with periods of tens to hundreds of seconds, leads to proto-neutron stars spinning at rates up to several hundred hertz shortly after formation. Recent population studies of young pulsars associated with supernova remnants, using hierarchical Bayesian inference on spin periods and derivatives, derive an initial spin period distribution that peaks at approximately 40–50 ms, with 90% of neutron stars born with periods less than about 0.5–0.8 s.^[46]^[47] For isolated neutron stars, the primary mechanism of rotational deceleration is the emission of magnetic dipole radiation from the aligned rotator, which applies a braking torque that extracts angular momentum. The vacuum dipole torque is expressed as N_{\rm vac} = -\frac{B^2 R^6 \Omega^3 \sin^2 \chi}{6 [c](/page/Speed_of_light)^3}, where B is the surface magnetic field strength, R the stellar radius, \Omega = 2\pi/P the angular frequency, \chi the obliquity angle between the magnetic and rotation axes, and [c](/page/Speed_of_light) the speed of light; this torque can be augmented by magnetospheric wind contributions. The resulting period derivative is typically \dot{P} \sim 10^{-15} s/s for standard parameters (B \sim 10^{12} G, R \sim 10 km).^[48]^[49] The characteristic spin-down age, derived assuming a constant braking torque and negligible birth period compared to the current one, is given by \tau_c = \frac{P}{2 \dot{P}}, yielding estimates of $10^4 to $10^7 years for the majority of radio pulsars and serving as a proxy for true age when evolutionary effects like magnetic field decay are minimal. In binary systems, however, accretion from a companion can counteract or reverse this spin-down by transferring angular momentum through a disc, with the accretion torque scaling as N_{\rm acc} \propto \dot{M} \sqrt{G M R_m}, where \dot{M} is the mass accretion rate, M the neutron star mass, and R_m the inner magnetospheric radius; this process recycles neutron stars into fast-spinning millisecond pulsars over gigayears.^[48]^[50] Measurements of proper motions for Galactic binary neutron stars, combined with distances and orbital tracing back to the Galactic plane, have provided kinematic ages in recent studies (2024), often ranging from tens to hundreds of megayears and revealing systematic discrepancies with spin-down ages that constrain birth velocities and initial spin states.^[51]

Glitches and anti-glitches

Neutron star glitches are sudden, discontinuous increases in the rotation frequency of the star, typically characterized by fractional changes in angular velocity of ΔΩ/Ω ≈ 10^{-9} to 10^{-6}, though some events reach up to 10^{-3}.^[52] These events are followed by a partial recovery phase over timescales of hours to days, where the spin frequency gradually decreases as the system relaxes.^[53] The Vela pulsar (PSR J0835−4510) is the most prominent example, exhibiting large glitches approximately every three years.^[52] The primary mechanism for glitches involves the transfer of angular momentum from a superfluid neutron component in the stellar interior to the rigidly rotating charged crust and magnetosphere.^[53] In this model, superfluid vortices are pinned to the nuclear lattice in the inner crust, allowing the superfluid to lag behind the decelerating crust due to external torques. When the pinning force is overcome—through vortex creep or sudden unpinning—a large number of vortices (on the order of billions) move outward, effectively speeding up the observed rotation of the crust.^[54] This process couples the superfluid reservoir, which may involve up to 10% of the star's moment of inertia, to the observable components.^[55] An alternative explanation posits starquakes, where accumulated strain in the elastic crust from spin-down leads to sudden cracking and release of stored elastic energy, imparting a spin-up.^[56] However, this mechanism is insufficient to explain large glitches, as the crust's moment of inertia is too small (less than 1% of the total) to provide the required angular momentum without invoking additional superfluid contributions.^[55] Anti-glitches, the rarer counterparts to glitches, manifest as abrupt spin-down events with fractional changes ΔΩ/Ω ≈ -10^{-7} to -10^{-4}, observed primarily in magnetars.^[57] Possible mechanisms include sudden outward motion of superfluid vortices due to reconnection or avalanches in the inner crust, transferring angular momentum from the crust to the superfluid.^[58] Magnetic reconnection in the highly twisted fields of magnetars may also drive these events by altering internal torques.^[59] In 2024, the first anti-glitch was detected in a rotation-powered pulsar (PSR B0540-69), suggesting internal reconfiguration akin to magnetar events and broadening the phenomenon beyond highly magnetized objects.^[57] Subsequent detections include an anti-glitch in PSR J1835-1106 reported in 2025.^[60] These occurrences often lack radiative counterparts, implying origins deep within the neutron star interior.^[61]

Radiation and Emission

Pulsar mechanisms

Pulsars are observed as rapidly varying sources of radiation due to the lighthouse model, in which a rotating neutron star emits beamed radiation from regions near its magnetic poles, sweeping across the observer's line of sight as the star spins. This model, proposed shortly after the discovery of pulsars, explains the periodic pulses as the result of a narrow beam, analogous to a lighthouse beacon, produced by co-rotating plasma in the magnetosphere along open magnetic field lines.^[62] The beam originates from the polar cap, the region at the magnetic pole where field lines are open to the light cylinder, and its visibility depends on the viewing angle relative to the magnetic and rotation axes.^[63] In the polar cap acceleration regions, strong electric fields parallel to the magnetic field lines accelerate charged particles to relativistic energies, enabling pair production and curvature radiation as key processes. Primary electrons, sourced from the neutron star surface or ionospheric plasma, are accelerated in vacuum gaps or partially screened regions above the polar cap, where the corotation electric field induces a potential drop of order 10^{12}-10^{14} volts. These electrons emit high-energy gamma rays via curvature radiation as they follow curved magnetic field lines, and in the strong magnetic fields (B ~ 10^{12} G), these photons convert into electron-positron pairs through magnetic pair production when their perpendicular energy exceeds ~2 m_e c^2 / sinθ, where θ is the angle to the field. The resulting secondary pairs further cascade, screening the electric field and forming sparks that structure the emission into sub-beams, consistent with observed pulse profiles. Radio emission from pulsars arises primarily from coherent plasma processes in the magnetosphere, where relativistic pair plasma generates low-frequency waves through instabilities. In the coherent curvature emission or plasma maser models, bunches of charged particles radiate in phase, achieving brightness temperatures exceeding 10^{25} K, far beyond incoherent synchrotron limits.^[64] These processes occur in the polar cap or slot gap regions, involving linear acceleration of pair plasma along field lines, leading to wave growth via two-stream or Weibel instabilities that convert into electromagnetic modes.^[65] The frequency-dependent pulse profiles, broader at low frequencies due to propagation effects in the plasma, support emission from altitudes of ~10-100 stellar radii, with conical beams evolving into fan-like structures at higher frequencies.^[64] For high-energy pulsars, gamma-ray emission originates from acceleration gaps in the outer magnetosphere, where unscreened electric fields in the slot or outer gap accelerate particles to produce GeV photons via curvature radiation. In the outer gap model, the gap extends from near the null surface to the light cylinder, with pair production limiting its size and enabling steady-state operation for young pulsars like the Crab. Primary electrons radiate gamma rays that pair-produce in the strong fields, creating a self-sustaining cascade, while inverse Compton scattering off soft photons can contribute at higher energies.^[66] This mechanism explains the phase-resolved spectra and off-pulse emission seen in Fermi LAT observations of over 200 gamma-ray pulsars. The geometry of pulsar emission depends on the alignment between the magnetic dipole axis and the rotation axis, with orthogonal rotators (near 90° inclination) producing wider beams and more complex profiles, while aligned rotators (near 0°) yield narrower, less detectable pulses. Electromagnetic torques from magnetic dipole radiation preferentially align the axes over the pulsar's lifetime, reducing the inclination angle and narrowing pulse widths with age, as evidenced by statistical studies of radio pulsar populations.^[67] This evolution, driven by the torque's dependence on sin²α (where α is the inclination), proceeds on timescales of 10^6-10^7 years for typical pulsars.^[68] Recent observations from the Imaging X-ray Polarimetry Explorer (IXPE) in 2024-2025 have provided polarization data that probe the magnetosphere geometry, revealing non-uniform polarization position angles across pulse phases for the Crab pulsar, consistent with emission from a striped wind or outer magnetosphere.^[69] These measurements, with polarization degrees up to 20% in the bridge pulse, support models of current sheets and field line tangling, offering insights into the transition from closed to open field regions.

Thermal and non-thermal spectra

The thermal emission from isolated neutron stars primarily originates from their cooling surfaces, manifesting as blackbody-like spectra with effective temperatures typically ranging from $10^5 to $10^6 K for young to middle-aged objects. These spectra are shaped by thin atmospheric layers, where models incorporating partially ionized hydrogen or helium reveal deviations from ideal blackbody radiation due to opacity effects and gravitational redshift.^[70] For instance, in hydrogen atmospheres, the emergent flux is enhanced at higher energies compared to a Planck function, allowing fits to observed X-ray data that constrain surface temperatures and radii more accurately than simple blackbody approximations.^[71] Atmosphere models distinguish between gaseous hydrogen layers and condensed matter surfaces, with the latter predicted for older neutron stars where heavy elements form solid phases. Spectral fitting of sources like the "Magnificent Seven" isolated neutron stars, discovered via ROSAT surveys, often favors thin hydrogen atmospheres over condensed surfaces, as the latter predict harder spectra with reduced low-energy flux. Absorption features in these spectra, such as broad edges around 0.2–0.5 keV, arise from photoionization in the atmosphere or cyclotron resonance in magnetic fields up to $10^{12} G, providing diagnostics of composition and magnetic geometry.^[72] Non-thermal emission in isolated neutron stars arises predominantly from synchrotron radiation in the magnetosphere, where relativistic electrons accelerated along open field lines produce power-law spectra extending from radio to gamma rays.^[73] These components, observed as tails in X-ray spectra with photon indices \Gamma \approx 1–2, dominate over thermal emission in younger pulsars and reveal particle acceleration mechanisms near the polar caps. In cooling curves derived from ROSAT and Chandra observations of the Magnificent Seven, non-thermal contributions are minimal, allowing clean thermal fits that trace neutrino and photon luminosity evolution over ages $10^4–10^6 years.

Binary Systems

X-ray binaries

Neutron star X-ray binaries are classified into low-mass X-ray binaries (LMXBs), where the companion star has a mass less than about 1 M_\odot, and high-mass X-ray binaries (HMXBs), where the companion exceeds 10 M_\odot. In LMXBs, the neutron star accretes material primarily through Roche-lobe overflow from the companion, leading to the formation of an accretion disk that channels matter toward the neutron star. In contrast, HMXBs typically involve accretion from the dense stellar wind of the massive companion, though Roche-lobe overflow can occur in some cases with Be-star companions. The accretion process in these systems is governed by the neutron star's strong magnetic field, which disrupts the inner accretion disk at the magnetospheric radius. In the propeller regime, particularly prominent in rapidly rotating neutron stars, infalling material is ejected by the rotating magnetosphere, preventing accretion when the magnetospheric radius exceeds the co-rotation radius. The maximum sustainable accretion rate is limited by the Eddington luminosity, corresponding to approximately $2 \times 10^{-8} \, M_\odot \, \mathrm{yr}^{-1} for a typical 1.4 M_\odot neutron star, beyond which radiation pressure halts further inflow. X-ray emission in these binaries arises from the interaction of accreting material with the neutron star's surface or magnetosphere. In LMXBs with weaker magnetic fields, X-rays are predominantly produced in the boundary layer, where the Keplerian disk matter spreads onto the neutron star surface and is heated to keV temperatures through frictional dissipation. In HMXBs with stronger fields, accretion funnels along magnetic field lines to the poles, forming shocks in the magnetosphere or accretion column that accelerate electrons and generate hard X-ray spectra through Comptonization. The spin evolution of neutron stars in LMXBs is characterized by "recycling," where sustained accretion transfers angular momentum, spinning up the star from initial periods of seconds to millisecond timescales over \sim 10^7 years. This process correlates with orbital periods, with faster rotators (up to 620 Hz) found in systems with longer orbits and higher accretion rates, eventually yielding active or quiescent millisecond pulsars after the accretion phase ends. In HMXBs, spin changes are more variable due to episodic wind accretion, often resulting in slower evolution without reaching millisecond periods. Observationally, LMXBs are subclassified using color-color diagrams, which plot soft versus hard X-ray colors to reveal track shapes during intensity variations. Atoll sources trace island-like or atoll-shaped paths, associated with moderate accretion rates and banana-to-island state transitions, while Z-sources follow a Z-like track at near-Eddington luminosities, reflecting higher mass transfer and distinct spectral states. These patterns, observed with instruments like the Rossi X-ray Timing Explorer, highlight differences in accretion disk properties and Comptonization, though some overlap suggests a continuum rather than strict separation. Recent multiwavelength campaigns, including X-ray observations from Chandra, XMM-Newton, and eROSITA alongside optical data, have identified new candidate transitional millisecond pulsars in sub-luminous disk states, such as 4FGL J0639.1-8009 and 4FGL J1824.2+1231, revealing variable power-law X-ray spectra and accretion signatures consistent with switching between rotation- and accretion-powered phases.

Mergers and nucleosynthesis

Neutron star binaries inspiral over billions of years due to energy loss via gravitational wave emission, with the early phases accurately described by the post-Newtonian approximation that expands general relativity in powers of the orbital velocity over the speed of light. As the stars approach merger, tidal interactions become significant, leading to deformation and eventual tidal disruption when the orbital separation reaches approximately 10-20 km, comparable to the neutron star radii of 10-15 km.^[74] This disruption ejects neutron-rich material and marks the transition from inspiral to the highly relativistic merger phase, where full numerical relativity simulations are required to model the dynamics. The gravitational wave signals from these mergers feature a characteristic "chirp" during inspiral, where the frequency and amplitude increase as the orbit tightens, allowing detectors like LIGO and Virgo to infer the chirp mass, defined as \mathcal{M} = \frac{(m_1 m_2)^{3/5}}{(m_1 + m_2)^{1/5}}, typically around 1.2 solar masses for observed events such as GW170817. Post-merger, the signal includes a ringdown phase dominated by quasi-normal modes of the remnant, providing probes of its structure and the equation of state. These signals have been detected multiple times, confirming the binary neutron star nature through the low chirp masses and absence of higher-mass black hole signatures. During the merger, dynamical ejecta with low electron fraction (Y_e \lesssim 0.25) is launched in neutron-rich outflows, enabling rapid neutron capture (r-process) nucleosynthesis that produces heavy elements beyond iron, including gold and uranium.^[75] This ejecta powers a kilonova transient, with the light curve showing a blue component from lanthanide-poor material (higher Y_e, more transparent, peaking in optical bands) and a red component from lanthanide-rich ejecta (lower Y_e, opaque, shifting to infrared). Observations of GW170817's counterpart, AT 2017gfo, revealed this two-component evolution, confirming r-process contributions to galactic enrichment. The merger remnant forms a hypermassive neutron star, with mass exceeding the Tolman-Oppenheimer-Volkoff limit for non-rotating stability (around 2-2.5 solar masses) but temporarily supported by differential rotation and thermal pressure.^[76] Depending on the total mass and equation of state, it may collapse to a black hole within milliseconds to seconds via gravitational instability, or stabilize as a massive neutron star if below the threshold.^[76] For GW170817, analyses favor a short-lived remnant collapsing promptly to a black hole, consistent with the lack of prolonged post-merger emission.^[77] As of November 2025, LIGO-Virgo-KAGRA has detected approximately 2–3 binary neutron star mergers during the O4 observing run, roughly doubling the confident detections since the O3 run.^[78] Advanced numerical simulations, now extending to approximately 1.5 seconds post-merger—the longest to date—reveal jet formation from magnetized accretion around the remnant black hole, driving potential gamma-ray bursts.^[79] These simulations also compute precise r-process yields, showing variations in heavy element production based on ejecta composition and neutron star equation of state, with lanthanide fractions influencing kilonova colors and galactic chemical evolution models.^[80]

Populations and Observations

Galactic distribution

Neutron stars are estimated to number between $10^8 and $10^9 in the Milky Way, based on models integrating the stellar birthrate, evolutionary lifetimes, and spatial distributions derived from observed samples.^[81] Of these, only a small fraction—approximately $10^5 to $10^6—are potentially observable as radio pulsars, limited by beaming effects, spin-down, and interstellar absorption that render most undetectable. The birthrate of neutron stars aligns closely with the core-collapse supernova rate in the Galaxy, estimated at 1 to 2 per century, primarily from massive stars exceeding 8 solar masses.^[82] The majority of neutron stars reside in the Galactic disk, with a vertical scale height of approximately 50–100 pc reflecting their birth in the thin disk amid star-forming regions, though natal kicks can disperse them further.^[81] A smaller halo population arises from binary neutron star mergers, where dynamical interactions and kicks eject systems from the disk, contributing to a diffuse distribution extending to several kiloparsecs above and below the plane.^[83] The age distribution of observable neutron stars is skewed toward younger objects, as spin-down reduces radio luminosity over time, making older isolated neutron stars harder to detect beyond characteristic ages of $10^6 to $10^7 years.^[81] However, recycled millisecond pulsars—spun up via accretion in binaries—represent an older cohort, with ages up to several billion years, providing insights into long-term evolution. Recent analyses using Gaia proper motions have backtracked pulsar trajectories to their birth locations, confirming associations with spiral arms such as the Perseus and Scutum-Centaurus structures.

Distances and kinematics

Determining the distances to neutron stars is crucial for understanding their luminosities, ages, and evolutionary contexts, with several complementary methods employed depending on the object's detectability across wavelengths. For radio pulsars, the primary technique involves measuring the dispersion measure (DM), which quantifies the integrated column density of free electrons along the line of sight, causing a frequency-dependent delay in radio pulse arrival times. Distances are then inferred by dividing the DM by models of the Galactic electron density distribution, such as the YMW16 or NE2001 models, yielding estimates accurate to within 20-50% for nearby sources.^[84] Astrometric parallax measurements provide geometric distances independent of interstellar medium assumptions; the European Space Agency's Gaia mission has delivered parallaxes for hundreds of pulsars via optical counterparts, while very long baseline interferometry (VLBI) achieves microarcsecond precision for radio-bright objects, as demonstrated by the PSRπ project using the Very Long Baseline Array (VLBA).^[85]^[86] For X-ray emitting neutron stars, interstellar absorption in spectra—modeled through photoelectric edges of elements like neon and magnesium—constrains the hydrogen column density N_H, which correlates with distance via empirical dust-to-gas ratios, particularly effective for isolated sources like RX J1856.5-3754.^[87]^[88] Neutron stars acquire high velocities, known as natal kicks, primarily during their formation in core-collapse supernovae, where asymmetries in the explosion impart momenta via hydrodynamic instabilities or neutrino-driven convection. These kicks typically range from 100 to 500 km/s, with evidence from pulsar proper motions and supernova remnant offsets supporting explosion anisotropies that accelerate the remnant while ejecting mass at ~10,000 km/s.^[89]^[90] In binary systems, kicks are generally lower, peaking at 40-50 km/s according to log-normal distributions derived from systemic velocities of double neutron star binaries, as these systems survive disruption only if the velocity impulse is moderated by the companion's gravitational binding.^[51] Proper motions, measured via long-term radio or optical astrometry, reveal the transverse components of these velocities, with many neutron stars exhibiting speeds exceeding 300 km/s and up to 1000 km/s or more. Runaway pulsars like PSR B1508+55, with a transverse velocity of ~1100 km/s inferred from VLBI observations, exemplify these high-speed ejections, implying kicks that unbound them from birth clusters or binaries.^[91]^[92] Such measurements, combined with radial velocity estimates from pulse profile changes or Doppler shifts, enable three-dimensional velocity vectors, highlighting a Maxwellian-like distribution biased toward higher speeds compared to progenitor stars. Trajectories of neutron stars often link them to their supernova remnant (SNR) origins, with proper motion projections backward in time aligning positions with remnant centers for young objects like those in the Crab or Vela SNRs. Bow shocks form when fast-moving neutron stars (~100-500 km/s) ram into the interstellar medium, creating asymmetric Hα or X-ray emitting structures, as observed in PSR J0437-4715 and the Guitar Nebula from PSR B2224+65, where shock apex offsets constrain space velocities and magnetic field strengths.^[93]^[94] In 2025, a 15-year monitoring campaign using the Very Large Array traced the proper motion of the young pulsar PSR J0538+2817, precisely mapping its post-kick trajectory and confirming a natal velocity of ~450 km/s consistent with asymmetric explosion models from its host SNR G180.1-1.8. Separately, analysis of a 2023 gamma-ray burst revealed the first "heartbeat" signals—rapid oscillations at 909 Hz—from a newborn, isolated millisecond magnetar, indicating spin-up from a proto-neutron star phase shortly after formation.^[95]^[96]

Subtypes

Radio pulsars

Radio pulsars represent the most common observed subtype of neutron stars, characterized by their periodic emission of radio pulses due to a rotating beam of radiation sweeping across the observer's line of sight, akin to a lighthouse effect. Approximately 4,000 radio pulsars have been discovered within the Milky Way galaxy as of November 2025, primarily through targeted radio surveys that detect these coherent emissions at frequencies typically between 100 MHz and several GHz.^[97] These objects are magnetized neutron stars with strong dipole magnetic fields that accelerate charged particles, producing the observed radio emission via curvature radiation and subsequent plasma processes. Key observational properties of radio pulsars include their spin periods, which range from about 30 milliseconds to 10 seconds, reflecting a broad evolutionary span from young to aged objects. The pulses arrive with a delay that increases at lower observing frequencies due to dispersive effects from free electrons in the interstellar medium; this is quantified by the dispersion measure (DM), defined as the integral of the electron density along the propagation path, typically expressed in units of pc cm^{-3}. Values of DM, often between 10 and 500 pc cm^{-3} for Galactic pulsars, enable distance estimates by modeling the electron distribution in the Galaxy, confirming that most radio pulsars reside within a few kiloparsecs of the Sun.^[98]^[98] The radio pulsar population divides into two primary classes based on spin and magnetic field characteristics. Normal pulsars, comprising the majority, are young objects with spin periods of 0.1–10 s and inferred surface magnetic fields of approximately 10^{12} G, derived from their spin-down rates assuming magnetic dipole braking. In contrast, millisecond pulsars form a distinct subset of older, rapidly rotating objects with periods under 30 ms and weaker fields around 10^{8}–10^{9} G; these are thought to originate from low-mass X-ray binaries where accretion torques "recycle" the neutron star, spinning it up before it resumes isolated radio emission. About 10–15% of known radio pulsars are millisecond types, offering exceptional rotational stability for precision timing applications.^[99]^[99] Ongoing radio surveys have dramatically expanded the known population, with major contributions from the Parkes Multibeam Pulsar Survey, which discovered over 800 pulsars in the 1990s–2000s using a 20-cm multibeam receiver; the Arecibo Pulsar ALFA (PALFA) survey, which has identified hundreds more through drift-scan observations at 1.4 GHz; and the Canadian Hydrogen Intensity Mapping Experiment (CHIME), a northern-hemisphere array operational since 2018 that detects low-frequency signals across a wide field. Collectively, recent efforts yield approximately 200–300 new radio pulsar discoveries per year as of 2025, enhancing statistical studies of the Galactic population.^[100]^[101]^[102] Radio pulsars evolve by losing rotational energy through magnetic dipole radiation, gradually lengthening their periods until they cross the "death line" in the spin period–period derivative (P–\dot{[P](/page/P′′)}) diagram. This boundary marks the minimum P and \dot{[P](/page/P′′)} values where the magnetosphere's accelerating electric potential drops below the threshold for sustaining electron–positron pair production cascades, which are essential for generating the coherent radio emission. Below this line, typically around P \approx 5–10 s for standard models, pulsars become radio-quiet, though some outliers challenge theoretical predictions.^[103]^[104]

Magnetars and exotic variants

Magnetars represent a subtype of neutron stars characterized by exceptionally strong magnetic fields exceeding $10^{14} gauss, which power their emission through magnetic field decay rather than rotational energy loss.^[105] These objects manifest observationally as soft gamma repeaters (SGRs), which exhibit recurrent gamma-ray bursts, and anomalous X-ray pulsars (AXPs), known for their steady X-ray luminosity punctuated by irregular flares.^[106] The decay of the internal magnetic field induces starquakes in the solid crust, releasing energy that heats the surface and drives outbursts, with typical burst energies reaching $10^{40} ergs.^[105] As of 2025, approximately 30 magnetars have been confirmed in the Milky Way, highlighting their rarity compared to the thousands of known radio pulsars.^[107] Exotic variants of neutron stars include hypothetical quark stars, proposed under the strange matter hypothesis where deconfined up, down, and strange quarks form a stable phase more compact than neutron matter.^[108] In this scenario, neutron stars could convert to strange quark stars if strange matter proves more stable at high densities, potentially explaining compact objects with radii below 10 km that challenge standard neutron star models.^[109] Isolated neutron stars like Calvera (1RXS J141256.0+792204) exemplify unusual properties, featuring a high surface temperature of about $4 \times 10^6 K and no detectable magnetic field or pulsar activity, suggesting it as a possible aged remnant with minimal field evolution.^[110] Central compact objects (CCOs) are another distinct class, typically found near the centers of young supernova remnants, emitting soft thermal X-rays from weak magnetic fields below $10^{11} gauss without radio or gamma-ray pulsations.^[111] These objects, numbering around a dozen, likely represent young neutron stars born with suppressed fields, possibly due to fossil configurations from progenitor stars, and serve as probes of early cooling phases.^[112] Recent observations in 2025 have identified candidates for transitional objects potentially involving quark phases, such as the gamma-ray burst GRB 240529A, whose multi-episode emission and X-ray afterglow suggest the collapse of a magnetar into a strange star.^[109] Laboratory analogs have confirmed vortex structures in superfluids mimicking neutron star interiors, providing indirect support for magnetar crust dynamics under extreme conditions. These exotics remain tentative, with ongoing multi-wavelength studies essential to distinguish them from standard neutron stars.

Historical Development

Theoretical foundations

In 1932, Lev Landau developed the theory of degenerate Fermi gases in stellar interiors, calculating the maximum mass limit for white dwarf stars supported by electron degeneracy pressure at approximately 1.5 solar masses, and suggesting that more massive configurations would require support from neutron degeneracy, implying the existence of denser "neutron-like" stars. This work laid a foundational conceptual framework for compact objects beyond white dwarfs, though it predated the neutron's discovery by James Chadwick earlier that year.^[113] Building on this, astronomers Walter Baade and Fritz Zwicky proposed in 1934 that core-collapse supernovae—distinguished from classical novae by their immense energy output—represent the transition of ordinary stars into compact remnants composed primarily of neutrons, forming stable "neutron cores" with radii around 10 kilometers and densities comparable to atomic nuclei. Their hypothesis linked these events to cosmic ray production and stellar evolution, predicting that such remnants could persist indefinitely due to neutron degeneracy pressure balancing gravity. In 1939, J. Robert Oppenheimer and George Volkoff provided the first rigorous theoretical models of these neutron stars by solving the general relativistic hydrostatic equilibrium equations—now known as the Tolman–Oppenheimer–Volkoff equation—for a degenerate neutron gas equation of state, demonstrating stable configurations with masses up to about 0.7 solar masses and radii of roughly 10 kilometers. These calculations confirmed the viability of neutron-supported stars against gravitational collapse under the assumed non-interacting Fermi gas model, though the precise equation of state remained approximate given the nascent understanding of nuclear forces. During the 1960s, John Archibald Wheeler and his collaborators advanced these ideas through pioneering numerical simulations of stellar collapse, exploring how iron-core implosions in massive stars could halt at neutron star densities or proceed to black hole formation depending on the core mass and nuclear physics. Wheeler's work emphasized the dynamical processes, including shock waves and neutrino emission, that might stabilize the nascent neutron star post-collapse.^[114] Pre-discovery theoretical discussions centered on the viability of degenerate neutron matter, with uncertainties in the high-density equation of state—particularly strong nuclear interactions and possible phase transitions—fueling debates over whether stable neutron stars could exist or if they would inevitably collapse further due to insufficient pressure support. These concerns persisted into the mid-1960s, as improved nuclear models sometimes suggested instability for certain mass ranges, tempering enthusiasm for the concept despite the earlier analytical successes.

Key discoveries and observations

The discovery of the first pulsar, designated CP 1919 (now PSR B1919+21), occurred on November 28, 1967, when graduate student Jocelyn Bell, under the supervision of Antony Hewish, detected regular radio pulses with a period of 1.337 seconds using a large radio telescope array at the Mullard Radio Astronomy Observatory in Cambridge, England. This breakthrough, published in 1968, revealed a new class of rapidly rotating neutron stars emitting beamed radio emission, fundamentally altering understandings of stellar remnants and earning Hewish and Martin Ryle the 1974 Nobel Prize in Physics (though Bell was not included).^[115] Subsequent surveys rapidly identified dozens more pulsars, confirming their galactic distribution and association with neutron stars. In 1974, Russell Hulse and Joseph Taylor discovered the first binary pulsar, PSR B1913+16 (the Hulse-Taylor pulsar), using the Arecibo Observatory, a system consisting of two neutron stars orbiting each other with a 7.75-hour period. Precise timing observations over subsequent years revealed an orbital decay rate matching general relativity's prediction for energy loss via gravitational wave emission to within 0.2%, providing the first indirect confirmation of gravitational waves and earning Hulse and Taylor the 1993 Nobel Prize in Physics. The 1982 discovery of the first millisecond pulsar, PSR B1937+21, with a rotation period of 1.5578 milliseconds, by Don Backer and colleagues at Arecibo, introduced a new population of rapidly spinning neutron stars and challenged models of isolated pulsar slowdown. This finding prompted M. Ali Alpar et al. to propose the "recycling" paradigm, in which old, slowly rotating neutron stars in binary systems accrete mass and angular momentum from companions, spinning up to millisecond periods before emerging as isolated fast rotators. Observations of over 100 such pulsars since have validated this mechanism, linking them to low-mass X-ray binaries. During the 1990s, soft gamma repeaters (SGRs) and anomalous X-ray pulsars (AXPs), first identified in the 1970s and 1980s as enigmatic bursting sources, were recognized as a unified class powered by decay of magnetic fields exceeding 10^14 gauss, termed magnetars. Key evidence came from Chryssa Kouveliotou et al.'s 1998-1999 observations with the Rossi X-ray Timing Explorer, which detected correlated radio afterglows, X-ray pulsations, and giant flares from sources like SGR 1900+14, linking their properties to neutron stars with extreme magnetism rather than rotation. This identification expanded neutron star subtypes and explained their high-energy outbursts without invoking accretion. The 2017 multimessenger event GW170817, detected by the LIGO and Virgo observatories on August 17, represented the merger of two neutron stars at a distance of 40 megaparsecs, accompanied by a short gamma-ray burst (GRB 170817A) observed by Fermi and INTEGRAL, followed by a kilonova AT 2017gfo visible across wavelengths from ultraviolet to infrared.^[116] This observation confirmed neutron star mergers as sites of rapid neutron capture (r-process) nucleosynthesis, producing heavy elements like gold, and provided the first direct evidence of such events driving gamma-ray bursts.^[116] Complementing this, NASA's Neutron star Interior Composition Explorer (NICER) mission yielded precise radius measurements from 2019 to 2023, constraining the radius of a 1.4 solar mass neutron star to 12.71^{+1.14}_{-1.08} km based on 2021 observations of the approximately 2.08 solar mass pulsar PSR J0740+6620 and refined values for PSR J0030+0451, constraining the neutron star equation of state and indicating radii around 12-13 km for typical masses. From 2023 to 2025, LIGO-Virgo-KAGRA detections of neutron star mergers roughly doubled the known sample beyond GW170817, including GW230529—a neutron star paired with a 2.5-4.5 solar mass mass-gap object—enhancing merger rate estimates to about 250 Gpc^{-3} yr^{-1}.^[117] In 2025, analysis of the long-duration gamma-ray burst GRB 230307A revealed quasi-periodic "heartbeat" signals oscillating at 909 Hz for 160 milliseconds, interpreted as pulsations from a newborn millisecond magnetar formed in the merger or collapse.^[96] Concurrently, James Webb Space Telescope observations in 2024 of the SN 1987A remnant detected compact infrared emission consistent with a dust-enshrouded neutron star, providing the first direct evidence of the compact object at its core after 37 years. In January 2025, the discovery of ASKAP J1832−0911, a radio-emitting neutron star with a rotation period of approximately 44 minutes—the longest known for such an object—challenged existing models of pulsar spin-down and emission mechanisms.^[118]

Notable Examples

Iconic pulsars

The Crab Pulsar (PSR B0531+21), located in the Crab Nebula supernova remnant, exemplifies a young, energetic neutron star with a rotation period of approximately 33 milliseconds and an age of about 970 years, derived from the historical supernova event of 1054 CE.^[119] It is notable for its multi-wavelength emission, including pulsed optical light detectable with ground-based telescopes and very high-energy (TeV) gamma rays observed by instruments like the Whipple Observatory, making it one of the few pulsars visible across such a broad spectrum. The Crab Pulsar is also renowned for its glitch activity, with over 24 glitches recorded since 1969, including sudden spin-ups that provide insights into the neutron star's internal superfluid dynamics.^[119] Its inferred surface magnetic field strength is around $10^{12} G, typical for young radio pulsars, powering the surrounding nebula through rotational energy loss. The Vela Pulsar (PSR B0833-45), associated with the Vela supernova remnant, rotates every 89 milliseconds and stands out for its high glitch rate, with 21 events documented since 1968, some as large as \Delta \nu / \nu \approx 10^{-6}, where \nu is the spin frequency.^[120] These glitches, occurring roughly every three years, are among the most frequent and well-studied, briefly referencing the broader phenomenon of rotational irregularities in neutron stars. It is particularly bright in gamma rays, with pulsed emission detected from MeV to GeV energies by satellites like Fermi-LAT, highlighting its role in understanding high-energy particle acceleration. Like the Crab, its magnetic field is estimated at \sim 10^{12} G, sustaining a luminous pulsar wind nebula.^[120] The Hulse-Taylor binary pulsar (PSR B1913+16), discovered in 1974, was the first observed binary neutron star system, featuring a 59-millisecond pulsar orbiting a companion neutron star every 7.75 hours in a highly eccentric orbit. Precise timing measurements revealed the orbital decay rate \dot{P_b} = -2.423 \times 10^{-12} s/s, matching general relativity's prediction for energy loss via gravitational wave emission to within 0.2%, providing the first indirect evidence of gravitational waves. The pulsars' masses are precisely determined at 1.438 M_\odot and 1.389 M_\odot, establishing key constraints on neutron star equations of state and binary evolution. Its magnetic field of about $1.4 \times 10^{12} G underscores the typical dipolar fields in recycled pulsars. In the globular cluster 47 Tucanae (NGC 104), at least 23 millisecond pulsars have been identified, representing a dense population recycled through accretion in the cluster's core. These include short-period examples like 47 Tuc A with a 4.3-millisecond spin, orbiting low-mass companions and exhibiting eclipsing binaries that probe cluster dynamics and binary formation.^[121] Their weak magnetic fields, typically below $10^{10} G, reflect spin-up from prolonged accretion, contrasting with isolated young pulsars and illustrating environmental influences on neutron star properties.^[122] Observations at radio and X-ray wavelengths reveal cluster-wide interactions, such as intrabinary shocks, tying these systems to broader studies of stellar evolution in dense environments.^[121]

Recent detections

In 2025, the LIGO-Virgo-KAGRA collaboration released the fourth gravitational-wave transient catalog (GWTC-4), incorporating 128 new detections from the ongoing observing runs, including several binary neutron star (BNS) mergers that provided unprecedented insights into compact object populations.^[123] Among these, GW250818k stood out as a candidate for a subsolar-mass BNS merger, detected in real-time with electromagnetic follow-up via the transient AT2025ulz, highlighting advances in multi-messenger astronomy.^[124] This event, combined with machine-learning algorithms enabling rapid inference of merger signals, marked the clearest post-merger remnants observed to date, where the hypermassive neutron star remnants exhibited prolonged stability before collapsing, offering constraints on the nuclear equation of state.^[125] These detections, totaling over a dozen BNS events in the catalog, underscored the increasing sensitivity of ground-based observatories to neutron star collisions at cosmological distances.^[126] Analysis published in September 2025 revealed the first 'heartbeat' signal possibly from a millisecond magnetar remnant, detected in the gamma-ray burst GRB 230307A (observed March 7, 2023) approximately 500 million light-years away.^[96] Detected by a combination of radio telescopes and gamma-ray monitors, the signal manifested as quasi-periodic bursts at 909 Hz, indicative of a magnetar formed in a binary neutron star merger. This rapid rotation and the signal's persistence for weeks before fading provided insights into post-merger dynamics and magnetar formation in extreme environments.^[96] Among isolated neutron stars, NASA's Chandra X-ray Observatory identified an eccentric source in May 2025, designated ASKAP J1832-0911, which exhibited bizarre 44-minute pulsations in both radio and X-ray bands, defying conventional models of neutron star behavior.^[127] Located near a supernova remnant in the Scutum constellation, this object—potentially a slowly rotating neutron star or exotic white dwarf—showed irregular intensity variations thousands of times longer than typical pulsar periods, suggesting an unconventional magnetic field geometry or orbital interaction that eludes standard evolutionary pathways.^[128] Complementing this, the X-ray Imaging and Spectroscopy Mission (XRISM) targeted the neutron star in the low-mass X-ray binary GX 13+1 in February 2024, with data analyzed in 2025 revealing an unexpectedly slow and dense wind outflow at just 1% of expected speeds for such systems.^[129] This 'cosmic fog' of ionized material, driven primarily by thermal rather than radiative forces, reached densities among the highest recorded, prompting revisions to accretion disk wind models in neutron star binaries.^[130] Kinematic studies advanced significantly in 2025 with a 15-year astrometric campaign culminating in August, tracking the proper motion of the neutron star in the supernova remnant G18.9–1.1 to precisely map its velocity vector and origin.^[95] Using data from NASA's Chandra X-ray Observatory, astronomers determined a transverse velocity of 264–474 km/s (distance-dependent) imparted during the supernova explosion of its progenitor, revealing asymmetric mass ejection as the dominant mechanism for such high-speed ejections in young neutron stars. This long-term observation not only confirmed the star's birthplace in a distant cluster but also refined models of supernova dynamics and neutron star launch mechanisms.^[95] Simulation efforts reached new milestones in 2025, with the longest-ever general relativistic magnetohydrodynamic simulations of BNS mergers, spanning over 100 milliseconds post-merger, incorporating realistic jet formation and black hole accretion.^[79] Conducted by an international team, these runs demonstrated how magnetic fields amplify relativistic jets in the remnant disk, providing the first detailed view of gamma-ray burst launching in neutron star systems and validating electromagnetic counterparts to gravitational-wave events.^[79] Separately, three-dimensional supernova simulations in February unveiled pathways to the lightest known neutron stars, with masses around 1.17 solar masses, arising from fallback accretion in asymmetric explosions of low-mass progenitors. These models, integrating advanced neutrino transport, explained the rarity of such objects and their role in binary evolution leading to mergers.^[131]

References

[1]
Star Types - NASA Science
Oct 22, 2024 · A neutron star forms when a main sequence star with between about eight and 20 times the Sun's mass runs out of hydrogen in its core. (Heavier ...
[2]
Neutron Stars - Imagine the Universe! - NASA
Neutron stars are formed when a massive star runs out of fuel and collapses. The very central region of the star – the core – collapses, crushing together ...Missing: properties | Show results with:properties
[3]
[PDF] The physics of neutron stars - arXiv
Jun 22, 2011 · Neutron stars, the most compact stars in the Universe, were given this name because their interior is largely composed of neutrons.
[4]
Neutrinos from Core-Collapse Supernova Explosions
In this article we give an overview of the supernova neutrinos emitted from core-collapse supernova explosions as well as the significance of their ...Abstract · Introduction · Scenario of core-collapse... · Observation of supernova...
[5]
Neutrino signal from proto-neutron star evolution: Effects of opacities ...
Feb 14, 2020 · We investigate the impact of charged-current–neutrino processes on the formation and evolution of neutrino spectra during the deleptonization of proto-neutron ...
[6]
Webb Finds Evidence for Neutron Star at Heart of Young Supernova ...
Feb 22, 2024 · The supernova, known as SN 1987A, was a core-collapse supernova, meaning the compacted remains at its core formed either a neutron star or a ...
[7]
[1602.05771] Do electron-capture supernovae make neutron stars ...
Feb 18, 2016 · The main aim is to determine whether these events are thermonuclear or core-collapse supernova explosions, and hence whether NSs are formed by ...
[8]
On the formation of neutron stars via accretion-induced collapse in ...
Feb 7, 2018 · Abstract:We investigate evolutionary pathways leading to neutron star formation through the collapse of oxygen-neon white dwarf (ONe WD) ...
[9]
Formation of neutron stars inside planetary nebulae via accretion ...
Jul 4, 2024 · Abstract:First, I demonstrate the possible formation channel for neutron stars (NSs) via the accretion-induced collapse (AIC) of ...
[10]
The formation of single neutron-stars from double white-dwarf ... - arXiv
Apr 7, 2020 · The merging of double white dwarfs (WDs) may produce the events of accretion-induced collapse (AIC) and form single neutron stars (NSs).
[11]
Evolution of double oxygen-neon white dwarf merger remnant - arXiv
Aug 21, 2023 · Double white dwarf (WD) merger process and their post-merger evolution are important in many fields of astronomy, such as supernovae, gamma-ray ...
[12]
[1410.5431] Super and massive AGB stars - IV. Final fates - arXiv
Oct 20, 2014 · We predict electron-capture supernova (EC-SN) rates for lower metallicities which are significantly lower than existing values from parametric ...
[13]
[2205.03989] Wide binary pulsars from electron-capture supernovae
May 9, 2022 · Neutron stars receive velocity kicks at birth in supernovae. Those formed in electron-capture supernovae from super asymptotic giant branch ...
[14]
Double Neutron Star Populations and Formation Channels
Jul 18, 2019 · Either future observations will fill in the gap in short-period DNSs for 0.3 < e < 0.6, or an alternative formation scenario is required.
[15]
Nucleosynthesis Conditions in Outflows of White Dwarfs Collapsing ...
Dec 3, 2024 · Accretion-induced collapse (AIC) or merger-induced collapse (MIC) of white dwarfs (WDs) in binary systems is an interesting path to neutron star (NS) and ...<|separator|>
[16]
https://arxiv.org/abs/2412.02756
[17]
[2407.06790] Constraining the dense matter equation of state ... - arXiv
Jul 9, 2024 · Pulse profile modeling of X-ray data from NICER is now enabling precision inference of neutron star mass and radius.Missing: 2025 | Show results with:2025
[18]
A More Precise Measurement of the Radius of PSR J0740+6620 ...
Oct 18, 2024 · Here we report an updated radius measurement, derived by fitting models of X-ray emission from the neutron star surface to NICER data accumulated through 2022 ...
[19]
The Neutron Star Mass-Radius Relation and the Equation of State of ...
May 31, 2012 · We show that the radius of a 1.4 solar mass neutron star lies between 10.4 and 12.9 km, independent of assumptions about the composition of the core.
[20]
[PDF] Implications of latest NICER data for the neutron star equation of state
(Received 16 December 2024; accepted 17 January 2025; published 6 February 2025). As an update to our previously performed Bayesian inference analyses of ...
[21]
The Maximum Mass of a Neutron Star - IOPscience
We obtain a secure upper bound on the neutron star mass equal to 2.9 M ☉. We also find that in order to reach the lowest possible upper bound of 2.2 M ☉,
[22]
Multimessenger constraints on the neutron-star equation of state ...
Dec 18, 2020 · We found that the radius of a 1.4 –solar mass neutron star is 11.75 − 0.81 + 0.86 km at 90% confidence and the Hubble constant is 66.2 − 4.2 + 4.4 at
[23]
The prospect of confining the equation of state of neutron stars with ...
Sep 3, 2025 · The radius of neutron stars obtained from the ensemble was confined within a 10.5–14.5 km bound. Among the NICER measurements, PSR J0437–4715 ...
[24]
Insights into Neutron Star Matter: EoS Models and Observations - arXiv
May 8, 2025 · This review classifies EoS models into hadronic matter, hybrid, and quark matter models, analyzing their assumptions, predictions, and constraints.
[25]
Constraining neutron-star matter with microscopic and macroscopic ...
Jun 8, 2022 · Our findings show that constraints from heavy-ion collision experiments show a remarkable consistency with multi-messenger observations and ...
[26]
[PDF] Evolution of neutron star magnetic fields - arXiv
Sep 12, 2021 · In this review we first briefly describe theoretical models of formation and evolution of magnetic field of neutron stars, paying special ...
[27]
[PDF] Magnetic Field Evolution in Neutron Star Crusts: Beyond the Hall Effect
Jan 20, 2022 · Abstract: Neutron stars host the strongest magnetic fields that we know of in the Universe. Their magnetic fields are the main means of ...
[28]
XRISM catches a pulsar's cosmic wind—and sees a surprising result
Oct 27, 2025 · The XRISM observatory is a joint mission for NASA, the European Space Agency (ESA), and the Japan Aerospace Exploration Agency (JAXA). The ...
[29]
[PDF] Neutron Star Cooling - arXiv
We review physical properties important for cooling such as neutrino emission processes and superfluidity in the stellar interior, surface envelopes of light ...
[30]
[PDF] Theory of cooling neutron stars versus observations - arXiv
Dec 12, 2007 · We review current state of neutron star cooling theory and discuss the prospects to constrain the equation of state, neutrino emission and ...<|control11|><|separator|>
[31]
Webb Finds Evidence for Neutron Star at Heart of Young Supernova ...
Feb 22, 2024 · NASA's James Webb Space Telescope has found the best evidence yet for emission from a neutron star at the site of a recently observed supernova.Missing: cooling | Show results with:cooling
[32]
None
Nothing is retrieved...<|control11|><|separator|>
[33]
[2407.11153] Neutron stars and the dense matter equation of state
Jul 15, 2024 · Here, we review theoretical approaches to calculate and model the neutron-star equation of state in various regimes of densities.
[34]
https://arxiv.org/abs/2407.11153
[35]
Equation of state for neutron stars with hyperons and quarks ... - arXiv
Jun 18, 2015 · We construct the equation of state (EoS) for neutron stars explicitly including hyperons and quarks. Using the quark-meson coupling model with ...
[36]
[2503.23413] Neutron stars and Constraints for the Equation of State ...
Mar 30, 2025 · In this proceeding, we review our current understanding of dense isospin symmetric and asymmetric matter and neutron star physics.
[37]
[1302.4499] I-Love-Q - arXiv
Feb 19, 2013 · We find universal relations between the moment of inertia, the Love number and the quadrupole moment that are independent of the neutron star's and quark star' ...
[38]
[2502.05513] The Equation of State of Neutron Stars: Theoretical ...
Feb 8, 2025 · This review explores observational and theoretical constraints on such transitions using multi-messenger astrophysics.
[39]
[PDF] Tolman–Oppenheimer–Volkoff (TOV) Stars
Dec 4, 2012 · The equation is presented as a somewhat intuitive extension of what was covered regarding stars in hydrostatic equilibrium. Furthermore, we.<|separator|>
[40]
[PDF] Relativity and Structure in Neutron Stars
Dec 9, 2011 · We explore the structure of neutron stars using the Tolman-Oppenheimer-Volkoff equation of hydrostatic equilibrium, the solution to which ...Missing: original | Show results with:original
[41]
Maximum mass of neutron stars | Phys. Rev. C
May 9, 2006 · Knowing the EOS, the equilibrium configurations of static neutron stars are obtained by solving the Tolman-Oppenheimer-Volkoff (TOV) equations
[42]
Phase transitions in dense matter and the maximum mass of neutron ...
We study the impact on the maximum possible neutron-star mass of an “exotic” core consisting of non-nucleonic matter.
[43]
Constraining the Generalized Tolman-Oppenheimer-Volkoff (GTOV ...
Aug 19, 2024 · In this work, we constrain the values of the parameters of the Generalized Tolman-Oppenheimer-Volkoff (GTOV) equation through Bayesian inference.Missing: applications | Show results with:applications
[44]
https://arxiv.org/abs/2408.10425
[45]
On the Initial Spin Period Distribution of Neutron Stars - IOPscience
We derive the initial spin period distribution of neutron stars by studying the population of young pulsars associated with supernova remnants.
[46]
The Spin Periods and Rotational Profiles of Neutron Stars at Birth
Aug 22, 2005 · Based on our calculated PNS spins, at ~ 200-300 milliseconds after bounce, and assuming angular momentum conservation, we estimate final neutron ...
[47]
Spin Evolution of Neutron Stars - MDPI
In this paper we review the basics of magneto-rotational properties of neutron stars focusing on spin-up/spin-down behavior at different evolutionary stages.
[48]
energy loss of a rotating magnetized neutron star - Oxford Academic
May 26, 2007 · The projection of the torque K on to the rotation axis KΩ is equal to W/Ω, ΩK=W. The conductivity of a neutron star crust, according to ...
[49]
[PDF] Spin evolution of neutron stars in binary systems
I review our understanding of the evolution of the spin peri- ods of neutron stars in binary stellar systems, from their birth as fast, spin-powered pulsars, ...
[50]
Kinematic constraints on the ages and kick velocities of Galactic ...
For this reason we are interested in estimating the ages of the Galactic BNSs, and we aim to give an overview of their spin-down ages based on braking in ...
[51]
Small glitches and other rotational irregularities of the Vela pulsar
Glitches are sudden increases in the rotation rate ν of neutron stars, which are thought to be driven by the neutron superfluid inside the star. The Vela pulsar ...
[52]
Post-glitch Recovery and the Neutron Star Structure: The Vela Pulsar
Jun 2, 2025 · The post-glitch recovery phase probes the coupling mechanisms between different stellar components, sheds light on the properties of neutron ...
[53]
Pinned vortex hopping in a neutron star crust - Oxford Academic
Pinning of vortices to ions can decouple the superfluid from the crust and create a reservoir of angular momentum. Sudden large-scale unpinning can lead to an ...
[54]
Pulsar Glitches: The Crust is not Enough | Phys. Rev. Lett.
In particular, many young neutron stars exhibit (more or less) regular glitches, where the observed spin rate suddenly increases (for a recent collection of ...Missing: starquakes insufficient
[55]
Neutron Starquakes and Pulsar Periods - Nature
Neutron star crust stress from shape changes can cause a sudden relaxation, a slight shape change, and a jump in angular velocity, similar to observed pulsar ...<|control11|><|separator|>
[56]
Discovery of the first anti-glitch event in the rotation-powered pulsar ...
Apr 28, 2024 · To date, observations of the few known anti-glitches have been made in magnetars or accreting pulsars. This is the first anti-glitch ...
[57]
Antiglitches in accreting pulsars from superfluid vortex avalanches
ABSTRACT. Three sudden spin-down events, termed 'antiglitches', were recently discovered in the accreting pulsar NGC 300 ULX-1 by the Neutron Star Interior.
[58]
simple mechanism for the anti-glitch observed in AXP 1E 2259+586
Based on this model, we predict that to achieve the observed sudden spin-down, a minimum magnetar-like value for the magnetic field strength is needed.
[59]
A Radiatively Quiet Glitch and Anti-glitch in the Magnetar 1E 2259+ ...
Hence, similar to spin-up glitches, anti-glitches can occur silently. This may suggest that these phenomena originate in the neutron star interior, and that ...
[60]
https://arxiv.org/abs/2505.14755
[61]
[PDF] Gold T. Rotating neutron stars as the origin of the pulsating radio ...
Enigmatic, rapidly pulsating radio sources were interpreted as rotating beacons sweeping over the Earth, radiating out from rapidly rotating,.
[62]
[1709.07179] The Nature of Coherent Radio Emission from Pulsars
Sep 21, 2017 · In this article we review the results from high quality observations in conjunction with theoretical models to unravel the nature of coherent ...
[63]
Coherent Emission Mechanisms in Radio Pulsars
A theory of pulsar radio emission generation, in which the observed waves are produced directly by the maser-type plasma instabilities.Missing: review | Show results with:review
[64]
A Three-dimensional Outer Magnetospheric Gap Model for Gamma ...
A three-dimensional pulsar magnetosphere model is used to study the geometry of outer magnetospheric gap accelerators, following seminal work of Romani and ...
[65]
Pulsar magnetic alignment and the pulsewidth–age relation
We test the hypothesis that pulsars evolve through a progressive narrowing of the emission cone combined with progressive alignment of the spin and magnetic ...MAGNETIC ALIGNMENT... · PULSEWIDTH EVOLUTION... · FITTING TO DATA
[66]
Pulsar Magnetic Axis Alignment and Counteralignment. - ADS
A model is proposed for the alignment or counteralignment of the magnetic axis of a pulsar with its symmetry axis by the motion of the magnetic distortion ...
[67]
Crab pulsar: IXPE observations reveal unified polarization ...
We present a phase-dependent analysis of the polarized emission from the Crab pulsar based on three sets of observations by the Imaging X-ray Polarimetry ...
[68]
[astro-ph/9604072] Model neutron star atmospheres with low ... - arXiv
Apr 13, 1996 · We show that the model atmosphere spectra are substantially different from the blackbody spectra. For light element atmospheres, the flux is ...
[69]
Measuring the basic parameters of neutron stars using model ... - arXiv
Oct 23, 2015 · Abstract:Model spectra of neutron star atmospheres are nowadays widely used to fit the observed thermal X-ray spectra of neutron stars.
[70]
[0905.3276] Models of magnetized neutron star atmospheres - arXiv
May 20, 2009 · Observed X-ray spectra of some isolated magnetized neutron stars display absorption features, sometimes interpreted as ion cyclotron lines.Missing: matter | Show results with:matter
[71]
[1406.7151] Synchrotron X-ray emission from old pulsars - arXiv
Jun 27, 2014 · We study the synchrotron radiation as the observed non-thermal X-ray emission from old pulsars (\gtrsim1-10Myr) to investigate the particle ...
[72]
Interpreting binary neutron star mergers
Mar 3, 2021 · We give an overview about general features of the dynamics and the gravitational wave signal during the binary neutron star coalescence.
[73]
r-process nucleosynthesis and kilonovae from hypermassive ...
Relatively neutron-rich ejecta that synthesize a substantial mass fraction of lanthanides, which are opaque to blue light, produce a 'red' kilonova, whereas the ...
[74]
The Aftermath of Binary Neutron Star Mergers: What Remains Behind?
Oct 19, 2021 · The fate of a remnant is dictated by the mass of the two merging neutron stars and the maximum mass a neutron star can support before it collapses to form a ...
[75]
Did rapid spin delay 2017 collapse of neutron stars into black hole?
Mar 1, 2022 · When two neutron stars merged in a 2017 cataclysm, the stars likely were spinning so fast it took about a second for them to collapse to a black hole.Missing: stable | Show results with:stable
[76]
[PDF] The LIGO-Virgo-KAGRA network and observations in the O4 run
Sep 20, 2024 · Probable NS-BH merger. • 230529 : single-detector event. Another possible NS-BH merger. (see next slide …) • 240615 : smallest localization ...
[77]
Breakthrough in simulating how neutron stars collide
May 30, 2025 · This video shows the longest and most complex self-consistent simulation to date of a binary neutron star merger with black hole and jet formation.
[78]
Jet from Binary Neutron Star Merger with Prompt Black Hole Formation
Aug 7, 2025 · We performed the longest numerical-relativity neutrino-radiation magnetohydrodynamics simulation for a binary neutron star merger that ...
[79]
Galactic neutron stars - Astronomy & Astrophysics
Neutron stars (NSs) produced in the Milky Way are supposedly ten to the eighth – ten to the ninth – of which only ∼2 × 103 are observed. Constraining the ...
[80]
[2009.03438] On the rate of core collapse supernovae in the Milky Way
Sep 3, 2020 · A conservative treatment of the errors yields a combined rate R=1.63 \pm 0.46 (100 yr)^{-1}; the corresponding time between core collapse ...
[81]
[2305.04955] Binary neutron star populations in the Milky Way - arXiv
May 8, 2023 · According to our fiducial model, the Square Kilometre Array will detect \sim 20 new BNSs in the Milky Way.Missing: estimate | Show results with:estimate
[82]
Seeking the nearest neutron stars using a new local electron density ...
Nov 26, 2024 · To that end, we revisit pulsar distance estimates extracted from the dispersion measure of pulsar radio waves scattering on electrons.
[83]
Binary Pulsar Distances and Velocities from Gaia Data Release 2
However, the Gaia parallax distance estimates cannot compete with the precision achieved with VLBI for PSR J0437−4715, and the Gaia parallax measurement ...
[84]
[PDF] Microarcsecond VLBI Pulsar Astrometry with PSRπ II. Parallax ...
Apr 22, 2019 · Abstract. We present the results of PSRπ, a large astrometric project targeting radio pulsars using the Very Long Baseline. Array (VLBA).
[85]
THE DISTANCE, MASS, AND RADIUS OF THE NEUTRON STAR IN ...
We determine the amount of the interstellar X-ray absorption by modeling the individual absorption edges of the elements Ne and Mg in the X-ray grating spectra ...
[86]
The role of interstellar absorption on the observed neutron star ...
Sep 11, 2006 · Searching for similar neutron stars with larger distances, one expects larger interstellar absorption resulting in harder X-ray counterparts.
[87]
Calibration of neutron star natal kick velocities to isolated pulsar ...
The associated supernova (SN) explosions eject matter at velocities of ∼ 10 000 km s−1 with a significant degree of asymmetry. Conservation of momentum ...2 Pulsar Velocity... · 4 Binary Evolution Model · 5.1 Bns Detection Rates
[88]
Natal Kick by Early-asymmetrical Pairs of Jets to the Neutron Star of ...
The two kick velocities are of about equal magnitude of ≃450 km s−1, which implies very energetic jets. Such jets can excite gravitational waves that present ...
[89]
Hyperfast pulsars as the remnants of massive stars ejected from ...
Recent proper motion and parallax measurements for the pulsar PSR B1508+55 indicate a transverse velocity of ∼1100 km s−1, which exceeds earlier measurements ...
[90]
Proper motions of six pulsars - Nature
Nov 20, 1975 · Tademaru and Harrison3 proposed a mechanism for the acceleration of pulsars, by asymmetric radiation, to velocities of 1,000 km s−1. Only rather ...
[91]
Probing the Low-velocity Regime of Nonradiative Shocks with ...
Nov 4, 2024 · While supernova remnants are the most well-studied example, neutron star (NS) bow shocks are also nonradiative and Balmer dominated. NS bow ...
[92]
The X-Ray Structure of the Pulsar Bow Shock G189.22 ... - NASA ADS
We present a deep observation with the Chandra X-Ray Observatory of the neutron star bow shock G189.22+2.90 in the supernova remnant (SNR) IC 443.
[93]
Researchers Waited 15 Years to Measure a Neutron Star's Journey
Aug 1, 2025 · With a measurement 15 years in the making, astronomers have pinned down the path of a neutron star launched by a collapsing star.
[94]
Researchers detect first 'heartbeat' of a newborn neutron star in ...
Sep 22, 2025 · The signal, spinning nearly a thousand times per second, points to the birth of a "millisecond magnetar"—a rapidly rotating neutron star with an ...Missing: young isolated
[95]
Chiming In on the Pulsar Search: Searching for Pulsars with CHIME
Sep 12, 2025 · While more than 3,700 pulsars have been discovered to date, increasing this number ... pulsars near the dense stellar population in the Milky ...
[96]
Chapter 6 Pulsars
E=2π2IP2≈2π2⋅1045gcm2(0.033s)2≈1.8×1049ergs. As magnetic dipole radiation extracts rotational energy, it slowly increases the period of a pulsar: ...Missing: torque | Show results with:torque<|separator|>
[97]
Binary and Millisecond Pulsars at the New Millennium - PMC
A very important additional difference between normal and millisecond pulsars is the presence of an orbiting companion. Orbital companions are much more ...
[98]
Parkes Multibeam Pulsar Survey – II. Discovery and timing of 120 ...
The survey aims to detect a large sample of pulsars for population studies (more than 600 have been discovered so far) as well as find individual objects of ...Missing: CHIME | Show results with:CHIME
[99]
Study of 72 Pulsars Discovered in the PALFA Survey - IOP Science
Jan 19, 2022 · We present new discoveries and results from long-term timing of 72 pulsars discovered in the Pulsar Arecibo L-band Feed Array (PALFA) survey.
[100]
First Discovery of New Pulsars and RRATs with CHIME/FRB
Nov 18, 2021 · We report the discovery of seven new Galactic pulsars with the Canadian Hydrogen Intensity Mapping Experiment's Fast Radio Burst (CHIME/FRB) backend.
[101]
Dependence of pulsar death line on the equation of state
Aug 17, 2017 · Pulsar death line can be defined in a $P\text{--}\dot{P}$ diagram. Traditionally, radio-loud pulsars are supposed to locate above the death
[102]
Regimes of Pulsar Pair Formation and Particle Energetics
The standard definition of a pulsar death line implies that pulsar radio emission turns off if the energetics of accelerated particles drops below the minimum ...Missing: Pdot | Show results with:Pdot
[103]
[1507.02924] Magnetars: the physics behind observations - arXiv
Jul 10, 2015 · Magnetars are the strongest magnets in the present universe and the combination of extreme magnetic field, gravity and density makes them unique laboratories ...
[104]
[1503.06313] Magnetars: Properties, Origin and Evolution - arXiv
Mar 21, 2015 · Magnetars are neutron stars in which a strong magnetic field is the main energy source. About two dozens of magnetars, plus several candidates, are currently ...
[105]
Hubble spots a magnetar zipping through the Milky Way - Phys.org
Apr 30, 2025 · Unraveling the track of the traveling magnetar. There are only about 30 known magnetars in the Milky Way galaxy. These dense balls of neutrons ...
[106]
Evidence for quark-matter cores in massive neutron stars - Nature
Jun 1, 2020 · Massive neutron stars are predicted to have sizable quark-matter cores. This finding has important implications for the phenomenology of neutron stars.<|separator|>
[107]
https://phys.org/news/2025-04-hubble-magnetar-milky.html
[108]
NICER Study of Pulsed Thermal X-Rays from Calvera: A Neutron ...
Calvera (1RXS J141256.0+792204) is an isolated neutron star detected only through its thermal X-ray emission. Its location at high Galactic latitude (b = +37°) ...<|separator|>
[109]
[1711.07210] Central Compact Objects in Supernova Remnants
Nov 20, 2017 · Central Compact Objects (CCOs) are a handful of sources located close to the geometrical center of young supernova remnants.
[110]
[2106.00700] A kinematic study of central compact objects and their ...
Jun 1, 2021 · Central compact objects (CCOs) are a peculiar class of neutron stars, primarily encountered close to the center of young supernova remnants (SNRs).
[111]
https://arxiv.org/abs/1711.07210
[112]
1966ARA&A...4..393W Page 393
No readable text found in the HTML.<|control11|><|separator|>
[113]
[PDF] Antony Hewish - Nobel Lecture
The first signals from CP 1919. But the source also exhibited apparent shifts of right ascension of up to 90 seconds which was evidence against a celestial ...
[114]
Multi-messenger Observations of a Binary Neutron Star Merger
Oct 16, 2017 · The Astrophysical Journal Letters, Volume 848, Number 2 Focus on the Electromagnetic Counterpart of the Neutron Star Binary Merger GW170817
[115]
LIGO-Virgo-KAGRA (LVK) Collaboration Detected a Remarkable ...
Apr 5, 2024 · The LIGO Livingston detector observed a gravitational-wave signal from the collision of what is most likely a neutron star with a compact object that is 2.5 to ...Missing: NS- NS
[116]
[1410.0886] 45 Years of Rotation of the Crab Pulsar - arXiv
Oct 3, 2014 · Since 1968, the rotation rate has decreased by about 0.5\,Hz, interrupted only by sporadic and small spin up events (glitches). 24 of these ...Missing: age optical TeV emission
[117]
First results of the glitching pulsar monitoring programme at the ...
Although erratic, this pulsar exhibits major glitches every 2–3 years. On ... 5(b). Timing analysis of Vela's glitch. Top: variations in the rotational frequency ...
[118]
Observations of 20 Millisecond Pulsars in 47 Tucanae ... - IOP Science
Five of the pulsars now known in 47 Tucanae have orbital periods of a few hours and implied companion masses of only ~0.03 M☉. Two of these are eclipsed at some ...
[119]
[1706.04908] Long-term observations of the pulsars in 47 Tucanae
Jun 15, 2017 · This paper is the second in a series where we report the results of the long-term timing of the millisecond pulsars (MSPs) in 47 Tucanae with the Parkes 64-m ...
[120]
GWTC-4.0: Updated Gravitational-Wave Catalog Released | LIGO Lab
Aug 26, 2025 · August 26th 2025 marks an important date for the international LIGO, Virgo, and KAGRA (LVK) gravitational-wave network. The collaborations are ...
[121]
https://iopscience.iop.org/article/10.1086/308859
[122]
Real-time inference for binary neutron star mergers using ... - Nature
Mar 5, 2025 · Mergers of binary neutron stars emit signals in both the gravitational-wave (GW) and electromagnetic spectra.
[123]
AI finds merging neutron stars in real time - Max-Planck-Gesellschaft
March 05, 2025. Astronomy Astrophysics Gravitational Waves. When two neutron stars merge, gravitational waves propagate into space. Shortly after this ...
[124]
Eccentric 'Star' Defies Easy Explanation, NASA's Chandra Finds
May 28, 2025 · Scientists have discovered a star behaving like no other seen before, giving fresh clues about the origin of a new class of mysterious objects.
[125]
Eccentric 'Star' Defies Easy Explanation, NASA's Chandra Finds ...
May 28, 2025 · This is thousands of times longer than pulsars, which are rapidly spinning neutron stars that have repeated variations multiple times a second.
[126]
XRISM observes windy neutron star, Chandra studies a luminous ...
Sep 21, 2025 · 25, 2024, it aimed at a neutron star known as GX13+1, located ... A few days before XRISM observations, GX13+1 suddenly got much brighter.Missing: magnetic field
[127]
Stratified wind from a super-Eddington X-ray binary is slower than ...
Sep 17, 2025 · Here we show the XRISM Resolve spectrum of the galactic neutron star X-ray binary, GX 13+1, which reveals one of the densest winds ever seen in ...
[128]
Breakthrough simulations shed light on the lightest neutron star ...
Feb 22, 2025 · Scientists have taken a major leap forward in understanding the formation of neutron stars, thanks to new 3D supernova simulations conducted ...