Next Article in Journal
Dynamics and Exact Traveling Wave Solutions of the Sharma–Tasso–Olver–Burgers Equation
Next Article in Special Issue
Correction: Ferrari, A.; Rompotis, N. Exploration of Extended Higgs Sectors with Run-2 Proton–Proton Collision Data at the LHC. Symmetry 2021, 13, 2144
Previous Article in Journal
Deep-Learning-Based Surface Texture Feature Simulation for Surface Defect Inspection
Previous Article in Special Issue
Searching for New Physics in Hadronic Final States with Run 2 Proton–Proton Collision Data at the LHC
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Searching for Pairs of Higgs Bosons in the LHC Run 2 Dataset

by
Elizabeth Brost
1,† and
Luca Cadamuro
2,3,*,†
1
Brookhaven National Laboratory, Department of Physics, Upton, NY 11973, USA
2
IJCLab, Université Paris-Saclay, CNRS/IN2P3, 91405 Orsay, France
3
Department of Physics, Previously at University of Florida, Gainesville, FL 32611, USA
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Symmetry 2022, 14(7), 1467; https://doi.org/10.3390/sym14071467
Submission received: 1 April 2022 / Revised: 1 July 2022 / Accepted: 4 July 2022 / Published: 18 July 2022
(This article belongs to the Special Issue Recent Advance in Accelerator and Particle Physics)

Abstract

:
The discovery of the Higgs boson confirms the existence of a scalar sector of the standard model, responsible for electroweak symmetry breaking, but the nature and properties of the potential at the origin of this mechanism are still unknown. By studying the production of pairs of Higgs bosons (HH), physicists can directly measure the coupling of the Higgs boson to itself and thus determine the shape of this potential, which has far-reaching implications on the origin and evolution of our Universe. Because of this deep connection to the foundations of electroweak symmetry breaking, HH production is also an ideal place to search for manifestations of yet-unknown physics, such as modifications of the strength of the self-coupling and of the interaction between pairs of vector bosons and Higgs bosons. In this review article, we summarize the current searches for HH production at ATLAS and CMS, using the LHC Run 2 dataset, discuss the implications of our current constraints on physics beyond the standard model, and briefly review prospect for future HH searches.

The study of the production of Higgs boson pairs (HH) is the focus of a broad program of research developed at the CERN Large Hadron Collider (LHC) by the ATLAS [1] and CMS Collaborations [2]. This review aims at providing an overview of the motivations and the status of these searches at the LHC and of the prospects at future colliders. We introduce the topic in Section 1 by discussing the role of HH and of the Higgs boson self-coupling in our understanding of fundamental physics, followed by an overview of the HH production modes and decay channels in Section 1.2. The latest LHC Run 2 results are discussed in Section 2, with a particular focus on the recent developments that have improved the experimental sensitivity. Finally, we discuss the prospects for the measurement of HH in Section 3, and conclude with a summary in Section 4. Throughout this review, we tried to streamline the discussion and avoid the details of individual searches, which can be found in the corresponding papers referenced in the text, and rather focus on the general philosophy of the analyses. We hope that this choice will be a useful introduction to readers that want to learn about HH searches, and a comprehensive overview of the current results for experts on the topic.

1. Introduction

Characterizing the properties of the standard model of particle physics (SM) to deepen our knowledge of the fundamental mechanisms of Nature is the scientific focus of the experiments at the LHC. After the discovery of the Higgs boson (H) by the ATLAS and CMS Collaborations in 2012 [3,4,5], the focus of the high energy physics community has shifted toward the precise characterization of the properties of this particle. As the observation of the Higgs boson confirms the existence of a scalar sector of the SM, testing the mechanism at the origin of the existence of the Higgs boson itself is of primary importance. It is in this context that the study of the production of HH should be regarded, with its potential to shed light on key properties of this fundamental mechanism.
In this section, we briefly discuss the theoretical motivations to study HH as a tool to verify the validity of the SM description, to search for the existence of new physics beyond it, and to consider the implications of these measurements in our understanding of nature.

1.1. The Higgs Boson Self-Coupling

The Brout–Englert–Higgs (BEH) mechanism [6,7,8], which explains the quantum–relativistic nature of the mass of fundamental particles, postulates the existence of a complex scalar doublet of fields, Φ , that is subject to a potential of the form
V ( Φ Φ ) = μ 2 Φ Φ + λ ( Φ Φ ) 2 with μ 2 , λ > 0
which is minimal for non-zero field values defined by the relation
| Φ 2 | = μ 2 2 λ v 2 2
We refer to the constant v as the vacuum expectation value. When the symmetry of this potential is spontaneously broken to a ground state, the degrees of freedom of the field Φ give rise to the longitudinal polarizations of the W and Z bosons, and hence to their masses, and to a physical field H whose quantum corresponds to the Higgs boson, H. In particular, the resulting Lagrangian contains the following terms:
L 1 2 m H 2 H 2 + λ v H 3 + λ 4 H 4
meaning that the strengths of the trilinear ( H 3 ) and quadrilinear ( H 4 ) self-interactions are governed by the value of λ and thus directly depend on the shape of the scalar potential of Equation (1).
The possibility to probe the shape of the scalar potential through the direct measure of λ represents the key physics goal of the experimental study of HH production. The SM predicts the following relation for value of λ :
λ = m H 2 2 v 2 0.13
implying that once the measurement of the Higgs boson mass of m H = 125.25 ± 0.17 GeV [9] is taken into account, a precise prediction for the value of λ exists. HH measurements can thus provide a test of the consistency of the SM and a verification that the BEH mechanism is realized in nature.
Studying this interaction is not a simple verification of an as-yet unmeasured property of the SM. The interaction of the Higgs boson with fermions connects to the Yukawa coupling structure, and the interactions with vector bosons to the breaking of electroweak symmetry and the local quantization around the minimum of the scalar potential. In contrast, the self-interaction of the Higgs boson depends on a global shape property of the potential that is deeply connected to the foundations of the BEH mechanism. Measuring the Higgs boson self-coupling would be a confirmation that the scalar potential of Equation (1) does exist in Nature and would reveal fundamental information on its properties.
Owing to this deep connection with the fundamentals of the SM, the study of the Higgs boson self-coupling has implications in our understanding of the evolution of our Universe. The value of λ runs with the energy with a function that depends on the values of the top quark (t) and H masses, and an instability of the BEH potential minimum can occur if its value becomes negative at a certain scale [10], opening up the possibility of quantum tunneling from the current state to a new minimum. Figure 1 shows the running of λ once the measured t and H masses are taken into account, and the tunneling probability as function of their values. These measurements thus imply that, in the SM, our Universe is not stable but is in a state of meta-stability, where the timescale associated with the tunneling probability is comparable to or longer than the currently estimated age of the Universe. While there is not a one-to-one correspondence between the vacuum instability and the trilinear self-coupling, the experimental determination of λ could help shed light on whether the long-term existence of the electroweak vacuum is truly challenged in nature, or whether a new mechanism exists to ensure the stability of the Universe.
Other examples of the cosmological implications of the BEH potential are the role of the Higgs boson as the possible inflaton in the primordial Universe [11], or its role in primordial baryogenesis [12]. The study of HH production at the LHC and the determination of the Higgs self-coupling can deepen our understanding on these topics.

1.2. HH Production in Effective Field Theories

Many experimental and theoretical indications hint at the possibility that as-yet-undiscovered physics beyond the SM (BSM) exists in nature. Owing to its deep connection to the scalar potential and thus to the foundations of the scalar sector, the study of HH might reveal hints of that BSM physics.
Direct searches for new particles that decay to HH, responsible for resonant HH production, are an important area of direct searches for BSM physics that we will not cover in this review. For a discussion on the status of the LHC Run 2 resonant HH searches, see ref. [13].
Another possibility is that BSM physics manifests in the modification of the strength of the couplings involved in HH production. In particular, if such BSM physics is associated to an energy scale that is not directly reachable at the LHC, we could still observe its low-energy effects as modifications of the properties of HH production. The theoretical framework to study this case is provided by effective field theories (EFTs), where we expand the Lagrangian of the SM with higher-order operators that represent the low-energy effects induced by a complete theory at a higher energy scale Λ . EFTs provide a convenient theoretical framework to interpret the experimental results, potentially serving as a bridge to reinterpret them in the context of specific UV-complete models. The choice of the criterion to perform the expansion of the SM Lagrangian, and the choice of the symmetries that are imposed on the new theory and its fields, result in the definition of different EFTs.
Two main EFTs are currently discussed in the context of HH searches. The first is called the standard model EFT (SMEFT), where the Higgs field is considered to be a doublet in the theory, as in the usual SM description, and the theory is formulated in a S U ( 2 ) × U ( 1 ) invariant way, where the electroweak symmetry is broken by the scalar potential minimization. The SM Lagrangian is expanded with higher-order operators, whose contribution is suppressed by powers of 1 / Λ , giving
L = L SM + i c i Λ O i 5 + i c i Λ 2 O i 6 +
The values c i , that represent the “strength” of a specific operator, are referred to as Wilson coefficients. There is a single dimension-5 operator, and it violates the lepton number [14]. Consequently, it is neglected in the discussion of the Higgs (and HH) phenomenology, and the focus is on the dimension-6 operators, which provide the leading contribution in the 1 / Λ expansion.
The second EFT considered in HH is the Higgs EFT (HEFT), where the Higgs doublet is not present. In this case, the theory is invariant under U ( 1 ) and the electroweak symmetry is assumed to have been already broken. For the HH case, the HEFT results in a more generic phenomenology at the cost of a larger number of operators to be studied and constrained. More information of HEFT in the context of HH can be found in refs. [15,16].
In general, EFTs are a broad and active area of research that can encompass many types of experimental measurements beyond HH. We present in Section 1.3.2 an overview of the EFT phenomenology that can be studied in HH, without giving further details on the theoretical construct. A theoretical overview of EFTs can be found in ref. [17].

1.3. HH Production Modes

In proton–proton collisions at the LHC, Higgs bosons pairs can be produced through the mechanisms that are listed, together with the corresponding cross section in the SM, in Table 1. The experimental investigation so far focuses on the two main production modes, the fusion of gluons (ggF) and of vector bosons (VBF), whose corresponding Feynman diagrams at the lowest order are shown in Figure 2.

1.3.1. ggF Production with Anomalous Self-Coupling

The simplest approach to model BSM physics effects in HH consists of studying the case of self-coupling values that differ from the SM prediction. Deviations are parametrized in terms of the strength of the coupling with respect to the SM prediction as κ λ = λ / λ SM , such that the SM prediction corresponds to the case κ λ = 1 . Similarly, deviations of the strength of the Yukawa interaction y t between the top quark and the Higgs boson are parametrized as κ t = y t / y t SM .
From the Feynman diagrams of Figure 2a, it can be seen that the amplitude of ggF HH production can be written as
A κ λ κ t T + κ t 2 B
and the above relation holds in all orders in the QCD perturbative expansion. The symbols T and B are the amplitudes associated to the triangle and box diagrams at leading order (respectively on the left and right sides of Figure 2a), or the amplitudes associated to groups of diagrams with the same powers of κ λ and κ t at higher QCD orders. As a consequence, by computing the square of the amplitude of Equation (6), the HH cross-section dependence on the couplings takes the form
σ ( g g HH ) κ λ 2 κ t 2 T + κ λ κ t 3 I + κ t 4 B
with T = | T | 2 , B = | B | 2 , and I = | T * B + B * T | . It should be noted that a global factor κ t 4 can be taken out from the formula above, so that the cross-section shape (but not its magnitude) is fully described as a function of a single variable, κ λ / κ t . Therefore, it will be assumed κ t = 1 (as in the SM) in the following discussion.
The dependence of the total HH cross section on κ λ is shown, for ggF and the other production modes, in Figure 3 [27]. The ggF quadratic dependence has a minimum around κ λ 2.5 . It should be noted that the figure reports the dependence at NLO FTapprox for ggF, but this cross section is now known at NNLO with finite top quark mass effects [24], and this value is used by the full Run 2 analyses discussed in this review.
The T, B, and I distributions as a function of the HH invariant mass, m HH , are shown in Figure 4a. We use m HH since this is a key variable in experimental analyses, and the most important variable to represent the properties of the HH system—at leading order, the properties of the HH system can be fully described with only m HH , and the angle of one Higgs boson with respect to the beam line, which is approximately flat [31]. Higher values of m HH are correlated with a higher p T of the Higgs bosons and harder kinematics of the events. Anomalous values of κ λ alter the relative importance of these three components and thus the shape of the m HH distribution. In particular, the distribution in m HH is the softest for κ λ 5 and has a typical two-peak structure with a tail at high values for κ λ 2.5 , with strong effects on the shape for κ λ values between 0 and 7. In the SM, the m HH distribution is wide with a broad peak at m HH 400 GeV . These effects are illustrated in Figure 4b.
Experimental analyses use both the total cross section and the differential distributions in m HH to constrain the value of the self-coupling. At the current sensitivity, the capability to constrain the self-coupling largely stems from the enhanced HH cross section. However, the differential m HH information is important to develop analyses that are optimal over a broad range of κ λ values, and to solve the degeneracy between κ λ values that result in the same total cross section values.

1.3.2. ggF Production in EFTs

In the previous section, we implicitly assumed that BSM physics might manifest solely through the modification of the Higgs boson self-coupling, and that the interpretations of experimental searches in terms of constraints on κ λ assume that all the other couplings are fixed to the SM predictions. While this is a valid operative assumption if we want to determine how precisely we can measure the Higgs boson self-coupling given the available data, a more generic study of BSM physics effects requires the framework of EFTs introduced in Section 1.2.
The addition of new operators in the SMEFT and HEFT formalisms results in the existence of new effective interactions, and the associated Feynman diagrams are shown in Figure 5. In the SMEFT formalism, the method employed to determine the expansion results in the existence of an additional chromomagnetic operator that implies a t t g H vertex, but this is absent in the HEFT formalism. To date, this diagram is not simulated with the NLO-accurate Monte Carlo tools discussed below in Section 1.3.4, and therefore, it is not shown and further discussed here. Further information can be found in [32].
In addition to the t t H ( κ t ) and H H H ( κ λ ) interactions, we obtain in this formalism the t t H H ( c 2 ), g g H ( c g ) and g g H H ( c 2 g ) interactions, where the symbols in parentheses denote the strength of the coupling.
In the SMEFT formalism, there is a relation that binds the value of κ t to c 2 , and the values of c g to c 2 g , so that there are, in practice, three free parameters involved in HH production. Considering that much tighter constraints on κ t and c g are expected to be set in single Higgs boson measurements, this means essentially that the only parameter that HH measurements can help to constrain is κ λ .
Conversely, in the HEFT formalism, all five interactions are independent, and their effects should be all considered simultaneously when interpreting the HH results. The five couplings κ λ , κ t , c 2 , c g , and c 2 g correspond, respectively, to the coefficients of the relevant HEFT operators c h h h , c t t h , c t t h h , c g g h , and c g g h h , as they are denoted in ref. [33]. Since the exploration of a five-dimensional parameter space is computationally very challenging, experiments have thus far taken the practical approach of studying “shape benchmarks” that represent arbitrary combinations of ( κ λ , κ t , c 2 , c g , c 2 g ) that are associated to characteristic distributions of m HH . While these do not represent, per se, points of specific theoretical interest, they capture the main kinematic features that result from combinations of the couplings, guiding the development of analyses and allowing for an evaluation of their sensitivity to generic EFT scenarios. Two sets of benchmarks have been defined to date, one based on leading-order (LO) modeling of the signal [31,34,35] and a more recent one derived from the next-to-leading-order (NLO) simulation [36]. The values of these benchmarks are shown in Table 2.

1.3.3. VBF Production with Anomalous Couplings

As illustrated by the Feynman diagrams of Figure 2a, VBF HH production involves the VV H and VV H H couplings (with V = W , Z ). The strengths of these two interactions with respect to their SM predictions are represented by the coupling modifiers κ V and κ 2 V , respectively. In the SM, both interactions originate from the spontaneous breaking of electroweak symmetry, and their values are thus related to the same fundamental mechanism. As a consequence, a cancellation occurs in the amplitude for longitudinally polarized vector bosons. For generic κ V and κ 2 V couplings, at energies that are lower compared to the scale of the new physics, the amplitude depends on the coupling and on the energy s ^ as [37]
A ( V L V L HH ) s ^ v 2 κ 2 V κ V 2
In the SM, where κ V = κ 2 V = 1 , this amplitude is suppressed. However, if the Higgs boson is the manifestation of new strong dynamics at the TeV scale, such as in models with a composite Higgs boson [38,39], anomalous values of κ V or κ 2 V can result in a large enhancements of the VBF HH production rate at high m HH , as illustrated in Figure 6. The dependence of the total VBF HH cross section on the value of the anomalous couplings is shown in Figure 7 for various combinations of parameters. The effects of the cancellation from Equation (8) are clearly visible in the dependence of the cross section on κ V and κ 2 V . It should be noted that the independent modification of κ V and κ 2 V that we illustrate here is possible in a HEFT framework but not in SMEFT.

1.3.4. Monte Carlo Generators for HH

The most accurate simulation of ggF HH production is at the NLO precision with finite top quark mass effects [33,41,42], and includes the simulation of anomalous values of κ λ and of the EFT couplings in the EWChL framework. All the analyses based on the full Run 2 data set performed by the ATLAS and CMS Collaborations and discussed in this review use this simulation.
The VBF HH signal is simulated at the leading order precision with MadGraph, including effects from anomalous κ V , κ 2 V , and κ λ couplings from [37].

1.4. HH Decay Channels

Owing to the large number of decays of the Higgs boson, HH production gives rise to a phenomenologically rich set of final states, which are illustrated in Figure 8. When choosing a final state to be investigated in experimental searches, one Higgs boson is usually required to decay to either b b or W W to ensure a sizable branching fraction, while the choice of the decay of the second Higgs boson represents a trade-off between the abundance of the final state and the background contamination. The HH decay channels that have been chosen for searches in the LHC Run 2 data set by the ATLAS and CMS Collaborations are listed below with a very short overview of the background processes and the main challenges. The references to the published analyses of each channel, a more detailed discussion of the analysis methods, and an overview of the results are presented in Section 2.
  • b b ¯ b b ¯ ( B 34 % ) is characterized by the largest branching fraction. The experimental challenge is the large multijet background from QCD and t t ¯ production, that requires a high performance in the identification of b jets, efficient online identification at trigger level, and dedicated strategies for the modeling of the background and its separation from the HH signal. The current searches at the LHC cover the topology, where all four b jets are separately reconstructed (“resolved” topology), and the topology where the H decays at high Lorentz boost resulting in collimated decay products and the HH system is reconstructed with two large-radius jets (“boosted” topology).
  • b b ¯ τ τ + ( B 7.3 % ) is usually approached in τ τ final states with at least one decay to hadrons ( τ h ) and a neutrino, τ h τ h , τ μ τ h , and τ e τ h , that together collect more than 88% of the τ τ decays. The main irreducible backgrounds are t t ¯ and Z plus b jets production, plus instrumental backgrounds where hadronic jets are misidentified as τ and b jets. The kinematic information of the reconstructed objects is combined in a multivariate discriminant to identify the presence of a signal.
  • b b ¯ W ± W ( B 25 % ) is examined in the decays of the W W system containing either one (≈30%) or two (≈5%) light leptons (electrons and/or muons, while decays to taus are only indirectly accounted for in their subsequent decays to an electron or muon plus neutrinos). Presently, the results based on the complete Run 2 LHC data set only cover the two-lepton final state. The main background is t t ¯ , that is suppressed with the usage of multivariate discriminants.
  • b b ¯ Z Z * ( B 3.1 % ) is currently investigated in the Z Z 4 decay channel ( = e , μ , about 0.45% of Z Z decays) with the full Run 2 data set, although previous searches also investigate the 4 q and 2 q 2 ν decay channels of the Z Z system. The main backgrounds of the 4 search are Z, Z Z , and H production in association with jets. The analysis relies on the excellent m 4 resolution to suppress the Z and Z Z backgrounds and uses a boosted decision tree to reject single Higgs boson events.
  • b b ¯ γ γ ( B 0.26 % ) is characterized by a tiny branching fraction, but benefits from the low background contamination and the percent-level experimental resolution on the photon pair invariant mass to identify the signal. The main backgrounds are the continuum production of photons and jets, and single Higgs boson production in association with jets (ggF, t t H ). The latter can have an important impact on the sensitivity of the analysis since it is characterized by the same H γ γ signature, and is suppressed with dedicated multivariate techniques. The distribution of the diphoton invariant mass, m γ γ , is used to search for the presence of the signal, optionally with the simultaneous fit of the broader m b b distribution.
  • Decays of the HH system to W W W W ( B 4.6 % ), W W τ τ ( B 2.7 % ), and τ τ τ τ ( B 0.39 % ), where W and τ decays contain two, three, or four reconstructed leptons (e, μ , and τ h ), are collectively referred to as “multilepton” final state. Categories are defined based on the number, flavor, and electrical charge (same or opposite sign) of the leptons, and the main backgrounds are diboson production and processes with misidentified leptons. Boosted decision trees are used to enhance the presence of the signal in the selected data samples.

1.5. Overview of Earlier HH Results

The investigation of HH production at the LHC has been marked by spectacular progress over the past years. If the study of HH was initially regarded as difficult, even at the HL-LHC, because of its low cross section, the efforts of the ATLAS and CMS analysis teams have resulted in a large improvement of the experimental sensitivities. Channels such as b b ¯ b b ¯ that were originally thought to be impossible because of the large backgrounds, provide some of the most stringent constraints on HH production today. Several analyses investigate both the ggF and VBF production modes to probe the strength of both the H H H and the VV H H vertices. As this review was prepared, the results based on the full LHC Run 2 data set (about 140 fb 1 ) in many channels have been made public, and the experiments are working toward their combination. We will cover these latest results and discuss the main drivers of analysis improvements in Section 2.
It is instructive to quickly recap the previous set of results that correspond to the combination of analyses based on the partial Run 2 data set collected in 2015 and 2016 (about 36 fb 1 ). A summary of the 95% confidence level (CL) limit on the signal strength, defined as the ratio of the excluded cross section with respect to the SM prediction, is shown in Figure 9 for the ATLAS and CMS HH combinations [43,44].
These results show that there is not a single “golden” channel for the investigation of HH, and the final result benefits from the combination of several analyses—in particular, from the combination of b b ¯ γ γ , b b ¯ τ τ + , and b b ¯ b b ¯ . A summary these results is presented in Table 3. The different hierarchy between the ATLAS and CMS results per channel was traced to different analysis methodologies, and the choice to focus more effort on the SM or BSM optimizations.
The upper limits on the HH cross section as function of κ λ are illustrated in Figure 10 for the two experiments. The shape of the limit curves directly relates to the modification of the interference pattern as discussed in Section 1.3.1. For very large | κ λ | values, the squared amplitude of the triangle diagram dominates, while in the region between κ λ ∼0 and 7, the changes induced in the shape of m HH result in the large changes of the experimental sensitivity. The sensitivity is maximal when the spectrum is hardest ( κ λ 2.5 ) and minimal when the spectrum is softer ( κ λ 5 ). Since most of the cross section enhancement effects appear in the region close to the kinematic production threshold of m HH = 250 GeV , channels with a good acceptance at low m HH , such as b b ¯ γ γ , have a smaller gap in the sensitivities between these two extreme regimes when compared to channels such as b b ¯ b b ¯ , where the higher thresholds imposed on the final state objects result in a lower acceptance at low m HH .
When comparing these early results with the most recent ones illustrated in Section 2, it should be noted that neither the contribution from the VBF HH signal nor the latest uncertainty related to the top quark mass scheme are included here.

2. Summary of LHC Run 2 HH Results

2.1. Updated HH Searches with Full LHC Run 2 Data Set

ATLAS and CMS are currently pursuing the full Run 2 HH analysis program. Several results have already been updated using the full Run 2 integrated luminosity. Updates on previous channels, including the ‘big three” b b ¯ b b ¯ , b b ¯ τ τ + , and b b ¯ γ γ that lead the sensitivity, were accompanied by the analysis of new HH channels, where branching ratios are smaller or where the event topology requires additional creativity to benefit of the Run 2 data set. We present in what follows an overview of the analysis techniques and discuss the results at the end of the Section.

2.1.1. HH→ b b ¯ b b ¯

The b b ¯ b b ¯ channel benefits from the largest branching ratio, but must contend with the large multijet backgrounds and combinatorial challenges when reconstructing the two Higgs bosons. The precise modeling and the efficient rejection of the background are the key elements to maximize the sensitivity of the analyses. The CMS Collaboration has performed analyses of the full Run 2 data set in the b b ¯ b b ¯ channel, investigating both the ggF and the VBF production modes, in resolved and boosted regimes. The resolved analysis covers the bulk of the HH signal, while the boosted focuses on a specific corner of the phase space, where the background contamination is minimized, and a powerful identification algorithm for b b large radius jets can be deployed. The ATLAS Collaboration has published to date the study of the VBF production mode in the resolved topology with the full Run 2 data set.

Resolved HH → b b ¯ b b ¯

Trigger represents a fundamental challenge in b b ¯ b b ¯ searches, and both the ATLAS and CMS analyses [55,56] rely on online algorithms to trigger on events containing at least two (ATLAS) or three (CMS) b jets, achieving thresholds as low as 30 or 40 GeV on the softest jet. Additional forward jets in opposite η regions of the ATLAS or CMS detector are required to tag the VBF topology. In the ATLAS analysis, the VBF candidate jets are chosen by requiring | Δ η j j | > 5.0 and m j j > 1000 GeV , while in the CMS analysis, no selection on m j j is applied, but a boosted decision tree (BDT) discriminant is trained to separate the ggF and the VBF production modes.
The reconstruction of the two H candidates is performed by both analyses by choosing the pairing that minimizes the distance from a line in the ( m b b 1 , m b b 2 ) plane that connects the origin to the expected peak position of the HH signal. In the CMS case, the choice that maximizes the p T of the H candidates is additionally performed when the jet resolution limits the performance of this decision. As this procedure effectively minimizes the invariant mass difference between the two pairings but avoids directly imposing the m H = 125 GeV hypothesis, the combinatoric challenge is solved without sculpting the multijet background to look like the HH signal. Optimal resolution on b jets is ensured through the use of a regression algorithm to estimate the neutrino energy loss from B hadronization and decays.
Signal regions are defined by the presence of 4 b jets and by requiring the masses of the two H candidates to be compatible with the Higgs boson. The CMS analysis performs an additional categorization of the events depending on their properties (low and high m HH for ggF, SM and BSM kinematics for VBF), providing maximal sensitivity to both the SM production and anomalous couplings.
The ATLAS analysis uses the m HH distribution to search for the presence of the signal, as shown in Figure 11. In the CMS analysis, either m HH or a counting experiment are used in the VBF categories, depending on the expected total number of events. For the ggF case, the output of a BDT that combines several kinematic variables is used, looking for the signal as an enhancement at high scores of this discriminant. An example of the distributions of these variables are shown in Figure 12.
The accurate modeling of the background calls for data-driven methods, and the multijet processes are estimated from background-enriched regions where the b tagging requirements and the H mass requirements are inverted. In the ATLAS VBF analysis, an iterative procedure is implemented to progressively correct the background template, following the method developed for the partial Run 2 data set analysis [45]. In the CMS analysis, a new method was deployed. It uses a BDT with a dedicated metric to learn the differences between 3 b and 4 b data in a mass sideband region and to correct the control region data, enabling the simultaneous modeling of the variables that are used in the final discriminant.

Boosted HH → b b ¯ b b ¯

The boosted search performed by the CMS Collaboration [57,58] studies both the ggF and VBF signal. For the VBF case, the focus of the analysis is on anomalous κ 2 V couplings, which results in the hard m HH distributions shown in Figure 6. The two Higgs bosons are reconstructed as large-radius jets with p T above 400 and 500 GeV , and two standard radius jets are additionally required to reconstruct the outgoing VBF partons. For the ggF search, where the signal is expected to be softer, lower thresholds of 300 GeV are used for the jets. The ParticleNet algorithm [59,60], based on a powerful graph neural network, is used to identify the H b b ¯ decays and to improve the resolution on the jet mass, benefiting the definition of signal regions based on m H . Categories for different signal purity are defined based on the ParticleNet output to maximize the sensitivity. A combination of simulation for hadronic t t ¯ and data-driven methods (for QCD multijet) are employed to estimate the background from sidebands defined by inverted ParticleNet and Higgs boson mass requirements. A summary of the expected number of events in the boosted searches is shown in Figure 13.
Although the orthogonality with the resolved b b ¯ b b ¯ CMS analysis is not yet enforced in these results, the overlap is expected to be small since the boosted selections define a small portion of the total phase space, which could be removed from the events selected in the resolved analysis with a small impact on the sensitivity. The two topologies are thus expected to contribute equally to the global sensitivity.

2.1.2. HH→ b b ¯ τ τ +

The b b ¯ τ τ + channel trades off between the large branching ratio for Higgs bosons decaying to bottom quarks, and the stronger background rejection when selecting two tau leptons. Both Higgs bosons can be reconstructed, and their masses can be used to discriminate between the signal and background.
Both ATLAS [61] and CMS [62] have updated their b b ¯ τ τ + searches to the full Run 2 data set. Events are categorized by the tau decays ( τ h τ h , τ μ τ h , and τ e τ h ), and by the triggers used to select events (ATLAS), or by the number of b-tagged jets and by resolved, boosted, or VBF topology (CMS). Machine learning plays a big role in these analyses, from the b jet and tau tagging, to H b b candidate tagging (CMS), to signal and background discrimination. The latter combines the kinematic properties of the reconstructed final state objects to enhance the presence of signal at high values of the discriminant. Background modeling is a key challenge, and MC simulation is supplemented with dedicated control regions to constrain the normalization of key backgrounds involving Z bosons and top quarks. Signal extraction is performed in the MVA output in each category. Figure 14 shows the BDT output score in the ATLAS τ h a d τ h a d category and aggregated DNN output score, sorted by the expected signal-to-background ratio, for all the CMS categories combined.

2.1.3. HH→ b b ¯ γ γ

The b b ¯ γ γ decay channel benefits from the excellent m γ γ resolution for the extraction of results, as well as the efficient and relatively low-threshold diphoton triggers used by both experiments. A large challenge in the last analysis round was the low statistics in data, once a simple signal selection was made. The size of the full Run 2 data set allowed for more creativity in analysis design for this round.
Both ATLAS [63] and CMS [64] updated their b b ¯ γ γ searches to the full Run 2 data set. Events with at least two photons with p T 1 > 0.35 ( 0.33 ) × m γ γ and p T 2 > 0.25 × m γ γ and exactly two (at least two) b-tagged jets are selected with p T > 25 GeV by ATLAS (CMS). Both searches use advanced machine learning techniques to separate the tiny HH→ b b ¯ γ γ signal from the large background and from the contamination of single Higgs processes (mostly t t H ) that are characterized by the same H γ γ signature. A modified four-body-mass variable, M*, defined as M* = m γ γ j j ( m γ γ m H ) ( m j j m H ) is used to reduce dependence on the jet and photon energy resolutions. This is true in the case that the photons and jets originate from the signal process. Both searches take advantage of the good acceptance for low-p T photons by dividing the M* spectrum into a number of categories (ATLAS divides the M* spectrum in half and then creates two additional categories based on BDT score, CMS creates two M* categories for the VBF selection, and four M* categories for the ggF selection, before splitting into 3 further categories based on MVA score), for a total of 4 ATLAS categories and 14 CMS categories. Then, the m γ γ spectrum is used to extract the final results: ATLAS fits the m γ γ spectrum in each category, while CMS performs a simultaneous fit in m γ γ × m j j in each category. The fit variables are shown in Figure 15, and the results are presented at the end of this section.

2.1.4. Final States with Leptons

The decay channels involving light leptons in the final state, including b b ¯ VV * ( ν ν ), b b ¯ Z Z * ( ), W ± W W ± W , W ± W τ τ + , and τ τ + τ τ + bring unique challenges for signal and background separation, as well as HH event reconstruction. These channels often benefit from low light lepton ( e , μ ) trigger thresholds.
Multilepton: The CMS multilepton search [65] targets HH→ W ± W W ± W , W ± W τ τ + , and τ τ + τ τ + decays. Together, these channels make up 7.7% of all HH decays, so this could be an important place to gain sensitivity to SM HH production. This analysis divides the signal region into seven categories based on the number of light leptons and/or hadronic taus. The categories and their sensitivity to SM HH production can be seen in Figure 16a. This plot demonstrates the power of applying a common analysis technique to selections that are not individually very powerful in the overall HH combination, but together, they can set an impressive limit.
b b ¯ Z Z * (4): The CMS bbZZ(4) search [66] is performed in events with at least four leptons, forming two pairs of opposite sign, same flavor leptons, and at least two jets with p T > 20 GeV . The leading lepton is required to have p T > 20 GeV , and the subleading lepton is required to have p T > 12 (10) GeV for electrons (muons). No other requirements are placed on the jets, and the two jets with highest b-tagging score are chosen as the H b b candidate jets. A BDT is trained separately in nine categories based on lepton flavor ( 4 e , 2 e 2 μ , 4 μ ) and data-taking year (2016–2018). The combined BDT output score for all categories is shown in Figure 16b.
bb ν ν : The ATLAS bb ν ν search [67] starts by selecting events with exactly two oppositely charged light leptons (e, μ ) and exactly two b-tagged jets. A DNN classifier is constructed using kinematic information from each of the two Higgs candidates and their relative position in the detector, and trained on the dominant HH→ b b ¯ W ± W process. The signal region is split into events with same flavor lepton pairs and those with different flavor, which allows for better separate between the signal and the background from Z boson production. Figure 16c shows the comparison between the data and the background modeling, for the DNN discriminant variable, d H H .

2.1.5. Results and Summary

Table 4 presents the limits set on σ / σ S M , a comparison to the earlier Run 2 limits, and the limits set on the self-coupling modifier, κ λ , in each analysis. Figure 17 summarizes the limits on the HH production cross section compared to its SM value, for each of the analyses presented in this section.
Limits on the VBF SM cross section for the various channels are summarized in Table 5. The excluded cross section as a function of κ 2 V is shown in Figure 18 for the boosted b b ¯ b b ¯ CMS analysis, that provides the strongest constraint. Results are shown under the variation of κ 2 V only, and with the simultaneous variation of κ V as well. The excluded contours connect to the dependence of the theoretical cross section on these couplings of Figure 7c.
Limits on the shape benchmarks discussed in Section 1.3.2 are illustrated in Figure 19 for the CMS b b ¯ γ γ search (LO benchmarks) and the ATLAS combination of b b ¯ γ γ and b b ¯ τ τ + (NLO benchmarks). The scatter of the excluded cross sections, that vary by more than one order of magnitude, illustrates the sensitivity of the analyses to the kinematic effects, and shows that constraints on the SM cannot be simply ported to generic EFT sensitivities. As we are approaching the SM sensitivity for HH production, a systematic study of the EFT effects and of the contributions of the operators beyond simple benchmarks will help in extracting the most information from the experimental measurements.

2.2. Preliminary Run 2 Combination

The ATLAS Collaboration performed a preliminary combination [69] of the recent b b ¯ γ γ [63] and b b ¯ τ τ + [61] HH results. The observed (expected) combined 95% CL upper limit is found to be 3.1 (3.1) times the standard model prediction, which represents a factor of 3.2 improvement over the previous ATLAS HH combination, even though this preliminary combination only uses two channels. The observed (expected) allowed range for the self-coupling modifier, κ λ , is found to be −1.0 < κ λ < 6.6 (−1.2 < κ λ < 7.2), at 95% CL.
Both the b b ¯ γ γ and the b b ¯ τ τ + channels are limited by statistical uncertainties, and their dominant uncertainties are mostly independent between the two analyses. As a consequence, the combined result is close to a simple sum in the quadrature of the individual sensitivities, as shown in Figure 20. Since these assumptions hold for all the HH channels summarized in Table 4 and Figure 17, we can perform the same sum for all the channels and estimate that a sensitivity of the order of 2–3 times the SM prediction will be achieved by each experiment at the LHC Run 2.

2.3. Lessons Learned from the Run 2 HH Programme

Huge progress has already been made toward the eventual measurement of the Higgs pair production cross section. This progress can be attributed to the improvements in online and offline event selection, object reconstruction, and creative analysis techniques. These improvements have allowed the ATLAS and CMS collaborations to produce results that outpace the improvements expected simply from the addition of more data. Two elements that proved to be central to achieve these results are as follows:
HH event kinematics: understanding and exploring the unique kinematics of HH events with varied anomalous couplings has allowed this round of analyses to achieve unprecedented limits. As shown in Figure 4, varying κ λ from its SM value results in softer event kinematics. Following early Run 2 analyses, we used this information systematically in all channels, using the m HH information in the analyses and designing the selections to also target low- m HH events for κ λ extraction. This in turn highlights the importance of being able to select low- p T objects in the trigger, which has traditionally been a challenge for hadronic signatures because of the large multi-jet background at the LHC. The latter point is one that we should keep in mind when preparing for Run 3 and HL-LHC operations.
Machine learning was used to great effect in all of the HH analyses presented in this review. It contributes to these HH searches at all levels, from object identification and selection to separating the small HH signal process from large backgrounds. The importance of high-performance object identification is illustrated by the b b ¯ τ τ + and b b ¯ b b ¯ analyses, where improved neural network-based algorithms have boosted the sensitivities.

3. Prospects for Future Measurements

3.1. LHC Run 3

The LHC plans to increase the center-of-mass energy from s = 13 TeV to s = 13.6 TeV for the Run 3, which will start in mid-2022. The larger center-of-mass energy will result in a ≈20% larger HH production cross section for gluon–gluon fusion (ggF). LHC Run 3 will produce approximately 400 fb 1 of data per experiment, which will be especially beneficial to the HH searches, which are currently statistically limited.
Both ATLAS and CMS have finished “Phase-I” upgrades to the detectors, which will improve the quality of data that can be recorded during Run 3. The ATLAS Phase-I upgrade involves the installation of new muon detectors [70], which will improve muon identification and trigger capabilities, an upgrade of the liquid argon calorimeter readout electronics, which will provide higher granularity information to the trigger [71], and and upgrade of the trigger and data acquisition system to make use of the upgraded detector capabilities [72]. The CMS Phase-I upgrades, partly already operating during Run 2, included an improved pixel detector [73], an upgraded hadron calorimeter to cope with the higher pileup environment [74], new endcap muon detectors [75] and a replacement of the calorimeter, muon, and global trigger electronics [76], that will continue to assure efficient online event selection during Run 3.
Following the considerations highlighted in Section 2.2, we can expect that the Run 2 combined limit will improve by at least 400 / 140 1.7 . A combination of the results of ATLAS and CMS should bring roughly a factor 2 of improvement, meaning that a Run 3 LHC sensitivity at the level of the SM predictions is within reach. If the impressive pace of improvements discussed in Section 2 is maintained, a hint or an evidence of HH production is in reach during LHC Run 3.

3.2. HH Prospects at the High-Luminosity LHC

The full LHC physics program, expected to extend until ∼2040, will culminate in the construction and operation of the high-luminosity LHC (HL-LHC), which is an ambitious upgrade to the accelerator complex [77]. The HL-LHC will collide protons at s = 14 TeV , and plans to provide 3000 fb 1 of data to each of the ATLAS and CMS experiments. The data-taking conditions at the HL-LHC will be challenging: the instantaneous luminosity will increase to 5– 7 × 10 34 cm 2 s 1 , and the number of simultaneous collisions will increase to 140–200 per event.
The Phase-II upgrades of the ATLAS and CMS detectors [78,79] will facilitate data taking and event reconstruction in the much higher pileup environment of the HL-LHC. In order to maintain or lower the trigger thresholds with respect to those used in Run 2, the trigger and data acquisition systems will be substantially upgraded or replaced. New tracker systems will be installed, extending the tracking coverage up to about | η | = 4. The calorimeters will be upgraded with new electronics or replaced with new technologies in the challenging endcap regions. Detectors to precisely measure the time of arrival of particles will help to face the challenges of pileup. Finally, the muon systems will be upgraded with new electronics and additional muon chambers.
Early HL-LHC projections [80], based on extrapolations from early Run 2 HH analyses or on dedicated analyses of HL-LHC conditions with a parametric simulation, suggest that ATLAS and CMS would achieve together a combined significance of 4 σ , including the b b ¯ γ γ , b b ¯ τ τ + , b b ¯ b b ¯ , b b ¯ Z Z * (4), and b b ¯ VV * ( ν ν ) channels. The projected precision on the value of the self-coupling is 50%. Figure 21 presents the negative log-likelihood scan as a function of κ λ , for each of the individual HH channels considered and their combination. This corresponds to an exclusion range of 0.1 < κ λ < 2.3 at 95% CL.
A small number of individual channel HL-LHC projections was recently updated, and they are summarized here:
  • HH b b ¯ γ γ + b b ¯ τ τ + : The recent preliminary combination of the b b ¯ γ γ and b b ¯ τ τ + channels by ATLAS (described in Section 2 is projected to the HL-LHC luminosity, and results in a combined significance of 3.2 σ , assuming the baseline systematic uncertainty scenario described in the reference [81]. This result also excludes values of κ λ outside the range 0.0 < κ λ < 2.7 .
  • HH in final states with photons: CMS prepared HH projection in three channels, in final states including photons. These studies were performed using dedicated simulations of the HL-LHC conditions (so they are not projections of earlier published work), and using the generic detector delphes [82]. In the b b ¯ γ γ channel [83], a boosted decision tree is used to separate the signal and background, and the extraction of the signal significance is performed using a simultaneous fit to the m γ γ and m b b distributions. The significance is found to be 2.16 σ (to be compared to 1.8 σ in a previous projection). In the W ± W γ γ and τ τ + γ γ channels [84], events are categorized by the number of light leptons or taus, and then the signal extraction is performed in each category using the m γ γ distribution. The resulting significance is 0.22 σ , when combining all categories.
  • CMS ttHH (4b): The cross section for ttHH production in the standard model is incredibly small, at just 0.948 fb [85]. In this search [86], events with exactly one light lepton and at least four jets are fed into a 2-step deep neural network to first separate ttHH events from events containing top quarks and Z bosons, and then to separate the SM ttHH signal from all background processes, for the various jet and b-jet multiplicities. The statistical analysis is performed using the simultaneously fitting the DNN discriminants, and an upper limit is set on SM ttHH production of 3.14 times the SM cross section.
We expect the exciting trend of improvements due to creative analysis techniques to continue, and look forward to seeing further HH projections for the HL-LHC as they are prepared.

3.3. Future Colliders

While we expect that improvements to data-collection, reconstruction, and analysis techniques will improve enough such that observation of SM Higgs pair production will be possible at the HL-LHC, it is unlikely that the Higgs self-coupling will be measured at comparable precision to other Higgs couplings that have been studied at the LHC. A precision of 5–10% would give us sensitivity to a broad class of new particles that could couple to the Higgs boson, and a precision of 1% would be sensitive to quantum corrections to the value of the self-coupling [87].
Projections of the sensitivity to the Higgs self-coupling were performed for several proposed future colliders:
  • Low-energy e + e colliders such as the FCC-ee, CEPC, or the proposed low-energy ILC runs, are able to make indirect measurements of the Higgs self-coupling through higher-order corrections to single-Higgs processes. These machines will be able to measure the self-coupling with a precision of 40% [88].
  • Higher-energy e + e colliders such as CLIC or higher-energy ILC runs would be able to produce pairs of Higgs bosons, and thus have direct access to the self-coupling, which would allow an O(20%) determination of the Higgs self-coupling [88].
  • A μ + μ collider operating at s = 10 ( 30 ) TeV could achieve 5.6 (2.0)% precision on the self-coupling [89], studying pairs of Higgs bosons produced via vector boson fusion.
  • Very high energy p p colliders, such as the proposed 100 TeV FCC-hh, would also have the ability to produce pairs of Higgs bosons. A study [90] combining the b b ¯ γ γ , b b ¯ τ τ + , and b b ¯ b b ¯ decay modes projects a precision of 3.4–7.8% on the value of the self-coupling.

4. Conclusions

The study of HH production connects deeply to the foundations of the SM. This process gives direct access to the determination of the Higgs boson self-coupling, and thus to the shape of the scalar potential postulated in the Brout–Englert–Higgs mechanism, with deep implications on our understanding on the foundations of particle physics and on the cosmological history of our Universe. The measurement of HH provides at the same time a test of the consistency of the SM and a probe of the possible existence of as-yet-undiscovered physics beyond it.
The ATLAS and CMS Collaborations at the CERN LHC are pursuing a broad program of investigation of HH production. As this manuscript is being prepared, most of the analyses based on the full Run 2 p p data set, recorded in 2016–2018 at a center-of-mass energy of 13 TeV and corresponding to an integrated luminosity L 140 fb 1 , have been completed, and their statistical combination is being prepared. The most sensitive individual channels achieve a sensitivity in the Run 2 data set ranging between 4 and 8 times the SM prediction, meaning that their combination is expected to probe signals that are 2–3 times the SM prediction. The analyses currently study both the ggF and the VBF HH production modes, and set constraints on the strength of the self-coupling, κ λ , and on the VV HH interaction, κ 2 V . The current best constraint of κ λ is between −1 and 6.6 from the combination of the b b ¯ γ γ and b b ¯ τ τ + analyses, and will further improve, as the full set of channels will be combined. The best constraints of κ 2 V from the b b ¯ b b ¯ channel in a search for high Lorentz boost Higgs bosons and corresponds to the range 0.6 < κ 2 V < 1.4 , and excludes the hypothesis that the VV HH interaction does not exist when we assume that only the SM diagrams participate in VBF HH production.
When compared to the previous set of results, based on the partial Run 2 data set with L 36 fb 1 , the latest analyses show impressive improvements that are factors larger than the simple improvement expected from the increase in the size of the data set that is roughly proportional to 1 / L , i.e., 140 / 36 2 times better. These results were enabled by the improvements in the analysis techniques and in the identification of physics objects, often leveraging on new machine learning techniques. The larger data set and the growing interest for HH also enabled the study of new channels and of new topologies of the final states that contribute to boost the total sensitivity.
Such considerations set the ground for excellent prospects for the LHC Run 3, where a further increase in the total data set size to about 400 fb 1 is expected after 3 years of p p run in 2023–2025 at a center-of-mass energy of 13.6 TeV . With the sole increase in the luminosity and the combination of the results of ATLAS and CMS, a global LHC sensitivity at the level of the SM prediction is within reach. Further progress in the object identification methods and analysis techniques, and a careful preparation of the Run 3 data taking with more efficient triggers may bring a first hint or an evidence of the existence of HH production in only a few years from now.
At the HL-LHC, where a data set of about 3000 fb 1 of s = 14 TeV proton–proton collisions is expected to be collected by ATLAS and CMS, the ultimate LHC sensitivity to HH and to the self-coupling will be achieved. The prospects for these measurements were comprehensively studied in the context of the update of the European Strategy for Particle Physics submitted in 2018. These results show that a combination of the ATLAS and CMS HH analyses is expected to yield a significance of 4 σ for the SM signal and to constrain the value of κ λ to a precision of about 50% at the 68% confidence level. It should be noted, however, that these projections are based on the extrapolation or on the reproduction with a simplified simulation of the early Run 2 analyses based on the 2015 and 2016 data sets, and therefore mostly account for the increase in the integrated luminosity and the performance of the upgraded detectors. The improvements already observed in the timescale of less than five years with the most recent results show that larger data sets enable more sophisticated and powerful analyses, and suggest that these projections may be conservative, as supported by recent projections in selected final states.
The study of HH at future e + e colliders will require high center-of-mass energies that are achievable only in the latest stages of the linear machines (CLIC, ILC), while at lower energies and at circular colliders (FCC-ee, CEPC), indirect constraints on the self-coupling can be set from the study of single Higgs bosons. New machines, such as a μ + μ collider, might achieve the energy for the study of HH. Finally, the 100 TeV p p collider, FCC-hh, that is being discussed as the long-term CERN collider project has the capability to measure κ λ with a precision of ∼5%, providing an accurate map of the shape of the scalar potential.
There is no doubt that the study of HH is an essential part of the short-, mid-, and long-term programs of the research in particle physics, and that its study might hold the key to answer many open questions of high energy physics. We are currently in a very favorable moment of this program: we are about to fully exploit the Run 2 data set, getting ready to unlock the potential of LHC Run 3, and preparing for the challenges and opportunities of the HL-LHC. We are looking forward to everything that the study of HH production will reveal.

Author Contributions

E.B. and L.C. contributed equally to this work. All authors have read and agreed to the published version of the manuscript.

Funding

The work of E. Brost is supported by the U.S. Department of Energy, Office of Science, Office of High Energy Physics.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank Ramona Gröber and Javier Mazzitelli for the helpful discussions on the theoretical aspects of HH discussed in this document.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. ATLAS Collaboration. The ATLAS Experiment at the CERN Large Hadron Collider. JINST 2008, 3, S08003. [Google Scholar]
  2. CMS Collaboration. The CMS Experiment at the CERN LHC. JINST 2008, 3, S08004. [Google Scholar]
  3. ATLAS Collaboration. Observation of a new particle in the search for the Standard Model Higgs boson with the ATLAS detector at the LHC. Phys. Lett. B 2012, 716, 1–29. [Google Scholar] [CrossRef]
  4. CMS Collaboration. Observation of a New Boson at a Mass of 125 GeV with the CMS Experiment at the LHC. Phys. Lett. B 2012, 716, 30–61. [Google Scholar] [CrossRef]
  5. CMS Collaboration. Observation of a New Boson with Mass Near 125 GeV in pp Collisions at s = 7 and 8 TeV. J. High Energy Phys. 2013, 6, 081. [Google Scholar]
  6. Englert, F.; Brout, R. Broken Symmetry and the Mass of Gauge Vector Mesons. Phys. Rev. Lett. 1964, 13, 321. [Google Scholar] [CrossRef] [Green Version]
  7. Higgs, P.W. Broken Symmetries and the Masses of Gauge Bosons. Phys. Rev. Lett. 1964, 13, 508. [Google Scholar] [CrossRef] [Green Version]
  8. Guralnik, G.S.; Hagen, C.R.; Kibble, T.W.B. Global Conservation Laws and Massless Particles. Phys. Rev. Lett. 1964, 13, 585. [Google Scholar] [CrossRef] [Green Version]
  9. Particle Data Group. Review of Particle Physics. PTEP 2020, 2020, 083C01. [Google Scholar]
  10. Elias-Miro, J.; Espinosa, J.R.; Giudice, G.F.; Isidori, G.; Riotto, A.; Strumia, A. Higgs mass implications on the stability of the electroweak vacuum. Phys. Lett. B 2012, 709, 222–228. [Google Scholar] [CrossRef] [Green Version]
  11. Bezrukov, F.L.; Shaposhnikov, M. The Standard Model Higgs boson as the inflaton. Phys. Lett. B 2008, 659, 703–706. [Google Scholar] [CrossRef] [Green Version]
  12. Morrissey, D.E.; Ramsey-Musolf, M.J. Electroweak baryogenesis. New J. Phys. 2012, 14, 125003. [Google Scholar] [CrossRef]
  13. Veatch, J. Searches for Resonant Scalar Boson Pair Production Using Run 2 LHC Proton-Proton Collision Data. Symmetry 2022, 14, 260. [Google Scholar] [CrossRef]
  14. Jenkins, E.E.; Manohar, A.V.; Stoffer, P. Low-Energy Effective Field Theory below the Electroweak Scale: Operators and Matching. J. High Energy Phys. 2018, 3, 016. [Google Scholar] [CrossRef] [Green Version]
  15. Grober, R.; Muhlleitner, M.; Spira, M.; Streicher, J. NLO QCD Corrections to Higgs Pair Production including Dimension-6 Operators. J. High Energy Phys. 2015, 9, 092. [Google Scholar] [CrossRef]
  16. Buchalla, G.; Capozi, M.; Celis, A.; Heinrich, G.; Scyboz, L. Higgs boson pair production in non-linear Effective Field Theory with full mt-dependence at NLO QCD. J. High Energy Phys. 2018, 9, 057. [Google Scholar] [CrossRef] [Green Version]
  17. Brivio, I.; Trott, M. The Standard Model as an Effective Field Theory. Phys. Rep. 2019, 793, 1–98. [Google Scholar] [CrossRef] [Green Version]
  18. Dawson, S.; Dittmaier, S.; Spira, M. Neutral Higgs boson pair production at hadron colliders: QCD corrections. Phys. Rev. D 1998, 58, 115012. [Google Scholar] [CrossRef] [Green Version]
  19. Borowka, S.; Greiner, N.; Heinrich, G.; Jones, S.P.; Kerner, M.; Schlenk, J.; Schubert, U.; Zirke, T. Higgs Boson Pair Production in Gluon Fusion at Next-to-Leading Order with Full Top-Quark Mass Dependence. Phys. Rev. Lett. 2016, 117, 012001. [Google Scholar] [CrossRef] [Green Version]
  20. Baglio, J.; Campanario, F.; Glaus, S.; Mühlleitner, M.; Spira, M.; Streicher, J. Gluon fusion into Higgs pairs at NLO QCD and the top mass scheme. Eur. Phys. J. C 2019, 79, 459. [Google Scholar] [CrossRef]
  21. de Florian, D.; Mazzitelli, J. Higgs Boson Pair Production at Next-to-Next-to-Leading Order in QCD. Phys. Rev. Lett. 2013, 111, 201801. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Shao, D.Y.; Li, C.S.; Li, H.T.; Wang, J. Threshold resummation effects in Higgs boson pair production at the LHC. J. High Energy Phys. 2013, 7, 169. [Google Scholar] [CrossRef] [Green Version]
  23. de Florian, D.; Mazzitelli, J. Higgs pair production at next-to-next-to-leading logarithmic accuracy at the LHC. J. High Energy Phys. 2015, 9, 053. [Google Scholar] [CrossRef] [Green Version]
  24. Grazzini, M.; Heinrich, G.; Jones, S.; Kallweit, S.; Kerner, M.; Lindert, J.M.; Mazzitelli, J. Higgs boson pair production at NNLO with top quark mass effects. J. High Energy Phys. 2018, 5, 059. [Google Scholar] [CrossRef] [Green Version]
  25. Baglio, J.; Campanario, F.; Glaus, S.; Mühlleitner, M.; Ronca, J.; Spira, M. ggHH: Combined uncertainties. Phys. Rev. D 2021, 103, 056002. [Google Scholar] [CrossRef]
  26. Baglio, J.; Djouadi, A.; Gröber, R.; Mühlleitner, M.M.; Quevillon, J.; Spira, M. The measurement of the Higgs self-coupling at the LHC: Theoretical status. J. High Energy Phys. 2013, 4, 151. [Google Scholar] [CrossRef] [Green Version]
  27. Frederix, R.; Frixione, S.; Hirschi, V.; Maltoni, F.; Mattelaer, O.; Torrielli, P.; Vryonidou, E.; Zaro, M. Higgs pair production at the LHC with NLO and parton-shower effects. Phys. Lett. B 2014, 732, 142–149. [Google Scholar] [CrossRef]
  28. Ling, L.S.; Zhang, R.Y.; Ma, W.G.; Guo, L.; Li, W.H.; Li, X.Z. NNLO QCD corrections to Higgs pair production via vector boson fusion at hadron colliders. Phys. Rev. D 2014, 89, 073001. [Google Scholar] [CrossRef] [Green Version]
  29. Dreyer, F.A.; Karlberg, A. Fully differential Vector-Boson Fusion Higgs Pair Production at Next-to-Next-to-Leading Order. Phys. Rev. D 2019, 99, 074028. [Google Scholar] [CrossRef] [Green Version]
  30. Dreyer, F.A.; Karlberg, A. Vector-Boson Fusion Higgs Pair Production at N3LO. Phys. Rev. D 2018, 98, 114016. [Google Scholar] [CrossRef] [Green Version]
  31. Carvalho, A.; Dall’Osso, M.; Dorigo, T.; Goertz, F.; Gottardo, C.A.; Tosi, M. Higgs Pair Production: Choosing Benchmarks With Cluster Analysis. J. High Energy Phys. 2016, 4, 126. [Google Scholar] [CrossRef] [Green Version]
  32. Maltoni, F.; Vryonidou, E.; Zhang, C. Higgs production in association with a top-antitop pair in the Standard Model Effective Field Theory at NLO in QCD. J. High Energy Phys. 2016, 10, 123. [Google Scholar] [CrossRef]
  33. Heinrich, G.; Jones, S.P.; Kerner, M.; Scyboz, L. A non-linear EFT description of ggHH at NLO interfaced to POWHEG. J. High Energy Phys. 2020, 10, 021. [Google Scholar] [CrossRef]
  34. Carvalho, A.; Dall’Osso, M.; De Castro Manzano, P.; Dorigo, T.; Goertz, F.; Gouzevich, M.; Tosi, M. Analytical Parametrization and Shape Classification of Anomalous HH Production in the EFT Approach. Available online: https://cds.cern.ch/record/2199287 (accessed on 1 April 2022).
  35. Carvalho, A.; Goertz, F.; Mimasu, K.; Gouzevitch, M.; Aggarwal, A. On the reinterpretation of non-resonant searches for Higgs boson pairs. J. High Energy Phys. 2021, 2, 049. [Google Scholar] [CrossRef]
  36. Capozi, M.; Heinrich, G. Exploring anomalous couplings in Higgs boson pair production through shape analysis. J. High Energy Phys. 2020, 3, 091. [Google Scholar] [CrossRef] [Green Version]
  37. Bishara, F.; Contino, R.; Rojo, J. Higgs pair production in vector-boson fusion at the LHC and beyond. Eur. Phys. J. C 2017, 77, 481. [Google Scholar] [CrossRef] [Green Version]
  38. Giudice, G.F.; Grojean, C.; Pomarol, A.; Rattazzi, R. The Strongly-Interacting Light Higgs. J. High Energy Phys. 2007, 6, 045. [Google Scholar] [CrossRef] [Green Version]
  39. Contino, R.; Grojean, C.; Moretti, M.; Piccinini, F.; Rattazzi, R. Strong Double Higgs Production at the LHC. J. High Energy Phys. 2010, 5, 089. [Google Scholar] [CrossRef] [Green Version]
  40. Guerrero Ibarra, D.F. Search for Non-Resonant Higgs Boson Pair Production in the Four Bottom Quark Decay Channel with the CMS Experiment. Ph.D. Thesis, School of the University of Florida, Gainesville, FL, USA, 2021. Available online: https://cds.cern.ch/record/2791685 (accessed on 1 April 2022).
  41. Heinrich, G.; Jones, S.P.; Kerner, M.; Luisoni, G.; Vryonidou, E. NLO predictions for Higgs boson pair production with full top quark mass dependence matched to parton showers. J. High Energy Phys. 2017, 8, 088. [Google Scholar] [CrossRef]
  42. Jones, S.; Kuttimalai, S. Parton Shower and NLO-Matching uncertainties in Higgs Boson Pair Production. J. High Energy Phys. 2018, 2, 176. [Google Scholar] [CrossRef] [Green Version]
  43. ATLAS Collaboration. Combination of searches for Higgs boson pairs in pp collisions at s = 13 TeV with the ATLAS detector. Phys. Lett. B 2020, 800, 135103. [Google Scholar] [CrossRef]
  44. CMS Collaboration. Combination of searches for Higgs boson pair production in proton-proton collisions at s = 13 TeV. Phys. Rev. Lett. 2019, 122, 121803. [Google Scholar] [CrossRef] [Green Version]
  45. ATLAS Collaboration. Search for pair production of Higgs bosons in the bb¯bb¯ final state using proton-proton collisions at s =13 TeV with the ATLAS detector. J. High Energy Phys. 2019, 1, 030. [Google Scholar]
  46. CMS Collaboration. Search for nonresonant Higgs boson pair production in the bb¯bb¯ final state at s = 13 TeV. J. High Energy Phys. 2019, 4, 112. [Google Scholar]
  47. ATLAS Collaboration. Search for resonant and non-resonant Higgs boson pair production in the b b ¯ τ+τ decay channel in pp collisions at s = 13 TeV with the ATLAS detector. Phys. Rev. Lett. 2018, 121, 191801. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. CMS Collaboration. Search for Higgs boson pair production in events with two bottom quarks and two tau leptons in proton–proton collisions at s = 13 TeV. Phys. Lett. B 2018, 778, 101–127. [Google Scholar] [CrossRef]
  49. ATLAS Collaboration. Search for Higgs boson pair production in the γγb b ¯ final state with 13 TeV pp collision data collected by the ATLAS experiment. J. High Energy Phys. 2018, 11, 040. [Google Scholar]
  50. CMS Collaboration. Search for Higgs boson pair production in the γγb b ¯ final state in pp collisions at s = 13 TeV. Phys. Lett. B 2019, 788, 7–36. [Google Scholar] [CrossRef]
  51. ATLAS Collaboration. Search for Higgs boson pair production in the b b ¯ WW* decay mode at s =13 TeV with the ATLAS detector. J. High Energy Phys. 2019, 4, 092. [Google Scholar]
  52. CMS Collaboration. Search for resonant and nonresonant Higgs boson pair production in the b b ¯ νν final state in proton-proton collisions at s =13 TeV. J. High Energy Phys. 2018, 1, 054. [Google Scholar]
  53. ATLAS Collaboration. Search for Higgs boson pair production in the γγWW* channel using pp collision data recorded at s = 13 TeV with the ATLAS detector. Eur. Phys. J. C 2018, 78, 1007. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. ATLAS Collaboration. Search for Higgs boson pair production in the WW(*)WW(*) decay channel using ATLAS data recorded at s = 13 TeV. J. High Energy Phys. 2019, 5, 124. [Google Scholar]
  55. CMS Collaboration. Search for Higgs boson pair production in the four b quark final state in proton-proton collisions at s = 13 TeV. In CMS Physics Analysis Summary; CERN: Geneva, Switzerland, 2022; Available online: https://cds.cern.ch/record/2801882 (accessed on 1 April 2022).
  56. ATLAS Collaboration. Search for the HHb b ¯ b b ¯ process via vector-boson fusion production using proton-proton collisions at s = 13 TeV with the ATLAS detector. J. High Energy Phys. 2020, 7, 108. [Google Scholar]
  57. CMS Collaboration. Search for Higgs boson pair production via vector boson fusion with highly Lorentz-boosted Higgs bosons in the four b quark final state at s = 13 TeV. In CMS Physics Analysis Summary; CERN: Geneva, Switzerland, 2021; Available online: https://cds.cern.ch/record/2776802 (accessed on 1 April 2022).
  58. CMS Collaboration. Nonresonant pair production of highly energetic Higgs bosons decaying to bottom quarks. In CMS Physics Analysis Summary; CERN: Geneva, Switzerland, 2022; Available online: https://cds.cern.ch/record/2803672 (accessed on 1 April 2022).
  59. Qu, H.; Gouskos, L. ParticleNet: Jet Tagging via Particle Clouds. Phys. Rev. D 2020, 101, 056019. [Google Scholar] [CrossRef] [Green Version]
  60. CMS Collaboration. Identification of highly Lorentz-Boosted Heavy Particles Using Graph Neural Networks and New Mass Decorrelation Techniques. In CMS Detector Perfomance Summary; CERN: Geneva, Switzerland, 2020; Available online: https://cds.cern.ch/record/2707946 (accessed on 1 April 2022).
  61. ATLAS Collaboration. Search for Resonant and Non-Resonant Higgs Boson Pair Production in the b b ¯ τ+τ Decay Channel Using 13 TeV pp Collision Data from the ATLAS Detector. In ATLAS Conference Note; CERN: Geneva, Switzerland, 2021; Available online: https://cds.cern.ch/record/2777236 (accessed on 1 April 2022).
  62. CMS Collaboration. Search for nonresonant Higgs boson pair production in final states with two bottom quarks and two tau leptons in proton-proton collisions at s = 13 TeV. In CMS Physics Analysis Summary; CERN: Geneva, Switzerland, 2022; Available online: https://cds.cern.ch/record/2803419 (accessed on 1 April 2022).
  63. ATLAS Collaboration. Search for Higgs boson pair production in the two bottom quarks plus two photons final state in pp collisions at s = 13 TeV with the ATLAS detector. In ATLAS Conference Note; CERN: Geneva, Switzerland, 2021; Available online: http://cds.cern.ch/record/2798691 (accessed on 1 April 2022).
  64. CMS Collaboration. Search for nonresonant Higgs boson pair production in final states with two bottom quarks and two photons in proton-proton collisions at s = 13 TeV. J. High Energy Phys. 2021, 3, 257. [Google Scholar]
  65. CMS Collaboration. Search for Higgs boson pairs decaying to WWWW, WWττ, and ττττ in proton-proton collisions at s = 13 TeV . In CMS Physics Analysis Summary; CERN: Geneva, Switzerland, 2022; Available online: https://cds.cern.ch/record/2799151 (accessed on 1 April 2022).
  66. CMS Collaboration. Search for nonresonant Higgs boson pair production in the 4 leptons plus 2 b jets final state in proton-proton collisions at s = 13 TeV. In CMS Physics Analysis Summary; CERN: Geneva, Switzerland, 2020; Available online: https://cds.cern.ch/record/2725691 (accessed on 1 April 2022).
  67. ATLAS Collaboration. Search for non-resonant Higgs boson pair production in the bbνν final state with the ATLAS detector in pp collisions at s = 13 TeV. Phys. Lett. B 2020, 801, 135145. [Google Scholar] [CrossRef]
  68. ATLAS Collaboration. HEFT interpretations of Higgs boson pair searches in b b ¯ γ γ and b b ¯ τ τ final states and of their combination in ATLAS. In ATLAS PUB Note; CERN: Geneva, Switzerland, 2022; Available online: https://cds.cern.ch/record/2806411 (accessed on 1 April 2022).
  69. ATLAS Collaboration. Combination of searches for non-resonant and resonant Higgs boson pair production in the b b ¯ γ γ , b b ¯ τ+τ and b b ¯ b b ¯ decay channels using p p collisions at s = 13 TeV with the ATLAS detector. In ATLAS Conference Note; CERN: Geneva, Switzerland, 2021; Available online: https://cds.cern.ch/record/2786865 (accessed on 1 April 2022).
  70. ATLAS Collaboration. New Small Wheel Technical Design Report. 2013. Available online: https://cds.cern.ch/record/1552862 (accessed on 1 April 2022).
  71. ATLAS Collaboration. ATLAS Liquid Argon Calorimeter Phase-I Upgrade: Technical Design Report. 2013. Available online: https://cds.cern.ch/record/1602230 (accessed on 1 April 2022).
  72. ATLAS Collaboration. Technical Design Report for the Phase-I Upgrade of the ATLAS TDAQ System. 2013. Available online: https://cds.cern.ch/record/1602235 (accessed on 1 April 2022).
  73. CMS Tracker Group. The CMS Phase-1 Pixel Detector Upgrade. J. Instrum. 2021, 16, P02027. [Google Scholar] [CrossRef]
  74. CMS Collaboration. CMS Technical Design Report for the Phase 1 Upgrade of the Hadron Calorimeter. 2012. Available online: https://cds.cern.ch/record/1481837 (accessed on 1 April 2022).
  75. CMS Collaboration. CMS Technical Design Report for the Muon Endcap GEM Upgrade. 2015. Available online: https://cds.cern.ch/record/2021453 (accessed on 1 April 2022).
  76. CMS Collaboration. CMS Technical Design Report for the Level-1 Trigger Upgrade. 2013. Available online: https://cds.cern.ch/record/1556311 (accessed on 1 April 2022).
  77. Apollinari, G.; Béjar Alonso, I.; Brüning, O.; Fessia, P.; Lamont, M.; Rossi, L.; Tavian, L. High-Luminosity Large Hadron Collider (HL-LHC): Technical Design Report V. 0.1; CERN Yellow Report; CERN: Geneva, Switzerland, 2017. [Google Scholar]
  78. ATLAS Collaboration. ATLAS Phase-II Upgrade Scoping Document; CERN: Geneva, Switzerland, 2015; Available online: http://cds.cern.ch/record/2055248 (accessed on 1 April 2022).
  79. CMS Collaboration. CMS Phase II Upgrade Scope Document; CERN: Geneva, Switzerland, 2015; Available online: https://cds.cern.ch/record/2055167 (accessed on 1 April 2022).
  80. Cepeda, M.; Gori, S.; Ilten, P.; Kado, M.; Riva, F.; Abdul Khalek, R.; Aboubrahim, A.; Alimena, J.; Alioli, S.; Lves, A.; et al. Report from Working Group 2: Higgs Physics at the HL-LHC and HE-LHC. In CERN Yellow Report; CERN: Geneva, Switzerland, 2019; Volume 7, pp. 221–584. [Google Scholar]
  81. ATLAS Collaboration. Projected sensitivity of Higgs boson pair production combining the b b ¯ γ γ and b b ¯ τ+τ final states with the ATLAS detector at the HL-LHC. In ATLAS PUB Note ATL-PHYS-PUB-2022-005; CERN: Geneva, Switzerland, 2022; Available online: http://cds.cern.ch/record/2802127 (accessed on 1 April 2022).
  82. Ovyn, S.; Rouby, X.; Lemaitre, V. DELPHES, a Framework for Fast Simulation of a Generic Collider Experiment. Available online: https://arXiv.org/abs/0903.2225 (accessed on 1 April 2022).
  83. CMS Collaboration. Prospects for non-resonant Higgs boson pair production measurement in bbγγ final states in proton-proton collisions at s = 14 TeV at the High-Luminosity LHC. In CMS Physics Analysis Summary; CERN: Geneva, Switzerland, 2022; Available online: https://cds.cern.ch/record/2803918 (accessed on 1 April 2022).
  84. CMS Collaboration. Prospects for HH measurements in the WW γ γ and τ τ γ γ final states in proton-proton collisions at s = 14 TeV at the High Luminosity-LHC. In CMS Physics Analysis Summary; CERN: Geneva, Switzerland, 2022; Available online: https://cds.cern.ch/record/2804003 (accessed on 1 April 2022).
  85. LHC Higgs Cross Section Working Group. Handbook of LHC Higgs Cross Sections: 4. Deciphering the Nature of the Higgs Sector. 2016. Available online: http://www.arXiv.org/abs/1610.07922 (accessed on 1 April 2022).
  86. CMS Collaboration. Search for the nonresonant t t ¯ HH production in the semileptonic decay of the top pair and the Higgs pair decay into b quarks at the HL-LHC. In CMS Physics Analysis Summary; CERN: Geneva, Switzerland, 2022; Available online: https://cds.cern.ch/record/2804085 (accessed on 1 April 2022).
  87. Di Micco, B.; Gouzevitch, M.; Mazzitelli, J.; Vernieri, C. Higgs boson potential at colliders: Status and perspectives. Rev. Phys. 2020, 5, 100045. [Google Scholar] [CrossRef]
  88. Di Vita, S.; Durieux, G.; Grojean, C.; Gu, J.; Liu, Z.; Panico, G.; Riembau, M.; Vantalon, T. A global view on the Higgs self-coupling at lepton colliders. J. High Energy Phys. 2018, 2, 178. [Google Scholar] [CrossRef] [Green Version]
  89. Han, T.; Liu, D.; Low, I.; Wang, X. Electroweak couplings of the Higgs boson at a multi-TeV muon collider. Phys. Rev. D 2021, 103, 013002. [Google Scholar] [CrossRef]
  90. Mangano, M.L.; Ortona, G.; Selvaggi, M. Measuring the Higgs self-coupling via Higgs-pair production at a 100 TeV p-p collider. Eur. Phys. J. C 2020, 80, 1030. [Google Scholar] [CrossRef]
Figure 1. Running of the value of λ as function of the energy scale (a) and regions of stability (green), meta-stability (yellow) and instability (red) as function of the H and top masses (b). The grey area denotes the values preferred by the experimental measurements [10].
Figure 1. Running of the value of λ as function of the energy scale (a) and regions of stability (green), meta-stability (yellow) and instability (red) as function of the H and top masses (b). The grey area denotes the values preferred by the experimental measurements [10].
Symmetry 14 01467 g001
Figure 2. Representative diagrams at the lowest order for the production of HH via gluon fusion (a) and vector boson fusion (b).
Figure 2. Representative diagrams at the lowest order for the production of HH via gluon fusion (a) and vector boson fusion (b).
Symmetry 14 01467 g002
Figure 3. Total HH production cross section as function of κ λ [27].
Figure 3. Total HH production cross section as function of κ λ [27].
Symmetry 14 01467 g003
Figure 4. Distributions of the box, interference, and triangle components that contribute to the SM signal (a). Since these have, respectively, no, linear, and quadratic dependence on κ λ , the m HH shapes largely differ for various coupling hypotheses (b).
Figure 4. Distributions of the box, interference, and triangle components that contribute to the SM signal (a). Since these have, respectively, no, linear, and quadratic dependence on κ λ , the m HH shapes largely differ for various coupling hypotheses (b).
Symmetry 14 01467 g004
Figure 5. Representative diagrams at the lowest order for ggF HH production in and EFT description. The diagram associated to the chromomagnetic operator, that appear only in the SMEFT formalism but not in the HEFT one, and is thus far not modeled in the experimental analyses, is not shown here.
Figure 5. Representative diagrams at the lowest order for ggF HH production in and EFT description. The diagram associated to the chromomagnetic operator, that appear only in the SMEFT formalism but not in the HEFT one, and is thus far not modeled in the experimental analyses, is not shown here.
Symmetry 14 01467 g005
Figure 6. Distribution of m HH for VBF HH events for the SM (red) and for values of κ 2 V corresponding to 0 (no VV H H interaction, orange) and 2 (blue) [40].
Figure 6. Distribution of m HH for VBF HH events for the SM (red) and for values of κ 2 V corresponding to 0 (no VV H H interaction, orange) and 2 (blue) [40].
Symmetry 14 01467 g006
Figure 7. Value of the 13 TeV VBF HH cross section as function of the κ 2 V coupling (a) and as a two-dimensional function of ( κ 2 V , κ λ ) (b) and ( κ 2 V , κ V ) (c). For all figures, all the other HH couplings not explicitly indicated on the axes are assumed to be fixed to their SM value.
Figure 7. Value of the 13 TeV VBF HH cross section as function of the κ 2 V coupling (a) and as a two-dimensional function of ( κ 2 V , κ λ ) (b) and ( κ 2 V , κ V ) (c). For all figures, all the other HH couplings not explicitly indicated on the axes are assumed to be fixed to their SM value.
Symmetry 14 01467 g007
Figure 8. Summary of the branching fractions for the main HH final states.
Figure 8. Summary of the branching fractions for the main HH final states.
Symmetry 14 01467 g008
Figure 9. Upper limits at the 95% CL on the signal strength of the SM HH signal from the combination of the results based on the partial Run 2 data set collected in 2016 (approximately 36 fb 1 ) by the ATLAS (a) and CMS (b) Collaborations [43,44].
Figure 9. Upper limits at the 95% CL on the signal strength of the SM HH signal from the combination of the results based on the partial Run 2 data set collected in 2016 (approximately 36 fb 1 ) by the ATLAS (a) and CMS (b) Collaborations [43,44].
Symmetry 14 01467 g009
Figure 10. Upper limits at the 95% CL on the HH production cross section as function of κ λ from combination of the results based on the partial Run 2 data set collected in 2016 (approximately 36 fb 1 ) by the ATLAS (a) and CMS (b) Collaborations. The intersection with the theoretical prediction for the production cross section defines the range of κ λ values that are excluded by these searches [43,44].
Figure 10. Upper limits at the 95% CL on the HH production cross section as function of κ λ from combination of the results based on the partial Run 2 data set collected in 2016 (approximately 36 fb 1 ) by the ATLAS (a) and CMS (b) Collaborations. The intersection with the theoretical prediction for the production cross section defines the range of κ λ values that are excluded by these searches [43,44].
Symmetry 14 01467 g010
Figure 11. Four-body-mass distribution, m b b b b , in the signal region of the ATLAS VBF b b ¯ b b ¯ search. Two BSM HH signal models are overlaid on the stacked backgrounds and data points [56].
Figure 11. Four-body-mass distribution, m b b b b , in the signal region of the ATLAS VBF b b ¯ b b ¯ search. Two BSM HH signal models are overlaid on the stacked backgrounds and data points [56].
Symmetry 14 01467 g011
Figure 12. Distributions of the events observed in the CMS resolved b b ¯ b b ¯ analysis. Distribution of the invariant mass of the HH system in the VBF selection (a), and of the BDT output score for one ggF category (b) [55].
Figure 12. Distributions of the events observed in the CMS resolved b b ¯ b b ¯ analysis. Distribution of the invariant mass of the HH system in the VBF selection (a), and of the BDT output score for one ggF category (b) [55].
Symmetry 14 01467 g012
Figure 13. Distributions of the events observed in the CMS boosted b b ¯ b b ¯ analysis. Expected number of events for the VBF categories (a), and distribution of the invariant mass of one H candidate in the ggF categories (b) [58].
Figure 13. Distributions of the events observed in the CMS boosted b b ¯ b b ¯ analysis. Expected number of events for the VBF categories (a), and distribution of the invariant mass of one H candidate in the ggF categories (b) [58].
Symmetry 14 01467 g013
Figure 14. Final fit distributions in the b b ¯ τ τ + searches. The ATLAS BDT output score for the (most sensitive) τ h a d τ h a d category is shown on the left. [61] The CMS aggregated DNN output score distribution, sorted by the log of the signal-to-background ratio, is shown on the right [62].
Figure 14. Final fit distributions in the b b ¯ τ τ + searches. The ATLAS BDT output score for the (most sensitive) τ h a d τ h a d category is shown on the left. [61] The CMS aggregated DNN output score distribution, sorted by the log of the signal-to-background ratio, is shown on the right [62].
Symmetry 14 01467 g014
Figure 15. Final fit distributions in the b b ¯ γ γ searches. The ATLAS m γ γ spectrum in the (most sensitive) high-mass, BDT tight category is shown in the top panel [63]. The CMS m γ γ and m j j spectra for all analysis categories are shown in the bottom panel [64].
Figure 15. Final fit distributions in the b b ¯ γ γ searches. The ATLAS m γ γ spectrum in the (most sensitive) high-mass, BDT tight category is shown in the top panel [63]. The CMS m γ γ and m j j spectra for all analysis categories are shown in the bottom panel [64].
Symmetry 14 01467 g015
Figure 16. Seven categories in CMS HH multilepton search. The observed and expected 95% CL upper limits on σ / σ S M for each category and their combination are shown (a) [65]. Inclusive BDT distribution in the CMS b b ¯ Z Z * ( 4 ) search. Data and signal and background processes are shown for a combination of all nine categories (b) [66]. The discriminant d H H for the ATLAS bb ν ν search. Data and background processes are shown, with the signal yield increased by a factor of 20 overlaid (c) [67].
Figure 16. Seven categories in CMS HH multilepton search. The observed and expected 95% CL upper limits on σ / σ S M for each category and their combination are shown (a) [65]. Inclusive BDT distribution in the CMS b b ¯ Z Z * ( 4 ) search. Data and signal and background processes are shown for a combination of all nine categories (b) [66]. The discriminant d H H for the ATLAS bb ν ν search. Data and background processes are shown, with the signal yield increased by a factor of 20 overlaid (c) [67].
Symmetry 14 01467 g016
Figure 17. The 95% CL limits on σ / σ S M , for each analysis presented in this section. ATLAS analyses are shown in blue/red, and CMS analyses are shown in green/yellow.
Figure 17. The 95% CL limits on σ / σ S M , for each analysis presented in this section. ATLAS analyses are shown in blue/red, and CMS analyses are shown in green/yellow.
Symmetry 14 01467 g017
Figure 18. 95% CL upper limit on the VBF HH cross section set by the boosted CMS b b ¯ b b ¯ analysis as function of the κ 2 V coupling (a) and negative log-likelihood contours as function of κ 2 V and κ V (b) [58].
Figure 18. 95% CL upper limit on the VBF HH cross section set by the boosted CMS b b ¯ b b ¯ analysis as function of the κ 2 V coupling (a) and negative log-likelihood contours as function of κ 2 V and κ V (b) [58].
Symmetry 14 01467 g018
Figure 19. The 95% CL upper limit on the HH EFT benchmark defined at LO ((a) CMS b b ¯ γ γ analysis [64]) and at NLO ((b) ATLAS b b ¯ γ γ + b b ¯ τ τ + combination [68]).
Figure 19. The 95% CL upper limit on the HH EFT benchmark defined at LO ((a) CMS b b ¯ γ γ analysis [64]) and at NLO ((b) ATLAS b b ¯ γ γ + b b ¯ τ τ + combination [68]).
Symmetry 14 01467 g019
Figure 20. 95% CL upper limits on the signal strength of the SM HH signal (a) and on the self-coupling modifier, κ λ (b), set by the ATLAS Collaboration, using 139 fb 1 of data collected during Run 2 [69].
Figure 20. 95% CL upper limits on the signal strength of the SM HH signal (a) and on the self-coupling modifier, κ λ (b), set by the ATLAS Collaboration, using 139 fb 1 of data collected during Run 2 [69].
Symmetry 14 01467 g020
Figure 21. Negative log-likelihood scan as a function of the self-coupling modifier κ λ , from a conditional signal+background fit to the background and the SM HH signal. The colored dashed lines correspond to the combined ATLAS and CMS results by channel, and the black line shows the ATLAS + CMS combination [80].
Figure 21. Negative log-likelihood scan as a function of the self-coupling modifier κ λ , from a conditional signal+background fit to the background and the SM HH signal. The colored dashed lines correspond to the combined ATLAS and CMS results by channel, and the black line shows the ATLAS + CMS combination [80].
Symmetry 14 01467 g021
Table 1. Summary of the cross section predicted by the SM for the various HH production modes at the LHC at s = 13 TeV . For ggF, FTapprox refers to the approximate full top mass approximation. In all cases, a reference value of m H = 125 GeV is used.
Table 1. Summary of the cross section predicted by the SM for the various HH production modes at the LHC at s = 13 TeV . For ggF, FTapprox refers to the approximate full top mass approximation. In all cases, a reference value of m H = 125 GeV is used.
Prod. Mode σ [ fb ]Scale Unc. m t Unc.PDF + α S Unc.Precision
ggF31.05 5.0 % + 2.2 % 23 % + 6 % ± 3 % NNLO FTapprox [18,19,20,21,22,23,24,25]
VBF1.726 0.04 % + 0.03 % - ± 2.1 % N 3 LO [26,27,28,29,30]
W + HH 0.329 0.41 % + 0.32 % - ± 2.2 % NNLO [26]
W HH 0.173 1.3 % + 1.2 % - ± 2.8 % NNLO [26]
Z HH 0.363 2.7 % + 3.4 % - ± 1.8 % NNLO [26]
ttHH0.775 4.3 % + 1.5 % - ± 3.2 % NLO [27]
tj HH 0.0289 3.6 % + 5.5 % - ± 4.7 % NLO [27]
Table 2. Definition of the LO [31,34,35] and NLO [36] EFT shape benchmarks. For the NLO benchmarks, please note that we used the conversion c g = 3 2 c g g h and c 2 g = 3 c g g h h , where c g g h and c g g h h are the coefficient of the associated HEFT operators. The conversion table can be found in Ref. [33].
Table 2. Definition of the LO [31,34,35] and NLO [36] EFT shape benchmarks. For the NLO benchmarks, please note that we used the conversion c g = 3 2 c g g h and c 2 g = 3 c g g h h , where c g g h and c g g h h are the coefficient of the associated HEFT operators. The conversion table can be found in Ref. [33].
κ λ κ t c 2 c g c 2 g
LO benchmarks
17.51.0−1.00.00.0
21.01.00.5−0.80.6
31.01.0−1.50.0−0.8
4−3.51.5−3.00.00.0
51.01.00.00.8−1.0
62.41.00.00.2−0.2
75.01.00.00.2−0.2
815.01.00.0−1.01.0
91.01.01.0−0.60.6
1010.01.5−1.00.00.0
112.41.00.01.0−1.0
1215.01.01.00.00.0
NLO benchmarks
13.940.94−1/33/4−1
26.840.611/301
32.211.05−1/30−3/2
42.790.611/31/4−1/2
53.951.17−1/3−3/43/2
65.680.831/31/2−1
7−0.10.941−1/41/2
Table 3. 95% CL upper limits on σ / σ S M for HH searches in the partial Run 2 data set.
Table 3. 95% CL upper limits on σ / σ S M for HH searches in the partial Run 2 data set.
σ / σ SM Obs. (Exp.) LimitRef.
b b ¯ b b ¯ ATLAS12.9 (20.7)[45]
CMS74.6 (36.9)[46]
b b ¯ τ τ + ATLAS12.6 (14.6)[47]
CMS31.4 (25.1)[48]
b b ¯ γ γ ATLAS20.4 (26.3)[49]
CMS23.6 (18.8)[50]
b b ¯ VV * ATLAS300 (300)[51]
CMS79 (89)[52]
W ± W γ γ ATLAS230 (160)[53]
W ± W W ± W ATLAS160 (120)[54]
Table 4. The 95% CL upper limits on σ / σ S M in the Run 2 HH searches presented in this section, improvement compared to previous searches (where applicable), and 95% CL exclusion range for the self-coupling modifer, κ λ .
Table 4. The 95% CL upper limits on σ / σ S M in the Run 2 HH searches presented in this section, improvement compared to previous searches (where applicable), and 95% CL exclusion range for the self-coupling modifer, κ λ .
σ / σ SM
Obs. (Exp.)
Improvement
w.r.t 36 fb 1
κ λ Exclusion
Obs. (Exp.)
Ref.
b b ¯ b b ¯
(boosted)
CMS3.9 (7.8)4.7×[−2.3,9.4] ([−5.0,12.0])[55]
CMS9.9 (5.1)30×[−9.9,16.9] ([−5.1,12.1])[58]
b b ¯ τ τ + ATLAS4.7 (3.9)3.8×[−2.4,9.2] ([−2.0,9.0])[61]
CMS3.3 (5.2)4.8×[−1.8,8.8] ([−3.0,9.9])[62]
b b ¯ γ γ ATLAS4.2 (5.7)4.6×[−1.5,6.7] ([−2.4,7.7])[63]
CMS7.7 (5.2)3.6×[−3.3,8.5] ([−2.5,8.2])[64]
multi-CMS21.8 (19.6)-[−7.0,11.7] ([−7.0,11.2])[65]
bbZZ(4)CMS30 (37)-[−9.0,14.0] ([−10.5,15.5])[66]
bb ν ν ATLAS40 (29)--[67]
Table 5. The 95% CL upper limits on σ V B F / σ S M V B F in the Run 2 VBF HH searches presented in this section and 95% CL exclusion range for the H H V V coupling modifier, κ 2 V .
Table 5. The 95% CL upper limits on σ V B F / σ S M V B F in the Run 2 VBF HH searches presented in this section and 95% CL exclusion range for the H H V V coupling modifier, κ 2 V .
σ VBF / σ SM VBF
Obs. (Exp.)
κ 2 V Observed
Exclusion
κ 2 V Expected
Exclusion
Ref.
b b ¯ b b ¯ (boosted)CMS-0.62 < κ 2 V < 1.410.66 < κ 2 V < 1.37[58]
b b ¯ b b ¯ ATLAS840 (550)−0.43 < κ 2 V < 2.56−0.55 < κ 2 V < 2.72[56]
b b ¯ b b ¯ CMS226 (412)−0.1 < κ 2 V < 2.2−0.4 < κ 2 V < 2.5[55]
b b ¯ τ τ + CMS124 (154)−0.4 < κ 2 V < 2.6−0.6 < κ 2 V < 2.8[62]
b b ¯ γ γ CMS225 (208)−1.3 < κ 2 V < 3.5−0.9 < κ 2 V < 3.0[64]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Brost, E.; Cadamuro, L. Searching for Pairs of Higgs Bosons in the LHC Run 2 Dataset. Symmetry 2022, 14, 1467. https://doi.org/10.3390/sym14071467

AMA Style

Brost E, Cadamuro L. Searching for Pairs of Higgs Bosons in the LHC Run 2 Dataset. Symmetry. 2022; 14(7):1467. https://doi.org/10.3390/sym14071467

Chicago/Turabian Style

Brost, Elizabeth, and Luca Cadamuro. 2022. "Searching for Pairs of Higgs Bosons in the LHC Run 2 Dataset" Symmetry 14, no. 7: 1467. https://doi.org/10.3390/sym14071467

APA Style

Brost, E., & Cadamuro, L. (2022). Searching for Pairs of Higgs Bosons in the LHC Run 2 Dataset. Symmetry, 14(7), 1467. https://doi.org/10.3390/sym14071467

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop