{"id":"f49053cd-5dcd-43ae-8b45-d49cbaf2dcdb","arxiv_id":"2412.17919","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Fermionic dark matter halos around supermassive black holes can form non-power-law spikes or even density depletions, depending on the fermion mass and degeneracy, instead of the universal r^-3/2 spike of power-law halos.","lead":"The paper computes how dark matter made of heavy fermions rearranges around a supermassive black hole, and finds that the resulting spike shape depends strongly on the fermion mass. In some cases the black hole depletes, rather than enhances, the surrounding dark matter, which would change how dark-matter signals and gravitational waves from galactic centers are interpreted.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Neglect of DM self-gravity in the final state (Sec. IV) undermines the novel light-fermion spikes, since M_DM ~ M_BH precisely in the regime where enhancement and mass-dependence are claimed.","rationale":"The reader's weakest assumption is exactly the one I would flag. The paper is transparent about this limitation, which is a point in its favor, but the limitation sits at the center of the most novel results. The calculation of the final density (Eq. 3) uses the Schwarzschild metric for the grown BH, and the action-conservation mapping (Appendix A) connects the initial RAR metric to this pure-Schwarzschild final metric. When the DM mass enclosed in the spike region is comparable to M_BH, the true final metric must include that mass; ignoring it changes both the effective potential and the radial action, so the computed ρ(r) is not the self-consistent adiabatic outcome. The authors identify the affected regime themselves: mc^2 ≲ 250 keV and M_BH ≲ 4×10^6 M_sun, which is where the new non-power-law, mass-dependent spikes appear. Thus the headline claim is not yet secure. A focused recalculation with a self-consistent final metric for one representative light-fermion model would settle whether the qualitative conclusions survive. If they do, the paper's conclusions are strengthened; if not, the claim should be restricted to heavy fermions or to cases where the DM mass is subdominant. My recommended verdict is unchanged from the reader's CONDITIONAL, because the concern is the same and the paper already acknowledges it; no new objection beyond that is needed.","tokens_in":19135,"tokens_out":10560,"duration_ms":106873,"concrete_test":"Recompute the m = 100 keV, M_BH = 4×10^6 M_sun model (Fig. 2, middle) including the DM self-gravity in the final state: (1) evaluate M_DM(r) inside the spike from the published density profile; (2) update the final metric to g00(r) = 1 − 2[M_BH + M_DM(r)]/r (or solve the TOV equation with the DM contribution) and redo the action-conservation mapping and Eq. (3); (3) iterate until M_DM(r) converges. If the peak density or the power-law slope changes by more than ~10% in the region where M_DM ~ M_BH, the Schwarzschild-only spike is not self-consistent and the central claim is not robust. A simple diagnostic: check whether M_DM(<r_peak)/M_BH > 0.1; if so, the test is triggered.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—that fermionic spikes are non-power-law and can deplete rather than enhance—rests on the adiabatic remapping of orbits in a final spacetime taken to be pure Schwarzschild (Sec. II A, Eq. 2; Appendix A). The authors explicitly state (Sec. IV) that for fermions with mc^2 ≲ 250 keV and M_BH ≲ 4×10^6 M_sun, the enclosed DM mass inside the BH sphere of influence is comparable to M_BH. This is precisely the parameter region where the paper reports the novel enhanced, mass-dependent spikes (e.g., m = 100 and 250 keV curves in Figs. 2 and 3). If the remaining DM mass is not negligible, the final metric is not Schwarzschild; the effective potential V_eff (Eq. A6), the radial action (Eq. A7), and the energy mapping (Eq. A9) all change, and Eq. (3) as evaluated is not a self-consistent solution for the spike. Thus the quantitative predictions of the light-fermion spikes—and the claimed mass/degeneracy dependence—are not established. The authors' suggestion that this favors m ≳ 300 keV does not rescue the general claim because the paper's stated novelty is the general mass-dependent behavior, not just heavy fermions.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper computes relativistic dark-matter (DM) spike profiles around a Schwarzschild supermassive black hole (SMBH) that grows in the center of a fermionic DM halo described by the Ruffini-Argüelles-Rueda (RAR) model. In Model I, a small baryonic BH seed grows adiabatically and the final DM distribution is obtained by conserving radial action and angular momentum between the initial RAR phase-space distribution and final bound orbits around the BH (Sec. II A and Appendix A). In Model II, a dense fermion core grows by accreting baryons to a critical mass and then collapses suddenly to an SMBH, with the remaining particles remapped according to the instantaneous-collapse prescription of Appendix B. The principal claims are that fermionic spikes generally do not follow a simple power law, that the spike structure depends on the fermion mass and degeneracy state, that a sufficiently massive BH or heavy fermion mass can suppress rather than enhance the central DM density, and that the usual r^{-3/2} spike is recovered only in the Boltzmannian regime.","tokens_in":19440,"tokens_out":9422,"duration_ms":101740,"significance":"The paper is a genuine forward calculation, not a fit: it starts from a specified RAR equilibrium distribution, applies radial-action conservation with published numerical methods, and recovers the standard Gondolo-Silk r^{-3/2} spike in the low-degeneracy limit. If the results are correct, they would be relevant for indirect DM searches, gravitational-wave signatures of DM spikes, stellar-orbit constraints near Sgr A*, and the DM-origin SMBH-seed scenario. The numerical implementation is described in enough detail to be reproducible. However, the authors themselves identify the key weakness: in the light-fermion regime that produces the novel non-power-law spikes, the enclosed DM mass is comparable to the BH mass, so the neglect of DM self-gravity in the final state is not justified. That limitation directly affects the main quantitative novelty of the paper.","major_comments":[{"comment":"The central claim of mass- and degeneracy-dependent, non-power-law spikes for semi-degenerate fermions is computed in a regime where the authors themselves state that the enclosed DM mass is comparable to the BH mass. For fermions with mc^2 ≲ 250 keV and M_BH ≲ 4×10^6 M_sun, the final gravitational potential is not the pure Schwarzschild potential used in Eq. (3), and the effective potential (A6), radial action (A7), and energy mapping (A9) are all evaluated with the DM self-gravity omitted. The 100 keV and 250 keV curves in Figs. 2 and 3, which exhibit the strongest mass dependence and the claimed depletion/enhancement transition, are therefore not self-consistent solutions. The suggestion in Sec. IV that the approximation favors m ≳ 300 keV does not rescue the general abstract claim, because the novel light-fermion behavior is precisely what is not established. The authors should either recompute the spike with a metric that includes the remaining DM mass self-consistently, or explicitly restrict all quantitative claims to parameters for which M_DM << M_BH is verified for every displayed profile.","section":"Sec. IV; Sec. II A, Eq. (3); Appendix A, Eqs. (A6)-(A9)"},{"comment":"The sudden-collapse mapping in Model II is not sufficiently validated. Equation (B2) assumes that the particle velocity is unchanged at the same coordinate radius during the instantaneous transition from the critical core to a Schwarzschild BH, but the collapse from R_crit = 9.3 M_crit to the BH horizon occurs on a timescale comparable to orbital periods of particles near that radius. The resulting final density therefore depends on this unverified approximation. Because Model II is presented as an independent scenario for SMBH formation of DM origin, the authors should either demonstrate robustness against a finite-time collapse prescription or explicitly label the Model II spike as schematic pending a dynamical treatment.","section":"Sec. III B; Appendix B, Eq. (B2)"},{"comment":"The statement that 'for fermion masses above ~300 keV, there is no enhancement of the density of the DM surrounding the BH for any BH mass, but a suppression' is presented as a general result, but it is derived from a single family of RAR models with fixed core mass (M_c = 3.5×10^6 M_sun) and fixed outer boundary conditions. The threshold mass and the depletion behavior could depend on the chosen halo parameters. The claim should be either demonstrated across the allowed parameter space of the RAR model or restated as a property of the specific configurations studied here.","section":"Sec. IV, item 3; Fig. 3"}],"minor_comments":[{"comment":"The title line contains a typo: 'black hol e seeds' should read 'black hole seeds'.","section":"Title"},{"comment":"The variable r in Eq. (6) is used interchangeably with the dimensionless ratio r/M_BH; please state explicitly that the integration limits and energy thresholds are given in units of the black-hole mass.","section":"Sec. II A, Eq. (6)"},{"comment":"The phrase 'critic mass' appears twice and should be corrected to 'critical mass'.","section":"Sec. II C"},{"comment":"The probability in Eq. (B1) is written as P(r_f|E_i,L), while Eq. (B3) uses P(r_f|r_i,v_i); the notation should be made consistent and the conditioning variables clearly defined.","section":"Appendix B, Eqs. (B1) and (B3)"},{"comment":"The relation between the initial core mass M_dm = 1.5×10^6 M_sun, the baryon-to-DM fraction χ = 0.5, and the resulting critical mass M_crit = 4×10^6 M_sun is not transparent. Please specify whether M_dm in Eq. (17) is the initial pure-DM core mass or the core mass at the time of collapse, and how M_crit is obtained from χ and M_dm.","section":"Sec. II C and Sec. III B"}],"recommendation":"major_revision","confidential_remarks":"The paper is a well-structured extension of the Sadeghian-Ferrer-Will method to RAR fermionic halos, and the authors are transparent about the main approximation. The decisive issue is whether the final-state neglect of DM self-gravity can be controlled or repaired; if the authors restrict their quantitative claims to m ≳ 300 keV or implement a self-consistent BH+DM metric, the paper could be publishable. I would also ask the referee handling the paper to verify that the imported RAR halo parameters and critical masses from the authors' earlier work are used within their stated domain of validity."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this is a serious, well-written application of the relativistic adiabatic spike machinery to RAR fermionic core-halo halos, and it shows something genuinely new—the r^-3/2 spike is not universal, and a central BH can deplete the surrounding DM instead of enhancing it. But the paper's general claim overreaches: the authors' own admitted neglect of DM self-gravity in the final state breaks the calculation for the light-fermion cases (mc^2 ≲ 250 keV) that anchor the mass-dependence, so the 100–250 keV curves should be treated as uncomputed, not as predictions.\n\nWhat is new: applying the Sadeghian et al. GR treatment to finite-temperature, self-gravitating fermion halos. The equations are laid out in enough detail to be checked. In the Boltzmannian limit they recover G&S r^-3/2, which is a good sanity check. The depletion threshold—for a given fermion mass there is a BH mass above which density is suppressed, and for m≳300 keV no enhancement for any BH mass—is a real qualitative result, at least for heavy fermions. Model II, the sudden collapse of a critical DM core, is an interesting new scenario and the authors attempt to model the orbital reshuffling.\n\nSoft spots: the self-gravity neglect is the load-bearing one. They write in Sec. IV that for m≲250 keV and M_BH≲4×10^6 M_sun, the enclosed DM mass is comparable to the BH mass. In that regime the final metric is not Schwarzschild; the effective potential (A6), radial action (A7), and energy relation (A9) all change, so Eq. (3) is not a self-consistent spike. The authors acknowledge this and conclude that m≳300 keV is favored—which is a retreat from the abstract's 'general case.' It means the novel non-power-law, mass-dependent spikes are only established for the heavy-fermion corner. Also, no code or data are released, so the numerical results are not independently reproducible. Model II leans on imported critical masses from prior group papers and a sudden-collapse mapping that is heuristic; that is defensible as prior work, but the mapping itself is less tested than the adiabatic one.\n\nThis paper deserves a serious referee. The right outcome is a conditional accept with the authors forced to either restrict the claims to the regime where self-gravity is negligible, or compute the spike with the full self-gravitating final potential. A code release would also help. I'd bring it to a reading group; the central idea trains a useful light on how DM particle physics changes spike predictions.","headline":"A serious relativistic treatment of dark-matter spikes for fermionic core-halo halos that shows depletion is real, but the light-fermion curves are not self-consistent because the final DM self-gravity is neglected.","tokens_in":19948,"tokens_out":3497,"would_cite":true,"duration_ms":33729,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["83C57","85A05"],"pacs":["95.35.+d","04.70.-s","98.62.Js"],"model":"deepseek-v4-flash","headline":"This paper claims that once dark matter is treated as massive fermions in equilibrium, the spike around a supermassive black hole is generally not a power law, and for fermion masses above roughly 300 keV the black hole depletes the…","keywords":["dark matter","fermionic dark matter","dark matter spikes","supermassive black holes","RAR model","general relativity","adiabatic growth","core-halo profiles"],"falsifier":"Measure the density profile within the sphere of influence of a supermassive black hole such as Sgr A*: a clean $r^{-3/2}$ power law extending to the innermost bound orbit would contradict the predicted depletion for fermion masses above roughly 300 keV, while a non-power-law profile with a suppressed central density would support the fermionic core-halo prediction.","tokens_in":18903,"feed_emoji":"🕳️","tokens_out":8615,"duration_ms":73420,"temperature":0.7,"pith_summary":"The paper asks what happens to the dark matter overdensity around a supermassive black hole when the dark matter is a self-gravitating gas of massive fermions, rather than a featureless power-law halo. It finds a mass- and degeneracy-dependent outcome: only dilute, Boltzmannian fermions produce the familiar $\\rho \\sim r^{-3/2}$ spike, while semi-degenerate core-halo fermions develop non-power-law spikes or even a density depletion. The significance is that dark-matter annihilation signals, gravitational-wave dephasing, and stellar-orbit fits all assume a universal power-law spike. The paper also compares two histories of the black-hole seed, a small baryonic seed that grows adiabatically and a heavy seed formed by the baryon-triggered collapse of a fermionic core, and shows that the spike carries information about both the particle mass and the seed's origin.","feed_headline":"Black holes can deplete dark matter, not only boost it","feed_subtitle":"For heavy fermions the central density falls instead of rising, so spike probes depend on the particle's mass.","key_machinery":"The machinery is the adiabatic-invariant method in Schwarzschild spacetime applied to initial states generated by the RAR model, a relativistic self-gravitating Fermi gas with finite temperature, tidal cutoff, and a dense degenerate core surrounded by a dilute halo. In the baryonic-seed scenario, radial action and angular momentum are conserved during slow black-hole growth, so the initial distribution function maps to the final one; in the dark-matter-seed scenario, the collapse of the critical fermion core is treated as an instantaneous change of the metric potential with angular momentum conserved. The load-bearing identity is the relation between initial and final energies defined by equality of the radial actions, which converts the spike calculation into a phase-space integral over bound orbits, with the fermion mass $m$ and degeneracy parameter $\\theta_0$ entering through the initial RAR halo.","core_discovery":"The central claim is that the standard $r^{-3/2}$ dark-matter spike is a special case, not the generic outcome. Starting from equilibrium fermionic halos of the RAR model and growing a Schwarzschild black hole adiabatically, the authors find that core-halo configurations redistribute into spike profiles that depend on fermion mass and central degeneracy and do not follow a simple power law. For a fixed fermion mass there is a black-hole mass above which the spike density is depleted relative to the initial halo, and for fermion masses above roughly 300 keV no black-hole mass yields any enhancement. In the alternative seed scenario, where the black hole forms by the collapse of a dense fermion core, the same qualitative behavior holds with an additional dependence on the baryon-to-dark-matter fraction $\\chi$ inside the collapsing core. Dilute Boltzmannian fermions reproduce the classical $r^{-3/2}$ spike, matching earlier treatments.","pith_inferences":["Extension: the same mass-dependent logic would predict that fermionic cores around stellar-mass black holes and neutron stars produce smaller or absent spikes, potentially altering dark-matter dynamical-friction and X-ray binary constraints.","Extension: applying the same phase-space machinery to bosonic or solitonic central cores would test whether non-power-law, mass-dependent spike shapes are generic to self-gravitating central cores or specific to fermionic degeneracy pressure.","Extension: if heavy fermions deplete the spike, gamma-ray constraints from the Galactic center would be weaker for that mass range; a detected non-power-law emission morphology could then be turned into a fermion-mass measurement.","Extension: the dependence of the final spike on the baryon fraction $\\chi$ suggests that observed spike profiles could encode the accretion history of the black-hole seed, not only the particle mass."],"forward_implications":["Annihilation-flux predictions, gravitational-wave estimates, and stellar-orbit constraints that assume a universal $r^{-3/2}$ spike need revision when the host halo is a fermionic core-halo configuration.","Below the enhancement threshold the spike peak height stays roughly constant while its width shrinks as the black hole grows; above the threshold the density falls below the initial halo value.","For fermion masses above roughly 300 keV, a supermassive black hole suppresses rather than enhances the surrounding dark matter density for every black-hole mass considered.","The two seed scenarios give nearly identical spike shapes for the same final black-hole mass, differing mainly by the flattened peak left when the collapsed core is removed.","Spike observables become probes of the seed's origin and of the halo's degeneracy state, not just of the black-hole mass."],"supporting_citations":[{"why":"Defines the adiabatic spike scenario for a seed black hole growing in a dark-matter halo, the baseline prediction this paper generalizes.","marker":"[1]"},{"why":"Supplies the fully relativistic treatment of spike formation around a Schwarzschild black hole that the paper follows.","marker":"[4]"},{"why":"Introduces the RAR core-halo model used to construct the initial fermionic halo profiles.","marker":"[38]"},{"why":"Provides Milky-Way halo parameters and constraints used to set the boundary conditions for the host halos.","marker":"[39]"},{"why":"Establishes the formation and stability of fermionic core-halo halos in a cosmological framework, the physical basis of the initial states.","marker":"[40]"},{"why":"Details the baryon-induced collapse of a fermion core into a supermassive black hole, the basis of the paper's model II.","marker":"[48]"},{"why":"Shows how the collapse-produced heavy seeds grow to supermassive black holes, providing the growth context for model II.","marker":"[47]"},{"why":"Supplies the NFW power-law halo used as the comparison baseline for the classical spike.","marker":"[62]"}],"fun_headline_variants":["Dark matter spikes can vanish around heavy black holes","Fermion mass decides if black holes boost or deplete dark matter","Black hole seeds: dark matter spikes depend on particle mass","New model: black holes can strip dark matter, not just pile it up","Dark matter spikes: not always a power law, sometimes deplete"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that once the black hole is present, the remaining dark matter inside the spike feels only the black hole's gravity, with the self-gravity of the dark matter itself neglected, a regime the authors note is questionable for fermion masses below about 250 keV and black-hole masses below about $4\\times10^6$ solar masses.","fun_headline_variants_meta":{"raw":{"variants":["Dark matter spikes can vanish around heavy black holes","Fermion mass decides if black holes boost or deplete dark matter","Black hole seeds: dark matter spikes depend on particle mass","New model: black holes can strip dark matter, not just pile it up","Dark matter spikes: not always a power law, sometimes deplete"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000226,"raw_usage":{"total_tokens":1541,"prompt_tokens":1088,"completion_tokens":453,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":704,"completion_tokens_details":{"reasoning_tokens":365}},"tokens_in":704,"tokens_out":453,"duration_ms":4211,"temperature":1.0,"reasoning_tokens":365,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T05:08:29.075816+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the density profile within the sphere of influence of a supermassive black hole such as Sgr A*: a clean $r^{-3/2}$ power law extending to the innermost bound orbit would contradict the predicted depletion for fermion masses above roughly 300 keV, while a non-power-law profile with a suppressed central density would support the fermionic core-halo prediction.","supporting_citations":[{"cited_title":"This is markedly different from the spike density profile from power-law DM density halos, which depends only on the initial power-law index γ","cited_arxiv_id":null,"evidence_quote":"Defines the adiabatic spike scenario for a seed black hole growing in a dark-matter halo, the baseline prediction this paper generalizes."},{"cited_title":"In the RAR case, when the central fermions are in a state of low degeneracy (θ0 ≪ − 1), we recover the resulting DM spike of ρ ∝ r−3/2 typical of cored halos as the NSIS (see Fig","cited_arxiv_id":null,"evidence_quote":"Supplies the fully relativistic treatment of spike formation around a Schwarzschild black hole that the paper follows."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Details the baryon-induced collapse of a fermion core into a supermassive black hole, the basis of the paper's model II."},{"cited_title":"Modelling the Track of the GD-1 Stellar Stream Inside a Host with a Fermionic Dark Matter Core-Halo Distribution","cited_arxiv_id":"2404.19102","evidence_quote":"Shows how the collapse-produced heavy seeds grow to supermassive black holes, providing the growth context for model II."},{"cited_title":"Lynden-Bell, Statistical mechanics of violent relax- ation in stellar systems, Monthly Notices of the Royal Astronomical Society, Vol","cited_arxiv_id":null,"evidence_quote":"Supplies the NFW power-law halo used as the comparison baseline for the classical spike."}],"review_version":1}