⚛️ general relativity

Modelling mountains on accreting magnetized neutron stars

This paper presents a novel model that simultaneously accounts for magnetic stresses, deep crustal heating, and elastic responses to predict quadrupolar deformations on accreting magnetized neutron stars, revealing a critical accretion rate threshold that determines the sign of deformation and yielding gravitational wave strain estimates ( $\varepsilon\sim 10^{-11}$ ) consistent with current non-detections but potentially accessible to next-generation detectors.

Original authors: T. Brusco, B. Haskell, M. Razzano, M. Bejger, J. L. Zdunik

Published 2026-08-19

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: T. Brusco, B. Haskell, M. Razzano, M. Bejger, J. L. Zdunik

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Deep in the cosmos, hidden within the crushing gravity of dead stars, lies a mystery that could help us hear the universe in a new way. These dead stars, known as neutron stars, are the incredibly dense remnants of massive stars that have exploded. They are so heavy that a single teaspoon of their material would weigh as much as a mountain on Earth. Some of these stars are born with powerful magnetic fields and are constantly fed by a nearby companion star, a process called accretion. As this new material falls onto the surface, it heats up and can create uneven bumps, or "mountains," on the star's crust. If a neutron star spins while carrying such a bump, it should ripple the fabric of space and time, sending out continuous gravitational waves. These waves are faint vibrations that stretch and squeeze the universe itself, and detecting them would allow scientists to peer inside the star's core, revealing how matter behaves under conditions impossible to recreate on Earth. However, despite years of searching, these specific signals have remained elusive, leaving astronomers to wonder if the mountains are too small to see or if their models of how these stars work are missing a crucial piece of the puzzle.

A team of researchers has now built a new, more complete model to solve this puzzle, combining three forces that were previously studied in isolation. They focused on how the magnetic field, the heat from falling material, and the stiffness of the star's outer shell interact to shape these mountains. In their simulation, they treated the star's core as a fluid, the outer layer as a solid elastic crust, and the thin ocean of material on top as a fluid that conducts heat. They calculated how the magnetic field pushes on the star, how the heat from accretion expands the crust, and how the crust's own strength resists these changes. By solving the equations that govern these forces, they were able to predict the exact shape and size of the deformations that would form on the surface of a spinning, accreting neutron star.

The researchers discovered a surprising threshold that determines the shape of these mountains. They found that the outcome depends heavily on how fast the star is eating material from its companion. When the rate of accretion is high, the magnetic field and the heat from the falling matter work against each other, creating a specific type of bump. But when the accretion rate drops below a certain point, the roles reverse: the magnetic field and the thermal effects switch their influence, leading to a different shape. This means that the same star could look very different depending on how much material it is currently swallowing. The study also revealed that the internal state of the star's core, specifically whether the particles inside are in a superfluid state, changes the exact point where this switch happens. This suggests that the internal physics of the star plays a direct role in how its surface deforms.

The size of the mountains predicted by this new model is incredibly small, yet significant. The researchers calculated that the deformations could reach a level of about one part in ten billion. While this sounds tiny, it is large enough to potentially be detected by the next generation of gravitational wave detectors, such as the Einstein Telescope or Cosmic Explorer, which are currently being planned. These future instruments will be far more sensitive than the current detectors like LIGO and Virgo, which have not yet found these signals. The fact that the predicted signals are just below the reach of today's machines explains why we have not heard them yet, but it also offers hope that they are within our grasp in the near future.

Crucially, the study also identified a limit to how large these mountains can grow. If the magnetic field is too strong or the accretion is too uneven, the stress on the star's crust becomes so great that the solid shell would crack or flow like plastic, destroying the mountain. The researchers found that for most realistic scenarios, the crust remains intact, but only if the unevenness of the falling material is kept within a narrow range. This finding helps astronomers understand why we might not see signals from every accreting neutron star; some might be too smooth, while others might have their crusts broken by the very forces trying to build mountains.

Ultimately, this work provides a more realistic map of what these cosmic objects look like. By showing that magnetic fields, heat, and elasticity must be considered together, the study offers a clearer path for future searches. It suggests that the silence from current detectors is not a failure, but a reflection of the subtle balance of forces at play. As we build better tools to listen to the universe, this model will help scientists know exactly what to listen for, turning the faintest whispers of space-time into a clear voice from the heart of a neutron star.

Technical Summary: Modelling Mountains on Accreting Magnetized Neutron Stars

Problem Statement
Continuous gravitational waves (GWs) from rotating, non-axisymmetric neutron stars (NSs) in Low Mass X-ray Binaries (LMXBs) are primary targets for current and next-generation detectors. While the "torque balance" scenario provides an upper limit on GW emission by equating spin-up accretion torques with spin-down GW torques, it does not predict the specific physical deformations ("mountains") required to generate such signals. Existing literature has treated magnetic deformations, thermal effects from accretion, and crustal elasticity in isolation. However, a self-consistent model accounting for the coupled effects of poloidal magnetic fields, deep crustal heating, and the elastic response of the crust has been lacking. Such a model is necessary to reliably predict signal amplitudes and understand the interplay of forces that may lead to subtle cancellations or shape changes in the star.

Methodology
The authors present a perturbative model to compute the quadrupolar deformation of a magnetized, accreting NS with an elastic crust. The approach involves solving a system of linearized deformation equations around a spherical background in hydrostatic and thermal equilibrium.

Background Model: The star is modeled with a total mass of $1.4 M_\odot$ and an accreted BSk21 equation of state (EoS). The structure includes a fluid core (treated as a barotropic fluid with a polytropic index $\gamma=2$ ), a solid elastic crust, and a fluid ocean. Thermal structure is determined by solving coupled heat transport equations for shallow crustal heating (SCH), deep crustal heating (DCH), and neutrino cooling, considering both normal ( $\Delta=0$ ) and superfluid ( $\Delta=1$ MeV) core scenarios.
Magnetic Field: The field is assumed to be purely poloidal and static. The magnetic function is derived by connecting analytical solutions of the Grad-Shafranov equation in the fluid core and a modified equation in the solid crust, ensuring continuity of the field and the Lorentz force density at the crust-core interface.
Perturbation Framework:
- Core/Ocean: Treated as fluids where density perturbations are driven by the Lorentz force and barotropic EoS.
- Crust: The crust is treated as a solid lattice. The authors solve for the displacement vector and traction vector by incorporating the full elastic response (shear modulus) and non-barotropic effects (temperature-dependent composition changes due to electron capture reactions).
- Thermal-Magnetic Coupling: The model accounts for the magnetic field's influence on thermal conductivity (via the magnetization parameter $x_m$ ) and the asymmetric deposition of accreted matter on polar caps, which creates temperature asymmetries at the crust-ocean interface.
Numerical Solution: The linear differential systems for temperature perturbations and crustal deformations are solved numerically using the Radau IIA method, subject to boundary conditions at the crust-core and crust-ocean interfaces. The validity of the results is constrained by the Von Mises strain criterion ( $\sigma < 0.1$ ) to ensure the crust remains in the elastic regime.

Key Contributions

Unified Framework: This is the first model to simultaneously incorporate the coupled effects of a poloidal magnetic field, deep crustal heating, and crustal elasticity in a self-consistent manner.
Threshold Accretion Rate: The study identifies a critical threshold accretion rate, $\dot{M}_{th}$ $\dot{M}_{t h}$ , which depends on crustal microphysics and the superfluid energy gap in the core.
- Above $\dot{M}_{th}$ : Magnetic stresses and asymmetric accretion drive deformations of opposite signs (magnetic fields favor oblate shapes, while accretion favors prolate shapes).
- Below $\dot{M}_{th}$ : The roles are reversed; magnetic stresses and asymmetric accretion drive deformations of opposite signs compared to the high-accretion regime. Specifically, the paper notes that while magnetic fields typically produce oblate deformations, at low accretion rates, the indirect effect of the magnetic field on temperature perturbations dominates over the direct Lorentz force, leading to prolate configurations under strong magnetic fields.
Inverted Behavior: The authors report an unexpected inversion where, at low accretion rates, the indirect effect of the magnetic field on temperature perturbations (via the source term in heat equations) appears to dominate over the direct Lorentz force, leading to prolate configurations under strong magnetic fields.

Results

Eccentricity Magnitudes: The predicted eccentricities ( $e$ ) reach magnitudes up to $\sim 10^{-11}$ without exceeding the crustal breaking strain.
Dependence on Parameters: The eccentricity is found to be linearly dependent on the surface temperature perturbation. The slope of this relationship changes sign at $\dot{M}_{th}$ . For a normal core, $\dot{M}_{th} \approx (1\text{--}2) \times 10^{-9} M_\odot \text{yr}^{-1}$ ; for a superfluid core, $\dot{M}_{th} \approx (2\text{--}6) \times 10^{-9} M_\odot \text{yr}^{-1}$ .
Gravitational Wave Strain: The corresponding characteristic GW strains are generally below the sensitivity of current detectors (LIGO, Virgo, KAGRA) but fall within the projected sensitivity range of third-generation detectors like the Einstein Telescope (ET) or Cosmic Explorer.
Crustal Stability: The model reveals that for high magnetic fields ( $B_p \ge 10^{11}$ G) and high accretion rates, the crust is prone to breaking unless the accretion asymmetry parameter $D$ is very small. Conversely, at lower accretion rates, the crust can sustain larger asymmetries.

Significance and Claims
The paper claims that these results are consistent with the non-detection of continuous GWs from LMXBs in recent observational campaigns, as the predicted strains are below current sensitivity thresholds. However, the authors emphasize that reliable models are essential for selecting relevant targets for future searches with next-generation detectors.

The study highlights that magnetic stresses, thermal asymmetries, and crustal elasticity cannot be treated independently. The existence of a threshold accretion rate where the dominant deformation mechanism inverts suggests a complex magneto-thermal interplay that must be understood to interpret future GW detections or place meaningful constraints on the internal physics of neutron stars. The authors note that while the current model provides a flexible framework, future work is required to couple these crustal calculations with magnetohydrodynamic models of the ocean to determine the accretion asymmetry parameter $D$ from first principles.

More like this