ALPtraum: ALP production in proton beam dump experiments
Abstract
With their high beam energy and intensity, existing and nearfuture proton beam dumps provide an excellent opportunity to search for new very weakly coupled particles in the MeV to GeV mass range. One particularly interesting example is a socalled axionlike particle (ALP), i.e. a pseudoscalar coupled to two photons. The challenge in proton beam dumps is to reliably calculate the production of the new particles from the interactions of two composite objects, the proton and the target atoms. In this work we argue that Primakoff production of ALPs proceeds in a momentum range where production rates and angular distributions can be determined to sufficient precision using simple electromagnetic form factors. Reanalysing past proton beam dump experiments for this production channel, we derive novel constraints on the parameter space for ALPs. We show that the NA62 experiment at CERN could probe unexplored parameter space by running in ‘dump mode’ for a few days and discuss opportunities for future experiments such as SHiP.
Keywords:
Mostly Weak Interactions: Beyond Standard Model; Non collider experiments with beams: Fixed target experiments; ExoticsCERNPHTH2015293, DESY 15237 CERN, 1211 Geneva 23, Switzerland Institut für Theoretische Physik, Universität Heidelberg, Philosophenweg 16, 69120 Heidelberg, Germany DESY, Notkestrasse 85, 22607 Hamburg, Germany DESY, Notkestrasse 85, 22607 Hamburg, Germany DESY, Notkestrasse 85, 22607 Hamburg, Germany
1 Introduction
Over the last few years Nature has been quite reluctant to reveal its secrets beyond the Standard Model — certainly not for a lack of trying on our side. While the second run of the LHC still promises significant discovery potential, it also seems sensible to carefully examine our search strategy and look for potential places we have missed as well as for new opportunities that can be explored in parallel. In this endeavour it is only prudent that we first look at technologies and apparatuses that are already at our disposal and then turn to more ambitious future projects.
In this note we want to discuss an explicit and particularly promising example of this strategy: we will show that without major modifications an existing proton fixed target experiment, NA62 Hahn:1404985 () at the CERN SPS, offers discovery potential for axionlike particles (ALPs) in the MeV to GeV range and can be competitive with existing experiments in a rather short runtime. A proposed dedicated search at the SPS, SHiP Anelli:2015pba (); Alekhin:2015byh (), could then probe even deeper into untested parameter space.
Fixed target experiments are particularly useful to search for new weakly coupled particles in the MeV to GeV range Bjorken:1988as (); Blumlein:1990ay (); Bjorken:2009mm (); Andreas:2012mt (), because they nicely combine a sufficiently energetic reaction to produce particles in this mass regime with a sufficiently high reaction rate to probe small couplings. The two most common types of beams are electron and proton beams. Indeed some of the best current limits on axionlike particles in the MeV to GeV mass range originate from electron beam fixedtarget experiments Konaka:1986cb (); Riordan:1987aw (); Bjorken:1988as (); Davier:1989wz () as well as electronpositron colliders Jaeckel:2015jla (). As demonstrated in Blumlein:1990ay (); Bergsma:1985qz (); Batell:2009di (), another promising option are proton fixedtarget experiments, in particular those making use of the high intensity of 400 GeV protons from the SPS at CERN. It is therefore a worthwhile task to determine the sensitivity of existing experiments, such as the NA62 experiment, as well as proposed experiments especially optimised for the search of longlived neutral particles, such as SHiP.
To fully exploit these experiments for the search of axionlike particles one needs a reliable calculation of the production rates and angular distributions. Due to the composite nature of both the proton and the nucleus this seems particularly challenging. However, we argue that for a highenergy proton beam one can reliably calculate the cross section for the fusion of two photons into one ALP (socalled Primakoff production Halprin:1966zz ()) for ALPs in the interesting MeV to GeV mass range using simple atomic form factors. The reason for this is as follows: Both the proton and the nucleus are surrounded by the virtual photons that make up the usual electric field of a charged particle. Due to the highly relativistic nature of the incoming protons, the photons from the proton now ‘see’ the photons from the nucleus as highly blue shifted, thereby being able to provide more energy/momentum. This effect allows to create relatively massive ALPs from photons that, in the rest frame of the proton/nucleus, are too soft to be affected by the substructure of the proton/nucleus. In other words, the momentum transfer in the respective rest frames is sufficiently small that we can use simple electromagnetic form factors to describe the photon interactions.
At the same time, for the masses of interest to us, the typical photons from the nucleus have sufficient momentum to be essentially unaffected by the shielding from the electron shell, leading to an enhancement of the cross section proportional to the square of the nucleus charge. Coherent Primakoff production in this Goldilocks zone can therefore dominate over incoherent production processes and lead to detectable signals in present and nearfuture experiments for a significant range of ALP parameter space that is not yet explored. Crucially, it is possible for us to reliably calculate not only the production rate, but also angular distributions, which are particularly important to obtain a realistic calculation of the geometric acceptances that determine the sensitivity of a real experiment.
This paper is structured as follows. In section 2 we review the nature of and motivation for ALPs and for searches in the lowenergy/highintensity regime. In section 3 we then provide a detailed calculation of the Primakoff production of ALPs from high energy protons impinging on a fixed target, paying particular attention to the angular distributions. In section 4 we proceed to discuss the general features of the production and decay relevant for the calculation of the experimental acceptances before we use them in section 5 in order to derive new constraints on the ALP parameter space from CHARM Bergsma:1985qz () and NuCal Blumlein:1990ay (). In section 6 we then estimate the sensitivity of a potential run at NA62, followed by a similar analysis for the proposed SHiP experiment in section 7. We conclude in section 8.
2 A brief review and motivation of ALPs
Progress in particle physics has been guided by the paradigm of renormalizable interactions with dimensionless couplings. This paradigm suggests that any new particle to be discovered should be heavy and reveal its presence either directly in highenergy colliders or indirectly by mediating higherdimensional interactions, which could then be probed in precision measurements such as muon experiments or flavour experiments. It has however become increasingly clear that there are other options. Even light particles could still remain to be discovered, provided they have sufficiently small interactions with Standard Model (SM) particles, and therefore with our experiments.
Popular examples for such PseudoGoldstone bosons are axionlike particles, which are loosely defined as (pseudo)scalar particles coupled to the SM particles by dimension5 couplings to two gauge bosons or derivative interactions to fermions (the socalled axion portal Nomura:2008ru (); Batell:2009di ()). Light pseudoscalars have received a considerable amount of interest recently in the context of dark matter modelbuilding, because they may act as a mediator for the interactions between dark matter and SM particles. In these models it is easily possible to reproduce the observed dark matter relic abundance from thermal freezeout while evading the strong constraints from direct and indirect detection experiments Boehm:2014hva (); Berlin:2015wwa ().
In the present work we focus on pseudoscalar ALPs whose dominant interaction is with photons^{1}^{1}1Other possibilities are discussed in section 8. and we are interested in masses and energies of the order of MeV to GeV, significantly below the scale of electroweak symmetry breaking. The relevant Lagrangian is then
(1) 
where denotes the photon ALP coupling.
The origin of such an effective coupling to photons can be motivated in analogy to the case of the axion and is generic for any (pseudo)Goldstone boson of an axial symmetry with nonvanishing electromagnetic anomaly. One naturally expects the effective coupling to be of order
(2) 
Probing a small value of therefore effectively allows us to probe a large scale of the underlying more fundamental theory from which the pseudoGoldstone boson arises. We emphasise that we take the ALP mass and the ALPphoton couplings to be independent parameters.
In terms of the effective ALPphoton coupling , the ALP decay width is^{2}^{2}2We assume that the total decay width is dominated by the ALPphoton coupling and that couplings to leptons and quarks — if present — give subleading contributions. We discuss the impact of additional decay channels in section 8.
(3) 
and the ALP lifetime is given by . For an ALP with energy in the laboratory frame, the typical decay length is then given by
(4) 
A given experiment will be most sensitive to ALPs with a decay length comparable to the distance between target and the detector. Particles with shorter decay length are likely to decay before they reach the decay volume and the decay products will be absorbed. Crucially, larger couplings imply shorter decay lengths and therefore lead to an exponential suppression of the expected number of events in a given experiment. It is therefore a great challenge to probe ALPphoton couplings in the range for ALP masses above (cf. figure 1).^{3}^{3}3Both smaller couplings and smaller ALP masses are in fact very strongly constrained by astrophysical and cosmological observations. Larger couplings, on the other hand, can be tested directly at colliders such as LEP or the LHC Jaeckel:2015jla (). While these couplings are large enough to produce a significant number of ALPs in the target of a beam dump experiment, the fraction of ALPs that reach the detector depends sensitively on the detector geometry and the beam energy. The higher the beam energy and the shorter the distance between target and detector, the larger ALPphoton couplings can be probed. The high beam energy of proton beam dump experiments is therefore suited for making progress in the large coupling window. This effect can be seen in figure 1 (cf. section 5 for details). The turquoise region from the proton beam dump experiment CHARM extends beyond the limit from the electron beam dump experiment SLAC 137, even though the former has a longer distance to the decay volume. Nevertheless, this can only partially compensate the geometric limitations. For better sensitivity an improved geometry as present in NA62 and even more so at SHiP is essential.
For small couplings and masses the decay length is longer than the distance from the ALP production point to the most downstream detectors in the experiment. Hence only a fraction of ALPs decay such that the resulting photons can be detected in the detector (cf. eq. (2)). This leads to the typical diagonal ellipse for the probed regions as shown in figure 1.
3 ALP production in proton beam dumps
Proton beam dumps can produce ALPs via coherent scattering (i.e. elastic scattering of the entire proton on the nucleus in the target), incoherent scattering (i.e. scattering of individual quarks or gluons from the proton and/or the nucleus) and nonperturbative processes (e.g. in the decay of hadronic resonances). The best example for coherent scattering is the Primakoff process shown in figure 2. As we will see, in contrast to incoherent scattering and nonperturbative processes, the Primakoff process remains perturbative at low energies and is not affected by hadronic uncertainties, so that predictions can be calculated with good accuracy. Moreover, the Primakoff process is particularly interesting for beam dump experiments for the following reasons:

The cross section for coherent scattering is typically proportional to the square of the charge of the target nucleus, , while incoherent scattering and nonperturbative processes are roughly proportional to the mass number of the target nucleus, . For heavy target materials, such as lead or tungsten, this can lead to a very significant enhancement of the rate of ALP production. But even for mediumweight target materials like copper, iron or molybdenum Primakoff production is still competitive with other production mechanisms Blumlein:1990ay ().

The typical momentum transfer in coherent scattering is small, such that the ALPs produced via the Primakoff process typically have very small transverse momenta (i.e. momenta perpendicular to the direction of the beam). As a result, cross sections are very strongly peaked in the forward direction, such that even a relatively small detector far away from the target can have a large geometric acceptance. This makes it possible to use the Primakoff process to search for ALPs with very large decay length.
Primakoff production is therefore not only the theoretically cleanest production mode, but can be expected to give the dominant contribution for ALPs that do not directly couple to quarks.^{4}^{4}4For ALPs with direct quark couplings stronger constraints are expected to arise from flavourchanging rare decays, such as Hiller:2004ii (); Freytsis:2009ct (); Andreas:2010ms (); Dolan:2014ska (). However, it is very difficult to accurately predict the rate and distribution of ALPs produced in this way. In the present work we will therefore focus on this production mechanism.
The Primakoff process can be studied in two different ways. In the laboratory frame, where the target nucleus is at rest, it makes sense to calculate the equivalent photon spectrum for the proton beam and then consider the process , i.e. calculate the probability for the photon to emit an ALP before being absorbed by the nucleus. In the centreofmass frame, on the other hand, both the proton and the target nucleus are moving, so one can calculate the equivalent photon spectrum for each and then consider the photonfusion process . The former approach was considered previously in Bjorken:1988as () for the analysis of electron beam dumps, the latter approach was employed by Blumlein:1990ay () for a proton beam dump experiment.
Of course, both approaches are physically equivalent and should lead to identical predictions for the ALP production cross section and distribution. However, as we will discuss in detail below, in order to accurately predict the angular distribution of the ALPs produced in the beam dump, it is essential to properly take into account the transverse momentum of the virtual photons (see e.g. Blumlein:2011mv (); Blumlein:2013cua ()). These transverse momenta have been neglected entirely in Blumlein:1990ay (), whereas in Bjorken:1988as () only the virtual photon connected to the target nucleus is allowed to have nonzero transverse momentum.^{5}^{5}5As we will see below, this is a valid approximation for electron beam dump experiments, because the transverse momentum of the virtual photon cannot significantly exceed the mass of the radiating particle. However, the approximation is not sufficient when considering the case of proton beam dump experiments. We therefore revisit both approaches, carefully taking into account all transverse momenta, and derive the angular distribution for ALPs produced in proton beam dumps.
3.1 The equivalent photon approximation
The equivalent photon approximation (EPA), also known as the WeizsäckerWilliams approximation, provides a convenient framework for studying processes involving photons emitted from fastmoving charges. The basic idea is to replace the charged particle(s) in the initial state by photons following a distribution and then consider directly the interactions of these photons. Since the transverse momentum of the photons is typically very small, it is common in the literature to integrate the distribution over and consider only the energy distribution , assuming that all photons travel in the same direction as the charged particle they originated from.^{6}^{6}6This is for example the case for the photonfromproton mode implemented in MadGraph 5, version 2.3.3 Alwall:2011uj (). For our purposes, however, it will be crucial to accurately predict angular distribution down to very small scattering angles, so that we need to include the distribution of transverse momenta. Our discussion follows ref. Budnev:1974de () (see ref. Martin:2014nqa () for a recent review).
Let us consider a particle of mass and energy and define the energy fraction of the photon as . The magnitude of the momentum transfer between particle and photon is given by , where
(5) 
Defining , the distribution of photons can then be written as
(6) 
where is the finestructure constant and we have dropped a term proportional to in the second line, assuming . We note in particular that transverse momenta are predicted to be of order . This observation makes clear why it is a good approximation to neglect transverse momenta for the equivalent photon spectra of electron beams.
The equation above is valid for a pointlike particle with unit charge. For protons or ions it needs to be modified to include the appropriate electromagnetic form factors:
(7) 
For protons, we follow Budnev:1974de () and take
(8) 
where
(9) 
with the magnetic moment of the proton and . To make a conservative estimate we set the both form factors to zero for .
For heavy ions with charge number and mass number , one finds
(10) 
where and are the charge and magnetic moment form factors, respectively. Since the contribution from the magnetic moment is not proportional to , it can be neglected to good approximation. For the charge form factor we use the Helm form factor, i.e.
(11) 
where is the first spherical Bessel function of the first kind and we employ a parametrization as provided in Lewin:1995rx (), with
(12) 
where . The Helm form factor vanishes for and we set it to zero for .
For scattering on neutral atoms one furthermore needs to take into account screening of the nuclear charge by electrons. These screening effects will become important for small , implying that one has to consider the electromagnetic form factor of the entire atom rather than the one of the nucleus. In principle, it is straightforward to extract the atomic form factor at smaller momentum transfer from measurements (see e.g. Henke:1993eda ()). At energies above the Kshell threshold we do not expect large deviations from , but even for smaller energies the deviations may be relatively small. Comparing to experimentally measured form factors in the literature nist (), we conclude that for molybdenum it is a good approximation to take for keV, while for copper and iron holds for keV. For smaller momenta we set the form factor equal to zero, which has a negligible effect on the predicted cross section.
3.2 Method 1: Photon fusion
Let us first consider the production of ALPs in the centreofmass frame. If the proton has momentum in the laboratory frame and we denote the mass of the target nucleus with , the centreofmass energy is approximately given by
(13) 
The momentum of the two particles in the centreofmass frame is therefore
(14) 
and the corresponding energies are given by . Denoting the respective equivalent photon spectra by and , where and , the total cross section for ALP production is then given by
(15) 
where the cross section for the process is calculated in appendix A to be
(16) 
The invariant mass of the twophoton system is given by
(17) 
Moreover, the momentum of the ALP along the direction of the beam is given by . Changing variables from and to and and performing the integration over we obtain
(18) 
It is straightforward from this expression to extract the momentum distribution of ALPs along the direction of the beam, . However, as discussed above, we are also interested in the momentum distribution perpendicular to the direction of the beam, . Calculating this distribution is complicated by the fact that the transverse momenta of the two photons do not have to be parallel and therefore need to be added vectorially . Nevertheless, assuming that the angle between the two photons is uniformly distributed between 0 and , it is possible to analytically average over the photonphoton angle. As shown in appendix B, we obtain
(19) 
where the area of integration is defined by . The integral can easily be evaluated numerically. The resulting distribution can then be transformed into the laboratory frame, in order to determine . We show two examples for the resulting distributions for different ALP masses in figure 3 at .
3.3 Method 2: Photon absorption
An alternative way to calculate the rate and distribution of ALPs produced in a proton beam dump experiment is to consider the equivalent photon spectrum of the proton beam in the laboratory frame and calculate the probability that a photon emits an ALP before being absorbed by the nucleus. In other words, we need to calculate the cross section for the process and multiply with the distribution of photon momenta , where now . As shown in appendix A, the relevant cross section is approximately given by
(20) 
where is the form factor for the target nucleus (see eq. (11)) and
(21) 
with the transverse photon momentum, the angle between the ALP momentum and the beam direction and the angle between the transverse momenta of ALP and photon.
The cross section for the process is then given by
(22) 
where we have made use of the fact that and taken . Compared to eq. (19), the advantage of eq. (22) is that it only requires integration over one transverse momentum. However, in the latter case it is not possible to analytically perform the averaging over the angle of the photon in the transverse direction, so that both methods require a comparable computational effort. Moreover, as discussed in detail in appendix A, eq. (22) is valid only for small angles and small , whereas no such approximation was necessary for the derivation of eq. (19). Nevertheless, as long as these approximations are valid, the two methods should be in agreement. We confirm in figure 4 that this is indeed the case.
To conclude this section, let us comment on the ALP angular distribution. First of all, we emphasise that in figure 3 we show , in order to reflect the larger solid angle at larger values of . We find that (for a beam energy of ) this distribution peaks at around , implying that an ideal detector should aim to cover this angle in order to capture the majority of the ALPs produced in the target. For comparison, we show in figure 4, which remains finite for . However, we observe that — at least for small ALP masses — even this distribution peaks at nonzero . The reason is that, as can be seen from eq. (3.1), a photon with energy has a typical transverse momentum of order and hence a typical angle of . For a beam energy of the maximum is therefore found around .^{7}^{7}7To produce an ALP with mass and energy requires a minimum momentum transfer of . If this quantity is comparable to (where is the size of the target nucleus), any additional transverse momentum will be strongly suppressed by the form factor of the target nucleus. As a result, the peak of the differential cross section for nonzero disappears for large ALP masses and small ALP energies.
4 Experimental event rates
Having calculated the rate and distribution of ALPs produced in the beam dump, we now proceed to calculating the probability that these ALPs lead to an observable signal in a given experiment. As the typical experimental setup we assume that the target is immediately followed by an absorber of length and a decay volume of length . At the far end of the decay volume there is a photon detector, which we take to be cylindrically symmetric with radius , such that it covers an angle from the interaction point. A concrete example of such an experiment is NA62, which is shown in figure 5 and discussed in detail in section 6.
The first ingredient for calculating the signal acceptance is the probability that the ALP passes through the absorber without decaying or interacting and then decays inside the decay volume. For an ALP with decay length this probability is given by
(23) 
For it is a good approximation to neglect the second term because almost all ALPs that reach the decay volume will decay before arriving at the detector. The simplest way to estimate the signal acceptance is then to assume that every ALP decaying in the decay region leads to an observable signal in the detector, provided the ALP angle satisfies . We will use this approach in the following section for the analysis of CHARM and NuCal.
For present and future proton beam dump experiments, a more refined analysis strategy will be necessary in order to suppress potential backgrounds and allow the reconstruction of ALP properties in the case of a discovery. Most importantly, we require that both photons produced in the ALP decay are detected, i.e. we require two electromagnetic showers in the calorimeter coincident in time. The signal acceptance then needs to include

the probability that both photons produced in the ALP decay reach the detector located at the far end of the decay volume;

the probability that both photons are detected and give a signal that can be distinguished from potential backgrounds.
The first part depends mostly on the angle of the ALP, , and its boost , which determines the typical opening angle between the two photons produced in the laboratory frame. However, there will also be some dependence on the ALP laboratoryframe decay length , because the separation between the two photons will be smaller the closer the ALP can get to the detector before decaying. The second part of the probability can depend on the photon position and energy. More importantly it will depend on whether the separation between the two photons is large enough to identify two separate photons which are coincident in time. Note that for a typical calorimeter the detection probability for photons is very close to unity across the entire energy range of interest. It is straightforward to determine the detection probability as a function of and from MonteCarlo simulations of ALP decays in a given detector geometry. Nevertheless, before presenting the results of these simulations, it is useful to develop an intuitive understanding of the most important contributions.
First of all, if the direction of the ALP is such that , it is impossible for both photons produced in the ALP decay to reach the detector, so that the detection probability vanishes: . For , on the other hand, the detection probability depends on the expected separation of the two photons. Since the decay of an ALP into photons is isotropic in the ALP rest frame, the distribution of the opening angle between the two photons in the laboratory frame is given by
(24) 
which is strongly peaked towards the minimum opening angle (assuming ). If the ALP lifetime is short compared to the length of the decay volume, such that most ALPs will decay at the beginning of the decay volume, the typical separation of the two photons at the position of the detector will hence be .
This consideration enables us to constrain the acceptance as a function of the boost factor . Clearly, it is a necessary requirement that . Consequently, we conclude . If is significantly different from zero, even larger boost factors will be required in order for both photons to reach the detector. On the other hand, we also require that the separation between the two photons is larger than a certain minimal distance to observe two individual photons. This requirement can be translated into an upper bound on the boost factor: .
In summary, for short ALP decay length, , we can roughly estimate the detection probability as
(25) 
It is clear that this estimate is rather simplistic. The easiest way to obtain a more realistic estimate of the detection probability is to perform a MonteCarlo simulation of ALP decays. The basic idea of such a simulation is to randomly determine the decay point and the direction of the emitted photons for a large number of ALPs with given properties , and and then determine the fraction of events where both photons hit the detector with the required separation.
Once the detection probability is known, we can calculate the number of predicted events by comparing the crosssection for ALP production with the crosssection for protonnucleus scattering . According to Carvalho:2003pza (), this cross section is largely independent of the beam energy and can be parametrized as
(26) 
where is the mass number of the nucleus. The number of ALPinduced events in a given experiment is then
(27) 
where is the number of protons on the target material and
(28) 
is the fiducial cross section, which implicitly depends on the ALP mass and decay constant. Above, we have implicitly made the assumption that the target material is thick enough to absorb all protons to a good approximation.
Before we conclude this section, let us comment on the influence of the target material. Since, according to the EPA, the photons probe the charge of the entire target nucleus, the ALP production crosssection is approximately proportional^{8}^{8}8In fact, the enhancement for heavy nuclei is somewhat suppressed due to their larger size. Smaller photon momenta are sufficient to resolve the substructure and there is thus an additional form factor suppression. Comparing copper and tungsten we find that the naive scaling law is shifted in favour of copper by about 1020%. to . According to eq. (27), the ALP production cross section needs to be normalized to the total interaction rate of the protons with the target material. Thus, using the approximate expression for from Carvalho:2003pza (), the predicted number of events for different target materials but equal acceptance and runtime can roughly be rescaled using a factor , which is equal to 30 for iron, 34 for copper, 52 for molybdenum, 99 for tungsten and 111 for lead. We note that for heavy ions one has an additional enhancement by the charge of the ion. However, typically the flux of ions is considerably lower at least partially eating up this advantage. Nevertheless it may be worthwhile to check for heavy ion experiments with an advantageous geometry.
5 Existing constraints from past experiments
In this section we derive constraints from past proton beam dump experiments using coherent production of ALPs via the Primakoff process and compare them with other exclusion bounds on the ALP parameter space from the literature. The two most constraining experiments are^{9}^{9}9We have also taken NOMAD into consideration for which an analysis for ALPs from decays has been made previously Altegoer:1998qta (). NOMAD has an onaxis detector located away from a Beryllium target and we estimate it to be less constraining than CHARM for the model we are considering. CHARM Bergsma:1985qz (), which has recently been used to derive constraints on light scalars and pseudoscalars with couplings to quarks Bezrukov:2009yw (); Clarke:2013aya (), and NuCal Blumlein:1990ay (), which has been reanalysed for the case of hidden photons Blumlein:2011mv (); Blumlein:2013cua ().
The CHARM experiment has performed a search for hidden particles decaying into photons based on a dataset of protons on a copper target. The detector was placed 480 metres away from the beamdump and was 35 metres long. The detector was in transverse dimension and was placed away from the beam axis, so that it covered about of the full circle around the beam axis. The CHARM experiment required a single electromagnetic shower in the detector and quotes a signal acceptance of . The detection probability was therefore given by (see section 4)
(29) 
Since CHARM observed 0 events, we set a bound at 90% confidence level of events Clarke:2013aya ().
The NuCal experiment made use of the U70 proton beam with a beam energy of . The lower beam energy is compensated by the much smaller distance between target and detector of only and the comparably large detector length of . We follow the analysis strategy of Blumlein:2013cua (), requiring a minimum ALP energy of . With this requirement, the detector acceptance is approximately constant and equal to . Since the detector has a radius of , we then obtain
(30) 
In a dataset of protons on an iron target, NuCal observed 1 event compared to an expectation of 0.3 events. At 90% confidence level, we can therefore exclude any parameter point predicting more than 3.6 events.
The parameter regions excluded by CHARM and NuCal are shown in figure 1. We also show a compilation of other constraints on the ALP parameter space (taken from Jaeckel:2015jla () and references therein), most notably the ones set by the electron beam dump searches SLAC 141 Riordan:1987aw () and SLAC 137 Bjorken:1988as (). As emphasised in section 2 the sensitivity of the different experiments depends decisively on the overall geometry of the different setups and therefore varies significantly between the different experiments under consideration. As a result, our reanalysis of the constraints from CHARM and NuCal improves significantly upon existing constraints. Let us now discuss what further improvements can be achieved with present and nearfuture experiments.
6 Projections for an ALP search at NA62
The fixedtarget experiment NA62 at the CERN SPS aims to measure the branching ratio at the 10 % level within the next few years and is currently taking data. As s constitute a major background for this measurement, overall hermetic detection of photons is crucial. Consequently, NA62 could be adapted to search for ALPs decaying to two photons, and we will elaborate on this possibility.
The NA62 experiment uses protons from the SPS impinging on a beryllium target. The positively charged hadron beam of 75 GeV is then selected via an achromat: Only the wanted momentum component can pass through a set of apertures within the socalled ‘Target Attenuator eXperimental areas’ (TAXes)’, see Chapter 2.1.2 of Hahn:1404985 (). The particles passing the apertures are then guided by magnetic elements to the decay volume downstream, cf. figure 5 and Hahn:1404985 () for details on the detector.
For a potential ALP search, two considerations for the primary beam are important compared to ‘regular’ data taking

To reach an acceptable background level, the experiment must be run in ‘dump mode’. Technically, this is possible by ‘closing’ the two TAXes, which are approximately behind the target. This means that the TAXes are positioned such that no free aperture is being made available for main beam and all secondaries. In this way, almost only muons and neutrinos are able to reach the decay volume (both TAXes are approximately thick).

As argued in section 4, a high fixed target is preferable in producing ALPs through photon fusion. Removing the beryllium target from the beampath allows all 400GeV protons to impinge directly onto the copper of the first TAX, offset by downwards from the straight beam path due to the bend magnet.^{10}^{10}10In fact, both TAXes are equipped with two holed inserts (the apertures for the appropriate component of the secondaries) and one diameter solid tungsten insert, respectively. Impinging on the full tungsten insert would be obviously be even more favourable than impinging on copper. The solid tungsten insert would be large enough in diameter to absorb the full proton beam (which is of a couple of mm in transverse size). However the implemented range of movement for the TAXes does not allow to impinge on the tungsten of the first TAX directly. Doing so would require modifications in the target area and we do not consider this situation further. We also emphasise that, even in the configuration with the target in place, a sizable fraction (40 %) of the protons go through the target without interacting and eventually impinge on the copper of the TAX.
After its production in the TAX, the ALP needs to reach the decay volume (behind the ‘CHarged Anti Counter’). We therefore require that the ALP propagate along a distance of m before decay. As distance over which ALP decays can give a signal in the calorimeter we use m.
As mentioned above, an excellent photon vetoing is crucial for the main physics goal of NA62. For our estimate, we conservatively assume usage of the Liquid Krypton Calorimeter (LKr) for photon detection only. The LKr was previously used in NA48 and is described in detail in Fanti:2007vi (). A hole in the centre of the calorimeter hosts the vacuum tube for the beam ( radius) and the calorimeter itself is segmented in the transverse plane into cells of approximately by . Additional calorimeters (IRC, SAC) are installed to cover the calorimeter hole, but as we will see below the presence of the LKr calorimeter hole has little influence on the ALP acceptance.
Let be the twovectors of the two photons on the calorimeter plane and be their radial distance with respect to the beam axis. Then, to summarize, the photons are required to be

within the calorimeter, thus ;

far enough from the central hole to give a fully contained shower ;

separated by at least to avoid shower overlap.
Note that the transverse shape of the LKr calorimeter is octagonal. For our estimates we implement the decay of ALPs in a toy Monte Carlo simulation unrelated to the NA62 software. In our Monte Carlo we approximate the LKr by a cylinder.
With the requirements above implemented, we can obtain a realistic estimate of the detection probability (an explicit dependence on and only arises after folding the acceptance with the production crosssection). While it would be possible to consider photon losses in the spectrometer or Ring Imaging Cherenchov detector as well as detection efficiencies smaller than unity, we expect these effects to be small and therefore do not include them in our MonteCarlo simulation. Figure 6 shows the resulting detection probability as a function of the ALP angle and its boost factor for different values of the decay length. The largest acceptances are found for offangle ALPs and , and and can reach up to .
In the left panel of figure 6, where the condition is satisfied, we also show for comparison the estimate for the acceptance region derived in section 6. We find good agreement between our simple estimate and the results of the more detailed MonteCarlo simulation. However, for longer decay lengths, as considered in the right panel, the acceptance region becomes much larger than predicted by our estimate, because a significant fraction of the ALPs decays closer to the detector. In addition to the features discussed in section 4 one can observe a loss of sensitivity for and very large boost factors due to one or both photons being lost in the calorimeter hole. In agreement with expectations, this effect becomes relevant if both and are small compared to , implying and . Since the ALP production cross section is very small for (see figure 3), the presence of the calorimeter hole does not significantly affect the sensitivity of NA62 to ALPs.
In fact, the angular distribution of ALPs produced via the Primakoff process peaks at somewhat larger angles and smaller boost factors than what can be detected by the NA62 detector, implying that an ideal detector should be somewhat larger and closer to the target. Nevertheless, highenergy ALPs are typically emitted with very small angle , such that they match the acceptance window of the NA62 experiment. We show in figure 7 the fiducial cross section defined in eq. (28). While the left panel shows a typical distribution as obtained in large regions of parameter space, the right panel corresponds to ALPs with a rather short lifetime, such that only particles produced with very large boost factor reach the decay volume.
To calculate the predicted number of events in NA62, we take as nominal intensity protons per effective second on target. The effective flat top of the proton spill is approximately 3 seconds, 2 times every supercycle of 42 seconds, thus we have to multiply by an effective factor of 0.15. Within 24 hours, therefore protons on target are available. Assuming that all backgrounds can be rejected, NA62 can probe all those points in the (, ) parameter plane where at least three events are expected from the production and decay of ALPs. Figure 8 shows the resulting expected sensitivity for a datataking period of 1 day and 1 effective^{11}^{11}11Here, we do not account for detector or beamdowntime. So the actual runtime could be larger than a month. month, respectively. As can be seen from the right panel, such a modest integrated intensity is fully sufficient to significantly improve upon existing upper limits from ALPsearches at beam dump experiments and probe presently unexplored regions of parameter space. We emphasize that doing so does not require any significant modification of the existing experimental setup.
7 Sensitivity at the proposed SHiP facility
While NA62 offers a unique opportunity to probe unexplored parameter space with an existing experiment, there are various ways in which to further improve the sensitivity of future experiments. As discussed above, the production rate is enhanced by a factor for heavy elements in the target. Furthermore, a shorter distance between the target and the beginning of the decay volume would allow to probe ALPs with shorter lifetimes, corresponding to larger couplings. Placing the detector closer to the target has the additional advantage that it is possible to cover somewhat larger ALP production angles and make use of the fact that the differential cross section typically peaks at around .
All of these potential improvements have been incorporated in the SHiP proposal Anelli:2015pba (), which aims to search for a variety of weaklycoupled lowmass states with very high intensity beams Alekhin:2015byh (). In this section we present an estimate of the expected sensitivity of an ALPsearch at SHiP.^{12}^{12}12Previous estimates of the sensitivity of SHiP for ALPs were based on incoherent production of ALPs, i.e. Alekhin:2015byh (). Including the coherent production via Primakoff processes significantly increases the parameter region that can be probed by SHiP. We assume a total of protons on a molybdenum target.
The decay volume is proposed to begin approximately at a distance of from the target and is long. We consider the detector to have a radius of , such that SHiP can cover production angles up to . As for NA62, we require that both photons from the ALP decay reach the detector and that the separation between them is at least . Moreover, we require that the combined energy of both photons satisfies .
An estimate of potential backgrounds is beyond the scope of this work. We therefore assume backgrounds to be negligible and show the parameter region that would lead to at least three events in the SHiP experiment in figure 9. The combination of an optimised design and a very large number of protons on target allows SHiP to explore large parameter regions inaccessible for other kinds of searches. The projected sensitivity extends all the way up to ALP masses of about and ALPphoton couplings as small as .
8 Conclusions
Proton fixed target experiments provide exciting opportunities in the search for new light and weakly coupled particles with masses in the MeV to GeV range. However, they confront us with the difficulty of colliding two composite particles, the proton and the even more complicated atoms of the target material. This raises the question whether one can provide a reliable calculation of the production rates. In this work we have addressed this problem by studying the Primakoff production of axionlike particles (ALPs) in proton fixed target experiments.
For a highly energetic proton beam of 400 GeV, as provided by the SPS, and ALP masses in the MeV to GeV mass range significant production is possible from photons with sufficiently high momentum to not be affected by the electron shell, and sufficiently low momentum that they do not resolve the proton and the nucleus. In this range a reliable calculation can be performed by taking into account simple electromagnetic form factors. As an additional benefit the scattering is then coherent over the whole nucleus and production cross sections are enhanced by a factor of the nuclear charge squared.
Importantly, we have provided predictions for the angular distribution of the ALPproduction crosssection: Although transversal momenta of the produced ALPs are typically small, the detector in a beam dump experiment is placed far away from the target and therefore covers only a tiny angle from the production point. A precise determination of the expected spatial distribution for ALPinduced events and an accurate estimate of the geometric acceptance is therefore mandatory in preparing and analysing a real experimental run. Taking all these effects into account, we have shown that even with a rather modest beamtime requirement, the currently operating NA62 experiment would have a sizeable discovery potential for ALPs in the mass range of MeV. The proposed facility SHiP could extend this reach over the course of its running period up to masses of .
In the present work we have focussed on pseudoscalar ALPs that couple dominantly to photons. It is however straightforward to generalise our results to scalar ALPs as well as ALPs with additional couplings to fermions as follows:

Writing the coupling between scalar ALPs and photons as , we obtain identical expressions for the ALP lifetime and the ALP production cross sections as for the case of the pseudoscalar. Our analysis therefore applies to this case as well.

Even for ALPs with relatively large (derivative) couplings to fermions, the decay into photons will typically give the dominant contribution to the ALP decay length for Dolan:2014ska (), which is the region of interest for NA62. For larger ALP masses, as potentially testable at SHiP, decays into muons (and, for scalar ALPs, mesons) can significantly reduce the ALP decay length and hence suppress event rates.

In the presence of couplings to fermions, new production modes become available, such as ALPstrahlung, ALPpion mixing or ALP production in flavourchanging rare decays. Depending on the details of the model under consideration these processes may significantly boost the ALP production rate and therefore enhance the sensitivity of beam dump experiments.
In conclusion, we have worked out a detailed example how physics beyond the SM at comparably low mass could be probed with existing setups with comparably low effort, providing a complementary window to the opportunities that the highestenergy accelerators are offering us.
We would like to thank Walter Bonivento, David d’Enterria, Javier Redondo and Pedro Schwaller for useful discussions, and Johannes Blümlein, Matthew Dolan and Tommaso Spadaro for insightful comments on the manuscript. BD acknowledges very helpful conversations within the NA62 collaboration, particularly with Brigitte BlochDevaux, Augusto Ceccucci, Luigi Di Lella, Niels Doble, Lau Gatignon, Evgueni Goudzovski, Matthew Moulson, Mathieu PerrinTerrin, Mauro Raggi, Giuseppe Ruggiero, Tommaso Spadaro and Paolo Valente. JJ would like to thank the IPPP for a very enjoyable visit. FK and KSH are grateful to the CERN Theory Division for hospitality. This work is supported by the German Science Foundation (DFG) under the Collaborative Research Center (SFB) 676 “Particles, Strings and the Early Universe”, the Transregio TR33 “The Dark Universe”, as well as the ERC Starting Grant ‘NewAve’ (638528).
Appendix A Cross sections
Let us first calculate the cross section for the process . The matrix element for the photonphotonALP vertex is given by
(31) 
where is the photon polarisation vector. Averaging over initial state polarisations gives
(32) 
We hence obtain
(33) 
Evaluating this expression in the centreofmass frame yields
(34) 
We now want to calculate the cross section for the process , allowing nonvanishing transverse momentum for both the incoming photon and the outgoing ALP. In the lab frame (where the incoming nucleus is at rest), we adopt the following choice of coordinates:
(35) 
The 4momentum of the outgoing nucleus is then fully determined by energy and momentum conservation. Note also, that we require and , i.e. we require all particles to be onshell.^{13}^{13}13In principle, the incoming photon can be offshell. However, the dominant contribution will result from incoming photons that are (almost) onshell, so the chosen approximation should be valid. The photon momentum is then fully determined by specifying the energy and direction of the outgoing ALP, the transverse momentum of the photon and the angle between the transverse momentum of the photon and the transverse momentum of the ALP.
Under the assumption that and , we obtain to order^{14}^{14}14Note that the term proportional to is typically tiny compared to the terms proportional to and . We keep these terms nevertheless in order to be able to compare our final result with the one obtained in Bjorken:1988as () for . , and
(36) 
It is then possible to calculate all products of fourvectors at this order. In particular, we find
(37) 
Finally, we obtain for the squared matrix element
(38) 
To obtain the scattering cross section, we need to multiply this expression with the nuclear form factor and a phase space factor
(39) 
Reassuringly, this expression reduces to the one from Bjorken:1988as () upon setting .
Appendix B Angular averaging
Let us consider two random 2dimensional vectors and . The probability distribution of the magnitude of these two vectors is given by the functions and . The directions and of the two vectors are evenly distributed (i.e. isotropic). We are now interested in determining the probability distribution for the magnitude of . As a first step we define
(40) 
The distribution can now be written as
(41) 
Now we define the new angular variables and , noting that . Clearly, is independent of and hence we can directly perform the integration:
(42) 
The integral over vanishes unless , and satisfy the triangle inequalities: . With the range of integration thus restricted, we can solve
(43) 
Now we need to differentiate the argument of the function with respect to and substitute this solution:
(44) 
Finally, we note that the argument of the function has two roots within , leading to an additional factor of 2. We thus obtain
(45) 
It is an easy but useful cross check to integrate the righthand side over to make sure that is correctly normalised.
In some applications it may be useful to write the range of integration in a slightly different way. Clearly, we can integrate separately over the two regions and . In the first case, the triangle inequalities imply and . The integration can hence be written as
(46) 
References
 (1) NA62 Collaboration, NA62: Technical Design Document, 2010. Tech. Rep. NA621007, https://cds.cern.ch/record/1404985.
 (2) SHiP Collaboration, M. Anelli et al., A facility to Search for Hidden Particles (SHiP) at the CERN SPS, 1504.04956.
 (3) S. Alekhin et al., A facility to Search for Hidden Particles at the CERN SPS: the SHiP physics case, 1504.04855.
 (4) J. D. Bjorken, S. Ecklund, W. R. Nelson, A. Abashian, C. Church, et al., Search for Neutral Metastable Penetrating Particles Produced in the SLAC Beam Dump, Phys. Rev. D38 (1988) 3375.
 (5) J. Blümlein et al., Limits on neutral light scalar and pseudoscalar particles in a proton beam dump experiment, Z. Phys. C51 (1991) 341–350.
 (6) J. D. Bjorken, R. Essig, P. Schuster, and N. Toro, New FixedTarget Experiments to Search for Dark Gauge Forces, Phys. Rev. D80 (2009) 075018, [0906.0580].
 (7) S. Andreas, C. Niebuhr, and A. Ringwald, New Limits on Hidden Photons from Past Electron Beam Dumps, Phys. Rev. D86 (2012) 095019, [1209.6083].
 (8) A. Konaka et al., Search for Neutral Particles in Electron Beam Dump Experiment, Phys. Rev. Lett. 57 (1986) 659.
 (9) E. M. Riordan et al., A Search for Short Lived Axions in an Electron Beam Dump Experiment, Phys. Rev. Lett. 59 (1987) 755.
 (10) M. Davier and H. Nguyen Ngoc, An Unambiguous Search for a Light Higgs Boson, Phys. Lett. B229 (1989) 150.
 (11) J. Jaeckel and M. Spannowsky, Probing MeV to 90 GeV axionlike particles with LEP and LHC, 1509.00476.
 (12) CHARM Collaboration, F. Bergsma et al., Search for Axion Like Particle Production in 400GeV Proton  Copper Interactions, Phys. Lett. B157 (1985) 458.
 (13) B. Batell, M. Pospelov, and A. Ritz, Exploring Portals to a Hidden Sector Through Fixed Targets, Phys. Rev. D80 (2009) 095024, [0906.5614].
 (14) A. Halprin, C. M. Andersen, and H. Primakoff, Photonic Decay Rates and NuclearCoulombField Coherent Production Processes, Phys. Rev. 152 (1966) 1295–1303.
 (15) Y. Nomura and J. Thaler, Dark Matter through the Axion Portal, Phys. Rev. D79 (2009) 075008, [0810.5397].
 (16) C. Boehm, M. J. Dolan, C. McCabe, M. Spannowsky, and C. J. Wallace, Extended gammaray emission from Coy Dark Matter, JCAP 1405 (2014) 009, [1401.6458].
 (17) A. Berlin, S. Gori, T. Lin, and L.T. Wang, Pseudoscalar Portal Dark Matter, Phys. Rev. D92 (2015) 015005, [1502.06000].
 (18) J. L. Hewett et al., Fundamental Physics at the Intensity Frontier, 1205.2671.
 (19) G. Hiller, B physics signals of the lightest CP odd Higgs in the NMSSM at large tan beta, Phys. Rev. D70 (2004) 034018, [hepph/0404220].
 (20) M. Freytsis, Z. Ligeti, and J. Thaler, Constraining the Axion Portal with , Phys. Rev. D81 (2010) 034001, [0911.5355].
 (21) S. Andreas, O. Lebedev, S. RamosSanchez, and A. Ringwald, Constraints on a very light CPodd Higgs of the NMSSM and other axionlike particles, JHEP 08 (2010) 003, [1005.3978].
 (22) M. J. Dolan, F. Kahlhoefer, C. McCabe, and K. SchmidtHoberg, A taste of dark matter: Flavour constraints on pseudoscalar mediators, JHEP 03 (2015) 171, [1412.5174]. [Erratum: JHEP07,103(2015)].
 (23) J. Blümlein and J. Brunner, New Exclusion Limits for Dark Gauge Forces from BeamDump Data, Phys. Lett. B701 (2011) 155–159, [1104.2747].
 (24) J. Blümlein and J. Brunner, New Exclusion Limits on Dark Gauge Forces from Proton Bremsstrahlung in BeamDump Data, Phys. Lett. B731 (2014) 320–326, [1311.3870].
 (25) J. Alwall, M. Herquet, F. Maltoni, O. Mattelaer, and T. Stelzer, MadGraph 5 : Going Beyond, JHEP 06 (2011) 128, [1106.0522].
 (26) 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–281.
 (27) A. D. Martin and M. G. Ryskin, The photon PDF of the proton, Eur. Phys. J. C74 (2014) 3040, [1406.2118].
 (28) 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–112.
 (29) B. L. Henke, E. M. Gullikson, and J. C. Davis, XRay Interactions: Photoabsorption, Scattering, Transmission, and Reflection at E = 5030,000 eV, Z = 192, Atom. Data Nucl. Data Tabl. 54 (1993), no. 2 181–342.
 (30) “NIST XRay Form Factor, Atten. Scatt. Tables Form Page.” http://physics.nist.gov/PhysRefData/FFast/html/form.html.
 (31) J. Carvalho, Compilation of cross sections for proton nucleus interactions at the HERA energy, Nucl. Phys. A725 (2003) 269–275.
 (32) NOMAD Collaboration, J. Altegoer et al., Search for a new gauge boson in pi0 decays, Phys. Lett. B428 (1998) 197–205, [hepex/9804003].
 (33) F. Bezrukov and D. Gorbunov, Light inflaton Hunter’s Guide, JHEP 05 (2010) 010, [0912.0390].
 (34) J. D. Clarke, R. Foot, and R. R. Volkas, Phenomenology of a very light scalar (100 MeV 10 GeV) mixing with the SM Higgs, JHEP 02 (2014) 123, [1310.8042].
 (35) NA62 Collaboration, NA62 Detector Schematics, 2010. NA62PHOCGI2010001, http://cds.cern.ch/record/1293104.
 (36) N. Doble, “Beam Documentation for NA48/NA62.” http://doble.web.cern.ch/doble/.
 (37) NA48 Collaboration, V. Fanti et al., The Beam and detector for the NA48 neutral kaon CP violations experiment at CERN, Nucl. Instrum. Meth. A574 (2007) 433–471.