Mountains on Neutron Stars: How Magnetic Fields and Accretion Shape the Terrain ✛

Original: "Modelling mountains on accreting magnetized neutron stars"
arXiv:2608.17508v1 · 2026-08-18 · CC BY 4.0 · 2 min
A model combines magnetic stresses, accretion heating, and crustal elasticity to predict deformations of neutron stars.

Context

This work belongs to the line of research on neutron star "mountains"—deformations that may explain why rapidly spinning pulsars in low-mass X-ray binaries do not spin up to their limit. Neighboring studies consider either magnetic deformations, thermal effects in the crust, or elasticity separately. Here, for the first time, all three factors are combined in a single self-consistent model.

Methods

The authors construct a neutron star model with a mass of 1.4 solar masses, consisting of a liquid core, an elastic crust, and an ocean. The magnetic field is assumed to be purely poloidal, and its effect is included via the Lorentz force. Crustal heating from accretion is calculated from heat transport equations, accounting for electron capture reactions. Deformations are found by solving the linearized elasticity equations in the crust with boundary conditions at the interfaces with liquid regions.

Results

A threshold accretion rate is found above which magnetic stresses and asymmetric heating produce deformations of opposite sign, and below which their roles are reversed. The maximum eccentricity without crustal failure reaches e ≈ 7×10⁻¹¹. The characteristic gravitational wave amplitudes lie below the sensitivity of LIGO, Virgo, and KAGRA, but may be accessible to the Einstein Telescope or Cosmic Explorer.

Implications

Unlike neighboring works where effects were considered in isolation, here it is shown that the magnetic field and thermal deformations can compete and even change the sign of the resulting deformation depending on the accretion rate. This is a first step toward a self-consistent description of "mountains" on accreting neutron stars.

Future development

The model can be extended by including more realistic magnetic field configurations (e.g., twisted-torus) and equations of state with "pasta" at the base of the crust. It is also necessary to couple the crustal calculations with the magnetohydrodynamics of the ocean to determine the accretion asymmetry parameter from first principles.

Impact

The results will aid in the search for continuous gravitational waves from neutron stars in low-mass X-ray binaries, especially in planning observations with third-generation detectors.

Next steps

Couple the crustal model with MHD models of the ocean and accretion flow to compute the parameter D without free assumptions. Also test the robustness of the results for different equations of state and magnetic field configurations.

Key open problems

This work is directly related to an unsolved problem in neutron star astrophysics: why the observed spin frequencies of accreting pulsars are limited to about 700 Hz, although theoretically they could spin up to kilohertz values. The model shows that gravitational radiation from "mountains" could be that brake, but only under certain conditions.

The knowledge machine's view

This section is addressed to the author of the paper and their colleagues. Neighbouring works were found by meaning-based search across our archive, not by word overlap. We do not review: we show the field around the work.

What the work does

The study presents a model of an accreting neutron star with a liquid core, an elastic crust, and an ocean, simultaneously incorporating the magnetic field, crustal heating from accretion, and elastic stresses. By solving the linearized equations of elasticity, the authors determine the deformations and the corresponding ellipticity, uncovering a threshold accretion rate at which the sign of the resulting deformation flips, and estimate the amplitudes of gravitational waves.

Where its strength is

The strength of the work lies in combining magnetic, thermal, and elastic effects in a single self-consistent model, which reveals their competition and the sign change of the deformation as a function of the accretion rate.

Where it could go next

Worth trying

  • Investigate how non-poloidal components of the magnetic field alter the deformations.
  • Assess the impact of crustal composition (e.g., the presence of impurities) on elastic properties and the threshold accretion rate.
  • Calculate the gravitational-wave signal waveform for the obtained deformations and compare it with the sensitivity of planned detectors.

Why this may matter

Understanding the mechanisms of mountain formation on accreting neutron stars is crucial for interpreting future gravitational-wave observations and for constraining the equation of state of matter at supra-nuclear densities.

Original abstract

as written by the authors arXiv:2608.17508v1

Continuous gravitational waves from accreting neutron stars in Low Mass X-ray Binaries are one of the main targets for current and next generation ground based detectors. In order to select the most promising astrophysical sources, however, reliable predictions for the signals are required, and it is therefore necessary to develop models that consistently account for the combined effects of magnetic stresses, accretion-induced heating, and the elastic response of the crust.}{We present a model for computing the quadrupolar deformation, incorporating for the first time the coupled effects of a poloidal magnetic field, deep crustal heating, and crustal elasticity. Perturbations to the star's structure driven by the Lorentz force density and by thermally-induced density variations are computed by solving a system of linearised deformation equations in the crust, for which we consider the full elastic response, while the ocean and core treated as barotropic fluids. We identify a threshold accretion rate whose value depends on crustal microphysics and the superfluid gaps in the core, above which magnetic stresses and asymmetric accretion drive deformations of opposite sign, while below this threshold their roles are reversed. The predicted eccentricities reach magnitudes up to $$\varepsilon\sim 10^{-11}$$, corresponding to characteristic gravitational-wave strains accessible to next-generation detectors such as the Einstein Telescope or Cosmic Explorer, but generally below the sensitivity of current LIGO, Virgo and KAGRA interferometers. These results are consistent with the non-detection of continuous gravitational waves from accreting neutron stars in Low Mass X-ray Binaries in recent observational campaigns, but highlight the need of reliable models to understand the impact of gravitational wave emission in these systems and select relevant targets for future searches.

The maximum deformation that a neutron star crust can sustain corresponds to a "mountain" less than a millimeter high on an object with a radius of 12 kilometers, like a bump on a perfectly polished billiard ball.

h_c = \frac{4\pi^2 G}{c^4 d} I_0 |e| \omega^2
The amplitude depends on the eccentricity, moment of inertia, and the square of the angular velocity

Key numbers

  • maximum eccentricity: 7×10⁻¹¹
  • threshold accretion rate: about 0.1 Eddington
  • stellar mass: 1.4 solar masses
  • stellar radius: 12.66 km
  • core-crust radius: 11.55 km
Scientists
Enrico FermiPaul DiracSubrahmanyan ChandrasekharHendrik LorentzIsidor Isaac RabiAlbert Einstein
Concepts
gravitational wave parameter estimation magnetar magnetar pulsations in grb accretion rate x ray binary neutron star mass limit x ray magnetic circular dichroism neutron detection techniques quadrupole radiation formula xray polarization binary stellar evolution gravitational wave detector nicer asteroseismology crust nanograv hubble constant measurement magnetorotational supernova xmm newton pulsar phenomena