pith. machine review for the scientific record. sign in

arxiv: 2604.13557 · v1 · submitted 2026-04-15 · 🌌 astro-ph.HE

Recognition: unknown

Relaxation of magnetically-confined mountains on accreting neutron stars through cross-field mass transport

Ryan Brunet , Andrew Melatos , Pedro Rossetto
Authors on Pith no claims yet

Pith reviewed 2026-05-10 12:49 UTC · model grok-4.3

classification 🌌 astro-ph.HE
keywords accreting neutron starsmagnetic mountainshydromagnetic instabilitiescross-field mass transportGrad-Shafranov equilibriamass quadrupole momentKulsrud-Sunyaev recipe
0
0 comments X

The pith

Cross-field mass transport prevents instabilities from demolishing magnetic mountains on accreting neutron stars and preserves a nonzero quadrupole moment indefinitely without ohmic dissipation.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper examines how hydromagnetic instabilities alter the structure of magnetically confined mountains once accreted mass exceeds a critical value on neutron stars. Ideal magnetohydrodynamics breaks down locally, allowing mass to diffuse across field lines according to the Kulsrud-Sunyaev recipe. A numerical scheme evolves the system through successive axisymmetric Grad-Shafranov equilibria that incorporate this transport. The resulting mass-flux adjustments nullify the instabilities locally instead of allowing them to destroy the mountain. This leaves a stable nonzero mass quadrupole moment that persists indefinitely in the absence of ohmic dissipation.

Core claim

Hydromagnetic instabilities modify the structure of a magnetically confined mountain on an accreting neutron star once the accreted mass exceeds a critical value. Ideal magnetohydrodynamics and flux freezing break down, and mass diffuses across magnetic field lines locally wherever instabilities are excited. A self-consistent iterative numerical scheme evolves an axisymmetric magnetic mountain through a quasistatic sequence of Grad-Shafranov equilibria modified by instability-driven cross-field mass transport obeying the semi-analytic Kulsrud-Sunyaev recipe. The results are compared to an artificially stabilised mountain in which flux freezing does not break down and there is no cross-field

What carries the argument

The iterative numerical evolution of axisymmetric Grad-Shafranov equilibria that incorporates local cross-field mass transport driven by the Kulsrud-Sunyaev instability recipe, which permits the mass-flux distribution to adjust and suppress further instability growth.

If this is right

  • The mountain structure self-adjusts locally so that instabilities are nullified rather than amplified.
  • A nonzero mass quadrupole moment is maintained indefinitely without ohmic dissipation.
  • The outcome differs from the artificially stabilised case that enforces flux freezing and forbids transport.
  • The mass-flux distribution evolves to a new equilibrium that remains stable against further hydromagnetic disruption.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Persistent quadrupole moments could sustain continuous gravitational-wave emission from accreting neutron stars over long accretion timescales.
  • The same local transport adjustment might operate in other systems where magnetic confinement meets plasma instabilities.
  • Ohmic dissipation, rather than instability-driven demolition, would then set the ultimate lifetime of such mountains.

Load-bearing premise

The semi-analytic Kulsrud-Sunyaev recipe accurately gives the rate and location of instability-driven cross-field mass transport and the mountain can be evolved as a quasistatic sequence of axisymmetric equilibria.

What would settle it

A calculation or observation showing that the mountain's mass quadrupole moment falls to zero once accreted mass exceeds the critical value, or that the structure is fully demolished by instabilities even when the Kulsrud-Sunyaev transport is included.

Figures

Figures reproduced from arXiv: 2604.13557 by Andrew Melatos, Pedro Rossetto, Ryan Brunet.

Figure 1
Figure 1. Figure 1: Cross-field mass transport: adjusting the initial mass-flux distri￾bution (black curve) by moving mass Δ𝑀 (red rectangle) across the unsta￾ble flux surface 𝜓 ′ = 𝜓, resulting in the adjusted distribution (red curve) given by Equation (21). The curvature 𝑑 2𝑀/𝑑𝜓′2 increases at the leading (𝜓 ′ = 𝜓 + Δ𝜓) and trailing (𝜓 ′ = 𝜓 − Δ𝜓) boundaries, resulting in the formation of new, smaller unstable regions, whic… view at source ↗
Figure 2
Figure 2. Figure 2: Unstable region as a function of the initial accreted mass, with 𝑀 (init) a = 10−8𝑀⊙, 𝜌max = 5.0 × 1014 kg m−3 (top left); 𝑀 (init) a = 10−7𝑀⊙, 𝜌max = 4.9 × 1015 kg m−3 (top right); 𝑀 (init) a = 10−6𝑀⊙, 𝜌max = 4.0 × 1016 kg m−3 (bottom left); and 𝑀 (init) a = 10−5𝑀⊙, 𝜌max = 1.9 × 1017 kg m−3 (bottom right). Stability is determined via (5). In each panel, unstable and stable points are coloured orange and b… view at source ↗
Figure 3
Figure 3. Figure 3: Unstable fraction of the total flux tube volume within ten scale heights of the stellar surface, as a function of the initial accreted mass 𝑀 (init) a . Mass is distributed according to the mass flux distribution (19), for accreted masses 10−7𝑀⊙ ≤ 𝑀a ≤ 10−5𝑀⊙. The two curves are for 𝑏 = 3 (blue) and 𝑏 = 10 (orange). In these one-shot experiments (see Section 4.1), the unstable region grows monotonically wi… view at source ↗
Figure 4
Figure 4. Figure 4: Nullifying the instability through cross-field mass transport: iteration of an unstable region for 𝑀a = 5 × 10−7𝑀⊙, 𝑏 = 10, according to the numerical scheme proposed by Kulsrud & Sunyaev (2020) and described in Section 3 and Appendix B. The iteration number corresponding to each snapshot is quoted in the legend at the top right of each panel. The instability is nullified after 218 iterations. The format o… view at source ↗
Figure 5
Figure 5. Figure 5: Iterative progress towards nullifying the instability: fraction of points Í 𝑖 𝑁𝑖/𝑁total that are unstable, as a function of iteration (left axis); and the fraction of the total accreted mass Í 𝑖 Δ𝑀𝑖/𝑀a that is transported, as a function of iteration (right axis). The subscript 𝑖 labels flux surfaces: 𝑁𝑖 is the number of unstable points on flux surface 𝜓𝑖 , and 𝑁total is the total number of points on all fl… view at source ↗
Figure 7
Figure 7. Figure 7: Quasistatic assembly of the mountain: mass-flux distribution 𝑑𝑀˜ /𝑑𝜓˜ at intermediate iterations, after mass transport stabilises the moun￾tain for that iteration and before accreting further mass at the next iteration. The top panel shows the mass-flux distributions prior to normalisation by the total accreted mass, and includes the initial PM04 profile (blue crosses) cor￾responding to Equation (19) with … view at source ↗
Figure 6
Figure 6. Figure 6: Path dependence of the mass-flux distribution: 𝑑𝑀/𝑑𝜓 for two marginally stable mountains with 𝑀a = 5 × 10−7𝑀⊙ after cross-field mass transport. One mountain (blue curve) is built up in one shot, with 𝑀 (init) a = 5 × 10−7𝑀⊙. The other mountain (orange curve) is built up quasistatically, by accreting a total of 𝑀a = 5 × 10−7𝑀⊙ in ten equal increments, each 5 × 10−8𝑀⊙, and nullifying the unstable region by c… view at source ↗
Figure 8
Figure 8. Figure 8: Magnetic flux surfaces of a representative mountain with 𝑀a = 1 × 10−5𝑀⊙ before (dashed curves) and after (solid curves) cross-field mass transport. The top panel is for 𝑏 = 3, the bottom for 𝑏 = 10. The contour colour corresponds to the value 𝜓/𝜓∗. surface, md decreases from its pre-accretion value mi at 𝑟 = 𝑅∗ to an asymptotic, reduced (i.e. screened) value at 𝑟 ≫ 𝑅∗ (e.g. at the edge of the simulation v… view at source ↗
Figure 9
Figure 9. Figure 9: Hydromagnetic configuration for a mountain with cross-field mass transport and 𝑏 = 10, 𝑀a = 10−5𝑀⊙. Panels display the density 𝜌 (top left); pressure 𝑝 (top right); magnetic field strength 𝐵 (bottom left); and current density |j| = |∇ × B|/𝜇0 (bottom right), with 𝐵∗ = 108 T, 𝑅∗ = 104 m. Contours of 𝜓 are overlaid. MNRAS 000, 1–23 (2026) [PITH_FULL_IMAGE:figures/full_fig_p013_9.png] view at source ↗
Figure 11
Figure 11. Figure 11: Mass ellipticity 𝜖 as a function of accreted mass 𝑀a. Analytic curves for 𝑏 = 3 (solid) and 𝑏 = 10 (dashed) from Payne & Melatos (2004). Numerical results using 𝑑𝑀/𝑑𝜓 from equation (19) for 𝑏 = 3 (triangles), 𝑏 = 10 (diamonds), and numerical results using equation (23) for 𝑏 = 3 (points) and 𝑏 = 10 (crosses). & Melatos (2004) observed 𝜖 increasing monotonically with 𝑀a, they found that 𝜖 for 𝑏 = 3 (green … view at source ↗
Figure 10
Figure 10. Figure 10: Reduction of the magnetic dipole moment due to accretion. Top panel plots the normalised dipole moment |md |/|mi | as a function of altitude above the stellar surface, 𝑥˜ = (𝑟 −𝑅∗ )/𝑥0. Solid curves are the reduced dipole moment for the functional fit equation (23), for 𝑀a/𝑀⊙ = 10−5 , 10−4 , 10−3 , 10−2 , and 5×10−2 , for 𝑏 = 10. Dashed curve corresponds to 𝑀a/𝑀⊙ = 10−5 from Payne & Melatos (2004), also f… view at source ↗
read the original abstract

Hydromagnetic instabilities modify the structure of a magnetically confined mountain on an accreting neutron star, once the accreted mass exceeds a critical value. Ideal magnetohydrodynamics and flux freezing break down, and mass diffuses across magnetic field lines locally, wherever instabilities are excited. Here a self-consistent, iterative, numerical scheme is used to evolve an axisymmetric magnetic mountain through a quasistatic sequence of Grad-Shafranov equilibria as a function of the accreted mass, $M_{\rm a}$, modified by instability-driven cross-field mass transport obeying the semi-analytic, Kulsrud-Sunyaev recipe. The results are compared to an artificially stabilised mountain, in which flux freezing does not break down, and there is no cross-field mass transport. It is shown that cross-field mass transport prevents instabilities from demolishing the mountain. Instead, the mass-flux distribution adjusts locally to nullify the instabilities and preserve a nonzero mass quadrupole moment indefinitely in the absence of ohmic dissipation.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 1 minor

Summary. The paper claims that cross-field mass transport, modeled via the semi-analytic Kulsrud-Sunyaev recipe, allows magnetically confined mountains on accreting neutron stars to relax through a quasistatic sequence of axisymmetric Grad-Shafranov equilibria. Instability-driven transport adjusts the local mass-flux distribution to suppress further instabilities, preserving a nonzero mass quadrupole moment indefinitely without ohmic dissipation. This is shown by comparing the evolving mountain to an artificially stabilized case in which flux freezing is enforced and no cross-field transport occurs.

Significance. If the Kulsrud-Sunyaev prescription is physically accurate in this regime, the work provides a self-consistent mechanism for sustaining magnetic mountains against hydromagnetic instabilities. It has implications for neutron-star magnetic-field evolution, spin-down, and gravitational-wave emission from accreting systems. The forward iterative scheme (rather than fitting to a target quadrupole) is a methodological strength.

major comments (3)
  1. [Abstract / Methods] The central claim—that transport self-adjusts the mass distribution to nullify instabilities and sustain a nonzero quadrupole—rests entirely on the accuracy of the semi-analytic Kulsrud-Sunyaev recipe for local mass fluxes. The manuscript provides no validation, error analysis, or comparison to full MHD simulations to confirm that the recipe remains reliable in the nonlinear, strongly curved-field mountain equilibria (abstract and methods description).
  2. [Abstract] No numerical resolution, convergence tests, or error estimates are reported for the iterative Grad-Shafranov scheme or the instability-detection step that triggers transport. This is load-bearing because the quasistatic sequence and local adjustment mechanism depend on accurate detection and transport at each accreted-mass increment.
  3. [Results (comparison section)] The comparison run (artificially stabilized mountain with transport forbidden) demonstrates only the effect of forbidding transport; it does not test whether the chosen Kulsrud-Sunyaev law supplies the correct rates and locations. If the recipe over- or under-estimates transport, the iterative nullification mechanism may fail.
minor comments (1)
  1. [Abstract] The abstract is concise but could explicitly cite the original Kulsrud-Sunyaev reference for the transport recipe and state the value (or range) adopted for the critical accreted mass that triggers instability onset.

Simulated Author's Rebuttal

3 responses · 1 unresolved

We thank the referee for their detailed and constructive report. We address each major comment below and indicate the revisions made to the manuscript.

read point-by-point responses

  1. Referee: [Abstract / Methods] The central claim—that transport self-adjusts the mass distribution to nullify instabilities and sustain a nonzero quadrupole—rests entirely on the accuracy of the semi-analytic Kulsrud-Sunyaev recipe for local mass fluxes. The manuscript provides no validation, error analysis, or comparison to full MHD simulations to confirm that the recipe remains reliable in the nonlinear, strongly curved-field mountain equilibria (abstract and methods description).

    Authors: We agree that the central results depend on the applicability of the Kulsrud-Sunyaev prescription in the nonlinear regime of curved-field equilibria, and that the manuscript does not include new validation against full MHD simulations. The paper adopts this established semi-analytic recipe (originally derived for interchange instabilities) to explore the consequences of instability-driven transport within a quasistatic Grad-Shafranov framework. In the revised manuscript we have added an explicit limitations paragraph in the Methods section that states the assumptions of the recipe, cites its original derivation and prior uses in related contexts, and notes that direct MHD validation in strongly curved geometries remains an open task for future work. revision: partial

  2. Referee: [Abstract] No numerical resolution, convergence tests, or error estimates are reported for the iterative Grad-Shafranov scheme or the instability-detection step that triggers transport. This is load-bearing because the quasistatic sequence and local adjustment mechanism depend on accurate detection and transport at each accreted-mass increment.

    Authors: We acknowledge the omission of these numerical details. The revised manuscript now includes a dedicated numerical-methods subsection that specifies the grid resolution (128 radial by 256 angular zones), the convergence criterion for the Grad-Shafranov iteration (maximum residual < 10^{-8}), and the instability-detection threshold. We have performed and reported convergence tests by repeating the evolution at halved and doubled resolution; the final quadrupole moment changes by less than 4 percent. Error estimates on the quadrupole are derived from the GS residual and from the sensitivity to the precise instability threshold, and these are now quoted in the Results section. revision: yes

  3. Referee: [Results (comparison section)] The comparison run (artificially stabilized mountain with transport forbidden) demonstrates only the effect of forbidding transport; it does not test whether the chosen Kulsrud-Sunyaev law supplies the correct rates and locations. If the recipe over- or under-estimates transport, the iterative nullification mechanism may fail.

    Authors: The comparison is designed to isolate the qualitative effect of permitting cross-field transport versus enforcing strict flux freezing. We accept that it does not independently verify the transport rates or spatial distribution supplied by the Kulsrud-Sunyaev law. In the revised Results section we have clarified that the observed stabilization and sustained quadrupole are predictions of the model under the adopted transport prescription, and we explicitly cross-reference the limitations paragraph added in Methods. revision: partial

standing simulated objections not resolved
  • Direct validation of the Kulsrud-Sunyaev recipe against full nonlinear MHD simulations in the relevant curved-field, high-beta regime.

Circularity Check

0 steps flagged

No significant circularity; forward numerical evolution produces emergent outcome

full rationale

The paper evolves a sequence of axisymmetric Grad-Shafranov equilibria numerically, at each step applying cross-field mass transport according to the cited Kulsrud-Sunyaev semi-analytic recipe wherever instabilities are present. The central result—that the mass-flux distribution self-adjusts to suppress further instabilities and sustain a nonzero quadrupole—is an output of this iterative process, not imposed by construction, parameter fitting to the target quadrupole, or redefinition of inputs. No quoted step in the abstract or described method reduces the claimed prediction to an input quantity by definition or self-citation chain. The recipe is treated as an external physical prescription rather than derived within the paper, so the derivation remains self-contained against the stated assumptions.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The model rests on the Kulsrud-Sunyaev transport recipe (treated as given) and the assumption that ideal MHD breaks down only locally at instabilities. No new entities are introduced.

free parameters (1)
  • critical accreted mass for instability onset
    Determines when cross-field transport begins; value not specified in abstract.
axioms (2)
  • domain assumption Axisymmetric Grad-Shafranov equilibria remain valid throughout the quasistatic evolution
    Invoked to evolve the mountain as a sequence of equilibria.
  • domain assumption Kulsrud-Sunyaev recipe gives the correct local cross-field mass flux
    Used to modify the equilibria when instabilities are present.

pith-pipeline@v0.9.0 · 5475 in / 1374 out tokens · 33068 ms · 2026-05-10T12:49:06.998683+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

3 extracted references · 1 canonical work pages

  1. [1]

    OurNon-StableUniverse

    Aasi J., et al., 2015, Phys. Rev. D, 91, 062008 Abbott B., et al., 2007, Phys. Rev. D, 76, 082001 Abbott B. P., et al., 2017a, Phys. Rev. D, 95, 122003 Abbott B. P., et al., 2017b, ApJ, 847, 47 Abbott B. P., et al., 2019, Phys. Rev. D, 100, 122002 Abbott R., et al., 2020, ApJ, 902, L21 Abbott R., et al., 2022, Phys. Rev. D, 105, 022002 AnderssonN.,Glamped...

    work page doi:10.1007/978-3-319-22309-4 2015

  2. [2]

    Compared to the exponential form of𝑑𝑀/𝑑𝜓(19) used by Payne & Melatos (2004), the intermediate𝑑𝑀/𝑑𝜓profiles in Figure 7 are shallower. We quantify this in Figure C2 by plotting the mass contained in the polar cap,𝑀(𝜓≤𝜓a), as a function of the accreted mass as the mountain is assembled and comparing it to𝑀(𝜓≤ 𝜓a)without mass transport. In the former case, w...

    2004

  3. [3]

    The curves are fits of the form (23); the points are simulation outputs fitted by the curves

    The third (green points and curves) corresponds to the highest mass reached during the assembly sequence. The curves are fits of the form (23); the points are simulation outputs fitted by the curves. APPENDIX D: MOUNTAIN ON A NEWBORN NEUTRON STAR In this paper, we assume accretion occurs over long time-scales 𝑀a/ ¤𝑀a ≳10 8 yr, characteristic of low mass X...

    2011