Search for invisible decays of Higgs bosons in the vector boson fusion and associated ZH production modes

A search for invisible decays of Higgs bosons is performed using the vector boson fusion and associated ZH production modes. In the ZH mode, the Z boson is required to decay to a pair of charged leptons or a b b-bar quark pair. The searches use the 8 TeV pp collision dataset collected by the CMS detector at the LHC, corresponding to an integrated luminosity of up to 19.7 inverse femtobarns. Certain channels include data from 7 TeV collisions corresponding to an integrated luminosity of 4.9 inverse femtobarns. The searches are sensitive to non-standard-model invisible decays of the recently observed Higgs boson, as well as additional Higgs bosons with similar production modes and large invisible branching fractions. In all channels, the observed data are consistent with the expected standard model backgrounds. Limits are set on the production cross section times invisible branching fraction, as a function of the Higgs boson mass, for the vector boson fusion and ZH production modes. By combining all channels, and assuming standard model Higgs boson cross sections and acceptances, the observed (expected) upper limit on the invisible branching fraction at m[H] = 125 GeV is found to be 0.58 (0.44) at 95% confidence level. We interpret this limit in terms of a Higgs-portal model of dark matter interactions.


Introduction
The discovery of a Higgs boson [1][2][3], together with the absence of any experimental hint of physics beyond the standard model (SM) at the Large Hadron Collider (LHC), have had a major impact on proposed theoretical models for new physics. All measurements of the recently observed 125 GeV boson to date indicate compatibility with the SM Higgs boson, but the associated uncertainties are large, and the possibility for non-SM properties remains. Moreover, although additional SM-like Higgs bosons have been excluded over a wide range of masses, additional Higgs bosons with exotic decay modes remain a possibility. * e-mail: cms-publication-committee-chair@cern.ch Invisible Higgs boson decays are possible in a wide range of models, for example through decays to neutralinos in supersymmetric models [4], or graviscalars in models with extra dimensions [5,6]. In general, interactions of the Higgs boson with the unknown dark matter (DM) sector may introduce invisible decay modes, and bounds on these decays can constrain DM models. In so-called "Higgs-portal" models of DM interactions [7][8][9], the Higgs boson takes the role of mediator between the SM particles and the DM particle. Recent theories proposing that the Higgs boson played a central role in the evolution of the early universe [10] provide further motivation to understand the relationship between the Higgs boson and DM.
Indirect constraints on non-SM decay modes of the recently observed Higgs boson have been inferred from the visible SM decay modes by including an additional non-SM partial width term in the combined fit to the data [3]. The resulting upper limit on the non-SM branching fraction is 0.89, at 95 % confidence level (CL). Direct searches for invisible Higgs boson decays, H(inv), are possible by requiring that the Higgs boson recoils against a visible system. Such searches were performed at LEP [11][12][13], using the ZH associated production mode. They excluded at 95 % CL an invisible Higgs boson of mass smaller than 105 GeV and produced with a cross section higher than 0.2 times the standard model ZH cross section. Phenomenological studies of hadron collider searches for H(inv) have considered all production mechanisms [14][15][16][17][18][19][20]. Recently, the ATLAS Collaboration reported a search for invisible decays of a Higgs boson produced in association with a Z boson that decays to leptons [21], placing an upper limit on the invisible Higgs boson branching fraction of 0.75 at 95 % CL for m H = 125.5 GeV. The ATLAS Collaboration also searched for an invisibly decaying Higgs boson in association with either a W or Z boson decaying hadronically [22].
Here we report searches for H(inv) in the ZH mode, where the Z boson decays to leptons or a bb quark pair, and the first search for H(inv) in the vector boson fusion (VBF) production mode, where the Higgs boson is produced in associa- Fig. 1 The Feynman diagrams for Higgs production in the VBF (left), Z( )H (center) and Z(bb)H (right) channels. The Higgs boson is assumed to decay invisibly tion with two quarks, as shown in Fig. 1 (left). Although the VBF signal benefits from a relatively large SM cross section, the final state of two jets plus missing transverse energy (E miss T ) suffers from large backgrounds. However, the backgrounds can be controlled by utilizing the distinct topology of the VBF process, in which the two jets are produced in a forward/backward configuration, with large invariant mass, and are well separated in rapidity. In addition, hadronic activity in the rapidity gap between the two scattered quarks is reduced, due to the absence of color flow in the VBF process. The ZH signal, shown in Fig. 1 (center) and (right), provides a complementary search to the VBF analysis. Despite a lower SM production cross section, the final state of a Z boson with large E miss T provides a clear topology with much lower backgrounds. We maximize the sensitivity of the search by including decays of the Z boson to leptons and bb quark pairs, which we refer to as Z( )H(inv), and Z(bb)H(inv), respectively, where represents either an electron or a muon. The Higgs boson production modes we consider here rely only on the Higgs boson coupling to the electroweak vector bosons. New physics that introduces invisible decays of the Higgs boson may also modify these couplings.
In the following sections of this article, we present a brief overview of the Compact Muon Solenoid (CMS) experimental apparatus, physics object reconstruction and datasets in Sects. 2 to 4, followed by a description of the event selection and background estimation for each of the three search channels in Sects. 5 to 7. We then present the results of the searches, and their combination, as upper limits on the production cross section times invisible branching fraction in Sect. 8. In Sect. 9 we interpret these cross section upper limits in terms of a Higgs-portal model of dark matter interactions, and we summarize our conclusions in Sect. 10.

The CMS apparatus
The central feature of the CMS apparatus is a superconducting solenoid of 6 m internal diameter, providing a magnetic field of 3.8 T. Within the volume of the superconducting solenoid are a silicon pixel and strip tracker, a lead tungstate crystal electromagnetic calorimeter (ECAL), and a brassscintillator hadron calorimeter, each composed by the barrel and endcap detectors. Muons are measured with detection planes made using three technologies: drift tubes, cathode strip chambers, and resistive-plate chambers, embedded in the steel flux-return yoke outside the solenoid. Extensive forward calorimetry complements the coverage provided by the barrel and endcap detectors. Data are selected online using a two-level trigger system. The first level, consisting of custom made hardware processors, selects events in less than 1 µs, while the high-level trigger processor farm further decreases the event rate from around 100 kHz to a few hundred Hz before data storage. The CMS experiment uses a righthanded coordinate system, with the origin at the nominal interaction point, the x axis pointing to the center of the LHC, the y axis pointing up (perpendicular to the LHC plane), and the z axis along the counterclockwise-beam direction. The polar angle θ is measured from the positive z axis and the azimuthal angle φ is measured in the x-y plane. The pseudorapidity, η, is defined as − ln[tan(θ/2)]. A more detailed description of the CMS apparatus can be found in Ref. [23].

Data samples and Monte Carlo simulation
The analyses presented here all use the 8 TeV data sample collected by the CMS Collaboration during 2012, corresponding to an integrated luminosity of 19.5 fb −1 in the VBF channel, 19.7 fb −1 in the Z( )H(inv) channel, and 18.9 fb −1 in the Z(bb)H(inv) channel. The Z( )H(inv) channel also uses the 7 TeV dataset collected during 2011, corresponding to 4.9 fb −1 . The uncertainty assigned to the luminosity measurement is 2.6 % (2.2 %) at √ s = 8 (7)TeV [24]. Backgrounds arising from sources other than pp collisions are suppressed using a set of filters that remove events due to anomalous calorimeter signals, beam halo identified in the muon endcaps, inoperable calorimeter cells, and tracking failure. We further require a well reconstructed vertex within the interaction region; |z| < 24 cm, r < 2 cm, where r = x 2 + y 2 .
The VBF signal is simulated using the powheg 2.0 event generator [25][26][27][28][29][30][31], while the Z( )H(inv) and Z(bb)H(inv) signals are simulated with pythia 6.4.26 [32]. The background processes are simulated using MadGraph 5.1.1 [33], with the exception of some minor backgrounds-specifically, the VH(bb) background to the Z(bb)H(inv) analysis is simulated with powheg 2.0, the diboson backgrounds in the VBF analysis are simulated with pythia 6.4.26, and the single-top-quark backgrounds in the VBF and Z( )H(inv) analyses use powheg 1.0. The QCD multijet background is simulated with pythia 6.4.26. All samples use the leadingorder CTEQ6L1 parton distribution functions (PDFs) [34], apart from the VH(bb) powheg samples, which use the nextto-leading-order (NLO) CTEQ6M PDFs [34]. Where yields are estimated directly from MC simulation, the PDF uncertainty is estimated using the PDF4LHC prescription [35,36]. For all Monte Carlo (MC) samples, the detector response is simulated using a detailed description of the CMS detector based on the Geant4 package [37]. Minimum bias events are superimposed on the generated events to simulate the effect of multiple pp interactions per bunch crossing (pileup). Simulated events are weighted such that the distribution of the number of pileup interactions reproduces that observed in data. The mean number of pileup interactions per bunch crossing was approximately 9 in 2011, and 21 in 2012. Additional weights are applied to simulated events to ensure trigger efficiency, lepton identification efficiency, and b-tagging efficiency match measurements from data.
The VBF and ZH production cross sections are taken from Refs. [38,39]. The ZH searches are performed in the boosted regime, where the Higgs boson has significant transverse momentum ( p T ), and thus, potential differences in the p T spectrum of the Z and Higgs bosons between data and MC generators could introduce systematic effects in the signal acceptance and efficiency estimates. Two sets of calculations are available that estimate the NLO electroweak corrections [40][41][42] and next-to-next-to-leading order (NNLO) QCD [43] corrections to vector boson plus Higgs boson production in the boosted regime. Both sets of corrections are applied to the signal MC samples. For VH production, the estimated uncertainty arising from the NLO electroweak corrections is 2 %, and from the NNLO QCD corrections is 5 %. In addition, we include NNLO electroweak corrections [44] to the ZZ and WZ background processes as a function of the p T of the Z boson.

Event reconstruction
The reconstructed interaction vertex with the largest value of i p T 2 i , where p Ti is the transverse momentum of the ith track associated with the vertex, is selected as the primary event vertex. This vertex is used as the reference vertex for all relevant objects in the event, which are reconstructed with a particle-flow algorithm [45,46]. The pileup interactions affect jet momentum reconstruction, missing transverse energy reconstruction, lepton isolation, and b-tagging efficiencies. To mitigate these effects, all charged-hadrons that do not originate from the primary interaction are identified by a particle-flow-based algorithm and removed from consideration in the event. In addition, following Ref. [47], the average neutral energy density from pileup interactions is evaluated on an event-by-event basis from particle-flow objects and used to compute a correction to the reconstructed jets in the event and to the summed energy in the isolation cones used for leptons.
Muons are reconstructed in the pseudorapidity range |η| < 2.4. Two muon reconstruction algorithms are used [48]: one in which tracks in the silicon tracker are matched to signals in the muon detectors, and another in which a global track fit is performed using hits in both the tracker and muon detectors. The muon candidates used in the analysis are required to be successfully reconstructed by both algorithms. The efficiency to reconstruct a muon of p T > 5 GeV is larger than 95 %, while the probability to misidentify a hadron as a muon is below 0.1 %. Further identification criteria are imposed on the muon candidates to reduce the fraction of tracks misidentified as muons. These include the number of measurements in the tracker and in the muon systems, the fit quality of the global muon track and its consistency with the primary vertex.
Electron reconstruction requires the matching of an energy cluster in the ECAL with a track in the silicon tracker [49]. Electron identification relies on a multivariate technique that combines observables sensitive to the amount of bremsstrahlung along the electron trajectory, the geometrical and momentum matching between the electron trajectory and associated clusters, as well as shower-shape observables. Additional requirements are imposed to remove electrons produced by photon conversions. In this analysis, electrons are considered in the pseudorapidity range |η| < 2.5, excluding the 1.44 < |η| < 1.57 transition region between the ECAL barrel and endcap, where electron reconstruction is suboptimal. Jets are reconstructed from particle-flow objects using the anti-k T clustering algorithm [50], with a distance parameter of 0.5, as implemented in the fastjet package [51,52]. Jets are found over the full calorimeter acceptance, |η| < 5. Jet energy corrections are applied as a function of the pseudorapidity and transverse momentum of the jet [53]. Jets resulting from pileup interactions are removed using a boosted decision tree (BDT), implemented in the TMVA package [54], with the following input variables: momentum and spatial distribution of the jet particles, charged-and neutral-particle multiplicities, and consistency of the charged hadrons within the jet with the primary vertex. The missing transverse momentum vector is calculated as the negative of the vectorial sum of the transverse momenta of all particle-flow objects identified in the event, and the magnitude of this vector is referred to as E miss T in the rest of this article. Jets that originate from the hadronization of b quarks are referred to as "b jets". The CSV b-tagging algorithm [55] is used to identify such jets. The algorithm combines the information about track impact parameters and secondary vertices within jets in a likelihood discriminant to provide separation between b jets and jets originating from light quarks, gluons, or charm quarks. The output of this CSV discriminant has values between zero and one; a jet with a CSV value above a certain threshold is referred to as being "b tagged". The efficiency to tag b jets and the rate of misidentification of non-b jets depend on the threshold chosen, and are typically parameterized as a function of the p T and η of the jets. These performance measurements are obtained directly from data in samples that can be enriched in b jets, such as tt and multijet events (where, for example, requiring the presence of a muon in the jets enhances the heavy-flavor content of the events). Several thresholds for the CSV output discriminant are used in this analysis. Depending on the threshold used, the efficiency to tag jets originating from b quarks is in the range 50-75 %, and the probability to incorrectly tag jets originating from c quarks, and light quarks or gluons as b jets are 5-25, and 0.15-3.0 %, respectively.

Search strategy
In the VBF mode, the Higgs boson is produced in association with two final-state quark jets separated by a large rapidity gap, and having high invariant mass. Loosely following the selection criteria discussed in Ref. [16], we select final states with two jets and large missing transverse energy and utilize the distinct topology of the VBF jets to discriminate the invisible Higgs boson signal from background.
The dominant backgrounds in this channel result from Z(νν)+jets, and W( ν)+jets, where the charged lepton is not identified. These backgrounds are estimated using control regions with a Z or W boson decaying to well identified charged leptons, in association with the same dijet topology used for the signal region. We then extrapolate from the control regions to the signal region using factors obtained from MC simulation. The background due to QCD multijet processes, where the E miss T arises from mismeasurement, is also estimated from data. Minor SM backgrounds, arising from tt, single-top, diboson, and Drell-Yan( )+jets processes are estimated from MC simulation.
We use the observed yield in the signal region, together with the estimated background, to perform a single-bin counting experiment.

Event selection
We use events collected with a trigger that requires E miss T > 65 GeV, in association with a pair of jets with p T j1 , p T j2 > 40 GeV, in a VBF-like topology. The jets are required to be in opposite forward/backward halves of the detector, well separated in pseudorapidity ( η jj = |η j1 − η j2 | > 3.5), and with high invariant mass (M jj > 800 GeV). For robustness against pileup, any pair of jets satisfying these criteria is accepted by the trigger. At the trigger level, the E miss T calculation does not include muons, allowing control samples of W(μν)+jets and Z(μμ)+jets events to be taken on the same trigger. The trigger efficiency is measured in events recorded on a singlemuon trigger, as a function of p T j2 (since the leading jet, j1, is effectively always above threshold for the regions considered), M jj , and E miss T , and the measured efficiency is applied to all MC samples.
The offline selection then proceeds as follows. We reject backgrounds from Z and W bosons by vetoing any event with an identified electron [49] or muon [56] with p T > 10 GeV. The VBF tag jet pair is then identified as the leading jet pair. This pair is required to pass tightened versions of the trigger selection, specifically p T j1 , p T j2 > 50 GeV, |η| < 4.7, η j1 , η j2 < 0, η jj > 4.2, and M jj > 1,100 GeV. The missing-energy requirement is E miss T > 130 GeV. Multijet backgrounds are reduced to a low level by requiring the azimuthal separation between the tag jets to be small, φ jj < 1.0 radians, since the background peaks at φ jj = π radians while the signal is roughly flat in φ jj . Finally, we apply a central-jet veto (CJV) to any event that has an additional jet with p T > 30 GeV and pseudorapidity between those of the two tag jets.
The lepton and central jet veto thresholds are set to low values at which reconstruction is known to be reliable, while the remaining thresholds are determined by optimizing the selection to give the best signal significance, calculated using a profile likelihood method that incorporates all systematic uncertainties, for a Higgs boson with m H = 125 GeV and 100 % invisible branching fraction. The thresholds on jet p T , M jj , and E miss T are constrained to be above the point where the trigger is 95 % efficient. This constraint effectively determines the jet p T and E miss T thresholds, since signal significance only worsens when these thresholds are raised above this point. Distributions of M jj , η jj , φ jj , and central jet p T in background and signal MC simulation are shown in Fig. 2, along with the thresholds applied after optimization of the selection.
After all selection requirements, an hypothetical signal equivalent to 125 GeV Higgs boson with B(H → inv) = 100 % and produced via the VBF process with SM couplings, is reconstructed with an efficiency of (6.8 ± 0.3) × 10 −3 , corresponding to a yield of 210 ± 29(syst) events. The requirements on the VBF tag jet p T and topology, M jj , and E miss T are all correlated and affect the signal efficiency by comparable amounts. A small signal yield from the gluonfusion process is also expected, where the VBF requirements may be satisfied by initial-state radiation. Based on powheg simulation, we estimate this to be 14 ± 10(syst) events.

Background estimation
The Z(νν)+jets background is estimated from data using observable Z(μμ) decays. We define a Z control region as for the signal region, with the following changes to the event selection: the lepton veto is replaced with a requirement of an oppositely charged pair of well reconstructed and isolated muons each with p T > 20 GeV, and invariant mass 60 < M μμ < 120 GeV, a veto is applied on any additional leptons with p T > 10 GeV, and the E miss T is recomputed after removing the muons from the Z boson decay. The number of Z(νν) events in the signal region is then predicted using: The ratio of cross sections, σ (Z → νν)/σ (Z/γ * → μμ) = 5.651 ± 0.023(syst) is calculated with mcfm [57] for m Z/γ * > 50 GeV, the mass range of the MC sample. The selection efficiencies in the signal region, ε s ZMC = (1.65 ± 0.27(syst)) × 10 −6 , and the control region, ε c ZMC = (1.11 ± 0.17(syst)) × 10 −6 , are estimated from DY( )+jets simulation, ignoring the muons when computing the efficiency in the signal region. The observed yield in the control region is N c μμobs = 12 events. The background in the control region-estimated from tt, diboson and single-top MC samples-is N c bkg = 0.23 ± 0.15(syst) events. The resulting estimate of the Z(νν) background in the signal region is 99 ± 29(stat) ± 25(syst) events. The source of systematic uncertainty in the background estimates will be described in Sect. 5.4. Figure 3 shows the E miss T and dijet invariant mass, M jj , distributions with a relaxed set of criteria for the Z control region, with M jj > 1,000 GeV and no requirements on η jj , φ jj , or CJV. In this figure, the simulated background is normalized to the data. It should be noted that our estimates of the dominant V+jets background are insensitive to the overall normalization of the simulation, which cancels in the ratio. φ jj , or CJV, and with the M jj requirement relaxed to 1,000 GeV. The simulated background from different processes is shown cumulatively, and normalized to the data, with its systematic uncertainty shown as a hatched region. The lower panels show the ratio of data to the simulated background, again with the systematic uncertainty shown as a hatched region The W(eν)+jets and W(μν)+jets backgrounds are estimated from single-lepton control samples. We define W(μν) and W(eν) control regions in a similar way to the Z boson background. In the W(μν) region, the lepton veto is replaced with a single μ requirement and a veto on any additional leptons, and the E miss T is recomputed after removing the muon from the W boson decay. The W(eν) region is defined similarly, with a single electron requirement and additional lepton veto, but here the E miss T is not recomputed, since the electron energy is already included in the E miss T at trigger level.

Events / 20 GeV
The number of W( ν) (where = e, μ) events in the signal region, N s is then estimated using: where N s WMC and N c WMC are the number of events in the signal and control regions in the W( ν)+jets MC simulation. The ratio N s WMC /N c WMC is equal to 0.347 ± 0.045(syst) for W(μν) and 1.08 ± 0.21(syst) for W(eν). In the W(μν) control region the observed yield is 223 events, with a background of 30.4 ± 7.0(syst) events. The observed yield in the W(eν) control region is 65 events, with a background of 7.1± 4.7(syst) events. The W(μν) background in the signal region is then estimated to be 66.8 ± 5.2(stat) ± 15.7(syst) events, and the W(eν) background to be 62.7±8.7(stat)±18.1(syst) events.
The background arising from W(τ ν)+jets, where the tau lepton decays hadronically (τ h ) is estimated using a slightly different method, since a tau lepton veto is not applied in the invisible Higgs boson signal selection. Hadronically decaying taus are reconstructed using the "hadron plus strips" algorithm [58]. This uses charged hadrons and neutral electromagnetic objects (photons) to reconstruct hadronic tau decay modes with one or three charged particles, in the range |η| < 2.3. A control region is defined, requiring one hadronic tau with p T > 20 GeV and |η| < 2.3, no additional leptons, and the remaining signal region selection. However, in the W(τ h ν) control region, the CJV is not applied in order to increase the yield. The number of W(τ h ν) events in the signal region, N s τ h , is then estimated from the control region in the same way as the W(μν) and W(eν) backgrounds. A yield of 32 events is observed in the control region, with the background estimated from the MC simulation to be 15.2 ± 3.6(syst) events, giving an estimate of the W(τ h ν) background in the signal region of 53±18(stat)±18(syst) events.
In order to cross check the backgrounds from V+jets processes (where V represents either a W or a Z boson), which dominate in the signal region, the W(μν) control region and MC simulation is used to compute yields in other control regions. For example, the yield in the Z(μμ) region is given by: Similar expressions are used to estimate yields in the W(eν) and W(τ h ν) control regions. In all cases, the predictions from data agree with the observed yield within the uncertainty.
The QCD multijet background in the signal region is estimated using the fractions of events passing the E miss T and CJV requirements. We define regions A, B, C, and D as follows, after the full remaining selection: We estimate the QCD multijet component in regions A, B, and C from data, after subtracting the electroweak backgrounds using estimations from simulation. The QCD multijet component in the signal region D can then be estimated using N D = N B N C /N A , where N i is the number of events in region i. This method is based on the assumption that the E miss T and the CJV are uncorrelated, which has been checked by comparing the E miss T distribution, below the 130 GeV threshold, in events passing and failing the CJV. The maximum difference in the E miss T distribution between these two samples is 40 %, which is assigned as a systematic uncertainty of the method. We predict the QCD background in the signal region to be 30.9 ± 4.8(stat) ± 23.0(syst) events. Furthermore, the method is tested on a high statistics sample with selections equivalent to those in the signal region, but dominated by QCD multijet events by changing the φ jj requirement to φ jj > 2.6 radians. In this sample, we observe 2,551 ± 57(stat) events in the pseudosignal region after subtraction of backgrounds, which are estimated from MC simulation. The QCD multijet component is predicted to be 2959 ± 58(stat), which is compatible with the observation within the systematic uncertainty. To give further confidence in this estimate, we perform a crosscheck using an ABCD method based on the E miss T and φ jj variables, which gives a prediction consistent with the main method.
The remaining SM backgrounds in the signal region-due to tt, single-top, VV and DY( )+jets-are estimated from MC simulation to be 20.0 +6.0 −8.2 (syst) events. The total expected background is 332 ± 36(stat) ± 45(syst). The background estimates are summarised in Table 1 along with the expected yield for a signal with m H = 125 GeV and B(H → inv) = 100 %.

Systematic uncertainty
The V+jets background estimates are affected by large statistical uncertainties, ranging from 5-30 %, due to control samples in data. The systematic uncertainty in the V+jet background estimates is dominated by the statistical uncertainty in the MC samples used to calculate the control-tosignal region translation factors. Additional important uncertainties arise due to jet and E miss T energy scale and resolution. These are estimated by varying the scales and resolutions associated with jets and unclustered energy within their uncertainties and recomputing the E miss T , resulting in a 13 % systematic uncertainty in the signal acceptance; 7-15 % in the V+jets background estimates; and 60 % uncer- scale uncertainties. We consider theoretical uncertainties in the vector boson fusion signal yield resulting from PDF uncertainties and factorization and renormalization scale uncertainties. The uncertainty in the gluon fusion signal yield is dominated by MC modelling of initial-state radiation, amongst other effects, and is estimated to be 60 % by comparing different MC generators. This has a modest overall effect since the gluon fusion yield is small. These uncertainties are summarized in Table 2, where they are quoted with respect to the total background or signal yield. The combined effect of all background uncertainties results in a relative increase of about 65 % in the expected upper limit on the B(H → inv).

Results
As shown in Table 1, we observe 390 events the signal region in data, compatible with the background only prediction. Figure 4 shows the E miss T and M jj distributions in data and simulated backgrounds in the signal region. The simulated V+jets backgrounds shown in this figure are normalized to the estimates from data given in Table 1.

Search strategy
The final state in the Z( )H(inv) channel consists of a pair of high-p T isolated leptons from the Z boson decay, high E miss T from the undetectable Higgs boson decay products, and limited jet activity. Since the signal cross section is orders of magnitude lower than those for inclusive DY+jets, W+jets, and tt, stringent requirements are needed to isolate the signal. We apply an event selection that is optimized for m H = 125 GeV while still being suitable for the other Higgs boson mass values considered. After this selection, the dominant backgrounds arise from ZZ and WZ processes, which are modelled using MC simulation. Smaller background contributions, from DY+jets, tt, WW, and W+jets, are modelled using control regions in data. For each value of the Higgs boson mass, the final background and signal yields used to calculate limits are obtained from a fit to the twodimensional distribution of the transverse mass, m T , of the dilepton-E miss T system, and the azimuthal separation of the two leptons.

Event selection
We use dielectron and dimuon triggers with p T > 17 GeV ( p T > 8 GeV) thresholds for the leading (subleading) lepton, together with single-muon triggers that allow recovery of some residual trigger inefficiencies. For data taken during periods when the instantaneous luminosity was low enough The E miss T (top) and M jj (bottom) distributions in data and MC after the full selection in the VBF search signal region. The simulated background from different processes is normalized to the estimates obtained from control samples in data, and shown cumulatively, with the total systematic uncertainty shown as a hatched region. Note that the QCD multijet background is not shown due to limited MC statistics, which results in a small apparent discrepancy between data and the backgrounds shown at low values of E miss T and M jj . The cumulative effect of a signal from a Higgs boson with SM VBF production cross section, m H = 125 GeV and B(H → inv) = 100 % is also shown to allow it, we also use a dimuon trigger with a p T > 7 GeV threshold for each muon.
The offline selection starts by requiring two well-identified, isolated leptons of the same flavor and opposite sign (e + e − or μ + μ − ), each with p T > 20 GeV. The invariant mass of the pair must be within ±15 GeV of the Z boson mass. To reduce the large potential background from DY( )+jets events, where the E miss T arises from mismeasurement, any event containing two or more jets with p T > 30 GeV is rejected. The remaining zero-and one-jet samples are treated separately in the analysis because of their significantly different signal-to-background ratios.
The top-quark background is further suppressed by rejecting events containing a bottom-quark decay identified by either the presence of a soft-muon or by the CSV b-tagging algorithm described in Sect. 2. The tagged b jet is required to have p T >20 GeV and to be reconstructed within the To reduce the WZ background in which both bosons decay leptonically, events containing additional electrons or muons with p T > 10 GeV are rejected. After all selection requirements, most of the remaining WZ background is from the decay mode (W → τ ν)(Z → ).
The remaining event selection uses three variables: E miss T , φ( , E miss T ), and |E miss T − p T |/ p T , where p T is the transverse momentum of the dilepton system. The last two variables effectively suppress reducible background processes like DY( )+jets and top-quark production. We optimized the selection criteria applied to these variables, in order to obtain the best expected exclusion limits at 95 % CL for m H = 125 GeV. For each possible set of selections, we repeat the full analysis, including the shape fits described in Sect. 6.5 below, the estimation of backgrounds from control data samples, and the systematic uncertainties. The final selection criteria obtained after optimization are: E miss T > 120 GeV, φ( , E miss T ) > 2.7 and |E miss T − p T |/ p T < 0.25. The efficiency of the full selection for the Z( )H(inv) signal at m H = 125 GeV is 5.6 %, estimated from MC simulation.

Background estimation
After the full selection, the dominant backgrounds arise from WZ and ZZ processes, which are modeled using MC simulation. The pre-fit normalization of these backgrounds is obtained from their respective NLO cross sections computed with mcfm.
The DY( )+jets background is modeled from an orthogonal control sample of events with a single isolated photon produced in association with jets (γ + jets). This choice has the advantage of providing a large statistics sample, which resembles Z boson production in all important aspects: production mechanism, underlying event conditions, pileup scenario, and hadronic recoil [59]. The kinematic distributions and overall normalization of the γ + jets events are matched to Z( ) + jets in data through event weights, determined as a function of the Z boson p T measured from data. This procedure takes into account the dependence of the E miss T on the associated hadronic activity.
Further discrepancies can arise due to differences in the pileup distribution of the γ + jets sample due to the fact that photon data was collected with triggers whose prescales varied as a function of photon threshold and data-taking period. These are taken into account by further weighting events in the control sample, according to the distribution of number of reconstructed vertices in the signal sample. The electroweak backgrounds to the control sample, involving photons and neutrinos, are subtracted using predictions from MC simulation.
This procedure yields an accurate model of the E miss T distribution in DY( )+jets events, as shown in Fig. 5 (left), which compares the E miss T distribution of the weighted γ +jets events, summed with other backgrounds, to the E miss T distribution of the dilepton events in data. Figure 5 also compares the distributions of (center) φ( , E miss T ) and (right) |E miss T − p T |/ p T obtained from this background model to the same distributions in the dilepton sample. The differ-ence between data and background predictions is less than 10 % in these distributions, which is negligible compared to the estimated systematic uncertainties after the final selection. The uncertainties in the electroweak background to the photon control sample yield a 100 % uncertainty in the normalization of the residual DY( )+jets background. However, since the Drell-Yan background after the full selection is very small, the large uncertainty has negligible impact on the final results.
The remaining background processes do not involve Z boson production, and are referred to as non-resonant backgrounds. Such backgrounds arise mainly from leptonic W boson decays in tt, tW decays and WW events. Also included in the estimate of non-resonant backgrounds are small contributions from single-top-quark events produced from s-and t-channel processes, W+jets production, and Z → τ τ events in which τ leptons produce electrons/muons and E miss T . We estimate these backgrounds using a control sample in data, consisting of events with opposite-charge differentflavor dilepton pairs (e ± μ ∓ ) that otherwise pass the full selection. The backgrounds in the e + e − and μ + μ − final states are then estimated by applying scale factors (α ee , α μμ ) to the number of events in the control sample, N eμ : We compute the two factors α ee and α μμ in the sidebands (SB) of the Z peak (40 < m < 70 GeV and 110 < m < 200 GeV) by using the following relations: where N SB ee , N SB μμ , and N SB eμ are the number of events in the Z sidebands counted in a top-quark-enriched sample of e + e − , μ + μ − , and e ± μ ∓ final states, respectively. The requirements for this sample are E miss The validity of the procedure for computing the scale factor is checked by closure tests on simulated samples. This method accounts for possible differences in probability for electrons and muons to pass the trigger and selection requirements. We also cross-check the methods by calculating α ee and α μμ from the Z peak region as follows: where N peak ee , N peak μμ , are the number of dielectron and dimuon events in a Z control sample. This method takes advantage of the equality between the production rates for Z → ee and Z → μμ and equates the ratio of observed dilepton counts to the square of the ratio of efficiencies. From the comparison of methods and the closure tests, we derive an uncertainty of 25 % on the normalization of the non-resonant background in addition to the contribution from the statistical uncertainties on the control samples. The background in the signal region, estimated using the methods described above, are shown in Table 3, along with the expected yield for a signal with m H = 125 GeV and 100 % invisible branching fraction. Table 4 lists the systematic uncertainties affecting this search. The most important uncertainties are those associated with theory, affecting both the signal acceptance and the dominant WZ and ZZ backgrounds. The uncertainties arising from missing higher-order QCD corrections are estimated by scaling the renormalization and factorization scales up and down by a factor of two, while those associated with PDFs are estimated using the PDF4LHC prescription [35,36].

Systematic uncertainty
The uncertainties related to jet and E miss T energy scale and resolution, lepton p T scale, and reconstruction efficiency affect the signal and all backgrounds, and are estimated as for the search in the VBF mode (see Sect. 5.4).
Uncertainties of approximately 100 %, which are derived from the data by comparing different estimation methods and conducting closure tests, are assigned to the non-resonant backgrounds. Due to the small size of the control samples, the relative uncertainties are large, but absolute contribution of these backgrounds is small.
The combined signal efficiency uncertainty is estimated to be ∼12 %, and the total uncertainty in the background estimations is about ∼15 %, dominated by the theoretical uncertainties mentioned above. The combined effect of all systematic uncertainties results in a relative increase of about 35 % in the expected upper limit on the B(H → inv). Table 3, the total number of observed events is 134 with an estimated background of about 138 events, while the expected signal yield is 35 events. The final limits on a signal are determined using a profile likelihood fit to the normalizations and the shapes of selected distributions in the signal region. For the 8 TeV data, we use the two-dimensional distribution of the azimuthal dilepton separation ( φ ) and the m T of the dilepton-E miss  one-dimensional fit to m T alone. The expected ratio of signal to background increases at high values of m T and low values of φ , giving the shape analysis greater sensitivity than a limit obtained from event counts alone. The transverse mass m T is given by the formula

As shown in
This definition of m T , which treats both the lepton pair and the recoiling undetected system as massless, is found to yield The two center-of-mass energies (7 and 8 TeV), two lepton flavors (e and μ), and two jet multiplicities (0 and 1), define eight disjoint samples that are treated separately in the likelihood calculation. The shapes and normalizations of the signal and of each background component are allowed to vary within their uncertainties, and correlations in the sources of systematic uncertainty are taken into account. The m T distribution in the 7 TeV data, and the φ distribution in the 8 TeV data, in the signal region are shown in Fig. 6 for illustration. As can be seen, the observed data are consistent with the predicted backgrounds.

Search strategy
The Z(bb)H(inv) search closely follows the strategy of the CMS search for SM Z(νν)H(bb) [60], sharing the same E miss T + bb final state, though the bb resonances have different masses. The event selection requires large E miss T , equivalent to the boost of the Higgs boson [61], and a jet pair consistent with a Z → bb decay. The signal yield after the final selection is estimated using a BDT trained on simulated background and signal MC samples, by fitting BDT output for background and signal to that obtained from data.
The backgrounds in this channel arise from production of W and Z bosons in association with jets (V+jets), tt, singletop-quark, diboson (VV), and QCD multijet production. The SM Higgs process, Z(νν)H(bb), has a negligible effect on this search, due to the different mass of the bb resonance and good di-jet mass resolution, which is about 10 %. The Z(νν)H(bb) process is therefore treated as an independent background process.
Since the VV production cross section is only a small factor larger than that of standard model VH, and given the nearly identical final state for VZ with Z(bb), the VV process has been used as a benchmark to validate the search strategy used here [60].

Trigger
A suite of four E miss T triggers is used for this search, due to the challenge of maintaining acceptance as the instantaneous luminosity increases. A trigger with E miss T > 150 GeV is used for the full 8 TeV data set. To increase acceptance at lower E miss T , we also use triggers requiring jets in addition to E miss T . For the early data-taking period, a trigger requiring E miss T > 80 GeV together with two jets with |η| < 2.5 and p T > 30 GeV was used. However, as the average instantaneous luminosity reached 3 × 10 33 cm −2 s −1 , this was replaced with a trigger requiring E miss T > 100 GeV, two jets with individual p T above 60 and 25 GeV respectively, the vector sum of the two jet p T to be above 100 GeV, and finally a veto on any jet with p T > 40 GeV and closer than 0.5 radians in φ to the E miss T direction. Finally, a trigger was used that requires E miss T > 80 GeV, together with two jets having |η| < 2.5 and p T > 20 GeV or p T > 30 GeV, depending on the luminosity conditions, and at least one of the jets tagged by the online CSV b-tagging algorithm [55].
For Z(bb)H(inv) events with E miss T > 170 GeV, the combined trigger efficiency is near 100 % with respect to the offline event reconstruction and selection, described in the next section. For events with E miss T between 130 and 170 GeV

Event selection
The event selection in this channel is designed to enhance heavy-flavor production and a Higgs boson with high Lorentz boost, with reasonable kinematic thresholds consistent with the trigger selection, and to provide sufficient statistics to perform the BDT training properly. The event selection is summarized in Table 5. Backgrounds to the signal are substantially reduced by a large E miss T requirement. In this regime, where the Higgs boson has substantial boost, the Z and Higgs bosons are separated by a large azimuthal opening angle, we therefore require φ(Z, H) > 2.0 radians. We define "low", "intermediate", significance, defined as the ratio of the E miss T and the square root of the scalar sum of transverse energy of all particle-flow objects, to be greater than three.
To reduce the tt and WZ backgrounds, events with isolated leptons with p T > 15 GeV are rejected.
The Z boson candidate is defined to be the pair of central (|η| < 2.5) jets, above minimum p T thresholds given in Table 5, that has the greatest vector sum of transverse momenta, p T jj . Each event is required to pass minimum requirements on p T jj as well as the invariant mass of the jet pair, M jj . In the low-E miss T category, events with two or more jets in addition to this pair are vetoed. Each jet in the Z boson pair are required to be tagged by the CSV algorithm. Separate thresholds are applied to the jets with higher (CSV max ), and lower (CSV min ), values of the CSV discriminator. The background from V+jets and VV processes is reduced significantly through b tagging, leaving the background in the signal region dominated by sub-processes where the two jets originate from genuine b quarks.
The Z boson mass resolution is improved by roughly 10 % by applying regression techniques similar to those used by the CDF Collaboration [62] and in the VH(bb) search by the CMS Collaboration [60]. This results in a resolution of approximately 10 %, after all event selection criteria are applied, with a few percent bias on the mass.
The selection is optimized to give the best signal significance, for a signal with m H = 125 GeV and B(H → inv) = 100 %. After all selection criteria, the efficiency for a signal with m H = 125 GeV and B(H → inv) = 100 % is 4.8 %, while for the most sensitive region of the BDT distribution, defined in Sect. 7.5, it is 1.75 %. The effect of the selection on signal and background can be seen in Fig. 7 which shows the M jj and CSV min distributions after all other selection requirements. The simulated background contributions are displayed cumulatively, and the uncertainty in the total background is shown as a hatched region. The arrows correspond to the cuts applied for the final selection as described in Table 5. The panels below both distributions show the ratio of observed data to expected background events As mentioned above, a BDT is used in the final stage of the analysis to discriminate signal from backgrounds. The BDT is trained using simulated samples for signal and all background processes after the full selection described above. This is performed separately for each Higgs boson mass hypothesis, which cover the range 105 < m H < 145 GeV in 10 GeV steps. The set of input variables to the BDT is chosen by iterative optimization from a larger number of potentially discriminating variables, and is listed in Table 6.

Background estimation
All backgrounds are modeled using MC simulation. Control regions in data are used to validate the simulated distributions used as input to the BDT. These control regions are also used to obtain scale factors to correct the pre-fit normalizations of the dominant Z+jets, W+jets and tt backgrounds. We use the same control regions as defined in Ref. [60] for the Z(νν)H(bb) search. For W backgrounds, the control region is defined using the same kinematic selection as the signal region apart from the lepton veto, which is inverted. For Z backgrounds we require a mass veto around the Higgs boson mass hypothesis. In addition we split the Z and W backgrounds into heavy-flavor enriched regions, by requiring the same b-tag as the signal region, and light-flavor enriched regions, by inverting the b-tag definition of the signal region. For the tt background, the control region is defined by inverting the lepton veto and additional jet criteria, with respect to the signal region definition.
To obtain the scale factors by which the simulated event yields are adjusted, a set of binned likelihood fits are performed to the CSV min distributions of events in the control regions. These fits are done simultaneously in all control regions, and the normalization of each background process is allowed to vary independently. Fits to several other variables are also performed, to verify consistency. The scale factors account not only for cross section discrepancies, but also residual differences in physics object selection. For the Z and W backgrounds, separate sets of scale factors are obtained for each process according to how many of the two jets selected in the Z boson reconstruction originate from a b quark. These are labelled: V+udscg for the case where none of the jets originates from a b-quark, V+b for the case where only one of the jets is from a b quark, and V+bb for the case where both jets originate from b quarks. The scale factors obtained are all close to and compatible with unity, except the V+b background where the scale factor is closer to 2, as seen in Ref. [60]. Table 7 shows the expected signal and background yields, estimated from MC simulation as described above. Figure 8 shows the distribution of CSV b-tag discriminant and dijet p T in the Z+bb and W+bb enriched regions, respectively. The high-E miss T category is shown, after the data/MC scale factors are applied. Table 8 lists the uncertainties considered in this channel. The values quoted are for the most sensitive region of the analysis (S/B > 3.5 %), which corresponds to requirements on the BDT output of >0.8, >0.7, and >0.2 in the low, intermediate, and high-E miss T categories, respectively.  Important theoretical uncertainties arise in the signal yield estimation from factorization and renormalization scales, as well as PDF uncertainties, and are estimated as for the Z( )H(inv) and VBF searches. In addition, uncertainties arising from the QCD NNLO and electroweak NLO corrections discussed in Sect. 3 are included.

Systematic uncertainty
The background estimates are unaffected by theoretical uncertainties, since they are corrected using data/MC scale factors, as discussed in Sect. 7.4. However, uncertainties in the background normalization arising from the scale factors themselves are accounted for, by propagating other system-atic uncertainties (jet energy scale, jet energy resolution, b tagging efficiency) to the control regions and repeating the fit procedure. Cross section uncertainties of 15 % each are assigned to the single-top-quark backgrounds in the t-and tW-channels, resulting in approximately 1 % uncertainty in the sum of all backgrounds. For the diboson backgrounds, a 7 % cross section uncertainty is assigned, consistent with the CMS measurement of this process [64], which results in an uncertainty of approximately 4 % in the total background.
As indicated in Table 8, uncertainties affecting the shape of the BDT output are also considered: trigger efficiency, jet These translate into yield uncertainties in the 3-5 % range, depending on the channel and the specific process. The shape of the BDT output distribution is also affected by the shape of the CSV distribution, and is therefore recomputed as the CSV distribution is varied within its uncertainties. The shape uncertainty due to MC modelling of backgrounds is estimated by comparing MadGraph and herwig++ results for the V+jets backgrounds, and comparing MadGraph with powheg for tt.
The combined effect of all systematic uncertainties results in a relative increase of about 20 % in the expected upper limit on the B(H → inv).

Results
The number of events observed in data are shown alongside the background estimates in Table 7, for the most sensitive regions of the analysis as defined in the previous section. The BDT output distributions of the three E miss T categories are shown in Fig. 9. In the Z(bb)H(inv) search, limits are determined using a fit to the BDT output distribution. This is performed separately for each Higgs boson mass hypothesis, every 10 GeV in the range 105-145 GeV. In the fit, the shape and normalization for signal and each background component are allowed to vary within the systematic and statistical uncertainties described in Sect. 7.5. These uncertainties are treated as nuisance parameters in the fit, with appropriate correlations taken into account. All nuisance parameters, including the scale factors described in Sect. 7.4 are adjusted by the fit.

Cross section limits
No evidence for a signal is observed in any of the three searches. We set 95 % CL upper limits on the Higgs boson production cross section times invisible branching fraction, B(H → inv), for the VBF and ZH production modes separately. Limits are calculated using a CL s method [65,66], based on asymptotic formulae from Ref. [67], following the standard CMS Higgs boson searches combination technique [3,68]. Systematic uncertainties are incorporated as nuisance parameters and treated according to the frequentist paradigm described in Ref. [68]. We also present 95 % CL limits on Higgs boson production cross section times invisible branching fraction normalised to the SM production cross section [38,39], which we will denote ξ = σ · B(H → inv)/σ SM . We present limits on ξ for the VBF and ZH modes separately and from the combination of all channels. It should be noted that the assumption of SM production cross sections is an arbitrary choice, as a sizeable invisible width would indicate physics beyond the SM, which may also modify the production cross-section. However, an alternative choice of Table 8 Summary of the uncertainties in the Z(bb)H(inv) channel. The numbers indicate the change in the total background estimate or in the total signal acceptance when each systematic effect is varied according to its uncertainties. Those uncertainties designated as "Norm." only affect the normalization of the contributions, while those designated "Shape" also affect the shapes of the BDT output. In the case of shape variations, the numbers indicate the range of changes across the bins of the distributions. Signal uncertainties are quoted for m H = 125 GeV and B(H → inv) = 100 %. Due to correlations, the total systematic uncertainty is less than the sum in quadrature of the individual uncertainties. The effect is evaluated in the most sensitive region of the BDT output Events / 0.1 Under the assumption of SM production cross sections and acceptances, we may interpret limits on ξ as limits on the invisible branching fraction of the 125 GeV Higgs boson. Figure 10 (top) shows the observed and median expected 95 % CL limits on the Higgs boson production cross section times invisible branching fraction, as a function of the Higgs boson mass, for the VBF production mode. Figure 10 (bottom) shows the corresponding limit on ξ . Assuming the SM VBF production cross section and acceptance, this corresponds to an observed (expected) upper limit on B(H → inv) of 0.65 (0.49) for m H = 125 GeV.
The 95 % CL observed and median expected upper limits on the Higgs boson production cross section times invisible branching fraction for the ZH production mode are shown  By assuming production cross sections as for the SM Higgs boson, the results of the three individual searches may be combined and interpreted as a limit on the invisible branching fraction of the 125 GeV Higgs boson. The statistical combination fully accounts for correlations between nuisance parameters in the individual searches. The most important correlations are unsurprisingly those associated with the signal uncertainty in the ZH searches, due to PDF and renormalization/factorization scale variation uncertainties. The most important correlated uncertainties are, in decreasing order of importance, the jet energy scale uncertainty, those associated with the signal uncertainty, due to PDF and renormalization/factorization scale variation uncer- tainties, the total integrated luminosity uncertainty, the lepton momentum scale uncertainties, the jet energy resolution uncertainty and the E miss T energy scale and resolution uncertainties. The resulting 95 % CL limit on ξ is shown in Fig. 12  and summarised in Table 9. Assuming the SM production cross section and acceptance, the 95 % CL observed upper limit on the invisible branching fraction for m H = 125 GeV is 0.58, with an expected limit of 0.44. The corresponding observed (expected) upper limit at 90 % CL is 0.51 (0.38). These limits significantly improve on the indirect 95 % CL limit of B(H → inv) < 0.89 obtained from visible decays [3].

Dark matter interactions
We now interpret the experimental upper limit on B(H → inv), under the assumption of SM production cross section, in the context of a Higgs-portal model of DM interactions [7][8][9]. In these models, a hidden sector can provide viable stable DM particles with direct renormalizable couplings to the Higgs sector of the SM. In direct detection experiments, the elastic interaction between DM and nuclei exchanged through the Higgs boson results in nuclear recoil which can be reinterpreted in terms of DM mass, M χ , and DM-nucleon cross section. If the DM candidate has a mass below m H /2, the invisible Higgs boson decay width, Γ inv , can be directly translated to the spin-independent DM-nucleon elastic cross section, as follows for scalar (S), vector (V), and fermionic (f) DM, respectively [8]: Here, m N represents the nucleon mass, taken as the average of proton and neutron masses, 0.939 GeV, while √ 2v is the Higgs vacuum expectation value of 246 GeV, and β = 1 − 4M 2 χ /m H 2 . The dimensionless quantity f N [8] parameterizes the Higgs-nucleon coupling; we take the central val-  Fig. 13 Upper limits on the spin-independent DM-nucleon cross section σ SI χ−N in Higgs-portal models, derived for m H = 125 GeV and B(H → inv) < 0.51 at 90 % CL, as a function of the DM mass. Limits are shown separately for scalar, vector and fermion DM. The solid lines represent the central value of the Higgs-nucleon coupling, which enters as a parameter, and is taken from a lattice calculation, while the dashed and dot-dashed lines represent lower and upper bounds on this parameter. Other experimental results are shown for comparison, from the CRESST [71], XENON10 [72], XENON100 [73], DAMA/LIBRA [74,75], CoGeNT [76], CDMS II [77], COUPP [78], LUX [79] Collaborations ues of f N = 0.326 from a lattice calculation [69], while we use results from the MILC Collaboration [70] for the minimum (0.260) and maximum (0.629) values. We convert the invisible branching fraction to the invisible width using B(H → inv) = Γ inv /( SM +Γ inv ), where SM = 4.07 MeV. Figure 13 shows upper limits at 90 % CL on the DMnucleon cross section as a function of the DM mass, derived from the experimental upper limit on B(H → inv) for m H = 125 GeV, in the scenarios where the DM candidate is a scalar, a vector, or a Majorana fermion.

Summary
A search for invisible decays of Higgs bosons has been performed, using the vector boson fusion and associated ZH production modes, with Z → or Z → bb. No evidence for a signal is observed in any channel. Using a CL s method, upper limits are placed on the Higgs boson production cross section times invisible branching fraction, for the VBF and ZH channels separately and combined. These results improve the exclusion in terms of σ · B(H → inv)/σ SM for m H > 113 GeV with respect to the limits obtained at LEP [11]. By assuming standard model production cross sections, and combining all channels, the upper limit on the invisible branching fraction of a Higgs boson for m H = 125 GeV, is found to be 0.58, with an expected limit of 0.44, at 95 % confidence level. These limits assume the signal acceptance of a SM Higgs boson. These constraints are more stringent than the indirect limits obtained from visible Higgs boson decays. Finally, the result is interpreted in a Higgs-portal model of dark matter [9]. Strong limits, beyond those from direct searches, are obtained on the dark matter nucleon cross section for light dark matter.