Differentially Rotating Neutron Stars

Updated 22 January 2026

Differentially rotating neutron stars are compact remnants with angular velocity varying across the star, arising in core-collapse supernovae and mergers.
They feature diverse equilibrium configurations such as spheroidal and quasi-toroidal forms that influence maximum mass support and stability.
Magnetohydrodynamic processes and microphysical effects shape their evolution, impacting gravitational-wave emission and electromagnetic counterparts.

A differentially rotating neutron star is a compact stellar remnant in which the angular velocity varies nontrivially as a function of position, in contrast to uniform (rigid) rotation. Differential rotation naturally arises in the aftermath of core-collapse supernovae and binary neutron star mergers, where hydrodynamical and magnetic processes have not yet redistributed angular momentum. The resulting structures—often labeled "hypermassive" or "übermassive" when supported above the mass limit for uniform rotation—exhibit unique equilibrium configurations, stability properties, and oscillation spectra, profoundly impacting their gravitational-wave and electromagnetic signatures, as well as their secular evolution and collapse behavior.

1. Rotation Laws and Equilibrium Families

The angular velocity profile Ω is typically encoded through a "rotation law" specifying how Ω falls off with distance from the axis. The standard Komatsu-Eriguchi-Hachisu (KEH) "j-constant" law prescribes the specific angular momentum as $u^t u_\phi = A^2 (\Omega_c - \Omega)$ , so $\Omega(\varpi) = \Omega_c / \left[1 + (\varpi/A)^2 \right]$ , where $\Omega_c$ is the central angular velocity, $\varpi$ the cylindrical radius, and $A$ the differential rotation length scale (Espino et al., 2019, Gondek-Rosinska et al., 2016, Szewczyk et al., 2023, Tootle et al., 8 Jan 2026). Realistic merger remnants are better described by more general multi-parameter laws, such as the Uryū-type rotation profile, which allows off-axis maxima and more accurate post-merger fits (Staykov et al., 2023, Cipriani et al., 5 Dec 2025). The degree of differential rotation is often quantified by $\hat{A}^{-1} = R_e / A$ , with values $\hat{A}^{-1} \sim 0.7-1$ typical for merger remnants (Szewczyk et al., 2023).

Solution space is classified into several equilibrium families (Espino et al., 2019, Gondek-Rosinska et al., 2016):

Type A: Spheroidal branch, density maximum at center, extends continuously from the nonrotating limit to mass shedding.
Type B/C: Quasi-toroidal configuration, exhibiting off-axis density maxima and, at high rotation, low-density "funnels" along the axis. Such configurations can support the largest masses.
Type D: Pinched/toroidal, double mass-shedding sequences; typically astrophysically marginal.

Inclusion of additional physics—such as hyperons, Δ-isobars, dark matter, or scalar fields—alters the detailed equilibrium structure but generally preserves this taxonomy (Farrell et al., 2023, Cipriani et al., 5 Dec 2025, Staykov et al., 2023).

2. Maximum Mass, Stability, and Dynamical Evolution

The presence of differential rotation can increase the maximum allowable rest (baryonic) mass $M_{0,\text{max}}$ well beyond the uniformly rotating (Kepler) or nonrotating (TOV) limit. For realistic nuclear equations of state (EOS), the ratio $M_{0,\text{max}} / M_{0,\text{TOV}}$ can reach 2–2.5 for quasi-toroidal Type B/C models, and up to ~4 for moderately stiff polytropes with tailored rotation profiles (Espino et al., 2019, Gondek-Rosinska et al., 2016, Szewczyk et al., 2023). The maximum is realized for moderate differential rotation; if the shear becomes too extreme, the mass support may decrease or configurations become unstable.

Stability is assessed via the turning-point method: along equilibrium sequences at constant rest mass or angular momentum, the onset of dynamical or secular instability corresponds to extrema in the mass–central density curve (Szewczyk et al., 2023). Numerical time evolutions confirm that quasi-toroidal hypermassive stars above the rigid-rotation mass limit are dynamically stable against axisymmetric collapse until they cross this threshold (Szewczyk et al., 2023, Tootle et al., 8 Jan 2026).

Key physical processes affecting the secular evolution include:

Magnetorotational instability (MRI) and MHD turbulence: Efficiently convert differential rotation into magnetic field energy and drive angular momentum transport on timescales $\sim 10^1-10^2$ ms (Miravet-Tenés et al., 8 Sep 2025, Kiuchi et al., 2012, Siegel et al., 2014). This leads to redistribution of angular momentum, eventual collapse, and set the lifetime of hypermassive remnants.
Viscous effects: Modeled via effective "α-viscosity" in axisymmetric simulations, also drive angular momentum transport, heat the envelope/torus, and power neutrino-rich ejecta (Shibata et al., 2017).
Cosmic censorship and collapse end-states: Stars can have equilibrium configurations with $J/M^2 > 1$ (supra-Kerr), but dynamical collapse always results in sub-Kerr black holes; naked singularities do not form generically (Giacomazzo et al., 2011).

3. Oscillation Spectrum and Gravitational-wave Emission

Differentially rotating neutron stars exhibit a rich spectrum of fluid and spacetime modes. Fundamental (quasi-radial $l=0$ ; quadrupolar $l=2$ ) modes are prominent GW sources in the post-merger phase (Jaraba et al., 15 Jan 2026, Yip et al., 2024, Chirenti et al., 2013). 2D axisymmetric and full 3D dynamical simulations and perturbative calculations show that:

The frequencies of the $F$ -mode ( $l=0$ ) and $^2f$ -mode ( $l=2$ ) scale nearly linearly with stellar compactness and the kinetic-to-binding energy ratio $T/|W|$ , for a range of differential rotation measures (Yip et al., 2024, Jaraba et al., 15 Jan 2026, Chirenti et al., 2013).
Differential rotation typically suppresses GW emission from $f$ -modes by 5–15% and shifts mode frequencies by a few percent, with larger shifts at higher shear (Chirenti et al., 2013, Jaraba et al., 15 Jan 2026). The CFS instability thresholds can be crossed in prograde $m=2$ modes for sufficiently high $T/|W|$ .
GW-driven non-axisymmetric (bar-mode, low- $T/|W|$ ) instabilities are possible in strongly differentially rotating configurations, enhancing GW luminosity and observability (Tootle et al., 8 Jan 2026, Giacomazzo et al., 2011).

Universal relations for oscillation frequencies have been established, enabling inference of stellar compactness and degree of differential rotation from GW data, though the presence of strong differential rotation must be explicitly accounted for due to significant systematic frequency shifts (Yip et al., 2024).

4. Microphysics and Exotic Physics

Differential rotation modifies the role of microphysical composition, but certain trends are robust:

The fractional increase in maximum mass afforded by differential rotation is nearly independent of detailed composition—whether the EOS contains hyperons, Δ-isobars, or just nucleons, the enhancement remains at ~16–17% for fixed rotation-law parameters (Farrell et al., 2023).
Inclusion of a dark-matter core reduces the maximum mass at zero or low rotation but this effect is mitigated at high differential rotation; strong rotation "restores" mass support lost to the dark core (Cipriani et al., 5 Dec 2025).
In scalar-tensor gravity, strong-field scalarization dramatically increases the maximum mass and angular momentum supportable by differentially rotating stars, delays the transition to quasi-toroidal structure, and introduces additional branches and stability features (e.g., a second turning point), with observable consequences for GW signals and post-merger morphology (Doneva et al., 2018, Staykov et al., 2023).

Thermal non-barotropicity is also relevant in newborn or post-merger neutron stars. New formalism enables the construction of non-barotropic, stationary differentially rotating models—an essential advance for realistic remnant modeling (Camelio et al., 2019).

5. Magnetic Topology, Outflows, and Electromagnetic Counterparts

The interplay between differential rotation and magnetic topology underpins outflow and EM emission phenomenology:

Magnetically driven winds: Differential-rotation-driven winding amplifies the toroidal field $B_\phi \sim t$ , launching quasi-isotropic, baryon-loaded winds once magnetic pressure accumulates near the surface (Siegel et al., 2014, Kiuchi et al., 2012). The stationary EM luminosity $L_{EM} \propto B^2 R^3 \Omega$ matches the observed timescales and X-ray plateau luminosities in short-GRB afterglows.
Magnetic topology: Open field lines enforce rigid corotation via Alfvénic coupling, while closed flux tubes allow differential shear layers long-lived on viscous timescales, as shown in full-MHD studies (Anzuini et al., 2020). This impacts pre-merger and post-merger angular momentum distribution, and may influence jet formation and observed EM counterparts.

Magnetic and turbulent transport ultimately erase differential rotation, setting observational windows for hypermassive remnants and modulating kilonova yields via viscous and MHD-driven mass ejection (Shibata et al., 2017, Kiuchi et al., 2012).

6. Numerical Implementation and Universal Relations

Modeling differentially rotating neutron stars requires robust, high-accuracy equilibrium and evolution solvers. Recent advances include:

Spectral initial-data codes such as FUKA/KADATH, supporting multiple coordinate systems and rotation laws, validated against RNS and dynamical-evolution tests (Tootle et al., 8 Jan 2026).
Two-fluid formulations for baryon–dark matter systems, extended RNS codes, and handling of tabulated, finite-temperature EOS (Cipriani et al., 5 Dec 2025, Staykov et al., 2023).
Well-balanced hydrodynamics and metric-primitive evolution to preserve equilibrium and enable precise mode extraction (Jaraba et al., 15 Jan 2026).

Universal "four-hair" relations, established in the weak-field limit, provide quasi–EoS-independent links between higher multipole moments and the first four (mass monopole, angular momentum, quadrupole, and current-octupole), robust also to mild differential rotation (Bretz et al., 2015). These greatly aid in constraining neutron star parameters from multimessenger data.

In summary, differentially rotating neutron stars represent a multi-faceted regime of relativistic astrophysics, combining intricate rotation profiles, microphysics, strong-field gravity, magnetohydrodynamics, and dynamic instability phenomena. Their modeling, supported by increasingly sophisticated numerical and theoretical frameworks, is central to interpreting multi-messenger signals from core-collapse and merger environments, probing the physics of dense matter, and testing the strong-field regime of gravity (Gondek-Rosinska et al., 2016, Espino et al., 2019, Szewczyk et al., 2023, Yip et al., 2024, Tootle et al., 8 Jan 2026, Siegel et al., 2014, Jaraba et al., 15 Jan 2026, Doneva et al., 2018).