Light in the beam dump
ALP production from decay photons in proton beam-dumps
The exploration of long-lived particles in the MeV-GeV region is a formidable task but it may provide us a unique access to dark sectors. Fixed-target facilities with sufficiently energetic and intense proton beams are an ideal tool for this challenge. In this work we show that the production rate of Axion-Like-Particles (ALPs) coupled pre-dominantly to photons receives a significant contribution from daughter-photons of secondary and mesons created in the proton shower. We carefully compare the PYTHIA simulated spectra of such secondaries to experimental literature, compute the ALP flux from the Primakoff conversion of these photons, and finally revisit existing limits on ALPs and update the prospects for a set of existing and future searches. Our results show that taking this production mechanism into account significantly enhances the sensitivity compared to previous studies based on coherent ALP production in primary proton-nucleus interactions.
At present, the search for new particles is driven by a number of pertinent observations that cannot be explained within the Standard Model (SM). Perhaps most prominent amongst them is Dark Matter (DM), but baryogenesis as well as the strong CP problem also ask for a solution. At the same time, we are challenged by the (so-far) non-detection of new particles. This could be explained by the new particles being extremely weakly coupled. If this is the case they can escape constraints from searches with colliders. But their small coupling also makes such particles long-lived.
In particular, very weakly coupled particles with masses in the MeV-GeV range provide exciting phenomenology and interesting connections to Dark Matter (cf., e.g., Boehm:2003hm (); ArkaniHamed:2008qn (); Freytsis:2010ne (); Dienes:2013xya (); Berlin:2015wwa (); Alekhin:2015byh (); Dolan:2017osp (); Hochberg:2018rjs ()) but also to baryogenesis Akhmedov:1998qx (); Asaka:2005pn (); Shaposhnikov:2008pf () and perhaps even to the strong CP problem Alves:2017avw (); Berezhiani:2000gh (); Agrawal:2017cmd (). Consequently, Heavy Neutral Leptons, Dark Photons and pseudo-scalars of such masses have received increasing attention over the past years (cf., e.g., Beacham:2019nyx () for a recent overview). The latter – pseudo-scalars – also known as as ‘Axion-like particles’ (ALPs) – are the subject of the present document.
ALPs have recently received a considerable amount of interest in the context of dark matter model-building. They may act as a mediator for the interactions between DM and SM particles and thereby allow reproducing the correct Dark Matter relic abundance via thermal freeze-out. At the same time this helps evading the strong constraints from direct and indirect detection experiments Freytsis:2010ne (); Dienes:2013xya (); Berlin:2015wwa (); Dolan:2017osp (). Other motivations for ALPs with masses above include their potential connection to explaining the observed value of the magnetic moment of the muon Marciano:2016yhf (). Also, via the so-called ‘relaxion’ mechanism, ALPs may play a crucial role in electroweak symmetry breaking Graham:2015cka () and in the solution of the hierarchy problem. A concrete implementation Flacke:2016szy () of such a relaxion may yield signatures that are observable with the methods described in this paper. Moreover, it is worth reiterating that recently models have been proposed that allow the QCD axion to live at the MeV-GeV scale Alves:2017avw (); Berezhiani:2000gh (); Agrawal:2017cmd (). Beyond these phenomenological considerations, ALPs are also motivated by top-down extensions of the SM such as string theory Svrcek:2006yi (); Arvanitaki:2009fg (); Acharya:2010zx (); Cicoli:2012sz (). The crucial feature is that a weakly broken or anomalous shift symmetry allows their mass to be much smaller than the fundamental scale, making them accessible to experimental tests. The same shift symmetry then also ensures that the interaction is suppressed by the fundamental scale and therefore very weak. In summary, new searches for ALPs are a well-motivated and timely task.
The main aim of the present paper is to study the production of ALPs in proton beam dumps from decay photons of secondary mesons. We therefore focus on pseudoscalar ALPs whose dominant interaction is with photons. We use the Lagrangian
where is the photon-ALP coupling and is the electromagnetic field strength.
To discover ALPs with masses on the order of MeV-GeV, proton-fixed target facilities (or rather proton beam-dump experiments) are well suited111Indeed there is a sizable number of such beams around the world, cf. e.g. Shiltsev:2013zma (); Shiltsev:2014jpa (); Shiltsev:2017mle (); Alemany:2019vsk () whose suitability for ALP searches has recently been discussed in Harland-Lang:2019zur ().. The strength of such an experimental setup is that it can provide sufficient energy to produce MeV-GeV scale particles, while ensuring that all of the protons in the beam ultimately interact. Moreover, decay volumes spanning tens of meters allow ALPs of various lifetimes to be detected. Overall this combination provides high sensitivity but also excellent complementarity to experiments at low-energy colliders, such as Belle-II Dolan:2017osp (), which can explore the region below a but relatively strong coupling, as well as to experiments at the LHC, which are sensitive mostly to masses above Jaeckel:2012yz (); Mimasu:2014nea (); Jaeckel:2015jla (); Mariotti:2017vtv ().
For ALPs coupled predominantly to photons, there are (at least) two important production mechanisms to be considered. In both cases, one of the photons is provided by a nucleus at rest, constituting the target/dump material. The second photon can be provided
by the charged proton itself. This is sometimes referred to as ‘photon-from-proton’ (PFP) mode, where the photon distribution around the proton is often computed in the Weizsaecker-Williams approximation222More precisely usually variants of the equivalent photon approximation Budnev:1974de () are used Dobrich:2015jyk (); Dolan:2017osp (). In particular in the region of low ALP masses this approximation is problematic and receives sizable corrections that lower the cross section Harland-Lang:2019zur (). This further increases the relative importance of the production mechanism that we will discuss in the present paper.
through a decay, notably from , but also from other neutral mesons.
As already noted, the second process is the main focus of our paper.
Although the importance of the inclusion of axion production from secondary pions, was, for example, already pointed out in Tsai:1986tx (), it is still somewhat under-appreciated. Indeed all the proton beam dump constraints shown in the overview Fig. 1 are calculated using only the PFP production mode Dobrich:2015jyk (); Dolan:2017osp (). To our knowledge, only333 A recent study performed for the FASER experiment Feng:2018pew () at LHC includes the estimate of ALP production directly from decays (which they find subdominant) as well as Primakoff-produced ALPs from decays, albeit in a different kinematic regime. a study put forward in Berlin:2018pwi () which determined prospects for SeaQuest (after its proposed ECAL upgrade) takes into account an estimate of the ALP yield stemming from the Primakoff-conversion of photons from decays in the dump444While also considering Primakoff production from real photons the PrimEx and GlueEx experiments considered in Aloni:2019ruo () are effectively a photon fixed-target experiment where the photons are produced via Bremsstrahlung off an electron beam and then shot onto the target.. However, some simplifying assumptions are made and no full Monte Carlo is set-up.
With our work, we want to close this gap in the literature and give improved estimates for the ALP production from meson-decay photons. The sensitivity improvement with respect to the case when only the PFP production is included will be discussed.
Our study is particularly timely since some of the experiments that can impact the parameter space have started taking data, notably NA62 NA62:2017rwk (), or are close to data-taking.
In practice calculating the photon flux inside the beam dump is far from trivial due to the non-perturbative nature of meson production. We therefore first carefully compare the yields for and other meson and their related angular distributions from PYTHIA Sjostrand:2014zea () simulations to data from past experiments and then use it to determine the photon flux inside the dump.
While here we are interested mostly in ALP production, the and other meson spectra in the dump are not only of relevance to the production of ALPs but can also be the source of Dark Photons and other exotic particles. Thus our work in comparing the yield from PYTHIA simulations to data is of more general interest.
Our paper is structured as follows. A first direct comparison of the PYTHIA simulation output with the experimental data is performed in Section 2: proton-proton and proton-beryllium interactions are separately analyzed. In Section 3 we review and discuss the computation of the ALP yield through production from the meson-decay photons. Finally, in Section 4 we re-evaluate existing experimental constraints and make estimates for future sensitivities taking the additional production mechanism into account. We discuss the conclusions in Section 5.
2 Neutral meson yields in proton beam dumps
The simulation of the production rates of secondary mesons is a challenging task, since the formation of mesons is complicated due to non-perturbative physics. In this section we therefore validate our PYTHIA simulations with experimental data. While no measurements of inclusive neutral meson production are available at exactly the desired energies and target materials employed in the experiments considered in Section 4, we nevertheless have data covering the energy range from 60 GeV to 450 GeV and different target materials, in particular hydrogen and beryllium. Furthermore we can compare to the production of different types of mesons. Putting this together allows us to have at least some confidence in the employed meson spectra. The total cross section for inclusive production agrees with the data within an uncertainty of 20% in the entire beam momentum range of interest. We also indicate kinematic regions where the results are more uncertain. Including and neglecting the contributions from these regions we provide an estimate of the uncertainty of the limits and sensitivities in Section 4.
2.1 400 GeV proton beam on a hydrogen target
Measurements of secondary meson production from a proton beam dumped onto a hydrogen target have been performed at the beginning of the 1990’s by the NA27 experiment operating at the LExan Liquid hydrogen Bubble Chamber (LEBC) with the European Hybrid Spectrometer (EHS). Results for and production from the LEBC-EHS AguilarBenitez:1991yy () allow a direct comparison of the proton-proton interaction expectation from the PYTHIA simulation program Sjostrand:2014zea () with the experimental results. The simulation includes elastic, inelastic non diffractive, and single-, double-diffractive processes. Parton densities for protons are defined using the CTEQ 5L set Lai:1999wy (), a widely-used leading-order QCD parametrization with .
The measurements from LEBC-EHS report that () are produced with an average multiplicity of () per incident proton. A total of of the produced s stem from the decay of secondary particles (mostly mesons). These figures can be compared with the output of the PYTHIA simulation. From it, the total () production multiplicity is (). A total of of the produced s stem from the decay of , , , or , in good agreement with the experimental data. For the total cross section, the SoftQCD set of PYTHIA version 8.2 yields 39.9 mb summing up single- and double-diffractive, non-diffractive, and elastic processes, in good agreement with the experimental data which slightly exceeds 40 mb ppPDG ().
In Figs. 2 and 3, the distributions measured at LEBC-EHS for the squared transverse momentum (), rapidity (), and Feynman variable have been compared to the PYTHIA output, after applying to the Monte Carlo the experimental selection criteria: for (), the condition is , where the Feynman variable is computed as , represents the beam axis direction, and is evaluated in the center of mass frame. To obtain the differential cross sections, the MC output is scaled according to the total proton-proton cross section, 39.14 mb for a proton beam. Four regions of are defined: a central region, for , where the MC overestimates the data by a factor less than two; an intermediate region, for , where data and MC agree; a fragmentation region, for where MC underestimates the data by a factor less than three, and a forward contribution, for , associated to the inelastic diffraction mechanism, where the MC largely underestimates the data. Inserting in the simulation the double-pomeron exchange with the Minimum-Bias Rockefeller model Ciesielski:2012mc () is seen to slightly improve the data-MC comparison for the central region, while not affecting the other regions. Including the double-pomeron exchange increases the total proton-proton cross section by 0.47 mb for a proton beam energy of 400 GeV.
In the following, the domain has been divided into the eight bins defined in Tab. 1. In general, the probability that photons from decays produce ALPs in the experimental acceptance increases with . As discussed in Section 3, we use our knowledge of the quality of our simulated meson spectra by making the following estimate for the uncertainty. As a baseline we take into account the full simulation results including all bins (0-7). This can then be compared to a conservative estimate that only includes bin 5, where the agreement between simulation and data is very good, and bin 6 where the simulation underestimates the data by a moderate factor of up to 3.
2.2 400 GeV proton beam on a beryllium target
We can now take the next step and allow for different target materials. First we note our expectation that 400 GeV protons interacting with a fixed proton target do not differ significantly from those of a neutron target, as far as the or production is concerned. This is confirmed by the PYTHIA Monte Carlo both for the scattering distributions and cross sections. Following from this we obtain the cross sections for larger target nuclei by an appropriate scaling with the geometric cross section , which in this subsection we compare to the data.
Measurements of meson yields from a 400 GeV proton beam on beryllium targets have been performed in Atherton:1980vj (), and data was taken for four values of the secondary particle momenta (60, 120 and 300 GeV) and two values of transverse-momentum (0 and 500 MeV) at different target lengths. To complement these measurements at a lower momentum range of secondary particles and in view of evaluation of neutrino fluxes for NOMAD and CHORUS, the NA56/SPY experiment Ambrosini:1999id () published yields in the range of secondary momentum from 7 to 135 GeV with a proton beam of 450 GeV. In order to make these experimental data useful for further applications, a very useful parametrization was developed in Bonesini:2001iz (), and sometimes is referred to as ‘BMPT’555In the context of the production of exotic particles, BMPT was, e.g., also employed in deNiverville:2016rqh () to predict yields of sub-GeV Dark Matter production in beam dump experiments.. As this parametrization was developed by the extrapolation of 400 GeV and 450 GeV data, it is suited to be employed for the use-case of NA62 and SHiP Alekhin:2015byh (); Anelli:2015pba (), while care has to be taken, when extrapolating to the NuCal Blumlein:1990ay (), and SeaQuest beam energies of 70 GeV and 120 GeV, respectively.
If the contribution of heavier meson and resonance decays is removed, the inclusive cross section is expected to be approximately equal to the average of the and inclusive cross sections. To validate our estimates for based on PYTHIA, we directly compare the MC output to the inclusive invariant cross section obtained in Ambrosini:1999id () for and . For completeness, we also consider the emission of and , protons and anti-protons.
Figure 4 shows the results of our comparison, which is quite good. To arrive at this comparison we have
accounted for the angular acceptance of NA56, corresponding to emission angles below approximately 0.7 mrad;
accounted for the NA56 target-efficiency factor: (where is the target length and is the appropriate proton interaction length)
scaled the MC results for the proton-beryllium total cross section, a factor of higher than that for proton-proton scattering.
Above, we have assumed a target length 100 mm and a proton interaction length 423 mm. In Figure 4, left-hand side, we observe a very satisfying agreement for all available data from NA56 except a slight underestimation of the yield of at large momenta. As NA56 data might include tertiary production up to a certain extent, an underestimation with respect to the simulation output might be expected. On the right-hand-side of Figure 4, for completeness we also show the proton and anti-proton inclusive invariant cross sections measured at NA56 compared to the expected output from PYTHIA.666The agreement observed on the secondary production would allow a reliable estimate of the tertiary production in a dump, which can be expected to be not negligible with respect to the primary production. This evaluation depends on the detailed structure of the dump and is beyond the scope of the present paper.
2.3 Proton beam energies below 400 GeV
To compare the inclusive production of light mesons with available literature, PYTHIA simulations of proton-proton interaction have been produced for proton beam energies of 70, 120, and 250 GeV: the first two values correspond to the beam energy of the NuCal and SeaQuest experiments, while the third value corresponds to the NA22 experiment Adamus:1986ta (), providing the most complete experimental data available below 400 GeV. The total cross section of production and the average number of emitted mesons are shown in Fig. 5 as a function of the beam energy: a general agreement (within 20% relative uncertainty) is observed between data sigmatotLiterature () and MC.
Data and MC differential cross sections have been compared at 250 GeV, as a function of the squared transverse momentum and of the Feynman variable (Fig. 6). Data and MC transverse momentum differential cross sections are seen to agree (even if the data range is limited), while MC underestimates the data for . Again, this will lead to conservative estimates of the expected ALP yields in the following. The correlation between average transverse momentum and observed in data is quite well reproduced by the MC, as shown in Fig. 7.
3 Production of ALPs from meson decay photons
Having validated our PYTHIA spectra for neutral secondaries produced in proton interactions against data in the previous section 2, we are now ready to compute their impact on ALP production. We will also compare it with the contribution of ALPs from the photon-from-proton-mode, which has been evaluated in Ref. Dobrich:2015jyk (). As both processes have different initial and final states there is no interference and both contributions can be be added together. In the following we will therefore concentrate on the new contribution from the mesons.
Taking the meson distributions as an input we set up a Monte-Carlo simulation that proceeds along the following steps:
Decay of the neutral mesons produced in the dump into photons.
Compute the cross-section of ALP production from these photons in the target nucleus.
Mimic appropriate experimental acceptances and cuts.
Evaluate a sensitivity prospect at fixed number of incident protons for the situation of zero background.
Let us start with the first step. The decay length of the neutral pion is given by
The decay length of the other neutral mesons is even smaller. Therefore, effectively the meson decay into photons is instantaneous, i.e. all mesons decay inside the target.
The meson decay therefore yields a distribution of real photons depending on the energy and the angle with respect to the beam axis. Due to the symmetry there is no dependence of this photon distribution with respect to rotations around the beam axis. The distribution is shown in Fig. 8 for proton beam energies of and .
In the second step the photons produce ALPs via the Primakoff conversion on the nucleus. We use laboratory-frame coordinates, where the nucleus is at rest. By a suitable rotation we can choose the photon to be moving in the -direction. We then have,
Note, however, that in this coordinate system the -direction is slightly rotated compared to the beam axis, with the angle arising from the angle of the meson production and subsequent decay into the desired photon. In the Monte-Carlo we then rotate back to the system with the -axis along the beam axis777In practice we do this by using a coordinate system as in Dobrich:2015jyk () that allows for an angle with respect to the beam axis right from the beginning..
We use the approximate form of the cross section as in Dobrich:2015jyk ()888We have checked by comparing to more complete expressions in Cadamuro:2010cz (); Aloni:2019ruo () that the approximation to the cross section is excellent in the regime giving relevant contributions within the experimental acceptance.. In this simple coordinate system the cross section, approximated for small , reads,
As in Dobrich:2015jyk () the electromagnetic form-factor , is taken to be of the Helm form999At very low momentum transfer (depending on the target material) the electron shell also shields the charge, reducing the form factor. However, as discussed in Dobrich:2015jyk () this region only gives a very small contribution to the signal which we neglect here. (cf. Woods:1954zz ()),
with the first spherical Bessel function of the first kind. For the nuclear radius we use Lewin:1995rx ()
To simplify the evaluation we set the form factor to zero for values , i.e. above the first zero of the Helm form factor (as in Dobrich:2015jyk ()).
Let us also emphasize again that in the present case we are dealing with a distribution of real, on-shell photons. This is in contrast to the photon from proton (PFP) mode where we have an effective parton distribution of virtual photons. Therefore, in this production mode we do not expect to be affected by the corrections to the equivalent photon approximation discussed in Harland-Lang:2019zur ().
The cross section is then determined by folding the (probability) distribution of the photons from the mesons with the cross section for an individual photon 101010Let us note at this point that formally the PFP production and the production from decay photons is at different order in . Meson production is a strong process and therefore essentially independent of whereas “radiating” a photon from the proton requires an extra electromagnetic interaction. This is accounted for by a factor of in the relevant photon distribution functions for the proton (cf. also Berlin:2018pwi ()).
As an example we show in Fig. 9 the resulting cross sections for several example energies and target materials relevant for the analysis in the next section (solid lines). This is then compared to the cross section for the PFP mode shown as dashed lines. For all considered energies , and , the PFP mode is sub-dominant. Also, for all energies, production from mesons is particularly favored since the spectrum is also considerably harder. The harder spectrum is useful for detecting ALPs with relatively large couplings since the higher -factor allows ALPs to decay outside dump and shielding regions.
For a better understanding of the behavior shown in Fig. 9 we note that the Primakoff cross section, Eq. (4), is peaked at small angles. More precisely, it is peaked at angles . For sufficiently small masses and high energies this angle is very small. The ALP then has the same angle with respect to the beam axis as the incoming decay photon. Therefore, to a large degree the energy and angular dependence is dominated by that of the initial photon distribution shown in Fig. 8. However, this is only an approximate statement since the cross section (4) also has a significant tail towards larger angles where it is ultimately cut off by the form factor.
To get a better quantitative feeling for the involved scattering angles let us note that the radii, according to Eq. (6), are in the range for the elements we are considering. With for small masses, and the form factor setting in at we find that the tail is typically starting to be cut off at angles . For energies , initial photon angles within a few milliradians of the desired ALP angle contribute. Stated differently, for energies and angles milliradians, the ALP angle is dominated by the angle of the decay photon given in Fig. 8.
Fig. 9 also demonstrates our procedure to estimate the uncertainty in our production rate. For we show two lines for the meson production. One (purple) in which we include all mesons in the full kinematic region and a more conservative estimate (black) where we include only the part of the meson spectrum that is most trustworthy. To be more precise we define two regions in . One covering all the bins, 0-7, in Tab. 1, whereas the other, more conservative one includes only the bins 5-6 where we find better agreement between Monte Carlo and data in our validation procedure. While this is very conservative, it does not fully appreciate the degree to which we underestimate the cross section. Indeed bin 6 gives a significant contribution to the cross section. From Fig. 3 we can see that in this region the simulated cross section underestimates the cross section by a factor which can easily be . Our conservative estimate therefore underestimates the cross section by perhaps . Validated spectra improving in this region could therefore allow for a significantly increased sensitivity.
In the third step we then take into account all the relevant experimental cuts,
The ALP decay has to happen outside the target and in the decay volume.
Crucially this gives the exponential dependence on the decay length, , where is the decay length, the distance from the target to the decay volume and the length of the decay volume.
The decay photons of the ALP have to reach the detector with a suitable minimum energy and other criteria are required for the photons to be detected.
We will describe the relevant details and approximations when we discuss the individual experiments. After all those cuts we obtain the fiducial cross section,
Before we take the final step, let us also note that we have employed the simplification that both the meson production and the photon conversion into ALPs happen at the beginning of the target. This neglects the finite distance traveled before the proton interacts inside the target as well as the additional distance traveled by the photon before it is converted into an ALP. The typical distances are of the order of the proton and photon radiation length, respectively. As both are of the order of cm in the relevant energy range, this should be a minor effect compared to the target sized of the order of m for the experiments we will consider in Section 4. For experiments with targets in the cm range this would have to be taken into account.
In the final step we now need to compare the fiducial cross section for the detectable ALP production from a photon with the cross section for the photon to be absorbed in the target material. We then get for the total number of events,
where is the total number of photons produced in the meson decays. is the total cross section for photons to be absorbed in the target material and can be determined from the radiation length ,
where is the mass of the nucleus.
We stress that this is different from the PFP mode, where the relevant cross section would be the cross section for proton-nucleus interactions.
For the photon cross section we use values for the radiation length from gammaxs (), whereas for the proton cross section we employ from Carvalho:2003pza (), where is the mass number of the nucleus. For example in the case of copper the photon cross section for energies is about for copper, for iron and for molybdenum. This is in contrast to the the proton cross section which is only , and , respectively.
While the photon absorption cross section is significantly larger we find that in many cases this is more than compensated by a number of other factors. In particular we produce on average more than 1 meson per proton and each meson gives two photons. Importantly, as we have seen in Fig. 9, the spectrum of ALPs from these photons is also harder than the one from the PFP mode, which is advantageous for the detection in particular at larger ALP masses.
4 Updated sensitivities for fixed-target experiments
In the following, we will first update the exclusion contours for the past fixed target experiments CHARM Bergsma:1985qz () and NuCal with ALPs produced from the decay of . We then project sensitivities for the existing NA62 NA62:2017rwk () and SeaQuest Berlin:2018pwi () set-ups. Finally, we give projections for the proposed SHiP Alekhin:2015byh (); Anelli:2015pba () facility.
4.1 Past experiments: NuCal and CHARM update
For CHARM Bergsma:1985qz (), we make use of the following parameters: The detector was located at a distance away from the proton dump, and the protons were dumped into a copper target. The detector (for CHARM this is identical to the decay volume) was in length and in transverse dimensions. The detector is off-set transversally by from the beam axis and this is accounted for in the MC. According to Ref. Bergsma:1985qz (), CHARM was sensitive to events with a single electromagnetic shower in acceptance. CHARM quotes a signal acceptance of 51%, which we include in our estimate. The number of protons on target (POT) is .
NuCal Blumlein:1990ay (), made use of the U70 proton beam facility with a beam energy of ‘only’ . However, NuCal profits from a comparably small distance between target and detector of only and a detector length of . We adapt the analysis strategy of Blumlein:2013cua (), and require a minimum ALP energy of and at least one photon detected. In this way, the acceptance is approximately constant and equal to . The detector has a radius of . In a dataset of protons on an iron target, NuCal observed 1 event compared to a background expectation of 0.3 events. At 90% confidence level, we can therefore exclude any point in the parameter space predicting more than 3.6 events.
The changes in the limits from CHARM and Nucal when including ALPs produced from decayed s, can be understood from Figure 9.
For NuCal, where the beam energy is 70 GeV, the plot in the right-hand side of the figure illustrates that, including the yield appreciably changes the existing limits. This can be seen in Figure 10 (l.h.s.), where the brown dashed line shows the NuCal limits in PFP mode only, while the yellow region is the reach considering the added yield for NuCal from PFP and a conservative (bins 5–6) range from decayed mesons. However, the position of the upper part of the exclusion contour (at large couplings) is mostly determined by the experiment’s geometry. Thus, even the revised NuCal limit does not drastically alter the untested parameter space that should be probed by new experiments, discussed in the next section.
Also for CHARM, including the ALP production channel via s improves the sensitivity with respect to the CHARM reach from PFP only considerably111111A recent study Dobrich:2018jyi () on ALPs coupled to fermions is another example of the importance of carefully recasting results from past experiments, particularly CHARM, using PYTHIA and MC simulations to appreciate their full impact.. Figure 10 (l.h.s.) shows a zoomed version of the existing limits with updated CHARM contour (solid red) compared to the previous curve (dashed magenta). However, the CHARM limits lie still within the limits of E137 and the new NuCal limit. Thus we summarize that the experimental future landscape, does not change significantly even after the CHARM and NuCal updates. However, as we will see is the next subsection, the inclusion of the meson contribution, drastically changes the prospect sensitivity for forthcoming searches.
4.2 Current and future set-ups
The NA62 experiment NA62:2017rwk () has been built to achieve a precise measurement of the ultra-rare decay . Besides its main goal, NA62 has a rich program to search for exotic particles, including long-lived particles that can be produced in the up-stream copper beam collimator, into which the primary SPS proton-beam is fully (in dump-mode) or partially (in standard, parasitic data-taking) dumped (see, e.g. Dobrich:2018ezn () for more details).
To be sensitive to a fully neutral final state, NA62 has to be run in beam-dump-mode. In our MC, we model NA62 using the following parameters: The distance between the beam-defining collimator (used to dump the beam) and the start of the fiducial volume is and the vacuum decay region before the Liquid Krypton Calorimeter (LKr) is long. In addition, we require the following acceptance conditions: Both photons produced in the ALP decay need to be detected at a minimum mutual distance of in the LKr. Moreover, these photons need to be away from the LKr central hole and their combined energy needs to be above . The target material for NA62 is copper and we show NA62 prospects for two different choices of POT.
In Figure 10 (l.h.s.) we show the sensitivity for NA62 at POT (corresponding to a one-day run), using (purple dashed) PFP only, using the full range (bins 0–7, blue dotted) and the central bins 5–6 (red), respectively. We also show the sensitivity for POT PFP to facilitate comparison with previous results Dobrich:2015jyk (). As shown, the contribution from the decay of ’s dominates the sensitivity prospect for a proton beam and constitutes one of the most relevant results of our study. Regarding the question, whether all or only a conservative number of bins should be chosen, we observe that the difference between the reach in both situations is small but visible at ALP masses of few tens of MeV and small couplings. This can be understood by looking at Figure 9. The central production corresponding to bin 4 contributes to low photon momenta and correspondingly low ALP momenta. As we can see in Fig. 10 this is more relevant for the sensitivity at low masses and couplings. The forward contribution to production corresponding to bin 7, is not relevant in increasing the cross section at high photon energies. Not including bins other than 5 and 6 has a minor impact on our projections and thus for all following computations we adopt this conservative condition.
In Figure 10 (r.h.s) we show again the prospect for NA62, this time however at POT (corresponding to a few months of data taking) while summing up contribution from production and PFP. Comparing this to the prospects shown, e.g. in Beacham:2019nyx (); Dobrich:2018ezn () underlines the importance of including the yield of ALPs produced by s in these estimates.
The sensitivity reach for SeaQuest to ALPs produced as a result of secondary decays was estimated previously in Berlin:2018pwi (). As outlined in Berlin:2018pwi (), the sensitivity to a di- final state requires the installation of an ECAL, potentially adapted from the PHENIX detector at BNL. Our analysis improves the study put forward in Berlin:2018pwi () in several key regards. Firstly, as outlined in Sect. 2, we have validated the differential cross sections obtained in PYTHIA against experimental literature in a wide energy range. Secondly, as for all considered set-ups, we have implemented a full Monte Carlo of ALP production and decay according to the geometry laid out in Berlin:2018pwi (). Lastly, we consider also the sub-dominant PFP contribution for ALP production in our estimate.
We use the geometric setup described in Berlin:2018pwi (). Also, as in Berlin:2018pwi (), we assume the need for 10 signal events to detect a signal beyond the background fluctuations. The calorimeter is placed between tracking stations 3 and 4 (at downstream of the target) We use a fiducial volume that has a length of , in between meters 7 and 8 of the experiment, and a geometric acceptance of in transverse directions. The target material of SeaQuest is iron and we assume the phase-I statistics of POT. Similarly to other set-ups, we require both photons to be detected at a minimum energy of each, a total energy of at least and a minimum mutual distance of 10 cm to avoid shower overlap given the photon shower Moliere radius.
The resulting prospects for SeaQuest are shown as the brown curve on the r.h.s. of Figure 10. Compared to the estimates of Berlin:2018pwi (), the expected sensitivity covers a somewhat larger area of parameter space.
Finally, we model the prospects for detection of ALPs in the SHiP Alekhin:2015byh (); Anelli:2015pba () calorimeter as follows121212The prospect geometry for SHiP has changed since the publication of Dobrich:2015jyk () and we follow the layout of shipCDS () for our estimate: The fiducial region is taken to be downstream of the production point. The calorimeter is positioned at after the beginning of the decay volume. We ask both photons to be in an acceptance area of 5 10 m. Both photons should have a minimum energy of 1 GeV and a combined energy of 3 GeV and be at least apart. The result can be seen in Fig 10 r.h.s as green curve. Target material is molybdenum and the POT are . Note that the envisaged SHiP calorimeter Bonivento:2018eqn () has the potential of reconstructing the photon direction, thereby allowing ALP mass reconstruction. Compared to the results shown in Beacham:2019nyx () the mass reach increases considerably, from to .
Proton beam dump experiments are a popular and versatile tool to explore the dark sector in the MeV to GeV range that may be connected to a number of open problems in particle physics, notably dark matter. Important examples of current and near future experiments are NA62 (in beam dump mode), SeaQuest, and SHiP.
New particles can be produced in primary interactions of the proton with the target material but also in the decay of secondary mesons. However, further important contributions to the production can arise from the interaction of secondary or even tertiary particles with the target material. While these production mechanisms have been noted and even occasionally used Berlin:2018pwi (); Feng:2018pew () they are still somewhat under-appreciated and many sensitivity calculations do not take them into account.
In this paper we have performed a detailed investigation of axion-like particle (ALP) production. More precisely we discussed the production of ALPs from the following process: Protons interact with the target nucleus and produce neutral mesons. The mesons (mostly ) decay into two which subsequently can interact with another target nucleus to produce an ALP via the Primakoff process. We show that this gives a significant contribution to the production of ALPs which is also kinematically well suited for detection in typical experimental setups. Indeed, for experiments with high beam energies such as NA62 or SHiP, this is the dominant contribution in the region of interest and significantly extends the mass reach, e.g. by a factor of in the case of the SHiP experiment.
A crucial input for the calculation of the production with secondary or even higher order particles are the spectra of these particles inside the target. In particular for mesons theoretical predictions are challenging due to the non-perturbative nature of the meson production processes. We have therefore validated our simulation results from PYTHIA 8.2 against a variety of measurements, thereby giving an estimate of the reliability of the simulation results and the impact this has on the sensitivity calculation for the experiments. We find that the impact of the uncertainty is moderate despite a relatively large uncertainty of our generated meson spectra in some regions of phase space. Further, very desirable improvement could come from two directions. First of all, our simulations only include mesons produced from the primary proton beam but, also for meson production, secondary interactions may play a sizable role. Including these secondaries will be an important next step. Second, the discrepancy of data and MC should be clarified by an extended study of existing data or, if needed, new measurements.
All in all, interactions of secondary particles in the beam dump are a powerful additional production mode for new very weakly coupled particles. Further studies that go beyond the example presented in this work are needed and in preparation.
We would like to thank M. Pospelov for giving the nudge that started this work as well as P. DeNiverville and F. Kahlhoefer for very useful discussions on the presented physics. We thank R. Wanke for clarifications on the SHiP prospectus setup. TS would like to thank P. Di Nezza for useful discussions on the validation of the MC simulation. BD acknowledges support through the European Research Council (ERC) under grant ERC-2018-StG-802836 (AxScale).
- (1) C. Boehm and P. Fayet, “Scalar dark matter candidates,” Nucl. Phys. B 683 (2004) 219 [hep-ph/0305261].
- (2) N. Arkani-Hamed, D. P. Finkbeiner, T. R. Slatyer and N. Weiner, “A Theory of Dark Matter,” Phys. Rev. D 79 (2009) 015014 [arXiv:0810.0713 [hep-ph]].
- (3) M. Freytsis and Z. Ligeti, “On dark matter models with uniquely spin-dependent detection possibilities,” Phys. Rev. D 83 (2011) 115009 [arXiv:1012.5317 [hep-ph]].
- (4) K. R. Dienes, J. Kumar, B. Thomas and D. Yaylali, “Overcoming Velocity Suppression in Dark-Matter Direct-Detection Experiments,” Phys. Rev. D 90 (2014) no.1, 015012 [arXiv:1312.7772 [hep-ph]].
- (5) A. Berlin, S. Gori, T. Lin and L. T. Wang, “Pseudoscalar Portal Dark Matter,” Phys. Rev. D 92 (2015) 015005 [arXiv:1502.06000 [hep-ph]].
- (6) S. Alekhin et al., “A facility to Search for Hidden Particles at the CERN SPS: the SHiP physics case,” Rept. Prog. Phys. 79 (2016) no.12, 124201 [arXiv:1504.04855 [hep-ph]].
- (7) M. J. Dolan, T. Ferber, C. Hearty, F. Kahlhoefer and K. Schmidt-Hoberg, “Revised constraints and Belle II sensitivity for visible and invisible axion-like particles,” JHEP 1712 (2017) 094 [arXiv:1709.00009 [hep-ph]].
- (8) Y. Hochberg, E. Kuflik, R. Mcgehee, H. Murayama and K. Schutz, “Strongly interacting massive particles through the axion portal,” Phys. Rev. D 98 (2018) no.11, 115031 [arXiv:1806.10139 [hep-ph]].
- (9) E. K. Akhmedov, V. A. Rubakov and A. Y. Smirnov, “Baryogenesis via neutrino oscillations,” Phys. Rev. Lett. 81 (1998) 1359 [hep-ph/9803255].
- (10) T. Asaka and M. Shaposhnikov, “The nuMSM, dark matter and baryon asymmetry of the universe,” Phys. Lett. B 620 (2005) 17 [hep-ph/0505013].
- (11) M. Shaposhnikov, “The nuMSM, leptonic asymmetries, and properties of singlet fermions,” JHEP 0808 (2008) 008 [arXiv:0804.4542 [hep-ph]].
- (12) D. S. M. Alves and N. Weiner, “A viable QCD axion in the MeV mass range,” JHEP 1807 (2018) 092 [arXiv:1710.03764 [hep-ph]].
- (13) Z. Berezhiani, L. Gianfagna and M. Giannotti, “Strong CP problem and mirror world: The Weinberg-Wilczek axion revisited,” Phys. Lett. B 500 (2001) 286 [hep-ph/0009290].
- (14) P. Agrawal, J. Fan, M. Reece and L. T. Wang, “Experimental Targets for Photon Couplings of the QCD Axion,” JHEP 1802 (2018) 006 [arXiv:1709.06085 [hep-ph]].
- (15) J. Beacham et al., “Physics Beyond Colliders at CERN: Beyond the Standard Model Working Group Report,” arXiv:1901.09966 [hep-ex].
- (16) W. J. Marciano, A. Masiero, P. Paradisi and M. Passera, “Contributions of axionlike particles to lepton dipole moments,” Phys. Rev. D 94 (2016) no.11, 115033 [arXiv:1607.01022 [hep-ph]].
- (17) P. W. Graham, D. E. Kaplan and S. Rajendran, “Cosmological Relaxation of the Electroweak Scale,” Phys. Rev. Lett. 115 (2015) no.22, 221801 [arXiv:1504.07551 [hep-ph]].
- (18) T. Flacke, C. Frugiuele, E. Fuchs, R. S. Gupta and G. Perez, “Phenomenology of relaxion-Higgs mixing,” JHEP 1706 (2017) 050 [arXiv:1610.02025 [hep-ph]].
- (19) P. Svrcek and E. Witten, “Axions In String Theory,” JHEP 0606 (2006) 051 [hep-th/0605206].
- (20) A. Arvanitaki, S. Dimopoulos, S. Dubovsky, N. Kaloper and J. March-Russell, “String Axiverse,” Phys. Rev. D 81 (2010) 123530 [arXiv:0905.4720 [hep-th]].
- (21) B. S. Acharya, K. Bobkov and P. Kumar, “An M Theory Solution to the Strong CP Problem and Constraints on the Axiverse,” JHEP 1011 (2010) 105 [arXiv:1004.5138 [hep-th]].
- (22) M. Cicoli, M. Goodsell and A. Ringwald, “The type IIB string axiverse and its low-energy phenomenology,” JHEP 1210 (2012) 146 [arXiv:1206.0819 [hep-th]].
- (23) V. Shiltsev, S. Henderson, P. Hurh, I. Kourbanis and V. Lebedev, “Issues and RD Required for the Intensity Frontier Accelerators,” arXiv:1305.6917 [physics.acc-ph].
- (24) V. Shiltsev, S. Henderson, P. Hurh, I. Kourbanis and V. Lebedev, “Issues and RD Required for the Intensity Frontier Accelerators,” arXiv:1409.5426 [physics.acc-ph].
- (25) V. Shiltsev, “Fermilab Proton Accelerator Complex Status and Improvement Plans,” Mod. Phys. Lett. A 32 (2017) no.16, 1730012 [arXiv:1705.03075 [physics.acc-ph]].
- (26) R. Alemany et al., “Summary Report of Physics Beyond Colliders at CERN,” arXiv:1902.00260 [hep-ex].
- (27) L. Harland-Lang, J. Jaeckel and M. Spannowsky, “A fresh look at ALP searches in fixed target experiments,” arXiv:1902.04878 [hep-ph].
- (28) J. Jaeckel, M. Jankowiak and M. Spannowsky, “LHC probes the hidden sector,” Phys. Dark Univ. 2 (2013) 111 [arXiv:1212.3620 [hep-ph]].
- (29) K. Mimasu and V. Sanz, “ALPs at Colliders,” JHEP 1506 (2015) 173 [arXiv:1409.4792 [hep-ph]].
- (30) J. Jaeckel and M. Spannowsky, “Probing MeV to 90 GeV axion-like particles with LEP and LHC,” Phys. Lett. B 753 (2016) 482 [arXiv:1509.00476 [hep-ph]].
- (31) A. Mariotti, D. Redigolo, F. Sala and K. Tobioka, “New LHC bound on low-mass diphoton resonances,” Phys. Lett. B 783 (2018) 13 [arXiv:1710.01743 [hep-ph]].
- (32) V. M. Budnev, I. F. Ginzburg, G. V. Meledin and V. G. Serbo, “The Two photon particle production mechanism. Physical problems. Applications. Equivalent photon approximation,” Phys. Rept. 15 (1975) 181.
- (33) B. Döbrich, J. Jaeckel, F. Kahlhoefer, A. Ringwald and K. Schmidt-Hoberg, “ALPtraum: ALP production in proton beam dump experiments,” JHEP 1602 (2016) 018 [JHEP 1602 (2016) 018] [arXiv:1512.03069 [hep-ph]].
- (34) Y. S. Tsai, “Axion Bremsstrahlung By An Electron Beam,” Phys. Rev. D 34 (1986) 1326.
- (35) J. L. Feng, I. Galon, F. Kling and S. Trojanowski, “Axionlike particles at FASER: The LHC as a photon beam dump,” Phys. Rev. D 98 (2018) no.5, 055021 [arXiv:1806.02348 [hep-ph]].
- (36) D. Aloni, C. Fanelli, Y. Soreq and M. Williams, “Photoproduction of axion-like particles,” arXiv:1903.03586 [hep-ph].
- (37) A. Berlin, S. Gori, P. Schuster and N. Toro, “Dark Sectors at the Fermilab SeaQuest Experiment,” Phys. Rev. D 98 (2018) no.3, 035011 [arXiv:1804.00661 [hep-ph]].
- (38) E. Cortina Gil et al. [NA62 Collaboration], “The Beam and detector of the NA62 experiment at CERN,” JINST 12 (2017) no.05, P05025 [arXiv:1703.08501 [physics.ins-det]].
- (39) T. Sjöstrand et al., “An Introduction to PYTHIA 8.2,” Comput. Phys. Commun. 191 (2015) 159 [arXiv:1410.3012 [hep-ph]].
- (40) H. L. Lai et al. [CTEQ Collaboration], “Global QCD analysis of parton structure of the nucleon: CTEQ5 parton distributions,” Eur. Phys. J. C 12 (2000) 375 [hep-ph/9903282].
- (41) M. Aguilar-Benitez et al., “Inclusive particle production in 400-GeV/c p p interactions,” Z. Phys. C 50 (1991) 405.
- (42) C. Patrignani et al. [Particle Data Group], “Review of Particle Physics,” Chin. Phys. C 40 (2016) no.10, 100001.
- (43) R. Ciesielski and K. Goulianos, “MBR Monte Carlo Simulation in PYTHIA8,” PoS ICHEP 2012 (2013) 301 [arXiv:1205.1446 [hep-ph]].
- (44) H. W. Atherton et al., “Precise Measurements of Particle Production by 400-GeV/ Protons on Beryllium Targets,” CERN-80-07, CERN-YELLOW-80-07.
- (45) G. Ambrosini et al. [NA56/SPY Collaboration], “Measurement of charged particle production from 450-GeV/c protons on beryllium,” Eur. Phys. J. C 10 (1999) 605.
- (46) M. Bonesini, A. Marchionni, F. Pietropaolo and T. Tabarelli de Fatis, “On Particle production for high-energy neutrino beams,” Eur. Phys. J. C 20 (2001) 13 [hep-ph/0101163].
- (47) P. deNiverville, C. Y. Chen, M. Pospelov and A. Ritz, “Light dark matter in neutrino beams: production modelling and scattering signatures at MiniBooNE, T2K and SHiP,” Phys. Rev. D 95 (2017) no.3, 035006 [arXiv:1609.01770 [hep-ph]].
- (48) I. V. Aiinenko et al. [NA22 Collaboration], “Inclusive Production in , and Interactions at 250-GeV/,” Z. Phys. C 35 (1987) 7 [Yad. Fiz. 47 (1988) 429] [Sov. J. Nucl. Phys. 47 (1988) 271].
- (49) K. Jaeger et al., “Inclusive and Production in 12.4-GeV/c p p Interactions,” Phys. Rev. D 11 (1975) 1756 and references 11–17 therein.
- (50) D. Cadamuro, S. Hannestad, G. Raffelt and J. Redondo, “Cosmological bounds on sub-MeV mass axions,” JCAP 1102 (2011) 003 [arXiv:1011.3694 [hep-ph]].
- (51) R. D. Woods and D. S. Saxon, “Diffuse Surface Optical Model for Nucleon-Nuclei Scattering,” Phys. Rev. 95 (1954) 577.
- (52) J. D. Lewin and P. F. Smith, “Review of mathematics, numerical factors, and corrections for dark matter experiments based on elastic nuclear recoil,” Astropart. Phys. 6 (1996) 87.
- (53) https://www.nist.gov/pml/xcom-photon-cross-sections-database
- (54) J. Carvalho, “Compilation of cross sections for proton nucleus interactions at the HERA energy,” Nucl. Phys. A 725 (2003) 269.
- (55) F. Bergsma et al. [CHARM Collaboration], “Search for Axion Like Particle Production in 400-GeV Proton - Copper Interactions,” Phys. Lett. 157B (1985) 458.
- (56) M. Anelli et al. [SHiP Collaboration], “A facility to Search for Hidden Particles (SHiP) at the CERN SPS,” arXiv:1504.04956 [physics.ins-det].
- (57) J. Blumlein et al., “Limits on neutral light scalar and pseudoscalar particles in a proton beam dump experiment,” Z. Phys. C 51 (1991) 341.
- (58) J. Blümlein and J. Brunner, “New Exclusion Limits on Dark Gauge Forces from Proton Bremsstrahlung in Beam-Dump Data,” Phys. Lett. B 731 (2014) 320 [arXiv:1311.3870 [hep-ph]].
- (59) B. Döbrich, F. Ertas, F. Kahlhoefer and T. Spadaro, “Model-independent bounds on light pseudoscalars from rare B-meson decays,” Phys. Lett. B 790 (2019) 537 [arXiv:1810.11336 [hep-ph]].
- (60) B. Döbrich [NA62 Collaboration], “Dark Sectors at fixed targets: The example of NA62,” [arXiv:1807.10170 [hep-ex]] or https://indico.cern.ch/event/523655/contributions/2246416/
- (61) https://cds.cern.ch/record/2654870?ln=en
- (62) W. M. Bonivento [SHiP Collaboration], “Studies for the electro-magnetic calorimeter SplitCal for the SHiP experiment at CERN with shower direction reconstruction capability,” JINST 13 (2018) no.02, C02041.