Impact of Evaporation Barriers on Solar-Captured Dark Matter Distribution and Evaporation Mass

arxiv: 2509.26192 · v3 · submitted 2025-09-30 · ✦ hep-ph

Impact of Evaporation Barriers on Solar-Captured Dark Matter Distribution and Evaporation Mass

Xuan Wen This is my paper

Pith reviewed 2026-05-18 12:09 UTC · model grok-4.3

classification ✦ hep-ph

keywords solar captured dark matterevaporation barrierdark matter distributionevaporation massnon-thermal distributionorbit-space calculationin-medium attractionphase-space structure

0 comments p. Extension

The pith

An evaporation barrier from solar medium attraction suppresses dark matter evaporation by shifting the distribution toward tightly bound orbits.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper presents an orbit-space calculation of the non-thermal distribution of solar-captured dark matter when an evaporation barrier is present. The barrier deepens the effective potential and deforms near-threshold phase space, moving particles away from weakly bound escape trajectories and toward core-crossing orbits. This suppresses evaporation overall and lowers the evaporation mass that sets the low-mass reach of searches. Although the bulk population stays near thermal equilibrium, the near-threshold tail acquires non-thermal features visible in the projected velocity spectrum. The calculation shows that this tail controls the low-mass sensitivity in the barrier regime.

Core claim

The barrier not only deepens the effective potential but also reshapes the near-threshold phase-space structure, displacing the equilibrium distribution away from weakly bound, escape-prone trajectories and toward more tightly bound core-crossing orbits, thereby suppressing evaporation and lowering the evaporation mass.

What carries the argument

orbit-space calculation of the non-thermal distribution in the presence of a smooth in-medium attraction barrier that deforms bound orbits

If this is right

The low-mass reach of solar dark matter searches extends further because fewer particles evaporate.
The projected velocity spectrum develops characteristic non-thermal features near threshold.
Particles on repeated core-crossing orbits are preferentially retained over escape-prone ones.
The near-threshold tail becomes essential for accurate predictions in the barrier regime.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

Similar medium-induced barriers could modify capture and retention calculations for dark matter in other stars or planets.
Incorporating this orbit deformation might change expected event rates in low-mass direct detection experiments that rely on solar capture.
Extending the calculation to time-dependent or spatially varying barriers could test robustness against more realistic plasma conditions.
The approach offers a controlled way to include medium effects without full microphysical scattering simulations.

Load-bearing premise

The in-medium attraction can be treated as a smooth evaporation barrier that deforms bound orbit space in a controlled way.

What would settle it

A direct simulation or observation of the near-threshold velocity spectrum that fails to show the predicted non-thermal structure or shows no reduction in evaporation mass when medium effects are included.

Figures

Figures reproduced from arXiv: 2509.26192 by Xuan Wen.

**Figure 1.** Figure 1: The equilibrium distribution fχ(E, L) for mχ = 2.5 GeV at β = 0 (left) and β = 1 (right). The energy E and angular momentum L are expressed in units of GM⊙/R⊙ and p GM⊙R⊙, respectively. Only the colored parameter region corresponds to bound states. solutions, we solve the linear Boltzmann equation in the two-dimensional orbit space spanned by the specific energy E and angular momentum L, using a Monte Carl… view at source ↗

**Figure 2.** Figure 2: The integrated distribution function in velocity vχ for β = 0 in orbit (E = −1.175, L = 0.1223), and for β = 1, 2, 3 in orbit (E = −1.567, L = 0.1223). 0.0 0.2 0.4 0.6 0.8 0 1 2 3 4 0.0 0.2 0.4 0.6 0.8 0 1 2 3 4 5 6 0.0 0.2 0.4 0.6 0.8 0 1 2 3 4 0.0 0.2 0.4 0.6 0.8 0 1 2 3 4 5 [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗

**Figure 3.** Figure 3: The integrated distribution function in radius r for β = 0 in orbit (E = −1.175, L = 0.1223), and for β = 1, 2, 3 in orbit (E = −1.567, L = 0.1223). 4 [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗

**Figure 4.** Figure 4: The ratio of the non-thermal to the approximated thermal distribution for β = 0 at mχ = 1.5, 2.0, 2.5, 3.0, 3.5, and 4.0 GeV. 6 [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗

**Figure 5.** Figure 5: The ratio of the non-thermal to the approximated thermal distribution for β = 1 at mχ = 0.5, 1.0, 1.5, 2.0, 2.5, 3.0, 3.5, and 4.0 GeV. 7 [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗

**Figure 6.** Figure 6: The ratio of the non-thermal to the approximated thermal distribution for β = 2 at mχ = 0.5, 1.0, 1.5, 2.0, 2.5, 3.0, 3.5, and 4.0 GeV. 8 [PITH_FULL_IMAGE:figures/full_fig_p008_6.png] view at source ↗

**Figure 7.** Figure 7: The ratio of the non-thermal to the approximated thermal distribution for β = 3 at mχ = 0.5, 1.0, 1.5, 2.0, 2.5, 3.0, 3.5, and 4.0 GeV. 9 [PITH_FULL_IMAGE:figures/full_fig_p009_7.png] view at source ↗

**Figure 8.** Figure 8: The parameter regions in the (mχ, σp) plane for four benchmarks. The near-equilibrium regions (tanh(t⊙/τe) ≃ 1) are shown as the darker-colored areas, while the lighter shades correspond to 0.9 ≤ tanh(t⊙/τe) ≤ 1. In the red (blue) regions, evaporation (annihilation) is subdominant in setting the solar DM population, while the purple band marks the crossover where the two effects are of similar importance. … view at source ↗

read the original abstract

Evaporation determines the low-mass reach of solar-captured dark matter because that reach is controlled by the small population of particles closest to the escape threshold. We present an orbit-space calculation of the non-thermal distribution of captured dark matter in the presence of an evaporation barrier generated by a smooth in-medium attraction sourced by the solar medium. We show that the barrier not only deepens the effective potential but also reshapes the near-threshold phase-space structure, displacing the equilibrium distribution away from weakly bound, escape-prone trajectories and toward more tightly bound core-crossing orbits, thereby suppressing evaporation and lowering the evaporation mass. Although the bulk population remains near thermal equilibrium, the near-threshold tail, as reflected in the projected velocity spectrum, acquires characteristic non-thermal structure because the barrier deforms the bound orbit space and preferentially retains particles that repeatedly traverse the hot solar core. The near-threshold tail is therefore essential for determining the low-mass reach of solar dark-matter searches in the barrier regime, and our orbit-space treatment captures the relevant physics in a controlled way.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper's orbit-space treatment shows how a smooth evaporation barrier displaces near-threshold dark matter orbits toward core-crossing ones, which could suppress evaporation and shift the low-mass cutoff for solar capture searches.

read the letter

The key takeaway is that this work treats the barrier as deforming bound orbit phase space rather than just adding a deeper potential well. That displacement moves the equilibrium tail away from escape-prone trajectories and toward more bound core-crossing paths, which the authors argue suppresses evaporation and lowers the evaporation mass. The bulk distribution stays close to thermal while the near-threshold part picks up non-thermal structure from the deformed density of states. This is the concrete new element relative to earlier thermal models of solar DM capture and evaporation. The calculation builds directly on standard capture assumptions without introducing obvious circularity, and the abstract frames the orbit-space step as a controlled approximation. That framing is reasonable on its face and gives a clear mechanism for why the tail matters for low-mass reach in neutrino and direct-detection limits. A clear limitation is the absence of any reported numbers, error bars, or direct comparison to Monte Carlo orbit simulations in the material I have. Without those, it is hard to judge how large the suppression actually is or how sensitive it is to the barrier depth scale. The stress-test concern about back-reaction on the scattering kernel is also worth checking in the full derivation: if the barrier alters relaxation rates, the assumed steady-state orbit distribution may not hold. Those points are fixable with additional plots or appendices rather than fatal. The paper is aimed at people who set limits on sub-GeV dark matter using solar capture signals. A reader working on evaporation cutoffs or on the mapping from capture rate to observable flux will get direct use from the orbit-space picture even if the quantitative shift turns out modest. It is coherent enough on its own terms to merit a serious referee who can check the orbit deformation math and ask for the missing benchmarks.

Referee Report

2 major / 2 minor

Summary. The paper claims that a smooth in-medium attractive potential in the Sun acts as an evaporation barrier that both deepens the effective potential and deforms the near-threshold bound-orbit phase space of captured dark matter. This deformation shifts the equilibrium distribution away from weakly bound, escape-prone trajectories toward more tightly bound core-crossing orbits, suppressing evaporation, lowering the evaporation mass, and imprinting non-thermal structure on the near-threshold velocity tail while leaving the bulk population near-thermal. The orbit-space treatment is presented as capturing the relevant physics in a controlled manner for determining the low-mass reach of solar DM searches.

Significance. If the central result holds, the work would adjust the low-mass cutoff for solar-captured dark matter searches by reducing the evaporation rate through phase-space reshaping. The emphasis on the non-thermal tail and the orbit-space formalism provides a potentially useful refinement over purely thermal models, particularly if the barrier effect can be shown to be robust against microphysical details.

major comments (2)

Abstract (paragraph on orbit-space calculation): the central claim that the barrier deforms bound orbit space and displaces the distribution toward core-crossing orbits while preserving the same relaxation processes requires demonstration that the barrier-induced force does not modify the scattering kernel or capture/thermalization rates; without this, the assumed orbit-space equilibrium may not be the correct steady state.
Abstract: no quantitative results, error estimates, or comparison to simulations are provided to support the magnitude of the evaporation suppression or the shift in evaporation mass, leaving the size of the effect on the low-mass reach unassessed.

minor comments (2)

Abstract: the description of the 'projected velocity spectrum' acquiring non-thermal structure would benefit from a brief definition or reference to how the projection is performed.
Abstract: clarify the origin and range of the barrier depth scale, which appears as an input parameter in the model.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for the constructive comments, which have helped us clarify key aspects of the orbit-space treatment. We respond to each major comment below.

read point-by-point responses

Referee: Abstract (paragraph on orbit-space calculation): the central claim that the barrier deforms bound orbit space and displaces the distribution toward core-crossing orbits while preserving the same relaxation processes requires demonstration that the barrier-induced force does not modify the scattering kernel or capture/thermalization rates; without this, the assumed orbit-space equilibrium may not be the correct steady state.

Authors: We agree that an explicit justification is needed. In the model the in-medium attraction is introduced as a smooth, long-range mean-field potential that augments the effective gravitational potential governing orbital motion. The scattering kernel itself is determined by short-range, microscopic interactions with individual solar nuclei; because these interactions occur on length scales much smaller than the spatial variation of the barrier, their rates and angular distributions remain unchanged. We have added a dedicated paragraph in the revised Section 2 that spells out this scale separation and confirms that the relaxation processes entering the orbit-space master equation are unaffected. The steady-state distribution is therefore still obtained by balancing the same capture and scattering rates against the modified escape probability. revision: yes
Referee: Abstract: no quantitative results, error estimates, or comparison to simulations are provided to support the magnitude of the evaporation suppression or the shift in evaporation mass, leaving the size of the effect on the low-mass reach unassessed.

Authors: The referee correctly observes that the abstract is qualitative. The body of the paper already contains the orbit-space integrals that yield the modified distribution and the resulting evaporation rate; however, to make the magnitude of the effect transparent we have inserted a new subsection (now Section 4.3) that reports numerical values for the evaporation-mass shift and the fractional suppression of the evaporation rate relative to the barrier-free case. These numbers are accompanied by a brief error budget arising from the assumed smoothness of the potential and from the truncation of the orbit-space basis. While a full Monte-Carlo simulation of the solar interior lies outside the present scope, we compare the orbit-space results directly to the standard thermal evaporation formula and discuss the implications for the low-mass cutoff of solar DM searches. revision: yes

Circularity Check

0 steps flagged

No circularity: orbit-space calculation is independent of fitted inputs or self-referential loops

full rationale

The paper introduces an orbit-space treatment of solar-captured DM in the presence of a smooth in-medium evaporation barrier. The abstract and provided text describe the barrier deepening the effective potential and deforming near-threshold phase space to shift the distribution toward core-crossing orbits, suppressing evaporation. No equations, parameter fits, or self-citations are exhibited that reduce this suppression effect or the non-thermal tail structure to a quantity defined or fitted within the same work. The derivation builds on standard solar DM capture assumptions without the central claim reducing by construction to those inputs or to a prior result by the same author. The calculation is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 1 invented entities

The work relies on standard solar DM capture assumptions and introduces the barrier as a new modeling element without external validation.

free parameters (1)

Barrier depth scale
Smooth in-medium attraction strength that sets the barrier height; appears required to produce the reported suppression.

axioms (2)

domain assumption Dark matter particles lose energy via scattering with solar nuclei and can be gravitationally captured.
Invoked to establish the initial captured population before barrier effects are applied.
ad hoc to paper The solar medium generates a smooth attractive potential that acts as an evaporation barrier.
Central modeling choice introduced to deform the orbit space.

invented entities (1)

Evaporation barrier no independent evidence
purpose: To reshape near-threshold phase space and suppress evaporation
New postulated effect from in-medium attraction; no independent falsifiable signature outside this calculation is provided.

pith-pipeline@v0.9.0 · 5704 in / 1490 out tokens · 41503 ms · 2026-05-18T12:09:01.548107+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

We solve the linear Boltzmann equation in the two-dimensional orbit space spanned by the specific energy E and angular momentum L, using a Monte Carlo scheme... Escape is implemented through the barrier-modified condition determined by β

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

25 extracted references · 25 canonical work pages · 1 internal anchor

[1]

J.321560

Gould A 1987Astrophys. J.321560

work page
[2]

Abbasi Ret al.(IceCube) 2022Phys. Rev. D105062004

work page
[3]

Bell N F, Dolan M J and Robles S 2021JCAP2021004

work page
[4]

Hochberg Y, Zhao Y and Zurek K M 2016Phys. Rev. Lett.116011301

work page
[5]

Schutz K and Zurek K M 2016Phys. Rev. Lett.117121302

work page
[6]

Hochberg Y, Pyle M, Zhao Y and Zurek K M 2016JHEP08057

work page
[7]

Hochberg Y, Charaev I, Nam S W, Verma V, Colangelo M and Berggren K K 2019Phys. Rev. Lett.123151802

work page
[8]

Chiles J, Charaev I, Lasenby R, Baryakhtar M, Huang J, Roshko A, Burton G, Colangelo M, Van Tilburg K, Arvanitaki A, Nam S W and Berggren K K 2022Phys. Rev. Lett.128231802

work page
[9]

Hochberg Y, Kahn Y, Kurinsky N, Lehmann B V, Yu T C and Berggren K K 2021Phys. Rev. Lett.127151802

work page
[10]

Astrophys.480847

Orlando E and Strong A W 2008Astron. Astrophys.480847

work page
[11]

Linden T, Beacom J F, Peter A H G, Buckman B J, Zhou B and Zhu G 2022Phys. Rev. D 105063013

work page
[12]

Bell N F, Dent J B and Sanderson I W 2021Phys. Rev. D104023024

work page
[13]

Niblaeus C, Beniwal A and Edsj¨ o J 2019JCAP2019011

work page
[14]

Feng J L, Smolinsky J and Tanedo P 2016Phys. Rev. D93115036

work page
[15]

Nguyen T T Q, Linden T, Carenza P and Widmark A 2025 arXiv:2501.14864

work page internal anchor Pith review Pith/arXiv arXiv 2025
[16]

Acevedo J F, Leane R K and Smirnov J 2024JCAP2024038 12 IOP PublishingJournalvv(yyyy) aaaaaa Authoret al

work page
[17]

Liang Z L, Wu Y L, Yang Z Q and Zhou Y F 2016JCAP2016018

work page
[18]

Liang Z L, Tang Y L and Yang Z Q 2018JCAP2018035

work page
[19]

J.294663–673

Spergel D N and Press W H 1985Astrophys. J.294663–673

work page
[20]

Serenelli A, Basu S, Ferguson J W and Asplund M 2009Astrophys. J. Lett.705L123–L127

work page
[21]

J.328919–939

Gould A 1988Astrophys. J.328919–939

work page
[22]

J.321571

Gould A 1987Astrophys. J.321571

work page
[23]

J.388338

Gould A 1992Astrophys. J.388338

work page
[24]

Busoni G, De Simone A and Huang W C 2013JCAP2013010

work page
[25]

Garani R and Palomares-Ruiz S 2017JCAP2017007 13

work page

[1] [1]

J.321560

Gould A 1987Astrophys. J.321560

work page

[2] [2]

Abbasi Ret al.(IceCube) 2022Phys. Rev. D105062004

work page

[3] [3]

Bell N F, Dolan M J and Robles S 2021JCAP2021004

work page

[4] [4]

Hochberg Y, Zhao Y and Zurek K M 2016Phys. Rev. Lett.116011301

work page

[5] [5]

Schutz K and Zurek K M 2016Phys. Rev. Lett.117121302

work page

[6] [6]

Hochberg Y, Pyle M, Zhao Y and Zurek K M 2016JHEP08057

work page

[7] [7]

Hochberg Y, Charaev I, Nam S W, Verma V, Colangelo M and Berggren K K 2019Phys. Rev. Lett.123151802

work page

[8] [8]

Chiles J, Charaev I, Lasenby R, Baryakhtar M, Huang J, Roshko A, Burton G, Colangelo M, Van Tilburg K, Arvanitaki A, Nam S W and Berggren K K 2022Phys. Rev. Lett.128231802

work page

[9] [9]

Hochberg Y, Kahn Y, Kurinsky N, Lehmann B V, Yu T C and Berggren K K 2021Phys. Rev. Lett.127151802

work page

[10] [10]

Astrophys.480847

Orlando E and Strong A W 2008Astron. Astrophys.480847

work page

[11] [11]

Linden T, Beacom J F, Peter A H G, Buckman B J, Zhou B and Zhu G 2022Phys. Rev. D 105063013

work page

[12] [12]

Bell N F, Dent J B and Sanderson I W 2021Phys. Rev. D104023024

work page

[13] [13]

Niblaeus C, Beniwal A and Edsj¨ o J 2019JCAP2019011

work page

[14] [14]

Feng J L, Smolinsky J and Tanedo P 2016Phys. Rev. D93115036

work page

[15] [15]

Nguyen T T Q, Linden T, Carenza P and Widmark A 2025 arXiv:2501.14864

work page internal anchor Pith review Pith/arXiv arXiv 2025

[16] [16]

Acevedo J F, Leane R K and Smirnov J 2024JCAP2024038 12 IOP PublishingJournalvv(yyyy) aaaaaa Authoret al

work page

[17] [17]

Liang Z L, Wu Y L, Yang Z Q and Zhou Y F 2016JCAP2016018

work page

[18] [18]

Liang Z L, Tang Y L and Yang Z Q 2018JCAP2018035

work page

[19] [19]

J.294663–673

Spergel D N and Press W H 1985Astrophys. J.294663–673

work page

[20] [20]

Serenelli A, Basu S, Ferguson J W and Asplund M 2009Astrophys. J. Lett.705L123–L127

work page

[21] [21]

J.328919–939

Gould A 1988Astrophys. J.328919–939

work page

[22] [22]

J.321571

Gould A 1987Astrophys. J.321571

work page

[23] [23]

J.388338

Gould A 1992Astrophys. J.388338

work page

[24] [24]

Busoni G, De Simone A and Huang W C 2013JCAP2013010

work page

[25] [25]

Garani R and Palomares-Ruiz S 2017JCAP2017007 13

work page