Probing the Dark Axion Portal with Muon Anomalous Magnetic Moment

We propose a new scenario of using the dark axion portal at one-loop level to explain the recently observed muon anomalous magnetic moment by the Fermilab Muon g-2 experiment. Both axion/axion-like particle (ALP) and dark photon are involved in the same vertex with photon. Although ALP or dark photon alone cannot explain muon $g-2$, since the former provides only negative contribution while the latter has very much constrained parameter space, dark axion portal can save the situation and significantly extend the allowed parameter space. The observed muon anomalous magnetic moment provides a robust probe of the dark axion portal scenario.

From the BNL result to the Fermilab one, both the central value and the uncertainty decreases. Huge amount of work has been done to match the unprecedented precision. The SM contribution to a µ contains four parts [8], The first two are the QED and electroweak (EW) predictions, respectively, while a HVP µ is the hadronic vacuum polarization (HVP) and a HLbL µ the hadronic light-by-light (HLbL) contribution. Although the biggest source of uncertainty comes from the hadronic part [9][10][11][12], the most recent calculations [13][14][15] have included the updated measurement of the hadronic contributions [16][17][18]. The latest theoretical calculation [21] gives where the uncertainty mainly comes from the hadronic vacuum polarization a HVP µ and the light-by-light part a HLbL µ .
In this letter, we explore the possibility that the dark axion portal [56] with coupling among ALP a, photon γ, and a massive dark photon γ can explain the observed muon anomalous magnetic moment at the Fermilab Muon g-2 experiment. Since this dimension-5 operator was proposed only recently and involves two invisible particles, its coupling C aγγ is not strongly constrained yet. With TeV scale new physics, C aγγ ∼ 3/TeV, a sizable parameter space is still available as we will elaborate in this letter.
The Dark Axion Portal Contribution -The dark axion portal [56] establishes the connection between the visible sector with the dark one via not just a single ALP or a single photon but both of them, where Through this dimension-5 operator, the CP-violating ALP a couples with the dual field strength of dark photon X µν ≡ 1 2 µναβ X αβ where X µν ≡ ∂ µ X ν − ∂ ν X µ . As shown in Fig. 1, the dark axion portal can contribute to the muon anomalous magnetic moment if the ALP and dark photon also couple with muon, L y µ a aμ(iγ 5 )µ − eμγ ν µX ν .
Here, y µ a is the Yukawa coupling with ALP while is the kinetic mixing between the dark photon and photon, 1  photons can also contribute by replacing the dark photon in Fig. 1 with a photon [34,57]. However, due to its stringent constraint, we omit this diagram for simplicity. The contribution of dark axion portal depicted in Fig. 1 is divergent. With the cut-off regularization, the result can be expressed in terms of the ultra-violet (UV) scale Λ, where the loop function G is as a function of the ALP mass m a , the dark photon mass m γ , and the muon mass m µ . The dark axion portal contribution (8a) has linear dependence on the Yukawa coupling y µ a . When y µ a is larger than the SM counterpart m µ /v where v ≈ 246 GeV is the Higgs vacuum expectation value, the dark axion portal contribution can be enhanced in comparison with the SM one. Similar feature has been observed and named as chiral enhancement in many models but with quadratic dependence [58][59][60].
Since the interaction in (6) is non-renormalizable, it is only valid up to some cut-off scale ∼ C −1 aγγ . In addition, the divergent loop integral can be assumingly regularized by a similar cut-off Λ with origin from the same UV physics. However, the predicted a µ in (8a) has mild dependence on Λ which only appears in a log function. Orders of variation in Λ can only change a µ by several times which can be easily compensated by tuning couplings. For comparison, the Yukawa coupling y µ a and the kinetic mixing parameter are dimensionless, hence cannot directly reflect the new physics scale.
The dark axion portal can also contribute to the muon anomalous magnetic moment at two-loop level via the photon vacuum polarization sub-diagram [61]. However, this contribution is always negative and numerically negligible. For example, with C aγγ = 3 TeV −1 , the two-loop contribution is roughly two orders of magnitude smaller than its one-loop counterpart. Therefore, we neglect their contribution and focus on the one-loop diagrams in this letter. Note that the dark axion portal coupling, C aγγ = 3 TeV −1 , adopted here satisfies existing experimental constraints [61].
In principle, the ALP can also couple with the SM Z boson, 1 2 C aγZ aF µνZ µν [62,63]. Then a similar contribution from the ALP-photon-Z vertex, by replacing the dark photon γ in Fig. 1 with Z. The analytical formula (8) still applies after replacing the dark photon mass m γ by the Z boson mass m Z , the coupling constants C aγγ by C aγZ and e by g V = g c W 1 4 − s 2 W ≈ 4.5 × 10 −2 e of the vector part of the Z-muon coupling while the axial-vector part does not contribute due to the mismatch of parity and charge conjugation properties. The good thing is that the Z coupling with muon and the Z boson mass have already been measured, hence reducing the number of parameters by two. However, the current bound on the coupling C aγZ 0.03 TeV −1 [64] from the anomalous Z decay Z → γa [65] is rather stringent. The contribution from the ALP-photon-Z vertex is negligibly small.
The Individual Contribution of ALP or Dark Photon -The ALP or dark photon alone can also contribute to the muon anomalous magnetic moment as shown in Fig. 2. We first consider the contribution from the ALP which is finite, Note that this ALP-only contribution is negative [34] since F a (η) 0 where η ≡ m µ /m a . The pseudoscalar case is completely different from the scalar scenario which can contribute a positive term. A pseudoscalar alone cannot explain why the observed a exp µ is larger than the SM prediction a SM µ unless the experimental measurement is smaller than the theoretical prediction.
with η ≡ m µ /m γ . Different from the loop function F a for the ALP-only contribution, F γ 0 always holds. It seems that the single contribution from the dark photon can explain the observed muon anomalous magnetic moment. However, the required parameter space in the m γ − plane to explain the observed ∆a µ has already been excluded by other experimental bounds [49] as we will discuss in detail below.
Parameter Space -It is instructive to compare the three distinct contributions: (8a) for dark axion portal, (9) for ALP, and (10) for dark photon, where C aγγ ∼ TeV −1 . The loop function ratios, G/F a and G/F γ , are dropped out since they are comparable with each other. For the dark axion portal contribution to dominate, a µ a a µ , a γ µ , the ALP and dark photon masses are bounded from below, m a y µ a / 10 GeV and m γ /y µ a √ 10 GeV. If the two dimensionless couplings are comparable with each other, ∼ y µ a , both the ALP mass m a and the dark photon mass m γ are around the GeV scale. Richer mass patterns can be realized by tuning the two couplings to change the mass limits, or allowing the ALP-only and dark photon-only contributions in Fig. 2 to be comparable with the dark axion portal one in Fig. 1. Below we explore the allowed 2)e (purple) [46]; 2) resonant production of dark photon at the BaBar experiment (green) [73]; 3) dark photon from pion decay at the NA48 experiment (red) [74]; 4) dark photon production from various mesons at LHCb [75,76] and Higgs at CMS (blue) [77]; 5) electroweak precision observables (yellow) [78]. The hashed band is the allowed region for explaining the muon g-2 with dark photon alone while the black contour is obtained together with dark axion portal at 95% C.L. parameter spaces in detail.
As argued above, the interesting ALP mass is around GeV to a few hundreds of GeV scale. In this range, the experimental constraints [66][67][68] mainly come from SN1987a, beam dump experiment NA62, low energy electron positron colliders such as BaBar, and collider searches at LHC. We summarize in Fig. 3 those constraints that can apply to the configuration in this letter, where the ALP only couples with muon rather than electron or tau.
The purple region at the left-bottom corner of Fig. 3 is excluded by the supernova (SN) cooling rate from the SN1987a observation [69]. The bound reaches y µ a ≈ 10 −6 and can exclude the mass region up to 0.2 GeV.
The red region in Fig. 3 is excluded by the NA62 experiment [66] using the search of rare kaon decay channel K → µνa. The probe of m a is limited by the kaon mass (∼ 494 MeV), explaining why the excluded region can only extend to 400 MeV. For comparison, we also show the projected sensitivity from the µ + N → µ + N + a process at the NA64 muon beam dump experiment (brown dashed line) [67].
The green region in Fig. 3 was obtained from the BaBar experiment with final-state radiation of the invisible Z which further decays into a µ + µ − pair, e + e − → µ − µ + Z → µ − µ + µ − µ + . Although this bound is not originally obtained for ALP, it can be easily converted [66,70]. The BaBar experiment is an electron-positron collider with center-of-mass energy around 10 GeV. Consequently, the sensitive mass region is m a ∈ [0.1, 4] GeV.
The similar situation happens for the CMS rare Z decay constraint [71] shown as the blue region in Fig. 3 which covers the range from 5 GeV up to 50 GeV. For comparison, we also show the projected sensitivities at Belle II (green dashed line) and HL-LHC [68] (blue dashed line).
It is evident that there is a very large region of the parameter space to be explored for masses above 0.2 GeV. All other constraints compiled in [66][67][68] involves Yukawa couplings with either electron or tau leptons and hence cannot apply to the configuration considered in this letter. The gap between the BaBar and CMS region can be covered by the future HL-LHC searches (blue dashed line) and Belle II can also increase the sensitivity below 4 GeV (green dashed lines) [68]. Even so, the available parameter space is still quite sizable.
The current bounds [49,72] on the dark photon mass m γ and its kinetic mixing with photon has been summarized in Fig. 4. These constraints covers the dark photon mass range from 1 MeV to 1 TeV. In contrast to the ALP case, since the kinetic mixing leads to universal coupling between the dark photon and all charged leptons, those experimental constraints involving electron can also apply here. It is interesting to see that the electron anomalous magnetic moment, (g − 2) e , excludes the very light dark photon scenario (purple region) [46]. The dark photon can be resonantly produced at the BaBar experiment. With 10 GeV center-of-mass energy, the excluded region (green) spans from around 20 MeV up to roughly 10 GeV [73]. In between, the NA48 searches for dark photon from pion decay cover the red region from 9 MeV to 100 MeV [74]. The decays of various mesons at LHCb [75,76] and Higgs at CMS [77] also give strong constraints shown as gaped blue regions. At CMS, the dark photon is produced from Higgs decay, h → Zγ and h → γ γ , and it further decays to a pair of muons. The yellow region comes from the electroweak (EW) precision observables [78]. It is interesting to observe that the EW precision observables fill the gap around Z boson mass.
Revival with Dark Axion Portal -With the available parameter space compiled in Fig. 3 for ALP and Fig. 4 for dark photon, we are ready to explore the allowed region for explaining the recently observed muon anomalous magnetic moment at Fermilab. This can be quantitatively done with χ 2 function, with central value ∆a µ and uncertainty σ(∆a µ ) taken from (5). The new physics prediction ∆a NP µ here contains the four parameters, m a and y µ a for ALP as well as m γ and for dark photon, in addition to the fixed coupling C aγγ = 3 TeV −1 and cut-off Λ = 1 TeV. To illustrate the allowed parameter space of ALP, for each point of Fig. 3 the values of m a and y µ a are fixed while the dark photon parameters m γ and are varied to ob-tain the smallest value χ 2 min (m a , y µ a ). The dark photon parameters scan in the range m γ ∈ [10 −3 , 350] GeV and ∈ [10 −5 , 1]. However, those points that fall inside the experimentally excluded regions of Fig. 4 are not included in the scan. The resulting χ 2 min (m a , y µ a ) is then a marginalized χ 2 function of just the two ALP parameters. Similar procedures can also produce a marginalized χ 2 min (m γ , ) after scanning the ALP parameters in the range of m a ∈ [10 −2 , 350] GeV and y µ a ∈ [10 −6 , 1] but deducting experimentally excluded regions.
So from Fig. 3 we can read off the values of m a and y µ a , but not the corresponding values of m γ and . The similar situation happens for Fig. 4. The black contours of m a and y µ a in Fig. 3 are obtained with χ 2 min (m a , y µ a ) < 5.99 at 95% C.L. and similarly for Fig. 4.
The black contours in Fig. 3 cover a large part of the remaining parameter space. As argued at the beginning of this section, the ALP mass is bounded from below, m a y µ a / 10 GeV, in order for the dark axion portal contribution to dominate. With smaller Yukawa coupling y µ a , the ALP mass can also be smaller and hence cover the whole mass range in Fig. 3. The ALP contribution is always negative [34] no matter what is the sign of the Yuakwa coupling y µ a as indicated by (9). This forbids the possibility of using only ALP to explain the observed positive ∆a µ . However, in the presence of dark axion portal, the ALP or pseudoscalar at large receives significant parameter space to explain the muon anomalous magnetic moment. Although the black contour is marginalized over the dark photon parameters, namely its mass m γ and kinetic mixing , we can still show the dependence on the dark photon mass by specifying the mass range of dark photon to be marginalized over. The solid black contour is obtained with m γ ≤ 200 GeV while the dashed one with 200 GeV ≤ m γ ≤ 350 GeV. It is interesting to see that with larger dark photon mass, the allowed ALP parameter space becomes larger with the Yukawa coupling y µ a touching down to as small as 10 −4 . The ALP solution can be readily saved by the dark axion portal.
For the dark photon parameter space illustrated in Fig. 4, almost all mass range below 200 GeV has been experimentally constrained to 10 −3 . Especially, the required parameter space, the hashed region in Fig. 4, for dark photon to explain the muon anomalous magnetic moment has been excluded by various observations including electron (g − 2) e , NA48/2, BaBar, and LHCb+CMS. It is very interesting to see that the dark axion portal coupling can also help to save the situation. Now the required parameter space significantly expands to the black contours, the solid one for m a ≤ 15 GeV and the dashed one for 15 GeV ≤ m a ≤ 350 GeV. The heavy mass region, m γ O(10) GeV, still has sizable space. Even the low mass region around m γ ≈ 10 MeV opens for 15 GeV ≤ m a ≤ 350 GeV.
In the large mass limit (m a , m γ m µ ), the total contribution shows decoupling features. To make it explicit, we fix the ALP parameters, m a = 110 GeV and y µ a = 0.082, as an example. The total contribution ∆a µ (m γ , ) is then a function of the two dark photon parameters m γ and . Fig. 5 shows the allowed region of ∆a µ (m γ , ) = (251 ± 59) × 10 −11 as grey area.
With larger dark photon mass, the dark photon coupling also needs to increase to maintain the prediction of ∆a µ (m γ , ). Otherwise, for a fixed , the predicted ∆a µ (m γ , ) would decrease with the dark photon mass. This decoupling behavior can be understood analytically with the approximate forms of (8), (9) and (10) in the large mass limit, For m γ m a , we can see that a µ ∝ ln Λ m γ and a γ µ ∝ 1/m 2 γ . Both terms decrease with m γ and hence have decoupling behavior. Fig. 6 shows the decoupling behavior in the ALP parameter space. Taking the two dark photon parameters, m γ = 100 GeV and = 0.003, the allowed region of ∆a µ (m a , y µ a ) = (251 ± 59) × 10 −11 in the m a − y µ a plane is the grey area. Since y µ a ∼ O(10 −1 ) is relatively large, both a a µ and a µ provide the dominant contributions, one has quadratic dependence on y µ a and the other linear. Consequently, for a given ALP mass m a , there are two possible solutions for the Yukawa coupling y µ a to match the experimental result. This is why the grey area of Fig. 6 is a circle rather than a line as in Fig. 5.
For m a m γ , we can also see decoupling features: with larger ALP mass the prediction ∆a µ (m a , y µ a ) be- comes smaller. This is because both a µ ∝ ln Λ ma and a a µ ≈ (m 2 µ /m 2 a ) ln(m a /m µ ) decrease with m a . With larger ALP mass, the Yukawa coupling y µ a should also increase in order to maintain the same prediction of ∆a µ (m a , y µ a ). However, with both dark photon and ALP parameters, the decoupling features are not transparent. For illustration, we take m γ = 100 GeV, m a = 110 GeV, and = 3 × 10 −3 (combined with our assumptions: C aγγ = 3 TeV −1 and Λ = 1 TeV), leading to, a µ ≈ 4.2 × 10 −8 y µ a , a a µ ≈ −1.4 × 10 −7 (y µ a ) 2 , (14) together with a γ µ ∆a µ . Then either y µ a = 0.082 or y µ a = 0.22 gives the observed discrepancy. The points (m a , y µ a ) = (110 GeV, 0.082 or 0.22) and (m γ , ) = (100 GeV, 3 × 10 −3 ) are inside our 95% C. L. band in Fig. 3 and Fig. 4, as expected, but y µ a and need not to be of the same order. In other words, the presence of both contributions can in fact enlarge the parameter space, including larger values of m a and m γ .
Conclusion -The latest measurement of the muon anomalous magnetic moment at the Fermilab Muon g-2 experiment further enhances the discrepancy with theoretical prediction from 3.7 σ to 4.2 σ. This clearly indicates that there is something new beyond the SM, although a decisive conclusion still awaits more data. On one hand, this discrepancy enhances theoretical exploration of possible solutions. But on the other, some solutions have been already excluded, including either the ALP or dark photon scenario. Even though the dark axion portal was originally motivated as a way to connect the visible world with the dark side, it can surprisingly save the ALP and dark photon for explaining the muon anomalous magnetic moment. Since dark matter contributes five times more energy density in the Universe than the ordinary matter and the latter already has rich particle spectrum, there is no reason to assume that the dark sector is composed of a single particle. In this sense, the dark axion portal provides a more interesting option than just ALP or dark photon. And muon anomalous magnetic moment can provide a robust probe of this new scenario.