Sensitivity of future collider facilities to WIMP pair production via effective operators and light mediators Ning Zhou, David Berge, LianTao Wang, Daniel Whiteson, and Tim Tait Department of Physics and Astronomy, University of California, Irvine, CA 9697 GRAPPA Institute, University of Amsterdam, Netherlands Department of Physics, University of Chicago, Chicago, IL We present extrapolations of the current mono-jet searches at the LHC to potential future hadron collider facilities: LHC, as well as pp colliders with s = or 0 TeV. We consider both the effective operator approach as well as one example of a light mediating particle. PACS numbers: arxiv:07.57v [hep-ex 9 Jul 0 Though the presence of dark matter in the universe has been well-established, little is known of its particle nature or its non-gravitational interactions. A vibrant experimental program is searching for a weakly interacting massive particle (WIMP), denoted as, and interactions with standard model particles via some as-yet-unknown mediator. One critical component of this program is the search for pair-production of WIMPs at particle colliders, specifically pp at the LHC via some unknown intermediate state. If the mediator is too heavy to be resolved, the interaction can be modeled as an effective field theory with a four-point interaction, otherwise an explicit model is needed for the heavy mediator. As the final state WIMPs are invisible to the detectors, the events can only be seen if there is associated initial-state radiation of a standard model particle [, see Fig, recoiling against the dark matter pair. The LHC collaborations have reported limits on the cross section of pp + X where X is a hadronic jet [, 5, photon [6, 7, and other searches have been repurposed to study the cases where X is a W [8 or Z boson [9,. In each case, limits are reported in terms of the mass scale M of the unknown interaction expressed in an effective field theory [, 9, though the limits from the mono-jet mode are the most powerful [0. In this paper, we study the sensitivity of possible future proton-proton colliders in various configurations (see Table I) to WIMP pair production using the mono-jet final state. We consider both effective operators and one example of a real, heavy Z -boson mediator. ANALYSIS TECHNIQUE The analysis of jet+ E T events uses a sample of events with one or two high p T jets and large E T, with angular cuts to suppress events with two back-to-back jets (multi-jet background). The dominant remaining background is Z ν ν in association with jets, which is indistinguishable from the signal process of +jets. The estimation of the background at large E T Is problematic in simulated samples, due to the difficulties of accurately modeling the many sources of E T. The experimental results, therefore, rely on data-driven background estimates, typically extrapolating the Z ν ν contribution from Z µµ events with large Z boson p T. In this study, we begin from experimentally reported values [, 5 of the background estimates and signal efficiencies (at s = 7 TeV, L = 5 fb, E T > 50 GeV), and use simulated samples to extrapolate to higher center-of-mass q g, γ, Z, or W q FIG. : Pair production of WIMPs ( ) in proton-proton collisions at the LHC via an unknown intermediate state, with initial-state radiation of a standard model particle.
energies, where no data is currently available. At the higher collision energies and instantaneous luminosities of the proposed facilities, the rate of multi-jet production will also be higher, requiring higher E T thresholds to cope with the background levels and the trigger rates. Simulated Samples We generate signal events as well as events for the dominant background Z ν ν+jets in madgraph5 [6, with showering and hadronization by pythia [7. In each case, we generate events with zero or one hard additional colored parton and use the MLM scheme to match the matrix-element calculations of madgraph5 to the parton shower evolution of pythia. Extrapolation In the case of each potential facility, we must choose a E T threshold for the analysis. For a given threshold at a specific facility, estimating the sensitivity of the jet+ E T analysis requires the dark-matter signal efficiency an estimate of the Z ν ν+jets background the uncertainty of the Z ν ν+jets background. For the signal efficiency, we modify the reported experimental efficiency at s = 7 TeV [, 5 to estimate the efficiency at a higher E T cut. The estimated signal yield N sig ( s, L, E T > X) for a facility with center-of-mass energy s and integrated luminosity L is N sig ( ɛ ET>X s, L, E T > X) = L ɛ 0 σ( s) ɛ ET>50 where ɛ 0 is the published signal efficiency. In each case, the efficiency ɛ ET>X of a E T threshold is measured at parton-level using the simulated samples, see Fig., and the cross sections σ( s) are leading order. In the case of the background estimate, we extrapolate from the reported background estimate, denoted s=7,l=5, ET>50 Nbg, by scaling with an extrapolation factor E b : where E b is N bg ( s=7,l=5, ET>50 s, L, E T > X) = E b Nbg E b ( s, L, E T > X) = L 5 fb ɛ E T>X σ( s) ɛ ET>50 σ( s = 7) accounting for the relative efficiency of a higher E T cut and larger background cross sections at increased center-ofmass energies. The relative uncertainty on the background N is scaled from the reported s=7,l=5, ET>50 relative uncertainty, ( N bg N bg ) s=7,l=5, E T>50, using the extrapolation factor E b as given above: bg N bg ( s, L, E T > X) = ( Nbg N bg Eb N bg ) s=7,l=5, E T>50 Together, the background estimate with uncertainties and the dark-matter signal efficiencies allow us to calculate the power of the jet+ E T analysis.
Fraction of Events Z νν+j D5 +j D8 +j D9 +j 5 Missing Transverse Energy FIG. : Missing transverse energy for signal ( +jets) and background (Z ν ν+jets) samples generated at s =7,, and 0 TeV with thresholds of 80, 500, 00 and 500 GeV, with m = GeV. TABLE I: Details of current and potential future pp colliders, including center-of-mass energy ( s), total integrated luminosity (L), the threshold in E T, and the estimated signal and background yields. s [TeV ET L [fb N D5 N bg 7 50.9 7. 970 ± 60 550 00 500 00 ± 80 0 000 00 760 ± 750 000 8. 870 ± 50 0 5500 000. 6 ± 90 RESULTS FOR EFFECTIVE FIELD THEORIES Given the expected background and uncertainties, we use the CLs method [, to calculate expected 90% confidence limits on contributions from new sources. Together with the estimated signal efficiencies, we calculate cross-section limits. As the predicted cross sections depend on M, we can therefore derive limits on M, see the top panel of Figs,, and 5. These are then translated in limits on the -nucleon cross section, see the right panel of Figs,, and 5. In addition, we study the luminosity dependence of the results at s = 0 TeV, see Figs 6, 7, and 8. In Fig 9, we map to WIMP pair annihilation cross-section limits. Our predictions are compared to Fermi-LAT limits from a stacking analysis of Dwarf galaxies [9, including a factor of two to convert the Fermi-LAT limit from Majorana to Dirac fermions, and to projected sensitivities of CTA [8. RESULTS FOR ON-SHELL MEDIATORS The EFT approach is useful when the current facility does not have the necessary center-of-mass energy to produce on-shell mediators. The next-generation facility, however, may have such power. In this section, we study the sensitivity of the proposed facilities to a model in which the heavy mediator is a Z which couples to as well as q q [ 5. We generate events as before, and measure the efficiency at parton level using simulated events. The coupling of the Z is a free parameter in this theory, but particularly interesting values are those which correspond to the limit of previous facilities on M. That is, an EFT model of the Z interaction has M = g Z M Z
M D5 SI n cross section [cm 6 7 8 9 0 D5 CDMS low energy CoGeNT 0 pp0, /ab pp, /ab LHC, /ab LHC, 00/fb 5 6 7 8 9 XENON0 0 XENONT LHC, 00/fb LHC, /ab pp, /ab pp0, /ab m m FIG. : Limits at 90% CL in M (left) and in the spin-independent WIMP-nucleon cross section (right) for different facilities using the D5 operator as a function of m. fixing the relationship between g Z and M Z. Figure shows the expected limits on the Z model at a facility with s = TeV and L = 00 fb, in terms of the cross section σ(pp Z j for E T > 550 GeV and in terms of g Z. The g expected limits can be compared to the curve with g Z = M Z M ; the cross-section limits can be compared to the predicted cross section assuming g Z = M Z M. Similar results for other facilities are shown in Figures,, and. Acknowledgements We acknowledge useful conversations with Roni Harnik and Patrick Fox. DW and NZ are supported by grants from the Department of Energy Office of Science and by the Alfred P. Sloan Foundation. The research of TMPT is supported in part by NSF grant PHY-09707 and by the University of California, Irvine through a Chancellor s fellowship. [ M. Beltran, D. Hooper, E. W. Kolb, Z. A. Krusberg, and T. M. Tait, JHEP 09 (0) 07, arxiv:0.7 [hep-ph. [ P. J. Fox, R. Harnik, J. Kopp, and Y. Tsai, Phys.Rev. D85 (0) 0560, arxiv:9.98 [hep-ph. [ J. Goodman, M. Ibe, A. Rajaraman, W. Shepherd, T. M. Tait, et al., Phys.Rev. D8 (0) 60, arxiv:08.78 [hep-ph. [ G. Aad et al. [ATLAS Collaboration, JHEP 0, 075 (0), arxiv:.9 [hep-ex. [5 S. Chatrchyan et al. [CMS Collaboration, JHEP 09, 09 (0), arxiv:06.566 [hep-ex. [6 G. Aad et al. [ATLAS Collaboration, Phys. Rev. Lett., 080 (0), arxiv:09.65 [hep-ex. [7 S. Chatrchyan et al. [CMS Collaboration, Phys. Rev. Lett. 8, 680 (0), arxiv:0.08 [hep-ex. [8 Y. Bai and T. M. P. Tait, Phys. Lett. B 7, 8 (0), arxiv:08.6 [hep-ph.
5 M D8 SD n cross section [cm 6 7 8 9 0 D8 SIMPLE 0 Picasso 0 COUPP 0 IceCube bb + IceCube W W XENONT pp0, /ab pp, /ab LHC, /ab LHC, 00/fb 5 6 7 8 9 LHC, 00/fb LHC, /ab pp, /ab pp0, /ab m m FIG. : Limits at 90% CL in M (left) and in the spin-dependent WIMP-nucleon cross section (right) for different facilities using the D8 operator as a function of m. [9 G. Aad et al. [ATLAS Collaboration, JHEP 0, 8 (0), arxiv:.6096 [hep-ex. [ L. M. Carpenter, A. Nelson, C. Shimmin, T. M. P. Tait and D. Whiteson, Phys. Rev. D 87, 07005 (0), arxiv:.5 [hep-ex. [ M. Beltran, D. Hooper, E. W. Kolb and Z. C. Krusberg, Phys. Rev. D 80, 0509 (009), arxiv:0808.8 [hep-ph. [ W. Shepherd, T. M. P. Tait and G. Zaharijas, Phys. Rev. D 79, 0550 (009), arxiv:090.5 [hep-ph. [ Q. -H. Cao, C. -R. Chen, C. S. Li and H. Zhang, JHEP 8, 08 (0), arxiv:09.5 [hep-ph. [ J. Goodman, M. Ibe, A. Rajaraman, W. Shepherd, T. M. P. Tait and H. -B. Yu, Phys. Lett. B 695, 85 (0), arxiv:05.86 [hep-ph. [5 Y. Bai, P. J. Fox and R. Harnik, JHEP, 08 (0), arxiv:05.797 [hep-ph. [6 A. Rajaraman, W. Shepherd, T. M. P. Tait and A. M. Wijangco, Phys. Rev. D 8, 0950 (0), arxiv:8.96 [hep-ph. [7 R. C. Cotta, J. L. Hewett, M. P. Le and T. G. Rizzo, arxiv:.055 [hep-ph. [8 F. J. Petriello, S. Quackenbush and K. M. Zurek, Phys. Rev. D 77, 500 (008), arxiv:080.005 [hep-ph. [9 Y. Gershtein, F. Petriello, S. Quackenbush and K. M. Zurek, Phys. Rev. D 78, 09500 (008), arxiv:0809.89 [hep-ph. [0 N. Zhou, D. Berge and D. Whiteson, Phys. Rev. D 87, 0950 (0), arxiv:0.69 [hep-ex. [ A. L. Read, J.Phys.G G8 (00) 69 70. [ T. Junk, Nucl. Instrum. Methods A, 5 (999), arxiv:990006 [hep-ex. [ H. An, X. Ji and L. -T. Wang, JHEP 07, 8 (0), arxiv:0.89 [hep-ph. [ M. T. Frandsen, F. Kahlhoefer, A. Preston, S. Sarkar and K. Schmidt-Hoberg, JHEP 07, (0) [arxiv:0.89 [hep-ph. [5 I. M. Shoemaker and L. Vecchi, Phys. Rev. D 86, 050 (0) [arxiv:.557 [hep-ph. [6 J. Alwall, M. Herquet, F. Maltoni, O. Mattelaer, and T. Stelzer, JHEP 6 (0) 8, arxiv:6.05 [hep-ph. [7 T. Sjostrand, S. Mrenna, and P. Z. Skands, JHEP 0605 (006) 06, arxiv:hep-ph/06075 [hep-ph. [8 M. Doro et al., Astropart.Phys. (0) 89, arxiv:08.556 [astro-ph. [9 Fermi-LAT Collaboration, M. Ackermann et al., Phys.Rev.Lett. 7 (0) 0, arxiv:8.56 [astro-ph.he.
6 M D9 SD n cross section [cm 7 8 9 0 6 D9 SIMPLE 0 Picasso 0 COUPP 0 IceCube bb + IceCube W W XENONT pp0, /ab pp, /ab LHC, /ab LHC, 00/fb 5 6 7 8 9 LHC, 00/fb LHC, /ab pp, /ab pp0, /ab m m FIG. 5: Limits at 90% CL in M (left) and in the spin-dependent WIMP-nucleon cross section (right) for different facilities using the D9 operator as a function of m. M D5 6 SI n cross section [cm 7 8 9 0 D5 CDMS low energy CoGeNT 0 pp0, e6/pb pp0, e6/pb pp0, e5/pb pp0, e/pb pp0, e/pb pp0, e/pb pp0, e/pb 5 6 7 8 9 XENON0 0 XENONT pp0, e/pb pp0, e/pb pp0, e/pb pp0, e/pb pp0, e5/pb pp0, e6/pb pp0, e6/pb m m FIG. 6: Limits at 90% CL in M (left) and in the spin-independent WIMP-nucleon cross section (right) for increasing luminosity using the D5 operator as a function of m.
7 M D8 6 SD n cross section [cm 7 8 9 0 D8 Picasso 0 COUPP 0 SIMPLE 0 IceCube bb + IceCube W W XENONT pp0, e6/pb pp0, e6/pb pp0, e5/pb pp0, e/pb pp0, e/pb pp0, e/pb pp0, e/pb 5 6 7 8 9 pp0, e/pb pp0, e/pb pp0, e/pb pp0, e/pb pp0, e5/pb pp0, e6/pb pp0, e6/pb m m FIG. 7: Limits at 90% CL in M (left) and in the spin-dependent WIMP-nucleon cross section (right) for increasing luminosity using the D8 operator as a function of m. M D9 6 SD n cross section [cm 7 8 9 0 D9 Picasso 0 COUPP 0 SIMPLE 0 IceCube bb + IceCube W W XENONT pp0, e6/pb pp0, e6/pb pp0, e5/pb pp0, e/pb pp0, e/pb pp0, e/pb pp0, e/pb 5 6 7 8 9 pp0, e/pb pp0, e/pb pp0, e/pb pp0, e/pb pp0, e5/pb pp0, e6/pb pp0, e6/pb m m FIG. 8: Limits at 90% CL in M (left) and in the spin-dependent WIMP-nucleon cross section (right) for increasing luminosity using the D9 operator as a function of m.
8 /s qq [cm 95% CL limit on <σ v> for 0 9 5 LHC, 00/fb LHC, /ab pp, /ab pp0, /ab x FermiLAT bb D5 CTA Halo CTA Fornax /s qq [cm 95% CL limit on <σ v> for 9 0 5 LHC, 00/fb LHC, /ab pp, /ab pp0, /ab x FermiLAT bb D8 CTA Halo CTA Fornax 6 value CTA Segue 6 value CTA Segue 7 7 8 8 9 m 9 m FIG. 9: Limits at 95% CL on WIMP pair annihilation for different facilities using the D5 (left) or D8 (right) operator as a function of m. TeV, 00/fb TeV, 00/fb Z g σ(pp Z ) [fb Limit m =0 Theory m =0 Limit m =00 Theory m =00 Limit, m =0 Fixed M^, m =0 Limit, m =00 Fixed M^, m =00 0 00 000 000 000 5000 6000 7000 0 00 000 000 000 5000 6000 7000 FIG. : Sensitivity at s = TeV, L = 00 fb to a dark matter pairs produced through a real Z mediator. Left, expected limits on the coupling g Z versus Z mass for two choices of m for events with E T > 550 GeV; also shown are the values of g Z which satisfy g /m Z = /M, where M are limits from s = 7 TeV, L = 5 fb. Right, production cross section as a function of Z mass, compared to expected limits, where g Z depends on m Z as in the left pane.
9 Z g TeV, 000/fb σ(pp Z ) [fb TeV, 000/fb Limit m =0 Theory m =0 Limit m =00 Theory m =00 Limit, m =0 Fixed M^, m =0 Limit, m =00 Fixed M^, m =00 0 00 000 000 000 5000 6000 7000 0 00 000 000 000 5000 6000 7000 FIG. : Sensitivity at s = TeV, L = 000 fb to a dark matter pairs produced through a real Z mediator. Left, expected limits on the coupling g Z versus Z mass for two choices of m for events with E T > 0 GeV; also shown are the values of g Z which satisfy g /m Z = /M, where M are limits from s = TeV, L = 00 fb. Right, production cross section as a function of Z mass, compared to expected limits, where g Z depends on m Z as in the left pane. Z g TeV, 000/fb σ(pp Z ) [fb TeV, 000/fb Limit m =0 Theory m =0 Limit m =00 Theory m =00 Limit, m =0 Fixed M^, m =0 Limit, m =00 Fixed M^, m =00 0 000 000 60008000000000000600080000000 0 000 000 60008000000000000600080000000 FIG. : Sensitivity at s = TeV, L = 000 fb to a dark matter pairs produced through a real Z mediator. Left, expected limits on the coupling g Z versus Z mass for two choices of m for events with E T > 750 GeV; also shown are the values of g Z which satisfy g /m Z = /M, where M are limits from s = TeV, L = 000 fb. Right, production cross section as a function of Z mass, compared to expected limits, where g Z depends on m Z as in the left pane.
Z g 0 TeV, 000/fb σ(pp Z ) [fb 0 TeV, 000/fb Limit m =0 Theory m =0 Limit m =00 Theory m =00 Limit, m =0 Fixed M^, m =0 Limit, m =00 Fixed M^, m =00 0 000 0000 0000 0000 50000 60000 70000 80000 0 000 0000 0000 0000 50000 60000 70000 80000 FIG. : Sensitivity at s = 0 TeV, L = 000 fb to a dark matter pairs produced through a real Z mediator. Top, expected limits on the coupling g Z versus Z mass for two choices of m for events with E T > 5500 GeV; also shown are the values of g Z which satisfy g /m Z = /M, where M are limits from s = TeV, L = 000 fb. Bottom, production cross section as a function of Z mass, compared to expected limits, where g Z depends on m Z as in the top pane.