Searching for the axion-like particle at the EIC

The axion-like particle (ALP) is a well motivated new particle candidate for beyond the standard model. In this study, we propose to probe the ALP via photon fusion scattering at the upcoming Electron-Ion Collider (EIC) with electron and proton energies of and , respectively. We can constrain the effective coupling strength between the ALP and photons to be at the confidence level with an integrated luminosity of for the mass range . Such bounds may be significantly improved if we consider the nucleus beam at the EIC. We also demonstrate that the limits from the EIC can be stronger than the off Z-pole measurement at the LEP and light-by-light scattering with pp collisions at the LHC.


I. INTRODUCTION
Axion-like particles (ALPs) are widely predicted in new physics beyond the standard model (BSM). A well known example is the Pseudo-Nambu-Goldstone boson (PNGB) from new global symmetry breaking, which was designed to solve the strong CP problem [1,2]. The ALP has received considerable attention in the particle physics and cosmology communities because it can solve the naturalness problems [3,4] and may be a compelling dark matter candidate for the universe [5][6][7]. The landscape of ALPs is rich, and the phenomenology is determined by their mass and couplings with SM particles. In general, ALPs may couple to gauge bosons, fermions, and Higgs bosons (see Refs. [8][9][10][11][12][13] for a general discussion). Many dedicated experiments have been proposed to search for ALPs based on their mass region and couplings. For example, the interactions between ALPs and fermions can be probed in rare decays [10,[14][15][16][17][18][19][20][21]. The bound of couplings to gauge bosons could be constrained by astrophysics and cosmology observations, for example, stellar evolution, big bang nucleosynthesis, and anisotropies in the cosmic microwave background for a light ALP [22][23][24][25], and by collider searches for heavy ALPs with light-by-light or vector-boson fusion type scattering .
BR(a → γγ) = 1 In this study, we consider the case in which the ALP a predominantly couples to photons, so that the branching ratio and the effective Lagrangian can be parameterized as where a is the ALP with mass , and is the (dual) field strength tensor of the photon. The coupling strength has been seriously constrained by electron/ proton beam-dump experiments [50,51], +invisible [30], inclusive [26], a photon-beam experiment [52], and at Belle II [33] with a mass range from MeV to 10 GeV. Above 10 GeV, the parameter space can be probed by light-by-light scattering in heavy-ion collisions at the Large Hadron Collider e + e − → γa(→ γγ) (LHC) energy [31,49], at the LEP [27], and measurements at the LHC [39][40][41][42][43]48] and future lepton colliders [44][45][46][47]. However, the bound from GeV to tens of GeV is less limited by current experiments compared to the MeV (or smaller) mass range. Such light ALPs are desirable for the community because it could be a natural product of new global symmetry breaking if it is a PNGB. In this paper, we propose to search for ALPs via photon fusion production at the upcoming Electron-Ion Collider (EIC) (see Fig. 1) 1) , which could be complementary to measurements at the LHC, and heavy-ion collisions at the LHC energy and lepton colliders when probing the ALP in this mass range. The ALP may be produced through s-channel and tchannel photon fusion scattering at the EIC 2) (see Fig. 1).

II. PRODUCTION OF ALPs
The production rate of the signal depends on the coupling strength , ALP mass , and branching ratio , that is, where denotes the cross section from the s-channel and t-channel, respectively. The narrow width approximation is applied for s-channel scattering because the decay width [54]. We check that the interference effects between the signal and background is negligible and can be ignored in our analysis. For simplicity, we take the branching ratio in this study.  Below, we utilize MadGraph5 [55] to calculate the signal cross section with electron and proton energies of and , respectively, at the leading order with NNPDF sets [56]. The factorization and renormalization scales in our calculation are chosen as the default scale in MadGraph5, that is, the transverse mass. To avoid the soft and collinear divergence, the following kinematic cuts for the final states are applied: where and with denote the transverse momentum and pseudorapidity of particle m, respectively.
The cone distance , where denotes the azimuthal angle of particle m. Figure 2 shows the production rates for the processes (red dashed line), (blue dashed line), and (green dashed line) without any kinematic cuts for the photons in the final state. It clearly shows that the cross section of the signal is dominantly determined by s-channel scattering and the contribution from the t-channel is negligible because of the suppression of the phase space. In addition, the narrow width approximation for s-channel production works well in this process. In the same figure, the solid cyan and black lines represent the cross sections from and , respectively, after we impose the kinematic cuts in Eq. (3) for the photons. Note that the kinematic cuts of the photons decrease the cross section significantly when the ALP mass . This arises because the invariant mass of a photon pair from an on-shell ALP decay is given by  . The solid and dashed lines correspond to the production rate with and without the kinematic cuts (see Eq. (3)) for the photons in the final states, respectively. We consider both the s-channel and t-channel for the red line, whereas only the contribution from the s-channel is considered for the green and blue lines.

Yandong Liu, Bin Yan
Chin. Phys. C 47, 043113 (2023) 1) The elastic scattering of two photons could also be used to search the ALP, but the signature is different from our scenario [53].

aγγ aZγ aZZ aZγ aZZ
2) The gauge-invariant operators can generate the effective couplings , and at the same time, but the contributions from and will be highly suppressed by the propagator of Z boson and can be ignored.

043113-2
, where denotes the transverse momentum of the photons. For a heavy ALP, the kinematic cuts could be satisfied automatically, whereas there is a phase space suppression effect for light ALPs, which is induced by the kinematic cuts of the photons (see Eq. (3)). Such an effect also generates a peak at approximatelt for the production rate distributions of the processes with the cuts on the photons (see the cyan and black lines).

III. COLLIDER SIMULATION
Next, we perform a detailed Monte Carlo simulation to explore the potential of probing the ALP at the EIC. The main irreducible backgrounds arise from the processes and . The cross sections after including the kinematic cuts in Eq. (3) are and . We also take into account the possibility that an electron or jet is misidentified as a photon in this study. The reducible backgrounds may originate from the processes and . The numbers shown inside the brackets denote the production cross section after imposing the cuts in Eq. (3). The other backgrounds involving multi-electrons and jets, for example, and , are negligible when we consider the basic cuts in Eq. (3) and mistag efficiencies. We generate both the signal and backgrounds using MadGraph5 [55] with the kinematic cuts in Eq. (3). Parton level events are passed to PYTHIA8 [57] for parton showering and hadronization, and the detector effects are simulated by Delphes [58]. In the detector simulation, we use the EIC delphes card, which was generated by Arratia and Sekula [59], based on the parameters in Ref. [60] and utilized in Ref. [61][62][63]. The detector parameters in Ref. [60] dictate the tracking momentum resolution, secondary-vertex resolutions, calorimeter energy resolutions, and particle identification performance at the EIC. Among them, the photon energy resolution plays a key role in our simulation, which is , with , , and . Now, we require the following set of preselection cuts on the reconstruction objects: m, n = γ, e, j where denote the photon, electron, and jet from the detector reconstruction, respectively. Note that the kinematic threshold for the reconstruction objects at the EIC may be looser compared to that at the LHC owing to the significantly lower collider energy (see the cuts of HERA as a Ref. [64]).
In the signal events, photons arise from the decay of m a the ALP and thus exhibit a peak at approximately in the invariant mass distribution. However, the photons of the backgrounds originate from the radiation of the electron and quark; thus, the peak position of the invariant mass distribution from the photon pair is determined by the transverse momentum of the photons; see Fig. 3 for the normalized invariant mass distributions of the photon pair from the signal and backgrounds. Because the typical decay width of the ALP is considerably smaller than the resolution of the di-photon invariant mass at the EIC [60], we further require the invariant mass of the first two leading photons within the mass window, We also check that other kinematic observables (for example, ) cannot significantly improve the cut efficiency in this process. This arises from the fact that both the signal and backgrounds share a similar topology to the Feynman diagram, and as a result, both should exhibit a similar distribution. However, the photons share the energy from the ALP in the signal process; thus, the distribution should peak at , whereas for the backgrounds tend to produce a soft spectrum because the cross sections are enhanced by the soft and/or collinear singularity in that phase space region. However, information may be correlated to ; therefore, the additional cut on cannot improve the sensitivity to probe the ALP. We show the cut efficiencies with several benchmark ALP masses after the kinematic cuts from Eqs. (4) and (5)   those of the irreducible backgrounds ( and ). However, the reducible backgrounds would still be dominated after considering the production rates of these processes. Because the kinematic cuts in Eq. (4) may be satisfied automatically for a heavy ALP, the cut efficiency for the signal tends to be constant when .
Equipped with the signal and background production cross sections and collider simulation efficiencies, the upper limit on the effective coupling at the confidence level can be obtained in terms of [65] √ where and are the numbers of signal and background events, respectively. Given the integrated luminosity of 300 fb , the upper limit on the effective coupling is presented in Fig. 4 (orange region). It shows that can be constrained to be TeV at the confidence level by assuming null results when directly searching for ALPs with at the EIC. This result may be further improved if we consider the nucleus beam at the EIC because the cross section would be enhanced by the atomic weight , where Z and N are the number of protons and neutrons in nucleus, respectively [53]. To roughly estimate the upper limit on from the nucleus beam, we use lead (Pb) as an example. The cross sections of the signal and backgrounds from the Pb beam could be obtained by properly rescaling the cross sections at electron-proton collisions, that is, , with for the Pb beam 1) . Although the mixing of protons and neutrons in the nucleus will change the total density for up and down quarks in the nuclear PDFs relative to those of the proton, the total effects will mildly modify the total cross sections and kinematic distributions [66], which do not significantly alter the conclusions of this study. Under this approximation, we find that the effective coupling from the Pb beam may be improved by several times compared to that of the proton beam; see the gray region of Fig. 4.

IV. SUMMARY AND DISCUSSION
As discussed in Sec. I, ALPs with masses of tens of GeV have been widely searched for by colliders. These searches include the measurements of tri-photons on and off the Z-pole at the LEP [67,68], searches for the same final states at the LHC [27], and light-bylight scattering in heavy-ion collisions at the LHC energy and pp collisions at the LHC [31,41,49]. However, the bounds of from several of the above measurements depend on the assumption of the interaction. This arises from the fact that both couplings may be generated from the dimension-5 operators and , where and are the field strength tensors of and , and and g are the corresponding gauge couplings. The effective coupling strengths of and are and , respectively. Therefore, the   We compare the investigation of at the EIC with other measurements in Fig. 4. It is evident that the Z-pole measurement (yellow region), compared with the other processes, yields the strongest constraint on the value of . However, this conclusion is only available when we consider one operator in the analysis. We also notice that the expected limit from the EIC with a proton beam (orange region) may be stronger than the off Z-pole measurement (green region) at the LEP and light-by-light scattering at the LHC (pink region), whereas the expected limit from the Pb beam at the EIC (gray region) may be comparable to the measurement in Pb-Pb collisions (blue region). Finally, we emphasize that photon fusion production at the EIC is complementary to other processes in the measurement of ALP-photon coupling.