REVIEW 3 major objections 4 minor 45 references
Critical and super-critical scatterings in baryogenesis and leptogenesis
T0 review · 3 major / 4 minor · reviewed 2026-08-06 · deepseek-v4-flash
Pith's one-line read A class of scatterings with rates that grow as the universe cools can generate the observed cosmic matter excess.
desk verdict New scaling route to asymmetry generation that deserves a referee, but the concrete model has two load-bearing gaps: the s/u topology proof and the Vµ IR cutoff. 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 the interaction topology forced by $U(1)_{B-L}$ charge assignments: $\rho$ carries charge $-2$, and $X_1$ and $X_2$ carry charge $-1$. For a non-relativistic $X$ and a relativistic $\rho$, each of the $s$- and $u$-channel fermion-exchange diagrams in $X_i \rho \to X_j \rho$ contributes $\delta M \propto 1/k_{\mathrm{CM}}$ at leading order, but if both channels are present their leading terms cancel destructively, leaving a momentum-independent amplitude. The charge assignment admits only one channel per flavour-changing process, so no cancellation occurs and the cross-section scales as $1/k_{\mathrm{CM}}^2$, i.e. as $1/T^2$ in the bath; with $n_\rho \propto T^3$ the rate $R \propto T$, while $H \propto T^2/M_{\mathrm{Pl}}$. This 'super-critical' scaling is what lets the scattering rate outlive the Hubble rate. The companion processes $X_i V_\mu \leftrightarrow \bar X_j \rho$ have one gauge vertex, which is momentum suppressed, giving $\sigma \propto 1/k_{\mathrm{CM}}$ and a critical rate $\propto H$; these together with $X_2 \to \bar X_1 \rho$ deplete the heavier species and freeze out the asymmetry.
What would settle it
Compute the tree-level amplitude for $X_i \rho \to X_j \rho$ with generic $U(1)$ charges, allowing both $s$- and $u$-channel fermion exchanges: if the $1/k_{\mathrm{CM}}$ terms cancel and the thermally averaged cross-section stops growing as $1/T^2$, the claimed super-critical rate does not occur. A second check is to vary parameters until the condition of Eq. (6) is violated; the resulting O(1) thermal mixing would rotate the mass eigenstates and quench the asymmetry below the observed baryon-to-entropy ratio.
Extended reading notes
Core claim
The central claim is that particle-number- or CP-violating scatterings with thermally averaged cross-sections scaling as $T^n$ with $n \le -1$ can drive asymmetry generation more efficiently than decays, because their rates track or overtake the Hubble rate at late times. For $n=-1$ the rate is Hubble-like ('critical'); for $n<-1$ it inevitably exceeds Hubble ('super-critical'), regardless of coupling strength. The paper realizes this in a two-flavour $U(1)_{B-L}$ model with heavy Dirac fermions $X_1, X_2$ at roughly 10 PeV and a complex scalar $\rho$: the charge assignments leave only one of the $s$- and $u$-channel fermion-exchange diagrams available for $X_i \rho \to X_j \rho$, so the leading $1/k_{\mathrm{CM}}$ terms do not cancel and $\sigma \propto 1/k_{\mathrm{CM}}^2$, or $1/T^2$ after thermal averaging. The CP-violating process $X_2 \rho \to X_1 \rho$ seeds a flavour asymmetry that becomes a global $B-L$ asymmetry, with freeze-out set by depletion of $X_2$; the benchmark numerical solution yields $|\Delta_\rho| \approx |\varepsilon| Y_1^{\mathrm{FO}}$, close to maximal efficiency, and matches the observed $Y_B$. Unlike super-critical dark-matter annihilation, the relativistic-target kinematics requires only $\gamma \ge 1$, consistent with unitarity at low momentum.
Load-bearing premise
The mechanism rests on the charge assignments under $U(1)_{B-L}$ leaving only one of the $s$- and $u$-channel fermion-exchange diagrams, so the leading $1/k_{\mathrm{CM}}$ terms in the amplitude do not cancel; if both channels contribute, or if thermal mixing restores the interfering topology, the cross-section becomes momentum-independent and the super-critical growth is lost.
Editorial extensions
If this is right
- Asymmetry generation can remain efficient down to scales well below the heavy fermion mass, because the scattering rate grows rather than decays as temperature drops.
- The required CP-violating phases and $B-L$-violating couplings can be much smaller than in ordinary decay-based leptogenesis, since the super-critical rate overtakes Hubble independently of the coupling strength.
- Freeze-out is fixed not by the expansion rate but by the depletion of the heavier fermion population through decays and critical scatterings.
- Models in which the CP-violating scalar mediator carries a gauge charge and participates in symmetry breaking should generically exhibit the same super-critical dynamics.
- The same mechanism can serve asymmetric dark matter and low-scale leptogenesis scenarios, where only one relativistic target species is needed.
Reading between the lines
- Beyond the paper: a purely group-theoretic diagnostic could classify $U(1)$ charge assignments by whether the $s$- and $u$-channel topologies interfere, turning the mechanism into a model-selection criterion.
- Beyond the paper: if the mechanism is generic, low-scale searches for $B-L$ gauge bosons and CP-odd scalars become telling probes, because the couplings need not be small enough to suppress the scattering rate.
- Beyond the paper: the same rate-scaling logic suggests that scatterings with $n<-1$ could dominate other rare processes at late times, so bounds on washout should be re-derived in models with scalar mediators carrying gauge charge.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper identifies a class of early-universe scattering processes in which the rate per non-relativistic particle can grow relative to the Hubble rate as the universe cools. Parametrizing the thermally averaged cross-section as <σv> ∝ T^n, the authors call n=-1 critical and n<-1 super-critical. The general argument is clean: a non-relativistic particle scattering off a relativistic bath with σ ∝ 1/k_CM^2 has Γ = n_T<σv> ∝ T, while H ∝ T^2. As a concrete realization, they construct a U(1)_{B-L} model with two heavy Dirac fermions X1,X2 (charge -1) and a complex scalar ρ (charge -2), in which CP-violating X2ρ→X1ρ scattering is claimed to have σ ∝ 1/k_CM^2, producing a rapidly growing asymmetry that freezes out when X2 is depleted. A Boltzmann calculation with a benchmark parameter set yields an asymmetry near the observed value. The paper also addresses thermal corrections to the X mass matrix and the IR-divergent t-channel V_μ exchange.
Significance. The general scaling argument is simple, correct, and potentially important: it shows that scattering processes need not become inefficient in an expanding universe if the target is relativistic and the cross-section grows fast enough, and it correctly notes that this evades the unitarity obstruction that applies to non-relativistic DM annihilation. If the model's topology claim holds, the mechanism offers a new and efficient route to baryogenesis/leptogenesis with possibly smaller couplings and lower scales. The manuscript is largely self-contained, with explicit appendices for decay rates, cross-sections, thermal mass corrections, and the Boltzmann decomposition. The main caveats are that the crucial topology assertion is not demonstrated and the numerical illustration is a benchmark fit rather than a predictive scan.
major comments (3)
- [Figs. 1-2 and the text between them] The entire super-critical behavior in the model rests on the statement that the U(1)_{B-L} charges of X_i (-1) and ρ (-2) leave only one of the s- and u-channel fermion-exchange diagrams for X2ρ→X1ρ, so that the two leading 1/k_CM contributions do not cancel and σ ∝ 1/k_CM^2. This claim is asserted but not demonstrated. Table III quotes cross-sections only in the Δm_X→0 limit, which obscures the propagator structure. Please provide an explicit charge-flow or amplitude-level demonstration that exactly one diagram contributes, and in particular rule out an antiparticle-exchange (crossed) diagram; if both diagrams contribute, the leading 1/k_CM terms cancel and the super-critical scaling disappears.
- [Appendix B, Eqs. (B5)-(B6)] The t-channel V_μ exchange divergences are regularized by introducing an angular cutoff θ_min and retaining only the finite parts (I1≈-2ln2, I2≈-1/3). No physical regulator (Debye screening, thermal width) is introduced, and no proof of scheme independence is given. Since the CP asymmetry coefficients in Table III contain terms proportional to g_{B-L}^2 ln2, the size of ε and hence the final asymmetry shown in Fig. 3 may depend on this ad hoc prescription. Please replace the cutoff by a physical regulator or demonstrate that the quoted results are regulator-independent.
- [Results, Table I and Fig. 3] The benchmark parameters in Table I are chosen so that the final asymmetry matches the observed value; the numerical agreement is therefore an illustration, not a prediction. The abstract and discussion claim that the mechanism is efficient 'even at low scales and with small couplings,' but the benchmark uses m1=30 PeV and no parameter scan is presented. A scan or at least a sensitivity estimate in (m1, Δm_X, α_ij) is needed to support the low-scale claim and to show that the freeze-out behavior in Fig. 3 is generic rather than selected.
minor comments (4)
- [Table III, definition of Y] The definition of Y just above Table III is garbled by the typesetting ('Y≡16· |y11y∗12+ y12y∗22|2'); it should read |y11 y12^* + y12 y22^*|^2 (or whichever combination is intended) so that the cross-sections can be checked.
- [Eq. (11)] Equation (11) quotes d ln ΔX/d ln x ≈ -α11 α_B-L M_Pl/m1 without a numerical coefficient or derivation; please provide the numerical factor or a reference.
- [Eq. (12)] The statement in Eq. (12) that |Δρ| ≃ |ε| Y_1^FO is presented as a numerical observation; it would be useful to state whether this is a fitted relation or a derived consequence of the Boltzmann system.
- [Discussion] There are several typographical artifacts in the text, including 'raison d'ˆetre' in the Discussion and the typesetting of T_{B-L} in the cosmology section; a careful proofreading pass would improve readability.
Circularity Check
No significant circularity: the super-critical scaling is derived from the model and its cross-sections; the benchmark is an illustration, not a fitted prediction.
full rationale
The paper's central derivation chain is self-contained: Eq. (3) fixes the interactions, the U(1)_{B-L} charge assignments are used to argue that only one fermion-exchange channel contributes to X2 rho -> X1 rho, Table III lists the resulting sigma proportional to 1/k_CM^2, and Eqs. (7)-(10) convert this into scattering rates and asymmetry evolution. None of these steps redefines an output as an input or fits a parameter and then calls it a prediction. Table I is explicitly a benchmark ('which yield an asymmetry near the observed value'), not a fitted prediction, and the observed Y_B is not claimed as a derived prediction. The load-bearing topology assertion, given in the paragraph beginning 'However, the interaction topology dictated by...' around Figs. 1 and 2, is not displayed in detail, and the non-perturbative CP contribution is deferred to companion paper [24]; these are completeness or support issues, not circular reductions. Self-citations [19, 20] are used only to contrast with DM annihilation and do not carry the derivation. Therefore no significant circularity is identified.
Assumptions & free parameters
free parameters (9)
- m1 =
30 PeV
- Delta m_X / m_X =
0.2
- alpha_B-L =
10^-5
- alpha_11 =
10^-4
- alpha_22 =
10^-2
- alpha_12 =
10^-2
- Theta (CP phase) =
pi/4
- w_ij =
0
- T_B-L =
~10 TeV (indicative)
assumptions (6)
- standard math The universe is radiation-dominated with H proportional to T^2/M_Pl during the relevant epoch.
- domain assumption The Sakharov conditions hold and sphalerons convert the lepton asymmetry to a baryon asymmetry.
- domain assumption Rho remains in chemical equilibrium with the plasma via rapid rho rho* annihilations.
- domain assumption The X fermions are highly non-relativistic and their mass eigenstates are approximately the zero-temperature states, with thermal mixing negligible under Eq. (6).
- ad hoc to paper Infrared-divergent angular integrals from t-channel V_mu exchange are regularized by a cutoff angle and only the finite part is retained.
- domain assumption The B-L phase transition at T ~ 10 TeV does not interfere with asymmetry generation at T ~ a few PeV.
invented entities (4)
-
U(1)B-L gauge boson V_mu
independent evidence
-
Complex scalar rho with B-L charge -2
-
Heavy Dirac fermions X1 and X2
-
Anomaly-canceling fermions chi
independent evidence
Cite this review
Pith. "Pith review of Critical and super-critical scatterings in baryogenesis and leptogenesis." pith.science (2026). https://pith.science/paper/A62ZQQF5
@misc{pith2026250711532,
author = {Pith},
title = {Pith review of: Critical and super-critical scatterings in baryogenesis and leptogenesis},
year = {2026},
howpublished = {\url{https://pith.science/paper/A62ZQQF5}},
note = {Machine review of arXiv:2507.11532}
}
read the original abstract
In many theories, matter-antimatter asymmetries originate from out-of-equilibrium decays and scatterings of heavy particles. While decays remain efficient, scattering rates typically drop below the Hubble rate as the universe expands. We point out the possibility of scatterings between non-relativistic particles and the relativistic bath whose cross-sections grow with decreasing temperature, leading to scattering rates that track or exceed the Hubble rate at late times. This results in soaring asymmetry generation, even at low scales and with small CP- or baryon/lepton-violating couplings.
Figures
Reference graph
Works this paper leans on
-
[1]
N. Aghanim et al. (Planck), Astron. Astrophys. 641, A6 (2020), [Erratum: Astron.Astrophys. 652, C4 (2021)], arXiv:1807.06209 [astro-ph.CO]
arXiv 2020
-
[2]
R. J. Cooke, M. Pettini, and C. C. Steidel, Astrophys. J. 855, 102 (2018), arXiv:1710.11129 [astro-ph.CO]
arXiv 2018
-
[3]
Nucleosynthesis Predictions and High-Precision Deuterium Measurements
S. Riemer-Sørensen and E. S. Jenssen, Universe 3, 44 6 (2017), arXiv:1705.03653 [astro-ph.CO]
work page Pith review arXiv 2017
-
[4]
A. D. Sakharov, Pisma Zh. Eksp. Teor. Fiz. 5, 32 (1967)
1967
-
[5]
K. Petraki and R. R. Volkas, Int. J. Mod. Phys. A 28, 1330028 (2013), arXiv:1305.4939 [hep-ph]
arXiv 2013
-
[6]
K. M. Zurek, Phys. Rept. 537, 91 (2014), arXiv:1308.0338 [hep-ph]
arXiv 2014
-
[7]
Weinberg, Phys
S. Weinberg, Phys. Rev. Lett. 42, 850 (1979)
1979
-
[8]
J. N. Fry, K. A. Olive, and M. S. Turner, Phys. Rev. Lett. 45, 2074 (1980)
work page 1980
Show all 45 references
-
[9]
Fukugita and T
M. Fukugita and T. Yanagida, Phys. Lett. B 174, 45 (1986)
1986
-
[10]
Yoshimura, Phys
M. Yoshimura, Phys. Rev. Lett. 41, 281 (1978), [Erra- tum: Phys.Rev.Lett. 42, 746 (1979)]
1978
-
[11]
E. W. Kolb and S. Wolfram, Nucl. Phys. B 172, 224 (1980), [Erratum: Nucl.Phys.B 195, 542 (1982)]
1980
-
[12]
Bento and Z
L. Bento and Z. Berezhiani, Phys. Rev. Lett. 87, 231304 (2001), arXiv:hep-ph/0107281
2001 arXiv
-
[13]
Buchmuller, R
W. Buchmuller, R. D. Peccei, and T. Yanagida, Ann. Rev. Nucl. Part. Sci. 55, 311 (2005), arXiv:hep- ph/0502169
2005
-
[14]
Davidson, E
S. Davidson, E. Nardi, and Y. Nir, Phys. Rept. 466, 105 (2008), arXiv:0802.2962 [hep-ph]
2008 arXiv
-
[15]
Langacker, Phys
P. Langacker, Phys. Rept. 72, 185 (1981)
1981
-
[16]
Falkowski, J
A. Falkowski, J. T. Ruderman, and T. Volansky, JHEP 05, 106 (2011), arXiv:1101.4936 [hep-ph]
2011 arXiv
- [17]
-
[18]
Binder, M
T. Binder, M. Garny, J. Heisig, S. Lederer, and K. Ur- ban, Phys. Rev. D 108, 095030 (2023), arXiv:2308.01336 [hep-ph]
2023 arXiv
-
[19]
M. M. Flores and K. Petraki, Phys. Lett. B 858, 139022 (2024), arXiv:2405.02222 [hep-ph]
2024 arXiv
-
[20]
Petraki, A
K. Petraki, A. Socha, and C. Vasilaki, (2025), arXiv:2505.20443 [hep-ph]
2025
- [21]
-
[22]
Chakrabortty, P
J. Chakrabortty, P. Konar, and T. Mondal, Phys. Rev. D 89, 095008 (2014), arXiv:1311.5666 [hep-ph]
2014 arXiv
-
[23]
V. A. Rubakov and M. E. Shaposhnikov, Usp. Fiz. Nauk 166, 493 (1996), arXiv:hep-ph/9603208
1996 arXiv
-
[24]
Flores, K
M. Flores, K. Petraki, and A. Socha, In preparation (2025)
2025
-
[25]
R. E. Cutkosky, J. Math. Phys. 1, 429 (1960)
1960
-
[26]
Gondolo and G
P. Gondolo and G. Gelmini, Nucl. Phys. B 360, 145 (1991)
1991
-
[27]
B´ elanger, F
G. B´ elanger, F. Boudjema, A. Goudelis, A. Pukhov, and B. Zaldivar, Comput. Phys. Commun. 231, 173 (2018), arXiv:1801.03509 [hep-ph]
2018 arXiv
- [28]
-
[29]
Carena, A
M. Carena, A. Daleo, B. A. Dobrescu, and T. M. P. Tait, Phys. Rev. D 70, 093009 (2004), arXiv:hep-ph/0408098
2004 arXiv
-
[30]
Cacciapaglia, C
G. Cacciapaglia, C. Csaki, G. Marandella, and A. Stru- mia, Phys. Rev. D 74, 033011 (2006), arXiv:hep- ph/0604111
2006
-
[31]
R. N. Mohapatra, Phys. Rev. Lett. 56, 561 (1986)
1986
-
[32]
R. N. Mohapatra and J. W. F. Valle, Phys. Rev. D 34, 1642 (1986)
1986
-
[33]
H. A. Weldon, Phys. Rev. D 26, 2789 (1982)
1982
-
[34]
E. J. Levinson and D. H. Boal, Phys. Rev. D 31, 3280 (1985)
1985
-
[35]
M. L. Bellac, Thermal Field Theory , Cambridge Mono- graphs on Mathematical Physics (Cambridge University Press, 2011)
2011
-
[36]
M. G. Mustafa, Eur. Phys. J. ST 232, 1369 (2023), arXiv:2207.00534 [hep-ph]. 7 APPENDICES Appendix A: Cosmological history
2023 arXiv
-
[37]
Phase transition and ϱ decays The B − L symmetry is restored at high temperatures by thermal corrections to the potential (see e.g., [28]), which give ρ the effective mass, m2 ρ(T ) ≈ g2 B−L + λρ 12 + λρH 6 − 3X k=1 |yχ k |2 4 ! T 2 − µ2 ρ, (A1) where we have included contribu...
-
[38]
The global charge carried by X1 must cascade into SM leptons before the EW phase transition, such that it is reprocessed into a baryon num- ber by sphalerons
Cascade of the asymmetry The X1 fermions decay into SM leptons with rate ΓX1→LlH ≃ |λX 1l|2m1/(16π). The global charge carried by X1 must cascade into SM leptons before the EW phase transition, such that it is reprocessed into a baryon num- ber by sphalerons. However, after th...
-
[39]
Type A:The χ fermions may couple to the ρ boson via Lρ,χ = − 3X k=1 yχ k 2 ρ†χc kχk + h.c., (A3) where we have diagonalized the yχ matrix
Stable relics In both Type A and B scenarios, summarized in Ta- ble II, the anomaly-canceling fermions χ are stable relics, but with distinct cosmological roles. Type A:The χ fermions may couple to the ρ boson via Lρ,χ = − 3X k=1 yχ k 2 ρ†χc kχk + h.c., (A3) where we have diag...
-
[40]
Neutrino masses The model can also account for the active neutrino masses via an X-induced inverse seesaw mechanism [31, 32]. After B − L and EW symmetry breaking, the one- generation mass mixing matrix, in the (νL, XL, Xc R)T ba- 8 sis, is MF = 0 0 λXvEW/ √ 2 0 yLvB−L/ √ ...
-
[41]
(A5) This result holds both in Type A and B scenarios, in neither of which do theχ fermions mix with νL and XL,R
The third eigenstate, consisting mainly of νL, has mass mν ≈ 1 2 √ 2 yL|λ X |2 vB−L v2 EW m2 X . (A5) This result holds both in Type A and B scenarios, in neither of which do theχ fermions mix with νL and XL,R. Appendix B: Interactions and rates
-
[42]
Decays The X2 fermion, being the heaviest field, decays into ¯X1 and ρ. For ∆ mX ≪ mX, ρ is emitted with energy Eρ ≃ ∆mX, and the corresponding decay rate is Γ2→ρ¯1 ≃ 2α12 ∆mX " 1 − mρ ∆mX 2#1/2 , (B1) where the temperature-dependent mass mρ is found in Eq. (A1). The thermally...
-
[43]
We consider s-wave contribu- tions only
Scatterings The tree-level cross-sections, σtree, for 2-to-2 processes — excluding the interactions Xi ¯Xj ↔ Xi′ ¯Xj′, which will be discussed in a companion paper [24] — are listed in Table III and grouped according to the totalB −L charge of the interacting states. We consid...
-
[44]
4, is given by Σ(p) = Z d4k (2π)4 i k2 i (yL∗ jk PR + yR∗ jk PL) × i (/p − /k + mj) (p − k)2 − m2 j i (yL ijPL + yR ijPR)
Self-energy of the X field The one-loop contribution to the self-energy of the X field, shown in Fig. 4, is given by Σ(p) = Z d4k (2π)4 i k2 i (yL∗ jk PR + yR∗ jk PL) × i (/p − /k + mj) (p − k)2 − m2 j i (yL ijPL + yR ijPR). (C7) In the imaginary time formalism, Σ( p) can be e...
-
[45]
γ + ¯γ 2 1 + Y eq 2 Y eq 1 − β + ¯β 2 1 − Y eq 2 Y eq 1 # Y eq ρ + (ϵ2ρ→1ρ + ϵ¯2ρ→¯1ρ)
Non-relativistic approximation We seek to find the X mass matrix at a period relevant for the asymmetry generation. To that end, we compute the structure functions in the non-relativistic limit p0 = q m2 i + |p|2 ≫ |p| ∼ p miT ≫ k0, |k| ∼T. (C14) Neglecting the mass difference...
Reviewed August 6, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.