{"id":"bf2dcf44-8e5d-468d-b3d9-42a76cb6be4e","arxiv_id":"2505.07948","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Dark matter heat transport can erase convective cores in solar-mass stars, yielding asteroseismic constraints on dark matter-nucleon scattering and a 4 sigma hint of dark matter-electron scattering in KIC 8228742 that conflicts with direct detection limits.","lead":"This paper simulates how captured dark matter can transport heat inside stars and erase their convective cores, then uses asteroseismic data from one subgiant star to constrain dark matter interactions. It reports a 4 sigma preference for dark matter-electron scattering around 1 GeV, while noting this region is excluded by direct detection experiments.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (2.20) is internally inconsistent as printed: for K much less than K0 it gives L proportional to sigma^3, not the stated L proportional to sigma^-1, so the Knudsen peak in Fig. 2 cannot follow.","rationale":"The paper does something genuinely useful: it applies a Monte Carlo-calibrated transport treatment to asteroseismology, re-calibrates stellar nuisance parameters per DM model, and gives credible constraints for SD DM-nucleon scattering. The DM-nucleon exclusion contour in Fig. 10 is a reasonable result within the model assumptions. The problem is the headline DM-electron hint. That hint lives near the Knudsen peak of the transport luminosity, and the equation that defines the luminosity interpolation is not self-consistent as printed. In the LTE limit, Eq. (2.20) predicts L proportional to sigma0^3, while the text and Fig. 2 require L proportional to sigma0^-1 and a peak; the only way to reconcile these is to replace the printed factor with its reciprocal. If the code follows the printed equation, the peak and the best-fit island disappear; if the code follows the text, the paper needs a correction and, more importantly, validation of K0 for electrons. Either way, the quantitative claim of a greater than 4 sigma preference cannot be accepted on the strength of this manuscript. The statistical issues noted by the reader, including the single-star selection, the grid scan, and the lack of a look-elsewhere penalty, are secondary but reinforce the conditional verdict. I therefore keep the reader's CONDITIONAL verdict: the method is promising, but the headline hint is not demonstrably robust until Eq. (2.20) and the electron-target K0 are fixed and validated.","tokens_in":28315,"tokens_out":13220,"duration_ms":120229,"concrete_test":"Inspect the Capt'n General or MESA implementation, or rerun the Monte Carlo solver of Refs. [28,82], to see whether Lcalib is coded as 0.5 [1 + (K0/K)^2]^-1 LSP or as the printed 0.5 [1 + (K0/K)^2] LSP, and reproduce Fig. 2. If the code matches the printed form, the LTE branch will rise as sigma0^3 instead of falling as sigma0^-1 and the Knudsen peak will disappear. Separately, run the Monte Carlo transport for 1 GeV DM-electron in a 1.25 solar mass model over sigma0 = 10^-36 to 10^-33 cm^2, fit K0 for the electron target, and recompute the chi^2_r02 grid for KIC 8228742. If K0_e differs from 0.4 by more than about 20 percent, or if the interpolation direction is inverted, the (1 GeV, 3 x 10^-34 cm^2) best-fit island and its 4 sigma preference must be re-evaluated.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In Sec. 2.4, Eq. (2.20) gives Lcalib = 0.5 (1 + (K0/K)^2) LSP, with K = lchi/rchi proportional to sigma0^-1. In the LTE limit K much less than K0, the bracket is approximately (K0/K)^2, so Lcalib is approximately 0.5 (K0/K)^2 LSP, which scales as sigma0^3. The immediately following text states that in this limit the factor suppresses the luminosity, grows as K^2, and gives L proportional to sigma0^-1; Fig. 2 also shows a Knudsen peak. These statements describe the reciprocal of the printed expression, i.e. Lcalib = 0.5 [1 + (K0/K)^2]^-1 LSP, which for K much less than K0 behaves as 0.5 K^2/K0^2 LSP and therefore scales as sigma0^-1. As printed, Eq. (2.20) cannot produce a maximum in the L-sigma curve at all. This is not a cosmetic typo: the DM-electron preference at mchi = 1 GeV and sigma0 = 3 x 10^-34 cm^2 is driven by models near the Knudsen peak, and every derived constraint and the claimed 4 sigma improvement passes through this luminosity calibration. The paper also states that K0 approximately 0.4 was found for constant DM-nucleon interactions and depends on interaction and target type, yet no Monte Carlo validation for DM-electron is shown here. A reader cannot verify the headline signal from the manuscript as written.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the effects of captured asymmetric dark matter on stellar structure and evolution, using a Spergel-Press analytic heat transport formula with a Monte-Carlo-calibrated Knudsen correction developed in the authors' earlier work. It implements these effects in MESA via the Capt'n General code, first for idealized 1.1-1.5 solar-mass main-sequence stars and then for a self-consistently recalibrated model of the subgiant KIC 8228742. From the r02 frequency-separation ratios of this star, the paper derives constraints on spin-dependent DM-nucleon scattering and reports a ~4 sigma preference for DM-electron scattering with m_chi = 1 GeV and sigma_0 ~ 3e-34 cm^2, a region already in strong tension with direct-detection limits.","tokens_in":28579,"tokens_out":6290,"duration_ms":61584,"significance":"If the transport calibration is correct, the paper demonstrates a new asteroseismological probe of light dark matter, with a genuine methodological improvement in that stellar nuisance parameters are re-calibrated for every DM model point. The inclusion of DM-electron scattering in stellar heat transport is new, and the paper makes falsifiable predictions for the r02 ratios of stars near the convective-core boundary. These strengths are substantial. However, the central detection claim rests on an unverified calibration step and a statistically fragile improvement, so the significance depends on whether those issues can be resolved.","major_comments":[{"comment":"As printed, Eq. (2.20) is internally inconsistent with the text that follows it. For K << K0, Lcalib = 0.5 (1 + (K0/K)^2) LSP ~ 0.5 (K0/K)^2 LSP, and since K is proportional to sigma_0^{-1} while LSP is proportional to sigma_0, this gives Lcalib proportional to sigma_0^3. The text states that in this limit the correction 'suppresses the luminosity and grows as ~K^2, and so L ~ sigma_0^{-1}'; that behavior, and the Knudsen peak in Fig. 2, follows from Lcalib = 0.5 [1 + (K0/K)^2]^{-1} LSP instead. This is a load-bearing error: the claimed DM-electron preference at (m_chi, sigma_0) = (1 GeV, 3 x 10^{-34} cm^2) sits near the peak, and every derived constraint and the 4 sigma improvement pass through this luminosity calibration.","section":"Sec. 2.4, Eq. (2.20)"},{"comment":"The calibration constant K0 ~ 0.4 is imported from the authors' prior work on constant DM-nucleon interactions, yet the text explicitly says that K0 depends on interaction and target type. No Monte Carlo validation is shown for DM-electron interactions or for the subgiant model used for KIC 8228742, and Ref. [82] is an unreviewed preprint. Because the transport luminosity controls the size of the convective core and hence the r02 signal, an unvalidated K0 for electrons could shift the Knudsen peak and change the location and strength of the claimed preference; the manuscript provides no way for a reader to assess this.","section":"Sec. 2.4 and Refs. [28, 82]"},{"comment":"The reported >4 sigma preference is the difference between the no-DM chi^2_r02 of 29.16 and the best grid point's chi^2_r02 of 5.99, interpreted with 2 degrees of freedom. This significance is unreliable for two reasons. First, the best point was selected from a grid over (m_chi, sigma_0), so the look-elsewhere effect from scanning the grid is not accounted for. Second, Model B has chi^2_star = 0.396, which is worse than the no-DM value of 0.082, so the improvement is not a global improvement of the full fit; a proper model comparison must include the full likelihood and the extra parameters rather than the r02 diagnostic alone.","section":"Sec. 4.1.1, Table 1, Fig. 11"}],"minor_comments":[{"comment":"Equation (4.3) is garbled as printed (\"22X 15 =\") and should be written as a sum over the radial orders shown in Fig. 12, presumably n = 15 to 22.","section":"Sec. 4.1, Eq. (4.3)"},{"comment":"Both captions say \"1.5 M_sun (left)\" for the third panel, but the panels are left, middle, and right; the last should read \"(right)\".","section":"Fig. 4 and Fig. 5 captions"},{"comment":"The second integral in Eq. (A.25) is labeled \"I_vchi>vT\" but should be labeled \"I_vchi<=vT\" (or similar) to match the split defined in Eq. (A.23).","section":"Appendix A, Eq. (A.25)"},{"comment":"Equation (2.4) writes the differential scattering rate as dsigma/dv, but later expressions in the same section use dsigma/d(cos theta); the relationship between the two notations is not stated.","section":"Sec. 2.1, Eq. (2.4)"}],"recommendation":"major_revision","confidential_remarks":"The paper advertises a >4 sigma hint in the abstract, which will attract attention. The editor should ask the authors to resolve the sign inconsistency in Eq. (2.20), to provide the Monte Carlo calibration for DM-electron interactions (or clearly restrict the claims to calibrated regimes), and to correct the statistical interpretation of the 4 sigma preference with a trials factor and a global likelihood. If the sign error is confined to the manuscript text and the code implements the reciprocal form, the central claims may survive, but the current version is not self-contained enough for the results to be verified."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: the paper is genuinely useful for the method—convective-core erasure as a probe of DM-electron scattering, with proper recalibration of stellar parameters per grid point—but you cannot trust the headline 4-sigma preference for mχ=1 GeV electrons. Eq (2.20) as printed is internally inconsistent: for K≪K0 it gives L∝σ^3, not σ^{-1}, so the Knudsen peak in Fig 2 cannot follow. The text describes the reciprocal of the printed expression. Since every constraint and the claimed signal passes through this luminosity calibration, the central result is unsupported as written.\n\nWhat is new and good: The finite-temperature capture treatment for electrons is a clear improvement; thermal motion genuinely matters and they show it. The per-grid-point recalibration of nuisance parameters is something prior work skipped, and it makes the DM-nucleon constraints more credible. The analytic v4/q4 capture and transport expressions in Appendix A, including a correction to earlier literature, look like a real service to the field. The MESA/Capt'n General setup is standard, reproducible, and the pipelines are described well enough to re-run.\n\nSoft spots, in order: (1) the Eq (2.20) reciprocal error above; (2) K0≈0.4 is imported from the same group's Monte Carlo calibration for constant DM-nucleon interactions, with no validation shown here for DM-electron scattering or for the specific mass/evolution stage of KIC 8228742; (3) the 4-sigma preference comes from a single star, has no look-elsewhere penalty, and the paper itself notes XENON1T's S2-only search excludes the preferred region; (4) the observational uncertainties on the calibration are large and the paper admits this. None of these invalidate the method, but together they mean the hint needs much more work before it should be quoted.\n\nWho it is for: anyone working on DM capture in stars or asteroseismic limits on sub-GeV DM. It deserves a serious referee. I'd send it to peer review with a clear request to fix Eq (2.20), validate or re-derive the electron transport calibration, and add a multi-star or look-elsewhere discussion before the hint is claimed.","headline":"A promising new asteroseismic probe of sub-GeV DM, but the headline 4-sigma hint rests on an inconsistent luminosity formula that needs fixing before the claim can be taken at face value.","tokens_in":29194,"tokens_out":2820,"would_cite":true,"duration_ms":26626,"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":"Oscillation data for the nearby subgiant KIC 8228742 favour dark-matter-electron scattering with a 1 GeV particle at more than 4 sigma, while dark-matter heat transport can erase the star's convective core.","keywords":["dark matter","asteroseismology","asymmetric dark matter","convective core","dark matter-electron scattering","spin-dependent scattering","stellar heat transport","subgiant star"],"falsifier":"Run a direct Monte Carlo simulation of heat transport by $\\sim1\\,\\mathrm{GeV}$ dark matter scattering on electrons in a KIC 8228742-like stellar model at $\\sigma_0\\simeq3\\times10^{-34}\\,\\mathrm{cm^2}$: if the resulting temperature and sound-speed profiles do not erase the convective core, or if the predicted $r_{02}$ ratios do not match the observed trend, the reported $4\\sigma$ preference disappears. A second nearby subgiant of similar mass whose $r_{02}$ ratios agree with the standard model would likewise disfavour the dark-matter-electron interpretation.","tokens_in":28015,"feed_emoji":"🔭","tokens_out":16071,"duration_ms":139674,"temperature":0.7,"pith_summary":"Stars just above the convective-core threshold are sensitive probes of asymmetric dark matter: if captured dark matter does not annihilate, it conducts heat out of the core, cools it, and can erase the convective core. The paper simulates this using Monte-Carlo-calibrated heat transport and compares the predicted oscillation-frequency ratios to measured ratios for the nearby subgiant KIC 8228742. It reports new constraints on spin-dependent dark-matter-nucleon scattering and a more than $4\\sigma$ preference for dark-matter-electron interactions at $m_\\chi\\simeq1\\,\\mathrm{GeV}$ and $\\sigma_0\\simeq3\\times10^{-34}\\,\\mathrm{cm^2}$, visible as an improvement in the $r_{02}$ frequency-ratio fit from $\\chi^2_{r02}\\simeq29.16$ to $5.99$. The preference sits in strong tension with Earth-based direct-detection limits, but the paper argues it is exactly the kind of signature asteroseismology is built to find.","feed_headline":"Pulsations of a nearby star hint at 4σ dark-matter-electron scattering","feed_subtitle":"Captured dark matter would cool the core of KIC 8228742 and erase its convective core.","key_machinery":"The load-bearing object is the calibrated heat-transport luminosity: the analytic isothermal, long-mean-free-path transport expression is rescaled by a Knudsen-number correction, $$L_{\\mathrm{calib}}=0.5\\left(1+(K_0/K)^2\\right)L_{\\mathrm{SP}},$$ with $K=\\ell_\\chi(0)/r_\\chi$ and $K_0\\simeq0.4$. This fixes how much energy the captured dark matter removes from the core and deposits in the outer layers at each radius, and therefore whether the convective core survives. The comparison observable is the frequency-separation ratio $r_{02}(n)=d_{02}(n)/\\Delta_1(n)$, a combination of small and large frequency separations that is insensitive to surface effects and sensitive to the sharp sound-speed discontinuity at a convective-core boundary. The analysis also recalibrates the stellar nuisance parameters separately for every dark-matter mass and cross-section point, rather than reusing the no-dark-matter stellar model.","core_discovery":"The central claim is that a captured population of asymmetric, non-annihilating dark matter can erase the convective core of a star near the $\\sim1.25\\,M_\\odot$ threshold, and that this structural change shows up in asteroseismic frequency-separation ratios. Using the subgiant KIC 8228742, the paper obtains constraints on spin-dependent dark-matter-nucleon scattering and finds a best-fit dark-matter-electron model with $m_\\chi=1\\,\\mathrm{GeV}$ and $\\sigma_0\\simeq3\\times10^{-34}\\,\\mathrm{cm^2}$ that is more than $4\\sigma$ better than the no-dark-matter model, improving $\\chi^2_{r02}$ from $29.16$ to $5.99$. The paper reads this as evidence for an extra heat-transporting component inside the star, while noting that the required cross section greatly exceeds current Earth-based direct-detection bounds, so a dark-matter interpretation would need an additional ingredient such as a local density enhancement or a velocity-dependent interaction.","pith_inferences":["The sharpest test would be to compute $K_0$ for electron scattering from first principles: a direct Monte Carlo transport simulation inside a KIC 8228742-like model would show whether the $4\\sigma$ best fit survives when the electron calibration is not borrowed from the nucleon case.","A sample of several $\\sim1.25\\,M_\\odot$ subgiants with measured $r_{02}$ ratios could serve as a population-level check; if only one star shows the convective-core-erasure signature, a dark-matter explanation would need something local, such as a dark-matter overdensity, rather than a universal particle property.","The best-fit electron cross section lies orders of magnitude above terrestrial electron-recoil limits, so a viable particle model would likely require a velocity- or momentum-dependent form factor that suppresses Earth-based rates while remaining efficient inside the star; the paper mentions this route but does not compute it for this star.","Because the transport calibration was validated for constant cross sections, the appendix's new analytic expressions for $v^{2n}$ and $q^{2n}$ interactions invite a re-fit of the KIC 8228742 data for non-constant cross sections, where the preferred mass and cross section could shift by factors of order unity."],"forward_implications":["Convective-core erasure becomes a generic, observable signature of asymmetric dark matter in stars just above $\\sim1.2\\,M_\\odot$, and dark-matter-electron scattering can produce it even where nucleon scattering is kinematically suppressed.","The single-star data from KIC 8228742 tighten spin-dependent dark-matter-nucleon constraints below $m_\\chi\\sim3\\,\\mathrm{GeV}$, reaching cross sections near $10^{-37}\\,\\mathrm{cm^2}$, in a low-mass region where evaporation may deplete the captured population.","If the preference for dark-matter-electron scattering is real, the star's core is cooler and denser than the standard model predicts, and the $r_{02}$ ratios should follow the near-linear trend the best-fit model reproduces.","Velocity- and momentum-dependent cross sections give qualitatively similar core erasure, so the reported effects are not tied to the constant-cross-section choice.","With more precise observations and additional subgiants of similar mass, the $4\\sigma$ preference should either be confirmed as a population-wide pattern or disappear."],"supporting_citations":[{"why":"Gives the analytic isothermal heat-transport expression that the paper rescales with the Knudsen correction.","marker":"[61]"},{"why":"Provides the capture and evaporation formalism used to set the dark-matter population and mark the low-mass evaporation region.","marker":"[42]"},{"why":"Supplies the finite-temperature capture treatment for electrons and the evaporation-mass estimates used for the low-mass region.","marker":"[25]"},{"why":"Supplies the Monte Carlo simulations and the calibrated heat-transport rescaling that the paper adopts.","marker":"[28]"},{"why":"Extends the Monte Carlo calibration to realistic stellar and planetary models and supports the K0 value used here.","marker":"[82]"},{"why":"The earlier asteroseismic study of KIC 8228742 whose no-dark-matter model and r02 fit the paper recalibrates and extends to electron scattering.","marker":"[78]"},{"why":"Source of the observed oscillation frequencies used to construct the r02 frequency-separation ratios.","marker":"[107]"},{"why":"Source of the spectroscopic and asteroseismic constraints used in the stellar calibration.","marker":"[106]"},{"why":"The direct-detection limit that is in strong tension with the paper's dark-matter-electron best-fit region.","marker":"[112]"}],"fun_headline_variants":["Asteroseismology hints at 4σ dark matter-electron scattering","Star's missing core hints at dark matter: 4σ, but tension","Dark matter may erase a star's convective core, pulsations say","Subgiant pulsations point to dark matter electron hits (4σ)"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"Everything rests on the calibration factor that converts the analytic dark-matter heat-transport luminosity into the actual transported luminosity; the paper takes the value fitted for constant dark-matter-nucleon scattering and applies it to electrons without re-deriving it for the electron case.","fun_headline_variants_meta":{"raw":{"variants":["Asteroseismology hints at 4σ dark matter-electron scattering","Star's missing core hints at dark matter: 4σ, but tension","Dark matter may erase a star's convective core, pulsations say","Subgiant pulsations point to dark matter electron hits (4σ)"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000236,"raw_usage":{"total_tokens":1518,"prompt_tokens":971,"completion_tokens":547,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":587,"completion_tokens_details":{"reasoning_tokens":468}},"tokens_in":587,"tokens_out":547,"duration_ms":5332,"temperature":1.0,"reasoning_tokens":468,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T22:06:05.167803+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a direct Monte Carlo simulation of heat transport by $\\sim1\\,\\mathrm{GeV}$ dark matter scattering on electrons in a KIC 8228742-like stellar model at $\\sigma_0\\simeq3\\times10^{-34}\\,\\mathrm{cm^2}$: if the resulting temperature and sound-speed profiles do not erase the convective core, or if the predicted $r_{02}$ ratios do not match the observed trend, the reported $4\\sigma$ preference disappears. A second nearby subgiant of similar mass whose $r_{02}$ ratios agree with the standard model would likewise disfavour the dark-matter-electron interpretation.","supporting_citations":[{"cited_title":"On asymmetric dark matter constraints from the asteroseismology of a subgiant star","cited_arxiv_id":"2109.12671","evidence_quote":"The earlier asteroseismic study of KIC 8228742 whose no-dark-matter model and r02 fit the paper recalibrates and extends to electron scattering."},{"cited_title":"Oscillation mode frequencies of 61 main sequence and subgiant stars observed by Kepler","cited_arxiv_id":"1204.3147","evidence_quote":"Source of the observed oscillation frequencies used to construct the r02 frequency-separation ratios."}],"review_version":1}