{"id":"b29fbce7-5392-47d2-8316-af56f8ed79b7","arxiv_id":"2411.14395","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Modern thermodynamically consistent crust models reproduce the observed cooling of MXB 1659-29 and IGR J17480-2446 only if shallow heating and a low-conductivity layer are added, just as traditional models require.","lead":"This paper simulates neutron star crust cooling in two soft X-ray transients using new accreted crust models with realistic nuclear compositions. It finds these modern models fit observations about as well as traditional ones, but only with extra unmodelled ingredients: shallow heating and reduced thermal conductivity.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The SGC-model comparison rests on the untested nHD neutron-diffusion assumption; if diffusion is kinetically suppressed on accretion timescales, the claimed new-model content and the comparison's meaning collapse.","rationale":"I read the paper as a model-comparison exercise: the authors simulate crust cooling with SGC and F+18 models and find both can match observations only with extra ingredients such as shallow heating and reduced conductivity. I agree with the reader's conditional verdict and with the identification of the nHD equilibrium assumption as the weakest load-bearing point. The paper offers no estimate of the neutron diffusion timescale, and the entire distinctiveness of the SGC models rests on that assumption. The K and RP variants already fail against MXB 1659-29, and IGR J17480-2446 requires an added amorphous layer, so the SGC family is not uniformly successful; the conclusion effectively rests on the SB variant plus auxiliary assumptions. I do not see an internal inconsistency or an overreading: the authors explicitly acknowledge the auxiliary ingredients, and the data tables are made available. However, the new-model claim is only as secure as the nHD assumption. Since the reader's conditional verdict already captures this external-validity risk, no change to the verdict is needed.","tokens_in":16136,"tokens_out":5233,"duration_ms":59186,"concrete_test":"Compute the neutron diffusion coefficient in the inner crust from a Boltzmann transport treatment using the SGC compositions and realistic neutron-nucleus and impurity scattering cross sections at T ~ 10^8 K and densities 10^11-10^13 g cm^-3, then integrate the diffusion time across a local pressure scale height. Compare this time with the accretion residence time of a mass layer for the accretion rates used in Section 3 (roughly 10^-11 to 10^-9 solar masses per year). If the diffusion time exceeds the residence time by an order of magnitude in any layer, the nHD equilibrium is not justified and the SGC-based cooling curves are not reliable; if it is shorter, the concern is settled in the authors' favor.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The strongest claim—that SGC and traditional accreted-crust models are similar in explanatory power—depends on the nHD equilibrium condition adopted from Gusakov and Chugunov (2020, 2021) and Shchechilin et al. (2021-2023). Under this condition, free neutrons in the inner crust diffuse rapidly enough to maintain hydrostatic and diffusion equilibrium while accreted matter is advected inward, and the SGC composition, impurity parameter Qimp, and heating profiles are computed on that basis. The present paper does not test this condition; it imports it. If neutron diffusion is actually slow compared with the local accretion residence time, the SGC tables describe an unphysical crust, and the simulation comparison in Section 3 is not a test of the models it claims to compare. The problem is load-bearing rather than cosmetic because the K and RP ash variants already fail against MXB 1659-29, and IGR J17480-2446 requires a hand-placed amorphous layer; the nHD assumption is what gives the SGC family its distinctive composition-dependent content and its claimed status as a 'new' alternative to the traditional picture in which free neutrons move with the nuclei. Without independent support for the diffusion timescale, the central conclusion is conditional on an assumption that is not established anywhere in the manuscript.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper simulates the thermal evolution of the neutron stars in two soft X-ray transients, MXB 1659-29 and IGR J17480-2446, using the new SGC family of accreted-crust models based on realistic thermonuclear ash compositions and the nHD equilibrium condition, and compares the resulting cooling curves with those of the traditional Fantina et al. (2018) (F+18) model and with observed post-outburst temperatures. The main finding is that the SGC models, in particular the Superburst (SB) ash variant, can reproduce the MXB 1659-29 cooling curves only after introducing a shallow-heating source, while all SGC variants require an additional low-conductivity amorphous layer to reproduce IGR J17480-2446. The authors conclude that the new and traditional crust models are similar in their explanatory power and that both need extra ingredients beyond the current deep-crustal-heating theory.","tokens_in":16595,"tokens_out":4569,"duration_ms":51521,"significance":"If the conclusions hold, the paper is a valuable calibration result: it shows that the thermodynamically consistent SGC crust models, despite their more realistic composition and heating profiles, do not remove the long-standing need for shallow heating and for a strongly reduced thermal conductivity in some sources. This is important for future attempts to use crust-cooling observations to constrain dense-matter microphysics, because it narrows the space of model predictions that can be tested against data. The paper also reduces the number of free parameters by taking the impurity parameter Qimp and composition directly from the SGC tables, and it makes the adopted crust tables publicly available, which is a useful contribution. The significance is limited, however, because the key comparisons depend on several fitted quantities (Esh, the outburst-II accretion rate, and the location of the amorphous layer) and because the nHD equilibrium condition underlying all SGC models is imported without an independent test in this manuscript.","major_comments":[{"comment":"The entire SGC family of crust models rests on the nHD equilibrium condition, which assumes that free neutrons in the inner crust diffuse on timescales shorter than the accretion residence time so that the composition, Qimp, and heating profile follow the Gusakov-Chugunov and Shchechilin et al. calculations. The present paper adopts this condition without evaluating its kinetic feasibility. If neutron diffusion is suppressed on accretion timescales, the SGC composition and heating tables describe an unphysical crust, and the comparison in Figs. 3, 5, and 6 would not be a test of the claimed new-model content. Please add a quantitative justification, either an order-of-magnitude estimate of the neutron diffusion timescale versus the advection timescale or a reference to a direct calculation, and state how the conclusions would change if diffusion is slow.","section":"Sections 1 and 2"},{"comment":"The paper states that Esh is adjusted so that each simulated light curve matches the first observation after outburst I, and that Mdot_II is chosen as the value providing the best fits after outburst II. Consequently, the agreement of the SB and F+18 curves with the data is not a parameter-free prediction of the models; it is a test of the shape of the cooling curve after normalizing to one point. This does not by itself invalidate the conclusion that shallow heating is required, but the comparison between model families is weaker than the text implies. Please quantify the sensitivity of the conclusions to the fitted values, for example by showing the range of Esh and Mdot_II that is consistent with the data, and indicate whether the relative ranking of the models is robust within those ranges.","section":"Section 3.1, Fig. 3"},{"comment":"The amorphous layer for IGR J17480-2446 is introduced with an assumed structure factor S=1 (Qimp=<Z>^2) and its location between about 3e10 and 3e11 g cm^-3 is chosen to minimize the discrepancy with the observed cooling curve. This is a genuinely new ad hoc ingredient, and its position is a fitted parameter. The conclusion that the SGC models require a strong conductivity reduction is therefore conditional on this placement and on the independent-scatterer model. Please either provide a physical argument for why an amorphous layer should appear preferentially at that density or show that the inferred need for reduced conductivity is insensitive to moving the layer within a plausible range.","section":"Section 3.2, Fig. 6"},{"comment":"The central claim that 'SGC models are similar to traditional models' is based mainly on the SB ash variant: the K and RP variants are found to be incompatible with the MXB 1659-29 observations (Fig. 3 and the discussion in Section 3.1). Since the three ash models are presented as representative rather than as a weighted set, the conclusion generalizes from one member of the SGC family. Please discuss whether the ash composition appropriate for MXB 1659-29 can be constrained independently, and whether the similarity conclusion would survive if K or RP were the correct ashes for this source.","section":"Section 3.1 and Conclusions"}],"minor_comments":[{"comment":"The adopted accretion rate during outburst II is given as Mdot_II = 1.8e-9 M_sun/yr, but no uncertainty is quoted; given that this value is fitted to the cooling data, an error estimate would help the reader judge the significance of the agreement in the right panel of Fig. 3.","section":"Section 3.1"},{"comment":"The legend uses subscripts on GC, K, RP, and SB to denote Poi values, but the subscript values are small and hard to read in a printed figure; please enlarge the legend fonts or list the Poi values explicitly in the caption.","section":"Figure 3 caption"},{"comment":"There is a typo: 'vertical doted lines' should read 'vertical dotted lines'.","section":"Figure 6 caption"},{"comment":"The phrase 'unlike the past analyses, which treated ... Qimp as free ... fitting parameter' is slightly overstated, because in the present analysis the location and extent of the amorphous layer for IGR J17480-2446 are still effectively free parameters; please rephrase to distinguish the treatment of Qimp from the new fitted layer parameters.","section":"Section 4, first paragraph"},{"comment":"The fitting formulas for P(rho), rho(P), and Ynf are useful, but the paper does not state the accuracy of the Ynf fit in the extrapolated region rho > rho_max; a sentence on the expected uncertainty of the extrapolation would be helpful.","section":"Appendix A"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a careful, modest observational-modeling paper by a group that has made central contributions to this area. The main results are likely to be of interest to the JHEAp readership. The reason for major revision, rather than acceptance, is that the central comparison relies on several fitted ingredients and on the unexamined nHD assumption; none of these issues seems impossible to fix within the scope of the paper, but they need to be addressed explicitly before the similarity claim can be regarded as established. I do not see grounds for rejection: the claims are proportionate, the data are real, and the available tables are a useful resource."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nThe short version: this is a careful, useful consolidation paper. It takes the new SGC accreted-crust models (Kepler, Extreme rp, Superburst ash compositions under the nHD equilibrium condition) and runs them through outburst/cooling simulations for MXB 1659-29 and IGR J17480-2446, comparing against the traditional Fantina+18 model. The main result is modest and credible: the new models do not outperform the old ones, and both need extra ingredients—shallow heating and, for IGR, a low-conductivity layer—to match the light curves. That is a genuinely useful negative result for the neutron-star thermal evolution community.\n\nWhat is new: this is the first application of the SGC models to full outburst cooling curves, and the paper ships analytical fits to the inner crust EoS and public tables for composition and heating. That is reproducible, and the fits are honestly characterized (1–2% pressure, 3.5% density inversion). The consistent treatment of MXB 1659-29 through two outbursts with fixed parameters is better than much of the literature.\n\nSoft spots, in proportion. The shallow heating energy Esh is tuned to the first cooling point for each model, the accretion rate for outburst II is chosen to give the best fit, and the amorphous layer in IGR is placed to minimize discrepancy. These are fitting choices, not free-floating inventions—the paper is transparent about them—but they mean the comparison is conditional on those choices. The larger caveat is the nHD assumption itself: the SGC models assume neutron diffusion keeps the inner crust in hydrostatic and diffusion equilibrium during accretion. That assumption is imported from the authors' prior work and is not tested here. If neutron diffusion is kinetically suppressed on accretion timescales, the SGC tables describe an unphysical crust and the comparison loses its new-model content. I don't think that kills the paper—it is clearly an application paper, and the authors flag the conditional nature of their models—but it means the central conclusion is 'given nHD, SGC and traditional models are similar,' not an unconditional statement.\n\nWho this is for: anyone working on neutron star crust cooling, SXT light curves, or accreted crust microphysics. It deserves a serious referee: the modeling is careful, the data products are useful, and the conclusion, while modest, is exactly the kind of consolidation that moves the field forward. I'd recommend engaging with it, not desk-rejecting it.","headline":"Careful consolidation: SGC crust models behave like traditional ones and still need shallow heating; the imported nHD assumption is the load-bearing caveat.","tokens_in":16998,"tokens_out":2415,"would_cite":true,"duration_ms":22965,"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 new nHD-equilibrium crust models fit observed transient cooling as well as the traditional model, and both require unexplained extras (shallow heating and reduced conductivity) to match the data.","keywords":["neutron stars","dense matter","soft X-ray transients","crust cooling","accreted neutron star crust","nHD equilibrium","shallow heating","thermal conductivity"],"falsifier":"One decisive test is a nuclear-physics estimate of the diffusion coefficient of free neutrons in the inner crust at densities of order $10^{11}$ to $10^{12}$ g/$cm^{3}$; if the resulting equilibration time is longer than a typical accretion outburst, the nHD equilibrium condition on which all the new models rest is not satisfied, and the paper's comparison would no longer test a physically realized crust.","tokens_in":15946,"feed_emoji":"⭐","tokens_out":12387,"duration_ms":105372,"temperature":0.7,"pith_summary":"Soft X-ray transients are neutron stars that are episodically heated by accretion and then cool in quiescence; their cooling curves are one of the few direct probes of the structure and composition of the neutron-star crust. This paper tests whether a new family of accreted-crust models, built on the assumption that free neutrons diffuse into hydrostatic and diffusion equilibrium during accretion, can reproduce the observed post-outburst cooling of the two best-monitored sources, MXB 1659-29 and IGR J17480-2446. The result is that the new models and the traditional 2018 model are similar in their ability to fit the data: the new models still need an extra shallow heating source to explain MXB 1659-29, and they still need a strongly suppressed thermal conductivity in a thin outer-crust layer to explain the slow cooling of IGR J17480-2446. If the paper is right, the cooling observations currently in hand do not distinguish between the two crust-model families, and the extra physics those observations demand is real and still unidentified.","feed_headline":"New crust models still need extra heating to fit neutron-star cooling","feed_subtitle":"Both new and traditional crust models need shallow heating and lower conductivity to match transients.","key_machinery":"The load-bearing construction is the nHD equilibrium condition: during accretion, free neutrons in the inner crust redistribute by diffusion until the crust reaches hydrostatic and diffusion equilibrium. This condition fixes the composition, the impurity parameter Qimp (the charge variance that scatters electrons and reduces thermal conductivity), and the heat released per accreted baryon, so these inputs are no longer free parameters but outputs of a nuclear reaction-network calculation. The paper combines three nHD-based crust models, corresponding to three representative thermonuclear ash compositions, with analytical fits for the inner-crust equation of state and a stellar thermal-evolution code, and compares the resulting cooling curves with the observed light curves of the two transients. The nHD assumption is the essential difference from the traditional model, in which free neutrons are assumed to move together with the nuclei.","core_discovery":"The paper's central finding is a comparison rather than a new mechanism: when the three new nHD-equilibrium crust models (the paper's K, RP, and SB ash versions) are used to simulate the same two transients, their cooling curves fall in the same range as those of the traditional 2018 model. For MXB 1659-29, the SB version with the lowest allowed outer-inner crust pressure matches the observed light curves if a shallow heating source of about 0.85 MeV per baryon is added; the K and RP versions cool too slowly because their impurity parameter is large, and thus do not fit. For IGR J17480-2446, none of the new models matches the observed slow cooling unless a layer near the bottom of the outer crust is treated as amorphous, with the thermal conductivity suppressed to the level of independent ion scattering. The total deep crustal heating in the new models is only about 0.2-0.6 MeV per accreted baryon, compared to about 1.5-2 MeV in the traditional model, yet the extra ingredients required to fit observations are the same for both families.","pith_inferences":["Beyond the paper, the fitted 'shallow heating' energy of about 0.85-1 MeV per baryon is comparable in size to the deep crustal heating itself, which suggests that whatever mechanism produces it, possibly nuclear reactions in the outer envelope rather than crust physics, is a first-order ingredient for all future transient-cooling fits.","The paper's need for a low-conductivity layer in one source but not the other hints that the conductivity suppression may be source-dependent, perhaps tied to accretion history or magnetic field; that is testable by extending the same fits to other transients with well-measured cooling curves.","A Bayesian model-comparison over a larger sample of transients could quantify how strongly the data prefer the new nHD-equilibrium models over the traditional model, something the two-source comparison here leaves open.","If the nHD equilibrium is correct, then the extra 'shallow heating' term may be absorbing a real deficit in the deep crustal heating of the new models; resolving that would require a microphysical calculation rather than further tuning of the heating parameter."],"forward_implications":["If the paper's result holds, the Superburst-ash nHD model can describe the MXB 1659-29 cooling data, but only with an added shallow heating of roughly 0.85 MeV per accreted baryon.","The K and RP versions of the new models are disfavored for MXB 1659-29 because their high impurity parameter makes the crust cool too slowly.","For IGR J17480-2446, the new models predict faster cooling than observed unless a layer near the bottom of the outer crust is assigned a strongly suppressed thermal conductivity, modeled as an amorphous solid.","The observed cooling therefore does not yet choose between the new nHD-equilibrium crust models and the traditional model; both need the same kinds of extra physics.","The much smaller deep crustal heating in the new models (about 0.2-0.6 MeV per baryon versus 1.5-2 MeV) does not, by itself, break the degeneracy, because the required shallow heating compensates in the quasi-steady quiescent state."],"supporting_citations":[{"why":"Supplies the nHD equilibrium condition and the universal heating-efficiency formula from which the new crust models are built.","marker":"Gusakov and Chugunov (2020, 2021)"},{"why":"Provides the detailed composition, equation of state, impurity parameter, and heat release for the three nHD-based crust models used here.","marker":"Shchechilin et al. (2021, 2022, 2023)"},{"why":"Defines the traditional accreted-crust model that serves as the baseline comparison.","marker":"Fantina et al. (2018)"},{"why":"The authors' earlier study that set up the simulation procedure and the shallow-heating treatment for these transients.","marker":"Potekhin et al. (2023)"},{"why":"Supplies the consistent observational light curves and outburst histories for MXB 1659-29.","marker":"Parikh et al. (2019)"},{"why":"Provides the observed cooling curve and the pre-outburst quasi-equilibrium temperature for IGR J17480-2446.","marker":"Ootes et al. (2019)"},{"why":"Fixes the input parameters, including accretion rate and stellar mass, for the MXB 1659-29 simulations.","marker":"Potekhin and Chabrier (2021)"},{"why":"Introduced the shallow-heating requirement that the paper still needs to fit observations.","marker":"Brown and Cumming (2009)"}],"fun_headline_variants":["Realistic crust models still require shallow heating to fit data","New crust physics doesn't eliminate need for shallow heating","Crust models: new and old alike need extra ingredients to match SXTs","Neutron star cooling: new crust models fall short without extra heat"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The new models all assume that free neutrons diffuse rapidly enough in the inner crust during accretion to reach hydrostatic and diffusion equilibrium; if that diffusion is too slow, the predicted compositions, impurity parameters, and heating profiles do not apply.","fun_headline_variants_meta":{"raw":{"variants":["Realistic crust models still require shallow heating to fit data","New crust physics doesn't eliminate need for shallow heating","Crust models: new and old alike need extra ingredients to match SXTs","Neutron star cooling: new crust models fall short without extra heat"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000338,"raw_usage":{"total_tokens":1885,"prompt_tokens":979,"completion_tokens":906,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":595,"completion_tokens_details":{"reasoning_tokens":832}},"tokens_in":595,"tokens_out":906,"duration_ms":9205,"temperature":1.0,"reasoning_tokens":832,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T15:13:09.763560+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"One decisive test is a nuclear-physics estimate of the diffusion coefficient of free neutrons in the inner crust at densities of order $10^{11}$ to $10^{12}$ g/$cm^{3}$; if the resulting equilibration time is longer than a typical accretion outburst, the nHD equilibrium condition on which all the new models rest is not satisfied, and the paper's comparison would no longer test a physically realized crust.","supporting_citations":[{"cited_title":"Crustal heating in accreting neutron stars from the nuclear energy-density functional theory","cited_arxiv_id":null,"evidence_quote":"Defines the traditional accreted-crust model that serves as the baseline comparison."},{"cited_title":"Thermodynamically consistent equation of state for an accreted neutron star crust","cited_arxiv_id":null,"evidence_quote":"Supplies the nHD equilibrium condition and the universal heating-efficiency formula from which the new crust models are built."},{"cited_title":"https://ui.adsabs.harvard.edu/abs/2021MNRAS.507.3860S/abstract , 507, 3860","cited_arxiv_id":null,"evidence_quote":"Provides the detailed composition, equation of state, impurity parameter, and heat release for the three nHD-based crust models used here."},{"cited_title":"Thermal evolution of neutron stars in soft X-ray transients with thermodynamically consistent models of the accreted crust","cited_arxiv_id":null,"evidence_quote":"The authors' earlier study that set up the simulation procedure and the shallow-heating treatment for these transients."},{"cited_title":"Consistent accretion-induced heating of the neutron-star crust in MXB 1659--29 during two different outbursts","cited_arxiv_id":null,"evidence_quote":"Supplies the consistent observational light curves and outburst histories for MXB 1659-29."},{"cited_title":"Continued cooling of the accretion-heated neutron star crust in the X-ray transient IGR J17480 - 2446 located in the globular cluster Terzan 5","cited_arxiv_id":null,"evidence_quote":"Provides the observed cooling curve and the pre-outburst quasi-equilibrium temperature for IGR J17480-2446."},{"cited_title":"Crust structure and thermal evolution of neutron stars in soft X-ray transients","cited_arxiv_id":null,"evidence_quote":"Fixes the input parameters, including accretion rate and stellar mass, for the MXB 1659-29 simulations."},{"cited_title":"Mapping crustal heating with the cooling light curves of quasi-persistent transients","cited_arxiv_id":null,"evidence_quote":"Introduced the shallow-heating requirement that the paper still needs to fit observations."}],"review_version":1}