REVIEW 4 major objections 6 minor 1 cited by
Constraining decaying dark matter models with gravitational lensing and cosmic voids
T0 review · 4 major / 6 minor · reviewed 2026-08-06 · deepseek-v4-flash
Pith's one-line read Decaying dark matter should leave a measurable imprint in the weak-lensing signal of cosmic voids.
desk verdict A useful exploratory prediction for void lensing as a DDM probe, but the signature rests on an ad hoc interaction term that needs microphysical justification before the forecast means much. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The machinery is a covariant two-fluid description of dark matter in a locally rotationally symmetric spacetime. One dust fluid (density $\rho$, 4-velocity $u^a$) represents the mother particles; the other (density $\lambda$, 4-velocity $v^a$) represents the decay products. Decay is encoded in the interaction term $I^a_{(1)} = -\Gamma\rho\,(u^a + w^a)$ with $w^a = (4/\pi^2)\, v_i\, \delta\, z^a$, where $\Gamma$ is the decay rate, $v_i$ is the injection velocity, and $z^a$ is the preferred spatial direction; this term produces an effective acceleration $A = -(4/\pi^2)\Gamma\delta\sqrt{2\epsilon}$ that drives daughter particles into voids. The evolution equations (5)--(11) are integrated from Lemaître--Tolman initial conditions set at redshift 150, and the resulting density profiles are projected onto the lens plane to compute convergence $\kappa = \Sigma/\Sigma_c$ and shear $\gamma = (\bar{\Sigma}-\Sigma)/\Sigma_c$.
What would settle it
Stack the weak-lensing convergence and shear of several thousand voids at lens redshift $z\approx0.25$ with sources near $z\approx0.5$ and compare the $\kappa$--$\gamma$ distribution to the no-decay prediction at better than 2% precision. If the observed distribution matches the no-decay case and shows no deficit of high-convergence voids in the 10--40 Mpc radius range, the DDM scenario with $\Gamma\gtrsim0.5H_0$ and injection velocities of 25--75 km s$^{-1}$ is ruled out; a measured deficit that scales with void size and lifetime would confirm the model.
Extended reading notes
Core claim
The central claim is a new observational route to dark-matter microphysics: the decay rate $\Gamma$ and the mass-splitting parameter $\epsilon$ (equivalently the injection velocity $v_i$) leave a distinguishable imprint on the gravitational lensing produced by cosmic voids. In the model, decay transfers particles from the comoving mother dust into a daughter fluid moving with velocity $v_i$; because matter density is higher near void edges, more decay products flood inward, reducing the central density contrast and the edge gradient. The result is that the maximal weak-lensing convergence $\kappa$ and shear $\gamma$ are suppressed relative to the no-decay case. The suppression is not visible for individual voids but appears statistically, and it survives Gaussian noise only at roughly 2% measurement precision; at 5--10% the signal is washed out.
Load-bearing premise
The entire predicted signal rests on the ad hoc interaction term of Eq. (2), which assumes that dark-matter decay produces a coherent bulk velocity field directed along the density gradient; if actual decays do not generate such a velocity field, the flooding of voids and the associated lensing deficit do not occur.
Editorial extensions
If this is right
- The void size function is not a useful probe: changes in injection velocity produce no statistically significant shift in the void number density as a function of radius.
- Weak lensing is the sensitive observable: at fixed shear, decaying dark matter creates an excess of low-convergence voids and a deficit of high-convergence voids, with the effect stronger for smaller voids and shorter lifetimes.
- Individual voids cannot be diagnosed one by one; the signature is statistical and requires stacking many voids, with a measurement precision of about 2% in convergence.
- The nonlinear treatment reaches injection velocities as low as 75 km s$^{-1}$ and lifetimes of order the Hubble time, extending the parameter space accessible relative to earlier linear modelling that required 90 km s$^{-1}$ and 5 Gyr lifetimes.
Reading between the lines
- If the interaction term is replaced by a particle-physics-derived decay kernel, the effective velocity $w^a$ may not align with the void density gradient; the flooding signature would then weaken, so the paper's detectability forecast is an upper bound on what generic DDM models can do.
- Because the signal enters through the combination $\Gamma\sqrt{\epsilon}$, void lensing alone cannot separate the decay rate from the mass splitting; combining it with Doppler magnification or void-galaxy correlations could break this degeneracy.
- A direct testable prediction: in future deep surveys, the stacked void lensing signal should show a deficit of high-convergence voids that grows as void radius shrinks from 30--40 Mpc to 10--20 Mpc, and that strengthens as the assumed lifetime drops from $2H_0^{-1}$ to $0.5 H_0^{-1}$.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies a phenomenological two-fluid model of decaying dark matter (DDM) in which a mother dust fluid decays into a massive daughter fluid with a nonzero injection velocity. The model is applied to spherically symmetric cosmic voids using covariant LRS evolution equations taken from the authors' earlier work. The authors first calibrate the initial void mass profile so that the no-decay LTB void reproduces the empirical HSW void profile, then evolve the same initial conditions with and without decay and compute the weak-lensing convergence and shear signals. They find that DDM 'floods' voids, making them shallower and reducing the lensing amplitude, and they conclude that a statistical detection via stacked void weak lensing would require percent-level precision, possibly achievable with LSST. The paper also shows that the void size function is not sensitive to the considered decay parameters.
Significance. If the predicted 'flooding' signature were robust, this would constitute a genuinely new astrophysical probe of dark matter decay, complementary to CMB and large-scale-structure constraints, and the paper is commendable for targeting cosmic voids, where baryonic contamination is minimized. The authors are explicit that the interaction term is phenomenological and they provide forward-modeled lensing distributions rather than a full survey simulator. The main value is therefore as a proof-of-concept: it identifies a possible observable signature and quantifies, in a preliminary way, the precision needed to see it. However, the significance is limited by the ad hoc nature of the decay interaction, the calibration of initial conditions to the HSW profile, and the absence of a formal statistical detection calculation. The manuscript is honest about many of these limitations, but the central claim as currently stated goes beyond what the analysis can support.
major comments (4)
- [Sec. II B, Eqs. (2)-(3)] The spatial part of the interaction term, w^a = (4/pi^2) v_i delta z^a, is 'adopted' rather than derived from a particle decay model. For an isotropic two-body decay in a spherically symmetric environment, the mean daughter velocity at a given point should be controlled by the gradient of the parent density (a net flux from overdense to underdense regions), not by the local density contrast itself. A term proportional to delta is maximum at the void center and vanishes at the edge, whereas a diffusion-type flux would vanish at the center and peak at the boundary. Because this term is the sole agent producing the 'flooding' effect and hence the predicted lensing deficit, the central observable claim is conditional on an unvalidated phenomenological input. Please either derive the form from a concrete decay kinematics model, or explicitly treat the model as a toy and test the sensitivity of the lensing predictions to alternative spatial dependences (for example, a term proportional to the density gradient).
- [Sec. II H and Fig. 4] The initial void parameters m0, r0, and Delta r are fine-tuned so that the no-decay LTB void reproduces the empirical HSW profile. This is a calibration step, not a derivation from cosmological initial conditions, and it means that the DDM prediction is a relative modification of an empirically matched profile rather than an absolute prediction from a power spectrum. The procedure is not circular in the sense that the DDM parameters are not fit to the target lensing observable, but the resulting lensing distributions inherit the HSW assumptions. The paper should state this explicitly and quantify how the DDM lensing signal changes when the calibration parameters are varied within the ranges quoted in Sec. II D.
- [Sec. III, Fig. 8] The claim that detection 'requires measurements of weak lensing signal with a precision of 2%' is inferred from a single noise-injection exercise on one panel (the lower-left panel of Fig. 5) with Gaussian scatter at 2%, 5%, and 10%. No formal significance calculation is presented: there is no test statistic, no number of voids, no survey volume or area, and no detection threshold. The visual separation between the DDM and no-DDM distributions in Figs. 5-7 is also not quantified. Please provide a proper detectability estimate that specifies the survey parameters, the number of stacked voids, and a statistical test (e.g., a Kolmogorov-Smirnov or Anderson-Darling test on the convergence-shear distribution) with a stated significance level.
- [Sec. II C] The evolution equations (5)-(10) are carried over from the authors' previous paper [26] without derivation or independent validation in this manuscript. Since the new physics enters precisely through the interaction term in Eq. (2) and the resulting acceleration A in Eq. (4), the reader cannot verify from this paper alone that the equations are consistent or that the numerical implementation conserves total energy-momentum. At minimum, please summarize the derivation of the system, state the assumptions under which it is valid, and provide a convergence test or a comparison with an independent method for at least one representative void.
minor comments (6)
- [Abstract] There is a typo: 'phenomological' should be 'phenomenological'.
- [Sec. II F] The section heading 'W eak gravitational lensing' contains an erroneous space; it should be 'Weak gravitational lensing'.
- [Sec. II F, Eq. (17)] The quantity \bar{\Sigma} is called the differential surface mass density but is not defined before its first use; please define it explicitly and distinguish it from the shear scalar \Sigma used earlier in Sec. II C.
- [Sec. II D and Sec. II H] The notation for the injection velocity is inconsistent: the text uses both v_in and v_i. Please use a single symbol throughout.
- [Sec. II H] The void redshift and source redshift are drawn from normal distributions, but the number of realizations used to produce Figs. 5-7 is not stated. Please specify the sample size and whether the same realizations are used for the DDM and no-DDM cases.
- [Fig. 3 caption] The caption says 'parameters increasing size and decay rate (from left to right)' but does not list the actual values; please state the void and source redshifts and the decay parameters for each panel.
Circularity Check
No significant circularity: the DDM parameters are scanned inputs, the lensing predictions are compared against a no-decay model with identical initial conditions, and the cited prior work is a derivation pointer rather than a load-bearing uniqueness claim.
full rationale
The paper's central claim is that decaying dark matter with massive daughters shallows cosmic voids and reduces weak-lensing convergence and shear. This is a derived consequence of the adopted interaction term in Eqs. (2)-(3), not a re-statement of the target observable. The decay parameters Gamma and v_i are scanned over a grid, and the resulting lensing distributions are compared with the no-decay LTB case evolved from the same randomized initial conditions. No parameter is fitted to the lensing signal it is then said to predict. The initial void parameters m0, Delta r, and r0 are calibrated to reproduce the empirical HSW void profile, which is an external N-body benchmark; this calibration affects the baseline no-decay model but does not force the DDM deficit, which is a model-dependent deviation from that baseline. The paper's self-citation of [26] is used to avoid re-deriving the propagation system (5)-(11), but the equations are displayed and the interaction term is introduced in this paper, so the central lensing claim does not reduce to the cited prior work. The main vulnerability is that the spatial interaction term w^a = (4/pi^2) v_i delta z^a is an ad hoc phenomenological choice, not derived from a particle model; if realistic DDM decay produces isotropic daughter velocities rather than a coherent radial bulk flow, the flooding effect and the predicted few-percent detectability would not follow. That is a correctness or robustness risk, not a circularity: the prediction is not equivalent to its input by construction. Overall the derivation chain is self-contained against an external benchmark, so no circular step is exhibited.
Assumptions & free parameters
free parameters (4)
- Dark matter decay rate Γ =
0.5 H0, 1.0 H0, 2.0 H0
- Injection velocity vi (mass splitting ϵ) =
25, 75, 100, 135 km/s (ϵ up to about 10^-7)
- Initial void mass profile parameters m0, r0, Δr =
m0 in [-0.0225,-0.015], r0 in [45,200] kpc, Δr in (0.425-0.5) r0
- HSW profile parameters α, β, δc, rs =
Empirically calibrated distributions from Hamaus et al. 2014, Fig. 2
assumptions (7)
- domain assumption Two dark fluids are perfect dust in their own frames and permeate each other without interfaces (T^ab_(1)=ρu^a u^b, T^ab_(2)=λv^a v^b).
- domain assumption The daughter particle is non-relativistic (γ≈1) and the radiative decay component is negligible.
- ad hoc to paper The interaction term has the form I^a_(1) = -Γρ(u^a + w^a) with w^a = (4/π^2)v_i δ z^a, producing the effective acceleration A = -(4/π^2)Γδ√(2ϵ).
- domain assumption The evolution equations (5)-(11) and the LRS covariant decomposition are taken from the authors' previous paper [26].
- domain assumption Voids are isolated, spherically symmetric, with uniform bang time tB=0 and mass profile Eq. (12) initialized at zi=150.
- domain assumption Linearized weak lensing with convergence κ=Σ/ΣC and shear γ=(Σ̄-Σ)/ΣC applies, and Doppler magnification is neglected.
- domain assumption The HSW profile, Eq. (18), is a universal description of cosmic void density profiles.
invented entities (1)
-
Massive daughter dark matter fluid with injection velocity V^a relative to the mother fluid
independent evidence
Cite this review
Pith. "Pith review of Constraining decaying dark matter models with gravitational lensing and cosmic voids." pith.science (2026). https://pith.science/paper/CNK2KLU7
@misc{pith2026250708275,
author = {Pith},
title = {Pith review of: Constraining decaying dark matter models with gravitational lensing and cosmic voids},
year = {2026},
howpublished = {\url{https://pith.science/paper/CNK2KLU7}},
note = {Machine review of arXiv:2507.08275}
}
read the original abstract
Despite overwhelming observational evidence for dark matter, we still have no evidence of direct detection. Consequently, our knowledge about dark matter is limited, for example, we do not know if dark matter is a stable particle or if it decays. Without a theoretical particle model, the parameter space of possible decay models is highly variable, and astrophysical or cosmological means of indirectly constraining the phenomenological models are required. This paper investigates a scenario in which a dark matter decays and the dark daughter particle moves with respect to the comoving mother particle. The model is parameterised by the decay rate and the injection velocity of the dark matter particles, which can be converted to the mass ratio. In previous work, a simpler model was used to investigate the evolution of cosmic voids typified as regions with low content of galaxies and non-baryonic matter. It was found that the growth of S-type voids is modified by the dark matter decay, leading to imprints at the present day. Here we extend our study and improve the method used to model the decay. We also study the gravitational weak-lensing signal that will be able to detect or constrain the parameter space of decaying dark matter. The results of this study suggest that future weak lensing surveys may provide unique probes of the phenomological parameters of dark matter.
Figures
Figures from the paper (5 more)
Forward citations
Cited by 1 Pith paper
-
Why Cosmic Voids Matter: Pristine Evolution
Cosmic voids traced by halos become stable at late times, and the matter around them evolves linearly, supporting their use as clean dark-energy probes.
Reference graph
Works this paper leans on
- [26]
-
[1]
P. Bull, Y. Akrami, J. Adamek, T. Baker, E. Bellini, J. B. Jimenez, E. Bentivegna, S. Camera, S. Clesse, J. H. Davis, et al., Physics of the Dark Universe 12, 56 (2016)
work page 2016
-
[2]
E. Di Valentino, O. Mena, S. Pan, L. Visinelli, W. Yang, A. Melchiorri, D. F. Mota, A. G. Riess, and J. Silk, Clas- sical and Quantum Gravity 38, 153001 (2021). 10 (a) 10 < Rv < 20 Mpc, vin = 25 kms−1 (b) 20 < Rv < 30 Mpc, vin = 25 kms−1 (c) 30 < Rv < 40 Mpc, vin = 25 kms−1 (d) 10 < Rv < 20 Mpc, vin = 75 kms−1 (e) 20 < Rv < 30 Mpc, vin = 75 kms−1 (f) 30 <...
work page 2021
-
[3]
Cen, The Astrophysical Journal 546, L77 (2001)
R. Cen, The Astrophysical Journal 546, L77 (2001)
work page 2001
-
[4]
K. L. Pandey, T. Karwal, and S. Das, Journal of Cosmol- ogy and Astroparticle Physics 2020, 026 (2020)
work page 2020
- [5]
-
[6]
K. Enqvist, S. Nadathur, T. Sekiguchi, and T. Takahashi, Journal of Cosmology and Astroparticle Physics 2015, 067 (2015)
work page 2015
- [7]
Show all 53 references
-
[8]
Doroshkevich and M
A. Doroshkevich and M. Y. Khlopov, Soviet Astronomy Letters 11, 236 (1985)
1985
-
[9]
M.-Y. Wang, A. H. Peter, L. E. Strigari, A. R. Zentner, B. Arant, S. Garrison-Kimmel, and M. Rocha, Monthly Notices of the Royal Astronomical Society 445, 614 (2014)
2014
-
[10]
Boehm, A
C. Boehm, A. Riazuelo, S. H. Hansen, and R. Schaeffer, Physical Review D 66, 083505 (2002)
2002
-
[11]
D. An, K. A. Meissner, P. Nurowski, and R. Penrose, M.N.R.A.S. 495, 3403 (2020), ISSN 0035-8711, URL https://doi.org/10.1093/mnras/staa1343
2020 doi
- [12]
-
[13]
Simon, G
T. Simon, G. F. Abell´ an, P. Du, V. Poulin, and Y. Tsai, arXiv preprint arXiv:2203.07440 (2022)
2022 arXiv
-
[14]
M.-Y. Wang, R. A. Croft, A. H. Peter, A. R. Zentner, and C. W. Purcell, Physical Review D 88, 123515 (2013)
2013
- [15]
-
[16]
E. B. Holm, T. Tram, and S. Hannestad, Journal of Cosmology and Astroparticle Physics 2022, 044 (2022)
2022
-
[17]
Lacroix, J
T. Lacroix, J. Silk, E. Moulin, and C. Boehm, Phys. Rev. D 94, 123008 (2016), [arXiv:1603.05977]
2016 arXiv
-
[18]
H. A. Clark, P. Scott, R. Trotta, and G. F. Lewis, JCAP 2018, 060 (2018), [arXiv:1612.01539]
2018 arXiv
-
[19]
Van de Weygaert and E
R. Van de Weygaert and E. Platen, Int. J. Mod. Phys. Conf. Ser. 1, 41–66 (2011)
2011
-
[20]
L. A. Thompson and S. A. Gregory, arXiv preprint arXiv:1109.1268 (2011)
2011 arXiv
-
[21]
Peebles, The Astrophysical Journal 557, 495 (2001)
P. Peebles, The Astrophysical Journal 557, 495 (2001)
2001
-
[22]
R. K. Sheth and R. Van De Weygaert, Monthly Notices of the Royal Astronomical Society 350, 517 (2004)
2004
-
[23]
Ceccarelli, D
L. Ceccarelli, D. Paz, M. Lares, N. Padilla, and D. G. Lambas, Monthly Notices of the Royal Astronomical So- ciety 434, 1435 (2013)
2013
-
[24]
Sutter, P
P. Sutter, P. Elahi, B. Falck, J. Onions, N. Hamaus, A. Knebe, C. Srisawat, and A. Schneider, Monthly No- tices of the Royal Astronomical Society445, 1235 (2014)
2014
-
[25]
Icke, Monthly Notices of the Royal Astronomical So- ciety 206, 1P (1984)
V. Icke, Monthly Notices of the Royal Astronomical So- ciety 206, 1P (1984)
1984
-
[27]
Amendola, J
L. Amendola, J. A. Frieman, and I. Waga, Monthly No- tices of the Royal Astronomical Society 309, 465 (1999)
1999
-
[28]
Krause, T.-C
E. Krause, T.-C. Chang, O. Dore, and K. Umetsu, The Astrophysical Journal Letters 762, L20 (2012)
2012
-
[29]
M. R. Buckley and A. H. Peter, Physics Reports 761, 1 (2018)
2018
-
[30]
E. S. Lester and L. K. Forbes, Journal of Engineering Mathematics 120, 29 (2020)
2020
-
[31]
A. R. King and G. F. R. Ellis, Communications in Math- ematical Physics 31, 209 (1973)
1973
-
[32]
G. F. Ellis, Cargese lectures in physics 6, 1 (1973)
1973
-
[33]
G. F. Ellis and H. Van Elst, in Theoretical and Observa- tional Cosmology (Springer, 1999), pp. 1–116
1999
-
[34]
C. G. Tsagas, A. Challinor, and R. Maartens, Physics Reports 465, 61 (2008)
2008
-
[35]
Contarini, G
S. Contarini, G. Verza, A. Pisani, N. Hamaus, M. Sahl´ en, C. Carbone, S. Dusini, F. Marulli, L. Moscardini, A. Renzi, et al., Astronomy & Astrophysics 667, A162 (2022)
2022
-
[36]
Contarini, A
S. Contarini, A. Pisani, N. Hamaus, F. Marulli, L. Moscardini, and M. Baldi, The Astrophysical Journal 953, 46 (2023)
2023
-
[37]
Blandford and R
R. Blandford and R. Narayan, Annual review of astron- omy and astrophysics 30, 311 (1992)
1992
-
[38]
Bartelmann, Classical and Quantum Gravity 27, 233001 (2010)
M. Bartelmann, Classical and Quantum Gravity 27, 233001 (2010)
2010
-
[39]
Bonvin, Physical Review D 78, 123530 (2008)
C. Bonvin, Physical Review D 78, 123530 (2008)
2008
-
[40]
Bartelmann and P
M. Bartelmann and P. Schneider, Physics Reports 340, 291 (2001)
2001
-
[41]
Sasaki, Progress of Theoretical Physics 90, 753 (1993)
M. Sasaki, Progress of Theoretical Physics 90, 753 (1993)
1993
-
[42]
Hamaus, P
N. Hamaus, P. Sutter, and B. D. Wandelt, Physical re- view letters 112, 251302 (2014)
2014
-
[43]
M. R. Hossen, S. A. Ema, K. Bolejko, and G. F. Lewis, Monthly Notices of the Royal Astronomical Society 509, 5142 (2022)
2022
-
[44]
M. R. Hossen, S. A. Ema, K. Bolejko, and G. F. Lewis, Monthly Notices of the Royal Astronomical Society 513, 5575 (2022)
2022
-
[45]
Wang and A
M.-Y. Wang and A. R. Zentner, Physical Review D 85, 043514 (2012)
2012
-
[46]
Ivezi´ c, S
ˇZ. Ivezi´ c, S. M. Kahn, J. A. Tyson, B. Abel, E. Acosta, R. Allsman, D. Alonso, Y. AlSayyad, S. F. Anderson, J. Andrew, et al., The Astrophysical Journal 873, 111 (2019)
2019
-
[47]
Mahony, M
C. Mahony, M. C. Fortuna, B. Joachimi, A. Korn, H. Hoekstra, S. J. Schmidt, and L. D. E. S. C. A. D. S. S. R. M. H. H. D. C. J. Harry, Monthly Notices of the Royal Astronomical Society 513, 1210 (2022)
2022
-
[48]
Bertone, N
G. Bertone, N. Bozorgnia, J. S. Kim, S. Liem, C. Mc- Cabe, S. Otten, and R. Ruiz de Austri, JCAP 2018, 026 (2018), [arXiv:1712.04793]
2018 arXiv
-
[49]
G. Aad, B. Abbott, D. C. Abbott, A. Abed Abud, K. Abeling, D. K. Abhayasinghe, S. H. Abidi, O. S. AbouZeid, N. L. Abraham, H. Abramowicz, et al. (ATLAS Collaboration), Phys. Rev. Lett. 126, 121802 (2021), URL https://link.aps.org/doi/10. 1103/PhysRevLett.126.121802
2021
-
[50]
Zhang, A
D. Zhang, A. Abdukerim, Z. Bo, W. Chen, X. Chen, Y. Chen, C. Cheng, Z. Cheng, X. Cui, Y. Fan, et al. (PandaX Collaboration), Phys. Rev. Lett. 129, 161804 (2022), URL https://link.aps.org/doi/10. 1103/PhysRevLett.129.161804
2022
-
[51]
Aprile, K
E. Aprile, K. Abe, F. Agostini, S. Ahmed Maouloud, L. Althueser, B. Andrieu, E. Angelino, J. R. Angevaare, V. C. Antochi, D. Ant´ on Martin, et al. (XENON Collaboration), Phys. Rev. Lett. 129, 161805 (2022), URL https://link.aps.org/doi/10. 1103/PhysRevLett.129.161805
2022
-
[52]
Bolejko, C
K. Bolejko, C. Clarkson, R. Maartens, D. Bacon, N. Meures, and E. Beynon, Physical Review Letters 110, 021302 (2013)
2013
-
[53]
D. J. Bacon, S. Andrianomena, C. Clarkson, K. Bole- jko, and R. Maartens, M.N.R.A.S. 443, 1900 (2014), [arXiv:1401.3694]
2014 arXiv
Reviewed August 6, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.