{"id":"455b6dd0-96e3-4bba-96d6-aefc4f2f039b","arxiv_id":"2502.03598","paper_version":2,"verdict":"REJECT","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"high","formal_verification":"none","parameter_count":1,"one_line_summary":"The paper predicts that primordial Higgs bosons maintained a constant chemical fugacity of 0.69 and a nonthermal, cold momentum distribution below T = 25 GeV.","lead":"A new calculation claims that Higgs bosons in the hot early universe stayed out of chemical equilibrium, with an abundance deficit of about 31 percent, and developed a cold momentum distribution below 25 GeV. The cause would be that Higgs decays into virtual W and Z bosons have no efficient inverse reaction, a detail that could matter for modeling the electroweak epoch and baryogenesis.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Omitted reverse of h -> WW*, ZZ* in Eq. (21) is the load-bearing step; detailed balance restores Upsilon_h = 1, so the chemical nonequilibrium claim does not follow.","rationale":"The reader's weakest_assumption identifies the same load-bearing concern: Eq. (21) omits the inverse of the virtual-decay channel, and the justification for that omission is not physically sound. In a thermal bath with all non-Higgs species in equilibrium, unitarity and detailed balance force the reverse rate to equal the forward decay rate at equilibrium. Including the reverse term changes the stationary fugacity from 0.69 to 1, which removes the paper's central chemical nonequilibrium result. This is a correctness risk in the rate-equation input, not merely a numerical uncertainty. The additional arithmetic problem in Eq. (29) strengthens the case that the quantitative claim is not supported by the equations as written. The kinetic nonequilibrium claim is also questionable because the main criterion used is comparison of scattering with production rather than with the Hubble rate or decay rate, but the chemical claim alone is sufficient to sustain the original REJECT verdict. I therefore see no reason to change the reader's verdict.","tokens_in":9784,"tokens_out":9709,"duration_ms":104264,"concrete_test":"Re-derive the stationary fugacity with the exact time-reversed process included: add the reverse rate R_rev = n_h^eq Gamma_{h->WW*,ZZ*} (valid when all decay products are at Upsilon_i = 1) to the production side of Eq. (22), so the virtual-decay channel contributes (1 - Upsilon_h) R_decay instead of -Upsilon_h R_decay. Solve Eq. (26) with dUpsilon_h/dt = 0. If the result is Upsilon_h = 1 rather than 0.69, as detailed balance requires, the central claim is unsupported. As a numerical cross-check, evaluate the thermally averaged 3->1 (or 4->1) reverse cross section from the same Glover-Ohnemus-Willenbrock matrix element at T = 100, 50, 20 GeV and confirm it is comparable to the forward decay rate.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central chemical claim depends entirely on omitting the reverse of h -> WW*, ZZ* from the population equation. In Eq. (21) the virtual-decay channel enters only as -Upsilon_h R_{h->WW*,ZZ*}; the authors state (Sec. 4 and Sec. 6) that the inverse is suppressed by extra powers of weak couplings and drop it entirely. That suppression argument is not valid: the time-reversed amplitude for, e.g., W f fbar -> h (or the four-fermion inverse of h -> WW* -> 4f) has the same power of g and g' as the forward decay, and integrated over a thermal bath with all non-Higgs species at fugacity 1, detailed balance requires a reverse rate per volume equal to n_h^eq Gamma_{h->WW*,ZZ*}. Adding that term turns the virtual-decay contribution into (1 - Upsilon_h) R_decay, so the stationary solution of Eq. (26) becomes Upsilon_h = (R_fusion + R_decay)/(R_fusion + R_decay) = 1, not 0.69. The claimed nonequilibrium is therefore an artifact of the omitted back-reaction, not a computed property. Eq. (29) is also internally inconsistent: Gamma_fusion/(Gamma_fusion + Gamma_decay) = 0.69 would imply Gamma_fusion/Gamma_decay ~ 2.23, not 0.69, and no independent derivation of the number is given.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":10102,"tokens_out":4306,"duration_ms":42211,"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":[{"comment":"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.","section":"Sec. 4, Eqs. (21)-(22)"},{"comment":"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.","section":"Sec. 4, Eq. (29)"},{"comment":"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.","section":"Sec. 5, Fig. 5"}],"minor_comments":[{"comment":"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.","section":"Abstract and Sec. 1"},{"comment":"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.","section":"Eq. (17)"},{"comment":"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.","section":"Sec. 4, Eqs. (24)-(28)"},{"comment":"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'.","section":"Fig. 4 caption"}],"recommendation":"reject","confidential_remarks":"This is a soundness rejection rather than a scope or novelty rejection. The central chemical claim depends on the omission of the reverse reaction in Eq. (21), and the text explicitly admits that omission without providing a valid suppression argument. The internal inconsistency of Eq. (29) is a separate, decisive problem because the headline value 0.69 is not derived. Even if the authors were to include the back-reaction and redo the stationary solution, the paper's main conclusion would likely be replaced by Upsilon_h = 1, so the issue cannot be repaired by local revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the central number in this paper, Υ_h = 0.69, is not supported. The population equation (21) drops the inverse of h → WW*, ZZ*, and the justification given in Secs. 4 and 6 is not valid. In a thermal bath with the final-state fermions in equilibrium, the reverse rate per volume is forced by detailed balance to be n_h^eq Γ_decay. Adding that term turns the loss channel into (1 − Υ_h)R_decay, and the stationary solution of Eq. (26) is Υ_h = 1, not 0.69. The chemical non-equilibrium claim is an artifact of the omitted back-reaction. Eq. (29) is also internally inconsistent: Γ_fusion/(Γ_fusion+Γ_decay) is not equal to Γ_fusion/Γ_decay, and no derivation of the value 0.69 is provided.\n\nWhat the paper does well: it assembles the relevant rates—bottom and top fusion, virtual W/Z decay, and Higgs–quark scattering—into a single temperature-dependent inventory, and the framing of the question is legitimate. The entropy-temperature treatment is standard and clean. The rate plots would be useful to someone starting a more rigorous study of Higgs thermalization in the early universe.\n\nThe kinetic claim is weaker than the abstract suggests. The authors compare the scattering rate to fusion and decay, but the conventional freeze-out criterion is scattering versus Hubble expansion; their own Fig. 4 shows both scattering and fusion exceed H by many orders of magnitude across the whole range. If they intend the lifetime argument (a Higgs decays before it scatters), they need to state that criterion explicitly and, ideally, solve the Boltzmann equation for the momentum distribution. As written, the 'cold' conclusion is not established.\n\nBottom line: this is a preliminary study with a load-bearing error in the chemical claim. It is not ready for publication. I would not invest referee time in the current version; a desk reject with an indication that the detailed-balance issue must be addressed is appropriate. The rate inventory in Figs. 3–5 could become the basis for a better paper if the authors revise the rate equations and either correct or remove the chemical non-equilibrium result.","headline":"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.","tokens_in":10637,"tokens_out":7821,"would_cite":false,"duration_ms":74086,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"In the primordial quark-gluon plasma, the Higgs boson was chemically out of equilibrium with fugacity 0.69 and kinetically cold below 25 GeV.","keywords":["Higgs boson","primordial quark-gluon plasma","chemical nonequilibrium","fugacity","kinetic nonequilibrium","detailed balance","virtual W/Z decay","early Universe"],"falsifier":"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.","tokens_in":9567,"feed_emoji":"⚛️","tokens_out":9411,"duration_ms":75954,"temperature":0.7,"pith_summary":"This paper claims that in the quark-gluon plasma of the early Universe, between temperatures of about 130 GeV and 10 GeV, the Higgs boson never reached chemical equilibrium: its fugacity, the factor that multiplies the equilibrium particle distribution, sat at $\\Upsilon_h = 0.69$, meaning a 31 percent deficit in abundance. The deficit arises because the Higgs is produced by two-particle fusion, mainly bottom-quark fusion, but disappears through a decay channel, $h \\to WW^*, ZZ^*$, in which one gauge boson is virtual and necessarily falls apart; the inverse three-particle recombination is weak-interaction suppressed and is dropped from the rate equation. The paper also claims that below 25 GeV the two-body scattering rate for Higgs–quark collisions falls below the production rate, so the Higgs momentum distribution is \"cold,\" shaped by production kinematics rather than by the plasma temperature. If correct, the Higgs is a concrete example of a species that fails full thermal equilibrium in the early Universe despite rapid strong-interaction rates during the expansion era.","feed_headline":"Primordial Higgs never reached equilibrium: 31% short","feed_subtitle":"Virtual W and Z decays break detailed balance, leaving the Higgs cold below 25 GeV.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Supplies the Higgs mass, total decay width, branching ratios (b bbar, W, Z), and the present baryon-to-photon ratio used to set the early-Universe abundances.","marker":"[1]"},{"why":"Establishes the h -> W W* (and analogously Z Z*) decay with one real and one virtual gauge boson, the one-way loss channel that breaks detailed balance.","marker":"[2]"},{"why":"Provides the relativistic kinetic-theory context for why decays through virtual particles are not balanced by inverse reactions, the premise for omitting the back-reaction.","marker":"[3]"},{"why":"Supplies the thermal two-body scattering rate formula R_12->34 used to compute Higgs-quark kinetic scattering rates.","marker":"[7]"},{"why":"Provides the thermal inverse-decay rate R_12->3 formula used for Higgs production by fermion fusion.","marker":"[14]"},{"why":"Together with [14], supports the thermal rate formalism for unstable-particle production applied to the Higgs fusion rates.","marker":"[15]"}],"fun_headline_variants":["Primordial Higgs off equilibrium by 31% in QGP","Virtual W and Z decays freeze Higgs abundance at 0.69","Higgs momentum went cold below 25 GeV in early plasma","First QGP showing kinetic and chemical nonequilibrium together"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Primordial Higgs off equilibrium by 31% in QGP","Virtual W and Z decays freeze Higgs abundance at 0.69","Higgs momentum went cold below 25 GeV in early plasma","First QGP showing kinetic and chemical nonequilibrium together"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000207,"raw_usage":{"total_tokens":1373,"prompt_tokens":890,"completion_tokens":483,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":506,"completion_tokens_details":{"reasoning_tokens":412}},"tokens_in":506,"tokens_out":483,"duration_ms":4263,"temperature":1.0,"reasoning_tokens":412,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-09T04:23:10.973822+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"Higgs Boson Decay to One Real and One Virtual W Boson,","cited_arxiv_id":null,"evidence_quote":"Establishes the h -> W W* (and analogously Z Z*) decay with one real and one virtual gauge boson, the one-way loss channel that breaks detailed balance."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the relativistic kinetic-theory context for why decays through virtual particles are not balanced by inverse reactions, the premise for omitting the back-reaction."},{"cited_title":"Unstable Hadrons in Hot Hadron Gas in Laboratory and in the Early Universe,","cited_arxiv_id":null,"evidence_quote":"Provides the thermal inverse-decay rate R_12->3 formula used for Higgs production by fermion fusion."},{"cited_title":"Pion and muon production in electron-positron photon plasma","cited_arxiv_id":"0803.1588","evidence_quote":"Together with [14], supports the thermal rate formalism for unstable-particle production applied to the Higgs fusion rates."}],"review_version":1}