Probing the Rare Four-Bottom Higgs Decay Hto bbar b bbar b at the HL-LHC and ILC
Pith reviewed 2026-05-20 04:25 UTC · model grok-4.3
The pith
The Standard Model predicts a branching ratio of order 1.6 times 10 to the minus 3 for the rare Higgs decay to four bottom quarks, with destructive interference among amplitudes, and shows this channel can reach observable significance at a
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The central claim is that the rare decay H to b bbar b bbar has a branching ratio of order 1.6 times 10 to the minus 3 when all leading amplitudes are included, with destructive interference playing a phenomenologically relevant role. This decay can be probed in associated production modes, yielding a statistical significance of approximately 3.5 in pp to W H to W b bbar b bbar at 14 TeV with the full HL-LHC dataset and a significance above 5 sigma in e plus e minus to Z H to Z b bbar b bbar at 250 GeV with only 300 inverse femtobarns using a boosted decision tree analysis that exploits four-b kinematic correlations.
What carries the argument
The boosted decision tree multivariate analysis that exploits kinematic correlations among the four-bottom final state in associated Higgs production to separate signal from backgrounds.
Load-bearing premise
The multivariate boosted decision tree analysis can reliably exploit kinematic correlations among the four-b final state without significant mismodeling of backgrounds or detector effects in the high-purity working points.
What would settle it
An observed event rate in the four-bottom final state from W H or Z H production that deviates by more than the combined statistical and systematic uncertainties from the rate expected for a branching ratio of 1.6 times 10 to the minus 3 would falsify the central claim.
Figures
read the original abstract
We propose the rare SM Higgs decay $H\to b\bar b b\bar b$ as a probe of the structure of Higgs interactions with bottom quarks and gauge bosons, and as a baseline for searches for new physics producing four-bottom final states in Higgs decays. We compute the leading contributions to this decay, including the dominant $H\to b\bar b g\to b\bar b b\bar b$ topology, the sizeable $H\to ZZ^\ast\to b\bar b b\bar b$ channel, and the loop-induced $H\to gg\to b\bar b b\bar b$ contribution. We find a branching ratio of order $1.6\times10^{-3}$ and show that destructive interference among the leading amplitudes is phenomenologically relevant. We demonstrate that this decay can be probed in associated Higgs production at both the HL-LHC and the ILC. For $pp\to WH\to Wb\bar b b\bar b$ at $\sqrt{s}=14$ TeV, we use a multivariate analysis based on boosted decision trees to exploit correlations among the four-$b$ kinematic observables. At $3000~{\rm fb}^{-1}$, the statistical significance reaches about $3.5$, while a tighter high-purity working point gives $S/B\simeq5\%$ with significance close to $3\sigma$. A combined high-luminosity LHC dataset could therefore make this rare decay observable. For $e^+e^-\to ZH\to Zb\bar b b\bar b$ at the ILC with $\sqrt{s}=250$ GeV, we demonstrate that the cleaner collider environment gives a high-purity signal sample. In the nominal setup, the multivariate analysis gives a significance above $5\sigma$ already at $300~{\rm fb}^{-1}$. At integrated luminosities of order $1-3~{\rm ab}^{-1}$, the branching ratio can be measured with several-percent precision.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper computes the leading amplitudes for the rare SM Higgs decay H→bbbb, including the dominant H→bbg→bbbb topology, H→ZZ*→bbbb, and loop-induced H→gg→bbbb, reporting a branching ratio of order 1.6×10^{-3} with phenomenologically relevant destructive interference. It then presents multivariate BDT analyses for associated production pp→WH→Wbbbb at the HL-LHC (√s=14 TeV) and e+e−→ZH→Zbbbb at the ILC (√s=250 GeV), projecting statistical significances of ~3.5σ at 3000 fb^{-1} (or ~3σ in a high-purity S/B≃5% working point) and >5σ already at 300 fb^{-1}, respectively.
Significance. If the BDT-based projections are robust, the work supplies a useful SM baseline for four-bottom final states and a probe of Higgs-bottom and Higgs-gauge couplings. The explicit inclusion of interference among amplitudes is a positive feature of the theoretical calculation.
major comments (2)
- [HL-LHC analysis] HL-LHC analysis section: The quoted 3.5σ significance (and the high-purity ~3σ point) is obtained from a BDT trained on simulated four-b kinematics; no systematic variations of jet energy scale, b-tagging, or irreducible backgrounds (tt̄, Z+bb) are reported, nor is any data-driven validation of the background modeling in the high-purity region described. This directly affects the load-bearing observability claim.
- [Theoretical calculation] Branching-ratio calculation: While destructive interference is stated to be relevant, the manuscript does not provide the explicit interference terms or the relative sizes of the H→bbg, H→ZZ*, and loop-induced contributions in a single equation or table, making it difficult to verify the net 1.6×10^{-3} result.
minor comments (2)
- [Abstract] The abstract states the ILC significance is “above 5σ already at 300 fb^{-1}” but does not specify the exact working point or luminosity scaling used for the 1–3 ab^{-1} precision projection.
- Notation for the four-b final state is occasionally inconsistent (b b-bar b b-bar vs. bbbb); a single convention should be adopted throughout.
Simulated Author's Rebuttal
We thank the referee for the thorough review and constructive feedback on our manuscript. We address the major comments below and have revised the manuscript accordingly to improve clarity and robustness.
read point-by-point responses
-
Referee: [HL-LHC analysis] HL-LHC analysis section: The quoted 3.5σ significance (and the high-purity ~3σ point) is obtained from a BDT trained on simulated four-b kinematics; no systematic variations of jet energy scale, b-tagging, or irreducible backgrounds (tt̄, Z+bb) are reported, nor is any data-driven validation of the background modeling in the high-purity region described. This directly affects the load-bearing observability claim.
Authors: We acknowledge the validity of this observation. Our analysis focuses on the potential statistical significance using a BDT on simulated signal and background samples to demonstrate the feasibility of observing this rare decay. However, we agree that a complete assessment requires consideration of systematic uncertainties. In the revised manuscript, we will add a paragraph discussing the expected systematic uncertainties based on current LHC performance (e.g., jet energy scale uncertainties of a few percent, b-tagging efficiencies), and note that the quoted significances are statistical. We will also clarify that a full experimental analysis would incorporate data-driven background estimation. This revision will better contextualize the observability claims without altering the core methodology. revision: yes
-
Referee: [Theoretical calculation] Branching-ratio calculation: While destructive interference is stated to be relevant, the manuscript does not provide the explicit interference terms or the relative sizes of the H→bbg, H→ZZ*, and loop-induced contributions in a single equation or table, making it difficult to verify the net 1.6×10^{-3} result.
Authors: We appreciate this suggestion for improving transparency. In the updated version of the paper, we will introduce a new table that lists the individual branching ratios for each contributing process (H→bbg, H→ZZ*, loop-induced H→gg), the pairwise interference terms, and the total after interference. This will explicitly demonstrate the destructive interference and allow readers to verify the net branching ratio of approximately 1.6×10^{-3}. We believe this addition will address the concern directly. revision: yes
Circularity Check
No significant circularity; SM amplitudes and MC-based significances are independent of target observables
full rationale
The paper computes the H→bbbb branching ratio from explicit leading-order SM amplitudes (H→bbg, H→ZZ*, loop-induced gg) with interference, using standard Feynman rules and external parameters. Collider projections (3.5σ at HL-LHC, >5σ at ILC) follow from BDT classification on simulated event kinematics for associated production, without fitting any parameter to the four-b final state or renaming fitted inputs as predictions. No self-citation chains, uniqueness theorems, or ansatze are invoked to force the central results. The derivation remains self-contained against SM benchmarks and standard simulation tools.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Standard Model Feynman rules and parton shower modeling for b-quark fragmentation are accurate for this process.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We compute the leading contributions to this decay, including the dominant H→b¯bg→b¯bb¯b topology, the sizeable H→ZZ∗→b¯bb¯b channel, and the loop-induced H→gg→b¯bb¯b contribution. We find a branching ratio of order 1.6×10−3 and show that destructive interference among the leading amplitudes is phenomenologically relevant.
-
IndisputableMonolith/Foundation/AlexanderDuality.leanalexander_duality_circle_linking unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
For pp→WH→Wb¯bb¯b at √s=14 TeV, we use a multivariate analysis based on boosted decision trees to exploit correlations among the four-b kinematic observables.
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
ForH→b ¯b→b ¯bg→b ¯bb¯b, the distribution is stretched horizontally, withMb1b2 concentrated in the 80–115 GeV range whileMb3b4 peaks at 10–20 GeV
-
[2]
ForH→ZZ ∗ →b ¯bb¯b, the distribution is vertical, withMb1b2 clustering around theZresonance andM b3b4 in the 10–35 GeV range
-
[3]
ForH→gg→b ¯bb¯b, thedistributionformsanarrowdiagonalstripextendingfrom(M b1b2, Mb3b4) = (115,10)GeV to(70,60)GeV. Altogether, the two-dimensional approach reveals distinctive horizontal, vertical, and diagonal patterns for the three subprocesses. These features, largely washed out in one-dimensional projections, provide clear discrimination and underscore...
-
[4]
the transverse momenta of the four selectedb-tagged jets,pT b1,p T b2,p T b3, andp T b4
-
[5]
the energies of the four selectedb-tagged jets in the reconstructed Higgs rest frame, defined as the centre-of-mass frame of the selected four-bsystem. These energies are ordered asE1 > E2 > E3 > E4, where the ordering refers to the reconstructed Higgs rest frame:E1,E 2,E 3, andE 4
-
[6]
the minimum and maximum invariant masses among all possible two-bcombinations,Mmin bibj and M max bibj
-
[7]
the invariant mass of the two most energeticbjets in the reconstructed Higgs rest frame,M12
-
[8]
the invariant mass of the two least energeticbjets in the reconstructed Higgs rest frame,M34
-
[9]
the invariant mass of the three least energeticbjets in the reconstructed Higgs rest frame,M234
-
[10]
the invariant mass of the selected four-bsystem,M4b
-
[11]
Altogether, the BDTG is trained using4 + 4 + 2 + 1 + 1 + 1 + 1 + 6 = 20input variables
the cosines of the angles between all pairs of the four selectedbjets in the reconstructed Higgs rest frame,cosθ ij, giving six variables. Altogether, the BDTG is trained using4 + 4 + 2 + 1 + 1 + 1 + 1 + 6 = 20input variables. The left panel of Fig. 15 shows the BDTG response for the full signal and full background samples. The two distributions are clear...
work page 2041
-
[12]
Georges Aad et al. Observation of a new particle in the search for the Standard Model Higgs boson with the ATLAS detector at the LHC.Phys. Lett., B716:1–29, 2012. doi:10.1016/j.physletb.2012.08.020
work page internal anchor Pith review doi:10.1016/j.physletb.2012.08.020 2012
-
[13]
Observation of a new boson at a mass of 125 GeV with the CMS experiment at the LHC
Serguei Chatrchyan et al. Observation of a new boson at a mass of 125 GeV with the CMS experiment at the LHC. Phys. Lett., B716:30–61, 2012. doi:10.1016/j.physletb.2012.08.021
-
[14]
Observation of the Diphoton Decay of the Higgs Boson and Measurement of Its Properties
Vardan Khachatryan et al. Observation of the Diphoton Decay of the Higgs Boson and Measurement of Its Properties. Eur. Phys. J., C74(10):3076, 2014. doi:10.1140/epjc/s10052-014-3076-z
-
[15]
Georges Aad et al. Measurement of Higgs boson production in the diphoton decay channel in pp collisions at center-of-mass energies of 7 and 8 TeV with the ATLAS detector.Phys. Rev., D90(11):112015, 2014. doi:10.1103/PhysRevD.90.112015
-
[16]
Measurement of the Properties of a Higgs Boson in the Four-Lepton Final State
Serguei Chatrchyan et al. Measurement of the Properties of a Higgs Boson in the Four-Lepton Final State. Phys. Rev., D89(9):092007, 2014. doi:10.1103/PhysRevD.89.092007
-
[17]
Georges Aad et al. Measurements of Higgs boson production and couplings in the four-lepton channel in pp collisions at center-of-mass energies of 7 and 8 TeV with the ATLAS detector.Phys. Rev., D91(1): 012006, 2015. doi:10.1103/PhysRevD.91.012006
-
[18]
Serguei Chatrchyan et al. Measurement of Higgs Boson Production and Properties in the WW Decay Channel with Leptonic Final States.JHEP, 01:096, 2014. doi:10.1007/JHEP01(2014)096
-
[19]
Observation and measurement of Higgs boson decays toWW ∗ with the ATLAS detector
Georges Aad et al. Observation and measurement of Higgs boson decays to WW∗ with the ATLAS detector. Phys. Rev., D92(1):012006, 2015. doi:10.1103/PhysRevD.92.012006
-
[20]
Georges Aad et al. Study of (W/Z)H production and Higgs boson couplings usingH→W W∗ decays with the ATLAS detector.JHEP, 08:137, 2015. doi:10.1007/JHEP08(2015)137. 32
-
[21]
Albert M. Sirunyan et al. Measurements of properties of the Higgs boson decaying to a W boson pair in pp collisions at√s=13 TeV. Phys. Lett., B791:96, 2019. doi:10.1016/j.physletb.2018.12.073
-
[22]
Evidence for the 125 GeV Higgs boson decaying to a pair ofτleptons.JHEP, 05:104, 2014
Serguei Chatrchyan et al. Evidence for the 125 GeV Higgs boson decaying to a pair ofτleptons.JHEP, 05:104, 2014. doi:10.1007/JHEP05(2014)104
-
[23]
Evidence for the Higgs-boson Yukawa coupling to tau leptons with the ATLAS detector
Georges Aad et al. Evidence for the Higgs-boson Yukawa coupling to tau leptons with the ATLAS detector. JHEP, 04:117, 2015. doi:10.1007/JHEP04(2015)117
-
[24]
Observation of the Higgs boson decay to a pair ofτleptons with the CMS detector
Albert M Sirunyan et al. Observation of the Higgs boson decay to a pair ofτleptons with the CMS detector. Phys. Lett., B779:283–316, 2018. doi:10.1016/j.physletb.2018.02.004
-
[25]
Observation ofH→b¯bdecays andV Hproduction with the ATLAS detector
Morad Aaboud et al. Observation ofH→b¯bdecays andV Hproduction with the ATLAS detector. Phys. Lett., B786:59–86, 2018. doi:10.1016/j.physletb.2018.09.013
-
[26]
A. M. Sirunyan et al. Observation of Higgs boson decay to bottom quarks.Phys. Rev. Lett., 121(12): 121801, 2018. doi:10.1103/PhysRevLett.121.121801
-
[27]
D. de Florian et al. Handbook of LHC Higgs Cross Sections: 4. Deciphering the Nature of the Higgs Sector. CERN YellowRep. Monogr., 2:1–869, 2017. doi:10.23731/CYRM-2017-002
-
[28]
David B. Kaplan. Flavor at SSC energies: A New mechanism for dynamically generated fermion masses. Nucl. Phys. B, 365:259–278, 1991. doi:10.1016/S0550-3213(05)80021-5
-
[29]
Inclusive search for a highly boosted Higgs boson decaying to a bottom quark- antiquark pair.Phys
Albert M Sirunyan et al. Inclusive search for a highly boosted Higgs boson decaying to a bottom quark- antiquark pair.Phys. Rev. Lett., 120(7):071802, 2018. doi:10.1103/PhysRevLett.120.071802
-
[30]
Aram Hayrapetyan et al. Measurement of the Higgs boson production via vector boson fusion and its decay into bottom quarks in proton-proton collisions at √s= 13 TeV. JHEP, 01:173, 2024. doi: 10.1007/JHEP01(2024)173
-
[31]
Higgs boson decay into four bottom quarks in the SM and beyond
Jun Gao. Higgs boson decay into four bottom quarks in the SM and beyond. JHEP, 08:174, 2019. doi:10.1007/JHEP08(2019)174
-
[32]
Christensen, and Alexander Pukhov
Alexander Belyaev, Neil D. Christensen, and Alexander Pukhov. CalcHEP 3.4 for collider physics within and beyond the Standard Model. Comput. Phys. Commun., 184:1729–1769, 2013. doi: 10.1016/j.cpc.2013.01.014
-
[33]
E. Boos, V. Bunichev, M. Dubinin, L. Dudko, V. Ilyin, A. Kryukov, V. Edneral, V. Savrin, A. Semenov, and A. Sherstnev. CompHEP 4.4: Automatic computations from Lagrangians to events.Nucl. Instrum. Meth. A, 534:250–259, 2004. doi:10.1016/j.nima.2004.07.096
-
[34]
Torbj˜A¶rn Sj˜A¶strand, Stefan Ask, Jesper R. Christiansen, Richard Corke, Nishita Desai, Philip Il- ten, Stephen Mrenna, Stefan Prestel, Christine O. Rasmussen, and Peter Z. Skands. An Introduction to PYTHIA 8.2.Comput. Phys. Commun., 191:159–177, 2015. doi:10.1016/j.cpc.2015.01.024
-
[35]
J. de Favereau, C. Delaere, P. Demin, A. Giammanco, V. Lema˜A®tre, A. Mertens, and M. Selvaggi. DELPHES 3, A modular framework for fast simulation of a generic collider experiment.JHEP, 02:057,
-
[36]
doi:10.1007/JHEP02(2014)057
work page internal anchor Pith review doi:10.1007/jhep02(2014)057 2014
-
[37]
AssociatedW Hproduction at hadron colliders: a fully exclusive QCD calculation at NNLO.Phys
Giancarlo Ferrera, Massimiliano Grazzini, and Francesco Tramontano. AssociatedW Hproduction at hadron colliders: a fully exclusive QCD calculation at NNLO.Phys. Rev. Lett., 107:152003, 2011. doi: 10.1103/PhysRevLett.107.152003
-
[38]
Higher-order QCD effects for asso- ciated WH production and decay at the LHC.JHEP, 04:039, 2014
Giancarlo Ferrera, Massimiliano Grazzini, and Francesco Tramontano. Higher-order QCD effects for asso- ciated WH production and decay at the LHC.JHEP, 04:039, 2014. doi:10.1007/JHEP04(2014)039
-
[39]
LHC Higgs Cross Section Working Group. Handbook of LHC Higgs Cross Sections: 14 TeV Report.https: //twiki.cern.ch/twiki/bin/view/LHCPhysics/CERNYellowReportPageAt14TeV. Accessed: 2025-01-20
work page 2025
-
[40]
Ansgar Denner, Stefan Dittmaier, Stefan Kallweit, and Alexander M¨ uck. HAWK 2.0: A Monte Carlo program for Higgs production in vector-boson fusion and Higgs strahlung at hadron colliders.Comput. Phys. Commun., 195:161–171, 2015. doi:10.1016/j.cpc.2015.04.021
-
[41]
Harlander, Anna Kulesza, Vincent Theeuwes, and Tom Zirke
Robert V. Harlander, Anna Kulesza, Vincent Theeuwes, and Tom Zirke. Soft gluon resummation for gluon-induced Higgs Strahlung.JHEP, 11:082, 2014. doi:10.1007/JHEP11(2014)082
-
[42]
Harlander, Heidi Rzehak, and Tom J
Lukas Altenkamp, Stefan Dittmaier, Robert V. Harlander, Heidi Rzehak, and Tom J. E. Zirke. Gluon-induced Higgs-strahlung at next-to-leading order QCD. JHEP, 02:078, 2013. doi: 10.1007/JHEP02(2013)078
-
[43]
Ansgar Denner, Stefan Dittmaier, Stefan Kallweit, and Alexander Muck. Electroweak corrections to Higgs-strahlung off W/Z bosons at the Tevatron and the LHC with HAWK.JHEP, 03:075, 2012. doi: 10.1007/JHEP03(2012)075
-
[44]
NNLO QCD corrections to the Higgs-strahlung processes at hadron colliders.Phys
Oliver Brein, Abdelhak Djouadi, and Robert Harlander. NNLO QCD corrections to the Higgs-strahlung processes at hadron colliders.Phys. Lett. B, 579:149–156, 2004. doi:10.1016/j.physletb.2003.10.112. 33
-
[45]
M. L. Ciccolini, S. Dittmaier, and M. Kramer. Electroweak radiative corrections to associated WH and ZH production at hadron colliders.Phys. Rev. D, 68:073003, 2003. doi:10.1103/PhysRevD.68.073003
-
[46]
John M. Campbell, R. Keith Ellis, and David L. Rainwater. Next-to-Leading Order QCD Predictions forW+ 2 Jet andZ+ 2 Jet Production at the CERN LHC. Phys. Rev. D, 68:094021, 2003. doi: 10.1103/PhysRevD.68.094021
-
[47]
Fernando Febres Cordero, Laura Reina, and Doreen Wackeroth.w- andz-boson production with a massive bottom-quark pair at the large hadron collider. Phys. Rev. D, 74:034007, 2006. doi: 10.1103/PhysRevD.74.034007
-
[48]
F. R. Anger, F. Febres Cordero, H. Ita, and V. Sotnikov. NLO QCD predictions forW b ¯bproduc- tion in association with up to three light jets at the LHC. Phys. Rev. D, 97(3):036018, 2018. doi: 10.1103/PhysRevD.97.036018
-
[49]
Associated production of a W boson and massive bottom quarks at next-to-next-to-leading order in QCD
Luca Buonocore, Simone Devoto, Stefan Kallweit, Javier Mazzitelli, Luca Rottoli, and Chiara Savoini. Associated production of a W boson and massive bottom quarks at next-to-next-to-leading order in QCD. Phys. Rev. D, 107(7):074032, 2023. doi:10.1103/PhysRevD.107.074032
-
[50]
J. Alwall, R. Frederix, S. Frixione, V. Hirschi, F. Maltoni, et al. The automated computation of tree- level and next-to-leading order differential cross sections, and their matching to parton shower simulations. JHEP, 1407:079, 2014. doi:10.1007/JHEP07(2014)079
-
[51]
TMVA: Toolkit for Multivariate Data Analysis
TMVA Collaboration. TMVA: Toolkit for Multivariate Data Analysis. 2007
work page 2007
-
[52]
R. Brun and F. Rademakers. ROOT: An object oriented data analysis framework.Nucl. Instrum. Meth. A, 389:81–86, 1997. doi:10.1016/S0168-9002(97)00048-X
-
[53]
M. Zerlauth. Hl-lhc operational scenarios and performance outlook. Presentation at Linear Collider Work- shop (LCWS 2025), 2025. URLhttps://agenda.linearcollider.org/event/10594/contributions/ 56953/attachments/40557/64535/LCWS-2025-(HL-)LHC-v0.2.pdf. Slides include integrated luminosity scenarios up to 4000 fb−1
work page 2025
-
[54]
Ties Behnke, James E. Brau, Brian Foster, Juan Fuster, Mike Harrison, James McEwan Paterson, Michael Peskin, Marcel Stanitzki, Nicholas Walker, and Hitoshi Yamamoto. The International Linear Collider Technical Design Report - Volume 1: Executive Summary. 2013
work page 2013
-
[55]
F. Garwood. Fiducial limits for the poisson distribution. Biometrika, 28(3/4):437–442, 1936. doi: 10.1093/biomet/28.3-4.437
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.