Transient discs around isolated accreting neutron stars

Anton V. Biryukov; Arthur G. Suvorov; Marina D. Afonina; Sergei B. Popov

REVIEW

Transient discs around isolated accreting neutron stars

Marina D. Afonina, Anton V. Biryukov, Sergei B. Popov, Arthur G. Suvorov

High Energy Astrophysical Phenomena · Submitted 2 Sep 2026

Pith's verdict: A transparent first census of disc-accreting isolated neutron stars whose headline numbers rest on a disc-formation assumption the authors themselves flag as unproven.

T0 · AI review

How this review was made: An AI referee (deepseek/deepseek-v4-flash) read the paper on September 8, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

T0 AI review An AI read the paper and answered a fixed list of questions.
T1 Audit Automatic checks that the references and files the paper relies on exist and match.
T2 Compute The paper's code or calculations re-run with the same results.
T3 Formal statement The claims written out in a formal math language a computer can read.
T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 3 major objections · 5 minor · Ask about this paper · The compiled record · Read the full review →

Challenge this review Comment (0) Read on arXiv ↗

Cite ▾

arXiv ID2609.02442 (version 2)

Pith IDZR26IOJO, Pith's permanent record of this paper and its review

CitationAnton V. Biryukov, Arthur G. Suvorov, Marina D. Afonina, Sergei B. Popov. "Transient discs around isolated accreting neutron stars." arXiv:2609.02442 (2026). https://pith.science/paper/ZR26IOJO

Cite this reviewBibTeX for Pith's review of this paper

Claim paper

Why Pith says this

Highlight any text to comment on it.

Editor's note

read the letter

You should know two things. First, this is the first population synthesis that counts disc-accreting isolated neutron stars in the Galaxy, and it is done transparently: 10^7 trajectories, two ISM models, two field-evolution models, four propeller models, and sampling errors from 95 subsamples. The new material is the torque-ratio criterion n_d = j_t/j_K (Eq. 16), the dAINS fractions and flux counts (Table 4, Fig. 2), and the ECME quenching and spectrum conditions (Eqs. 35 and 45). The qualitative result — dAINSs are a small, low-velocity, near-plane minority of old accretors — is robust across every propeller model that yields accretors at all. Credit where due: the LPT speculation is cleanly separated from the synthesis, the caveats are stated rather than buried, and nothing is fitted to the absent X-ray detections, so the central derivation is not circular.

The soft spots are real, and the paper names most of them. The stress-test concern lands: in Sec. 2 they note that the cooling time is of order the torque variation time and that convection in settling accretion can carry away angular momentum, then proceed to assume j_t > j_K is sufficient for a disc. If a thin disc cannot actually form, the dAINS population could be near zero even with Kolmogorov turbulence; the MHD alternative already cuts N_D by one to two orders of magnitude and eliminates it entirely for constant fields. So the up-to-2% figure is an upper envelope, not a robust prediction. A second distinct weakness: the abstract's headline few hundred sources at F_X >= 10^-14 erg/s/cm2 is never checked against the null X-ray searches cited in their own introduction, which is the obvious test of the model. Minor points: the R_A prefactor uncertainty is mentioned but not propagated, disc luminosity is excluded from the flux estimates, and the pitch-angle exponents l and m in the LPT spectrum are chosen so that l >= 3 reproduces GLEAM-X J1627 rather than predicted.

Bottom line: the quantitative claims are plausible but conditional on an unvalidated disc-formation criterion, and the paper would be stronger if it framed N_D as an upper limit until that criterion is tested. It still deserves a serious referee — the structure, honesty, and analytic results justify referee time. I would bring it to reading group and would cite it, with the caveat attached.

Referee report

Significance

If the central population numbers are robust, the paper would identify a new, potentially numerous class of hidden isolated neutron stars concentrated near the Galactic plane and in the low-velocity tail, and it would offer a quantitative route to explaining some long-period radio transients without invoking binaries or exotic spin-down. The study has real strengths: 10^7 trajectories, 95-subsample error estimates in Table 4, explicit comparison of constant and decaying fields, one- and two-phase ISM models, and a refreshingly honest discussion of the model limitations. The main weakness is that the entire population census is conditional on a disc-formation criterion that the manuscript itself acknowledges as unvalidated, so the predicted numbers should be read as scenario-dependent rather than as a firm count. The falsifiable spatial and velocity predictions are useful regardless of the exact normalization.

Major objections (3)

[Sec. 2, Eq. (16), Table 4] The disc-formation criterion n_d = j_t/j_K > 1 is load-bearing for every population number in Table 4 and Figures 1–2, but it is adopted despite the manuscript's own caveats in Sec. 2: the cooling time can be of the same order as the torque-variation time, and convection at the settling-accretion stage can carry away captured angular momentum and prevent disc formation. The sentence stating that disc formation is 'quite realistic' for the lowest-velocity objects accreting from dense ISM is not derived from any calculation and does not directly address the parameter combinations that dominate N_D. As written, the central prediction of up to ~2% of Galactic neutron stars as dAINSs is an upper-limit scenario rather than a demonstrated result; the authors should either supply a quantitative cooling and angular-momentum argument, or explicitly frame the population numbers as conditional on an idealized thin-disc assumption.
[Sec. 2 and Sec. 4, Table 5] Equation (16) relies on the extrapolation v_t(r) = 10 km/s (r/2×10^20 cm)^(1/3) from measured scales of about 100 AU and above down to Bondi radii ~0.1–10 AU, where the paper states there are no robust turbulence measurements. The manuscript's own MHD alternative, Eq. (23), reduces n_d by an order of magnitude and, according to Table 5, lowers N_D by factors of 10–100 and to zero for constant fields. The sampling uncertainties quoted in Table 4 are therefore not the dominant uncertainty in the central claim, and the abstract's 'few hundred' bright sources at F_X ≥ 10^-14 erg/s/cm^2 come from the optimistic branch. The authors should state explicitly that the Kolmogorov-normalization numbers are an upper bound unless independent calibration of the turbulent velocity at Bondi scales is provided.
[Sec. 2, Eqs. (19)–(22), Fig. 2] The X-ray fluxes and polar-cap temperatures used to count bright dAINSs are computed from the Bondi rate Mdot in Eqs. (19)–(22), but Sec. 2 states that in the settling-accretion regime 'the accretion rate onto the surface is sufficiently reduced,' and Sec. 6.1 repeats that settling accretion may influence disc properties. The manuscript never specifies whether disc formation removes the settling suppression and restores Bondi-rate surface accretion. If it does not, the surface luminosity and the numbers of sources above F_X = 10^-14 erg/s/cm^2 in Fig. 2 would be overestimated, potentially by a large factor. This point needs to be stated explicitly and, ideally, quantified.

Minor objections (5)

[Sec. 5.3, Eqs. (43)–(46)] The spectral prediction S(nu) ∝ nu^-(ell+3)/2 at nu >> nu_A is matched to GLEAM-X J1627 by allowing ell >= 3 and then choosing ell = 3, m = 2 for the near-cutoff slope; because ell and m are free, this is a consistency check rather than a prediction, and the text should label it as such.
[Sec. 5.2, Eqs. (26) and (29)] The cross-sectional area A(r) and the quenching condition assume a dipolar flux tube anchored at the Alfven radius, but the parameter h is never linked to the disc scale height or to the cooling arguments of Sec. 2; a sentence defining H and its expected range would improve reproducibility.
[Sec. 4, Fig. 3] The cumulative distribution is built from only seven LPT sources, and the text notes the poor statistics, but the figure would be more reproducible if the individual source latitudes and uncertainties were given in a table.
[Sec. 6.1] The sentence 'The question of low accretion rate the second group of questions' is missing a verb, and the Ferrario (2026) bibliographic entry appears incomplete.
[Abstract, Sec. 4] The abstract's 'up to ~2%' value is tied to the one-phase exponentially-decaying-field model, while the two-phase constant-field model gives 0.028%; a sentence indicating that the range spans more than an order of magnitude depending on the model would help readers.

How Pith checks papers → · Full audit of this paper →

The paper

The paper claims

Up to 6 million isolated neutron stars in the Milky Way may host transient accretion discs, with a few hundred bright enough to detect in X-rays.

The reading

This paper uses population synthesis to estimate how many isolated neutron stars can form an accretion disc from interstellar gas because of turbulence. The key finding is that disc formation requires very low relative velocities (below about 40 km/s) and weak magnetic fields, conditions satisfied by only a small fraction of the Galactic neutron star population. In the most optimistic scenario, with decaying magnetic fields and Kolmogorov turbulence, up to 2% of all neutron stars (about six million objects) may host a disc at some point. However, most are too dim to detect; only a few hundred are predicted to have X-ray fluxes above 10⁻¹⁴ erg/s/cm². The paper also suggests that these disc-accreting isolated neutron stars could be the engines behind some long-period radio transients through electron cyclotron maser emission.

What carries the argument

The central condition for disc formation is n_d = j_t / j_K > 1, where j_t is the specific turbulent angular momentum of captured ISM material at the Bondi radius and j_K is the Keplerian specific angular momentum at the magnetospheric (Alfvén) radius. This ratio depends on NS velocity, magnetic field, and ISM density as n_d ≈ 3.9 (v/10 km/s)^(-65/21) (B/10¹² G)^(-2/7) (n/1 cm⁻³)^(1/7) (Eq. 16). The turbulent velocity at the Bondi scale is extrapolated from larger scales using a Kolmogorov (α=1/3) or MHD (α=1/2) power law, which carries the most uncertainty.

Extended reading notes

Core claim

Through a large population synthesis tracking spin evolution, magnetic field decay, ISM density maps, and turbulence, the authors calculate that disc-accreting isolated neutron stars (dAINSs) can number up to ∼6×10⁶ in the Galaxy, corresponding to ∼2% of all isolated neutron stars, but only if the magnetic field decays exponentially (τ≈1.5 Gyr) and the turbulent velocity follows Kolmogorov scaling. Under constant-field models the number drops by a factor of ∼10, and under MHD turbulence (α=1/2) it drops by another factor of 10–100. The dAINSs are confined to low velocities (≲40 km/s for decaying fields, ≲20 km/s for constant fields) and lie within a few hundred parsecs of the Galactic plane. Their effective polar‑cap temperatures lie in the range 4×10⁶–10⁷ K, and a few hundred are expected at X‑ray fluxes >10⁻¹⁴ erg/s/cm². The paper further shows that the accretion rates and magnetic fields of these objects are compatible with the quenching condition for electron‑cyclotron maser emission, making them plausible progenitors of a subset of long‑period radio transients.

Load-bearing premise

The turbulent velocity of the ISM at scales comparable to the Bondi radius (0.1–10 AU) is extrapolated from measurements at scales ≥100 AU using a Kolmogorov power law, with no direct observational verification at sub‑AU scales; the MHD alternative (α=1/2) reduces disc numbers by factors of 10–100 or more.

Editorial extensions

If this is right

If dAINSs exist, they contribute a population of soft X‑ray sources (kT ∼ 0.3–0.9 keV) that could be detected by surveys like eROSITA, with a few hundred at fluxes above 10⁻¹⁴ erg/s/cm².
Disc accretion boosts the accretion efficiency compared to spherical settling, making dAINSs more luminous than previously considered AINSs and thus more likely to be observed.
The transient nature of the disc (disc lifetime comparable to the torque variation time ∼R_G/v) implies that individual sources may appear and disappear on timescales of months to years.
The spin periods of dAINSs can reach thousands of seconds, matching the range observed in long‑period radio transients (LPTs).
The radio emission mechanism (electron cyclotron maser) predicts circular polarisation at high frequencies and linear polarisation at low frequencies, a distinctive signature that can be tested with current radio telescopes.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

If the true ISM turbulent velocity at Bondi radii (0.1–10 AU) is an order of magnitude lower than the Kolmogorov extrapolation used here, then the number of dAINSs would be negligible, consistent with the complete absence of confirmed AINS candidates to date.
The paper's strongest quantitative prediction—the number of dAINSs at a given X‑ray flux—can be directly confronted with future deep X‑ray surveys (e.g., Athena, AXIS), providing a crisp observational test.
The LPT emission model developed in §5 suggests that dAINSs with decaying magnetic fields and thick discs (h∼1) are the most promising LPT progenitors; future multi‑frequency radio monitoring of candidate LPTs could distinguish between this scenario and white‑dwarf binary models.
The population synthesis framework could be extended to include a full treatment of cooling and convection in the settling‑accretion regime, which might further reduce the disc formation fraction and tighten the parameter space.

Original abstract

Mature, isolated neutron stars can accrete from the interstellar medium. Due to turbulence, the accreted gas can have a substantial angular momentum and form a disc around the compact object. In this paper, we perform a population synthesis of isolated neutron stars in the Milky Way that specifically tracks the possibility of disc formation, which typically requires a low spatial velocity of the compact object ($\lesssim 40$ km s$^{-1}$). In general, weak magnetic fields favour disc formation as the magnetosphere occupies a smaller volume. Still, even in the case of fields decaying exponentially over a characteristic timescale of $\sim 1.5$ Gyr, we find that for several realistic models of propeller spin-down, only a small fraction of accretors (down to $\sim 0.02$ %) attain discs. However, disc-accreting isolated neutron stars are relatively numerous among the brightest sources. We estimate that their number can reach a few hundred at X-ray fluxes of $\gtrsim 10^{-14}$ erg s$^{-1}$ cm$^{-2}$. We speculate that disc-accreting isolated neutron stars can manifest as long-period radio transient sources via electron cyclotron maser emission, predicting spectra and critical conditions for quenching.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13.0 via TypeSafe API, question bank abstract-p1-20260923, 2026-09-23: 1 request, 6,330 input tokens.

What the abstract itself states, as probabilities: claims first 0.05 · states limitation 0.62 · human participants 0.01 · controlled comparison 0.06 · randomized assignment 0.01 · states quantitative result 0.98 · compares with prior results 0.11 · mentions code or data availability 0.02.

The abstract, sentence by sentence

arXiv listing abstract v1. Results and the main claim are underlined; show every role.

Mature, isolated neutron stars can accrete from the interstellar medium. Due to turbulence, the accreted gas can have a substantial angular momentum and form a disc around the compact object. In this paper, we perform a population synthesis of isolated neutron stars in the Milky Way that specifically tracks the possibility of disc formation, which typically requires a low spatial velocity of the compact object ($\lesssim 40$ km s$^{-1}$). In general, weak magnetic fields favour disc formation as the magnetosphere occupies a smaller volume. Still, even in the case of fields decaying exponentially over a characteristic timescale of $\sim 1.5$ Gyr, we find that for several realistic models of propeller spin-down, only a small fraction of accretors (down to $\sim 0.02$ %) attain discs. However, disc-accreting isolated neutron stars are relatively numerous among the brightest sources. We estimate that their number can reach a few hundred at X-ray fluxes of $\gtrsim 10^{-14}$ erg s$^{-1}$ cm$^{-2}$. We speculate that disc-accreting isolated neutron stars can manifest as long-period radio transient sources via electron cyclotron maser emission, predicting spectra and critical conditions for quenching.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#	sentence	role	the authors' own work	states the main conclusion	strength
s0	Mature, isolated neutron stars can accrete from the interstellar medium.	background 1.00	0.04	0.02	asserted 0.75
s1	Due to turbulence, the accreted gas can have a substantial angular momentum and form a disc around the compact object.	background 1.00	0.08	0.03	hedged possibility 0.79
s2	In this paper, we perform a population synthesis of isolated neutron stars in the Milky Way that specifically tracks the possibility of disc formation, which typically requires a low spatial velocity of the compact object ($\lesssim 40$ km s$^{-1}$).	contribution 0.63	0.95	0.07	asserted 0.94
s3	In general, weak magnetic fields favour disc formation as the magnetosphere occupies a smaller volume.	background 0.65	0.17	0.04	asserted 0.82
s4	Still, even in the case of fields decaying exponentially over a characteristic timescale of $\sim 1.5$ Gyr, we find that for several realistic models of propeller spin-down, only a small fraction of accretors (down to $\sim 0.02$ %) attain discs.	result 1.00	0.95	0.30	asserted 1.00
s5	However, disc-accreting isolated neutron stars are relatively numerous among the brightest sources.	result 0.84	0.63	0.23	asserted 0.99
s6	We estimate that their number can reach a few hundred at X-ray fluxes of $\gtrsim 10^{-14}$ erg s$^{-1}$ cm$^{-2}$.	result 1.00	0.91	0.17	hedged possibility 0.77
s7	We speculate that disc-accreting isolated neutron stars can manifest as long-period radio transient sources via electron cyclotron maser emission, predicting spectra and critical conditions for quenching.	interpretation 0.78	0.59	0.20	speculative 1.00

What kind of paper the classifier reports this is: simulation study 0.95 · evidence from simulation or computation 1.00.

What was examined, and what was not

sources read

arXiv listing abstract v1

not examined

full text
figures and tables
references and citing sentences
relations to other papers

levels run

P abstract

not run

S sections
T body sentences
C citations
R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures

**Figure 1.** Distribution of accreting NSs with discs over six parameters: characteristic velocity

**Figure 2.** Number of dAINSs in the Milky Way with fluxes exceeding

**Figure 3.** Normalised cumulative density distribution of several observed LPTs

**Figure 4.** Maximum accretion rate from expression (35) for canonical choices

**Figure 5.** Fraction of accreting dAINSs such that the accretion rate lies below

**Figure 6.** Normalised frequency spectrum,

Figures reproduced from arXiv:2609.02442 by the authors.

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.