Probing the explanation of the muon (g-2) anomaly and thermal light dark matter with the semi-visible dark photon channel.

We report the results of a search for a new vector boson ( A ' ) decaying into two dark matter particles χ 1 χ 2 of different mass. The heavier χ 2 particle subsequently decays to χ 1 and an off-shell Dark Photon A ' ∗ → e + e - . For a sufficiently large mass splitting, this model can explain in terms of new physics the recently confirmed discrepancy observed in the muon anomalous magnetic moment at Fermilab. Remarkably, it also predicts the observed yield of thermal dark matter relic abundance. A detailed Monte-Carlo simulation was used to determine the signal yield and detection efficiency for this channel in the NA64 setup. The results were obtained re-analyzing the previous NA64 searches for an invisible decay A ' → χ χ ¯ and axion-like or pseudo-scalar particles a → γ γ . With this method, we exclude a significant portion of the parameter space justifying the muon g-2 anomaly and being compatible with the observed dark matter relic density for A ' masses from 2 m e up to 390 MeV and mixing parameter ε between 3 × 10 - 5 and 2 × 10 - 2 .

Despite its great success, the Standard Model (SM) does not provide a complete description of nature. For example it cannot explain the origin of dark matter, the neutrino masses and the baryon asymmetry problem. Furthermore, interesting discrepancies between some SM predictions and measurements have been observed. These include the LHCb results challenging lepton universality [1] and the long standing discrepancy of the muon anomalous magnetic moment [2] which was recently confirmed [3]. The combination of the Brookhaven and Fermilab muon g-2 results compared to the latest theoretical calculations using dispersion relations [4] leads to a discrepancy of :It should be noted that when compared to the latest QCD lattice calculations this is reduced to about 1.4 [5]. In order to help elucidate the origin of this difference, a new experiment aiming to measure the contribution of hadronic corrections is being prepared at CERN [6].

In terms of new physics, among many interesting proposals, a way to explain this discrepancy is to introduce a 1-loop correction involving a dark sector massive gauge field to the QED 3-point vertex [7]. The Dark Photon can couple to both charged Dark Matter (DM) fields with coupling strength , and to SM leptons via kinetic mixing with the SM photon field A. Considering models with a diagonal coupling to DM and SM fields, two decay modes are possible, (visible mode) [8] and (invisible mode) [9]. These were excluded as explanations of the anomaly by the combined results of NA64 and BaBar [9-11] (for ) and NA48 [12] (for ). In addition, the prospected sensitivity for NA64 running in muon mode to probe the anomaly was estimated recently [13] in the anomaly free gauge extension of the SM.

In this work we consider an alternative model for decay. This was initially proposed to recover an explanation of the discrepancy still within the Dark Photon paradigm [14, 15]. This model is characterised by the decaying into a heavier and a lighter Dark Matter states. While is a non-interacting stable state which determines the DM relic abundance, is unstable and de-excites by emitting a and an off-shell , that subsequently decays into an electron-positron pair. Thus, a new semi-visible decay mode combining characteristics of both the visible and invisible decay signatures emerges. It is very remarkable that such a model can potentially explain both and the Dark Matter thermal relic abundance for 300 MeV 1 GeV, thus making it of great phenomenological interest [14, 15].

In this study, we focus on a Dark Matter model that extends the SM symmetry group with a dark sector , which is spontaneously broken by a dark Higgs field [15]. The gauge mixing with the SM photon via the term , governed by the parameter , generates a massive Dark Photon after spontaneous symmetry breaking.

The diagonalisation of the mixed gauge interaction [16] allows the removal of from the pure mixed-gauge Lagrangian. However, a new coupling appears between and the SM electromagnetic current, with interaction strength . The novel interaction term is the main phenomenological feature of the model: the possibility to have SM final states produced via kinetic mixing of with A. Nevertheless, the unique feature of the semi-visible model, that distinguishes it from the invisible and visible channels, is found in the Yukawa dark sector. In the unbroken theory, we start from a pseudo-Dirac field charged under . The chiral projections of can then be coupled with to produce Dirac and Majorana mass terms, after spontaneous symmetry breaking [15].

In this pseudo-Dirac scenario, the right and left Majorana masses are the same and are strongly suppressed relative to the Dirac mass . The diagonalisation of the mass matrix results in two eigenstates: a lighter stable and a heavier DM particle, with mass difference [17].

The effective Lagrangian for the semi-visible model is:In Eq. 1, the coupling to muons is neglected since the di-muon production threshold is not relevant for the sub-GeV mass range of probed in this study. The absence of elastic diagonal interaction terms derives from the choice of a pseudo-Dirac field, where only off-diagonal terms are permitted. This allows us to neglect the invisible and visible decay modes, where conversely only diagonal terms are present. A sketch of the dominating decay chain is shown in Fig. 1. The Dark Photon decays promptly in a lighter and a heavier via the mentioned inelastic interaction, followed by the subsequent decay .

Fig. 1

Production of and subsequent semi-visible decay chain of a Dark Photon,

The width of the process was calculated at leading order. A numerical approach was used to compute the 3-body decay phase-space, implemented in the module MadWidth of MadGraph5aMC@NLO [18]. Thus, a correction is obtained to the previous analytical approximations (valid for ) of the formula from Refs. [15, 19]. The newly attained formula readswhere is the correction factor estimated using both Madgraph and CalcHep [20], found to be in good agreement. The quoted uncertainty on the K factor takes into account the difference between Madgraph and CalcHep and the slight dependence of K on . The K-factor was found to be basically insensitive to the other parameters of the model.

An upper bound for the dark sector coupling constant can be found by requiring the absence of a Landau pole for the effective coupling constant up to an energy scale  TeV: [21, 22]. In this study, a benchmark value of is used, compatible with other literature [14, 15]. Nevertheless, a discussion on the implications of different choices in our results will be provided. Furthermore, the resonance of the thermal averaged non-relativistic co-annihilation DM cross section present at [17, 19] can be avoided by setting as benchmark such as in [15, 17, 19]. Finally, the parameter has kinematic limits and . A relatively large mass splitting is chosen in this study, as strong bounds for lower already exist as explanation of by BABAR and E137 [10, 14, 15]. A complete discussion of the achievable range, up to the limit , is provided below.

In this work, we present a direct search for the semi-visible signature using the NA64 experiment located at CERN SPS. The Dark Photons are produced in the process as 100 GeV electrons coming from the H4 beamline scatter inside the NA64 electromagnetic calorimeter (ECAL). The production mechanism is identical to the one described in Ref. [9]. The setup is schematically shown in Fig. 2. The experiment uses a set of scintillator and veto counters, a magnet spectrometer consisting of two dipole magnets, and a set of tracking detectors (six micromegas chambers [23], three straw detectors [24] and two GEMs [25]) to define the incoming beam. A synchrotron radiation detector (SRD) is used to suppress the hadron contamination in the beam. The electrons are absorbed in a lead-scintillator sandwich Shashlick-type ECAL target of 40 radiation lengths. Downstream from the ECAL, a large high-efficiency VETO counter and three iron hadronic calorimeters (HCALs) complete the setup. An additional HCAL module is placed along the unbent beam path to further suppress background from upstream interactions. Further details about the setup can be found in Refs. [9, 26, 27].

Fig. 2

A schematic view of the NA64 invisible mode setup used in 2018

The is produced in the target via Dark-Bremsstrahlung [28] and decays promptly into . The long-lived travels through the VETO and HCAL1, which acts as a veto to reject particles leaking from the target, before decaying through . The result of this decay would be observed in the experiment through two characteristic signatures. The first (S1), is identified by the presence a significant energy deposition in HCAL2 or HCAL3, as the result of the decay to (see Fig. 3). In the second signature (S2), decays outside the fiducial volume limited by the last HCAL module, traveling a distance m. All energy of the produced is undetected and the signature is effectively equivalent to the one used for the search (invisible mode). These two signatures closely resemble the ones used in the previous NA64 analysis searching for pseudo-scalar and axion-like particles (see more details in Refs. [9, 26]). Since the data has already been unblinded, we chose a conservative approach and re-cast this analysis using the pseudo-Dirac model instead of performing a new one. We checked that the optimization of the selection criteria will enhance the sensitivity up to around 10%. This will be used in the future to perform a blinded analysis when new data will be available. In the following, we provide a brief description of the method.

Fig. 3

A schematic view of an event from a produced after a 100 GeV scatters off in the active dump, . The particle decaying within HCAL2 corresponds to the S1 signature (see text for more details)

Impurities from the beam, consisting mostly in (1%) are suppressed using the SRD to a level of [29]. The tracking system measures the incoming energy of the electrons which is required to be within 3 GeV from the nominal 100 GeV beam energy. The ECAL serves as an active dump measuring the energy deposition of the incoming particles. We require no activity in the VETO and the first HCAL module to reject any initial beam particles that penetrate the ECAL. For S2, no activity is required in all HCAL modules, since is assumed to decay outside the fiducial volume. Finally, for S1, a cut is applied on the variable , defined as the ratio between the energy deposit in the periphery and the total energy deposited in the HCAL2 and HCAL3 modules.

The leading background of S1 is caused by traveling undetected through the first HCAL module and deposit energy in either HCAL2 or HCAL3. The S2 signature can instead be misidentified due to large missing energy produced by electron-hadron interactions along the beamline. The total expected background for the two signatures is detailed in the previously published analyses, where it was determined that 0.19±0.07 events are expected in the signal box of S1 [26] and 0.53 ± 0.17 are expected for S2 [9, 11]. A full discussion of the uncertainties is also found in these references. In particular, the two leading contributions are for the parametrisation of the form factor in the cross-section and an additional coming from the data-MC discrepancy in the dimuon-yield (events where the interaction is detected) [11, 26]. The difference in the estimate of the width, between MadGraph and CalcHEP, leads to a negligible uncertainty . All uncertainties, summed in quadrature, do not exceed 20%.

The background for this new search is the same as the previous searches for and . Thus, the applied cuts are already optimized for the best coverage of . An exception is the cut applied on the variable R. A larger tail for high values of R is expected in this model, for two reasons: (i) smaller energy due to the emission in the original prompt decay, and (ii) the three-body decay of which increases the phase space of the final decay products . To be conservative, the cut of was not changed, and the task of its optimization is left for a future analysis on new data. Instead, its effect was calculated with our detailed MC simulation. We found that on average this cut has an efficiency of 52%, weakly dependent on the parameter of the model, with a minimum value of 45%.

The signal yield was calculated using a full MC simulation based on the Geant 4 toolkit [30]. The framework used for the previous NA64 analysis of the 2016–2018 data [9, 26] was extended to include the new model containing the semi-visible decay. Both and particles are assumed to have no interaction with the detectors. For the S1 signature, is forced to decay inside the fiducial volume, i.e., the space between the beginning of HCAL2 and the end of HCAL3. The event is then weighted by the probability of such decay to take place. In the case of S2, we assume that the full energy of is lost, and a weight corresponding to the probability for the Dark Photon to decay beyond all NA64 subdetectors is assigned to the event. This simulation is performed for a grid on the () plane to estimate the expected number of events for different masses and mixing strengths. We use , and a mass splitting as benchmark for these simulations.

The exclusion limit was calculated using the multibin limit setting technique with the modified frequentist approach (the code based on the ROOSTATS package [31]), using the profiled likelihood as a test statistic [32, 33]. The corresponding 90% exclusion limit was obtained using Eq. 4 of [26] to compute the expected signal yield. The results are summarized in Fig. 4 in the 2D plane (;), where the relevant estimated bound of E137 and Babar are also shown [14, 15] together with the projection of Nucal [34, 35] and CHARM [36, 37]. Using the benchmark values discussed above, our data exclude the explanation of the semi-visible model up to a mediator mass GeV including a so-far uncovered area close to this boundary. Even though this might look only like as a slight improvement of the current bounds, we would like to stress that these are the first experimental limits in this region of parameter space obtained with a full analysis of the data including all efficiencies and uncertainties. In Fig. 5, the region of the parameter in the central band of the anomaly is also shown in the two planes (;) and (;). To project this space we use the same assumption as in Fig. 4 for the parameters not considered, and we set as the epsilon compatible with the central band of , a convention used in previous studies of this model [15, 37]. Our results exclude the unexplored area for GeV in the band, leaving space for models in which is larger than 0.4 GeV, or with large mass splitting that can still explain the anomaly. The largest limitation to probe the missing region comes from the increasingly short decay time of , which makes the chance of detection vanishingly small. This is particularly relevant for the mass splitting (see Fig. 5 bottom), since the scaling of the decay width adds a large suppression to the signal yield. For , our limits become rapidly weak.

Fig. 4

The NA64 90% exclusion limit for a new vector boson with a coupling to electrons with decay mode . The limits were derived in the (, ) assuming , and a mass splitting . The red band shows the region of parameter space within two sigma from the world average of [3]. The blue band shows the same region before the results at Fermilab were published. Constraints from BABAR and E137 are also shown following the recasting done in Refs. [14, 15], together with the bounds of NuCal and CHARM [37]. A thick black line shows the combination of parameters compatible with a DM thermal relic scenario

Fig. 5

The NA64 90% exclusion limit for a new vector boson with a coupling to electrons with decay mode . The limits were derived in the (;) (top) and (;) (bottom) plane assuming , and a , where is the value in the central band of the anomaly. Constraints from Babar and E137 are also shown following the recasting done in Refs. [14, 15], together with the bounds of NuCal and CHARM [37]. A thick black line shows the combination of parameters compatible with a DM thermal relic scenario. The projected limit for EOT using a compact HCAL1 are drawn in the (;) plane

In order to increase the sensitivity for short-lived , the HCAL1 length should be shortened as much as possible while keeping the residual background under control. The use of an absorber with very short nuclear interaction length (Tungsten or Uranium compared to the stainless steel used in the current HCAL1 ) and an optimization of the ratio between active and passive material could bring to an HCAL1 shorter by a factor of three with respect to the current length (50 cm versus 150 cm). A similar study is planned for optimizing the design of the Tungsten target used for searching for into visible final states [38]. With a shorter HCAL1 the interesting parameter space that can account for the discrepancy and saturate the thermal relic density would be completely covered in all the dimensions of the model. A more compact HCAL1 would be also beneficial to improve the NA64 sensitivity of the search for axion-like particles (ALPS) into two photons, , for larger ALP-photon couplings."

In this work, we analysed the data collected by the NA64 experiment during three different runs in the “invisible-mode” configuration considering a new pseudo-Dirac scenario characterized by the decay as signal candidate. In this model, the decay of the mediator results in both SM and DM particles in the final states, for an effective signature that combines features of both invisible and visible mode. This scenario can provide an explanation to the anomaly, recently confirmed at Fermilab [3], and at the same time is compatible with a freeze-out scenario capable to explain the observed DM relic-density. The previous limits on this model were improved by this analysis, excluding GeV at 90% confidence level (C.L.), assuming a DM coupling and a mass splitting . A large region of parameter space characterized by short living remains unexplored as an exciting prospect for future searches.