NNLO QCD⊕QED corrections to Higgs production in bottom quark annihilation

We present next-to-next-to leading order (NNLO) quantum electrodynamics (QED) corrections to the production of the Higgs boson in bottom quark annihilation at the Large Hadron Collider (LHC) in the five flavor scheme. We have systematically included the NNLO corrections resulting from the interference of quantum chromodynamics (QCD) and QED interactions. We have investigated the infrared (IR) structure of the bottom quark form factor up to two-loop level in QED and in $\mathrm{QCD}\ifmmode\times\else\texttimes\fi{}\mathrm{QED}$ using $\mathrm{K}+\mathrm{G}$ equation. We find that the IR poles in the form factor are controlled by the universal cusp, collinear and soft anomalous dimensions. In addition, we derive the QED as well as $\mathrm{QCD}\ifmmode\times\else\texttimes\fi{}\mathrm{QED}$ contributions to soft distribution function as well as to the ultraviolet renormalization constant of the bottom Yukawa coupling up to second order in strong coupling and fine structure constant. Finally, we report our findings on the numerical impact of the NNLO results from QED and $\mathrm{QCD}\ifmmode\times\else\texttimes\fi{}\mathrm{QED}$ at the LHC energies taking into account the dominant NNLO QCD corrections.