{"id":"d1e9f778-9312-44a4-a4dd-a010ae89f607","arxiv_id":"2506.09841","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Strong poloidal magnetic fields deform neutron stars and reduce central density, suppressing the direct Urca process and slowing thermal cooling.","lead":"Neutron star cooling simulations show that extremely strong internal magnetic fields can slow down cooling by shutting off the fast direct Urca neutrino process. The effect appears even though the magnetic field is too weak to change the local equation of state, because the field changes the star's geometry and central density.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The fast-to-slow cooling contrast hinges on the 1.4 Msun CMF star being above the direct Urca threshold and on omitting pairing; if either fails, the magnetic-field effect vanishes.","rationale":"The reader identified precisely the same load-bearing assumption: the result depends on the CMF EoS opening direct Urca at 1.4 Msun, and on the absence of pairing. My independent reading confirms this is the weakest link. The paper provides a plausible mechanism and internally consistent calculations, but the headline effect is conditional on two microphysical inputs that are neither demonstrated nor varied. Since the reader's verdict is already CONDITIONAL with this exact caveat, my stress-test does not change the verdict; it reinforces the need for the authors to verify the DU threshold in their EoS and to test at least one pairing scenario. The concern is not an internal inconsistency; it is a model-dependence that the paper acknowledges in its limitation statement, and the proposed test would settle it directly.","tokens_in":9503,"tokens_out":8594,"duration_ms":109764,"concrete_test":"Run the identical 2D cooling calculation for: (a) the same 1.4 Msun CMF star but with a standard nucleon pairing model (e.g., S0 neutron and P2 proton gaps); and (b) an alternative EoS (e.g., SLy4 or GM1 with and without hyperons) that has DU closed at 1.4 Msun. If including pairing, or switching EoS, removes the fast-to-slow transition with increasing magnetic field, the central claim is confirmed to be model- and pairing-dependent. Additionally, publish the proton fraction profile versus density for the zero-field CMF star to verify the DU threshold crossing quantitatively.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim rests on the assertion that a 1.4 solar-mass star built with the CMF EoS permits the direct Urca (DU) process in the absence of a magnetic field, so that the field-induced reduction of the DU region produces the observed transition from fast to slow cooling. This is stated in Sec. 3 but not demonstrated with a composition profile, and it is a single-model, single-mass point. The same section explicitly defers nucleon pairing to future work, which matters because in realistic neutron stars pairing suppresses DU neutrino emission for much of the cooling epoch. If the CMF EoS does not actually cross the DU threshold at 1.4 Msun, or if pairing is included, the zero-field star would already cool slowly and the magnetic field would no longer cause a fast-to-slow transition. The paper's own limitation statement ('This is usually alleviated by the inclusion of appropriate pairing among nucleons') confirms the centrality of this assumption. Without a demonstration of DU being open in the unpaired, unmagnetized CMF star and a check of how pairing modifies the result, the quantitative cooling curves and the non-linear relaxation-time increase are not robust beyond this specific model.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the thermal evolution of non-rotating, axisymmetric neutron stars with strong poloidal magnetic fields, using the CMF equation of state for the microphysics and the Astreus code to solve the coupled Einstein-Maxwell equations for stellar structure and a 2D finite-difference scheme for the cooling equations. For a canonical 1.4 solar-mass star, the authors construct configurations with central fields up to about 4.7e17 G, solve the cooling evolution including standard neutrino emission processes, and report that stars with the strongest fields cool more slowly than their unmagnetized counterpart. They attribute this to the magnetic-field-induced reduction in the volume where the direct Urca process operates, and they report a nonlinear increase in the thermal relaxation time with increasing field strength, fitting it with a four-parameter rational function. The central claim is that magnetic fields can alter cooling behavior even when too weak to modify the Fermi distributions directly, through the effect of curvature on the stellar composition and the DU region.","tokens_in":9766,"tokens_out":2962,"duration_ms":41224,"significance":"If the result is robust, it is a novel and interesting proof-of-principle: magnetic fields of the order of 1e17 G could change the thermal evolution of neutron stars by modifying the active direct Urca region through general-relativistic deformation, rather than through a direct effect on the equation of state or on single-particle phase space. The work combines a nontrivial 2D general-relativistic magnetized equilibrium solver with a 2D cooling solver, which is a technical achievement. The paper also provides a concrete falsifiable prediction: for a fixed gravitational mass, stronger magnetic fields produce higher surface temperatures at ages of 1e2-1e6 yr and a longer thermal relaxation time. However, the claim currently rests on a single equation of state, a single mass, an unpaired treatment of nucleons, and a model-specific DU threshold, so its quantitative predictions should be treated cautiously until those dependencies are characterized.","major_comments":[{"comment":"The central mechanism requires that the unmagnetized 1.4 solar-mass CMF star permit the direct Urca process, but the paper does not demonstrate this. The sentence \"our microscopic model allows for the DU process in stars with 1.4 M_sun\" is asserted without showing the beta-equilibrium proton fraction (or any composition profile) versus density and its relation to the DU threshold. Figure 4 shows DU regions only for magnetized configurations, not for f0 = 0. Without a field-free composition profile demonstrating that the DU threshold is actually crossed, the identification of the fast-cooling branch with direct Urca is unsupported and the central fast-to-slow transition could be an artifact of the particular EoS or the chosen mass.","section":"Sec. 3"},{"comment":"The omission of nucleon pairing is not a minor simplification here because the paper's stated effect is the suppression of DU: the text even notes that \"This is usually alleviated by the inclusion of appropriate pairing among nucleons.\" In realistic neutron stars, pairing gaps suppress DU neutrino emission over much of the cooling epoch, so including pairing could substantially reduce or even eliminate the contrast between low-field and high-field cooling curves shown in Fig. 2. The authors should either include a pairing model in at least one representative run, or quantitatively argue (e.g., via the relevant temperature range and expected gap magnitudes) that pairing would not erase the reported fast-to-slow transition.","section":"Sec. 3, Fig. 2"},{"comment":"The manuscript reports quantitative cooling curves and relaxation times but provides no numerical convergence tests, resolution studies, or uncertainty estimates. The claim of a nonlinear increase in relaxation time, and the fitted parameters of Eq. (7), are only meaningful if the 2D cooling solutions are converged with respect to grid spacing and time step. The paper should show at least one convergence test for a representative strong-field configuration and, if possible, estimates of the numerical error in the extracted relaxation times.","section":"Sec. 3, Figs. 2-3"}],"minor_comments":[{"comment":"The definition of the relaxation time t_w = max |d ln T_s / d ln t| is not explained: as written the expression is dimensionless, while t_w is presented in years. Please clarify how this quantity is extracted from the cooling curves and what physical condition it represents.","section":"Eq. (6) and Fig. 3"},{"comment":"The four-parameter fit to six data points is presented as \"a good fit\" with a 95% confidence band in Fig. 3, but no residuals, reduced chi-squared, or parameter uncertainties are given. Since this fit is illustrative rather than load-bearing, a brief statement of the fit quality would suffice.","section":"Sec. 3, Eq. (7)"},{"comment":"The caption says the magnetic moment is given \"in Gaussians, where 1 Gaussian = 1e-3 A m^2\"; the unit and conversion should be defined consistently, and the table would benefit from listing the gravitational mass explicitly for each configuration.","section":"Table 1"},{"comment":"The reference \"J Zapata, R Negreiros, T. S., & Jaikumar, P. 2022\" has garbled author formatting and should be corrected.","section":"References"},{"comment":"The introduction refers to conclusions in Sec. 5, but the conclusions section is numbered 4; renumber or adjust the cross-reference.","section":"Sec. 4 vs. text"},{"comment":"The DU-active regions would be much easier to interpret if the figure also showed the f0 = 0 case and if the color scale for the yellow region were accompanied by a quantitative definition of the DU criterion (e.g., the proton-fraction threshold as a function of density).","section":"Fig. 4"}],"recommendation":"major_revision","confidential_remarks":"The paper presents an interesting mechanism and a credible technical setup, but the central claim currently depends on the unpaired CMF model allowing DU in the field-free 1.4 Msun star, which is asserted rather than demonstrated. The pairing caveat is acknowledged in the text but its potential to erase the reported effect is not quantified. I would be comfortable with publication after the authors show the composition profile, include a pairing sensitivity check, and add convergence information; these are load-bearing but fixable within the manuscript's scope."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new thing here is that the authors actually do it: full 2D general-relativistic magnetized structure matched to 2D cooling evolution. They show a clean mechanism—the magnetic field adds to the gravitational mass, so a star with fixed gravitational mass has lower baryonic density and a smaller direct Urca (DU) region. That makes a 1.4 solar-mass star switch from fast to slow cooling and produces a non-linear increase in relaxation time. That calculation did not exist before, and the paper explains it clearly.\n\nCredit where due: the equations in Sec. 2 are standard but correctly assembled; the cooling code is the authors' own and used consistently; Fig. 4 is a direct, convincing illustration of the DU region shrinking with field strength. The four-parameter fit to the relaxation time is presented as an illustration and does not feed back into the cooling calculation, so I don't hold that against them.\n\nThe soft spots are real but not fatal. The load-bearing assumption is that the CMF equation of state at 1.4 solar masses has DU open in the unmagnetized star. The paper says it does, but it never shows the proton fraction vs. density curve or where the DU threshold is crossed. That is a one-figure check the referee should request. Second, pairing is explicitly deferred to future work. Their comment \"Regardless of pairing or lack thereof\" is too quick: if pairing suppresses DU in the field-free star, the zero-field star would already cool slowly and the magnetic field would no longer cause a fast-to-slow transition. So the quantitative cooling curves are not robust beyond this unpaired, single-EoS, single-mass model. Third, the surface fields of 7-8 x 10^16 G are above most observational estimates for magnetars; citing a 10^16 G central estimate and calling the fields \"not far from observed values\" overstates the case. I also would have liked convergence tests or a code release, though that is a minor issue for a first study.\n\nThis is a useful proof-of-principle for the neutron-star cooling community, which still mostly works in 1D. It deserves a serious referee rather than a desk rejection, but the referee should ask for the DU threshold demonstration and a sensitivity check with pairing. I'd sent it to review.","headline":"First 2D GR cooling of magnetized neutron stars with a clear DU-suppression mechanism, but the central fast-to-slow contrast rests on an undemonstrated, single-EoS, unpaired assumption.","tokens_in":10226,"tokens_out":2199,"would_cite":true,"duration_ms":27723,"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":"Magnetic fields near 10^17 G can shut down the fastest neutrino-cooling channel in neutron stars.","keywords":["neutron star cooling","magnetic fields","direct Urca process","axisymmetric general relativity","thermal relaxation time","equation of state","neutrino emission"],"falsifier":"A decisive check is to run the same two-dimensional cooling calculation with an equation of state for which a 1.4-solar-mass star never reaches the direct Urca proton fraction: if the cooling curves then barely change with magnetic field strength, the paper's central claim is falsified. Alternatively, observing a 1.4-solar-mass neutron star with an inferred surface field of $7$–$8\\times10^{16}$ G that cools as fast as unmagnetized stars would contradict the prediction.","tokens_in":9323,"feed_emoji":"🧲","tokens_out":7780,"duration_ms":67479,"temperature":0.7,"pith_summary":"Neutron stars of the same mass can cool very differently depending on how strong their internal magnetic field is. The paper shows that central fields of $3$–$4\\times10^{17}$ G — about ten times stronger than typical surface fields inferred from observations — are enough to shut down the direct Urca process, the fastest neutrino cooling channel, even though such fields are too weak to change the equation of state or the Fermi distribution of particles. The effect works because the magnetic field adds to the gravitational mass, so a star with a fixed observed mass and a stronger field has a lower baryonic density and therefore a smaller core region meeting the conditions for direct Urca. The result matters because it turns cooling curves into a possible probe of the interior magnetic field of neutron stars.","feed_headline":"Magnetic fields can flip neutron-star cooling from fast to slow","feed_subtitle":"Same-mass stars could stay hotter for centuries, turning cooling curves into a probe of internal fields.","key_machinery":"The central mechanism is the coupling between Einstein's and Maxwell's equations for an axisymmetric poloidal magnetic field, solved together with a thermal evolution equation for temperature and heat flux. The named pieces are the chiral mean field (CMF) equation of state, which fixes the proton fraction and therefore whether direct Urca — the fast neutrino-emitting process requiring a proton fraction above a threshold — can run, and the direct Urca threshold itself. The workhorse relation is that the gravitational mass of a magnetized star receives a contribution from the electromagnetic energy, so two stars with the same gravitational mass but different field strengths have different baryonic densities; a stronger field lowers the central density and shrinks the direct Urca region. Thermal relaxation time, defined from the maximum slope of the cooling curve, then grows non-linearly with the field strength, well fitted by a rational function of the current function $f_0$.","core_discovery":"Cooling a 1.4-solar-mass neutron star with a poloidal magnetic field in full general relativity, the authors find that stars with central fields around $3$–$4\\times10^{17}$ G (surface fields $7$–$8\\times10^{16}$ G) remain significantly hotter than their unmagnetized counterparts. The reason is geometric rather than microscopic: the electromagnetic field contributes to curvature and hence to the gravitational mass, so a star with fixed gravitational mass and larger field has lower baryonic content and lower central baryon density. This shrinks the region where the proton fraction exceeds the direct Urca threshold, converting a fast-cooling star into a slow-cooling one and increasing the thermal relaxation time non-linearly with field strength. The transition is seen as a change in cooling regime between current-function values $f_0 = 2.0$ and $f_0 = 2.5$, with the direct-Urca-active volume losing its spherical shape and shrinking as the field grows.","pith_inferences":["Inference: the same mechanism should be mass-selective; the cooling contrast induced by magnetic fields should be largest for stars whose field-free proton fraction sits near the direct Urca threshold, and smaller for stars far above or below it.","Inference: if nucleon pairing is added, absolute temperatures will change, but the geometry-driven shrinkage of the direct Urca region should survive because it is tied to baryon density, not to the pairing gap.","Inference: the ellipsoidal direct Urca region implies anisotropic neutrino emission that may also produce a small aspherical momentum kick, a testable connection to neutron star natal kicks that the paper does not pursue.","Inference: a straightforward extension would be to map the same calculation at several masses and compare the resulting cooling curves against a sample of thermally emitting neutron stars with well-measured masses; a clean separation by inferred surface field would test the claim observationally."],"forward_implications":["Magnetized neutron stars with central fields above roughly $3\\times10^{17}$ G will cool slowly and stay hotter for longer than field-free stars of the same mass, producing distinct cooling curves after about 100 years of age.","The thermal relaxation time increases non-linearly with magnetic field strength, with a fast-to-slow cooling transition around central fields of $3$–$4\\times10^{17}$ G.","The direct Urca active region shrinks and becomes ellipsoidal as the field grows, so neutrino emission becomes spatially anisotropic inside the star.","Observed surface temperatures of neutron stars, combined with an independent mass measurement, can in principle distinguish stars with strong internal fields from those without.","Because the field effect operates through baryon density rather than through particle microphysics, it persists at field strengths too low to alter the equation of state."],"supporting_citations":[{"why":"Supplies the chiral mean field equation of state that fixes the proton fraction and determines where direct Urca can run in the 1.4-solar-mass star.","marker":"Dexheimer & Schramm 2008"},{"why":"Provides the general-relativistic structure equations for magnetized neutron stars, including the hydrostatic equilibrium equation with the Lorentz force used in the stellar models.","marker":"Cardall et al. 2001b"},{"why":"Establishes that poloidal magnetic fields in neutron stars intensify from the surface toward the core, motivating the field strengths considered here.","marker":"Bocquet et al. 1995a"},{"why":"Introduces the Green's function expansion method used by the Astreus code to solve the field equations with flat-space boundary conditions.","marker":"Komatsu et al. 1989"},{"why":"Expands the Green's function technique for relativistic stellar models, underpinning the numerical construction of the magnetized stars.","marker":"Cook et al. 1992"},{"why":"Derives the axisymmetric thermal evolution equations and the alternating-direction implicit integration method adopted in the cooling calculation.","marker":"Negreiros et al. 2012"},{"why":"Catalogs the neutrino emission processes, including direct and modified Urca, used in the cooling calculation.","marker":"Yakovlev & Pethick 2004"},{"why":"Establishes the connection between the non-linear increase of the relaxation time and the transition from fast to slow cooling.","marker":"Sales et al. 2020"},{"why":"Shows that magnetic fields around $10^{17}$ G do not alter the equation of state, supporting the claim that the cooling effect here is geometric rather than microscopic.","marker":"Chatterjee et al. 2015"},{"why":"Shows magnetic fields affect crust thickness, one of the channels through which the field changes the relaxation time.","marker":"Franzon et al. 2017"}],"fun_headline_variants":["Magnetic fields suppress fast cooling in neutron stars","Strong magnetic fields make neutron stars cool slower","Magnetic fields flip neutron-star cooling regime","Neutron stars with strong fields stay hotter longer"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that a 1.4-solar-mass neutron star built from the chosen equation of state would, without a magnetic field, be just dense enough in its core to run direct Urca; if that threshold were not crossed, or if pairing suppressed the process, the magnetic field could not cause the reported fast-to-slow cooling transition.","fun_headline_variants_meta":{"raw":{"variants":["Magnetic fields suppress fast cooling in neutron stars","Strong magnetic fields make neutron stars cool slower","Magnetic fields flip neutron-star cooling regime","Neutron stars with strong fields stay hotter longer"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00016,"raw_usage":{"total_tokens":1206,"prompt_tokens":895,"completion_tokens":311,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":511,"completion_tokens_details":{"reasoning_tokens":253}},"tokens_in":511,"tokens_out":311,"duration_ms":3498,"temperature":1.0,"reasoning_tokens":253,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T04:39:24.944046+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A decisive check is to run the same two-dimensional cooling calculation with an equation of state for which a 1.4-solar-mass star never reaches the direct Urca proton fraction: if the cooling curves then barely change with magnetic field strength, the paper's central claim is falsified. Alternatively, observing a 1.4-solar-mass neutron star with an inferred surface field of $7$–$8\\times10^{16}$ G that cools as fast as unmagnetized stars would contradict the prediction.","supporting_citations":[{"cited_title":"2012, Physical Review D—Particles, Fields, Gravitation, and Cosmology, 85, 104019","cited_arxiv_id":null,"evidence_quote":"Derives the axisymmetric thermal evolution equations and the alternating-direction implicit integration method adopted in the cooling calculation."},{"cited_title":"2020, Astronomy & Astrophysics, 642, A42","cited_arxiv_id":null,"evidence_quote":"Establishes the connection between the non-linear increase of the relaxation time and the transition from fast to slow cooling."},{"cited_title":"2017, Physical Review D, 96, 123005","cited_arxiv_id":null,"evidence_quote":"Shows magnetic fields affect crust thickness, one of the channels through which the field changes the relaxation time."}],"review_version":1}