Taming a leading theoretical uncertainty in HH measurements via accurate simulations for bbH production

We present a new simulation for Higgs boson production in association with bottom quarks ( bbH ) at next-to-leading order (NLO) accuracy matched to parton showers in hadronic collisions. Both contributions, the standard one proportional to the bottom-quark Yukawa coupling and the loop-induced one proportional to the top-quark Yukawa coupling from the gluon-fusion process, are taken into account in a scheme with massive bottom quarks. Therefore, we provide the full simulation of the bbH final state in the Standard Model, which constitutes also a crucial background to measurements for Higgs-boson pair ( HH ) production at the Large Hadron Collider when at least one of the Higgs bosons decays to bottom quarks. So far, the modeling of the bbH final state induced one of the dominant theoretical uncertainties to HH measurements, as the gluon-fusion component was described only at the leading order (LO) with uncertainties of 𝒪 (100