Crust cooling of the neutron star in Aql X-1: Different depth and magnitude of shallow heating during similar accretion outbursts
Pith reviewed 2026-05-24 20:33 UTC · model grok-4.3
The pith
Shallow heating in the crust of Aql X-1 was deeper and stronger during the 2016 outburst than during the similar 2013 outburst.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Observations after the 2016 outburst reveal a neutron star that was initially cooler and began cooling later than after the 2013 outburst. Thermal simulations of the crust require the shallow heating to have occurred at greater depth and with larger magnitude in 2016 to match the data, even though the outbursts were comparable in duration, fluence, and morphology. This implies that constant neutron star parameters play only a weak role in controlling shallow heating and that similar accretion events can drive different heating outcomes.
What carries the argument
Shallow heating, an unidentified heat source located in the outer layers of the neutron star crust, whose depth and total energy input are constrained by matching observed post-outburst cooling curves to thermal evolution simulations.
If this is right
- Basic neutron star parameters that remain fixed between outbursts do not strongly determine the properties of shallow heating.
- Outbursts with similar accretion morphology can produce markedly different shallow heating in the crust.
- The origin of shallow heating is likely linked to variable aspects of the accretion flow or changes in the crust state rather than to fixed source properties.
- Alternative explanations for the different quiescent evolution, such as changes in the crust between outbursts, remain possible but require further data to test.
Where Pith is reading between the lines
- If shallow heating varies with outburst details not captured by current models, repeated observations of the same source after multiple events could map which accretion parameters control the heating.
- This variability suggests that population studies of neutron star cooling in LMXBs may need to account for outburst-to-outburst differences rather than assuming uniform heating per source.
- Testing whether the same pattern appears in other frequently outbursting sources would clarify whether the effect is common or specific to Aql X-1.
Load-bearing premise
The two outbursts were similar enough in duration, total mass accreted, and accreted material composition that differences in cooling can be attributed solely to changes in the shallow heating parameters.
What would settle it
A measurement demonstrating that the 2013 and 2016 outbursts differed substantially in total accreted mass or duration, or a simulation showing that the observed cooling difference can be reproduced without altering shallow heating parameters.
read the original abstract
The structure and composition of the crust of neutron stars plays an important role in their thermal and magnetic evolution, hence in setting their observational properties. One way to study the crust properties is to measure how it cools after it has been heated during an accretion outburst in a low-mass X-ray binary (LMXB). Such studies have shown that there is a tantalizing source of heat, of currently unknown origin, that is located in the outer layers of the crust and has a strength that varies between different sources and different outbursts. With the aim of understanding the mechanism behind this "shallow heating", we present Chandra and Swift observations of the neutron star LMXB Aql X-1, obtained after its bright 2016 outburst. We find that the neutron star temperature was initially much lower, and started to decrease at much later time, than observed after the 2013 outburst of the source, despite the fact that the properties of the two outbursts were very similar. Comparing our data to thermal evolution simulations, we infer that the depth and magnitude of shallow heating must have been much larger during the 2016 outburst than during the 2013 one. This implies that basic neutron star parameters that do not change between outbursts, do not play a strong role in shallow heating. Furthermore, it suggests that outbursts with a similar accretion morphology can give rise to very different shallow heating. We also discuss alternative explanations for the difference in quiescent evolution after the 2016 outburst.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper presents new Chandra and Swift X-ray observations of Aql X-1 after its 2016 outburst and compares the neutron-star cooling curve to that observed after the 2013 outburst. Despite asserting that the two outbursts had very similar properties, the post-outburst temperature evolution differs markedly; thermal-evolution simulations are used to infer that both the depth and the magnitude of shallow heating were substantially larger in 2016 than in 2013. Alternative explanations are briefly discussed.
Significance. If the central inference is robust, the result shows that shallow heating can change dramatically between accretion episodes whose macroscopic properties appear comparable, thereby ruling out explanations that depend only on unchanging neutron-star parameters such as mass, radius or core temperature. This strengthens the empirical case that shallow heating is tied to variable aspects of the accretion flow or crust composition and supplies a concrete observational constraint for theoretical models of the heating mechanism.
major comments (2)
- [Introduction and observations sections] The claim that the 2013 and 2016 outbursts are sufficiently similar to isolate changes in shallow heating (abstract and introduction) is load-bearing for the central conclusion, yet the manuscript supplies no quantitative comparison of total fluence, integrated mass-accretion rate, outburst duration, or recurrence-time constraints. Without these metrics it is impossible to assess whether residual differences in accretion history could produce the observed cooling difference.
- [Modeling section] The thermal-evolution modeling (modeling section) that yields the new shallow-heating parameters is presented without tabulated best-fit values, uncertainties, or goodness-of-fit statistics for the 2016 data set; the abstract likewise omits these quantities. This prevents evaluation of whether the inferred change in heating depth and strength is statistically required by the data.
minor comments (2)
- The notation used for the shallow-heating parameters (Q_sh and depth) should be defined explicitly the first time it appears and kept consistent with the simulation code output.
- Figure captions should state the energy bands and instruments used for each data point so that the reader can immediately assess the contribution of Swift versus Chandra measurements.
Simulated Author's Rebuttal
We thank the referee for their thoughtful and constructive report. The two major comments identify areas where the manuscript can be strengthened with additional quantitative detail. We address each point below and will revise the paper accordingly.
read point-by-point responses
-
Referee: [Introduction and observations sections] The claim that the 2013 and 2016 outbursts are sufficiently similar to isolate changes in shallow heating (abstract and introduction) is load-bearing for the central conclusion, yet the manuscript supplies no quantitative comparison of total fluence, integrated mass-accretion rate, outburst duration, or recurrence-time constraints. Without these metrics it is impossible to assess whether residual differences in accretion history could produce the observed cooling difference.
Authors: We agree that explicit quantitative metrics would allow readers to evaluate the similarity claim more rigorously. In the revised manuscript we will add a dedicated subsection (or table) in the observations section that directly compares the two outbursts using the available all-sky monitoring data. Quantities will include total outburst fluence, integrated mass-accretion rate (estimated from bolometric flux), outburst duration above a fixed flux threshold, peak flux, and the time since the previous outburst. These values will be derived from the same public light curves used in the original analysis, with references to the data sources. This addition will make the similarity assessment transparent and will not alter the central conclusion. revision: yes
-
Referee: [Modeling section] The thermal-evolution modeling (modeling section) that yields the new shallow-heating parameters is presented without tabulated best-fit values, uncertainties, or goodness-of-fit statistics for the 2016 data set; the abstract likewise omits these quantities. This prevents evaluation of whether the inferred change in heating depth and strength is statistically required by the data.
Authors: We acknowledge the omission of tabulated fit results. In the revised version we will include a new table in the modeling section that reports, for both the 2013 and 2016 data sets: the best-fit shallow-heating depth and strength, their 1-sigma uncertainties, the reduced chi-squared (or equivalent goodness-of-fit metric), and the number of degrees of freedom. We will also update the abstract to quote the key quantitative results (e.g., the factor by which depth and magnitude differ). These additions will allow direct assessment of whether the change in shallow heating is statistically required. revision: yes
Circularity Check
No significant circularity; inference relies on external simulations and new data
full rationale
The paper's central inference—that shallow heating depth and magnitude differed between the 2013 and 2016 outbursts—is obtained by fitting standard neutron-star thermal evolution models (from the literature) to new Chandra/Swift observations of post-2016 cooling. No equations or steps in the provided text reduce the inferred heating parameters to quantities already fitted from the same dataset, nor do they rely on self-citations, ansatzes smuggled via prior work, or renaming of known results. The assumption that the outbursts were sufficiently similar is an external modeling choice whose validity can be tested independently; it does not create a definitional or fitted-input loop within the derivation itself. The work is therefore self-contained against external benchmarks.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.