the standard model of elementary particles ( sm ) has been amazingly successful , and currently the large hadron collider ( lhc ) experiment is searching for higgs boson , which is the only particle that has not been observed yet in the framework of the sm . so far , some interesting hints of the higgs boson with a mass of about 125 gev have been reported at the lhc  @xcite and its discovery may not be in a far future .
the higgs boson with a mass of about 100 gev is totally consistent with the electroweak ( ew ) precision measurements .
therefore , the sm well - describes the nature up to the ew scale .
we , however , think that the sm is not the final theory of elementary particles .
the hierarchy problem has been a strong driving force to think of the physics beyond the sm , and many ideas , such as supersymmetry , extra - dimension , and technicolor , have been proposed to solve the problem .
the lhc experiment is an ideal place to probe the models for the hierarchy problem , and the searches are currently going on .
the first round of the lhc searches , however , does not show any serious deviations from the sm prediction .
it seems that this negative search result starts creating a tension between new physics models and the requirements for the hierarchy problem
. therefore it will be good time to consider other approaches to think about new physics beyond the sm .
there is a different approach to consider the physics beyond the sm , based on a consideration for unexplained experimental results .
the muon anomalous magnetic moment ( muon g-2 ) is one of the most precisely measured observables  @xcite .
the theoretical prediction from the sm by several groups has suggested that there is a discrepancy between the experimental result and the sm prediction  @xcite :    @xmath5    where @xmath6 is a discrepancy between the experimental result @xmath7 and the sm prediction @xmath8 .
there are many discussions  @xcite on uncertainties of the sm prediction to understand whether this anomaly is real or not .
if the discrepancy can not be explained by the sm , this would be an evidence of the physics beyond the sm . at present
, there seems no satisfactory explanation for this anomaly within the sm .
therefore , it is worth while considering the physics beyond the sm seriously in order to explain the anomaly .
it is interesting to note that the anomaly of muon g-2 is the same size as 1-loop contributions induced by the ew gauge bosons in the sm .
this suggests that in order to explain the anomaly , new particles with masses of ew scale are required if the interactions are of order of ew gauge couplings .
although the anomaly of muon g-2 has been discussed in the context of new physics models ( for example , see ref .
@xcite in mssm , ref .
@xcite in extra dimension model , and ref .
@xcite in little higgs model ) motivated mainly by the hierarchy problem , we take a `` bottom - up '' style approach . starting from introducing new interactions with muon to generate new contribution to muon g-2 , we try to capture important features of new physics models , although we can only discuss a part of complete models . in this paper
, we consider models where muon has new yukawa type interactions with new particles to generate new contributions to muon g-2 .
we find that in a model where the chirarity flip of muon occurs only in external line of muon g-2 diagrams , it is necessary to introduce the relatively large yukawa coupling and the ew scale new particles . on the other hand , in a model where the chirality flip can happen in the internal lines of muon g-2 diagrams
, we can obtain favorable g-2 contributions without large yukawa couplings , and they are also consistent with the precision ew observables .
we also expect these particles may be observed directly and/or indirectly at the lhc experiment .
this paper is organized as follows . in the next section ,
we consider new physics models where muon has new yukawa couplings to explain the anomaly of muon g-2 .
one model is that only right - handed muon has new yukawa coupling .
another one is that both right- and left- handed muons have new yukawa interactions to enhance the contributions to the muon g-2 .
we show the consistent region of parameter space . in sec .  3 , we discuss effects of the new physics on ew observables and whether the parameter space which is consistent with the muon g-2 is also consistent with the precision measurements . in sec .  4 , we discuss phenomenology of these new physics models at the lhc .
especially we show the possible effect on the higgs boson decay of @xmath1 , and the direct production cross sections of new particles at the lhc .
is devoted to summary .
in order to explain the deviation from the sm prediction of muon g-2 by new physics , muon has to have new interactions with some charged particles . since the effective operator for muon g-2 couples to photon , the loop effects via the new interactions of muon with charged particles could induce the extra contribution to the muon g-2 .    in this paper , as such new interactions
, we consider new yukawa - type couplings with right - handed or / and left - handed muon .
we consider two cases : ( 1 ) model where only right - handed muon has new yukawa interaction and ( 2 ) model where both right - handed and left - handed muons have new yukawa interactions .
we consider the following yukawa interaction where the right - handed muon couples to new su@xmath11 singlet fermion @xmath10 and singlet scalar @xmath9 : @xmath12 here @xmath13 is the right - handed muon .
qed charges of new fermion @xmath10 and scalar @xmath9 are @xmath14 and @xmath15 , respectively , in order to have the gauge invariant yukawa interaction .
the masses of @xmath10 and @xmath9 are denoted by @xmath16 and @xmath17 , respectively . in order to simplify the model ,
we have assumed a @xmath18 parity , under which the sm particles are even and new particles @xmath9 and @xmath10 are odd .
this type of the parity may be interesting for the dark matter .
we also note that even if we impose the @xmath18 parity , the right - handed electron , for example , can have similar yukawa interaction with @xmath9 and @xmath10 .
this will cause severe flavor mixing problem .
we will not discuss the complete model here , however , we implicitly assume that @xmath10 or @xmath9 has approximate muon flavor number so that the flavor mixing is strongly suppressed .    ,
width=566 ]    since the right - handed muon couples to @xmath10 and @xmath9 , their radiative corrections induce the muon g-2 , @xmath19 as shown in fig .
[ diagramg2 ] .
the contribution @xmath19 is given by @xmath20 , \label{muong2}\end{aligned}\ ] ] where @xmath21 and @xmath22 are momenta of external muons and @xmath23 is a photon momentum and a limit @xmath24 is taken . here
@xmath25 are so called passarino - veltman functions @xcite , and their definitions in this paper are shown in appendix .
the explicit forms are @xmath26 where @xmath27 , @xmath28 and the higher order terms of @xmath29 are neglected . the first term in eq .
( [ muong2 ] ) comes from fig .
[ diagramg2](a ) and the second is from fig .  [
diagramg2](b ) .    ) as a function of @xmath14 and @xmath17 . here
we take @xmath30 and @xmath31 gev .
we show contours for @xmath32 and @xmath33 , from right to left , corresponding to @xmath34 and @xmath35 deviation from the measured value , respectively.,width=453 ]    in fig .
[ fig_amu ] , we show the contribution to @xmath19 as a function of @xmath14 and @xmath17 . here
we take @xmath30 and @xmath31 gev .
we show contours for @xmath32 and @xmath33 , from right to left , corresponding to @xmath34 and @xmath35 deviation from the measured value , respectively .
we note that @xmath36 dependence of @xmath19 is trivial , that is , @xmath19 is proportional to @xmath37 , as shown in eq .
( [ muong2 ] ) .
for example , if one takes @xmath38 , the values of @xmath19 in fig .  [ fig_amu ]
reduce by a factor @xmath39 .
one can see that the region with @xmath40 ( which corresponds to @xmath41 ) is disfavored by the data of muon g-2 .
it is interesting to notice that the neutral fermion ( @xmath42 , corresponding to @xmath43 ) is difficult to explain the anomaly in any values of @xmath17 , and on the other hand , the neutral scalar ( @xmath44 , corresponding to @xmath45 ) can potentially accommodate the anomaly if the scalar @xmath9 is not so heavy .
it is also interesting to note that the multi - charged fermion and scalar ( such as @xmath46 , corresponding to @xmath47 ) are also favored by the muon g-2 anomaly if @xmath17 is a right value for the anomaly .
therefore , the anomaly of muon g-2 constrains the qed charges as well as the mass scale of new particles .    2 ( @xmath19 ) as a function of @xmath16 and @xmath17 for ( a ) @xmath45 and and @xmath30 ( b ) @xmath48 and @xmath49 , respectively .
contours for @xmath32 and @xmath33 ( corresponding to @xmath34 and @xmath35 deviation from the measured value ) are shown , similar to fig .
[ fig_amu ] .
, title="fig:",width=321 ] 2 ( @xmath19 ) as a function of @xmath16 and @xmath17 for ( a ) @xmath45 and and @xmath30 ( b ) @xmath48 and @xmath49 , respectively .
contours for @xmath32 and @xmath33 ( corresponding to @xmath34 and @xmath35 deviation from the measured value ) are shown , similar to fig .
[ fig_amu ] .
, title="fig:",width=321 ]    in fig .
[ fig_amu_2 ] , we show @xmath19 as a function of @xmath16 and @xmath17 for ( a ) @xmath45 @xmath50 and @xmath30 and ( b ) @xmath48 @xmath51 and @xmath49 .
contours for @xmath32 and @xmath33 ( corresponding to @xmath34 and @xmath35 deviation from the measured value ) are shown , similar to fig .
[ fig_amu ] .
, one should simply multiply the value of @xmath19 in fig .
[ fig_amu_2 ] ( a ) and ( b ) by a factor @xmath52 and @xmath53 , respectively . ]
as one can see , new particles with their masses of about @xmath54 gev and the coupling @xmath55 are strongly favored by the anomaly of muon g-2 when @xmath45 . when the qed charge of @xmath10 decreases further negatively ( from @xmath56 to @xmath57 )
, the muon g-2 gets larger even if the yukawa coupling @xmath36 gets smaller , as shown in fig .
[ fig_amu_2 ] ( b ) .
since the anomaly of muon g-2 requires relatively light new particles ( @xmath58 ) and/or relatively large yukawa coupling ( @xmath59 ) in this scenario , the ew precision observables can be affected by these new particles and new interactions .
therefore , we will check the effects of these particles on the ew observables in a next section .
we would also like to point out that the anomaly of muon g-2 suggests the existence of relatively light new particles and/or multi - charged particles .
therefore , it is very interesting to know whether these particles can be found directly or indirectly at the lhc .
we will study the direct productions of these particles and the effects on higgs decay to @xmath60 at the lhc in a later section .      unlike the previous case , here we consider that both right- and left - handed muons have new yukawa interactions : @xmath61 where @xmath62 , @xmath63 and @xmath9 are the second generation su(2 ) doublet lepton , su(2 ) doublet and singlet scalars , respectively , and @xmath10 is su(2 ) singlet fermion , whose mass is @xmath16 .
here we have assumed the @xmath18 parity again in order to simplify the model . under the @xmath18 parity ,
the sm particles are even , and new particles @xmath9 , @xmath64 and @xmath10 are odd . we also implicitly assume the approximate muon flavor symmetry , under which @xmath10 or ( @xmath65 ) has muon flavor number , so that the flavor mixing of this type of new yukawa couplings are strongly suppressed .
the model we discuss in this paper may be a part of the complete model .
however , we think this part of the complete model is crucial for the muon g-2 .
the qed charges of new particles are represented in term of qed charge of @xmath10 field ( @xmath14 ) as follows : @xmath66 since @xmath67 , @xmath9 and @xmath68 can mix each other .
for example , the following gauge invariant term induces the @xmath69 mixing mass terms after the electroweak symmetry is spontaneously broken : @xmath70 here we parameterize the mass terms for @xmath9 and @xmath68 as follows : , mass term such as @xmath71 is also possible .
here we neglect such a mass term for simplicity .
even if we include such a mass term , our result does not change qualitatively . ]
@xmath72 thus diagonalizing this mass matrix , we define the mass eigenstates @xmath73 as @xmath74 here @xmath75 is a mixing unitary matrix , which diagonalizes the mass matrix shown in eq .
( [ mass_matrix ] ) .
the mass eigen values are taken to be @xmath76 .
this type of mixing terms are important to induce large contribution to muon g-2 .
the contributions to muon g-2 are summarized by @xmath77 where @xmath21 and @xmath22 are momenta of external muons , and @xmath23 is a photon momentum and a limit @xmath24 is taken .
here we show the explicit expressions of above passarino - veltman functions : @xmath78 where @xmath79 and we neglected the higher order of @xmath80 .    the effective operator which expresses the muon g-2 is written by @xmath81 where @xmath82 is a vacuum expectation value of higgs boson , and @xmath83 is a field strength of photon field , and @xmath84 is a typical scale related to new physics .
as one see , the chirality of muon has to flip in this interaction . in the case where only right - handed muon has the new yukawa interaction , the chirality flipping of muon happens in the external muon line in the loop diagram . on the other hands , in the present case where both right- and left - handed muons have the new yukawa interactions ,
the chirality flipping can occur in the internal fermion line , which is proportional to the mass of the fermion @xmath10 .
that is why there are terms which are proportional to @xmath85 in eq .
( [ muong2_model2 ] ) .
this contribution enhances the effects of muon g-2 .
therefore , it is very important to explain the anomaly of muon g-2 .    ) as a function of @xmath16 and @xmath86 in case of ( a ) @xmath45 and ( b ) @xmath42 .
contours for @xmath32 and @xmath33 ( corresponding to @xmath34 and @xmath35 deviation from the measured value ) are shown . here
we assume that @xmath87 , @xmath88 .
we also take ( a ) @xmath89 and ( b ) @xmath90 . , title="fig:",width=321 ] ) as a function of @xmath16 and @xmath86 in case of ( a ) @xmath45 and ( b ) @xmath42 .
contours for @xmath32 and @xmath33 ( corresponding to @xmath34 and @xmath35 deviation from the measured value ) are shown .
here we assume that @xmath87 , @xmath88 .
we also take ( a ) @xmath89 and ( b ) @xmath90 . ,
title="fig:",width=321 ]    in fig .
[ muong2_model2 - 1 ] , we show @xmath19 as a function of @xmath16 and @xmath86 ( which is a mass of lighter scalar state ) in case of ( a ) @xmath91 and ( b ) @xmath42 .
contours for @xmath92 , 10.1 , 18.1 , 26.1 , 34.1 , 42.1 and 50.1 ( corresponding to @xmath34 and @xmath35 deviation from the measured value ) are shown .
here we took @xmath93 , @xmath94 , and ( a ) @xmath89 ( b ) @xmath90 .
note that in order to generate the positive contribution to muon g-2 in fig .
[ muong2_model2 - 1 ] ( a ) , the sign of new yukawa couplings had to be taken as @xmath95 .
on the other hand , in fig .
[ muong2_model2 - 1 ] ( b ) , the sign of @xmath96 should be positive . as shown in fig .
[ muong2_model2 - 1 ] , compared to fig .
[ fig_amu_2 ] , even smaller yukawa couplings and heavier new particles can accommodate the anomaly of muon g-2 because of the enhancement mentioned above . as can be seen from fig .
[ muong2_model2 - 1 ] , when @xmath97 , the ew scale new particles are expected in order to explain the anomaly of muon g-2 . therefore , we will analyze the effects on the ew precision observables and show the consistent region of parameters in a next section .
in the previous section , we showed that the relatively light new particles are required in order to explain the anomaly of muon g-2 .
in addition , in the case where the only right - handed muon has new yukawa interaction , the new yukawa coupling should be relatively large for the anomaly of muon g-2 .
therefore , we should check if the relatively light particles and the relatively large yukawa coupling are consistent with the ew precision measurements . although models which we consider in this paper may be only a part of complete models
, we should know how these new particles affect ew observables .    in this section
, we adopt the formalism in refs.@xcite in order to include the oblique corrections as well as the vertex corrections in the ew observables .
first we briefly summarize the formalism  @xcite .    in the presence of ew to tev scale physics
, it is well known that the oblique corrections in gauge boson self - energy are important .
they are parameterized by peskin - takeuchi s s , t , u parameters  @xcite : @xmath98 here we use the notation @xmath99 and @xmath100 to refer to the cosine and sine of the weak mixing angle and @xmath101 .
in addition , in refs.@xcite , @xmath102 and @xmath103 parameters are introduced in order to account for the smaller corrections : @xmath104 furthermore , new particles have effects on the running qed coupling constant @xmath105  @xcite : @xmath106 here @xmath107 , @xmath108 , @xmath109 , and @xmath110 are the leptonic , the five - flavor hadronic , the top quark and the new physics contributions to the running of the qed coupling constant respectively .
the new physics contributions to the running of the qed coupling constant @xmath110 are defined as @xmath111 these oblique parameters contribute to ew observables , as shown in refs .
@xcite .
if muon has @xmath112 new yukawa couplings , extra yukawa couplings , for simplicity . ]
it induces non - universal vertex corrections to the @xmath113 coupling .
the standard model coupling of @xmath113 is given by @xmath114,\end{aligned}\ ] ] where @xmath115 are the standard model couplings , and their tree level contributions are @xmath116 and @xmath117 . when we take into account of corrections to @xmath113 vertex via the new yukawa coupling ( including the corresponding wave function renormalization ) ,
the @xmath113 coupling is modified by @xmath118.\end{aligned}\ ] ] here the vertex corrections generated by new particles are parameterized by @xmath119 .
the corrections to @xmath120 vertex generates the corrections to the fermi constant @xmath121 and therefore , we parameterize it as @xmath122 .
the @xmath122 for the new physics contributions from the vertex and box diagrams to the @xmath123-decay process is defined as @xmath124 here @xmath125 is the muon decay constant including effects of sm radiative corrections and new physics oblique corrections except the vertex and box corrections from new physics .
.experimental data and theoretical predictions of electroweak observables .
the predicted values of the sm and sample model discussed in sec .
3.2 are shown . [ cols="^,^,^,^,^,^",options="header " , ]      of the ew observables listed in table .
[ ewobs ] as a function of the higgs boson mass @xmath126 within the sm .
( b ) @xmath127 contours shown in dashed lines in s - t plane , assuming that other oblique corrections ( @xmath128 , @xmath103 , @xmath102 ) and vertex corrections are zero .
reference higgs boson mass and top quark mass are taken to be @xmath129 gev and 173.2 gev , respectively .
we also show the predicted s - t values ( shown in solid lines ) in the sm , varying the higgs boson mass and top quark mass from the reference values .
, title="fig:",width=302 ]   of the ew observables listed in table .
[ ewobs ] as a function of the higgs boson mass @xmath126 within the sm .
( b ) @xmath127 contours shown in dashed lines in s - t plane , assuming that other oblique corrections ( @xmath128 , @xmath103 , @xmath102 ) and vertex corrections are zero .
reference higgs boson mass and top quark mass are taken to be @xmath129 gev and 173.2 gev , respectively .
we also show the predicted s - t values ( shown in solid lines ) in the sm , varying the higgs boson mass and top quark mass from the reference values .
, title="fig:",width=302 ]    using the formalism in ref.@xcite , one can calculate ew observables listed in table  [ ewobs ] without assuming lepton universality .
first , we show a result of the sm fit of the ew observables in fig .
[ smfit ] ( a ) as a function of higgs boson mass @xmath126 .
the best fit point is at @xmath130 gev .
this is consistent with the result shown by the lep electroweak working group  @xcite . in table
[ ewobs ] , we show the sm fit when higgs boson mass is assumed to be 125 gev , which has been suggested by the latest lhc data . in fig .
[ smfit](b ) , @xmath127 contours are shown in dashed lines in s - t plane . here
we assume that other oblique corrections ( @xmath128 , @xmath103 , @xmath102 ) and vertex corrections are zero .
reference higgs boson mass and top quark mass are taken to be 125 gev , 173.2 gev , respectively .
the predicted values of s - t parameters are also shown by varying the higgs boson mass and top quark mass from the reference values .
we note that the sm fit with lighter higgs boson is good except for the muon g-2 , and small positive s and t @xmath131 ) can decrease the @xmath127 further .
therefore , we need new physics which largely contributes to muon g-2 , but whose effects on other ew observables are small . in this section ,
we analyze the ew observables in models discussed in the previous section . in this analysis , we assume @xmath132 gev .
since the new scalar @xmath9 and new fermion @xmath10 are su(2 ) singlet , they do not couple to w boson .
but if it has a qed charge , it can couple to photon and z boson . here
we list the corrections to self - energy functions of gauge bosons in this model .
the singlet fermion ( @xmath10 ) contributions are given by @xmath133 , \nonumber\\ \pi_{\gamma \gamma}^{(\chi)}(p^2)&= & -\frac{e^2}{4\pi^2}q_\chi^2\left [ m^2_{\chi}b_0(\chi,\chi)-p^2\left\{b_1(\chi,\chi ) + b_{21}(\chi,\chi)\right\ } -2(1-\epsilon)b_{22}(\chi,\chi ) \right ] , \nonumber \\
\pi_{\gamma z}^{(\chi)}(p^2 ) & = & \frac{geq_\chi^2s_w^2}{4\pi^2c_w } \left [ m_\chi^2b_0(\chi,\chi)-p^2\left\ { b_1(\chi,\chi)+b_{21}(\chi,\chi ) \right\ } -2(1-\epsilon)b_{22}(\chi,\chi ) \right ] .
\nonumber\\ \label{singlet_chi}\end{aligned}\ ] ] here @xmath134 , which are passarino - veltman functions and are shown explicitly in appendix .
we use the dimensional regularization in space - time dimension @xmath135 .
the singlet scalar ( @xmath9 ) contributions are given by @xmath136 here passarino - veltman functions @xmath137 and @xmath138 are shown in appendix .
we can easily show that peskin - takeuchi stu parameters in this case vanish ( @xmath139 ) , since the new particles do not have su(2)@xmath140 interactions .
thus , the leading contributions to the oblique corrections are @xmath102 parameter and @xmath110 : @xmath141\right\}\right ] \nonumber \\ & & + \frac{s_w^4 q_\phi^2}{3\pi}\left[1-\frac{12m_\phi^2}{m_z^2 } \left\{1-\sqrt{\frac{4m_\phi^2}{m_z^2}-1 } \tan^{-1}\left[\frac{1}{\sqrt{\frac{4m_\phi^2}{m_z^2}-1}}\right]\right\}\right],\end{aligned}\ ] ]    @xmath142\right ] \nonumber \\ & & -\frac{2\alpha q_\phi^2}{9\pi}\left [ 1-\frac{3m_\phi^2}{m_z^2 } + \frac{3}{4}\sqrt{\frac{4m_\phi^2}{m_z^2}-1 } \tan^{-1}\left[\frac{1}{\sqrt{\frac{4m_\phi^2}{m_z^2}-1}}\right]\right],\end{aligned}\ ] ]    where we assume @xmath143 and @xmath144 .    as discussed in the previous section , in order to explain the anomaly of muon g-2 , the new yukawa coupling in this model should be relatively large .
in this case , the potentially large vertex corrections may be expected , as shown in fig.[vertex_fig ] .
coupling.,width=566 ]    the results are expressed by @xmath145,\end{aligned}\ ] ] where @xmath21 , @xmath146 and @xmath23 are muon , anti - muon and z - boson momenta , respectively . here
passarino - veltman functions @xmath147 and @xmath148 are explicitly shown in appendix .    .
@xmath149 are shown as a function of @xmath16 and @xmath17 in case of @xmath45 and @xmath30 , where @xmath150 includes all corrections , and @xmath151 contains all corrections except for the vertex corrections .
( b ) @xmath127 as a function of @xmath16 and @xmath17 in case of @xmath45 and @xmath30.,title="fig:",width=302 ] .
@xmath149 are shown as a function of @xmath16 and @xmath17 in case of @xmath45 and @xmath30 , where @xmath150 includes all corrections , and @xmath151 contains all corrections except for the vertex corrections .
( b ) @xmath127 as a function of @xmath16 and @xmath17 in case of @xmath45 and @xmath30.,title="fig:",width=302 ]    .
@xmath149 is shown as a function of @xmath16 and @xmath17 in case of @xmath48 and @xmath49 , similar to fig .
[ vertex](a ) .
( b ) @xmath127 as a function of @xmath16 and @xmath17 in case of @xmath48 and @xmath49 .
, title="fig:",width=302 ] .
@xmath149 is shown as a function of @xmath16 and @xmath17 in case of @xmath48 and @xmath49 , similar to fig .
[ vertex](a ) .
( b ) @xmath127 as a function of @xmath16 and @xmath17 in case of @xmath48 and @xmath49 .
, title="fig:",width=302 ]    the vertex corrections affect ew observables such as @xmath150 , @xmath152 , and @xmath153 . among these observables , we find that the effect on the @xmath150 is the most important in the fit because @xmath150 is most precisely measured . in fig.[vertex ] ( a ) , we show the effect of vertex corrections in @xmath150 , that is , @xmath154 as a function of @xmath16 and @xmath17 in case of @xmath45 and @xmath30 . here
@xmath150 is a theoretical prediction in this model containing all corrections , on the other hand , @xmath151 contains all corrections except for the vertex corrections .
the difference @xmath155 shows the effect of the vertex corrections .
note that the size of experimental error of @xmath150 is @xmath156 for @xmath157 , as shown in table .
[ ewobs ] . as one can see in fig .
[ vertex](a ) , the vertex correction can change the prediction of @xmath150 about @xmath157 in the region where the muon g-2 can be explained by the new physics contributions . since the sm prediction of @xmath150 is smaller than that of the experimental result ( shown in table .
[ ewobs ] ) , the vertex corrections do not help to improve the prediction of @xmath150 . therefore , the region of small @xmath10 mass is very constrained .    in fig .
[ vertex ] ( b ) , @xmath127 is shown as a function of @xmath16 and @xmath17 in case of @xmath45 and @xmath30 .
as one can see from the figure , @xmath127 is large in the region of small @xmath16 because of the large vertex corrections , discussed above . in the region of larger @xmath16 and @xmath17 , it is difficult to explain the anomaly of muon g-2 because new particles are too heavy , and hence the @xmath127 gets larger .
the minimum of @xmath127 is around @xmath158 gev and @xmath159 gev , and we see that the light new particles are favored by the ew observables including muon g-2 .    in fig .
[ vertex2 ] ( a ) , the effect of vertex corrections on @xmath150 , @xmath155 is shown in case of @xmath48 and @xmath49 , similar to fig .
[ vertex ] ( a ) .
the vertex corrections can be as large as @xmath157 of the experimental error of @xmath150 , and they increase the @xmath127 . in fig .
[ vertex2 ] ( b ) , the @xmath127 is shown as a function of @xmath16 and @xmath17 .
as one can see , the region of small @xmath16 is disfavored by the vertex corrections of @xmath150 , as shown in fig .
[ vertex2 ] ( a ) .
the minimum of the @xmath127 is around @xmath160 gev and @xmath159 gev .
therefore , the ew precision measurements ( including muon g-2 ) prefer the relatively light new particles in this case .         contours ( without including muon g-2 result ) in s - t plane ( @xmath175 to 34 in dashed lines ) and s - t values ( solid lines with points ) in this model . here
we take @xmath89 , @xmath45 ( corresponding to @xmath176 and @xmath177 ) , and @xmath178 .
two solid lines correspond to two different values of @xmath179 , @xmath88 and @xmath180 , varying @xmath181 from 200 gev to 400 gev at 50 gev step .
, width=491 ]    in fig .
[ st_para1 ] , we show numerical result of @xmath127 contours in s - t plane ( without including muon g-2 result ) , shown in dashed lines . to draw the @xmath127 contours in s - t plane ,
we assume that other oblique corrections ( @xmath128 , @xmath102 , and @xmath103 ) and vertex corrections are zero . as one can see from the figure , slightly positive s and t ( @xmath182 and @xmath183 )
are favored by ew observables . in fig .
[ st_para1 ] , we also show the predicted values of s - t parameters in this model , shown in solid lines .
here we took @xmath89 , @xmath45 ( corresponding to @xmath176 and @xmath177 ) , and @xmath178 .
two solid lines correspond to two different values of @xmath179 , @xmath88 and @xmath180 .
the points on the solid lines represent the predicted values of s and t , varying @xmath181 from 200 gev to 400 gev at 50 gev step .
when @xmath178 and @xmath88 , @xmath184 gev and @xmath185 gev . as roughly shown in eqs .
( [ st ] ) , when @xmath186 gev and @xmath177 , @xmath187 and @xmath188 , and hence @xmath189 and @xmath190 .
as @xmath181 becomes larger and about 300 gev , both @xmath191 and @xmath192 get closer to zero because the doublet scalars are almost degenerate . then when @xmath181 is increased further , @xmath191 gets negative because @xmath193 , but @xmath190 . in fig .
[ st_para1 ] , one can see this behavior .
even when we increase @xmath179 , we can see that this behavior is almost same .
except for @xmath90 and @xmath42 ( corresponding to @xmath194 and @xmath195 ) .
, width=491 ]    in fig .
[ st_para2 ] , the same figure is shown except for @xmath90 and @xmath42 ( corresponding to @xmath194 and @xmath195 ) . when @xmath42 ( @xmath195 ) and @xmath186 gev , @xmath196 and @xmath187 , so that @xmath191 is negative .
as @xmath181 gets larger , @xmath191 gets larger . for @xmath197 gev ,
@xmath198 because the doublet scalars are almost degenerate . for @xmath199 gev
, @xmath191 becomes positive .
this behavior is different from the previous case with @xmath177 since the sign of @xmath174 is different .
note that @xmath192 is always positive unless su(2 ) scalar doublet is degenerate .
as one can see in both cases , the small non - degeneracy of su(2 ) scalar doublet can make @xmath127 better .    in our numerical analysis
, we also include the vertex corrections in @xmath113 , @xmath200 and @xmath201 vertices even though the yukawa couplings @xmath202 and @xmath203 can be small in this model .
the expression of the vertex corrections is summarized as follows : @xmath204,\\ \delta g_r^\mu & = & \frac{y_r^2}{16\pi^2}\left [ 2 \sum_{ij } \left\{\left(-\frac{1}{2}-q_\phi s_w^2\right)v_{2i}^ * v_{2j}-q_\phi     s_w^2 v_{1i}^ * v_{1j}\right\ } v_{1i}v_{1j}^ * c_{24}(s_i,\chi , s_j ; p , q - p)\right.\nonumber \\ & & + \sum_i q_\chi s_w^2    \frac{1}{2}-2c_{24}-m_z^2 ( c_{12}+c_{23})+m_\chi^2 c_0\right\}(\chi , s_i,\chi;q - p , p )   \nonumber \\ & & \left .
-\sum_i s_w^2|v_{1i}|^2(b_0+b_1)(s_i,\chi;p ) \right],\\ \delta g_l^{\nu_\mu } & = & \frac{y_n^2}{16\pi^2}\left [ 2\left(\frac{1}{2 } -s_w^2 q_{\phi_1 } \right ) c_{24}(\phi_1,\chi,\phi_1;q - p , p ) \right.\nonumber\\ & & + q_\chi s_w^2\left\ { \frac{1}{2}-2 c_{24}-m_z^2(c_{12}+c_{23})+m_\chi^2 c_0 \right\}(\chi,\phi_1,\chi;q - p , p ) \nonumber \\ & & \left .
-\frac{1}{2}(b_0+b_1)(\phi_1,\chi;p ) \right],\end{aligned}\ ] ] where @xmath21 and @xmath23 are momenta for muon in @xmath205(muon neutrino in @xmath206 ) and z - boson , respectively .
correction to @xmath207-vertex is written by @xmath208,\end{aligned}\ ] ] where @xmath21 , @xmath22 and @xmath23 are momenta for muon , muon - neutrino , and w - boson , respectively .
here we use approximation @xmath28 .
contours as a function of @xmath16 and @xmath86 in case of ( a ) @xmath45 and @xmath89 and ( b ) @xmath42 and @xmath90 . here
we take @xmath87 and @xmath88.,title="fig:",width=302 ]   contours as a function of @xmath16 and @xmath86 in case of ( a ) @xmath45 and @xmath89 and ( b ) @xmath42 and @xmath90 . here
we take @xmath87 and @xmath88.,title="fig:",width=302 ]    in fig .
[ chi2_model2 ] , we show @xmath127 contours as a function of @xmath16 and @xmath86 in case of ( a ) @xmath45 and @xmath89 and ( b ) @xmath42 and @xmath90 . here
we have taken @xmath88 , and @xmath87 .
@xmath181 is chosen so that the @xmath127 becomes the smallest , provided that @xmath209 for @xmath45 and @xmath210 for @xmath42 respectively . in both cases shown in fig .
[ chi2_model2 ] , the @xmath127 can be as small as 22 , which is smaller than those in the previous case where only right - handed muon has new yukawa coupling . as seen from fig .
[ chi2_model2 ] , the relatively light @xmath10 and @xmath211 are favored by the ew observables . in table .  1
, we show the predicted values of the ew observables in case of @xmath42 , @xmath212 gev , @xmath213 gev , @xmath214 gev@xmath215 , and @xmath216 gev@xmath215 . as shown in table .  1
, the fit is quite good .
therefore , it will be important to know what kind of effects we expect from these particles at the lhc . in a next section
, we will discuss the possible effects of these new particles at the lhc .
in the previous section , we show that in order to explain the anomaly of muon g-2 , it is strongly suggested that there should be the ew scale new particles . it is also suggested that multi - charged particles may be favored by the anomaly of muon g-2 . therefore , these particles may be reachable directly and/or indirectly at the lhc . here
we discuss the effects of these new particles in the higgs decay @xmath1 and their direct productions at the lhc .
the su(2 ) singlet and doublet scalars can couple to the higgs boson through the following scalar interactions : @xmath217 where @xmath218 is a standard model higgs boson doublet . in case where only su(2 ) singlet scalar @xmath9 exists , only coupling @xmath219 is non - zero , but in case where both su(2 ) singlet @xmath9 and doublet @xmath64 scalars are there , all couplings are allowed ( @xmath220 ) is zero , terms such as @xmath221 and @xmath222 are also possible . in our analysis , we assume that such terms are small , for simplicity .
even if we include such terms , our result does not change qualitatively . ] .    in the case where only right - handed muon has new yukawa interaction ,
there is su(2 ) singlet @xmath9 , and its interaction with higgs boson is @xmath223 where @xmath224 in this model and @xmath82 is a vacuum expectation value of higgs field , @xmath225 gev .
when @xmath220 , which is the qed charge of @xmath9 , is non - zero , @xmath9 can contribute to the higgs decay @xmath1 .
the decay width @xmath226 is given by @xmath227 where the amplitude @xmath228 at leading order is written by @xmath229 here the first , the second , and the third terms are top , w , and @xmath9 contributions , respectively , and @xmath230 .
the function @xmath231 are given by @xmath232,\\ f_w(\tau)&=&2 + 3\tau^{-1}+3\tau^{-1}(2-\tau^{-1})f(\tau),\\ f_\phi(\tau_\phi ) & = & \tau^{-1}\left[-1+\tau^{-1}f(\tau)\right],\end{aligned}\ ] ] where @xmath233 for @xmath234 .
branching ratio of @xmath1 is approximately given by @xmath235 because the total decay width does not change much , even if @xmath226 changes . here
@xmath236 and @xmath237 are the sm predictions at the leading order .
normalized by the sm prediction @xmath238 as a function of @xmath17 and @xmath239 .
here we take @xmath132 gev . ,
width=377 ]    notice that this ratio is a function of @xmath17 and @xmath240 when higgs boson mass @xmath126 is fixed . as we discussed earlier , in order to explain the anomaly of muon g-2 , larger @xmath220 may be favored . as can be seen in eq .
( [ hgg ] ) , larger @xmath220 induces larger effect on @xmath1 .    in fig .
[ higgs_lhc ] , the ratio @xmath241 is shown as a function of @xmath17 and @xmath239 .
here we took @xmath132 gev .
as one can see from the figure , the effect can be significant if @xmath239 is large and @xmath17 is @xmath242 gev .
for example , as we discussed , when @xmath243 , @xmath244 and @xmath245 gev , @xmath246 .
unfortunately , we can not predict the branching ratio @xmath247 in this model because @xmath248 is unknown .
however , it is not surprising that @xmath247 is very different from the sm prediction .    in the case where both right- and left - handed muons have new yukawa interactions , the su(2 ) singlet and doublet scalars have interactions with higgs boson , as shown in eq .
( [ higgs_int ] ) and they generates the following interactions : @xmath249 where @xmath250
they contribute to the decay width of @xmath1 ,    @xmath251    when @xmath45 , which corresponds to @xmath176 and @xmath177 ( as an example we discussed above ) , only @xmath252 can contribute to @xmath1 , because @xmath211 does not couple to photon .
the result is same as one in fig .
[ higgs_lhc ] for @xmath253 and @xmath254 . in this case
, the effect can be significant if @xmath255 is large .
for example , when @xmath256 gev and @xmath257 and @xmath56 , @xmath258 and 1.05 , respectively .
when @xmath42 , which corresponds to @xmath194 and @xmath195 , scalars @xmath73 can contribute to @xmath1 .
similarly to the previous case , if the higgs couplings @xmath259 and @xmath260 are large , the effect on @xmath1 can be significant . for example , when @xmath87 and @xmath261 , we obtain @xmath262 and @xmath263 for @xmath264 , on the other hand , @xmath265 and @xmath266 for @xmath267 . as a result , for @xmath178 , @xmath268 and @xmath269 for @xmath270 and @xmath56 ( 0.71 and 1.56 for @xmath271 and @xmath272 ) , respectively .
although it is difficult to predict the decay branching ratio of @xmath1 in this model because it depends on many unknown parameters in scalar interactions , it is possible to have large effect on @xmath1 process .      as we have discussed in the previous sections , the anomaly of muon g-2 strongly suggests the existence of ew scale new particles .
models we discussed in this paper may not be complete , namely , there may be other new particles in the complete models . therefore , in this paper , we can not discuss the well - defined signatures at the lhc . however , we think that it is important to know how many events of new particles can be generated at the lhc . here
we present the production cross sections of new particles at the lhc .     in case of @xmath45 :   @xmath273 , and ( b ) @xmath274 in case of @xmath195 :   @xmath275 at the lhc ( tevatron ) , as a function of @xmath16 and @xmath276 , respectively .
here we take @xmath277 tev and 14 tev for lhc and 1.96 tev for tevatron .
, title="fig:",width=302 ]   in case of @xmath45 :   @xmath273 , and ( b ) @xmath274 in case of @xmath195 :   @xmath275 at the lhc ( tevatron ) , as a function of @xmath16 and @xmath276 , respectively . here
we take @xmath277 tev and 14 tev for lhc and 1.96 tev for tevatron . , title="fig:",width=302 ]    in the case where the right - handed muon has new yukawa coupling , singlet scalar ( @xmath9 ) and singlet fermion ( @xmath10 ) exist in the model . when @xmath45 , only @xmath10 couples to the standard model gauge bosons ( @xmath278 and @xmath279 ) because @xmath44
. therefore , the fermion @xmath10 can be pair - produced via ew gauge interactions .
we calculate the pair production cross section of @xmath10 at the lhc ( @xmath280 ) and at tevatron ( @xmath281 tev ) as a function of @xmath16 , as shown in fig .
[ fig_lhc](a ) . to calculate the production cross sections , we have used the madgraph  @xcite .
in order to explain the anomaly of muon g-2 , the mass of the fermion @xmath10 should be of order of about 100 gev . from fig .
[ fig_lhc](a ) , the production cross section of @xmath10 is in the range of 0.5 - 0.0005 pb ( 1 - 0.001 pb ) for @xmath282 gev and @xmath277 tev ( @xmath283 tev ) at the lhc and also 0.1 - 0.001 pb for @xmath284 gev and @xmath281 tev at the tevatron . if we consider a luminosity of 10 fb@xmath285 for @xmath277 tev and 14 tev , the number of events produced at the lhc is of order of 5 - 5000 and 10 - 10000 , respectively .
therefore , number of signal events can be significant . in this case , if @xmath9 is lighter than @xmath10 , @xmath10 can decay to @xmath123 and @xmath9 .
if @xmath9 is a stable neutral particle , the final signature may be a @xmath286 + missing energy .
this will suffer from the standard model background such as @xmath287 production whose production rate is about 0.29 ( 0.51 ) pb for @xmath277 tev ( 14 tev ) at the lhc .
therefore , the detail study will be very important to discuss the discovery potential of such a signal event .
however , the detail signature will depend on the particle spectrum in the complete model .
thus we will leave the detail study .
if we consider multi - charged @xmath10 , the cross section increases by a factor @xmath288 , and the cross section further becomes larger , and hence the constraints from the lhc will be very important .
pair production is also not negligible in this case . ]    in the second model where one has a @xmath289 doublet scalar in addition to @xmath9 and @xmath10 , the scalars can be pair - produced via ew gauge interaction .
for example , fig .
[ fig_lhc ] ( b ) shows the production cross sections for the charged @xmath290 scalars as a function of @xmath276 in case of @xmath42 , which corresponds to @xmath194 and @xmath195 .
we show the predicted values for @xmath277 tev and 14 tev at the lhc and @xmath281 tev at the tevatron .
as can be seen from the figure , the cross sections are in the range of 0.1 - 0.0001 pb ( 0.2 - 0.0005 pb ) for @xmath291 gev at @xmath277 tev ( @xmath283 tev ) and also 0.02 - 0.0001 pb for for @xmath292 gev at the tevatron . for simplicity
, we have assumed that there are no @xmath9-@xmath68 scalar mixing . in this case , @xmath293 are also produced via @xmath294-boson .
if masses of @xmath289 doublet scalar are degenerated , we have checked that the cross section is similar to one in fig .
[ fig_lhc ] ( b ) .
in general , the signal events have the form @xmath295 where @xmath296 and @xmath297 .
therefore , fig .
[ fig_lhc ] shows the typical values of the cross section in this model .
although we can not discuss the complete signatures because our sample models are not complete , we stress that the searches for the production processes of these new particles via ew interactions are very important in models where the anomaly of muon g-2 can be explained by these new particles .
the current lhc experiment does not show any deviation from the sm predictions .
this result puts severe constraints on new physics models .
especially , colored particles are strongly constrained since the lhc as a hadron collider can easily produce the colored particles via qcd interactions .
for example , in the constrained minimal supersymmetric standard model ( cmssm ) , gluino as well as the first and the second generation squarks are strongly constrained and their typical lower mass limits are about 1 tev  @xcite .
therefore , models motivated by the hierarchy problem are getting more constrained , and especially cmssm as a solution to the hierarchy problem is strongly disfavored .
if we motivate new physics models from the reported anomaly of muon g-2 , our naive expectation is the existence of ew scale new particles because the size of anomaly of muon g-2 is as large as one induced by ew gauge bosons in the sm . in this paper
, we analyzed two examples , where muon has new yukawa type interactions with new particles .
one example was a model where right - handed muon has new yukawa interaction with su(2 ) singlet scalar and fermion . in this case , in order to explain the anomaly of muon g-2 as well as to satisfy the precision ew measurements , we found that the relatively large new yukawa coupling and ew scale new particles are strongly favored .
we also noted that the multi - charged particles ( whose qed charges are large ) are favored .
another example we analyzed was a model where both right- and left - handed muons have new yukawa interactions with su(2 ) doublet scalar and singlet scalar and fermion , so that the chirality flip of muon can occur in the internal fermion line of 1-loop muon g-2 diagrams . in this case ,
new contributions to muon g-2 are enhanced , and hence even smaller new yukawa couplings can explain the anomaly of muon g-2 .
we also showed that the ew observables are also well - fitted by the new physics contributions . in this case , the ew scale new particles are favored by the data .
although these two examples may be only a part of complete models , the muon g-2 as well as the ew observables can constrain the qed charge of new particles as well as the scale of new particles and the size of new yukawa couplings .
since the ew scale new particles can explain the anomaly of muon g-2 and they can be consistent with the precision ew measurements , it is important to study the possible phenomenology at the lhc .    since some of new particles have qed charges as well as couplings with higgs boson , the branching ratio of @xmath298 can be effected . especially , as we showed , multi - charged particles favored by the muon g-2 can significantly contribute to @xmath298 , and the effects can be sizable .
therefore , higgs physics will have an impact on models we discussed here .
we also computed direct production cross sections of new particles motivated by the muon g-2 .
since models discussed here may not be complete , we could not discuss the complete signatures for the models .
however , since the number of signal events are not negligible , there may be possibility to reveal new physics models for muon g-2 at the lhc .
since we do not need the colored particles to explain the anomaly of muon g-2 , the search for new particles which only have ew interactions will be very important .    finally , we comment on the possible completions of models discussed in this paper .
the structure of the model where both right- and left - handed muons have new yukawa type interactions is quite similar to that of the minimal supersymmetric standard model .
when @xmath195 , the scalar @xmath68 in su(2 ) doublet scalar @xmath64 and su(2 ) singlet scalar @xmath9 correspond to left- and right - handed smuons , respectively .
the su(2 ) singlet fermion @xmath10 is a bino , for example .
a part of a model discussed in ref .
@xcite , which is motivated by neutrino mass and dark matter is also similar to the model we discussed here .
therefore , our approach can capture crucial features of new physics models for the anomaly of muon g-2 and hopefully create an interesting way to find the physics beyond the standard model .
_ note added _    after we completed this work , we learned that atlas and cms collaboration at lhc observed a new particle consistent with the sm higgs boson at a mass of about 125 gev  @xcite . since we have taken the higgs boson mass to be 125 gev in our analysis , our conclusion does not change .
this work is supported in part by grant - in - aid for scientific research from the ministry of education , science , sports , and culture ( mext ) , japan , no .
22540273 ( kt ) .
this work is also supported in part by grant - in - aid for nagoya university global coe program , `` quest for fundamental principles in the universe : from particles to the solar system and the cosmos '' , from the mext , japan .
in this paper , passarino - veltman functions  @xcite are defined by @xmath299 \left[(k+p)^2-m_b^2+i \epsilon\right]},\nonumber \\
p^\mu b_1(a , b;p ) & = & 16\pi^2 \mu^{2\epsilon } \int \frac{d^n k}{i(2\pi)^n}\frac{k^\mu } { \left[k^2-m_a^2+i\epsilon\right]\left[(k+p)^2-m_b^2+i\epsilon\right ] } , \nonumber \\
p^\mu p^\nu b_{21}(a , b;p)&+&g^{\mu \nu}b_{22}(a , b;p ) \nonumber \\ & = & 16\pi^2 \mu^{2\epsilon } \int \frac{d^n k}{i(2\pi)^n } \frac{k^\mu k^\nu } { \left[k^2-m_a^2+i\epsilon\right]\left[(k+p)^2-m_b^2+i\epsilon\right ] } , \ ] ] @xmath300[(k+p_1)^2-m_b^2+i\epsilon][(k+p_1+p_2)^2-m_c^2+i\epsilon ] } , \nonumber   \\ & & \left ( p_1^\mu c_{11}+p_2^\mu c_{12 } \right)(a , b , c;p_1,p_2 ) \nonumber \\ & = & 16\pi^2 \mu^{2\epsilon } \int \frac{d^n k}{i(2\pi)^n}\frac{k^\mu } { [ k^2-m_a^2+i\epsilon][(k+p_1)^2-m_b^2+i\epsilon][(k+p_1+p_2)^2-m_c^2+i\epsilon ] } , \nonumber   \\ & & \left\{(p_1^\mu p_1^\nu c_{21 } + p_2^\mu p_2^\nu c_{22}+(p_1^\mu p_2^\nu+p_1^\nu p_2^\mu)c_{23 } + g^{\mu \nu}c_{24 } \right\ } ( a , b , c;p_1,p_2 ) \nonumber \\ & = & 16\pi^2 \mu^{2\epsilon } \int \frac{d^n k}{i(2\pi)^n}\frac{k^\mu k^\nu } { [ k^2-m_a^2+i\epsilon][(k+p_1)^2-m_b^2+i\epsilon][(k+p_1+p_2)^2-m_c^2+i\epsilon ] } , \nonumber   \\\end{aligned}\ ] ] where we use dimensional regularization in @xmath301 dimensions .
g.  aad _ et al . _ [ atlas collaboration ] , phys .
b * 710 * , 49 ( 2012 ) [ arxiv:1202.1408 ] ; s.  chatrchyan _ et al .
_ [ cms collaboration ] , arxiv:1202.1488 .
see also recent tevatron result ; [ tevnph ( tevatron new phenomina and higgs working group ) and cdf and d0 collaborations ] , arxiv:1203.3774 .
k.  nakamura _ et al .
_ [ particle data group collaboration ] , j.  phys .
g g * 37 * , 075021 ( 2010 ) ; b.  l.  roberts , chin .
c * 34 * , 741 ( 2010 ) [ arxiv:1001.2898 ] .
k.  hagiwara _ et al .
_ , j.  phys .
g g * 38 * , 085003 ( 2011 ) [ arxiv:1105.3149 ] .
see also f.  jegerlehner and a.  nyffeler , phys .
rept .   * 477 * , 1 ( 2009 ) [ arxiv:0902.3360 ] ; m.  davier _ et al .
_ , eur .
j.  c * 71 * , 1515 ( 2011 ) [ erratum - ibid .
c * 72 * , 1874 ( 2012 ) ] [ arxiv:1010.4180 ] ; f.  jegerlehner and r.  szafron , eur .
j.  c * 71 * , 1632 ( 2011 ) [ arxiv:1101.2872 ] ; t.  teubner _
et al . _ , aip conf .
proc .   * 1343 * , 340 ( 2011 ) ; nucl .
phys .  proc .
* 218 * , 225 ( 2011 ) .
t.  aoyama _ et al .
_ , arxiv:1205.5370 ; arxiv:1205.5368 ; r.  boughezal and k.  melnikov , phys .
b * 704 * , 193 ( 2011 ) [ arxiv:1104.4510 ] ; a.  e.  dorokhov , a.  e.  radzhabov and a.  s.  zhevlakov , eur .
j.  c * 71 * , 1702 ( 2011 ) [ arxiv:1103.2042 ] ; t.  goecke , c.  s.  fischer and r.  williams , phys .
d * 83 * , 094006 ( 2011 ) [ arxiv:1012.3886 ] ; eur .
j.  a * 47 * , 28 ( 2011 ) [ arxiv:1009.5297 ] ; a.  nyffeler , chin .  phys .
c * 34 * , 705 ( 2010 ) [ arxiv:1001.3970 ] ; a.  nyffeler , phys .
d * 79 * , 073012 ( 2009 ) [ arxiv:0901.1172 ] ; d.  k.  hong and d.  kim , phys .
b * 680 * , 480 ( 2009 ) [ arxiv:0904.4042 ] ; j.  prades , e.  de rafael and a.  vainshtein , [ arxiv:0901.0306 ] ; a.  e.  dorokhov and w.  broniowski , phys .
d * 78 * , 073011 ( 2008 ) [ arxiv:0805.0760 ] ; j.  bijnens and j.  prades , mod .
a * 22 * , 767 ( 2007 ) [ hep - ph/0702170 ] ; m.  hayakawa _ et al .
_ , pos lat * 2005 * , 353 ( 2006 ) [ hep - lat/0509016 ] ; k.  melnikov and a.  vainshtein , phys .  rev .
d * 70 * , 113006 ( 2004 ) [ hep - ph/0312226 ] ; j.  h.  kuhn , a.  i.  onishchenko , a.  a.  pivovarov and o.  l.  veretin , phys .
d * 68 * , 033018 ( 2003 ) [ hep - ph/0301151 ] ; t.  kinoshita and m.  nio , phys .  rev .
lett .   * 90 * , 021803 ( 2003 ) [ hep - ph/0210322 ] ; m.  j.  ramsey - musolf and m.  b.  wise , phys .
lett .   * 89 * , 041601 ( 2002 ) [ hep - ph/0201297 ] ; j.  bijnens , e.  pallante and j.  prades , nucl .
b * 626 * , 410 ( 2002 ) [ hep - ph/0112255 ] ; nucl .
b * 474 * , 379 ( 1996 ) [ hep - ph/9511388 ] ; j.  bijnens , e.  pallante and j.  prades , phys .
* 75 * , 1447 ( 1995 ) [ erratum - ibid .
* 75 * , 3781 ( 1995 ) ] [ hep - ph/9505251 ] ; i.  r.  blokland , a.  czarnecki and k.  melnikov , phys .
lett .   * 88 * , 071803 ( 2002 ) [ hep - ph/0112117 ] ; m.  knecht , a.  nyffeler , m.  perrottet and e.  de rafael , phys .  rev .
lett .   * 88 * , 071802 ( 2002 ) [ hep - ph/0111059 ] ; m.  knecht and a.  nyffeler , phys .
d * 65 * , 073034 ( 2002 ) [ hep - ph/0111058 ] ; z.  bern , a.  de freitas , l.  j.  dixon , a.  ghinculov and h.  l.  wong , jhep * 0111 * , 031 ( 2001 ) [ hep - ph/0109079 ] ; m.  hayakawa and t.  kinoshita , phys .
d * 57 * , 465 ( 1998 ) [ erratum - ibid .
d * 66 * , 019902 ( 2002 ) ] [ hep - ph/9708227 ] ; m.  hayakawa , t.  kinoshita and a.  i.  sanda , phys .
d * 54 * , 3137 ( 1996 ) [ hep - ph/9601310 ] ; m.  hayakawa , t.  kinoshita and a.  i.  sanda , phys .
lett .   * 75 * , 790 ( 1995 ) [ hep - ph/9503463 ] .
t.  moroi , phys .  rev .
d * 53 * , 6565 ( 1996 ) [ erratum - ibid .
d * 56 * , 4424 ( 1997 ) ] [ hep - ph/9512396 ] .
g.  c.  mclaughlin and j.  n.  ng , phys .
b * 493 * , 88 ( 2000 ) [ hep - ph/0008209 ] .
s.  r.  choudhury , a.  s.  cornell , a.  deandrea , n.  gaur and a.  goyal , phys .
d * 75 * , 055011 ( 2007 ) [ hep - ph/0612327 ] .
k.  kannike , m.  raidal , d.  m.  straub and a.  strumia , jhep * 1202 * , 106 ( 2012 ) [ arxiv:1111.2551 ] ; h.  davoudiasl , h.  -s .
lee and w.  j.  marciano , arxiv:1205.2709 .
g.  t hooft and m.  j.  g.  veltman , nucl .
b * 153 * , 365 ( 1979 ) ; g.  passarino and m.  j.  g.  veltman , nucl .
b * 160 * , 151 ( 1979 ) .
g.  -c .
cho _ et al .
_ , jhep * 1111 * , 068 ( 2011 ) [ arxiv:1104.1769 ] .
g.  -c .
cho and k.  hagiwara , nucl .
b * 574 * , 623 ( 2000 ) [ hep - ph/9912260 ] .
k.  hagiwara _ et al .
_ , z.  phys .
c * 64 * , 559 ( 1994 ) [ erratum - ibid .
c * 68 * , 352 ( 1995 ) ] [ hep - ph/9409380 ] .
m.  e.  peskin and t.  takeuchi , phys .
lett .   * 65 * , 964 ( 1990 ) ; phys .
d * 46 * , 381 ( 1992 ) .
g.  altarelli and r.  barbieri , phys .
b * 253 * , 161 ( 1991 ) ; g.  altarelli , r.  barbieri and s.  jadach , nucl .
b * 369 * , 3 ( 1992 ) [ erratum - ibid .
b * 376 * , 444 ( 1992 ) ] ; g.  altarelli , r.  barbieri and f.  caravaglios , phys .
b * 349 * , 145 ( 1995 ) .
i.  maksymyk , c.  p.  burgess and d.  london , phys .
d * 50 * , 529 ( 1994 ) [ hep - ph/9306267 ] .
tevatron electroweak working group , f.  t.  c.  collaboration and d.  collaboration , arxiv:1204.0042 .
tevatron electroweak working group and for the cdf and d0 collaborations , arxiv:1107.5255 .
see web page : http://lepewwg.web.cern.ch/lepewwg/    j.  alwall _ et al .
_ , jhep * 0709 * , 028 ( 2007 ) [ arxiv:0706.2334 ] .
for example , see g.  aad _ et al .
_ [ atlas collaboration ] , arxiv:1206.1760 ; s.  chatrchyan _ et al . _
[ cms collaboration ] , arxiv:1206.3949 .
t.  hambye , k.  kannike , e.  ma and m.  raidal , phys .
d * 75 * , 095003 ( 2007 ) [ hep - ph/0609228 ] .
g.  aad _ et al . _ [ atlas collaboration ] , [ arxiv:1207.7214 ] ; s.  chatrchyan _ et al . _
[ cms collaboration ] , phys .
b [ arxiv:1207.7235 ] .