{"id":"e05f8251-d980-4635-bad2-be8f3a04c1c6","arxiv_id":"2412.15464","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"A population synthesis model with detailed neutron star spin evolution reproduces the observed Galactic binary neutron star population and explains why heavy systems like GW190425 are rarely seen in radio surveys.","lead":"Radio surveys of our Galaxy find neutron star pairs with total masses near 2.6-2.7 solar masses, while gravitational wave detectors saw one much heavier merger. This modeling paper proposes that the missing heavy systems are simply radio-quiet: their neutron stars keep strong magnetic fields and stop broadcasting quickly, so radio and gravitational wave surveys naturally see different mass distributions.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Mass-dependent radio-lifetime result hinges on unvaried §2.2.3 assumptions on magnetar initial fields and accretion-driven decay; relaxing them may erase the predicted radio/GW mass split.","rationale":"The reader's weakest assumption is the same one I regard as load-bearing: the assignment of extremely strong initial fields to NSs from massive progenitors and the accretion-driven decay timescale in §2.2.3. My analysis sharpens the failure mode: the important direction is not 'accretion fails to bury the field' (which would leave higher B and hence shorter radio lifetimes, reinforcing the conclusion), but rather (a) lower B_ini for massive progenitors or (b) a shorter td_acc, both of which would leave the recycled first-born NS in heavy BNSs with B≲1e11-1e12 G, pushing P below P_death and extending radio lifetimes to ≳100 Myr. Under those variants the ~7.5% radio-detectable heavy fraction in the fiducial model would rise by a large factor, so the model's ability to explain the complete absence of Mtot>3 Msun Galactic BNSs would fail. The paper's sensitivity analysis (Section 4) is substantial — wind slope β, subsonic propeller vs settling, and EOS effects are all checked — but none of these touches B_ini or td_acc, the two parameters that set the mass dependence. This is a genuine gap because the magnetar fraction among core-collapse SNe is not observationally established as 100% for M≳20 Msun, and the accretion field decay timescale is known only to within an order of magnitude. For that reason, the central claim should remain CONDITIONAL until the proposed reruns are performed. The paper deserves credit for a detailed spin-evolution implementation, explicit death-line treatment, and a falsifiable prediction of 19-22% heavy merging fraction; the concern is not an internal inconsistency but an uncalibrated sensitivity at the pivot of the argument. My recommended verdict is therefore UNCHANGED (CONDITIONAL), with the condition that the B-field branch be tested. Agreement with the reader: agree, with the clarification that the 'not bury' phrase should read as 'too-fast decay' in the dangerous direction.","tokens_in":28697,"tokens_out":10063,"duration_ms":85513,"concrete_test":"Re-run the fiducial population synthesis with: (A) B_ini drawn for all NSs from the same log-normal (mean log10 B=13, sc=0.55) instead of 1e14-1e15 G for M_ZAMS≥20 Msun; (B) fiducial B_ini but accretion-decay timescale td=1e5 yr and 1e7 yr instead of 1e6 yr. For each, report (1) the fraction of survived radio-detectable BNSs with Mtot≥3 Msun (fiducial 7.5%; observed 0/19) and (2) the z=0 merging BNS fraction with Mtot≥3 Msun (fiducial 19-22%). If (1) rises above ~10% in any variant, or (2) changes by more than a factor of two, the claimed radio/GW mass split is contingent on those §2.2.3 assumptions.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central selection-effect claim requires that the first-born NSs in high-Mtot systems end their HMXB stage with B ≳ 1e12 G, so their dipole spin-down lifetime is short (cases B and D; Figs. 3 and 5). That outcome is produced by two coupled assumptions in §2.2.3: (i) NSs from ZAMS stars with M≥20 Msun are born with B_ini uniform in [1e14,1e15] G, rather than the log-normal (mean 1e13 G) used for lower-mass stars; and (ii) during accretion the magnetic field decays with td=1e6 yr, so in the short (~4 Myr) HMXB phase of massive binaries B drops only to ~1e13 G. If (i) is relaxed to the same log-normal for all NSs, the same td gives B_final~2e11 G after 4 Myr of accretion; the equilibrium spin period from RLOF is then ~0.1 s and the death-line crossing time exceeds ~100 Myr, making heavy BNSs radio-detectable for a substantial fraction of a Hubble time. Then the observed absence of Galactic BNSs with Mtot>3 Msun in radio surveys would not be reproduced by the model, and the claimed 19-22% merging fraction of heavy BNSs would be decoupled from the radio selection argument. The paper's Section 4 sensitivity tests vary wind prescription, subsonic accretion model, and EOS, but do not vary B_ini for massive progenitors or td_acc, even though these are the parameters that set the mass dependence of the radio lifetime.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents a binary population synthesis model, based on a modified BSE code, that follows the spin evolution of the first-born neutron star through the high-mass X-ray binary phase and uses the resulting radio lifetimes to construct the present-day population of Galactic binary neutron stars with pulsar components. The authors argue that low-total-mass BNSs (Mtot around 2.6–2.7 \\(M_\\odot\\)) are typically recycled into long-lived millisecond pulsars, whereas high-mass BNSs (Mtot above about 3 \\(M_\\odot\\), such as GW190425) retain strong magnetic fields and have short radio lifetimes, so they are largely absent from radio surveys. The model is compared with the observed P–Pdot, Porb–e, and total-mass distributions of the roughly 19 measured Galactic BNSs and is claimed to match them well, and it predicts that 19–22% of BNSs merging at redshift zero have Mtot > 3 \\(M_\\odot\\). The paper also discusses sensitivity to the wind velocity law, the subsonic accretion model, and the neutron star equation of state.","tokens_in":1725,"tokens_out":1749,"duration_ms":89471,"significance":"The paper addresses a genuine and timely puzzle: the mass discrepancy between radio-detected Galactic BNSs and the gravitational-wave events GW170817 and GW190425. Its strengths include a detailed treatment of multiple accretion and spin-down phases, the use of observationally motivated star formation and metallicity enrichment histories for the Milky Way, and explicit sensitivity tests for the wind prescription, the subsonic accretion mechanism, and the EOS. The prediction that about 19–22% of z=0 BNS mergers have Mtot > 3 \\(M_\\odot\\) is falsifiable with future gravitational-wave observations, and the proposed mechanism (progenitor-mass-dependent initial fields plus accretion-induced field decay) is physically plausible. However, the central selection-effect conclusion currently rests on a small number of unvaried assumptions, and the claimed agreement with observation is not quantified statistically, so the result should be viewed as promising rather than firmly established.","major_comments":[{"comment":"The central claim that high-Mtot BNSs are radio-quiet depends on two assumptions that are not varied in the robustness tests of Section 4: (i) neutron stars from ZAMS stars with M >= 20 \\(M_\\odot\\) are assigned a uniform initial magnetic field of 1e14–1e15 G rather than the log-normal distribution (mean 1e13 G) used for lower-mass progenitors, and (ii) the magnetic-field decay timescale during accretion is td = 1e6 yr. A concrete test shows the issue: if the same log-normal initial field is used for all neutron stars, Eq. (30) gives B ~ 2e11 G after 4 Myr of accretion, Eq. (29) gives an equilibrium spin period near 0.1 s, and the death-line crossing time from Eq. (11) is on the order of 100 Myr, so heavy BNSs would be radio-visible for a substantial fraction of a Hubble time and the observed absence of Mtot > 3 \\(M_\\odot\\) radio BNSs would not be reproduced. Since the radio/GW mass split is the paper's main result, the authors should repeat the population synthesis with these alternatives, or provide a quantitative justification for why the bimodal initial-field and fast-decay assumptions are necessary.","section":"Section 2.2.3, Eq. (30), and Section 4"},{"comment":"The manuscript repeatedly states that the model 'can well match' the observed P–Pdot, Porb–e, and total-mass distributions, but this is based on visual inspection with only ~19 observed systems and no statistical test (e.g., a two-sample Kolmogorov–Smirnov or Anderson–Darling test) is reported. Because the model produces thousands of synthetic pulsar-bearing BNSs, such tests on the relevant marginals are feasible and would give a quantitative measure of agreement. Without them, the 'well match' language is stronger than the evidence supports, especially given the small observed sample and the acknowledged neglect of radio selection effects.","section":"Section 3.2, Figs. 6 and 7"},{"comment":"The paper explicitly sets aside beaming, radio luminosity, and survey selection, treating every neutron star with P < P_death (Eq. (11)) as equally observable. For the central argument that the radio-selected sample is biased against high-Mtot BNSs, the relevant quantity is not radio lifetime alone but lifetime multiplied by beaming fraction and by the luminosity/selection probability. If beaming or luminosity correlate with spin period or magnetic field—as is observed for recycled versus young pulsars—the relative detectability of low- and high-mass systems could change, potentially altering the conclusion. The authors should include a simple beaming/luminosity prescription as a sensitivity test, or demonstrate explicitly that the lifetime effect dominates over plausible variations of these factors.","section":"Section 3.2 and Section 4"}],"minor_comments":[{"comment":"The first paragraph of the conclusions refers to the 'P – e diagram,' but this should be the 'P_orb – e diagram.'","section":"Section 5"},{"comment":"The phrase 'in the the cases C and D' contains a duplicated article and should read 'in the cases C and D.'","section":"Section 4"},{"comment":"The author name in 'Os lowski et al. (2011)' has an erroneous space and should be 'Osłowski et al. (2011).'","section":"Section 2.2.3"},{"comment":"The abbreviation 'HXMB' appears in several places (e.g., 'the HXMB stage') and should be 'HMXB.'","section":"Section 3.1"},{"comment":"The definitions of R_mag,1 and R_mag,2 are written with ambiguous parentheses; for example, 'R_mag,1 = (µ^2/2 ẌM √2GMNS)^(2/7)' should be typeset as [µ^2 / (2 ẌM √(2 G M_NS))]^(2/7).","section":"Section 2.2"},{"comment":"In the sentence citing the star formation and metallicity enrichment history, the citation to Licquia & Newman (2015) is misplaced because that work provides the total stellar mass of the Milky Way, not the SFH; please adjust the citation or the attributions.","section":"Section 4"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and makes a useful contribution to a debated problem. The main barrier to acceptance is that the central selection-effect conclusion is not yet robust to the two key magnetic-field assumptions described in Section 2.2.3, and the paper's robustness section does not vary them. I would support publication after the authors add the requested sensitivity runs and a quantitative statistical comparison with the observed BNS sample."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nThe paper is worth a careful read: it builds a detailed spin-evolution model for BNSs and argues that the narrow total-mass peak of radio-detected Galactic BNSs and the heavier GW190425 select different populations because high-mass BNSs are radio-quiet. That is a real and testable claim, and it goes beyond the earlier COMPAS work by Chattopadhyay et al. (2020), which the authors cite and compare against.\n\nWhat is genuinely good: the machinery is transparent, and the paper documents its recipes thoroughly, including accretion phases, field decay, kicks, and the star formation history. The forward model uses independent inputs for the most part—isolated pulsar field distributions, HMXB torque models, observed SFH—so the P–Pdot match is not a fit to the very data it explains. The paper also admits what it leaves out: beaming, luminosity, and survey selection, and it flags the small-number statistics of the observed sample. Those are honest caveats, not hidden ones.\n\nThe soft spots are proportional. The headline comparisons in the P–Pdot, Porb–e, and total-mass diagrams are visual; there are no statistical tests, and with only ~19 observed systems the discrimination power is modest. More importantly, the mass dependence of the radio lifetime rests on two coupled assumptions in Section 2.2.3: neutron stars from progenitors above 20 Msun are assigned uniform 1e14–1e15 G initial fields, and the accretion-driven field decay timescale is 1e6 yr. The sensitivity tests in Section 4 vary wind prescription, accretion model, and EOS, but not these two parameters. The stress-test arithmetic is right: if the high-mass progenitors draw from the same log-normal as the low-mass ones, the post-HMXB fields drop to ~2e11 G and heavy systems stay radio-bright for more than 100 Myr, which would erase the predicted radio/GW mass split. That is a load-bearing joint, and it has not been stress-tested.\n\nI would still send this to a referee. The model is coherent, the prediction about the merging fraction (19–22% above 3 Msun) is falsifiable with more GW events, and the paper is honest about its own uncertainty budget. The referee should push on the field assumptions and ask for either a variation of B_ini and t_d across the full population or a physical argument for why the magnetar-field assignment is robust. No code or data are provided, which makes the detailed numbers hard to reproduce; that is a separate but real cost.\n\nFor your reading group: maybe, if people want to see how far population synthesis can go with spin evolution. I would not cite it as the last word, but I would cite it as the clearest statement of the mass-dependent radio-lifetime selection effect.\n\nBest,\n[You]","headline":"Careful spin-evolution population synthesis with a genuinely new selection-effect claim, but the mass-dependent radio-lifetime result hangs on unvaried magnetar-field assumptions that a referee should push hard on.","tokens_in":29598,"tokens_out":2211,"would_cite":true,"duration_ms":20681,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper argues that the radio-silent majority of heavy binary neutron stars explains why the Milky Way's radio-selected systems cluster near 2.6–2.7 $M_\\odot$ while GW190425 weighs about 3.4 $M_\\odot$.","keywords":["binary neutron stars","neutron star spin evolution","pulsar recycling","radio selection effects","gravitational-wave sources","population synthesis","neutron star magnetic fields","magnetars"],"falsifier":"Find one Galactic binary neutron star with a recycled, radio-emitting pulsar and a total mass at or above 3 $M_\\odot$, well outside the observed 2.6–2.7 $M_\\odot$ peak. Under the model's magnetic-field assumptions such a system should be extremely rare; a single secure detection, or a radio survey that finds a population of them, would falsify the selection-effect explanation. A complementary test is to measure the initial spin-down-inferred fields of young pulsars whose progenitors' masses are constrained to exceed 20 $M_\\odot$; if those fields fall below $10^{14}$ G, the predicted radio silence of heavy binaries disappears.","tokens_in":28438,"feed_emoji":"🔭","tokens_out":7991,"duration_ms":59725,"temperature":0.7,"pith_summary":"This paper argues that the mass difference between radio-detected binary neutron stars in the Milky Way (total masses clustered near 2.6–2.7 $M_\\odot$) and the gravitational-wave event GW190425 (about 3.4 $M_\\odot$) is a selection effect rather than evidence for different formation channels. Using an upgraded binary population synthesis code that tracks neutron star spin through every accretion phase, the authors show that low-mass systems recycle the first-born neutron star into a long-lived millisecond pulsar, so they remain radio-bright for up to a Hubble time. Most high-mass systems, by contrast, keep strong magnetic fields and fade from the radio band within roughly a million years, making them nearly invisible to radio surveys. The model reproduces the observed $P{-}\\dot{P}$, $P_{\\rm orb}{-}e$, and total-mass distributions of Galactic binary neutron stars, and predicts that 19–22 percent of merging binaries at redshift zero weigh more than 3 $M_\\odot$, consistent with the gravitational-wave detections.","feed_headline":"Heavy neutron-star pairs hide from radio telescopes","feed_subtitle":"Spin evolution keeps low-mass binaries shining for eons while heavy ones fade in a million years.","key_machinery":"The load-bearing object is the neutron star spin evolution model grafted onto the binary star evolution code, which advances the first-born neutron star through an ejector phase, rapid rotator, propeller, super- and sub-Keplerian magnetic inhibition, subsonic settling accretion, Bondi-Hoyle-Littleton accretion, and finally Roche-lobe overflow spin-up. Each phase has an analytic torque or spin-down rate, with a death line separating radio-loud from radio-quiet pulsars, and an exponential magnetic field decay with a shorter timescale during accretion ($10^6$ yr) than otherwise ($10^9$ yr). The crucial input is the initial field: neutron stars from progenitors below 20 $M_\\odot$ draw log-normal fields around $10^{13}$ G, while those from more massive stars are assigned uniform fields of $10^{14}$–$10^{15}$ G, producing the short radio lifetimes of heavy binaries.","core_discovery":"The paper's central claim is that the paucity of massive binary neutron stars in the radio-selected Galactic sample is a consequence of neutron star spin and magnetic field evolution during the high-mass X-ray binary phase. First-born neutron stars in low total-mass binaries experience prolonged Roche-lobe overflow, which recycles them into millisecond pulsars whose radio lifetime can exceed the Hubble time; the observed narrow mass peak near 2.6–2.7 $M_\\odot$ therefore reflects what survives as radio-bright. In massive binaries the donor star is heavier than about 20 $M_\\odot$, so the newborn neutron star is assumed to be a magnetar with an initial field of $10^{14}$–$10^{15}$ G, which shortens the radio lifetime to about $10^5$ yr and leaves most heavy systems radio-quiet. With the Milky Way's star formation and metallicity history folded in, the model matches the observed $P{-}\\dot{P}$ and $P_{\\rm orb}{-}e$ diagrams and the total-mass histogram, and yields a merger population at $z\\sim0$ with 19–22 percent of systems above 3 $M_\\odot$, which the authors take to be compatible with GW170817 and GW190425.","pith_inferences":["If the model is right, third-generation gravitational-wave detectors should see a high-mass tail in the merger mass distribution that radio-selected Galactic samples will never fully reproduce, making the two observables complementary probes of the same population.","The same mass-dependent recycling logic may apply to other compact binaries, so any radio-versus-gravitational-wave mass discrepancy in future data could be a sign of selection rather than of new formation physics.","A targeted test is to measure magnetic fields of young pulsars whose natal progenitors can be shown, through kinematics or association, to have been more massive than 20 $M_\\odot$; a field distribution that does not extend to $10^{14}$–$10^{15}$ G would weaken the heavy-binary radio invisibility."],"forward_implications":["Radio surveys of the Milky Way should keep seeing recycled-pulsar binaries only below about 3 $M_\\odot$; heavy systems like GW190425 should stay radio-silent even if they are intrinsically common.","The intrinsic merger population is heavier than the radio-selected Galactic sample suggests: 19–22 percent of binary neutron star mergers at $z\\sim0$ exceed 3 $M_\\odot$.","A sample of roughly 60–100 additional Galactic binary neutron stars should sharpen the bimodal mass distribution and directly test whether the selection effect explains the missing heavy systems.","Radio-detectable Galactic binary neutron stars should be dominated by single recycled pulsars paired with radio-quiet neutron stars, with pulsar-pulsar binaries making up only about 3.7 percent of the detectable population."],"supporting_citations":[{"why":"supplies the binary star evolution code that the spin-evolution model is grafted onto.","marker":"Hurley et al. (2000, 2002)"},{"why":"provides the wind-fed accretion phases and torque formulas used for the first-born neutron star.","marker":"Karino (2020)"},{"why":"defines the death line and the accretion-recycling framework for pulsar spin evolution.","marker":"Bhattacharya & van den Heuvel (1991)"},{"why":"constrains the initial magnetic field distribution adopted for neutron stars from lower-mass progenitors.","marker":"Faucher-Giguère & Kaspi (2006)"},{"why":"adds constraints on the initial field distribution for the same low-mass progenitor population.","marker":"Igoshev & Popov (2013)"},{"why":"supplies the magnetar sample motivating the strong initial fields assigned to neutron stars from massive progenitors.","marker":"Olausen & Kaspi (2014)"},{"why":"provides the remnant mass function used to set neutron star masses from core-collapse and ultra-stripped supernovae.","marker":"Fryer et al. (2012)"},{"why":"gives the observed bimodal Galactic BNS mass distribution that the model must reproduce.","marker":"Huang et al. (2018)"},{"why":"defines case BB mass transfer and ultra-stripped supernovae, the channel that recycles low-mass first-born neutron stars.","marker":"Tauris et al. (2013, 2015)"},{"why":"provides an independent population synthesis whose $P{-}\\dot{P}$ distribution is compared with the present model.","marker":"Chattopadhyay et al. (2020)"}],"fun_headline_variants":["Why heavy neutron-star pairs are radio silent","Magnetars doom heavy neutron-star pairs to radio silence","Spin evolution hides massive neutron-star pairs from radio","Light neutron-star pairs outshine heavy ones on radio"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The explanation relies on neutron stars born from stars heavier than 20 solar masses being born with magnetic fields near magnetar strength ($10^{14}$–$10^{15}$ G) that do not decay away before the binary becomes radio-quiet; if those fields are weaker or decay faster, heavy binaries would remain radio-bright and the mass discrepancy would need another explanation.","fun_headline_variants_meta":{"raw":{"variants":["Why heavy neutron-star pairs are radio silent","Magnetars doom heavy neutron-star pairs to radio silence","Spin evolution hides massive neutron-star pairs from radio","Light neutron-star pairs outshine heavy ones on radio"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000803,"raw_usage":{"total_tokens":3621,"prompt_tokens":1127,"completion_tokens":2494,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":743,"completion_tokens_details":{"reasoning_tokens":2433}},"tokens_in":743,"tokens_out":2494,"duration_ms":16295,"temperature":1.0,"reasoning_tokens":2433,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T11:24:47.099722+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Find one Galactic binary neutron star with a recycled, radio-emitting pulsar and a total mass at or above 3 $M_\\odot$, well outside the observed 2.6–2.7 $M_\\odot$ peak. Under the model's magnetic-field assumptions such a system should be extremely rare; a single secure detection, or a radio survey that finds a population of them, would falsify the selection-effect explanation. A complementary test is to measure the initial spin-down-inferred fields of young pulsars whose progenitors' masses are constrained to exceed 20 $M_\\odot$; if those fields fall below $10^{14}$ G, the predicted radio silence of heavy binaries disappears.","supporting_citations":[{"cited_title":"2020, , 72, 95, 10.1093/pasj/psaa087","cited_arxiv_id":null,"evidence_quote":"provides the wind-fed accretion phases and torque formulas used for the first-born neutron star."}],"review_version":1}