{"id":"102473ec-4fc9-45a3-be00-920d57421720","arxiv_id":"2412.07953","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Updated black-hole formation bounds show Population III stars can constrain bosonic asymmetric dark matter at lower masses than neutron stars when a Bose-Einstein condensate forms.","lead":"The authors calculate when captured asymmetric dark matter would collapse neutron stars or Population III stars into black holes, and use the stars' survival to rule out dark matter masses and cross sections. The new result is that the first stars, observed by JWST, could extend these constraints to very light bosonic dark matter if it forms a Bose-Einstein condensate.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Thermalization is the load-bearing step for the new Pop III exclusion: Eq. (3.5) is a degenerate-NS formula, and the paper's own caveat lifts the bounds if the 10^6 yr Pop III stars have not thermalized their captured DM.","rationale":"I agree with the reader that thermalization is the weakest link. The reader's conditional verdict and my concern coincide: the Pop III low-mass exclusion is the paper's headline contribution, and it depends on thermalization within 10^6 yr. The manuscript itself contains the caveat that bounds must be lifted if thermalization fails, and §5 notes that much parameter space is excluded by the thermalization requirement, yet the paper does not provide a dedicated non-degenerate calculation for Pop III stars. This is an internal consistency risk, not a disagreement with consensus: the same machinery may be fine for NSs, but the new Pop III reach is not independently verified. The abstract's 10^-15 GeV claim adds a further BEC assumption that also requires a thermalized population. Because the reader already conditioned acceptance on this and related gaps, my stress test does not move the verdict; if the check passes, the conditional acceptance can be upgraded, and if it fails, the central claim should be rejected.","tokens_in":17282,"tokens_out":9534,"duration_ms":101429,"concrete_test":"Recompute the Pop III thermalization time with a non-degenerate target formalism appropriate to a hydrogen polytrope at T_c=2×10^7 K (e.g., Gould 1987; Garani & Palomares-Ruiz 2022), including relativistic kinematics if mχ≲T_c, and evaluate it at (mχ=10^-8 GeV, σ=10^-40 cm^2) and at the abstract's BEC endpoint. Compare with the 10^6 yr star age and overlay the resulting no-thermalization region on Fig. 5; if either advertised endpoint falls in that region, the central Pop III claim is not supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The new Pop III constraints—the advertised reach below mχ≈10^-8 GeV, and the abstract's 10^-15 GeV BEC line—are only valid if captured DM thermalizes within the 10^6 yr stellar lifetime. The paper says so itself in §3: 'we assume that the DM particles have thermalized. If this is not the case, we can no longer guarantee that a BH would form, so we must lift our bounds.' The only thermalization timescale quoted, Eq. (3.5), is the degenerate-neutron-star expression from [19], with p_F≈0.575 GeV in the numerator. The Pop III star is a non-degenerate hydrogen polytrope with T_c=2×10^7 K, so the target momentum is ∼(3 m_n T_c)^(1/2)≈2×10^-3 GeV and the scattering kinematics are not those of a Fermi sea; applying Eq. (3.5) as-is is not justified. The claimed endpoint (mχ=10^-8 GeV, σ=10^-40 cm^2) is exactly the regime where the thermalization requirement is most restrictive, and no independent non-degenerate (or relativistic, since mχ≪T_c) thermalization calculation is shown. If t_th exceeds 10^6 yr there, the headline Pop III exclusion must be lifted, and the BEC-based 10^-15 GeV reach fails as well, because BEC formation presupposes a thermalized DM population. This is a load-bearing assumption that the manuscript flags but does not resolve.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies constraints on asymmetric dark matter (ADM) via black hole formation in neutron stars (NSs) and Population III (Pop III) stars. For NSs, the authors extend the earlier analysis of McDermott, Yu, and Zurek by including multiscatter capture, a nucleon form-factor suppression, and a complete treatment of evaporation, including evaporation from a Bose-Einstein condensate (BEC). For Pop III stars, they apply the same formalism to a 1000 solar-mass, 10^6-year-old star in a high-density DM environment, claiming that observation of such a star can probe bosonic ADM down to m_chi ~ 10^-8 GeV at sigma ~ 10^-40 cm^2, and the abstract states a reach below 10^-15 GeV if a BEC forms. They also present single-integral approximations for the evaporation rate from arbitrary polytropic objects (Appendix A) and from a BEC state (Appendix B), and check the uniform-density radius approximation to within 1%.","tokens_in":17654,"tokens_out":8271,"duration_ms":74628,"significance":"If the derived bounds are correct, this work would extend ADM constraints to masses many orders of magnitude below direct-detection limits, using a class of objects (Pop III stars) that may be observed by JWST. The paper contains analytic derivations that are mostly self-consistent, and it explicitly checks the uniform-density approximation used for the DM radius. The evaporation approximations in the appendices, especially the BEC evaporation rate, are a useful technical contribution. However, the two load-bearing points—the thermalization of captured DM in Pop III stars and the consistency of the abstract's headline mass reach with the body's calculation—currently prevent the central claim from being fully established.","major_comments":[{"comment":"The thermalization timescale used to justify the Pop III exclusion regions is the degenerate-neutron-star expression from [19], with p_F ≈ 0.575 GeV in the numerator. The Pop III star is modeled as a non-degenerate n=3 polytrope with T_c = 2×10^7 K, for which the relevant target momentum is ∼(3 m_N T_c)^{1/2} ≈ 2×10^{-3} GeV. No justification is given for applying the degenerate formula to this environment, and the paper's own caveat in Section 3—'we assume that the DM particles have thermalized. If this is not the case, we can no longer guarantee that a BH would form, so we must lift our bounds'—means that if the thermalization time is underestimated, the headline exclusion at m_chi ≈ 10^{-8} GeV, sigma ≈ 10^{-40} cm^2 (Section 5) would be invalid. This is exactly the low-m_chi, high-sigma regime where the thermalization requirement is most restrictive. Please provide a non-degenerate thermalization calculation for Pop III stars or demonstrate that the excluded region is robust to the correct timescale.","section":"Section 3, Eq. (3.5); Section 5"},{"comment":"The abstract states that Pop III stars maintain 'efficacy below m_chi = 10^{-15} GeV' when a BEC forms, but Section 5 states explicitly that 'there must be a value of m_chi at which N_chi < N_Cha, which turns out to be sigma-independent and equal to approximately 10^{-8} GeV,' and the vertical line in the right panels of Figs. 4 and 5 is at m_chi ≈ 10^{-8} GeV. No curve in the body extends to 10^{-15} GeV. This is a direct contradiction between the headline claim and the body's own calculation; the abstract must be revised or the underlying calculation must be changed to support the claimed reach.","section":"Abstract vs. Section 5"}],"minor_comments":[{"comment":"The statement that the form factor 'suppresses the boundary of BH formation by approximately two orders of magnitude when m_chi > m_n' is ambiguous: the form factor suppresses the capture rate, which weakens (raises) the boundary; please rephrase to state that it weakens the bound.","section":"Section 2.1, after Eq. (2.5)"},{"comment":"The notation '0F1(;1+2/3 \\hat{\\phi}(r); \\tau(r))' should be typeset as _0F_1(; b; z) with the semicolon inside the argument; as written it appears to have an empty first argument, which is confusing.","section":"Section 2.2, Eq. (2.13)"},{"comment":"The MESA-derived stellar parameters (T_c = 2×10^7 K, R = 1.11 R_sun, M = 1000 M_sun) are stated without a reference or simulation details; please cite the MESA model or provide the relevant input file.","section":"Section 5"},{"comment":"The resulting 'closed-form approximations' still contain an unevaluated integral over xi (Xi(xi) in A.14 and the integral in A.18). Calling these 'closed-form' overstates the reduction; they are more accurately described as single-integral reductions of the original triple integral.","section":"Appendix A, Eqs. (A.14) and (A.18)"}],"recommendation":"major_revision","confidential_remarks":"The thermalization issue is the main technical risk: without a non-degenerate thermalization timescale for Pop III stars, the exclusion regions are conditional on an unjustified extrapolation. The abstract inconsistency is also problematic and should be corrected before the paper is reconsidered. The paper is well within the scope of JCAP and contains useful derivations; I see no need for rejection if these points are addressed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a credible extension of the scalar-ADM black-hole-formation program, and the Population III application is genuinely new. But the advertised low-mass reach rests on a thermalization timescale formula that was derived for degenerate neutron stars and applied to a non-degenerate hydrogen polytrope without re-derivation, and the abstract's 1e-15 GeV claim sits uneasily beside the body's sigma-independent 1e-8 GeV floor.\n\nWhat's actually good: the paper consistently combines multiscatter capture, the nucleon form-factor suppression, and evaporation, including a first approximate treatment of evaporation from a DM BEC. Appendix A derives closed-form evaporation rates for arbitrary polytropes, and the authors verify that the uniform-density radius approximation used for r_X is good to 1% over the parameter space. They also flag the thermalization assumption and mark where it forces the bounds to be lifted, which is honest.\n\nNow the soft spots, in rough order. First and most important: Eq. (3.5) is the degenerate-neutron-star thermalization formula from [19], with the Fermi momentum p_F in the numerator. Population III stars are non-degenerate; the target momentum scale is sqrt(m_n T_c), not p_F. The kinematics of DM-nucleon scattering in a thermal hydrogen plasma are not the same as scattering off a Fermi sea. The paper uses Eq. (3.5) for Pop III stars without justification. Since the claimed endpoint (mχ ≈ 1e-8 GeV, σ ≈ 1e-40 cm²) sits right where thermalization is most restrictive, and since BEC formation presupposes a thermalized population, this is load-bearing, not a footnote. The authors themselves say the bounds must be lifted if thermalization fails; they just haven't done the Pop III calculation.\n\nSecond, the abstract says 'efficacy below mχ=10^-15 GeV (assuming a BEC forms)', but Section 5 and Fig. 5 say the BEC-case floor is sigma-independent at about 10^-8 GeV. Maybe the 1e-15 number comes from an ultra-high-sigma regime where the evaporation suppression shuts off BEC evaporation, but the body doesn't connect those dots. As written, it reads as an internal inconsistency.\n\nThird, no code or data is shipped. Most of the machinery is imported from the authors' own previous papers, which is fine, but it means an independent check of the BEC evaporation and Knudsen-stitching approximations requires re-implementing a long chain.\n\nBottom line: the central idea is worth taking seriously and the paper deserves a proper referee. It should not be desk-rejected. But it needs a revision that either re-derives the thermalization timescale for non-degenerate stars or softens the Pop III claims, and reconciles the abstract with the body.","headline":"A credible extension of the scalar-ADM collapse program with a genuinely new Pop III application, but the headline low-mass reach rests on an unadapted thermalization formula and an abstract/body inconsistency.","tokens_in":18145,"tokens_out":3735,"would_cite":true,"duration_ms":36357,"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":"Observing a single intact Population III star can exclude bosonic asymmetric dark matter down to $10^{-8}$ GeV, or below $10^{-15}$ GeV if the captured dark matter forms a Bose-Einstein condensate.","keywords":["asymmetric dark matter","bosonic dark matter","black hole formation","neutron stars","Population III stars","dark matter capture","dark matter evaporation","Bose-Einstein condensate"],"falsifier":"Calculate the thermalization time from Eq. (3.5) for a $1000\\,M_\\odot$ Population III star of age $10^6$ years at the paper's deepest claimed exclusion points, for example $m_\\chi = 10^{-8}$ GeV with $\\sigma = 10^{-40}\\,\\mathrm{cm}^2$ and the BEC case below $10^{-15}$ GeV; if $t_{\\rm th}$ exceeds the stellar age anywhere in the excluded region, that part of the bound is not valid.","tokens_in":17048,"feed_emoji":"⭐","tokens_out":11922,"duration_ms":103102,"temperature":0.7,"pith_summary":"The paper argues that simply observing an intact star can rule out large regions of the allowed mass and interaction strength for bosonic asymmetric dark matter. The reason is that captured dark matter cannot annihilate, so it accumulates; if enough of it collects in the star's core it collapses into a black hole that swallows the star. The paper's new case for Population III stars is that a single high-mass first-generation star, the type easiest for JWST to see, should be destroyed if the dark matter has mass down to about $10^{-8}$ GeV at a scattering cross section of $10^{-40}\\,\\mathrm{cm}^2$, and down to $10^{-15}$ GeV if the captured dark matter forms a Bose-Einstein condensate. For neutron stars in very dense dark matter environments near the Galactic center, the paper shows that multiscatter capture and evaporation shift and reshape the previously derived exclusion boundary. If correct, these bounds extend dark matter searches into mass ranges far below what terrestrial direct detection can reach, using telescopes that already exist or are coming online.","feed_headline":"One first-generation star could rule out dark matter to 10^-8 GeV","feed_subtitle":"Captured bosonic dark matter would collapse the star into a black hole; its survival is the bound.","key_machinery":"The machinery is the population equation $dN_\\chi/dt = C - E N_\\chi$ together with two thresholds for black hole formation: $N_\\chi > N_{\\rm self}$ without a condensate, and $N_\\chi^0 > N_{\\rm Cha}^{\\rm boson}$ if a Bose-Einstein condensate forms. Capture is computed as a multiscatter sum over $N$ collisions, suppressed by a degenerate-neutron factor $\\xi_s$ at low mass and a nucleon form factor $F(\\delta p^2)$ at high mass; in four regions of parameter space the full sum reduces to closed-form rates used to locate analytic bounds. Evaporation is computed from the upscattering integral with a suppression factor $s_{\\rm evap}$, with isothermal and local-thermal-equilibrium dark matter distributions stitched by Knudsen number, plus a separate closed-form evaporation rate for particles in the BEC ground state. Polytropic Lane-Emden profiles provide the stellar density and temperature structure for both object classes.","core_discovery":"The central claim is that non-destruction of a neutron star or a Population III star excludes bosonic asymmetric dark matter whose parameters lie above the black-hole-formation boundary in the $\\sigma$--$m_\\chi$ plane. For a reference neutron star of $1.44\\,M_\\odot$, age $10^{10}$ years, in ambient dark matter density $\\rho_\\chi = 10^{13}\\,\\mathrm{GeV}\\,\\mathrm{cm}^{-3}$, the paper obtains the boundary including multiscatter capture, evaporation, and a finite-nucleon-size form factor; the main new qualitative feature is that the boundary enters multiscatter region I, so capture does not saturate until much higher $\\sigma$, allowing deeper $m_\\chi$ reach than previous single-scatter treatments. For a $1000\\,M_\\odot$ Population III star observed $10^6$ years after entering the main sequence, the paper claims bounds reaching $m_\\chi \\sim 10^{-8}$ GeV for $\\sigma = 10^{-40}\\,\\mathrm{cm}^2$, and below $10^{-15}$ GeV if a Bose-Einstein condensate forms, because the geometric capture limit is pushed to far higher $\\sigma$ than in neutron stars. Throughout, the paper treats dark matter evaporation self-consistently, including from the condensed ground state, and derives closed-form evaporation approximations for arbitrary polytropic objects and for a BEC. All of these bounds are explicitly lifted where the dark matter has not thermalized within the object's lifetime.","pith_inferences":["We infer that the same capture-and-evaporation machinery could be applied to white dwarfs and brown dwarfs, where the lower escape velocity makes evaporation stronger; that analysis would sharpen or relax the sub-GeV bounds those objects already provide, depending on how well the dark matter thermalizes.","We infer that if BEC-forming bosonic ADM exists near the lower end of the Population III reach, the first deep searches for Population III stars in dark-matter-rich minihalos should see an absence of such stars; a surviving star there would disfavor that candidate.","We infer that the thermalization requirement is the quiet bottleneck: future work should map the excluded region after imposing $t_{\\rm th} < t_\\ast$ exactly, since the closed-form estimate of Eq. (3.5) may be the main source of systematic uncertainty in the low-mass, high-$\\sigma$ corner."],"forward_implications":["If a single Population III star of about $1000\\,M_\\odot$ is observed at age $\\sim 10^6$ years, the absence of collapse excludes bosonic ADM with $\\sigma = 10^{-40}\\,\\mathrm{cm}^2$ down to $m_\\chi \\sim 10^{-8}$ GeV, and below $10^{-15}$ GeV if a Bose-Einstein condensate forms.","With many Population III stars observed, non-destruction in high-density early-universe minihalos would strengthen the exclusion beyond the single-object limit, subject to the thermalization condition.","Neutron star bounds in dense environments ($\\rho_\\chi \\gtrsim 10^9$ GeV cm$^{-3}$) must include multiscatter capture and evaporation; single-scatter-only treatments can misplace or overstate parts of the excluded region.","The closed-form evaporation rates for arbitrary polytropes and for BEC states let these bounds be recast quickly for other celestial objects without expensive numerical integration.","In the mass range probed, these stellar bounds are complementary to terrestrial direct searches, covering sub-keV masses below the neutrino-floor-limited reach of current detectors."],"supporting_citations":[{"why":"Supplies the baseline method for bosonic ADM collapse in neutron stars via $N_{\\rm self}$ and $N_{\\rm Cha}^{\\rm boson}$, which this paper extends.","marker":"[19]"},{"why":"Provides the multiscatter capture formalism in Eq. (2.3), summing capture over $N$ collisions.","marker":"[13]"},{"why":"Provides the form-factor suppression $F(\\delta p^2)$ for finite nucleon size used when $m_\\chi > m_n$.","marker":"[20]"},{"why":"Gives the per-particle evaporation rate $E$ and the evaporation mass scale about 1.4 keV for neutron stars.","marker":"[24]"},{"why":"Provides the evaporation suppression factor $s_{\\rm evap}$ and the Knudsen-number interpolation between isothermal and LTE dark matter distributions.","marker":"[25]"},{"why":"Sets the Population III star parameters, dark matter temperature, and first-star capture analysis that the Pop III bounds build on.","marker":"[9]"},{"why":"Gives the region-by-region closed-form multiscatter capture rates (regions I-IV) used to derive analytic bounds.","marker":"[7]"},{"why":"Source of the high Galactic-center dark matter density $\\rho_\\chi \\sim 10^{13}$ GeV cm$^{-3}$ scenario for neutron stars.","marker":"[28]"},{"why":"Provides the upscattering rate and evaporation integral formalism used throughout the evaporation calculation.","marker":"[22]"}],"fun_headline_variants":["Black hole formation in stars bounds asymmetric dark matter","Pop III stars could probe DM down to 10^-15 GeV","Neutron stars and Pop III stars as dark matter probes","Star survival excludes asymmetric dark matter"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The bounds rest on captured dark matter actually thermalizing inside the star within its lifetime; if a dark matter particle is captured but does not sink to the core, the paper cannot guarantee a black hole forms, and the exclusion region must be lifted.","fun_headline_variants_meta":{"raw":{"variants":["Black hole formation in stars bounds asymmetric dark matter","Pop III stars could probe DM down to 10^-15 GeV","Neutron stars and Pop III stars as dark matter probes","Star survival excludes asymmetric dark matter"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00076,"raw_usage":{"total_tokens":3430,"prompt_tokens":1054,"completion_tokens":2376,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":670,"completion_tokens_details":{"reasoning_tokens":2313}},"tokens_in":670,"tokens_out":2376,"duration_ms":17653,"temperature":1.0,"reasoning_tokens":2313,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T18:22:18.675294+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Calculate the thermalization time from Eq. (3.5) for a $1000\\,M_\\odot$ Population III star of age $10^6$ years at the paper's deepest claimed exclusion points, for example $m_\\chi = 10^{-8}$ GeV with $\\sigma = 10^{-40}\\,\\mathrm{cm}^2$ and the BEC case below $10^{-15}$ GeV; if $t_{\\rm th}$ exceeds the stellar age anywhere in the excluded region, that part of the bound is not valid.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the baseline method for bosonic ADM collapse in neutron stars via $N_{\\rm self}$ and $N_{\\rm Cha}^{\\rm boson}$, which this paper extends."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the form-factor suppression $F(\\delta p^2)$ for finite nucleon size used when $m_\\chi > m_n$."},{"cited_title":"Garani and S","cited_arxiv_id":null,"evidence_quote":"Gives the per-particle evaporation rate $E$ and the evaporation mass scale about 1.4 keV for neutron stars."},{"cited_title":"Garani and S","cited_arxiv_id":null,"evidence_quote":"Provides the evaporation suppression factor $s_{\\rm evap}$ and the Knudsen-number interpolation between isothermal and LTE dark matter distributions."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Sets the Population III star parameters, dark matter temperature, and first-star capture analysis that the Pop III bounds build on."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Source of the high Galactic-center dark matter density $\\rho_\\chi \\sim 10^{13}$ GeV cm$^{-3}$ scenario for neutron stars."},{"cited_title":"Gould,Weakly Interacting Massive Particle Distribution in and Evaporation from the Sun, ApJ321(1987) 560","cited_arxiv_id":null,"evidence_quote":"Provides the upscattering rate and evaporation integral formalism used throughout the evaporation calculation."}],"review_version":1}