Research Article | Open Access

Fayez Abu-Ajamieh, "Probing Scalar and Pseudoscalar Solutions to the Anomaly", *Advances in High Energy Physics*, vol. 2020, Article ID 1751534, 13 pages, 2020. https://doi.org/10.1155/2020/1751534

# Probing Scalar and Pseudoscalar Solutions to the Anomaly

**Academic Editor:**Antonio J. Accioly

#### Abstract

I investigate a class of models with scalar and pseudoscalar solutions to the anomaly for both the muon and the electron over the mass range of perturbativity ( GeV), with Yukawa couplings proportional to the lepton’s mass. In particular, I investigate the constraints from BaBar, beam dump experiments, decay measured quantities, LEP mono- searches, searches, and solar and horizontal branch (HB) star bounds. For a pseudoscalar, I find that no region in the parameter space can simultaneously provide a solution for both the electron and the muon anomalies while maintaining the required form of the couplings, and therefore, the pseudoscalar solution is disfavored. On the other hand, I find for the scalar case that there is an open window above ~30 MeV in the allowed region, but with significant tension with experiment for the region . In addition, there is a smaller window between ~350 KeV and 1 MeV that is not ruled out by cosmological observations. Part of the first open window is expected to be covered by the proposed NA64 experiment. Similar analysis can be readily applied to other proposed solutions to the anomaly, such as solutions with or with the dark photon.

#### 1. Introduction

An exciting piece of evidence for the existence of physics beyond the standard model (BSM) is the discrepancy between the predicted and the measured values of the muon anomalous magnetic moment . The current measured value [1–4] shows a discrepancy compared with the SM prediction [5–7]:

A similar less significant discrepancy of about was also observed for the electron [8]:

Although both discrepancies fall short of the limit required to confirm their existence, they nonetheless pose tantalizing hints for physics BSM. In addition, current experiments at Fermilab [9, 10] and at the J-PARC E34 collaboration [11, 12] are expected to yield improved experimental results in the near future.

New physics explanations of this anomaly include (see [13] for a comprehensive review) supersymmetry (see [14] for a review), a light boson [15–23] (also see [24] for a review), a scalar contribution within the framework of the 2 Higgs doublet model (2HDM) [25–32], additional fermions [33], leptoquarks [34, 35], and the dark photon [36].

Recently, there have been proposed solutions to this anomaly through a scalar [5] or a pseudoscalar axion-like particle (ALP) [6] in a general framework. In this short paper, I will investigate the viability of these solutions, explore the relevant experimental limits, and highlight the experimental probes for their discovery for the mass range of their validity.

For the case of a pseudoscalar, the effective interaction with photons and fermions can be parameterized by (This interaction can be viewed as an effective theory of a UV-complete model. One possible UV completion that is consistent with the SM EW theory was introduced in [5]. Such a model could lead to lepton-flavor violation through terms like . We ignore this possibility in this paper as it will not affect the interaction in Equation (3) and as it was studied in detail in [5]. The interested reader is instructed to refer to [5] for detailed analysis.) where is dimensionful coupling, is dimensionless Yukawa coupling, and , are the magnetic field strength tensor and its dual, respectively. For a scalar, is replaced with and there is no in the second term. Since (see Equations (5) and (6) below), and we can see from Equations (1) and (2) that within the allowed range of uncertainties, we could obtain ; (Notice that from Equations (5) and (6), , where . The assumption in Equation (4) makes , and thus .) then, we are motivated to define the Yukawa couplings to be proportional to the lepton mass: where is some model-dependent energy scale that is universal for all leptons, such as the axion decay constant or the radion constant. The coupling of the form given in Equation (4) has the additional advantage in that it arises in many UV completions ok pseudo Nambu-Goldstone bosons (PNGBs), such as axion models [37], left-right twin Higgs models [38], and dark matter models with a scalar portal to the dark sector [39] (also see [40]). PNGBs have an approximate shift symmetry, and their interactions are proportional to some universal symmetry breaking scale. In addition, such models have been widely discussed in literature as an effective theory for solving the anomaly (see for example [5, 6, 41]). I will focus on this form of Yukawa couplings throughout this paper.

It was shown in [5] that the discrepancy in can be explained by a scalar with , while in [6], it was shown that an ALP pseudoscalar can explain both of the electron and the muon anomalies by considering the NLO contributions.

The LO and NLO contributions to the are shown in Figure 1. The LO contribution for the scalar, as well as for the pseudoscalar, was calculated in [5]: where . On the other hand, the NLO contribution includes the Barr-Zee (BZ) contribution (top right diagram in Figure 1), the two-loop light-by-light (LBL) contribution (bottom left diagram in Figure 1), and the vacuum Polarization (VP) contribution (bottom right diagram in Figure 1). These contributions are the same for both the scalar and the pseudoscalar cases and are given by [6] where is either a scalar or a pseudoscalar, is the dimensionful coupling of to photons, and is some UV cutoff scale that is assumed to be much larger than . I will set the cutoff scale TeV throughout this paper. Notice that since the Lagrangian in Equation (3) is CP conserving, there will be no contribution to the lepton’s electric dipole moment (EDM). I will ignore the more general case where CP-violating terms are present.

#### 2. Favored Region

In this section, I will investigate the parameter space and try to establish the favored region for both the scalar and the pseudoscalar cases. Inspecting Equations (5)–(9), we notice the following: (i)For a scalar, the LO contribution is always positive, while for a pseudoscalar it is negative. As for the NLO contributions, we can see that the LBL and the VP are always positive, while the BZ contribution depends on the sign of (ii)Since the central measured anomaly for the muon is positive, while for the electron it is negative, a scalar solution can easily accommodate the muon anomaly. However, yielding the central measured electron anomaly would require the assumption that so that the BZ contribution can offset all other (positive) contributions. This would require large (nonperturbative) Yukawa couplings for the electron, and therefore, it is disfavored. Thus, I will assume that all couplings are positive for the scalar case. This means that the central measured electron anomaly cannot be produced. However, it is possible to show that the contribution is within (See added note at the end of this paper.)(iii)As the LO contribution for the pseudoscalar case is negative, it can easily accommodate the measured electron anomaly; however, in order to yield the (positive) measured muon anomaly, one needs somewhat large couplings to photons while keeping the Yukawa coupling somewhat small in order for the NLO contributions to dominate over the LO. Nonetheless, it is possible to find such solutions while maintaining perturbative couplings as we shall see below(iv)For the pseudoscalar case, if we assume that , then while could be either positive or negative. On the other hand, assuming that yields exactly the same solution but with opposite signs.

##### 2.1. Pseudoscalar

For concreteness, I will assume that . Figure 2 shows the allowed regions for and with . Notice that there is no overlap between the two regions even for small masses, which means that there is no region in the parameter space where Equation (4) is valid. As a matter of fact, there is no pseudoscalar solution in the whole parameter space where Equation (4) is true. In addition, for both Yukawa couplings to have the same sign, one needs , which is excluded by cosmological observations [42]; therefore, the pseudoscalar solution is disfavored. One can avoid these constraints by assuming that and ; however, one needs to justify this assumption. I will disregard the pseudoscalar solution in the remainder of this paper.

##### 2.2. Scalar

With the assumption that all couplings are positive, I attempt at finding the favored region in the parameter space for the scalar case. Here, I will only focus on solutions of the form in Equation (4) and discard all other possibilities.

Notice that we have three parameters, namely , , and . I will fix and keep and as free parameters. I will select appropriate benchmark points for by minimizing the of the electron and muon measurements:

More concretely, I fix the value of , then I use Equation (10) to find and that minimize , together with the value of at the minimum. Then, I will scan through a wide range of , and then set the benchmark points where is the smallest.

Figure 3 shows for several values of . As the plot shows, yields the lowest values of . Notice that becomes almost constant for smaller couplings. This is reasonable as when the coupling to photons becomes very small, the NLO contributions become negligible and the LO contribution is dominant. This high-level analysis seems to favor smaller couplings to photons, suggesting that the coupling to leptons is the dominant coupling. This is consistent with the cosmological constraints on ALPs (see for instance [42]). I will focus on this scenario and I will choose and as two benchmark points. Notice that to a good level of accuracy, the second benchmark point is representative of the entire region of , with very similar favored regions for the predicted mass and coupling.

Now, we can use Equations (6)–(9) in order to find the allowed region in the parameter space corresponding to a deviation from the central values. In order to set an upper limit on , we demand that all Yukawa couplings remain perturbative. Obviously, the most stringent bound comes from as it has the largest value. Requiring that , we obtain an upper bound on of GeV for . Figure 4 shows the allowed region corresponding to the two benchmark points. The plots show the bands for the Yukawa couplings to electrons and muons assuming that Equation (4) holds. In addition, the plots also show the region where the contribution to is within of the measured value assuming that the Yukawa couplings to leptons are independent of one another. Notice that the brown region corresponds to since we are assuming positive couplings. Thus, as noted earlier, there is no point in the parameter space that can yield the central value of the measured .

We can carry the analysis further to find the favored region in the parameters space for the benchmark points. Figure 5 shows the scalar mass and coupling to leptons that minimize , together with the and confidence level contours. Notice that this region is a subset of the allowed region in Figure 4. For , we find a predicted scalar mass of , with . For this point, one finds a predicted anomaly . On the other hand, for the second benchmark point (), the predicted mass and coupling are MeV and , respectively, with . This is consistent with the results found in [5].

Notice that the predicted values of the electron anomaly for these benchmark points are much smaller that the (absolute) measured central value of (see Equation (2)), albeit they are still within the limit. This is logical as we can see from Equation (6), if the Yukawa coupling has the form given in Equation (4). This is an important prediction to test this model. That is, we claim that if indeed this model is correct, then more accurate measurements of electron anomaly should yield . If future measurements of the electron anomaly are inconsistent with this, then the assumption in Equation (4) would be ruled out and other explanations would be needed. Another important prediction of this model is the tau anomaly. Given that Equation (4) predicts the anomaly to be proportional to the mass, then should be large enough to be measured. Specifically, we predict for .

#### 3. Experimental Probes and Limits

For the mass range 50 GeV, the most relevant constraints come from the BaBar experiment, beam dump experiments, the decay, the LEP mono- searches, and searches, and from the bounds on solar emission and the on emission from HB stars. Although there is an overlap between the results of this chapter and [43], where they discuss a similar effective model and a UV completion through the 2HDM, there are several novel features in this work, including bounds from the decay, LEP searches, and the solar and HB stars. In addition, here, we attempt at explaining both the electron and the muon anomalies and we extend the range of the mass over the entire range of validity. Thus, our work should be viewed as complementary to theirs.

##### 3.1. BaBar

Recent results from the BaBar experiment [44] searching for the process can be important for constraining the parameters space. The results can be used to extract the constraints on the process . ([45] extracts the constraints from BaBar’s results for the pseudoscalar case.) The tree-level Feynman diagrams that contribute to this process are shown in Figure 6, where I have neglected the diagrams where the scalar is radiated by the initial state particles or by the intermediate particle, since the coupling to muons dominates over the coupling to electrons and photons, and I am assuming a subleading coupling to the . Figure 4 in [44] sets an upper limit on the process . Therefore, we can extract the limits on the parameter space by requiring:

The excluded part of the parameter space is shown in the light gray region in Figure 7. Notice that it does not constrain the favored region corresponding to a deviation, although larger deviations would be constrained. This result also shows us that any new searches from BaBar (or similar experiments like Belle II) need to be at least an order of magnitude better in order to explore the favored region in the parameter space, or that searches should be made for a similar process with taus instead of muons, as the former has a larger coupling. Thus, there is a good motivation to search for processes like .

Notice that in calculating the branching fraction, we are assuming that the scalar can only decay to leptons or photons. If the scalar is allowed decay to other SM particles, the limits will become even weaker. Also notice that the upper limit in BaBar’s results found in Figure 4 of [44] is spiky, therefore, the upper limit used in Equation (11) contains significant uncertainty of about a factor of a few. In extracting , I was conservative and used the smallest cross-section. A less conservative estimate would relax the limits further.

##### 3.2. Decay ([46, 47] Also Discuss the 1-Loop Corrections to the Decay in the Context of R-Parity Violating Extensions to SUSY and in the Context of Type-II 2HDM)

The excellent measurements of the decay width and branching fractions present us with a potentially suitable tool for probing the parameter space. In particular, we can explore the limits associated with the scalar loop correction to the decay to a pair of leptons.

The scalar loop corrections can significantly affect the leptonic decay width of the boson. The NLO corrections to are shown in Figure 8, where the coupling of to is assumed to be subdominant compared with the coupling to leptons. Notice here that UV divergences in the leg corrections cancel that in the vertex correction and that for a massive the result is free of IR divergences.

Dropping the lepton mass in the loops, and keeping only as an IR regulator, which is justified for where this bound is moslty relevant, the NLO correction is approximately given by where is the LO decay width given by:

We can compare this correction with the branching fractions of the leptonic decays, which are given by [48]:

As it turns out, the decay to provides the most stringent constraints in spite of the larger uncertainty in its branching fraction. This is due to its larger coupling given the assumption in Equation (4).

The excluded region of the parameter space at a level is shown in red in Figure 7. As can be seen from the plots, there is some tension between the loop decay and the allowed region near GeV, although the alowed region is not fully excluded. Notice here that since we are dropping the mass of the tau, the bound will be less reliable for the range ; however, this is unimportant as in that region the decay constraint is far from the favored region and therefore our results are unaffected. On the other hand, in the region of interest where is larger than a few GeV, dropping the mass of the tau is justified and will not impact the results significantly. The same argument applies for the constraint obtained from we discuss below.

##### 3.3. Constraints from

The analysis conducted in the previous section can be extended to the searches conducted by the LEP experiment [49–64], the KEK collaboration [65–68], DESY-PETRA collaboration [69–74], and the SLAC-PEP experiment [75].

Here, the process of proceeds through the s-channel with a photon or a propagator, and the NLO loop correction due to will be identical to the case of the decay shown in Figure 8. So we can set the bound as where and is the tree level cross-section. The bound is shown in Figure 7 in light green. As the plot shows, and similar to the case of the decay, there is a significant tension with experimental results for GeV, although this region is not fully excluded.

##### 3.4. Muon Beam Dump Experiments

Muon beam dump experiments provide a powerful tool for probing the mass range MeV. The relevant constraints come from Orsay [76] and the E137 experiment at SLAC [77] (also see [41] for a summary).

In the Orsay beam dump experiment, searches for the light Higgs boson in the 2HDM were conducted through looking for the process . In [76], the coupling of the lighter Higgs to electrons is assumed to be where , the ratio of the two doublets’ VEVs. The results show the excluded region in the parameter space. Therefore, they can readily be extrapolated to this model by setting .

On the other hand, [77] presents the results of the E137 beam dump searches for axions produced via bremsstrahlung followed by the subsequent decay to . The results were extracted for the case of a scalar in [41], so I will just use their results.

I show these constraints in magenta (Orsay) and orange (E137) in Figure 7. As can be seen from the plots, muon beam dump experiments exclude the region between MeV and 30 MeV. On the other hand, the window between MeV is still open. This window is projected to be explored by the proposed NA64 project at CERN [78, 79]. The NA64 experiment is a fixed-target experiment that can run in the muon mode with a beam energy of 160 GeV and is designed for searching for missing energy GeV. This experiment can help probe this open window. The projected region in the parameter space is shown by the dashed line in Figure 7.

Another proposed experiment is Fermilab’s displaced decay search with a muon beam energy of 3 GeV [80]. However, the projected sensitivity of this experiment covers only a part of the projected sensitivity of the NA64 experiment, therefore I will not plot it here.

##### 3.5. Constraints from LEP Mono- and Searches

The LEP mono- searches were conducted to set limits on the number of neutrinos via studying the process . The results can be used to set limits on the parameter space by considering the process for the range . The tree-level process proceeds through the and channels; however, due to the smallness of the scalar’s coupling to the electron, the constraints are weak. On the other hand, stronger constraints can be obtained through the triangle diagram shown in Figure 9 where the runs in the loop.

Similar to the case of the Higgs, we can write the effective Lagrangian as where . The effective couplings can be extracted from the results of the Higgs. For example, we can use the result of the Higgs decay to where instead of the top running in the loop, we have the tau. This gives where and

Similarly, we find from the Higgs decay to . We obtain: where and

Here, , and . Armed with this, we can find the cross-section of the triangle diagram in Figure 9: where are the sine and cosine of the Weinberg angel, respectively. The results from ALEPH [81, 82], L3 [83–85], OPAL [86], and DELPHI collaborations [87] can be used in order set constraints on the scalar by requiring that at a level: for . Notice that for , could decay to ; however, for both benchmark points, the decay length is orders of magnitude larger than the dimensions of the detector, so appears as missing energy. On the other hand, for , where , , ,; can decay to a lepton pair, so we can use LEP searches for that region. The cross-section for can be readily obtained by multiplying the cross-section in Equation (29) by the appropriate branching fraction for the mass range. The constraints are shown in brown in Figure 7.

##### 3.6. Constraints from Solar Emission, HB Stars, and SN1987A

Supernova 1987 (SN1987A), HB stars, and solar emission constraints can impose stringent constraints on ALPs for masses (see [42, 88] for instance). However, such constraints are only relevant if the ALP’s dominant coupling is to photons. Since the coupling to photons in the type of models we are considering in this paper is favored to be much less than that to leptons (at least for the scalar case), those limits need to be revisited.

If is produced in the sun or in HB stars, then it could affect the measured energy loss rate of the star when it streams out, in addition to affecting the star’s evolution. Therefore one can obtain a bound on and by requiring that the amount of energy carried away by the scalar be less than the observed limits.

In this model where the scalar’s couplings to leptons are dominant, will be mainly produced in stars via its interaction with electrons. Assuming that the electrons are nonrelativistic and nondegenerate, which is a good assumption for the solar and HB stars’ mediums, the Compton-like scattering dominates over both Bremsstrahlung and electron-positron annihilation. For , we can neglect the recoil energy of the electron and to a good approximation the energy loss rate per unit volume is given by [89]: where is the energy of the photon and is the number density of the electrons, which in terms of the electron fraction in the star , the mass density of the star , and the atomic mass unit , is given by

This calculation was done by Grifols and Massó in [90] with neglected. Here, I will keep the mass of the scalar but assume it is less than and keep the assumption that the recoil energy of the electron is small. The cross-section is given by where , , and is some complicated function that I relegate to the appendix. The bound on the solar (HB) emission rate is [89]

Given this bound, we can use Equation (33) to solve Equation (31) numerically in order to find the excluded region in the parameter space. Here, I set the average temperature to be for the solar (HB) medium, the average density to be for the solar (HB) medium and use for both. The yellow region in Figure 7 shows the excluded part of the parameter space by the solar constraints, while the dark green region shows the excluded part by HB stars. As the plot shows, all masses below are excluded. However, we must note that for close to the electron mass, the bound is less rigorous as our assumption of a small electron’s recoil energy becomes less valid. However, I checked numerically that including the recoild energy of the final state enelctron does not significantly impact the bound. Also, notice that the HB constraints are much more stringent than the solar ones in spite of the weaker bound in Equation (34). This is due to the higher temperature of HB stars compared with the sun, which yields a larger Boltzmann factor in Equation (31).

Of course, this analysis is valid only if streams freely out of the star, i.e., if it does not get trapped inside the medium of the star. In order to verify the validity of this assumption, we can calculate the mean free path of and compare it with the radius of the star. The mean free path is given by

Assuming that the average photons energy is given by one can easily check that for the entire mass range of interest, the mean free path is orders of magnitude larger than the radius of the sun or the typical radius of an HB star, thereby justifying the free streaming assumption.

A similar argument applies for scalars produced in supernovas. However, it was shown in [40] that for a scalar produced in the supernova core, the mean free path is given by which means that for larger than GeV, the scalar gets trapped in the core and never streams out. Thus, SN1987A does not costrain the favored region and therefore we ignore it.

##### 3.7. Discussion

We have shown that the constraints from BaBar, beam dump experiments, the NLO correction to the decay, LEP mono- searches, searches, and solar and HB star constraint, exclude a significant part of the parameter space. Figure 7 shows that all masses below are excluded by solar and HB constraints; beam dump experiments exclude the region above 1 MeV up to , while masses above are in significant tension with the NLO corrections to the decay to and the NLO corrections to , although this region is not entirely excluded. This leaves a large open window between and that is most favored to be explored.

Part of this region is projected to be explored by the proposed NA64 experiment, which is projected to cover the mass range from 1 MeV and up to the dimuon mass. The first benchmark point corresponding to lies in this region. On the other hand, the second benchmark point lies outside this region but could be explored by the Belle II experiment [91]. The Belle II experiment has recently started collecting data and it is expected that by 2025, it would have reached a total integrated luminosity of . This projects it to be more sensitive than the BaBar experiment and therefore might help explore more of the parameter space above the dimuon threshold which is not covered by NA64 and up to . This leaves two windows, one above and the other is between and 1 MeV.

The international linear collider (ILC) [92] (if built) might help explore both regions. For the former, the ILC’s ultraprecision measurement of the production and subsequent decay can help improve the bounds from the decay by lowering the measured uncertainties, while for the latter region, the channel can help investigate the hypothetical process .

A final point to mention is that due to the assumption made in Equation (4), the dominant coupling will be to the tau lepton, and therefore collider and beam dump experiments limits are more stringent. If we assume a different type of coupling to leptons, say by assuming a suppressed coupling to the tau compared to the muon, then the constraints will be alleviated, and more of the parameter space will open. We will very briefly discuss one model where this can be achieved in the next section.

#### 4. A Radion Solution for the Anomaly?

In the Randall-Sundrum (RS) model [93], the radion is the scalar field that parameterizes the fluctuations of the extra dimension around its potential minimum. The radion could pose an interesting possibility for solving the anomaly due to its unique couplings to matter. More specifically, the radion’s coupling to matter is highly model-dependent and varies according to the localization of the matter fields on either of the branes or in the bulk. The coupling to brane-localized matter is given by [94]: where is the 4D radion field, is the radion constant on the UV and IR branes, respectively, and is the stress-energy tensor. Since the UV scale is typically many orders of magnitude larger than the IR scale, it is possible to suppress the coupling to the tau lepton compared to the muon by assuming that the former is localized on the UV brane, while assuming that the latter is localized on the IR brane. This way, one could alleviate all of the constraints (except for beam dump experiments) as the couplings will be rescaled by .

Although the typical mass of the radion is comparable to the electroweak (EW) scale (), much lighter masses can be achieved through the Contino-Pomarol-Rattzzi (CPR) mechanism [95] as was demonstrated in [96, 97].

Another interesting aspect of a radion solution is that the radion could couple to nucleons and pions through quarks and gluons [96], which presents additional experimental probes. Focusing on the coupling to pions, we can write the effective Lagrangian as where is the effective radion coupling to pions with dimension . If the radion is heavy enough, it could decay to :

The decay width in Equation (40) could be small for typical values of but could still be measurable. For instance, the decay width for a 400 MeV radion with GeV would be KeV.

For lighter masses, searches for the rare pion decay could provide an interesting search option. If we assume that dominates over , then the branching fraction of this hypothetical decay would be given by

For instance, GeV would yield a branching fraction comparable to the observed rare decay .

#### 5. Conclusions

The anomaly remains one of the best ways to search for physics BSM. In this paper, I investigated a class of models with a scalar/pseudoscalar that has a coupling to leptons proportional to the lepton’s mass.

We saw in this paper that for the case of a pseudoscalar solution, there is no region in the parameter space that could simultaneously solve both the electron and the muon anomalies with Yukawa couplings of the form in Equation (4). However, if no assumption is made regarding the form of Yukawa’s, it is possible to have a pseudoscalar solution for both of the anomalies. Nevertheless, this solution is not very attractive since it would require tuning the scalar’s coupling to photons to be somewhat large in order for both of the Yukawa couplings to have the same sign. Such a large coupling to photons would be disfavored by cosmological observations.

On the other hand, a scalar can simultaneously provide a solution for both anomalies while having the required form of the Yukawa couplings. In such a case, we demonstrated that such a solution favors smaller couplings to photons GeV, and we established the corresponding favored region in the parameter space that corresponds to two representative benchmark points. We investigated the experimental constraints from the BaBar experiment, beam dump experiments, the NLO corrections to the decay and to , LEP mono-, and searches, and from solar and HB emission bounds, and we saw that a significant part of the parameter space is excluded. In particular, we found only two open windows for between and 1 MeV, and between and , with the region above being in significant tension with experiment.

The tools used in this paper can be used to constrain other solutions to the anomaly, such as solutions that adopt the or the dark photon to explain it. I expect that the limits on these solutions would not be too different from the scalar case for the same range of masses and couplings, they are nonetheless worthwhile investigating.

Future experiments, such as NA64, Belle II, and the ILC can help explore significant regions of the parameter space, and one hopes that in the near future, enough data would be collected to shed more light on the remaining open windows, and thus help explore the viability of this solution.

#### Appendix

#### A. Explicit Expression for

Added note: After finalizing this paper, a new result [98] was published claiming a discrepancy in the measured electron anomaly:

If this yet uncorroborated result is true, then our results cannot explain both anomalies simultaneously at the level, although the results remain valid at the level. This is because the electron anomaly will be negative at the level and a novel explanation would be needed to explain the sign discrepancy between the two anomalies. In this case, the electron anomaly in this work would serve more as a constraint on the allowed region of muon anomaly, and the result in Figure 7 would be relevant for the muon anomaly only. Either way, a corroborating result is still needed to confirm the value and the sign of the electron anomaly.

#### Data Availability

All data are available online.

#### Conflicts of Interest

The author declares that he/she have no conflicts of interest.

#### Acknowledgments

I would like to thank John Terning for his valuable insight. I would also like to thank John Conway, Lloyd Knox, and Max Chertok for answering my questions.

#### References

- P. J. Mohr, B. N. Taylor, and D. B. Newell, “CODATA recommended values of the fundamental physical constants: 2010,”
*Reviews of Modern Physics*, vol. 84, no. 4, pp. 1527–1605, 2012. View at: Publisher Site | Google Scholar - Muon (g - 2) Collaboration, “Measurement of the Positive Muon Anomalous Magnetic Moment to 0.7 ppm,”
*Physical Review Letters*, vol. 89, no. 10, article 101804, 2002. View at: Publisher Site | Google Scholar - Muon (g - 2) Collaboration), “Measurement of the Negative Muon Anomalous Magnetic Moment to 0.7 ppm,”
*Physical Review Letters*, vol. 92, no. 16, article 161802, 2004. View at: Publisher Site | Google Scholar - Muon g-2 Collaboration, “Final report of the E821 muon anomalous magnetic moment measurement at BNL,”
*Physical Review D*, vol. 73, no. 7, article 072003, 2006. View at: Publisher Site | Google Scholar - C. Y. Chen, H. Davoudiasl, W. J. Marciano, and C. Zhang, “Implications of a light “dark Higgs” solution to theg
*μ*−2discrepancy,”*Physical Review D*, vol. 93, no. 3, article 035006, 2016. View at: Publisher Site | Google Scholar - W. J. Marciano, A. Masiero, P. Paradisi, and M. Passera, “Contributions of axionlike particles to lepton dipole moments,”
*Physical Review D*, vol. 94, no. 11, article 115033, 2016. View at: Publisher Site | Google Scholar - P. deNiverville, H. S. Lee, and M. S. Seo, “Implications of the dark axion portal for the muon
*g*−2,*B*factories, fixed target neutrino experiments, and beam dumps,”*Physical Review D*, vol. 98, no. 11, article 115011, 2018. View at: Publisher Site | Google Scholar - T. Aoyama, M. Hayakawa, T. Kinoshita, and M. Nio, “Tenth-order electron anomalous magnetic moment: contribution of diagrams without closed lepton loops,”
*Physical Review D*, vol. 91, no. 3, article 033006, 2015. View at: Publisher Site | Google Scholar - J. Grange, V. Guarino, P. Winter et al., “Muon (g-2) Technical Design Report,” http://arxiv.org/abs/1501.06858. View at: Google Scholar
- G. Venanzoni, “The New Muon g−2 experiment at Fermilab,”
*Nuclear and Particle Physics Proceedings*, vol. 273-275, pp. 584–588, 2016. View at: Publisher Site | Google Scholar - N. Saito and J-PARC g−2/EDM Collaboration, “A novel precision measurement of muon g - 2 and EDM at J-PARC,”
*AIP Conference Proceedings*, vol. 1467, p. 45, 2012. View at: Publisher Site | Google Scholar - T. Mibe, “Measurement of muon and EDM with an ultra-cold muon beam at J-PARC,”
*Nuclear Physics B - Proceedings Supplements*, vol. 218, no. 1, pp. 242–246, 2011. View at: Publisher Site | Google Scholar - M. Lindner, M. Platscher, and F. S. Queiroz, “A call for new physics: The muon anomalous magnetic moment and lepton flavor violation,”
*Physics Reports*, vol. 731, pp. 1–82, 2018. View at: Publisher Site | Google Scholar - D. Stockinger, “The muon magnetic moment and supersymmetry,”
*Journal of Physics G: Nuclear and Particle Physics*, vol. 34, no. 2, pp. R45–R91, 2007. View at: Publisher Site | Google Scholar - P. Langacker, “The physics of heavy gauge bosons,”
*Reviews of Modern Physics*, vol. 81, no. 3, pp. 1199–1228, 2009. View at: Publisher Site | Google Scholar - S. Baek, N. G. Deshpande, X. G. He, and P. Ko, “Muon anomalous g−2 and gauged L
*μ*−L*τ*models,”*Physical Review D*, vol. 64, no. 5, article 055006, 2001. View at: Publisher Site | Google Scholar - E. Ma, D. P. Roy, and S. Roy, “Gauged L
*μ*−L*τ*with large muon anomalous magnetic moment and the bimaximal mixing of neutrinos,”*Physics Letters B*, vol. 525, no. 1-2, pp. 101–106, 2002. View at: Publisher Site | Google Scholar - S. N. Gninenko and N. V. Krasnikov, “The muon anomalous magnetic moment and a new light gauge boson,”
*Physics Letters B*, vol. 513, no. 1-2, pp. 119–122, 2001. View at: Publisher Site | Google Scholar - M. Pospelov, “Secluded U(1) below the weak scale,”
*Physical Review D*, vol. 80, no. 9, article 095002, 2009. View at: Publisher Site | Google Scholar - J. Heeck and W. Rodejohann, “Gauged L
_{μ}−L_{τ}symmetry at the electroweak scale,”*Physical Review D*, vol. 84, no. 7, article 075007, 2011. View at: Publisher Site | Google Scholar - K. Harigaya, T. Igari, M. M. Nojiri, M. Takeuchi, and K. Tobe, “Muon g-2 and LHC phenomenology in the L
_{μ}− L_{τ}gauge symmetric model,”*Journal of High Energy Physics*, vol. 2014, no. 3, p. 105, 2014. View at: Publisher Site | Google Scholar - W. Altmannshofer, S. Gori, M. Pospelov, and I. Yavin, “Neutrino trident production: a powerful probe of new physics with neutrino beams,”
*Physical Review Letters*, vol. 113, no. 9, article 091801, 2014. View at: Publisher Site | Google Scholar - W. Altmannshofer, C. Y. Chen, P. S. Bhupal Dev, and A. Soni, “Lepton flavor violating explanation of the muon anomalous magnetic moment,”
*Physics Letters B*, vol. 762, pp. 389–398, 2016. View at: Publisher Site | Google Scholar - N. V. Krasnikov, “The muon (g-2) anomaly and a new light vector boson,” http://arxiv.org/abs/1702.04596. View at: Google Scholar
- E. O. Iltan and H. Sundu, “Probing Scalar and Pseudoscalar Solutions of the g-2 Anomaly,” http://arxiv.org/abs/hep-ph/010310. View at: Google Scholar
- Y. Omura, E. Senaha, and K. Tobe, “Lepton-flavor-violating Higgs decay h ⟶
*μτ*and muon anomalous magnetic moment in a general two Higgs doublet model,”*Journal of High Energy Physics*, vol. 2015, no. 5, p. 28, 2015. View at: Publisher Site | Google Scholar - A. Broggio, E. J. Chun, M. Passera, K. M. Patel, and S. K. Vempati, “Limiting two-Higgs-doublet models,”
*Journal of High Energy Physics*, vol. 2014, no. 11, p. 58, 2014. View at: Publisher Site | Google Scholar - L. Wang and X. F. Han, “A light pseudoscalar of 2HDM confronted with muon g-2 and experimental constraints,”
*Journal of High Energy Physics*, vol. 2015, no. 5, p. 39, 2015. View at: Publisher Site | Google Scholar - T. Abe, R. Sato, and K. Yagyu, “Lepton-specific two Higgs doublet model as a solution of muon g − 2 anomaly,”
*Journal of High Energy Physics*, vol. 2015, no. 7, p. 064, 2015. View at: Publisher Site | Google Scholar - J. Cao, P. Wan, L. Wu, and J. M. Yang, “Lepton-specific two-Higgs doublet model: experimental constraints and implication on Higgs phenomenology,”
*Physical Review D*, vol. 80, no. 7, article 071701, 2009. View at: Publisher Site | Google Scholar - B. Batell, M. Pospelov, and A. Ritz, “Multilepton signatures of a hidden sector in rareBdecays,”
*Physical Review D*, vol. 83, no. 5, article 054005, 2011. View at: Publisher Site | Google Scholar - S. Dutta, A. Goyal, and M. P. Singh, “Lepto-philic 2-HDM + singlet scalar portal induced fermionic dark matter,”
*Journal of High Energy Physics*, vol. 2019, no. 7, p. 76, 2019. View at: Publisher Site | Google Scholar - A. Freitas, J. Lykken, S. Kell, and S. Westhoff, “Testing the muon g-2 anomaly at the LHC,”
*Journal of High Energy Physics*, vol. 2014, no. 5, p. 145, 2014. View at: Publisher Site | Google Scholar - D. Chakraverty, D. Choudhury, and A. Datta, “A nonsupersymmetric resolution of the anomalous muon magnetic moment,”
*Physics Letters B*, vol. 506, no. 1-2, pp. 103–108, 2001. View at: Publisher Site | Google Scholar - K. Cheung, “Muon anomalous magnetic moment and leptoquark solutions,”
*Physical Review D*, vol. 64, no. 3, article 033001, 2001. View at: Publisher Site | Google Scholar - H. Davoudiasl, H. S. Lee, and W. J. Marciano, “Muong−2, rare kaon decays, and parity violation from dark bosons,”
*Physical Review D*, vol. 89, no. 9, article 095006, 2014. View at: Publisher Site | Google Scholar - C. Hagmann, H. Murayama, G. G. Raffelt, L. J. Rosenberg, and K. V. Bibber, 2008, http://pdg.lbl.gov/2008/reviews/rpp2008-rev-axions.pdf.
- G. L. Liu and Q. G. Zeng, “Muon
*g*-2 anomaly confronted with the higgs global data in the left-right twin Higgs models,”*European Physical Journal C: Particles and Fields*, vol. 79, no. 7, p. 612, 2019. View at: Publisher Site | Google Scholar - Y. Nomura and J. Thaler, “Dark matter through the axion portal,”
*Physical Review D*, vol. 79, no. 7, article 075008, 2009. View at: Publisher Site | Google Scholar - R. Essig, R. Harnik, J. Kaplan, and N. Toro, “Discovering new light states at neutrino experiments,”
*Physical Review D*, vol. 82, no. 11, article 113008, 2010. View at: Publisher Site | Google Scholar - C. Y. Chen, M. Pospelov, and Y. M. Zhong, “Muon beam experiments to probe the dark sector,”
*Physical Review D*, vol. 95, no. 11, article 115005, 2017. View at: Publisher Site | Google Scholar - J. Jaeckel and M. Spannowsky, “Probing MeV to 90 GeV axion-like particles with LEP and LHC,”
*Physics Letters B*, vol. 753, pp. 482–487, 2016. View at: Publisher Site | Google Scholar - B. Batell, N. Lange, D. McKeen, M. Pospelov, and A. Ritz, “Muon anomalous magnetic moment through the leptonic Higgs portal,”
*Physical Review D*, vol. 95, no. 7, article 075003, 2017. View at: Publisher Site | Google Scholar - BaBar Collaboration, “Search for a muonic dark force atBaBar,”
*Physical Review D*, vol. 94, no. 1, article 011102, 2016. View at: Publisher Site | Google Scholar - M. Bauer, M. Neubert, and A. Thamm, “Collider probes of axion-like particles,”
*Journal of High Energy Physics*, vol. 2017, no. 12, p. 44, 2017. View at: Publisher Site | Google Scholar - O. Lebedev, W. Loinaz, and T. Takeuchi, “Constraints on R-parity violating couplings from lepton universality,”
*Physical Review D*, vol. 61, no. 11, article 115005, 2000. View at: Publisher Site | Google Scholar - O. Lebedev, W. Loinaz, and T. Takeuchi, “Constraints on two-Higgs-doublet models at large tan
*β*from W and Z decays,”*Physical Review D*, vol. 62, no. 5, article 055014, 2000. View at: Publisher Site | Google Scholar - Particle Data Group, “Review of Particle Physics,”
*Chinese Physics C*, vol. 40, no. 10, article 100001, 2016. View at: Publisher Site | Google Scholar - The OPAL Collaboration, “Tests of the standard model and constraints on new physics from measurements of fermion-pair production at 189-209 GeV at LEP,”
*European Physical Journal C: Particles and Fields*, vol. 33, no. 2, pp. 173–212, 2004. View at: Publisher Site | Google Scholar - The OPAL Collaboration, “Search for Yukawa Production of a light neutral Higgs Boson at LEP,”
*European Physical Journal C: Particles and Fields*, vol. 23, no. 3, pp. 397–407, 2002. View at: Publisher Site | Google Scholar - The OPAL Collaboration, “Tests of the standard model and constraints on new physics from measurements of fermion-pair production at 189 GeV at LEP,”
*European Physical Journal C: Particles and Fields*, vol. 13, no. 4, pp. 553–572, 2000. View at: Publisher Site | Google Scholar - L3 Collaboration, “Measurement of cross sections and leptonic forward-backward asymmetries at theZ pole and determination of electroweak parameters,”
*Zeitschrift für Physik C Particles and Fields*, vol. 62, no. 4, pp. 551–573, 1994. View at: Publisher Site | Google Scholar - M. Z. Akrawy, G. Alexander, J. Allison et al., “Analysis of Z0 couplings to charged leptons,”
*Physics Letters B*, vol. 247, no. 2-3, pp. 458–472, 1990. View at: Publisher Site | Google Scholar - OPAL Collaboration, “Precision measurements of the neutral current from hadron and lepton production at LEP,”
*Zeitschrift für Physik C Particles and Fields*, vol. 58, no. 2, pp. 219–237, 1993. View at: Publisher Site | Google Scholar - M. Z. Akrawy, G. Alexander, J. Allison et al., “Measurements of the decay of the Z0 into lepton pairs,”
*Physics Letters B*, vol. 235, no. 3-4, pp. 379–388, 1990. View at: Publisher Site | Google Scholar - The DELPHI Collaboration, “Cross-sections and leptonic forward-backward asymmetries from the Z
^{0}running of LEP,”*European Physical Journal C: Particles and Fields*, vol. 16, no. 3, pp. 371–405, 2000. View at: Publisher Site | Google Scholar - M. Z. Akrawy, G. Alexander, J. Allison et al., “A combined analysis of the hadronic and leptonic decays of the Z0,”
*Physics Letters B*, vol. 240, no. 3-4, pp. 497–512, 1990. View at: Publisher Site | Google Scholar - P. Aarnio, P. Abreu, W. Adam et al., “Study of the leptonic decays of the Z0 boson,”
*Physics Letters B*, vol. 241, no. 3, pp. 425–434, 1990. View at: Publisher Site | Google Scholar - P. Abreu, W. Adam, F. Adami et al., “Determination of Z0 resonance parameters and couplings from its hadronic and leptonic decays,”
*Nuclear Physics B*, vol. 367, no. 3, pp. 511–574, 1991. View at: Publisher Site | Google Scholar - ALEPH Collaboration, “Update of electroweak parameters fromZ decays,”
*Zeitschrift für Physik C Particles and Fields*, vol. 60, no. 1, pp. 71–81, 1993. View at: Publisher Site | Google Scholar - ALEPH Collaboration, “Measurement of electroweak parameters fromZ decays into fermion pairs,”
*Zeitschrift fur Physik C: Particle Fields*, vol. 48, no. 3, pp. 365–391, 1990. View at: Publisher Site | Google Scholar - L3 Collaboration, “Search for a at the Z resonance,”
*Physics Letters B*, vol. 306, no. 1-2, pp. 187–196, 1993. View at: Publisher Site | Google Scholar - ALEPH Collaboration, “Improved measurements of electroweak parameters fromZ decays into fermion pairs,”
*Zeitschrift für Physik Physik C Particles and Fields*, vol. 53, no. 1, pp. 1–20, 1992. View at: Publisher Site | Google Scholar - ALEPH Collaboration, “Measurement of hadron and lepton-pair production from e
^{+}e^{−}annihilation at centre-of-mass energies of 130 and 136 GeV,”*Physics Letters B*, vol. 378, no. 1-4, pp. 373–384, 1996. View at: Publisher Site | Google Scholar - TOPAZ Collaboration, “Charge asymmetry measurement in e
^{+}e‑⟶*μ*^{+}*μ*^{−}and*τ*^{+}*τ*^{−}reactions at and 55 GeV,”*Physics Letters B*, vol. 208, no. 2, pp. 319–323, 1988. View at: Publisher Site | Google Scholar - MAC Collaboration, “A measurement of the forward-backward charge asymmetry at √s=29 GeV,”
*Physics Letters B*, vol. 218, no. 3, pp. 369–373, 1989. View at: Publisher Site | Google Scholar - AMY Collaboration, “Measurements of cross section and charge asymmetry for e
^{+}e^{−}⟶*μ*^{+}*μ*^{−}and e^{+}e^{−}⟶*τ*^{+}*τ*^{−}at = 57.8,”*Physics Letters B*, vol. 331, no. 1-2, pp. 227–235, 1994. View at: Publisher Site | Google Scholar - TOPAZ Collaboration, “Cross sections and charge asymmetries for e
^{+}e^{−}⟶*μ*^{+}*μ*^{−}and e^{+}e^{−}⟶*τ*^{+}*τ*^{‑}reactions at,”*Physics Letters B*, vol. 291, no. 1-2, pp. 206–212, 1992. View at: Publisher Site | Google Scholar - Cello Collaboration, “Measurement of e
^{+}e^{−}⟶*τ*^{+}*τ*^{−}at high energy and properties of the*τ*lepton,”*Physics Letters B*, vol. 114, no. 4, pp. 282–286, 1982. View at: Publisher Site | Google Scholar - TASSO Collaboration, “A measurement of electroweak effects in the reactione
^{+}e^{−}⟶*τ*^{+}*τ*^{−}at 35.0 and 42.4 GeV,”*Zeitschrift für Physik C Particles and Fields*, vol. 43, no. 4, pp. 549–555, 1989. View at: Publisher Site | Google Scholar - JADE Collaboration, “Tau-lepton production and decay at Petra energies,”
*Physics Letters B*, vol. 161, no. 1-3, pp. 188–196, 1985. View at: Publisher Site | Google Scholar - JADE Collaboration, “Final results on muon and tau pair production by the JADE Collaboration at PETRA,”
*Zeitschrift für Physik C Particles and Fields*, vol. 46, no. 4, pp. 547–553, 1990. View at: Publisher Site | Google Scholar - Cello Collaboration, “
*τ*Production and decay with the cello detector at petra,”*Physics Letters B*, vol. 222, no. 1, pp. 163–172, 1989. View at: Publisher Site | Google Scholar - JADE Collaboration, “Tests of the standard model in leptonic reactions at PETRA energies,”
*Zeitschrift für Physik C Particles and Fields*, vol. 30, no. 3, pp. 371–380, 1986. View at: Publisher Site | Google Scholar - K. K. Gan, I. Beltrami, B. G. Bylsma et al., “Measurement of the reaction e
^{+}e^{−}⟶*τ*^{+}*τ*^{−}at √s = 29 GeV,”*Physics Letters B*, vol. 153, no. 1-2, pp. 116–120, 1985. View at: Publisher Site | Google Scholar - M. Davier and H. Nguyen Ngoc, “An unambiguous search for a light Higgs boson,”
*Physics Letters B*, vol. 229, no. 1-2, pp. 150–155, 1989. View at: Publisher Site | Google Scholar - J. D. Bjorken, S. Ecklund, W. R. Nelson et al., “Search for neutral metastable penetrating particles produced in the SLAC beam dump,”
*Physical Review D*, vol. 38, no. 11, pp. 3375–3386, 1988. View at: Publisher Site | Google Scholar - S. Andreas, S. V. Donskov, P. Crivelli et al., “Proposal for an Experiment to Search for Light Dark Matter at the SPS,” http://arxiv.org/abs/1312.3309. View at: Google Scholar
- S. N. Gninenko, N. V. Krasnikov, and V. A. Matveev, “Muong−2 and searches for a new leptophobic sub-GeV dark boson in a missing-energy experiment at CERN,”
*Physical Review D*, vol. 91, no. 9, article 095015, 2015. View at: Publisher Site | Google Scholar - A. Chapelain and Muon g-2 Collaboration, “The Muon g-2 experiment at Fermilab,”
*EPJ Web of Conferences*, vol. 137, article 08001, 2017. View at: Publisher Site | Google Scholar - The ALEPH Collaboration, “Single- and multi-photon production in ee collisions at up to 209 GeV,”
*European Physical Journal C: Particles and Fields*, vol. 28, no. 1, pp. 1–13, 2003. View at: Publisher Site | Google Scholar - ALEPH Collaboration, “Single- and multi-photon production in e
^{+}e^{−}collisions at a centre-of- mass energy of 183 GeV,”*Physics Letters B*, vol. 429, no. 1-2, pp. 201–214, 1998. View at: Publisher Site | Google Scholar - L3 Collaboration, “Determination of the number of light neutrino species from single photon production at LEP,”
*Physics Letters B*, vol. 431, no. 1-2, pp. 199–208, 1998. View at: Publisher Site | Google Scholar - L3 Collaboration, “Single and multi-photon events with missing energy in e
^{+}e^{−}collisions at 161GeV<s<172Gev,”*Physics Letters B*, vol. 415, no. 3, pp. 299–310, 1997. View at: Publisher Site | Google Scholar - L3 Collaboration, “Single- and multi-photon events with missing energy in e
^{+}e^{−}collisions at LEP,”*Physics Letters B*, vol. 587, no. 1-2, pp. 16–32, 2004. View at: Publisher Site | Google Scholar - OPAL Collaboration, “Multiphoton production in e
^{+}e^{-}collisions at s =183-GeV,”*Physics Letters B*, vol. 438, p. 379, 1998. View at: Publisher Site | Google Scholar - The DELPHI Collaboration, “Photon events with missing energy in e
^{+}e^{-}collisions at =130-GeV to 209-GeV,”*European Physical Journal C: Particles and Fields*, vol. 38, p. 395, 2005. View at: Google Scholar - D. Cadamuro and J. Redondo, “Cosmological bounds on pseudo Nambu-Goldstone bosons,”
*Journal of Cosmology and Astroparticle Physics*, vol. 1202, p. 32, 2012. View at: Google Scholar - G. G. Raffelt,
*Stars as laboratories for fundamental physics: The astrophysics of neutrinos, axions, and other weakly interacting particles*, The University of Chicago Press, 1996. - J. A. Grifols and E. Masso, “Constraints on finite-range baryonic and leptonic forces from stellar evolution,”
*Physics Letters B*, vol. 173, no. 3, pp. 237–240, 1986. View at: Publisher Site | Google Scholar - Belle-II Collaboration, “Belle II technical design report,” http://arxiv.org/abs/1011.0352. View at: Google Scholar
- T. Behnke, J. E. Brau, B. Foster et al., “The international linear collider technical design report - volume 1: executive summary,” Tech. Rep., http://arxiv.org/abs/1306.6327. View at: Google Scholar
- L. Randall and R. Sundrum, “Large Mass Hierarchy from a Small Extra Dimension,”
*Physical Review Letters*, vol. 83, no. 17, pp. 3370–3373, 1999. View at: Publisher Site | Google Scholar - C. Csaki, M. L. Graesser, and G. D. Kribs, “Radion dynamics and electroweak physics,”
*Physical Review D*, vol. 63, no. 6, article 065002, 2001. View at: Publisher Site | Google Scholar - R. Contino, A. Pomarol, and R. Rattazzi,
*talk by R. Rattazzi at Planck 2010, CERN*, Madrid Christmas Workshop, 2010, talk by A. Pomarol, http://indico.cern.ch/getFile.py/access?contribId=163&resId=0&materialId=slides&confId=75810. - F. Abu-Ajamieh, J. S. Lee, and J. Terning, “The light radion window,”
*Journal of High Energy Physics*, vol. 2018, no. 10, p. 50, 2018. View at: Publisher Site | Google Scholar - F. Abu-Ajamieh, “International Journal of Modern Physics A,” vol. 33, no. 24, Article ID 1850144, 2018. View at: Publisher Site | Google Scholar
- T. Aoyama, T. Kinoshita, and M. Nio, “Theory of the anomalous magnetic moment of the electron,”
*Atoms*, vol. 7, no. 1, p. 28, 2019. View at: Publisher Site | Google Scholar

#### Copyright

Copyright © 2020 Fayez Abu-Ajamieh. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The publication of this article was funded by SCOAP^{3}.