{"id":"27cfec5d-d881-4724-8353-2fcd60625ce6","arxiv_id":"2505.18733","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Two-fluid magnetothermal simulations find that ambipolar diffusion cannot alone produce the persistent X-ray luminosity of magnetars.","lead":"This paper models how magnetic fields decay in neutron star cores and tests whether the released heat can explain magnetar X-ray glow. The simulations show that ambipolar diffusion alone is not enough, even at 5e15 gauss, re-focusing magnetar energy-source models.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Central conclusion hinges on unverified stability of the Grad-Shafranov equilibrium against non-axisymmetric instabilities; if unstable, ambipolar heating could be substantially stronger than modeled.","rationale":"The reader correctly identified the stability of the Grad-Shafranov equilibrium against non-axisymmetric modes as the weakest assumption in the argument. This is the single most load-bearing concern because it directly determines whether the magnetic field reaches a quiescent terminal state or continues to dissipate energy dynamically. The paper explicitly flags this as unresolved in Sec. VI.3, but the abstract and main conclusions are stated without this caveat, so the conclusion is conditional on the stability assumption. I considered other possible concerns: (i) the time-reparametrization procedure (Sec. V.A) appears internally consistent, since the magnetic evolution equations are invariant under the scaling of Eqs. (86)-(89) and the thermal equation (90) properly accounts for the T-dependence of gamma_np; (ii) the neglect of non-equilibrium Urca reactions is justified by the check that xi <~ 1 for the modeled fields (Fig. 7b), so it is not an internal inconsistency; (iii) the vacuum boundary condition at the core-crust interface is acknowledged to overestimate the evolution rate, which makes the heating estimate an upper bound rather than a lower one, reinforcing the paper's conclusion rather than threatening it; and (iv) superfluid/superconducting effects are outside the stated scope of the model. Thus none of these challenge the central claim as directly as the 3D stability question. The proposed concrete test—a linear stability analysis or 3D simulation of the specific equilibria—would settle the matter. Since the paper already receives a CONDITIONAL verdict due to this and related limitations, my assessment does not change the reader's verdict; the manuscript remains acceptable conditionally on the stability concern being addressed or explicitly placed in the conclusion.","tokens_in":31655,"tokens_out":10326,"duration_ms":93243,"concrete_test":"Perform a 3D linear stability analysis (or direct 3D MHD simulation) of the GS equilibria obtained in Sec. IV from models A, B, C0, and C04/C06/C08 after t ~ 0.1 tad. Add non-axisymmetric perturbations with azimuthal orders m=1 and m=2 and evolve for at least 10^4 yr (or compute the growth rates using the local Tayler criterion on the toroidal component). If the equilibria are stable with perturbation growth times longer than tad, the axisymmetric assumption is validated. If unstable modes grow on timescales < 10^3 yr, the heating rate and duration in Sec. V would be underestimated, and the central conclusion would need to be revised.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim—that ambipolar diffusion alone cannot explain magnetar quiescent luminosities—depends on the axisymmetric terminal state being the true endpoint of magnetic evolution. In Sec. IV the authors find that the field relaxes to a Grad-Shafranov equilibrium (Eq. 85), and in Sec. V the thermal balance is computed from the heating rate of this axisymmetric evolution. However, Sec. VI.3 acknowledges that purely poloidal and purely toroidal fields are unstable in 3D, and that barotropic stars may not support stable magnetic equilibria. If the GS equilibrium obtained here is unstable to non-axisymmetric perturbations on timescales comparable to or shorter than the ambipolar timescale tad (Eq. 65), the field would not become quiescent; instead, continued instability-driven dissipation would add magnetic energy release beyond what the axisymmetric runs capture. The comparison in Sec. V.C shows the model's peak surface luminosity at B=5e15 G (Eq. 92, L_s ~ 6.7e34 erg/s) is only a factor of a few below the brightest observed magnetars (e.g., SGR 1806-20, L_s ~ 1e35-3e35 erg/s), so a modest increase in dissipation efficiency or duration from 3D instabilities could close the gap. The paper does not perform any stability analysis of its specific equilibria, leaving this as a true load-bearing uncertainty rather than a conservative approximation.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper studies the long-term magnetothermal evolution of a neutron star core in the weak-coupling regime of ambipolar diffusion. The authors use an axisymmetric two-fluid model (neutrons and a combined charged fluid) with an artificial friction method, solve for the velocity and chemical-potential perturbations, evolve the magnetic field via the induction equation, and treat the thermal evolution through a time reparametrization that accounts for the temperature dependence of the neutron-proton drag. They find relaxation to a Grad–Shafranov equilibrium at constant temperature; with cooling included, only fields B ≳ 5×10^15 G produce a heating-cooling balance, lasting ~10^3 yr, and the corresponding surface luminosities are insufficient to explain the brightest magnetars.","tokens_in":32053,"tokens_out":27398,"duration_ms":215391,"significance":"If the axisymmetric evolution is representative, this is an important quantitative step in modeling core ambipolar diffusion with a two-fluid, GR-informed thermal treatment. The numerical scheme is validated against semi-analytical results in Ref. [44], and energy conservation is explicitly checked (Fig. 5). The scaling relations (91)–(93) provide compact, falsifiable predictions for how the balance temperature, surface luminosity, and equilibration time depend on field strength and length scales. The conclusion—that ambipolar diffusion alone cannot explain magnetar quiescent X-ray luminosities—is a concrete, testable statement that will inform both cooling calculations and future 3D models. The main caveat is the unverified stability of the axisymmetric GS equilibrium, which the authors explicitly acknowledge.","major_comments":[{"comment":"The central conclusion that ambipolar diffusion heating is insufficient to explain magnetar quiescent luminosities rests on the assumption that the axisymmetric Grad–Shafranov equilibrium reached in the simulations is the true terminal state. The paper acknowledges in Sec. VI.3 that purely poloidal and purely toroidal fields are unstable in 3D and that barotropic stars may not support stable equilibria, but it does not test the stability of the specific equilibria found here. This is load-bearing because Eq. (92) predicts L_s ≈ 6.7×10^34 erg/s at B = 5×10^15 G, only a factor of a few below the brightest sources such as SGR 1806-20 (L_s ≈ 1–3×10^35 erg/s); additional dissipation from non-axisymmetric instabilities could close this gap and invalidate the claim that ambipolar heating alone cannot explain the observed luminosities. The authors should either provide at least a linear stability analysis of their GS solutions or explicitly qualify the abstract and Sec. V.C conclusion as applying only to stable axisymmetric equilibria.","section":"Sec. VI.3, Sec. V.C, Eq. (92)"}],"minor_comments":[{"comment":"The numerical prefactor 5.6×10^3 yr in Eq. (65) appears inconsistent with the code units of Table I and with the later Eq. (93). Using the stated parameters (B = 10^15 G, T = 10^9 K, ℓ_B = 2 km, ℓ_c = 10 km) and the formula t_B ≈ 4πγ_np n_n n_c ℓ_B^4/(ℓ_c^2 B^2) gives ≈1.4×10^5 yr, not 5.6×10^3 yr. The scaling in Eq. (93) is consistent with the simulations, so Eq. (65) seems to contain a typographical or unit error that should be corrected.","section":"Sec. II.G, Eq. (65)"},{"comment":"The thermal reparametrization is valid for ξ ≲ 1, and Fig. 7(b) shows that ξ approaches unity at late times for B = 5×10^15 G. Please state explicitly the maximum value of ξ during the heating-cooling balance phase for the cases used in Fig. 9 and confirm that the derived L_s is unaffected by the incipient breakdown of the ξ ≲ 1 condition.","section":"Sec. V.B, Fig. 7(b)"},{"comment":"The initial temperature T∞_init = 10^9 K is above the nominal strong-to-weak coupling transition T_trans ≈ 5×10^8 K. The text justifies this by the negligible field evolution in the strong-coupling regime, but a sentence clarifying the time from which the weak-coupling evolution is applied would help avoid confusion.","section":"Sec. V.B"},{"comment":"Several entries in Table IV are upper limits (e.g., PSR J0007+7303, Swift J1818–1607); the plotting convention for these limits in Fig. 9 and the treatment of the associated statistical errors should be described.","section":"Appendix B, Table IV"},{"comment":"The term 'weak-coupling' may be misread as referring to weak interactions; a brief definition at first use (reduced collisional coupling between particles) would improve clarity.","section":"Abstract and Sec. I"},{"comment":"The phrase 'indicates, form left to right' should read 'indicates, from left to right'.","section":"Captions of Figs. 4 and 5"}],"recommendation":"major_revision","confidential_remarks":"The paper is well written and technically sound within its stated scope. The single load-bearing issue is the stability of the GS equilibrium; the authors should treat it head-on, either by a stability analysis or by carefully qualifying the abstract and conclusions. The numerical inconsistency in Eq. (65) should be corrected. I believe the paper is publishable after these revisions."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the paper does something genuinely new — it adds thermal evolution to the two-fluid ambipolar diffusion model via a time reparametrization and uses it to make a quantitative claim: ambipolar heating alone cannot explain magnetar quiescent X-ray luminosities. That conclusion is supported within the axisymmetric model, and the numerical work earns credit. The code is validated against semi-analytic results, energy conservation is checked, and the scaling laws are transparent. The novelty over Castillo et al. 2020 is the thermal coupling and the observational comparison, not the constant-temperature equilibria.\n\nThe soft spots are real but for the most part openly stated. The main one is the stability of the Grad-Shafranov equilibrium. The authors acknowledge in Sec. VI.3 that purely poloidal and toroidal fields are unstable in 3D and that barotropic stars may not support stable equilibria, but they do not test stability of their specific configurations. If the axisymmetric endpoint is unstable, continued magnetic dissipation would be stronger and the 'insufficient heating' conclusion could fail. This is not a minor technicality: the predicted luminosity at B=5e15 G (about 6.7e34 erg/s) is only a factor of a few below the brightest observed magnetars, so a modest increase in dissipation could close the gap. The abstract states the conclusion without this caveat, which overstates the current evidence.\n\nOther issues are minor: no code release, no propagated model uncertainties, and the microphysical inputs (e.g., gamma_np) are uncertain by factors of a few. The authors note the gamma_np overestimate and argue it does not change the conclusion; I think that is right because the conclusion is about insufficiency, not a precise luminosity. No circularity: the envelope fits are independent, not calibrated to magnetar data.\n\nWho benefits: magnetar modelers and anyone studying core field evolution. The paper is honest, well-structured, and the central claim is testable. It deserves a serious referee. I would send it to review and ask the authors to either strengthen the stability discussion or soften the abstract accordingly.","headline":"Solid axisymmetric two-fluid simulations that quantify ambipolar heating and show it falls short for magnetars, but the 3D stability caveat means the headline conclusion should be softer.","tokens_in":32576,"tokens_out":3283,"would_cite":true,"duration_ms":25306,"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":"The paper argues that ambipolar diffusion, the slow drift of charged particles through neutrons in a neutron star core, cannot by itself power the persistent X-ray luminosity of magnetars, and shows this with two-fluid axisymmetric…","keywords":["magnetars","ambipolar diffusion","neutron star cores","magnetothermal evolution","quiescent X-ray luminosity","Grad-Shafranov equilibrium","two-fluid magnetohydrodynamics","magnetic field decay"],"falsifier":"Run the same two-fluid model in three dimensions starting from the paper's initial configurations; if the Grad–Shafranov state is destroyed by non-axisymmetric instabilities within $\\sim10^3$ yr and the extra dissipation raises the surface luminosity above roughly $6.7\\times10^{34}\\,\\mathrm{erg\\,s^{-1}}$ for $B=5\\times10^{15}$ G, the conclusion that ambipolar heating alone is insufficient would be overturned.","tokens_in":31494,"feed_emoji":"🌠","tokens_out":6210,"duration_ms":47105,"temperature":0.7,"pith_summary":"The paper asks whether ambipolar diffusion—the slow drift of charged particles through neutrons in a neutron star core—can convert enough magnetic energy into heat to sustain the X-rays observed from magnetars. The authors simulate the long-term magnetic and thermal evolution of a normal, non-superfluid core in axial symmetry, treating it as two coupled fluids and using a time-reparametrization to follow the temperature as it changes. They find that only very strong fields, $B \\gtrsim 5 \\times 10^{15}\\,\\mathrm{G}$, can briefly balance ambipolar heating against neutrino cooling, holding the core hot for about $10^3$ years. Even then, the predicted surface X-ray luminosity stays below the brightest observed magnetars, so the paper concludes that ambipolar diffusion on its own cannot explain their persistent emission.","feed_headline":"Ambipolar diffusion alone cannot explain magnetar X-rays","feed_subtitle":"Even 5×10^15 G fields keep a heating-cooling balance for only ~1,000 years, short of the brightest sources.","key_machinery":"The central machinery is the two-fluid weak-coupling model, in which a neutron fluid and a single charged-particle fluid (protons and electrons) are coupled by collisions while the magnetic field is frozen into the charged fluid; the ambipolar velocity is $v_\\mathrm{ad} = \\mu \\nabla \\chi_n / (\\gamma_{np} n_c)$, with $\\chi_n$ the neutron chemical-potential perturbation. An artificial friction force filters out fast Alfvén waves and enforces non-penetration at the crust–core boundary, and a time-reparametrization $dt = [\\gamma_{np}(t)/\\gamma'_{np}]\\,dt'$ maps constant-temperature simulations onto an evolving temperature because the collision coefficient scales as $\\gamma_{np} \\propto T^2$. This combination lets one simulation at fixed temperature be rescaled to any field strength and thermal history, which is what allows the paper to scan magnetar parameters and compare with observations.","core_discovery":"On its own terms, the paper establishes that in the weak-coupling regime ambipolar diffusion drives a neutron star core toward a Grad–Shafranov equilibrium in which neutrons are in diffusive equilibrium and the Lorentz force is balanced by chemical-potential gradients in the charged-particle fluid. When the thermal evolution is included, the same evolution implies a heating–cooling balance only for $B \\gtrsim 5 \\times 10^{15}\\,\\mathrm{G}$, lasting roughly $10^3\\,[B/(5\\times10^{15}\\,\\mathrm{G})]^{-6/5}\\,\\mathrm{yr}$, with a corresponding surface luminosity $L_\\mathrm{s}^\\infty \\sim 6.7\\times10^{34}\\,(B/5\\times10^{15}\\,\\mathrm{G})^{0.75}\\,\\mathrm{erg\\,s^{-1}}$. Comparing these tracks with a sample of observed cooling neutron stars and magnetars, the paper finds that the brightest quiescent X-ray luminosities cannot be reproduced by ambipolar heating alone.","pith_inferences":["If the Grad–Shafranov equilibrium is unstable to non-axisymmetric perturbations, as purely poloidal and toroidal fields are known to be, magnetic dissipation would be stronger and shorter-lived, so the paper's negative conclusion hinges on that stability.","The same time-reparametrization strategy could be extended to cores with muons, where the charged component is non-barotropic and the final equilibrium is not a Grad–Shafranov state, possibly changing the heating luminosity.","Coupling the core to a realistic crust with finite conductivity would slow the core's magnetic evolution relative to the vacuum-boundary case, making it even harder for ambipolar diffusion alone to match observed luminosities."],"forward_implications":["If the paper is right, ambipolar diffusion in a normal core cannot be the main heating source for magnetar quiescent X-ray emission.","For $B \\gtrsim 5\\times10^{15}\\,\\mathrm{G}$, ambipolar heating delays cooling for about a thousand years before the field settles into Grad–Shafranov equilibrium and heating fades.","Two-fluid models that let neutrons move evolve faster than one-fluid models with fixed neutrons, so published timescales based on the one-fluid assumption need revision.","The neglect of non-equilibrium Urca reactions is justified only up to $B \\sim 5\\times10^{15}\\,\\mathrm{G}$; stronger fields would require a coupled treatment.","Observed luminous magnetars would need crustal dissipation, three-dimensional instabilities, or superfluid and superconducting effects to supply the missing luminosity."],"supporting_citations":[{"why":"Introduced ambipolar diffusion as a heating mechanism for magnetar quiescent luminosity, the scenario this paper tests.","marker":"[3]"},{"why":"Developed the energy-budget argument that strong-field evolution can power magnetar emission.","marker":"[4]"},{"why":"Supplies the numerical scheme, anelastic approximation, and time-reparametrization method extended here from the strong-coupling to the weak-coupling regime.","marker":"[30]"},{"why":"Analytical two-fluid treatment showing neutron bulk motion shortens the evolution timescale and predicting the Grad–Shafranov final state.","marker":"[34]"},{"why":"Earlier axisymmetric two-fluid ambipolar diffusion simulations whose magnetic evolution results are extended here to include thermal evolution.","marker":"[40]"},{"why":"Validates the artificial friction method by comparison with semi-analytical solutions in the limit of negligible friction.","marker":"[44]"},{"why":"One-dimensional thermal-evolution model with ambipolar heating that this two-fluid model improves upon.","marker":"[5]"},{"why":"Compilation of magnetar quiescent luminosities and ages used for the observational comparison.","marker":"[7]"}],"fun_headline_variants":["Ambipolar diffusion falls short for brightest magnetars","Magnetar X-rays exceed ambipolar heating predictions","Core ambipolar flow can't explain magnetar luminosity","Ambipolar diffusion can't heat brightest magnetars enough","Weak-coupling ambipolar diffusion can't power magnetars"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that the axisymmetric magnetic equilibrium reached in the simulations is stable in three dimensions; if that equilibrium breaks apart, dissipation would be stronger and the conclusion could flip.","fun_headline_variants_meta":{"raw":{"variants":["Ambipolar diffusion falls short for brightest magnetars","Magnetar X-rays exceed ambipolar heating predictions","Core ambipolar flow can't explain magnetar luminosity","Ambipolar diffusion can't heat brightest magnetars enough","Weak-coupling ambipolar diffusion can't power magnetars"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000565,"raw_usage":{"total_tokens":2757,"prompt_tokens":1104,"completion_tokens":1653,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":720,"completion_tokens_details":{"reasoning_tokens":1577}},"tokens_in":720,"tokens_out":1653,"duration_ms":9601,"temperature":1.0,"reasoning_tokens":1577,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T14:26:38.494150+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same two-fluid model in three dimensions starting from the paper's initial configurations; if the Grad–Shafranov state is destroyed by non-axisymmetric instabilities within $\\sim10^3$ yr and the extra dissipation raises the surface luminosity above roughly $6.7\\times10^{34}\\,\\mathrm{erg\\,s^{-1}}$ for $B=5\\times10^{15}$ G, the conclusion that ambipolar heating alone is insufficient would be overturned.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Earlier axisymmetric two-fluid ambipolar diffusion simulations whose magnetic evolution results are extended here to include thermal evolution."},{"cited_title":"Reisenegger, A & A 499, 557 (2009)","cited_arxiv_id":null,"evidence_quote":"Validates the artificial friction method by comparison with semi-analytical solutions in the limit of negligible friction."}],"review_version":1}