Search for Higgs boson pair production in the $\mathrm{b}\overline{\mathrm{b}} \mathrm{l}\nu \mathrm{l}\nu$ final state at $\sqrt{s} = 13~\mathrm{TeV}$

A search for pair-produced Higgs bosons decaying respectively into $\mathrm{b}\overline{\mathrm{b}}$ and $\mathrm{V}\mathrm{V}$ (with $\mathrm{V}$ either a $\mathrm{W}$ or a $\mathrm{Z}$ boson), with subsequent $\mathrm{V}\mathrm{V}$ decays into two leptons and two neutrinos, is presented. The analysis is based on a sample of proton-proton collisions at $\sqrt{s} = 13~\mathrm{TeV}$ at the LHC corresponding to an integrated luminosity of $2.30~\mathrm{fb}^{-1}$. The search signature is a resonance in the invariant mass distribution of the $\mathrm{b}$-jet pair at the Higgs boson mass in combination with high scores of a boosted decision tree discriminant based on kinematic information. Data and predictions from the standard model (SM) are in agreement within uncertainties. For the SM $\mathrm{h}\mathrm{h}$ hypothesis, the data are observed (expected) to exclude a production cross-section times branching ratio of $166.7$ ($92.8^{+59.9}_{-33.4}$) fb, corresponding to approximately $400$ times the SM cross section. Lack of deviation from the SM predictions in the observations is used to place constraints on different scenarios considering anomalous couplings which could affect the rate and kinematics of $\mathrm{h}\mathrm{h}$ production.