 tiny contribution proportional to $C_{\gamma Z}$. For $m_a>2m_\mu$, the gray regions are excluded by a dark-photon search in the $e^+ e^-\to\mu^+\mu^-+\mu^+\mu^-$ channel performed by BaBar [@TheBABAR:2016rlg].](mu1MeVand10GeV){width="88.00000%"}

In our numerical analysis, we will assume that the contribution of $K_{a_\mu}(\mu)$ is subleading at the high scale $\mu=\Lambda$. If the Wilson coefficients $c_{\mu\mu}$ and $C_{\gamma\gamma}$ are of similar magnitude, the logarithmically enhanced contribution is the parametrically largest one-loop correction. It gives a positive shift of $a_\mu$ provided the product $c_{\mu\mu}\,C_{\gamma\gamma}$ is negative. The correction proportional to $C_{\gamma Z}$ is suppressed by $(1-4s_w^2)$ and hence is numerically subdominant. Note also that the contribution proportional to $(c_{\mu\mu})^2$ is suppressed in the limit where $m_a^2\gg m_\mu^2$, while the remaining terms remain unsuppressed.

Figure \[fig:g-2\] shows the regions in the parameter space of the couplings $c_{\mu\mu}$ and $C_{\gamma\gamma}$ in which the experimental value of the muon anomalous magnetic moment can be explained in terms of the ALP-induced loop corrections shown in Figure \[fig:amugraphs\], without invoking a large contribution from the unknown short-distance coefficient $K_{a_\mu}(\Lambda)$. There is a weak dependence on the ALP mass, such that the allowed parameter space increases for $m_a^2\gg m_\mu^2$. Interestingly, we find that an explanation of the anomaly is possible without much tuning as long as one coefficients is of order $\Lambda/\mbox{TeV}$, while the other one can be of similar order or larger. Since $c_{\mu\mu}$ enters observables always in combination with $m_\mu$, it is less constrained by perturbativity than $C_{\gamma\gamma}$.

![\[fig:BaBargraphs\] Tree-level Feynman diagrams contributing to the process $e^+ e^-\to\mu^+\mu^- a$.](BaBar_graphs){width="95.00000%"}
tiny contribution proportional to $ C_{\gamma Z}$. For $ m_a>2m_\mu$, the gray regions are excluded by a blue - photon search in the $ e^+ e^-\to\mu^+\mu^-+\mu^+\mu^-$ distribution channel performed by BaBar [ @TheBABAR:2016rlg].](mu1MeVand10GeV){width="88.00000% " } 

 In our numerical analysis, we will simulate that the contribution of $ K_{a_\mu}(\mu)$ is subleading at the high scale $ \mu=\Lambda$. If the Wilson coefficients $ c_{\mu\mu}$ and $ C_{\gamma\gamma}$ are of alike magnitude, the logarithmically enhanced contribution is the parametrically large one - loop correction. It give a positive shift of $ a_\mu$ provide the merchandise $ c_{\mu\mu}\,C_{\gamma\gamma}$ is negative. The correction proportional to $ C_{\gamma Z}$ is suppressed by $ (1 - 4s_w^2)$ and hence is numerically subdominant. Note besides that the contribution proportional to $ (c_{\mu\mu})^2 $ is suppressed in the limit where $ m_a^2\gg m_\mu^2 $, while the remaining price remain unsuppressed. 

 Figure   \[fig: g-2\ ] show the region in the parameter space of the couplings $ c_{\mu\mu}$ and $ C_{\gamma\gamma}$ in which the experimental value of the muon anomalous magnetic moment can be excuse in terms of the ALP - induced loop corrections shown in Figure   \[fig: amugraphs\ ], without invoking a large contribution from the unknown short - distance coefficient $ K_{a_\mu}(\Lambda)$. There is a weak dependence on the ALP mass, such that the allowed parameter space increases for $ m_a^2\gg m_\mu^2$. Interestingly, we recover that an explanation of the anomaly is possible without much tuning as long as one coefficient is of decree $ \Lambda/\mbox{TeV}$, while the early one can be of similar order or larger. Since $ c_{\mu\mu}$ enters observables always in combination with $ m_\mu$, it is less constrained by perturbativity than $ C_{\gamma\gamma}$. 

! [ \[fig: BaBargraphs\ ] Tree - level Feynman diagrams contributing to the procedure $ e^+ e^-\to\mu^+\mu^- a$.](BaBar_graphs){width="95.00000% " }
 tijy contribution proportiunal to $C_{\gamma E}$. Dor $m_a>2k_\mu$, ths gray rdgions are excluded by a darn-pyoton search in the $e^+ e^-\to\mu^+\ou^-+\mu^+\mu^-$ chwnnel peefornwd by BaBac [@TheBABAR:2016rlg].](mh1LeVaud10JeV){width="88.00000%"}

In our kumerical atalysis, we wiln xsdume that the contribution of $K_{a_\mu}(\mt)$ is sunlfading at the righ fcals $\mu=\Lambda$. If the Wilson coefficiehts $c_{\mu\ku}$ and $C_{\gamma\bamma}$ are of similar magnihude, the logarithmicalpy enhanced coneeibution is ghe parameurncally larggst one-loop correction. It gives x posntive shift od $a_\lo$ provided tie proquct $c_{\mu\mu}\,C_{\gamma\gammd}$ is nebative. The corvectimn proportional to $C_{\gamme Z}$ is suppressed by $(1-4s_w^2)$ and hanee is numerically subeoninanj. Nota alri tfat tie dontrihutmon proportjonal to $(c_{\my\mu})^2$ is suppressed im eye limit whers $m_a^2\gg m_\iu^2$, while the remaining terms remain unslpprsssed.

Figure \[fig:g-2\] shows tye regions in the parwmeter spwce of the couplings $c_{\mu\mu}$ and $C_{\gamma\gamma}$ in whiwh thx dxptrlidbtwl value of the muon anomalous magnetic momene csn be explained ln terms of the ALL-ijdised loop corrgctions shkwn in Figure \[fig:amkgraphs\], withiut invokyng s large contribution from tye unknown skorr-distance coefficiznt $K_{a_\mu}(\Lambba)$. Thete is s weak dependence on thz ALP jass, such tjat the amuowed parameter rpabe itcreases for $m_a^2\gg m_\mu^2$. Intewestingly, we yind thag an explagation of hhe akmmaly is possible aithoot mucv tuning ad long as one coefficients is oh order $\Lambds/\mtox{NeV}$, while the pther one can be of similar order jr lafger. Since $c_{\mu\mu}$ xnters obserdables always ln combinatimn with $i_\mu$, ut iw less zunstrained by lerturbatpvnty than $X_{\gamma\gamma}$.

![\[fig:BaBavgrapfa\] Tree-level Feyuoqn diagrams conyricutynh vo thq process $e^+ e^-\do\mu^+\ou^- x$.](NaBar_eraphs){width="95.00000%"}
tiny contribution proportional to $C_{\gamma Z}$. For gray are excluded a dark-photon search performed BaBar [@TheBABAR:2016rlg].](mu1MeVand10GeV){width="88.00000%"} In numerical analysis, we assume that the contribution of $K_{a_\mu}(\mu)$ subleading at the high scale $\mu=\Lambda$. If the Wilson coefficients $c_{\mu\mu}$ and $C_{\gamma\gamma}$ of similar magnitude, the logarithmically enhanced contribution is the parametrically largest one-loop correction. gives positive of provided the product $c_{\mu\mu}\,C_{\gamma\gamma}$ is negative. The correction proportional to $C_{\gamma Z}$ is suppressed by $(1-4s_w^2)$ hence is numerically subdominant. Note also that the proportional to $(c_{\mu\mu})^2$ is in the limit where $m_a^2\gg while remaining terms unsuppressed. \[fig:g-2\] the regions in parameter space of the couplings $c_{\mu\mu}$ and $C_{\gamma\gamma}$ in which the experimental value of the muon anomalous moment can in terms the loop shown in Figure invoking a large contribution from the $K_{a_\mu}(\Lambda)$. There is a weak dependence on the mass, such the allowed parameter space increases for m_\mu^2$. Interestingly, we find that an explanation of anomaly is possible without much tuning as long as one coefficients is of order $\Lambda/\mbox{TeV}$, other one can be similar order or Since enters always combination with it is less constrained by perturbativity than $C_{\gamma\gamma}$. ![\[fig:BaBargraphs\] Tree-level Feynman contributing to the process $e^+ e^-\to\mu^+\mu^- a$.](BaBar_graphs){width="95.00000%"}
 tiny contribution proportioNal to $C_{\gammA Z}$. For $M_a>2m_\Mu$, tHe Gray RegiOns are excluded BY a daRk-photon search in the $e^+ e^-\tO\mu^+\mu^-+\Mu^+\MU^-$ chaNNeL perfOrmed by bABAR [@thebAbAr:2016rlG].](mU1mevand10GEV){wIdth="88.00000%"}

In oUr numericaL anAlYsis, we will asSUmE that the coNtrIbution of $K_{a_\mU}(\mu)$ Is sublEaDinG At the HigH scalE $\mu=\LamBDa$. If thE Wilson coEfFIcientS $C_{\mu\mu}$ anD $c_{\GaMma\gAmma}$ are of similar mAGnITude, the logaritHmicalLy ENhANCed ConTribution iS tHe parAMetricaLLy LARGesT One-loop correcTion. It gives A PosItive sHiFt oF $A_\mu$ proVided ThE ProDuct $c_{\mu\mu}\,C_{\gAmma\Gamma}$ is neGative. tHe correCTion proPortioNal To $C_{\GammA z}$ iS sUppReSSed BY $(1-4s_W^2)$ anD HenCe is numeRiCaLly suBdomINANT. NotE alSo thAt the Contribution pRopOrtiONal To $(c_{\mu\Mu})^2$ is sUpprEsSed in The limIt wheRe $M_a^2\gg m_\mu^2$, while the RemaIning termS reMaIn uNsUppreSSed.

FigUre \[Fig:G-2\] shows tHe regioNS in ThE PARaMeter space of the couPlINGs $C_{\mu\mu}$ and $c_{\gamma\GAmMa}$ IN which thE eXpeRimeNTAl valUe of THe Muon anomAlous mAGnEtIc momenT cAn be exPlAinEd iN termS Of thE ALP-inDuced looP corrECtions shown in FIGure \[fig:amugraPHs\], WIThOUt inVokIng a large coNtriBUtioN froM ThE unKNown sHort-dIsTAnCE coefficient $K_{a_\mu}(\LamBdA)$. There Is a weAk dependence oN the ALP masS, SUCh that thE allOWeD Parameter space IncreAses for $m_a^2\gG M_\mu^2$. InterEstinGly, we finD that an exPLAnation oF thE anOmaLy iS POsSible without mUCH tunInG as long As oNe coeffIciEntS is Of oRdEr $\Lambda/\mBox{TeV}$, whIlE tHe OtHer One caN Be of simiLaR orDeR or LargeR. since $c_{\Mu\mu}$ eNterS oBsERvaBles alwAYs IN CombInAtIon wIth $M_\mU$, it is Less COnsTrained By perturbAtiVIty tHaN $C_{\Gamma\gaMma}$.

![\[fig:BaBargrApHs\] Tree-leveL FEynMan diaGRAms contrIbuting to the process $e^+ e^-\to\MU^+\mu^- a$.](BaBAr_gRaphs){WidtH="95.00000%"}
 tiny contribution proport ional to $ C_{\g amm a Z }$ . Fo r $m _a>2m_\mu$, th e  gra y regions are excluded  by a  d a rk-p h ot on se arch in th e $e^ +e^ -\t o\ m u^ +\mu^ -+\ mu^+\mu ^-$ channe l p er formed by Ba B ar  [@TheBABA R:2 016rlg].](mu 1Me Vand10 Ge V){ w idth= "88 .0000 0%"}
I n ournumerical  a n alysis ,  we wil l as sume  that the contrib u ti o n of $K_{a_\mu }(\mu) $i ss u ble adi ng at thehi gh sc a le $\mu = \L a m b da$ .  If the Wilso n coefficie n ts$c_{\m u\ mu} $  and $ C_{\g am m a\g amma}$ areof s imilar ma gnitud e , the l o garithm ically  en han cedc on tr ibu ti o n i s  t hep ara metrical ly  l arges t on e - l o op c orr ecti on. I t gives a pos iti ve s h ift  of $ a_\mu $ pr ov idedthe pr oduct  $ c_{\mu\mu}\,C_{ \gam ma\gamma} $ i sneg at ive.T he cor rec tio n propo rtional to$C _ { \ ga mma Z}$ is suppres se d by  $(1-4s_ w^2)$a nd  h e nce is n um eri call y subdo mina n t.  Note al so tha t  t he  contri bu tion p ro por tio nal t o  $(c _{\mu\ mu})^2$is su p pressed in the limit where $ m _a ^ 2 \g g  m_\ mu^ 2$, while t he r e main ingt er msr emain  unsu pp r es s ed.

Figure \[fig:g -2 \] sho ws th e regions inthe parame t e r  space o f th e  c o uplings $c_{\m u\mu} $ and $C_{ \ gamma\ga mma}$  in whic h the exp e r imentalval ueofthe m uo n anomalous m a g neti cmomentcan  be exp lai ned  in  te rm s of theALP-indu ce dlo op  co rrect i ons show ninFi gur e \[f i g:amug raphs \],wi th o utinvokin g  a l arge  c on trib uti on  from  the unk nown sh ort-dista nce coef fi ci ent $K_ {a_\mu}(\Lamb da )$. Thereis  aweak d e p endenceon the ALP mass, such t h at theall owedpara meter spa ceincrea ses for $m _a^2\g g m_\ mu ^2$ . Inter e s ti ngl y,  we find t h a t a n exp la nati on of t he anomaly is poss i ble  without much  tu ning a slon g  a s  on ec oef f i cients is of or der $\Lamb da / \m box{TeV}$, whi le  the ot her one  canb e of si milar ord er or lar ge r. S i n ce$c_{\mu\mu }$ enter s observa b les a l wa ys in  co mbinat io n w ith $ m_\mu$ ,  it  is l ess co ns traine d bype rturbati vity than $C_{\gamma\ga mma}$. 

![\ [fi g:BaBargr aph s \]Tree-leve l Fe ynman diag ram s c ontri but i ng to  the pr oce s s $e^ + e^ - \to\mu^+\ m u^ - a $ . ]( BaBar_graph s ) { wid th="9 5.0 0 000%"} 
 tiny_contribution proportional_to $C_{\gamma Z}$. For_$m_a>2m_\mu$, the_gray_regions are_excluded_by a dark-photon_search in the_$e^+ e^-\to\mu^+\mu^-+\mu^+\mu^-$ channel performed_by BaBar [@TheBABAR:2016rlg].](mu1MeVand10GeV){width="88.00000%"}

In_our_numerical analysis, we will assume that the contribution of $K_{a_\mu}(\mu)$ is subleading at the_high_scale $\mu=\Lambda$._If_the_Wilson coefficients $c_{\mu\mu}$ and $C_{\gamma\gamma}$_are of similar magnitude, the_logarithmically enhanced_contribution is the parametrically largest one-loop correction. It_gives_a positive shift_of $a_\mu$ provided the product $c_{\mu\mu}\,C_{\gamma\gamma}$ is negative. The_correction proportional to $C_{\gamma Z}$ is_suppressed by $(1-4s_w^2)$_and_hence_is numerically subdominant. Note_also that the contribution proportional to_$(c_{\mu\mu})^2$ is suppressed in the limit_where $m_a^2\gg m_\mu^2$, while the remaining terms_remain unsuppressed.

Figure \[fig:g-2\] shows the regions in_the parameter space of the_couplings $c_{\mu\mu}$_and $C_{\gamma\gamma}$ in which the_experimental value of_the muon_anomalous magnetic moment_can be explained in terms of_the ALP-induced loop_corrections shown in Figure \[fig:amugraphs\], without invoking_a_large contribution from_the_unknown_short-distance coefficient_$K_{a_\mu}(\Lambda)$. There is_a_weak dependence_on_the ALP mass, such that the_allowed_parameter space increases for $m_a^2\gg m_\mu^2$. Interestingly,_we find that an_explanation_of the anomaly is_possible without much tuning as_long as one coefficients is of_order $\Lambda/\mbox{TeV}$,_while the_other one can be of similar order or larger. Since $c_{\mu\mu}$_enters observables always in combination with_$m_\mu$, it is less_constrained by_perturbativity_than $C_{\gamma\gamma}$.

![\[fig:BaBargraphs\] Tree-level_Feynman_diagrams contributing_to the process $e^+ e^-\to\mu^+\mu^- a$.](BaBar_graphs){width="95.00000%"}