{"id":"3a631b37-a9ea-4ecb-aa5d-177230e4e3de","arxiv_id":"2608.09395","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"In 2D GRMHD simulations, magnetised neutron stars hyperaccreting inside massive envelopes can halt accretion above B_surf ~2.3e13 G and launch ~1e46 erg/s precursor jets that still cannot unbind the envelope.","lead":"Simulations of a neutron star being swallowed by a massive companion show that strong magnetic fields can briefly stop the star's feeding and fire weak jets, but cannot blow the star apart. These precursor flares are a possible new X-ray signal from common-envelope mergers.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The factor-of-ten amplification from ordinary pulsar fields to the jet threshold (B_surf ≥ 2.3×10^13 G) is asserted via the alpha-omega dynamo, but the axisymmetric model suppresses that dynamo; the paper's own runs at B_surf = 4.64×10^12 G do not launch a jet, so the central astrophysical…","rationale":"The paper's central quantitative and observational payoff is that a hyperaccreting NS inside a common envelope can launch low-power jets before BH formation, and that this is not restricted to exotic magnetars because ordinary pulsar fields need only a factor-of-ten amplification. I read the simulations as honestly computed: the GRMHD code Gmunu is established, the microphysics (two-moment neutrino transport, 13-isotope network) is state of the art, and the dependence of the outcome on initial B_surf is clearly documented. For B_surf = 4.64e13 G, the simulations do show a magnetically dominated funnel and a ~1e46 erg/s polar outflow; for B_surf = 2.32e13 G, the jets are intermittent and depend on eta. The weakness is in the bridge from these initial-condition experiments to the astrophysical population. The paper's own Sec. IV E concedes that the 2D axisymmetric setup suppresses the alpha dynamo and that whether 3D field amplification can halt accretion before collapse is open. Yet Sec. IV B asserts that a modest amplification factor of about 10, 'easily and rapidly achieved' by the alpha-omega dynamo, suffices. Since the applied B_surf is the control parameter that separates jet-forming from non-jet-forming models, and since the simulations cannot regenerate poloidal flux, the factor-of-ten claim is not a result of this paper—it is an unsimulated assumption. The reader's CONDITIONAL verdict already captures this risk; my stress-test does not overturn the paper but sharpens the condition: the conclusions apply to initially strong fields, and the ordinary-pulsar extension is a hypothesis, not a finding. The proposed 3D run directly tests whether the missing alpha effect can close the gap; until then, the factor-of-ten bridge remains the most load-bearing unverified step.","tokens_in":25279,"tokens_out":9886,"duration_ms":108525,"concrete_test":"Run one 3D GRMHD simulation of the B_surf = 4.64e12 G, eta = 0.5 initial condition with the same Gmunu setup but sufficient resolution to resolve the fastest-growing MRI mode in the inner disk (e.g., Q_MRI ≥ 10), and track the poloidal field strength at the NS surface and the jet power at 500 km over the first ~0.3 s. If the poloidal field does not grow by roughly a factor of 10 to reach ~2.3e13 G before collapse, or if no jet with L_jet ≳ 1e46 erg/s breaks out, then the factor-of-ten amplification claim in Sec. IV B is unsupported and the precursor-flare scenario should be restricted to pre-existing magnetar-strength fields.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing step is not the conditional statement 'if B_surf ≥ 2.3e13 G then a jet can form,' which the simulations do show, but the astrophysical claim in Sec. IV B that this threshold is reachable from a canonical pulsar (B_surf ≲ 1e12 G) by an amplification factor of only ~10. That amplification is attributed to the alpha-omega dynamo and magnetic winding, yet the entire simulation suite is axisymmetric, in which Cowling's theorem forbids the alpha effect; the paper's own Fig. 7 shows poloidal field energy only decaying after a transient and never regenerating. The closest model to a pulsar field, B_surf = 4.64e12 G (with eta = 0.3 and 0.5), never launches a sustained jet. Sec. IV B asserts the factor-of-ten increase is 'easily and rapidly achieved,' while Sec. IV E states it 'remains an open question whether fully 3D magnetic field amplification could generate sufficient magnetic pressure to halt the hypercritical accretion and launch a successful jet before collapse occurs.' These two statements are in direct tension, and no simulation or detailed model in the paper supplies the missing amplification. Therefore the paper's broader conclusion that ordinary pulsars can produce precursor flares is not established by its own evidence, and the central claim should be read strictly as an upper-bound conditional for pre-existing strong fields.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents axisymmetric GRMHD simulations of a neutron star embedded in a rotating 15-solar-mass envelope, using grey two-moment neutrino transport and a 13-isotope nuclear network. The authors vary the envelope rotation parameter (eta = 0.1, 0.3, 0.5), the NS spin, and the initial surface magnetic field (B_surf from 4.64e10 to 4.64e13 G). They find that envelope rotation forms a centrifugal barrier and a thick disk, suppressing accretion and neutrino luminosity; that magnetic winding amplifies toroidal fields and drives magnetic towers; and that for B_surf >= 2.32e13 G accretion onto the NS can be halted, evacuating a polar funnel and launching low-power jets of order 1e46 erg/s. The diagnostic energy is found to be insufficient to unbind the massive envelope. The paper also argues that ordinary pulsar fields (B_surf ~1e12 G) need only a factor ~10 amplification to reach the jet-launching threshold, via the alpha-Omega dynamo or magnetic winding.","tokens_in":25588,"tokens_out":6844,"duration_ms":73287,"significance":"If the conditional strong-field results hold, this is a significant step: the study couples state-of-the-art microphysics to a systematic parameter survey of an underexplored common-envelope accretion regime, provides concrete precursor X-ray signatures, and sharpens the argument that a hyperaccreting NS alone cannot power the most energetic CEJSN events. The authors are also unusually explicit about their limitations, and the core forward-modeling framework is appropriate. However, the quantitative claims are restricted by axisymmetry, under-resolved MRI modes, a fixed metric, and a frozen NS core, and the paper's broad observational relevance currently rests on an amplification claim that is not demonstrated by the simulations.","major_comments":[{"comment":"The broadest astrophysical claim, that ordinary pulsar fields (B_surf ~1e12 G) can reach the jet-launching threshold B_surf ~2.32e13 G by only a factor ~10 amplification via the alpha-Omega dynamo or magnetic winding, is not supported by the simulations. The runs are axisymmetric, which by Cowling's theorem suppresses the alpha effect, and Fig. 7 shows that E_Bpol only decays after a transient, with no regeneration; the closest models, B_surf = 4.64e12 G, never launch a sustained jet. This is in direct tension with Sec. IV E, where the authors state that it remains an open question whether fully 3D amplification could generate sufficient magnetic pressure. Please either remove this claim, explicitly mark it as a speculative requirement supplied by physics outside the simulations, or provide a concrete dynamo/amplification model with quantitative support.","section":"IV B, IV E, Fig. 7"},{"comment":"The metric is held fixed and the MHD/radiation variables for r < 8 km are frozen during the evolution. Consequently, the simulations cannot track growth of the NS rest mass or approach to the maximum-mass limit, so the statements that the magnetic barrier 'delays prompt BH formation' and that the NS will eventually 'exceed its maximum mass limit and collapse into a BH' are extrapolations beyond the model. Please qualify these statements as inferences from the suppression of the accretion rate rather than directly simulated results, or extend the model to permit mass growth and collapse.","section":"II B, IV A, IV C"},{"comment":"The quantitative threshold B_surf ~2.32e13 G and the jet power ~1e46 erg/s depend on the MRI being at least partially captured, yet the authors state that Q_MRI < 8 can trigger refinement but the maximum resolution is capped and the fastest-growing MRI modes are not fully resolved. Because no convergence study is reported, it is unclear whether the threshold and jet power are robust or artifacts of limited resolution and axisymmetry. Please add resolution tests (or a dynamo closure) and, in the absence of such tests, present these quantities with explicit error bars or as order-of-magnitude limits rather than precise thresholds.","section":"III B, IV E"},{"comment":"The claim that the polar jet is 'very likely to drill through the polar overburden' is not directly supported by Fig. 14, where the diagnostic energy at t = 100 ms appears to fall well below E_bind,cone ~1e49 erg. The argument requires continuous energy injection over 1-2 s, which is not simulated; please recast this statement as a speculative extrapolation or provide a time-integrated breakout estimate.","section":"IV D, Fig. 14"}],"minor_comments":[{"comment":"In the abstract, 'and evacuates' should be 'and evacuating', and 'may launche' should be 'may launch'.","section":"Abstract"},{"comment":"The legends of Figs. 5 and 7 list the weakest magnetised model as B_surf = 4.64e11 G, while the text in Secs. II A 1 and III B uses B_surf = 4.64e10 G for B_c = 1e11 G. Please make the labels consistent across the figures and text.","section":"Figs. 5 and 7"},{"comment":"Equation (1) is typographically ambiguous; please write A_phi = B_c r_c^3 sin(theta) / [2 (r^3 + r_c^3)] and state explicitly that B_c is the central field strength.","section":"Eq. (1)"},{"comment":"Section IV D cites 'Abrahams et al., in prep' without a reference; please add a citation or remove the attribution.","section":"IV D"},{"comment":"The phrase 'fully coupled' in the abstract should be clarified, since the metric is fixed and the NS interior is frozen; 'fully coupled' refers to the microphysical modules, not to full spacetime dynamics.","section":"Abstract, II B"}],"recommendation":"major_revision","confidential_remarks":"The core conditional simulation results are credible and worth publishing after revision. The main risk is the unsupported ordinary-pulsar amplification claim in Sec. IV B, which overstates the astrophysical reach of the paper and should be removed or clearly demarcated as speculative. The reader's conditional verdict is appropriate; I do not see grounds for rejection, but the present version needs substantive revision of the discussion and of the quantitative presentation."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a good, heavy numerical paper that shows, within its 2D axisymmetric framework, that a hyperaccreting NS inside a massive envelope cannot power the main CEJSN explosion. The jets it produces are weak (~1e46 erg/s), only when the initial surface field is already above ~2e13 G, and the energy budget is far below what's needed to unbind the envelope. That negative result is credible and worth publishing. The soft spot is not the simulations but what the authors claim them to mean: they say ordinary pulsar fields of ~1e12 G can be amplified by a factor of ~10 to reach the jet threshold 'easily and rapidly' via an alpha-omega dynamo. Their own axisymmetric setup forbids the alpha effect, their poloidal field energy only decays after a transient, and their run at 4.64e12 G never launches a jet. That claim is in direct tension with their own Section IV E, which says whether 3D amplification could launch a jet remains open. This overreach should be fixed.\n\nWhat's genuinely new: this is the first study I know of that couples GRMHD, two-moment neutrino transport, and a 13-isotope network in the embedded-NS problem with rotation and initial magnetic fields. The parameter sweep (eta, B_surf, NS spin) is clean. The centrifugal barrier and magnetic tower diagnostics are clear. The conclusion that the NS phase is a precursor and the BH engine must do the main work is a useful, concrete target for the CEJSN community.\n\nWeaker spots, in order: (1) the pulsar-amplification claim above; (2) the MRI is under-resolved and there is no convergence study—the authors admit this, but it means the quantitative jet powers, especially for the 2.32e13 G borderline case, could shift in 3D; (3) the metric is fixed and the NS core frozen, which is fine for a few hundred ms but less so when claiming a magnetospheric barrier delays BH collapse over longer times; (4) they invoke a pre-ejected envelope of 50–90% mass loss to make polar breakout plausible, but that's external to the simulation.\n\nThe reader's take is right. I'd send this to a serious referee. The paper deserves time despite the overreach because the core result—the NS engine is weak—is probably robust, and the field needs that benchmark. The authors should be asked to either support the factor-of-ten with a concrete model or present it explicitly as speculation.\n\nWho is this for? CE/CEJSN modelers, transient theorists thinking about X-ray precursors to BH-driven explosions.","headline":"Solid, useful negative result on the NS engine in CEJSN; the paper's own claim that ordinary pulsar fields can reach the jet threshold via an alpha-omega dynamo is unsupported and contradicts its own axisymmetric results.","tokens_in":26156,"tokens_out":2654,"would_cite":true,"duration_ms":28927,"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":"Strong magnetic fields at a neutron star's surface can completely halt hypercritical accretion and launch low-power precursor jets, while delaying black hole collapse.","keywords":["common envelope evolution","neutron star accretion","hypercritical accretion","GRMHD","magnetic tower","jet launching","precursor flare","black hole formation"],"falsifier":"Run the same hyperaccretion setup in full 3D with resolution that captures the fastest MRI mode and the same initial surface fields; if no magnetospheric barrier forms at $B_{\\rm surf} \\sim 2.3 \\times 10^{13}$ G and no ${\\sim}10^{46}$ erg/s polar outflow appears before the neutron star collapses, the 2D barrier is an artifact. On the observational side, a campaign that finds no X-ray flares coincident with the sudden disappearance of red supergiants would rule out the precursor-flare channel.","tokens_in":25059,"feed_emoji":"🌠","tokens_out":5844,"duration_ms":61240,"temperature":0.7,"pith_summary":"This paper asks what happens when a neutron star is engulfed by a massive companion's envelope and accretes at hypercritical rates. It uses axisymmetric general-relativistic magnetohydrodynamic simulations with neutrino transport and nuclear burning to show that envelope rotation alone forms a centrifugal barrier and a thick accretion disk, suppressing accretion and lowering neutrino luminosity. Adding magnetic fields, it finds that for initial surface fields above about $2.3 \\times 10^{13}$ G, magnetic pressure can completely halt accretion at the neutron star surface, evacuate a low-density polar funnel, and launch low-power precursor jets with powers up to about $10^{46}$ erg/s. These jets delay the neutron star's collapse to a black hole and could appear as X-ray flares, but they lack the energy to unbind the massive envelope, so the subsequent black-hole-driven explosion remains the dominant outcome.","feed_headline":"Strong fields stall neutron star accretion and launch weak jets","feed_subtitle":"Simulations show surface fields above ~2e13 G delay black hole formation and power ~1e46 erg/s precursor flares.","key_machinery":"The central object is the magnetically dominated polar funnel and the magnetic tower that creates it. Differential rotation of the hypercritical accretion flow amplifies the toroidal magnetic field through the $\\Omega$-effect, giving linear field growth with time and toroidal energy growing as $E_{B,\\rm tor} \\propto t^2$; when the initial surface field exceeds ${\\sim}2.3 \\times 10^{13}$ G, the accumulated magnetic pressure halts equatorial accretion at the neutron star surface and evacuates a funnel where the magnetisation satisfies $\\sigma_{\\rm mag} > 1$. The jet power is then measured by integrating the energy flux only within that highly magnetised funnel.","core_discovery":"The central discovery is a magnetic barrier mechanism: a sufficiently magnetised neutron star with polar surface field $B_{\\rm surf} \\gtrsim 2.3 \\times 10^{13}$ G, embedded in a rotating hypercritical accretion flow, can have its surface accretion completely halted by magnetic pressure. Differential rotation winds the poloidal field into a toroidal magnetic tower whose pressure evacuates a low-density polar funnel, enabling a mildly relativistic, collimated outflow with jet power of order $10^{46}$ erg/s. Because that power falls far short of the envelope's binding energy (${\\sim}10^{51}$ erg total, ${\\sim}10^{49}$ erg within the polar cone), the jet acts as a precursor that clears a funnel and delays black hole formation without unbinding the star; the eventual black-hole-driven accretion still controls the final explosion. The paper presents this as an extension of earlier non-rotating, unmagnetised hyperaccretion simulations to the rotating, magnetised regime.","pith_inferences":["Since the paper itself leaves open whether full 3D $\\alpha$-$\\Omega$ dynamo action amplifies fields beyond the 2D $\\Omega$-effect alone, a natural extension is that the $B_{\\rm surf} \\sim 2.3 \\times 10^{13}$ G threshold could be reached from more ordinary pulsar fields, making precursor jets more frequent than the 2D runs suggest.","The choked-jet channel implies a hot, high-entropy cocoon is deposited inside the star before collapse; a subsequent black hole jet would inherit this cavity, potentially producing a delayed, bright transient rather than an immediate one.","A direct observational test is to search for red supergiant 'disappearances' accompanied by a ${\\sim}10^{46}$ erg/s X-ray flare lasting tens to hundreds of milliseconds; such events should be rarer than ordinary supernovae.","The very low $^{56}$Ni yield ($<10^{-9}\\,M_\\odot$) and pristine $Y_e \\sim 0.5$ composition of the unbound ejecta imply that any optical counterpart of the precursor phase would be dim and red, distinguishable from a normal supernova."],"forward_implications":["In models with $B_{\\rm surf} \\gtrsim 2.3 \\times 10^{13}$ G, the accretion rate at 15 km can drop to zero, temporarily decoupling the neutron star from the infall and delaying prompt black hole formation.","The strongest magnetised models sustain jet powers of order $10^{46}$ erg/s at 500 km, sufficient for low-luminosity gamma-ray bursts or X-ray precursors but not for unbinding a red supergiant envelope.","Faster envelope rotation ($\\eta=0.5$) widens the low-density polar funnel and lets the jet break out, while slower rotation ($\\eta=0.3$) leaves the jet choked or stuttering.","Weakly magnetised or slowly rotating engines follow a quiet-collapse track in which the neutron star is overwhelmed by accretion and becomes a black hole, appearing as a failed supernova with little electromagnetic signal.","The magnetically cleared funnel left behind by the precursor jet provides a low-density channel that a later black-hole-driven engine could exploit to power a classical long gamma-ray burst."],"supporting_citations":[{"why":"Establishes the non-rotating, unmagnetised hyperaccretion baseline and the numerical setup that this paper extends with rotation and magnetic fields.","marker":"[19]"},{"why":"Provides the closest prior 3D GRMHD treatment of magnetised accretion onto an unmagnetised neutron star, with simplified neutrino cooling, which this work builds beyond.","marker":"[26]"},{"why":"Introduces the common-envelope jets supernova scenario that the paper tests by asking whether the neutron star engine can launch feedback-driving jets.","marker":"[20]"},{"why":"Supplies the early framework for rapid infall onto neutron stars that motivates the hypercritical accretion and collapse question.","marker":"[18]"},{"why":"Provides the shock-breakout model used to predict the X-ray luminosity of the choked jets and compare with detector sensitivities.","marker":"[62]"},{"why":"Supplies the high-density nuclear equation of state used for the neutron star and hot accretion flow.","marker":"[33]"},{"why":"Provides the code used to generate the rotating neutron star initial models with the adopted spin profiles.","marker":"[27]"}],"fun_headline_variants":["Magnetic barrier halts neutron star infall, powers precursor jets","Magnetized neutron stars stall accretion, launch low-power jets","Strong fields cap neutron star accretion, produce X-ray flares","High-field neutron stars delay black hole death with weak jets","Magnetic pressure clears polar funnel, neutron star jets fall short"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The conclusion rests on the assumption that a 2D axisymmetric, MRI-limited simulation captures the dominant magnetic feedback; if full 3D dynamo action or unresolved MRI turbulence changes the accretion state, the predicted magnetic barrier and precursor jets could fail or look different.","fun_headline_variants_meta":{"raw":{"variants":["Magnetic barrier halts neutron star infall, powers precursor jets","Magnetized neutron stars stall accretion, launch low-power jets","Strong fields cap neutron star accretion, produce X-ray flares","High-field neutron stars delay black hole death with weak jets","Magnetic pressure clears polar funnel, neutron star jets fall short"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000243,"raw_usage":{"total_tokens":1624,"prompt_tokens":1133,"completion_tokens":491,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":749,"completion_tokens_details":{"reasoning_tokens":406}},"tokens_in":749,"tokens_out":491,"duration_ms":5361,"temperature":1.0,"reasoning_tokens":406,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T18:04:16.082645+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same hyperaccretion setup in full 3D with resolution that captures the fastest MRI mode and the same initial surface fields; if no magnetospheric barrier forms at $B_{\\rm surf} \\sim 2.3 \\times 10^{13}$ G and no ${\\sim}10^{46}$ erg/s polar outflow appears before the neutron star collapses, the 2D barrier is an artifact. On the observational side, a campaign that finds no X-ray flares coincident with the sudden disappearance of red supergiants would rule out the precursor-flare channel.","supporting_citations":[{"cited_title":"Effect of Neutron Star Jets on Common Envelope Evolution","cited_arxiv_id":"2607.10267","evidence_quote":"Provides the closest prior 3D GRMHD treatment of magnetised accretion onto an unmagnetised neutron star, with simplified neutrino cooling, which this work builds beyond."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the common-envelope jets supernova scenario that the paper tests by asking whether the neutron star engine can launch feedback-driving jets."},{"cited_title":"Simulating the convection in red super-giant stars: wobbling jets in common envelope evolution","cited_arxiv_id":"2607.13023","evidence_quote":"Provides the code used to generate the rotating neutron star initial models with the adopted spin profiles."}],"review_version":1}