REVIEW 3 major objections 4 minor 15 references
Higgs Thermal Nonequilibrium in Primordial QGP
T0 review · 3 major / 4 minor · reviewed 2026-08-09 · deepseek-v4-flash
Pith's one-line read In the primordial quark-gluon plasma, the Higgs boson was chemically out of equilibrium with fugacity 0.69 and kinetically cold below 25 GeV.
desk verdict The 0.69 fugacity is an artifact of omitting the reverse decay; the rate inventory is useful but the central chemical and kinetic claims are not established. 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 central object is the fugacity $\Upsilon_h$ that parametrizes chemical equilibrium in the phase-space distribution $f = 1/(\Upsilon^{-1} e^{E/T} \pm 1)$; $\Upsilon = 1$ is full abundance equilibrium, and $\Upsilon < 1$ is a deficit. The argument is carried by the population equation $d\Upsilon_h/dt = (1 - \Upsilon_h)\Gamma_{\mathrm{fusion}} - \Upsilon_h \Gamma_{h \to WW^*, ZZ^*}$, where $\Gamma_{\mathrm{fusion}}$ sums the two-particle fusion rates ($b\bar b$, $c\bar c$, $\tau\bar\tau$, $gg$ into $h$) and $\Gamma_{h \to WW^*, ZZ^*}$ is the one-way virtual-decay loss. Setting $d\Upsilon_h/dt = 0$ yields the 0.69 result. The kinetic claim is carried by the two-body scattering rate $\Gamma_{\mathrm{scattering}} = (R_{hb \to hb} + R_{ht \to ht})/n_h^{\mathrm{th}}$ computed from tree-level amplitudes; its crossing with $\Gamma_{\mathrm{fusion}}$ at $T = 25$ GeV marks the onset of a "cold" Higgs distribution.
What would settle it
Compute the full thermal rate for the inverse process $f + \bar f + W \to h$ (and its $Z$ analog) using on-mass-shell particles in the 10 to 130 GeV temperature range; if that rate times the relevant densities is not more than an order of magnitude below $\Gamma_{h \to WW^*, ZZ^*}$, the stationary solution of the population equation moves close to $\Upsilon_h = 1$, directly contradicting the paper's 0.69.
Extended reading notes
Core claim
On the paper's own terms, the discovery is a quantitative breach of detailed balance for the Higgs in the primordial QGP. Because $m_h = 124$ GeV lies below the $WW$ and $ZZ$ pair thresholds, the dominant depletion channel $h \to WW^*, ZZ^*$ produces at least one virtual gauge boson that decays with unit probability, while the reverse $3 \to 1$ process would need extra weak-interaction vertices; the paper therefore omits it from the population equation. Solving the stationary balance $(1 - \Upsilon_h)\Gamma_{\mathrm{fusion}} = \Upsilon_h \Gamma_{h \to WW^*, ZZ^*}$ gives $\Upsilon_h = 0.69$ across the whole epoch, so the Higgs abundance is always about 31 percent below the equilibrium yield. Separately, comparing the momentum-exchanging scattering rate $\Gamma_{hq \to hq}$ with the fusion rate $\Gamma_{\mathrm{fusion}}$ shows the two curves cross at $T = 25$ GeV; below that temperature scattering cannot keep up with production, and the Higgs momentum distribution is governed by the fusion process rather than by the ambient temperature. The paper presents this as the first known setting in which kinetic nonequilibrium coexists with chemical equilibrium in a relativistic plasma.
Load-bearing premise
The calculation's load-bearing premise, stated around Eq. (21), is that recombination of three on-shell particles into a Higgs is so weak-interaction-suppressed that it can be omitted; if a future calculation finds that reverse rate non-negligible, the stationary fugacity moves back toward 1 and the chemical nonequilibrium result collapses.
Editorial extensions
If this is right
- Any early-Universe computation that assumes full chemical equilibrium for the Higgs overestimates its number density by roughly 31 percent throughout the 130 GeV to 10 GeV QGP epoch.
- Below 25 GeV, the Higgs population was kinetically cold: individual Higgs bosons decayed before enough scattering events occurred to imprint the plasma temperature on their momenta.
- The persistent nonequilibrium spans the electroweak phase-transition region, so the paper argues the transition did not have to end rapidly, leaving more room for nonequilibrium electroweak baryogenesis.
- The Higgs-to-baryon density ratio, even with the 0.69 fugacity, is enormous (about $10^5$ at 10 GeV), so the out-of-equilibrium Higgs population vastly outnumbers the matter-antimatter asymmetry.
Reading between the lines
- The same detailed-balance-breaking logic should apply to any near-threshold resonance whose mass is just below the pair-production threshold of a strongly coupled heavy partner; the paper does not explore that generalization.
- A full computation of the inverse $3 \to 1$ recombination rate, which the paper drops as weak-suppressed, is the natural test: it would determine whether 0.69 is robust or whether a small reverse rate already moves the fugacity back toward 1.
- If the cold-Higgs effect is real, it should also affect the energy flow between the Higgs sector and the rest of the plasma below 25 GeV, which could feed back into the timing of electroweak processes; this is an implication the paper leaves for future kinetic-theory work.
- The same rate comparison could be run for laboratory heavy-ion collisions, where the temperature and expansion time scales differ; a detectable out-of-equilibrium Higgs or near-threshold scalar signal would test the mechanism outside cosmology.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript studies the chemical and kinetic equilibration of the Higgs boson in the early-universe quark-gluon plasma over the temperature range 130 GeV > T > 10 GeV. It introduces a fugacity parameter for the Higgs, writes a rate equation that includes 2-to-1 fusion reactions from bottom, charm, gluon, and tau channels and the decay h -> WW*, ZZ*, and deliberately omits the inverse of the virtual-boson decay channel. Setting the time derivative of the fugacity to zero yields the claimed value Upsilon_h = 0.69. The paper separately compares the total Higgs scattering rate with the fusion and decay rates and concludes that the Higgs momentum distribution is 'cold' for T < 25 GeV. Both central claims rest on the rate equation in Sec. 4 and on the rate comparison in Sec. 5.
Significance. If correct, the paper's claims would imply a persistent 31% chemical underabundance of Higgs bosons in the primordial QGP and a kinetically under-thermalized Higgs population below 25 GeV, with possible implications for electroweak baryogenesis scenarios. The paper is transparent about its main assumption: it explicitly states in Secs. 4 and 6 that the inverse of h -> WW*, ZZ* is omitted because it is supposedly higher order in weak couplings. I do not see a circularity in the definition of Upsilon_h, and there is no curve fitting to data. However, the paper does not provide machine-checked proofs, reproducible code, or an independent numerical derivation of the headline value; the number 0.69 is asserted rather than derived, and the central rate equation is truncated on a physical assumption that is not demonstrated.
major comments (3)
- [Sec. 4, Eqs. (21)-(22)] The omission of the inverse reaction to h -> WW*, ZZ* is the load-bearing step and is not justified in the manuscript. The time-reversed process, for example W f fbar -> h or the four-fermion inverse of h -> WW* -> 4f, has the same powers of the weak couplings g and g' as the forward decay when all initial particles are on shell, and in a thermal bath with all non-Higgs species in equilibrium detailed balance requires a reverse rate per volume equal to n_h^eq Gamma_{h->WW*,ZZ*}. Adding that term changes Eq. (22) to (1/V) dN_h/dt = (1 - Upsilon_h)(R_fusion + R_decay), whose stationary solution is Upsilon_h = 1, not 0.69. The text itself says in Sec. 4 that 'we have omitted the back-reaction process entirely,' and in Sec. 6 that the decay 'does not have a back reaction,' but this is exactly the point that needs proof. The assertion that the 3-to-1 inverse is suppressed by powers of g^2 or g'^2 is not a substitute for a phase-space and amplitude computation, and the chemical nonequilibrium claim does not follow without it.
- [Sec. 4, Eq. (29)] Equation (29) is internally inconsistent. The stationary solution of Eq. (26) is Upsilon_h = Gamma_fusion/(Gamma_fusion + Gamma_H->WW*,ZZ*), which equals Gamma_fusion/Gamma_decay only in the unphysical limit Gamma_decay = 0. The manuscript writes both expressions and claims both equal 0.69; this would require Gamma_fusion/Gamma_decay = 2.23, not 0.69. No calculation, table, or figure is given from which Gamma_fusion and Gamma_decay can be read off, so the abstract's central number 0.69 is not derived anywhere in the text, even under the authors' own truncated equation.
- [Sec. 5, Fig. 5] The kinetic nonequilibrium conclusion is not established by the comparison of total rates in Fig. 5. Whether the Higgs momentum distribution is 'cold' depends on the momentum-transfer rate from h+b/t scattering relative to the Hubble expansion rate and relative to the production and decay rates, and on the actual shape of the distribution produced by the fusion process. The manuscript solves no Boltzmann equation for f_h(p) and displays no computed momentum distribution. The statement in Sec. 5 that the produced distribution 'is not informed about ambient temperature' therefore goes beyond what the presented rate comparison can support.
minor comments (4)
- [Abstract and Sec. 1] The Abstract contains grammatical errors, including 'the Higgs bosons is always out of chemical abundance equilibrium'; the rest of the text also has typos such as 'elctro-weak' in Sec. 1 and 'the have gauge boson pairs' in Sec. 6.
- [Eq. (17)] The notation '+/-1' in Eq. (17) is not defined; the text should state that the upper sign applies to bosons and the lower sign to fermions, and should specify the statistical factors entering Phi(p3) more carefully.
- [Sec. 4, Eqs. (24)-(28)] The discussion conflates a stationary approximation with solving a time-dependent equation. Equations (24) and (26) are differential equations, but the paper sets dUpsilon_h/dt = 0 and never integrates the time evolution; the phrase 'we solve the population equation' in Sec. 6 should be replaced by a statement that the stationary fixed point is analyzed.
- [Fig. 4 caption] The Fig. 4 caption contains a typo: 'op quark scattering' should read 'top quark scattering', and the Fig. 2 caption has the doubled article 'the the number density'.
Circularity Check
No definitional circularity; the central 0.69 value is an unsupported model output rather than a recycled input.
full rationale
The paper's derivation chain is not circular in the sense used by this analysis. The fugacity value Υ_h = 0.69 is presented as the stationary solution of the population equation (26), not as a fitted parameter or as a restatement of a definition. The model assumes that h → WW*, ZZ* has no effective back-reaction and writes Eq. (21) with only a loss term; solving that equation then necessarily gives Υ_h < 1. Whether that assumption is physically justified is a soundness question, not a circularity question: the omission of the inverse reaction is an input hypothesis, and the claim of chemical nonequilibrium is a consequence of that hypothesis, but the hypothesis is not derived from the conclusion. The self-citations in the paper (e.g., refs. [10], [11], [14], [15]) supply the standard fugacity formalism and the inverse-decay rate formulas; these are tools and are not the source of the target value 0.69. The main substantive problems are non-circular: Eq. (29) is internally inconsistent as written, since Υ_h = Γ_fusion/(Γ_fusion + Γ_decay) cannot equal Γ_fusion/Γ_decay = 0.69, and no independent computation of the rates leading to 0.69 is shown in the text. These are correctness/completeness defects, not a reduction of the result to its inputs. A score of 2 reflects the presence of several self-citations in the formalism while the central claim still has independent algebraic content, albeit incompletely and inconsistently presented.
Assumptions & free parameters
free parameters (1)
- Higgs fugacity Υ_h =
0.69
assumptions (4)
- domain assumption The inverse 3-to-1 reactions to h -> WW*, ZZ* are negligible and can be omitted from Eq. (21).
- domain assumption All QGP species other than the Higgs are in chemical equilibrium with fugacity Υ_i = 1.
- domain assumption The stationary condition dΥ_h/dt = 0 with Hubble terms dropped is valid.
- standard math Boltzmann limit m_i/T >> 1 for heavy particles and nonrelativistic rate formulas.
Cite this review
Pith. "Pith review of Higgs Thermal Nonequilibrium in Primordial QGP." pith.science (2026). https://pith.science/paper/OU5ZFITJ
@misc{pith2026250203598,
author = {Pith},
title = {Pith review of: Higgs Thermal Nonequilibrium in Primordial QGP},
year = {2026},
howpublished = {\url{https://pith.science/paper/OU5ZFITJ}},
note = {Machine review of arXiv:2502.03598}
}
abstract
In this work we investigate the chemical and kinetic nonequilibrium dynamics of the Higgs boson during the primordial Universe QGP (quark-gluon plasma) epoch $130\mathrm{\,GeV}>T>10\mathrm{\,GeV}$. We show that the Higgs bosons is always out of chemical abundance equilibrium with a fugacity $\Upsilon_h = 0.69$ due to virtual decay channels. Additionally, Higgs momentum distribution is found to be ``cold'' for $T<25$\,GeV, since the scattering rate drops below the production rate.
Figures
Figures from the paper (2 more)
Reference graph
Works this paper leans on
-
[1]
Review of Particle Physics 2022,
R. L. Workman et al. [Particle Data Group], “Review of Particle Physics 2022,” PTEP 2022, 083C01 (2022) doi:10.1093/ptep/ptac097
-
[2]
Higgs Boson Decay to One Real and One Virtual W Boson,
E. W. N. Glover, J. Ohnemus and S. S. D. Willenbrock, “Higgs Boson Decay to One Real and One Virtual W Boson,” Phys. Rev. D 37 (1988), 3193 doi:10.1103/PhysRevD.37.3193
-
[3]
G. V. Vereshchagin and A. G. Aksenov, Relativistic Kinetic Theory with Applications in Astrophysics and Cosmology, Cambridge University Press (2017), Online ISBN: 9781107261365, doi:10.1017/9781107261365
-
[4]
Matter and Antimatter in the Universe,
L. Canetti, M. Drewes and M. Shaposhnikov, “Matter and Antimatter in the Universe,” New J. Phys. 14, 095012 (2012) doi:10.1088/1367-2630/14/9/095012 [arXiv:1204.4186 [hep-ph]]
arXiv 2012
-
[5]
D. E. Morrissey and M. J. Ramsey-Musolf, “Electroweak baryogenesis,” New J. Phys. 14, 125003 (2012) doi:10.1088/1367-2630/14/12/125003 [arXiv:1206.2942 [hep-ph]]
arXiv 2012
-
[6]
Quarks to Cosmos: Particles and Plasma in Cosmo- logical evolution,
J. Rafelski, J. Birrell, C. Grayson, A. Steinmetz and C. T. Yang, “Quarks to Cosmos: Particles and Plasma in Cosmo- logical evolution,” [arXiv:2409.19031 [hep-ph]]
-
[7]
Hadrons and Quark–Gluon Plasma,
J. Letessier and J. Rafelski, “Hadrons and Quark–Gluon Plasma,” Cambridge University Press, 2002, ISBN 978-1-009- 29075-3, 978-0-511-53499-7, 978-1-009-29070-8, 978-1-009-29073-9 doi:10.1017/9781009290753
-
[8]
Quarkochemistry in Relativistic Heavy Ion Collisions,
T. S. Bir´ o and J. Zim´ anyi, “Quarkochemistry in Relativistic Heavy Ion Collisions,” Phys. Lett. B 113, 6 (1982) doi:10.1016/0370-2693(82)90097-1; Report number: KFKI-1981-69 (Budapest)
Show all 15 references
-
[9]
Strangeness Production in the Quark-Gluon Plasma,
J. Rafelski and B. M¨ uller, “Strangeness Production in the Quark-Gluon Plasma,” Phys. Rev. Lett. 48, 1066 (1982) [erratum: Phys. Rev. Lett. 56, 2334 (1986)] doi:10.1103/PhysRevLett.48.1066; Report number: Print-82-0048 (Frankfurt)
1982 doi
-
[10]
Strangeness in Relativistic Heavy Ion Collisions,
P. Koch, B. Muller and J. Rafelski, “Strangeness in Relativistic Heavy Ion Collisions,” Phys. Rept. 142, 167-262 (1986) doi:10.1016/0370-1573(86)90096-7
1986 doi
-
[11]
Gluon production, cooling and entropy in nuclear collisions,
J. Letessier, J. Rafelski and A. Tounsi, “Gluon production, cooling and entropy in nuclear collisions,” Phys. Rev. C 50, 406-409 (1994) doi:10.1103/PhysRevC.50.406
1994 doi
-
[12]
Relic neutrinos: Physically consistent treatment of effective number of neutrinos and neutrino mass,
J. Birrell, C. T. Yang, P. Chen and J. Rafelski, “Relic neutrinos: Physically consistent treatment of effective number of neutrinos and neutrino mass,” Phys. Rev. D 89 (2014), 023008 doi:10.1103/PhysRevD.89.023008 [arXiv:1212.6943 [astro-ph.CO]]
2014 arXiv
-
[13]
From Quark-Gluon Universe to Neutrino Decoupling: 200 < T < 2MeV,
M. J. Fromerth, I. Kuznetsova, L. Labun, J. Letessier and J. Rafelski, “From Quark-Gluon Universe to Neutrino Decoupling: 200 < T < 2MeV,” Acta Phys. Polon. B 43, no.12, 2261-2284 (2012) doi:10.5506/APhysPolB.43.2261 [arXiv:1211.4297 [nucl-th]]
2012 arXiv
-
[14]
Unstable Hadrons in Hot Hadron Gas in Laboratory and in the Early Universe,
I. Kuznetsova and J. Rafelski, “Unstable Hadrons in Hot Hadron Gas in Laboratory and in the Early Universe,” Phys. Rev. C 82, 035203 (2010) doi:10.1103/PhysRevC.82.035203
2010 doi
-
[15]
Pion and muon production in e-, e+, gamma plasma,
I. Kuznetsova, D. Habs and J. Rafelski, “Pion and muon production in e-, e+, gamma plasma,” Phys. Rev. D 78 (2008), 014027 doi:10.1103/PhysRevD.78.014027 [arXiv:0803.1588 [hep-ph]]
2008 arXiv
Reviewed August 9, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.