{"id":"fd401d7c-3cd2-4986-aa3d-5e0d60921fa7","arxiv_id":"2412.15521","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Isolated binary evolution at solar metallicity produces essentially no millisecond pulsar-black hole binaries because the neutron star never accretes after formation.","lead":"This paper simulates how often neutron star-black hole binaries form in the Milky Way with the neutron star appearing first. It finds these systems never experience mass transfer after the neutron star forms, so their pulsars cannot spin up, meaning millisecond pulsar-black hole binaries likely cannot form from isolated pairs.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The categorical 'cannot form' conclusion rests on POSYDON's rotation-limited non-degenerate accretion; if that prescription is even partially wrong, the COMPAS channel reopens, so the null result is not yet supported as a general barrier.","rationale":"The reader's weakest-assumption analysis already identifies the rotation-limited accretion prescription as the load-bearing point, and my stress-test converges on the same concern. The paper is a careful, transparent population synthesis study; the disagreement between POSYDON and COMPAS is real and is explicitly discussed in Section 4.1. However, the abstract overstates the result by saying MSP-BH binaries 'cannot form' when the only reason the isolated channel is closed is a single uncertain microphysical prescription for accretion onto non-degenerate stars. This does not invalidate the paper's modeling or its useful catalog of NSBH formation channels, but it means the central claim should remain conditional. The reader's CONDITIONAL verdict already captures this correctly, so no verdict adjustment is needed. I do not see a more serious internal inconsistency or a reason to reject the paper; the main limitation is the untested sensitivity of the headline conclusion to the accretion treatment.","tokens_in":22007,"tokens_out":4200,"duration_ms":45365,"concrete_test":"Re-run the S16/fiducial population through POSYDON with a parameterized non-degenerate accretion treatment that decouples mass retention from the accretor's critical rotation (e.g., constant retention fractions beta = 0.5 and 1.0), keeping all other assumptions fixed, and check whether any NSBH binary experiences stable RLO or CE after the first SN. A complementary test is to take the recycled-NSBH progenitor binaries from Chattopadhyay et al. (2021) that COMPAS evolves into MSP-BH systems and evolve those exact initial conditions through POSYDON's HMS-HMS MESA grids with the rotationally enhanced wind term disabled. If post-first-SN mass transfer appears in either test, the null result is driven by the accretion prescription rather than by a general formation barrier.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim is that no NSBH binary in its populations undergoes mass transfer after the first supernova, hence MSP-BH binaries cannot form in isolation. The paper itself, in Section 4.1, shows that this null result is a direct consequence of POSYDON's treatment of mass accretion onto non-degenerate companions: accretion is limited by the accretor's critical rotation rate, and rotationally enhanced winds eject any additional transferred mass (Section 2.2). This restricts mass-ratio reversal during HMS-HMS evolution to binaries with near-unity ZAMS mass ratios, so the secondary evolves on a similar timescale and either triggers a double CE before the first SN (Channel II) or remains detached afterward (Channel I). In contrast, COMPAS assumes nearly conservative accretion, which permits mass-ratio reversal for more unequal binaries and enables a CE phase with accretion onto the NS after the first SN. The paper offers no independent calibration or test of the rotation-limited accretion assumption; it even notes that star-disk interactions could alter it (Section 2.2). Section 4.2 acknowledges uncertainties in supernova prescriptions, double CE onset, and grid resolution, but not the accretion assumption that underpins the no-MT result. Therefore the abstract's categorical statement is not robust to a plausible, code-dependent modeling choice, and the conclusion should be presented as conditional on this accretion physics rather than as a universal formation barrier.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper uses the POSYDON v2 population-synthesis code to evolve 10^7 solar-metallicity binaries and isolate the subpopulation of NS\\,–\\,BH binaries in which the neutron star forms first (NSBH). It reports Galactic NSBH birth rates below about 1.6\\,Myr$^{-1}$, two to three orders of magnitude below BH-first systems, and identifies two formation channels: fully detached evolution after the hydrogen main sequence (Channel I) and a double common-envelope phase before the first supernova (Channel II). The central result is that no NSBH binary in any of the tested populations undergoes mass transfer, stable or unstable, after the first supernova, so the neutron star cannot be recycled into a millisecond pulsar. The authors conclude that isolated binary evolution cannot form MSP\\,–\\,BH binaries and that dynamical environments may be the only remaining channel, placing this result in explicit contrast with the COMPAS-based conclusions of Chattopadhyay et al. (2021).","tokens_in":22269,"tokens_out":4136,"duration_ms":41802,"significance":"If the central null result holds, the paper would settle an important disputed question in binary population synthesis by closing the isolated formation channel for MSP\\,–\\,BH binaries at solar metallicity, with direct consequences for pulsar-survey targeting. The study's strengths are that the result emerges from detailed MESA-based binary grids rather than from fitted outcomes, and that the authors test multiple supernova remnant prescriptions, kick settings, common-envelope efficiencies, and an alternative interpolation scheme, while transparently listing many model limitations. The main caveat, identified in the paper itself in Section 4.1, is that the absence of post-supernova mass transfer is a direct consequence of POSYDON's rotation-limited accretion prescription for non-degenerate companions. Because that prescription is not independently calibrated or varied here, the categorical statement that MSP\\,–\\,BH binaries 'cannot form' is stronger than the modeling actually supports.","major_comments":[{"comment":"The central null result is load-bearing on a single code-specific assumption: POSYDON's treatment of mass accretion onto non-degenerate stars, which is limited by the accretor's critical rotation rate and by rotationally enhanced winds. The paper itself states in §4.1 that COMPAS's more conservative accretion allows mass-ratio reversal for more unequal ZAMS binaries, leading to a common-envelope phase and accretion onto the neutron star after the first supernova, whereas POSYDON restricts mass-ratio reversal to near-unity mass ratios. Section 2.2 also notes that star\\,–\\,disk interactions could change the accretion and spin-up assumptions. No test, calibration, or quantitative comparison of this assumption is provided. Because the abstract's conclusion \"MSP\\,–\\,BH binaries cannot form\" rests on this prescription, the claim is not yet robust. I would ask the authors to add a model variation that relaxes the rotation-limited accretion assumption (for example, allowing a fixed fraction of the transferred mass to be accreted regardless of spin, or adopting a COMPAS-like conservative accretion limit) and to report whether any post-supernova mass-transfer channel appears; alternatively, the conclusion should be explicitly conditional on the accretion treatment rather than stated as a universal formation barrier.","section":"§4.1, §2.2"},{"comment":"The model-uncertainties section does not include the rotation-limited accretion assumption among the factors that could change the main result. Section 4.2 discusses uncertainties in the core-collapse prescription, double-CE onset, stellar winds, interpolation, and grid resolution, but the discussion of the comparison with COMPAS in §4.1 makes clear that the absence of post-supernova mass transfer is controlled by the accretion physics. Since the abstract draws a categorical conclusion from this null result, the limitation section should explicitly state that the no-recycling result is conditional on POSYDON's treatment of non-degenerate accretion, and should discuss what would happen if that treatment were modified. This is not a cosmetic omission; it is the difference between a robust population-synthesis result and a code-dependent artifact.","section":"§4.2"},{"comment":"The conclusion that \"MSP\\,–\\,BH binaries cannot be formed\" overstates the scope of the simulations. The models cover only solar metallicity, only isolated binaries, and only the specific set of physical prescriptions adopted in POSYDON. The abstract itself contains the qualifier \"in isolation,\" but the final sentence, \"Thus, dynamical environments and processes may provide the only formation channels,\" presents the result as an absolute barrier. A more accurate statement would be that no MSP\\,–\\,BH binaries form through the isolated-binary channels modeled here at solar metallicity, given the tested prescriptions. The distinction matters for the paper's survey-strategy implications, and the authors should either soften the claim or provide evidence that the tested parameter space is representative of the physical possibilities.","section":"Abstract and §5"}],"minor_comments":[{"comment":"In the caption, \"potted\" should be \"plotted\": \"Binaries with initial mass ratios ≤ 0.95 are potted on the q = 0.9 grid slice.\"","section":"Figure 4 caption"},{"comment":"The caption begins \"T able 1.\" with a stray space; it should read \"Table 1.\"","section":"Table 1 caption"},{"comment":"The acronym \"F AST\" in the list of surveys should be \"FAST\".","section":"Introduction"},{"comment":"The phrase \"supported by the the Swiss National Science Foundation\" contains a duplicated \"the\" and appears twice.","section":"Acknowledgments"},{"comment":"The sentence describing the rotationally enhanced wind model says the wind keeps the stellar rotation rate \"always below its critical threshold,\" but the same paragraph discusses critically rotating stars; consider writing \"at or below\" for consistency with the earlier description.","section":"Section 2.2"},{"comment":"The legend labels \"Stable reverse mass-transfer phase\" and \"Unstable reverse mass-transfer phase\" are rendered in different colors in the figure, but the caption text does not define all color/symbol combinations unambiguously; a compact legend description in the caption would improve readability.","section":"Figure 4"}],"recommendation":"major_revision","confidential_remarks":"The paper is scientifically honest and clearly written, and the disagreement with the COMPAS-based study is handled fairly. My main concern is that the headline result, the impossibility of forming MSP\\,–\\,BH binaries in isolation, is explicitly tied to POSYDON's rotation-limited accretion prescription in Section 4.1, yet the abstract and conclusions present it as a universal barrier. This is fixable either by adding a test that varies the accretion efficiency or by conditioning the conclusion on the accretion treatment. I do not see any issue with the paper's fit to the journal or with novelty; the concern is about the strength of the claim relative to the evidence."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a careful paper that deserves referee time, but the headline claim needs a qualifier. The novel piece is not the qualitative conclusion—Kruckow et al. 2018 already got no recycled NSBHs from inefficient mass transfer—but the quantitative act: POSYDON v2 rates, the 100% mass ratio reversal, the 64% double CE fraction in the S16 model, and a clean mechanistic explanation for why COMPAS (Chattopadhyay et al. 2021) and ComBinE (Kruckow et al. 2018) disagree. The explanation—rotation-limited accretion onto non-degenerate companions prevents mass ratio reversal except for near-unity ZAMS mass ratios, so no post-first-SN mass transfer occurs—is both the paper's strength and its soft spot.\n\nWhat it does well: the model variations are sensible (two remnant prescriptions, two kick settings, three CE efficiencies, two interpolation schemes), and the null result holds across all of them. The authors are open about limitations: double CE onset is uncertain, grid resolution matters, future SN prescriptions could open new channels. The comparison with earlier codes is honest and detailed rather than dismissive. That is how population synthesis disagreements should be handled.\n\nThe soft spot is the load-bearing one: the no-mass-transfer outcome is a direct consequence of the rotation-limited accretion treatment, and the paper does not test that assumption. It even notes in Section 2.2 that star-disk interactions could change the spin-up and accretion behaviour, and in Section 4.1 that COMPAS's nearly conservative MT is what allows the recycled channel. If real accretion onto MS companions is more efficient than POSYDON assumes, the isolated MSP-BH channel reopens. The abstract's \"cannot form\" is therefore stronger than the evidence supports; Section 4.2 acknowledges future CC prescriptions as a caveat but not the accretion physics. Also, birth rates are quoted without uncertainties and the simulated populations are not released, which limits reproducibility for a paper whose main unit is a model population.\n\nOverall: this is a solid, useful study for pulsar survey targeting and for the NSBH population synthesis community. It does not settle the isolated channel question, but it sharpens it considerably and provides a concrete mechanism to test. A serious referee should ask the authors to soften the abstract, add an accretion-physics variation or an explicit calibration argument, and publish the population data. I would take it to journal review as is, with those requests.","headline":"Solid, transparent population synthesis study whose main null result is only as strong as POSYDON's rotation-limited accretion prescription.","tokens_in":22888,"tokens_out":2032,"would_cite":true,"duration_ms":18614,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper argues that isolated binary evolution at solar metallicity cannot form millisecond pulsar–black hole binaries, because no mass transfer onto the neutron star occurs after the first supernova.","keywords":["millisecond pulsars","neutron star-black hole binaries","binary population synthesis","mass transfer","common envelope evolution","pulsar recycling","supernova kicks","globular clusters"],"falsifier":"A confirmed detection of a millisecond pulsar in a binary with a black hole in the Galactic field, with no evidence of a dynamical exchange history, would falsify the claim that isolated evolution cannot form such systems. A more direct computational test would be to rerun the same populations with an accretion treatment that is not capped by the accretor's critical rotation and check whether any post-first-supernova mass transfer onto the neutron star appears.","tokens_in":21779,"feed_emoji":"🔭","tokens_out":6320,"duration_ms":48582,"temperature":0.7,"pith_summary":"This paper models the formation of neutron star–black hole binaries in the Galactic field at solar metallicity, concentrating on systems where the neutron star forms first and could in principle be spun up to millisecond periods. It finds such systems are rare, with birth rates below one per Myr in a Milky Way-like galaxy, and 2–3 orders of magnitude rarer than binaries where the black hole forms first. The paper reports that none of its simulated neutron-star-first binaries experience any mass transfer, stable or unstable, after the first supernova, which removes the only known mechanism for spinning the neutron star up to millisecond rotation. The conclusion, if correct, is that millisecond pulsar–black hole binaries cannot form from isolated binaries, and pulsar surveys should target dynamically active environments instead.","feed_headline":"Isolated binaries cannot make millisecond pulsar-black holes","feed_subtitle":"Population synthesis finds the neutron star never accretes after the first supernova, so only dense star clusters remain as birth sites.","key_machinery":"The load-bearing mechanism is the code's treatment of rotationally-limited accretion onto non-degenerate companions: a star can only accept mass until it reaches critical rotation, after which boosted winds eject further transferred material. This makes mass ratio reversal possible only for binaries with zero-age main-sequence mass ratios near unity, and it prevents the companion from later filling its Roche lobe and transferring mass onto the neutron star after the first supernova. The other central element is the double common envelope treatment, which lets near-equal-mass binaries tighten their orbits before the first supernova but still leaves them detached afterward.","core_discovery":"The paper claims that, within the physics implemented in its population synthesis code, every neutron-star-first black hole binary progenitor undergoes a mass ratio reversal before the first supernova, through mass transfer during hydrogen main-sequence evolution. These progenitors then either evolve fully detached (Channel I) or pass through a double common envelope phase before the first supernova (Channel II), but in both channels no mass transfer occurs after the first supernova. Because accretion is required to recycle a neutron star to millisecond periods, the authors conclude that millisecond pulsar–black hole binaries have no isolated formation channel in these models, leaving dynamical formation as the only remaining route.","pith_inferences":["A testable extension left implicit by the paper: if isolated millisecond pulsar–black hole binaries do form, they would most likely appear as wide, detached, unrecycled pulsar–black hole systems, so a field discovery of a recycled pulsar with a black hole companion would directly contradict this channel.","The gap between this result and an earlier population synthesis study that found recycled neutron stars appears to be driven mainly by the rotation cap on accretion; a systematic comparison of accretion physics, rather than supernova physics, would decide which population's neutron-star-first binaries are closer to reality.","Future supernova remnant prescriptions that produce different remnant masses, or a revised physical treatment of double common envelope onset, could reopen a narrow isolated formation channel even within this code.","The rotation-limited accretion assumption could be tested observationally by comparing simulated mass ratios of double neutron star and X-ray binary populations against measured systems, since those constrain how much mass non-degenerate accretors actually retain."],"forward_implications":["If correct, current and upcoming pulsar surveys should expect no Galactic-field millisecond pulsar–black hole detections from isolated binary evolution.","Dynamically active environments, such as globular clusters and nuclear clusters, become the only plausible birth sites for these systems.","The neutron-star-first black hole birth rate below 1 Myr$^{-1}$ means these systems contribute negligibly to gravitational-wave merger rates compared with black-hole-first systems.","The properties of these binaries, including which formation channel dominates, depend strongly on the supernova remnant prescription, natal kick assumptions, and common-envelope efficiency.","The one detected millisecond pulsar with a compact companion in NGC 1851 is consistent with a dynamical origin, matching the paper's conclusions."],"supporting_citations":[{"why":"Describes the population synthesis code's pre-computed grids, mass-transfer classification, and rotationally-limited accretion scheme that produce the null result.","marker":"Fragos et al. 2023"},{"why":"Introduces the code version used here, including reverse mass transfer and the double common envelope treatment in the binary model grids.","marker":"Andrews et al. 2024"},{"why":"The earlier study that reaches the opposite conclusion, finding recycled neutron stars in its populations; it is the central counterpoint the paper must explain.","marker":"Chattopadhyay et al. 2021"},{"why":"Earlier population synthesis that also finds no common envelope or stable mass transfer after the first supernova, supporting the paper's main qualitative result.","marker":"Kruckow et al. 2018"},{"why":"Provides the delayed remnant-mass prescription used for one set of the paper's models, which affects which binaries form neutron-star-first black holes.","marker":"Fryer et al. 2012"},{"why":"Provides the alternative remnant-mass prescription whose stochastic explosion behavior widens the parameter space for neutron-star-first black hole formation.","marker":"Sukhbold et al. 2016"},{"why":"Reports the discovered millisecond pulsar with a compact companion in NGC 1851, the observational anchor for the paper's dynamical-formation conclusion.","marker":"Barr et al. 2024"}],"fun_headline_variants":["No isolated path to pulsar-black hole binaries","Isolated binaries can't make MSP-BH pairs","Accretion never happens: no isolated MSP-BH","Pulsar spin-up impossible in isolation","Only dense clusters can make MSP-BH"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The null result rests on the code's assumption that a non-degenerate star can only accrete until it reaches critical rotation, so only near-equal-mass binaries reverse their mass ratios and no companion can later dump mass onto the neutron star.","fun_headline_variants_meta":{"raw":{"variants":["No isolated path to pulsar-black hole binaries","Isolated binaries can't make MSP-BH pairs","Accretion never happens: no isolated MSP-BH","Pulsar spin-up impossible in isolation","Only dense clusters can make MSP-BH"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001024,"raw_usage":{"total_tokens":4322,"prompt_tokens":951,"completion_tokens":3371,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":567,"completion_tokens_details":{"reasoning_tokens":3299}},"tokens_in":567,"tokens_out":3371,"duration_ms":18447,"temperature":1.0,"reasoning_tokens":3299,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T11:20:53.478497+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A confirmed detection of a millisecond pulsar in a binary with a black hole in the Galactic field, with no evidence of a dynamical exchange history, would falsify the claim that isolated evolution cannot form such systems. A more direct computational test would be to rerun the same populations with an accretion treatment that is not capped by the accretor's critical rotation and check whether any post-first-supernova mass transfer onto the neutron star appears.","supporting_citations":[],"review_version":1}