Abstract

Highly energetic proton/electron beam fixed-target experiments extend an opportunity to probe the sub-GeV dark matter and associated interactions. In this work, we have explored the sensitivity of DUNE (Deep Underground Neutrino Experiment) for sub-GeV leptophobic dark matter, i.e., this dark matter barely couples with the leptons. Baryon number gauge theory can predict the existence of leptophobic cold dark matter particle candidates. In our work, the dark matter candidate is considered to be scalar whose mass is defined by the symmetry breaking of new baryonic gauge group . In this scenario, a light scalar dark matter couples with the standard model candidates via vector boson mediator which belongs to the baryonic gauge group . This leptophobic dark matter dominantly couples to the quarks. Under this scenario, new parameter space for is explored by DUNE for leptophobic dark matter candidates. This new parameter space allowed to get a lower value than the present exiting constraint value of , i.e., .

1. Introduction

Existence of dark matter (DM) particle is supported by many gravitational phenomena in astrophysics and cosmology, but the nature of the dark matter is still to be discovered. These gravitational evidences motivate the physicist to know about the dark matter (DM) particle candidates and their interactions with ordinary matter if any. None of the running experiments is able to probe the DM candidates; hence, different theoretical assumptions or models are floating in the scientific community to describe their nature. DM candidates are assumed to be nonbaryonic, weakly interacting, and stable.

On the basis of DM mass and their velocities, DM candidates can be broadly classified into two categories: hot and cold dark matter. All relativistic and superrelativistic DM particles are called hot DM, whereas all nonrelativistic DM particles are known as cold DM. Most part of the DM is cold because they move with nonrelativistic speed, and they are capable of forming small structure galaxies rather than massive galaxies. While hot DM cannot form small structure of galaxies by the free streaming process. As we observe small structure galaxies, hence, it is assumed that DM is a mixture of small amount of hot DM with predominant cold DM. Different DM physics models need different experimental environment for the detection of DM. The experimental techniques proposed for dark matter detection can be broadly classified into three categories.

1.1. Direct Detection

In the direct detection technique, we suggest that the DM candidates can be detected by the measurement of recoil energy of nucleons, but this detection technique becomes less sensitive for the detection of sub-GeV dark matter because in this mass range of DM the recoil energy of nucleons decreases below the detector threshold value of the detector.

1.2. Indirect Detection

In indirect detection, the standard model (SM) particles produced via annihilation of DM particles are studied.

1.3. Collider Technique

Beam collider technique looks for the missing transverse energy in an event, for the DM detection.

In this work, we have focused on the MeV to GeV mass range of DM because a lot of work has been done for mass range 1 GeV to 10 TeV (WIMP (weakly interactive massive particles) as a DM candidate) at direct detection experiments, but these experiments have estimated null results for WIMP [15]. Hence, sub-GeV mass DM particles can be potential candidates of the DM family and the sub-GeV DM candidates can be explored by direct detection, in neutrino laboratories. To improve the sensitivity of DM candidates, fixed-target approach is considered in direct detection experiment, but one major drawback of this approach is the presence of huge neutrino background. These particles (neutrinos) can mimic the DM signatures hence for the detection of DM neutrino background reduction becomes necessary. To mitigate the neutrino background in our DM detection work, we have considered DUNE running in beam dump mode. For the first time, neutrino experiment in the beam dump mode was carried by MiniBooNE in 2014 [6, 7]. The fixed-target experiment in beam dump mode opens the door for the study of DM and hidden sector physics via different DM portals.

A large number of DM studies [6, 814] have focused on the kinetic mixing scenario, in which DM of hidden sector couples with the SM particles via kinetic mixing parameter , which kinetically mixes the dark photon mediator of hidden sector and ordinary photon. This is one of the possible ways of interaction of DM with SM particles, but there are other large number of scenarios of DM interaction with SM particles. These scenarios need to be probed for imposing better constraints on DM parameter space. In one of these scenarios, DM candidates of hidden sector interacts with the SM particles via a new massive vector gauge boson mediator of new baryonic gauge group . This vector gauge boson dominantly couples with the quarks and is called leptophobic dark matter.

In our work, using fixed-target experiment in beam dump mode, we have estimated the sensitivity of DUNE (Deep Underground Neutrino Experiment) [15] for capturing the signatures of leptophobic dark matter. A highly energetic proton beam of energy 120 GeV is used to produce a boosted dark matter beam through vector boson mediator . Here, three modes of DM production are taken into consideration in our analysis. These modes are pseudomeson decay, bremsstrahlung, and parton-level production channels. The DM signatures are captured in DUNE near detector by looking at DM-nucleon elastic scattering.

This paper is organized as follows. In Section 2, we have described the leptophobic DM model considered in our analysis which is followed by a note on DUNE near detector. In Section 3, we have discussed different modes of DM production. Section 4 takes up the expression for scattering cross-section of leptophobic DM. We have briefly summarized all the essential parameters that are used in our analysis and discussed about the expression of DM signal rates in Section 5. Section 6 illustrates several experimental and theoretical constraints which are relevant for leptophobic sub-GeV DM. In Section 7, we have briefly discussed the BdNMC simulation tool which is used for DM simulation in our analysis, and sensitivity results of DUNE for leptophobic sub-GeV DM are discussed. At the end, we have concluded our work, in Section 8.

2. Leptophobic Dark Matter Model and DUNE near Detector

The model considered for sub-GeV leptophobic scalar DM candidates, which are charged under a new baryonic gauge group , couples with the SM particles through a new vector boson mediator which dominantly interacts with quarks. We have used generation independent coupling of with quarks to make the model simple as it avoids the tree level flavor changing neutral current interactions. In this consideration, to allow the renormalizable Yukawa couplings of quarks to the Higgs boson of SM, coupling of with right and left-handed quarks should be the same. Our considered model, due to the addition of new gauge group to the standard model gauge group, suffers from gauge anomalies [16, 17] and can be considered as nonrenormalizable effective field theory with a cutoff [18]. New states, either at or below this cutoff (), must be introduced for the theory to remain consistent. The simplest choice for the addition of new states is new chiral heavy fermion which helps in the cancellation of the anomaly.

Here, in this work, our focus is on GeV scale phenomenology; hence, the exact details of the UV completion can be ignored, and we can focus our attention on low-energy effective field theory of a local symmetry under which the DM is charged. The Lagrangian of leptophobic scalar DM for low-energy effective field theory can be expressed as where is the field of SM, is the field strength of , , is the baryonic gauge coupling, and is the baryonic charge of . The represents the sum of baryonic current over all quark species, i.e., . Above Lagrangian includes baryonic coupling () scenario as well as kinetic mixing () interaction scenario. Two different scenarios at a time are possible and may give some interesting results, but for most cases, either the baryonic coupling scenario or the kinetic mixing scenario will dominate. Therefore, for leptophobic DM scenario, we have set the value of kinetic mixing parameter as we want to check constraint on DM parameter space in the presence of baryonic coupling.

The search of MeV scale dark matter in beam dump mode was first performed at the near detector of MiniBooNE, and this search developed the hope for DM searches at fixed-target experiments. The MiniBooNE experiment used 8 GeV accelerated proton beam to produce the DM beam, but the DM signal decreased after the mediator mass became greater than 1 GeV. The upcoming DUNE will use 120 GeV proton beam produced at Fermilab main injector for physics studies of neutrino oscillation hence providing higher energies for the study of dark matter. The far detector of DUNE is proposed to study precise neutrino oscillation physics, whereas the near detector will be capable of probing physics beyond the standard model, sterile neutrinos, nonstandard interaction, low mass DM, etc.

In this manuscript, we have focused on the possibility of upcoming DUNE to probe the dark matter sensitivity. The 120 GeV proton beam after hitting the Beryllium target produces charged mesons and neutral mesons. The neutral mesons further decay into a pair of DM via baryonic mediator . In beam dump mode, charged mesons get absorbed into thick steel pipe before they can decay into neutrinos. The produced dark matter will hit the 1 ton near detector made up of with the ratio of (: ) which will be placed at a distance of 574 m downstream from the center of the Beryllium target. This detector is a 5-m-long cylinder with 5 m diameter [19]. We have considered efficiency of DUNE near detector for sub-GeV dark matter candidates [20].

3. Production of Leptophobic Dark Matter

We have focused on three production channels of leptophobic dark matter for DUNE. The considered production modes are as follows: (i)Pseudoscalar meson decay(ii)Bremsstrahlung process(iii)Parton-level production

3.1. Pseudoscalar Meson Decay

This production channel dominates overall production modes of dark matter at lower masses of vector boson mediator . In this mode, vector boson mediator is produced via radiative decay of pseudoscalar mesons [8, 9]. These secondary mesons are produced from the primary interaction of and .

The pseudoscalar mesons coupling with vector boson takes place under the gauged Wess-Zumino-Witten (WZW) lagrangian [2123]. If the mass of the secondary meson is greater than the mass of the vector boson mediator , i.e., or if the mass of DM particle produced is such that , then will be produced on-shell, and further, it will decay into a pair of DM candidates. Branching ratio of pseudoscalar meson decay to dark matter particles is calculated using narrow width approximation [9] which is equal to the product of mesons decay to and decay of to dark matter. Where is the baryonic coupling constant, is the electromagnetic coupling constant, is the mass of the pseudoscalar mesons, and the value of is different for different mesons, i.e., and [24]. We have used the BMPT (Beryllium Material Proton Target) distribution fits [25] in rejection sampling to simulate the momentum and angular distribution of mesons. For the detailed calculation, check the reference [26].

3.2. Bremsstrahlung Process

For intermediate masses of vector boson mediator , resonant vector mesons decay by the bremsstrahlung process dominates over other decay processes. Here, the mass of is close to the mass of resonant vector mesons.

In this process, vector mesons mix with vector boson mediator . This process generates a nearly collimated beam of vector boson . The four momentum of incident proton of mass is where . The four momentum of outgoing vector boson mediator of mass is where , , and is a fraction of proton beam momentum carried away by outgoing vector boson in the direction of proton beam. Here, and are transverse and longitudinal components of the momenta of .

By Weizscker-Williams approximation, the rate of production of vector boson production per proton is as follows [14, 27, 28]: here, , is a function of atomic number , and is a splitting weight-function of photon which relates before and after differential scattering cross-section [27] here, .

Since radiative has time-like momentum and time-like form factor expresses off-shell mixing of vector bosons with vector mesons in appropriate kinematic region. Baryonic vector portal considers both proton form factor and neutron form factor [29]. These incorporate only isoscalar Breit-Wigner components [29] -like in the spacelike regime, and this form factor is not completely resolved for . Above 1 GeV, the form factors suppresses the rate of production of virtual bosonic mediator; hence, above this energy, direct parton-level production dominates over other channels.

To calculate the dark photon production rate, equation (5) must be integrated over and in a range that satisfies some kinematic conditions [27] expressed as

Using the above equation, and are selected for DUNE in the present work.

3.3. Parton-Level Production of Dark Matter

Above 1 GeV of vector boson mediator mass, this channel becomes significant. This process works under the narrow width approximation via and can be written as

Dark matter pair production cross-section at parton-level can be expressed as where cross-section for the production of vector boson and can be written as where , is the hadron-level center of mass energy, and is quark charge in the unit of positron electric charge. To calculate the cross-section of the DM production, we have used CTEQ6.6 PDFs [30] and have set which is allowed to vary from to . In above equation, is the parton distribution function (PDF) which gives the probability of extraction of quarks and gluons with longitudinal momentum fraction from a proton (neutron). Details of the cross-section calculation are discussed in the references [10, 14, 26].

4. Scattering Cross-Section of Leptophobic Dark Matter

In the considered leptophobic dark matter model, dark matter dominantly interacts with the quarks, whereas it does not couple with the leptons. In our work, we have focused on the neutral current elastic scattering of dark matter with nucleons present in the DUNE near detector. The differential cross-section of neutral current DM-nucleon elastic scattering which is similar to the neutrino-nucleon neutral current scattering [31, 32] can be expressed as where and represent energy of the incoming and outgoing DM, is the mass of nucleons , and is the momentum transfer. The and are monopole and dipole form factors [26], and the values of kinematic functions for complex scalar DM are listed as

5. Signal Rates

To simulate the DM event rates, it is essential to incorporate the relevant cuts of detector geometry and energy resolution in the simulation tool. In this work, the detector limitations are listed below [14]:

The expression for the DM events produced by the pseudoscalar mesons and vector mesons can be stated as where is the atomic number density of detector material, is the detector efficiency, is the total number of mesons produced in the target, is the total number of DM trajectories produced by relevant production channels, and is the length of DM trajectory in the detector. The DM-nucleon elastic scattering cross-section is defined as where represents the energy of incoming dark matter, represents the energy of outgoing dark matter, and is the minimum and maximum energy of outgoing DM which is calculated by the relevant experimental cuts that are derived by the experimental data of nucleon recoil momentum (). In baryonic vector portal, we take ( is atomic number) for the elastic or quasielastic scattering of DM. For parton-level production channel, we have substituted (total number of produced vector bosons) in place of in the above equation (13). Total DM events are evaluated by adding the DM events produced via three different channels considered in this work.

6. Constraints on the Leptophobic Dark Matter

In this paper, we are checking the sensitivity of DUNE fixed-target experiment in beam dump mode for leptophobic DM. Several essential constraints that need to be considered for checking the DUNE sensitivity for leptophobic DM are listed below in brief.

6.1. Thermal Relic Dark Matter

In large class of models where DM is considered to be a thermal relic, its cosmological abundance is calculated by the correct thermal relic density of DM (% of the energy density of the universe) which can be obtained by the measurement of s-wave annihilation cross-section during freeze-out. The s-wave annihilation of DM into charged SM particles, in particular, rules out DM masses below 10 GeV, but by adding the p-wave process into the annihilation cross-section of DM, the DM of masses less than 10 GeV can be produced. For the correct reproduction of DM annihilation cross-section with only s-wave process, the most viable annihilation mode of DM is its annihilation into neutrino like states. In this annihilation process, because of the weak interaction of baryonic neutrinos with matter, they are unable to ionize the hydrogen and helium gases; hence, the problem of energy injection around redshift is completely ignored, during and after the recombination. In this scenario, the required annihilation cross-section of thermal relic DM (of order) can be achieved. The required value of annihilation cross-section imposes constraints on the on-shell production mode and off-shell production mode of DM candidates. In off-shell production mode (; is the mass of vector boson mediator and is the DM mass), annihilation process can achieve thermal relic DM annihilation cross-section by imposing bound. Whereas for the on-shell production mode (), baryonic fine structure constant would require a slightly larger value. In both production modes, baryonic constant would necessarily be greater than the weak Fermi constant , [33, 34].

6.2. Direct Detection

Direct detection experiments probe the cross-section of DM-nucleon elastic scattering. The CRESST-II [2] experiment provides the best limit on the recoil energy of nucleons. The CRESST-II experiment can explore the sensitivity of DM masses below 0.5 GeV with detection threshold of nuclear recoil 307 eV. The DM-nucleon scattering cross-section for baryonic current can be expressed as where is the reduced mass of DM and nucleons.

6.3. Constraints on Masses of the Vector Boson Mediator

The CDF (Collider Detector at Fermilab) collaboration imposes a rigorous constraints on monojet, jet + missing energy. The limit imposed on the quarks coupling and is largely independent of vector boson mass [35, 36].

6.4.

The Brookhaven alternating gradient synchrotron [37] imposes limit on the branching ratio .

6.5.

Dimension-4 coupling of new light vector to the SM state leads to enhance the production rates; however, the current they couple to is conserved. These rare processes allow us to probe new and better constraints on the coupling of such vectors. These new limits arise from range of processes such as -decays, -decays, rare -decays, and changing meson decays. -decays constraints will be very useful at -factories and LHC experiment. For our work, we have considered K-decays only [3841]. The -decay imposes a limits on the branching ratio of for =1.8 MeV and for  MeV [42, 43].

6.6.

BES (Beijing Spectrometer) collaboration imposes a constraint on the branching ratio for larger value of [44].

6.7. Angular Dependence in Neutron Scattering

The constraints imposed by neutron scattering on the baryonic fine structure constant () which couples the DM candidates and SM particles for mediator vector boson mass  MeV is expressed below [45, 46]:

Existing constraints on the leptophobic DM model are shown in Figures 1 and 2. The plot (4) shows the constraints on the model in the parameter space, while the plot (5) shows the constraints on the model in the parameter space for 10 MeV DM mass and . The other constraints plotted in Figures 1 and 2 are taken from , , Monojet(CDF), Neutron Scattering, and experiments.

7. Simulation and Results

We have used BdNMC (Beam dump Neutrino Monte Carlo) [47] simulation tool to compute the sensitivity of DUNE near detector for leptophobic DM. This simulation tool works on the direct detection technique to probe the sensitivity of the considered experiment running in beam dump mode.

The proton beam of DUNE produces charged and neutral mesons after hitting the beryllium target. These mesons are allowed to propagate through the decay pipe which is made up of steel. The neutral mesons quickly decays in the decay pipe, while charged mesons propagate further and get absorbed in the decay pipe before they decay. The neutral mesons decay into vector mediator which further decays into a pair of DM candidates in the center of mass frame of . The total DM particles produced via distinct production channels (visit Section 3) reach the DUNE near detector which is a cylinder of length 5 m, diameter 5 m, and 1 ton () fiducial mass [19]. Using BdNMC simulation tool, the DM event rate is estimated by imposing production distribution function , specific geometric cuts [19], and recoil energy cuts on nucleons in equation (13). The trajectory of each DM particle which intersects the fiducial mass of near detector is recorded in the form of the energy of DM candidates.

In Figure 3, we have explored the signatures of DM-nucleon scattering events with varying mass of the vector boson mediator for a fixed value of DM mass, i.e., 10 MeV. From this piece of work, we have observed that the DM signatures coming from the DM produced via mesons decay channel, parton-level production channel, and bremsstrahlung are , , and , respectively, in the mass range from 10 MeV to 300 MeV of . For intermediate mass range of , i.e., from 300 MeV to 1 GeV, the DM events arising by DM produced via resonance decay channel are approximately . Whereas the DM signatures coming from the DM produced via parton-level production channel and mesons decay channel are negligible. For GeV, the main contribution of DM events arises from parton-level production channel. In this plot, we notice a sharp peak around  MeV which can be attributed to the resonance production via bremsstrahlung process at this point.

In Figure 1, we have presented the sensitivity of DUNE for leptophobic DM scenario in parameter space where is allowed to vary between 0.03 GeV and 2 GeV, and baryonic fine structure constant has been varied from to . This study is performed for the DM mass 10 MeV and kinetic mixing parameter . The threshold values for recoil nucleons and other essential parameters that have been used in simulation are mentioned in the Table 1. The DM-nucleon elastic scattering event contours for 1, 10, and 1000 events are shown in Figure 1 along with other experimental constraints which are shown by different line colors as mentioned in legends. From Figure 1, we have observed that the allowed parameter space between two sensitivity contours 1 and 10 events is [], the allowed parameter space between two sensitivity contours 10 and 100 events is [], and the allowed parameter space for DM events greater than 1000 events is []. A sharp peak around MeV which can be attributed to the resonance production via bremsstrahlung process at this point. This plot shows that DUNE can probe DM for baryonic fine structure coupling below the existing bound on the coupling value ; these lower coupling values are present for mass of vector boson mediator less than 200 MeV.

Figure 2 shows the contour plot of DM-nucleon elastic scattering events in plane where the variation of is from 0.01 GeV to 1 GeV and is varied from to . This study is performed for kinetic mixing parameter and vector boson mass . The threshold values for recoil nucleons and other essential parameters that are used in simulation for probing DM events are mentioned in Table 1. The three contours of Figure 2 represent the 1, 10, and 1000 events in the parameter space of . From Figure 1, we have observed that the allowed parameter space between two sensitivity contours 1 and 10 events is [], the allowed parameter space between two sensitivity contours 10 and 100 events is [], and the allowed parameter space for DM events greater than 1000 events is []. A sharp peak around  MeV which can be attributed to the resonance production via bremsstrahlung process at this point. In this study, we observe that the vector boson mediator shows weak coupling (below the existing limit of ) for DM masses less than 250 MeV.

Figures 1 and 2 shows that like DUNE, and can also explore the DM signatures at lower coupling values, but DUNE can explore DM signatures for higher values of as compared to the stated experiments. The MiniBooNE sensitivity for leptophobic DM at C.L. is also plotted for comparative study.

8. Conclusion

In our work, we have explored the sensitivity of DUNE in beam dump mode for leptophobic DM. In the considered DM model, DM candidates couple with the SM particles via vector boson mediator of a new baryonic gauge group with coupling strength . The leptophobic DM signatures are observed via DM-nucleon elastic scattering in the parameter space of and . The analysis of Figures 1 and 2 illuminates that the DUNE sensitivity for leptophobic light DM is better as the constraint imposed on the coupling value by DUNE is lower than the present constraint value of the coupling value . These lower coupling values can be explored by DUNE for vector boson mediator mass MeV. Therefore, DUNE in beam dump mode will be able to provide the new results for leptophobic DM. DUNE potential to probe dark matter for lower dark matter coupling values is better than the MiniBooNE potential. These results will help us to understand the nature of DM and its interactions.

Data Availability

All my data are self generated by the dark matter simulation tool “BdNMC”.

Conflicts of Interest

The authors declare that they have no conflicts of interest.