{"id":"e7ea1ba2-3da4-49b3-812b-dcbf414f74f1","arxiv_id":"2412.14342","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Monte Carlo simulations show the corrected Spergel and Press formalism reproduces dark matter heat transport in realistic Sun, brown dwarf, and Earth models, but with per-scenario fitted parameters.","lead":"This paper uses Monte Carlo simulations to test how captured dark matter conducts heat inside the Sun, a brown dwarf, and the Earth, using realistic density and temperature profiles. It finds that a corrected Spergel and Press approximation describes the simulated heat flow, and that previous dark matter evaporation rates for the Sun appear robust.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Single-trajectory steady-state convergence is asserted, not demonstrated; admitted finite-time bias can be absorbed into fitted A and K0, so Eq. (7)'s claimed accuracy is not yet established.","rationale":"I agree with the reader's weakest assumption; it is the most load-bearing prerequisite for the central claim. The reader also noted the in-sample nature of the fit as a separate weakness, but the more fundamental issue is whether the Monte Carlo ground truth itself is unbiased. The paper's own caveats make this concern live: finite-time suppression of evaporation is admitted, and convergence could not be obtained at higher masses. If the fitted parameters vary with run length, the central validation collapses because the agreement in Fig. 3 would partly reflect a systematic bias absorbed into A and K0. The proposed convergence test directly probes whether the finite-time bias is small compared to the quoted statistical errors at representative calibration points. If the parameters are stable with run length, the concern is retired; if not, the paper must either restrict its claim to converged regimes or provide a bias-correction procedure. This sharpens the reader's conditional verdict but does not change it, since the paper currently lacks the convergence evidence needed to fully support the abstract's claim.","tokens_in":15637,"tokens_out":9575,"duration_ms":90803,"concrete_test":"Run cosmion at a representative LTE-regime point (e.g., Sun SD, mχ = 3 GeV, σ0 = 10^-34 cm^2) for N, 2N, and 4N collisions, and re-fit A and K0 to each dataset. If A or K0 shifts by more than the quoted statistical error between 2N and 4N, the recorded heat transport is not steady-state and the Eq. (7) calibration is unreliable.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that Eq. (7), with A and K0 fitted to simulation, accurately describes DM heat transport in realistic potentials. The load-bearing condition is that the single-particle Monte Carlo has reached steady state: the text assumes 'By ergodicity, the particle traces out the full steady-state phase space distribution' (Sec. II), but no convergence diagnostics are shown. The paper itself reports in Sec. IV A that 'converged numerical values could not be obtained' for higher DM masses and in Sec. IV B that finite Monte Carlo time 'artificially suppress[es] the recorded rate.' Since the fitted A and K0 in Table I are derived from these same runs, any residual finite-time bias is absorbed into the fitted parameters and into the smooth curves of Fig. 3. The agreement between Eq. (7) and the simulations is therefore not independent evidence for the formalism's accuracy until steady-state convergence is established at the calibration points. In particular, the LTE side of the Knudsen transition, where the paper notes exponentially more collisions are needed, is exactly the region sampled when determining K0.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper extends the Monte Carlo framework of Banks et al. (2022) to simulate heat transport by asymmetric dark matter in realistic gravitational potentials and density/temperature/composition profiles of the Sun, a 0.01 solar-mass brown dwarf, and the Earth. The authors include both spin-dependent scattering on hydrogen and spin-independent scattering on multiple nuclear species, and they compare the simulated luminosity profiles and maximum luminosities with the corrected Spergel and Press formula, Eq. (7), using values of A and K0 fitted to the simulations and reported in Table I. They also compute solar dark-matter evaporation rates and compare them with analytic estimates from the literature, concluding that previous evaporation calculations appear robust. The Cosmion code is made publicly available.","tokens_in":15851,"tokens_out":4182,"duration_ms":39737,"significance":"If the central claim holds, this work is a valuable step for the dark-matter-in-stars community: it provides an open-source Monte Carlo tool, extends validation of heat-transport formalisms to realistic potentials and to spin-independent multi-isotope scattering, and gives concrete values for the effective parameters A and K0. The comparison in Fig. 2 of full radial luminosity profiles not directly used in the fitting provides some independent shape evidence for the corrected Spergel and Press form. However, because the maximum-luminosity agreement in Fig. 3 is obtained by fitting both parameters to the same simulation points, and because the paper explicitly reports convergence difficulties in parts of the parameter space, the headline claim that Eq. (7) 'remains accurate across all celestial objects considered' is not yet established as strongly as the text suggests.","major_comments":[{"comment":"The central accuracy claim rests on a two-parameter fit of A and K0 to the same Monte Carlo maximum luminosities shown in Fig. 3. The text states in Sec. IV A that A is left as a free parameter 'that we fit based on simulation results', and Table I reports these fitted values, so the agreement between the solid curves and the data points in Fig. 3 is in-sample rather than an independent test of Eq. (7). The radial profiles in Fig. 2 are not directly used in the fitting and provide some independent shape support, but the conclusion in Sec. V that Eq. (7) 'provides a good parametrization in all the regimes that we have tested' requires either an explicit out-of-sample test (for example, reserving a subset of mass/cross-section grid points for validation) or a goodness-of-fit statistic that accounts for the two fitted parameters.","section":"Sec. IV A, Eq. (7), Table I, Fig. 3"},{"comment":"The steady-state assumption is load-bearing but not demonstrated. Section II asserts that a single trajectory traces out the full steady-state phase-space distribution by ergodicity, yet Sec. IV A reports that converged numerical values could not be obtained at higher DM masses and that the LTE regime requires exponentially more collisions. Since Table I and Fig. 3 are calibrated on runs with 10^6 to 10^9 collisions, any residual finite-time bias can be absorbed into A and K0. I ask for explicit convergence diagnostics at representative calibration points, e.g., the maximum luminosity as a function of the number of collisions or a comparison of independent random seeds, and for a statement of the resulting systematic uncertainty on A and K0.","section":"Sec. II and Sec. IV A"},{"comment":"The evaporation comparison is weakened by the acknowledged finite-time suppression. The text says that particles caught in long orbits 'will artificially extend the time recorded during which no evaporation has taken place, and thus artificially suppress the recorded rate', and that the quoted error bars only reflect sqrt(N_evap)/t_sim. The concluding sentence in Sec. V that previous evaporation rates are 'consistent with our simulation results' is therefore stronger than the evidence shown. Please quantify the suppression, for example with longer runs at one or two representative masses or with an analytic bound on the one-sided bias, or soften the conclusion to a consistency check that explicitly accounts for this systematic.","section":"Sec. IV B and Fig. 5"},{"comment":"The claim in Sec. IV A that K0 is 'robustly found to be between 0.4 and 0.5' is not supported by the table values for the Sun, where K0 ranges from 0.271(12) at 20 GeV spin-dependent to 0.480(2) at 3 GeV spin-dependent, and the 3 GeV spin-independent entry is 0.271(2). Because K0 controls the position of the Knudsen transition and is used for the interpolated maps in Fig. 4, the mass and interaction-type dependence should either be discussed and explained or the 'robustly' wording should be revised.","section":"Table I"}],"minor_comments":[{"comment":"There are typos in the abstract ('asteroseismoloigcal') and in the Introduction ('upmost'), which should read 'asteroseismological' and 'utmost', respectively.","section":"Abstract and Introduction"},{"comment":"The Introduction lists a red giant star among the astrophysical bodies to be studied, but no red-giant results appear in the paper; either add such a simulation or remove the mention.","section":"Section I"},{"comment":"The caption refers to shaded regions from the Monte Carlo simulations, but the legend does not explicitly label the shaded band; please clarify the legend so the reader knows which entry corresponds to the shaded region.","section":"Fig. 1 caption"},{"comment":"The footnote about the missing square root in Algorithm 1 of Ref. [11] is useful, but it should specify which equation or algorithm line in that reference is being corrected so that readers can verify the implementation.","section":"Footnote 1"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a solid technical contribution with a useful public code, but the central validation is in-sample and the convergence of the simulations is asserted rather than demonstrated. I would encourage the editor to invite a revision in which the authors perform explicit out-of-sample checks and convergence diagnostics; with those additions, the paper could make a strong case for its main claim."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The real news here is that the authors took their earlier SHO-potential Monte Carlo code, generalized it to arbitrary gravitational potentials and to spin-independent scattering off multiple nuclear species, and are shipping the code. That is genuinely useful. The paper is also honest: it explicitly says converged runs could not be obtained for high DM masses and that finite simulation time artificially suppresses evaporation rates. Those admissions are in the text, not hidden.\n\nWhat the paper does well is show that the functional form L(r) = A / (1 + (K0/K)^2) L_SP(r), with A and K0 fitted per scenario, tracks the simulated luminosity curves in the Sun, a brown dwarf, and the Earth. The radial profile shape is not directly fitted—only the max luminosity is used to set A and K0—so the agreement in Fig. 2 is partly independent support. That matters.\n\nBut the stress-test concern lands. The two parameters are fitted to the same maximum luminosities that are then compared in Fig. 3, so the agreement is in-sample. The paper does not show convergence diagnostics for the single-particle phase-space distribution, and the LTE side of the Knudsen transition is exactly where the simulation struggles. Since K0 is determined near that transition, a finite-time bias could be absorbed into the fitted parameters. The evaporation section is explicitly limited: the error bars ignore the known suppression, and the concluding sentence that literature evaporation rates are \"trustworthy\" goes beyond what this method can currently demonstrate. The authors themselves note the LTE-regime convergence was not achieved.\n\nSo the central claim—that Eq. (7) is accurate across all objects and regimes—is not yet established. It is plausible, and the code gives others a way to check. What is missing is an out-of-sample test: fix A and K0 from one mass (or one object) and predict another, or show convergence of the phase-space distribution before fitting. That is a reasonable revision, not a fatal flaw.\n\nWho is this for? Anyone doing dark-matter heat transport in stars or planets. The paper deserves a serious referee: the software is public, the method is reproducible, and the question of which heat-transport formalism to trust is consequential for asteroseismology and fusion-rate predictions. Send it out, but ask for out-of-sample validation and a clearer statement of what is predictive versus what is interpolation.","headline":"Useful extension of the authors' own MC framework to realistic potentials and multi-species SI scattering, but the central validation is weakened by fitting A and K0 to the same data; deserves a serious referee with demands for out-of-sample tests.","tokens_in":16361,"tokens_out":957,"would_cite":true,"duration_ms":11770,"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":"This paper shows that a two-parameter corrected Spergel-Press formula describes dark-matter heat transport in the Sun, a brown dwarf, and the Earth.","keywords":["dark matter","heat transport","Spergel-Press formalism","Monte Carlo simulation","Knudsen regime","evaporation","spin-independent scattering","cosmion"],"falsifier":"Run the public cosmion code on a red-giant model at the same dark matter masses and cross sections used for the Sun, and compare the fitted A and K0 in Eq. (7) with Table I; if either parameter moves outside the quoted uncertainties, the claimed universality of the calibrated form is falsified. A cheaper check within the present data is to verify convergence of the 20 GeV solar spin-independent point, since the paper states that converged values could not be obtained at higher masses.","tokens_in":1956,"feed_emoji":"☀️","tokens_out":8128,"duration_ms":106073,"temperature":0.7,"pith_summary":"This paper claims that a two-parameter corrected version of the Spergel-Press formula describes the heat transported by dark matter in three very different astrophysical bodies: the Sun, a brown dwarf, and the Earth. The correction rescales the Spergel-Press luminosity by A divided by (1 plus (K0/K)^2), with A near 0.5 and K0 around 0.4, and it holds for both spin-dependent scattering on hydrogen and spin-independent scattering on multiple nuclear species. The paper validates this form against Monte Carlo simulations that track a dark matter particle through millions of collisions in each body's realistic gravitational potential, temperature, and density. It also reports that simulated evaporation rates from the Sun match earlier analytic estimates. The public code cosmion is released with the paper.","feed_headline":"Dark matter heat transport confirmed in Sun, dwarf, Earth","feed_subtitle":"Realistic-potential runs back the corrected Spergel-Press form in both spin channels for all three bodies.","key_machinery":"The central object is Eq. (7): L(r) = A/(1+(K0/K)^2) * L_SP(r), a two-parameter interpolation between the isothermal (Spergel-Press, K >> 1) and local-thermal-equilibrium regimes. Here L_SP is the Spergel-Press luminosity from Eqs. (2) and (5), K = l_chi(r=0)/r_chi is the Knudsen number, A tunes the overall normalization (the Spergel-Press form overpredicts by about a factor of two, so A approx 1/2), and K0 sets the cross section where transport switches regimes. The supporting machinery is the Monte Carlo random walk itself: trajectories integrated in the real potential $\\varphi$(r) with an RKF45 solver using optical depth as the dependent variable, collision rates summed over nuclear isotopes, and a species-selection truncation that keeps spin-independent simulations tractable.","core_discovery":"The central claim is that the calibrated Spergel-Press formalism, Eq. (7), reproduces the dark-matter heat-transport luminosity obtained from direct Monte Carlo integration of the Boltzmann collision equation in realistic gravitational potentials of the Sun, a 0.01 solar-mass brown dwarf, and the Earth. The two fitted parameters—the prefactor A, around 0.43 for spin-dependent scattering and up to 0.53 for spin-independent, and the Knudsen transition location K0, between roughly 0.27 and 0.48—capture departures from the idealized isothermal and local-thermal-equilibrium limits. The fit targets the maximum luminosity across roughly five orders of magnitude in cross section, covering the transition from the long-mean-free-path Knudsen regime to the local-thermal-equilibrium regime, for dark matter masses from 1 to 200 GeV depending on the body. The paper argues that the commonly used Gould-Raffelt LTE approach mismodels the shape of the heat transport profile near that transition, and that the corrected Spergel-Press form is the more reliable parameterization.","pith_inferences":["Inference: because A and K0 show only weak mass and species dependence, Eq. (7) may serve as a universal interpolant for any spherically symmetric body; a testable extension is to run cosmion on a red giant or white dwarf and check whether the fitted parameters stay in the same ranges.","Inference: the paper's finite-time suppression of recorded evaporation rates implies that the true evaporation floor for low-mass dark matter in the Sun could be higher than the simulation records; rare-event resampling or longer runs would give a sharper comparison to analytic evaporation rates.","Inference: the explicit Keplerian treatment of bound orbits outside the star is the time-reverse of halo capture, so the same code could be extended to compute capture rates from unbound orbits and verify them against standard capture formalisms, which the paper leaves to future work.","Inference: the claimed accuracy applies only to mass ranges where converged simulations were obtained; for the Sun these stop near 20 GeV, so the universality statement should not be read as covering the higher-mass regime where the paper reports non-convergence."],"forward_implications":["The two-parameter form of Eq. (7) can replace the Gould-Raffelt LTE treatment in stellar and planetary modeling, because the paper finds that the LTE form mismodels the heat transport shape near the Knudsen transition.","Predictions of fusion-rate modifications and asteroseismological signatures in the Sun can be computed with the corrected Spergel-Press form using the fitted A and K0 values from Table I.","Spin-independent scattering with multiple nuclear species is validated in the same framework, so heat transport in hydrogen-poor bodies such as the Earth is covered by the same calibrated form.","Evaporation rates from the Sun computed with the interpolation of Ref. [29] are consistent with the simulations, supporting the use of those rates in capture-evaporation equilibrium calculations.","The public cosmion code allows the same test to be run on other astrophysical bodies, extending the verified domain beyond the three objects considered."],"supporting_citations":[{"why":"Supplies the corrected Spergel-Press parameterization and the original Monte Carlo method that this work extends from a harmonic-oscillator potential to realistic potentials.","marker":"[1]"},{"why":"Defines the original isothermal (long-mean-free-path) heat transport formalism whose corrected version is being tested.","marker":"[3]"},{"why":"Defines the local-thermal-equilibrium formalism that the paper compares against and provides the simulation-based finding that Spergel-Press overpredicts luminosities by about a factor of two.","marker":"[4]"},{"why":"Provides the original Monte Carlo simulations in a simple harmonic oscillator potential that serve as the baseline for the present extension to realistic potentials.","marker":"[8]"},{"why":"Proposes the interpolation between isothermal and LTE radial distributions that the paper includes as a comparison in Fig. 1.","marker":"[9]"},{"why":"Provides the B16 AGSS09met standard solar model used for the Sun's realistic temperature, density, and composition profiles.","marker":"[12]"},{"why":"Supplies the compiled Earth temperature, density, and composition profiles used for the planet simulations.","marker":"[23]"},{"why":"Supports the factor-of-two overprediction of Spergel-Press transport that motivates the A approx 1/2 normalization.","marker":"[27]"},{"why":"Provides the analytic evaporation rates and the isothermal, LTE, and interpolated treatments that the paper's evaporation simulations are compared against.","marker":"[29]"}],"fun_headline_variants":["Monte Carlo shows dark matter heat transport works in real bodies","Dark matter heat transport formula validated in Sun, Earth, brown dwarf","Simulations back dark matter heat transport in realistic stars and planets","Asymmetric dark matter heat flow confirmed in three celestial bodies","New code verifies dark matter heat transport in Sun, dwarf, Earth"],"cache_read_input_tokens":18560,"weakest_assumption_plain":"The single-particle Monte Carlo tracks one dark matter particle through $10^{6}$ to $10^{9}$ collisions and relies on ergodicity to represent the steady-state ensemble; finite simulation time is known to suppress recorded evaporation rates and prevented converged results at higher dark matter masses in the paper.","fun_headline_variants_meta":{"raw":{"variants":["Monte Carlo shows dark matter heat transport works in real bodies","Dark matter heat transport formula validated in Sun, Earth, brown dwarf","Simulations back dark matter heat transport in realistic stars and planets","Asymmetric dark matter heat flow confirmed in three celestial bodies","New code verifies dark matter heat transport in Sun, dwarf, Earth"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000196,"raw_usage":{"total_tokens":1336,"prompt_tokens":896,"completion_tokens":440,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":512,"completion_tokens_details":{"reasoning_tokens":352}},"tokens_in":512,"tokens_out":440,"duration_ms":4258,"temperature":1.0,"reasoning_tokens":352,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T12:19:01.447345+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the public cosmion code on a red-giant model at the same dark matter masses and cross sections used for the Sun, and compare the fitted A and K0 in Eq. (7) with Table I; if either parameter moves outside the quoted uncertainties, the claimed universality of the calibrated form is falsified. A cheaper check within the present data is to verify convergence of the 20 GeV solar spin-independent point, since the paper states that converged values could not be obtained at higher masses.","supporting_citations":[{"cited_title":"The DM temperature Tχ in this case is the average of the nucleon temperatures with which the DM is in thermal contact with, weighted by the local inter- action rate","cited_arxiv_id":null,"evidence_quote":"Supplies the corrected Spergel-Press parameterization and the original Monte Carlo method that this work extends from a harmonic-oscillator potential to realistic potentials."},{"cited_title":"[4, 8] ( mχ = 1 kg, R⋆ = 2.4 m), and a realistic solar model using temperature, density and composition data from the AGS05 Standard Solar Model","cited_arxiv_id":null,"evidence_quote":"Defines the original isothermal (long-mean-free-path) heat transport formalism whose corrected version is being tested."},{"cited_title":"Po- sitions are drawn randomly from the radial distri- bution in Eq","cited_arxiv_id":null,"evidence_quote":"Defines the local-thermal-equilibrium formalism that the paper compares against and provides the simulation-based finding that Spergel-Press overpredicts luminosities by about a factor of two."},{"cited_title":"(19) To do this we implement the sampling algorithm outlined in Ref","cited_arxiv_id":null,"evidence_quote":"Provides the original Monte Carlo simulations in a simple harmonic oscillator potential that serve as the baseline for the present extension to realistic potentials."},{"cited_title":"1 Noting a missing square root in the conditional expression of Algorithm 1: ξ6 < p y2 + z − 2xyµ/(y + x)","cited_arxiv_id":null,"evidence_quote":"Proposes the interpolation between isothermal and LTE radial distributions that the paper includes as a comparison in Fig. 1."},{"cited_title":"Effect of hypo- thetical, weakly interacting, massive particles on energy transport in the solar interior,","cited_arxiv_id":null,"evidence_quote":"Provides the B16 AGSS09met standard solar model used for the Sun's realistic temperature, density, and composition profiles."},{"cited_title":"7.06 - temperatures, heat, and energy in the mantle of the Earth,","cited_arxiv_id":null,"evidence_quote":"Provides the analytic evaporation rates and the isothermal, LTE, and interpolated treatments that the paper's evaporation simulations are compared against."}],"review_version":1}