electronic mach - zehnder interferometer ( mzi ) @xcite can be realized with the edge states of the quantum hall liquids ( qhls ) .
together with the quantum antidots @xcite , mzis in the regime of the fractional quantum hall effect ( fqhe ) are expected @xcite to be useful for observation of the fractional statistics of fqhe quasiparticles .
in contrast to fractional quasiparticle charge , which has been confirmed in several experiments @xcite , there is no commonly accepted observation of anyonic statistics of the quasiparticles , which remains a challenging experimental problem .
currently , this problem attracts interest in the context of solid - state quantum computation , since individual manipulation of anyonic quasiparticles involving their braiding provides an interesting possible basis for implementation of the quantum information processing @xcite .
however , in typical interferometer - based experimental set - up the quasiparticles emerge as a continuum of gapless edge excitations that should be described by a 1d field theory @xcite .
individual quasiparticles can be realized in this theory only asymptotically in a special limit . in the fractional edge - states mzi ,
such a limit occurs at _ large _ voltages , when the system is characterized by the hamiltonian dual to the hamiltonian of the initial electron tunneling model of the mzi .
the latter is perturbative in electron tunneling at low voltages , and is much better defined , since weak electron tunneling is probably the most basic process in solid - state physics .    in the dual description derived from the electronic model by the instanton technique @xcite , mzi acquires @xmath7 different quantum states which differ by the the effective flux @xmath8 through it .
this flux contains , in addition to the external flux @xmath9 , a statistical contribution earlier also found @xcite for the antidots .
the tunneling of each quasiparticle changes this effective flux by @xmath10 and therefore switches the mzi from one flux state into another . since quasiparticle also carries the charge @xmath5 , the change of flux by @xmath10 results in the change of the interference phase for the quasiparticles by @xmath11 and the corresponding change of the rate of coherent tunneling through the interferometer .
summation over the @xmath2 flux states in calculation of the physical quantities restores @xmath12 periodicity of their @xmath9 dependence . at low voltages ,
this periodicity is guaranteed by @xmath12 periodicity of the electron tunneling amplitude .    in this work ,
we follow the approach to the mzi that allows us to obtain the uniform description of its transport properties in both regimes of electron and quasiparticle tunneling .
we consider the standard mzi geometry ( fig .
1 ) with two tunneling contacts between the two effectively parallel edges of qhls , but allow for , in general , different filling factors , @xmath0 , with @xmath1 , of these edges . in the symmetric case of equal propagation times along the two edges between the contacts , the corresponding 1d field theory permits an exact bethe ansatz solution @xcite . making use of this solution
, we calculate the zero - temperature full counting statistics @xcite of the charge transferred between the two edges forming mzi .
the transferred charge distribution is shown to reflect the anyonic braiding statistics of the tunneling quasiparticles , and decomposition of electrons into quasiparticles with increase of voltage .
( 3.,1.4 ) ( 0.6,-0.05)=1.9 in    the paper is organized as follows . in section 2
we introduce the model of symmetric mzi and its bethe ansatz solution . in section 3
, we use this solution to derive the general expression for the generating function @xmath13 of the distribution of the transferred charge .
the method employed in this calculation generalizes to the mzi transport through two point contacts the method developed earlier by saleur and weiss @xcite for a single point contact .
the charge transfer statistics found for one contact @xcite demonstrated fractionalization of electron charges with increasing bias voltage across the contact . in the case of mzi , one can expect that in addition to charge fractionalization , braiding properties of anyonic quasiparticles should emerge in the charge transfer statistics . to see this ,
we analyze our result for the transfer statistics in different regimes .
we find that the logarithm of the generating function separates into two parts , which for very different absolute values of the electron tunneling amplitudes in the two contacts , @xmath14 , become identical with the cumulant - generating functions of the two separate contacts . in general , however , each `` single - contact '' term accounts also for the interference between the contacts . at low voltages ,
these single - contact terms are combined in such a way that the charge transfer occurs only in electron units , with the low - voltage asymptotic describing the regular poisson distribution of tunneling electrons . in section 4
, we develop analogous qualitative interpretation of the large - voltage asymptotics of the bethe ansatz result for the generating function .
the interpretation is based on the @xmath2-state model of the quasiparticle tunneling described above . following the method of ref .
, the cumulant - generating function is calculated from the kinetic equation governing the quasiparticle transitions in the basis of the @xmath2 flux states .
direct comparison of this generating function with the leading asymptotics of the bethe ansatz result reveals their coincidence under a special choice of @xmath2 quasiparticle interference phases which are found as @xmath15 , where @xmath16 numbers the flux states of the interferometer .
the @xmath17-dependent , statistical , part of these phases agrees with the expected anyonic statistics of the quasiparticles , while the common phase is given by the electron interference phase @xmath18 with an additional phase shift @xmath19 for even @xmath2 .
the equivalence of the two distributions in the large - voltage limit proves that the bethe ansatz construction we implemented indeed describes the statistical transmutation of the effective flux through the mzi . in section 5 ,
we use the obtained generating function to find the first three cumulants of the transferred charge distribution .
we study , in particular , the voltage dependence of the first and the second cumulants proportional to the average current and the shot noise , respectively .
the ratio of the two , the fano factor , is an experimentally observable @xcite characteristics of the electron - quasiparticles decomposition .
finally , we consider the third cumulant , which determines the asymmetry of the transferred charge distribution around average . in section 6 , we look at the special case @xmath20 which allows to obtain these cumulants in terms of elementary functions for arbitrary voltages .
the results of this work are summarized in the conclusion .
we start our discussion with the electronic model of mzi ( fig .
1 ) formed by two single - mode edges with filling factors @xmath21 , @xmath22 .
electron operator @xmath23 of the edge @xmath17 is expressed using the standard bosonization approach @xcite as @xmath24 } .\ ] ] here @xmath25 are the two chiral bosonic modes propagating in the same direction ( to the right in fig .  1 ) , which satisfy the usual equal - time commutation relations
@xmath26=i \pi \mbox{sgn}(x)\ , \delta_{lp}$ ] .
the majorana fermions @xmath27 account for the mutual statistics of electrons in different edges , and @xmath28 is a common energy cut - off of the edge modes .
the fermi momenta @xmath29 correspond to the average electron density in the edges , while the operators of the density fluctuations are : @xmath30 .    in the symmetric case of equal times of excitation propagation between the contacts along the two edges ,
the two combinations of the bosonic operators @xmath25 that enter the electron tunneling terms of the two contacts can be expressed as the values at points @xmath31 of the same right - propagating chiral bosonic field : @xmath32 the langrangian describing electron tunneling in the two contacts can then be written as : @xmath33 \ , , \label{e2}\ ] ] where @xmath34 and @xmath35 are the absolute values and the phases of the dimensionless tunneling amplitudes .
the products of the majorana fermions @xmath36 were omitted from the lagrangian ( [ e2 ] ) , since they cancel out in each perturbative order due to charge conservation .
the phases @xmath35 include contributions from the external magnetic flux @xmath9 and from the average electron numbers @xmath37 on the two sides of the interferometer : @xmath38 + \mbox{const } \equiv -\kappa .\ ] ] the factor @xmath39 in the lagrangian ( [ e2 ] ) follows from the normalization of the bosonic field @xmath40 , which in the absence of tunneling is a free right - propagating chiral field .
this field undergoes successive scattering at the two contacts by the tunneling terms of the lagrangian .
the scattering breaks the charge conservation and therefore creates tunneling current .
the applied voltage can be introduced at first as a shift of the incoming field of one of the edges : @xmath41 . as one can see from eqs .
( [ phi- ] ) , such a shift translates into the shift of the tunneling field @xmath42 .
the thermodynamic bethe ansatz solution of the tunneling model was developed @xcite for a single - point tunneling contact with @xmath43 by application of one - particle boundary @xmath44-matrices @xcite to a distribution of the bosonic field excitations ( kinks , antikinks , and breathers ) introduced through the massless limit of the `` bulk '' sine - gordon model .
this solution was generalized @xcite to the two tunneling contacts relevant for the mzi problem by successive application of two boundary @xmath44-matrices to the same distribution of the excitations of the model . for the charge transport , only the kink - antikink ( and vice verse ) transitions are important , and their boundary @xmath44-matrices are written as @xmath45 here the standard energy scales @xmath46 , @xmath47 , are used to characterize the tunneling strength at the individual contacts , and @xmath48 ) \over \nu \gamma(\nu/[2(1-\nu ) ] ) } .\ ] ] the explicit relation between the energy scales @xmath46 and electron tunneling amplitudes is given below .
at zero temperature , dynamics of the liquid should be described with only one type of quasiparticles , e.g. , kinks , which fill out all available states with the `` bulk '' distribution @xmath49 , with @xmath50 being the quasiparticle momentum , up to some limiting momentum @xmath51 defined by the applied voltage .
each momentum-@xmath50 quasiparticle undergoes successive scattering at the two tunneling contacts independently of other quasiparticles .
the overall scattering process is described by the product of the two boundary @xmath44-matrices given by eq .
( [ e27 ] ) . our goal is to find the cumulant - generating function @xmath52 of the charge transfer between the two branches of the interferometer , which is defined , as usual , as the logarithm of the fourier transform of the probability distribution function of the transferred charge .
we measure the charge in units of the elementary electron charge by setting @xmath53 .
the independence of the scattering events of quasiparticles with different momenta implies then that @xmath52 can be found as a sum of logarithm of the generating functions of individual momentum states , and its long - time asymptotics is @xmath54 the generating function of one state with momentum @xmath50 @xmath55 is defined by the total transition probability @xmath56 of the momentum-@xmath50 kink into antikink . taking the product of the scattering matrices of the two contacts to find the total scattering matrix @xmath57 , and using the parametrization of the contact tunneling strengths as @xmath58
, we can write the probability @xmath56 from the corresponding matrix element of @xmath57 as @xmath59 here @xmath60 are the transition probabilities in the individual point contacts : @xmath61^{-1 } , \label{tau1}\ ] ] and the factor @xmath62 characterizes interference between the two contacts : @xmath63 without a loss of generality , we assume below a specific relation between the tunneling strength parameters of the two contacts , @xmath64 , and write them as @xmath65 , with @xmath66 .    the aim of our subsequent derivation is to find the cumulant - generating function @xmath52 in terms of the two generating functions @xmath67 for charge transfer in an individual point contact that was found from the bethe ansatz solution by saleur and weiss @xcite . this derivation does not need the explicit expressions for @xmath49 and @xmath51 which can be found in @xcite . following the approach for one contact , we first relate @xmath52 in eq .
( [ lnp1 ] ) to the effective tunneling current . to do this ,
we introduce the generalized tunneling probability @xmath68 : @xmath69^{-1 } \ , .
\label{tauc2}\ ] ] which is the solution of the following differential equation in the new parameter @xmath70 : @xmath71 that satisfies the initial condition @xmath72 .
one can extend eq .
( [ p1 ] ) for @xmath73 to include the parameter @xmath70 through the substitution @xmath74 .
equation ( [ tauc2 ] ) shows then that the logarithm of @xmath73 extended this way can be expressed as @xmath75\big{|}_{z=0}^{z = i\xi } \ , .
\label{der1}\ ] ] calculating the derivatives of eq .
( [ der1 ] ) with respect to @xmath70 and @xmath76 , and using eqs .
( [ tauc2 ] ) , one can see that @xmath77 combining eq .
( [ der2 ] ) with ( [ lnp1 ] ) in which @xmath73 is extended to include the parameter @xmath70 , one sees that the cumulant - generating function satisfies the following relation @xmath78 which expresses it through the total tunneling current @xmath79 in the two contacts that is defined by the generalized tunneling probability @xmath80 .    as the next step , one substitutes eq .
( [ tauc1 ] ) into ( [ tauc2 ] ) and casts the total tunneling probability @xmath80 into the following form @xmath81 where @xmath82 is defined by the conditions that @xmath83 and @xmath84 .
differentiation of eq .
( [ theta1 ] ) shows that the coefficient in eq .
( [ tauc4 ] ) in front of the sum can be written as the derivative of @xmath85 : @xmath86 using eqs .
( [ tauc4 ] ) and ( [ theta2 ] ) in the definition of the tunneling current @xmath79 ( [ lnp2 ] ) we obtain an important relation @xmath87 which expresses the derivative of the cumulant - generating function ( [ lnp2 ] ) for the charge transfer in a symmetric interferometer in terms of the tunneling current @xmath88 in one point contact between the two edges with the effective filling factor @xmath89 .
the current in one contact has been calculated @xcite from the bethe ansatz solution , and its tunneling conductance @xmath90=g_{1/m}[(v / a ) e^{-\theta/[2(m-1)]}]=i_{1/m}(\theta , v ) /v\ ] ] is given at zero temperature by a universal scaling function expressed in the form of the low- and high - voltage expansion series .
integrating eq .
( [ ic ] ) over @xmath70 , and using the result in eq .
( [ lnp2 ] ) , we express the generating function @xmath91 in the following form : @xmath92 } ] \ , .\label{lnp3 } \end{aligned}\ ] ] the explicit expansion series for @xmath93 in eq .
( [ lnp3 ] ) allow integration in each order .
the integration transforms the generating function ( [ lnp3 ] ) into the sum of the two generating functions @xmath67 for charge transfer in individual contacts and gives our key result : @xmath94 the single - contact generating function @xmath67 is known in terms of the low- and high - voltage expansion series @xcite : @xmath95 where @xmath96 and @xmath97 is the conductance quantum .
equation ( [ lnps1 ] ) representing the total generating function @xmath98 as the sum of the two single - contact generating functions @xmath67 seems to suggest that the charge transfer in the mzi is divided into two independent processes associated with the two point contacts of the mzi .
indeed , one manifestation of this division is that dependence of each of these processes on the bias voltage @xmath99 ( through the function @xmath100 ) is determined , as in eq .
( [ lnps2 ] ) , only by the characteristics energy scale @xmath46 of the corresponding junction .
such a division , however , is not complete .
the total generating function @xmath98 depends also on the charge dynamics in the interferometer as a whole , since each individual charge transfer process in one contact triggers multiple charge transfers involving interference between both contacts of the interferometer .
information about this charge dynamics enters eq .
( [ lnps1 ] ) through the function @xmath101 determined by eq .
( [ theta1 ] ) and sensitive to both interferometer contacts .
nevertheless , the interference can become irrelevant if the two contacts are strongly asymmetric .
for @xmath102 , the generating function @xmath98 is well approximated at low voltages by the single - contact @xmath103 defined by @xmath104 , i.e. by the strongest electron tunneling amplitude @xmath105  see eq .
( [ w1 ] ) below . at large voltages , the generating function @xmath98 is approximated by the single - contact @xmath67 defined by @xmath106 , i.e. , by the weakest electron tunneling amplitude @xmath107 .      at small voltages , when @xmath108 for both @xmath47 , the low - voltage expansions of the single - contact generating function @xmath67 ( [ lnps2 ] ) for both terms in eq .
( [ lnps1 ] ) can be combined as follows : @xmath109 \ , .\label{lnp4}\end{aligned}\ ] ] since @xmath110 is a polynomial of @xmath111 , and the latter , according to eq .
( [ theta1 ] ) , is a linear function of @xmath112 , this expansion of the mzi generating function shows that at low voltages , the charge transfer between the edges of the mzi is quantized in units of electron charge .    more explicitly , using the standard expansion of @xmath113 ( see eq .
1.331.4 in ref .  ) and the relations that follow from eqs .
( [ b1 ] ) and ( [ theta1 ] ) : @xmath114 \ , , \nonumber\\ r\equiv b\tanh\delta\theta_0=   { |t_{1b}^{m-1}+t_{2b}^{m-1}e^{i\kappa}|^2 \over t_{1b}^{2(m-1)}+t_{2b}^{2(m-1 ) } } \ , , \label{r } \\
\cosh\delta\theta_0= { t_{1b}^{2(m-1)}+t_{2b}^{2(m-1 ) } \over 2(t_{1b}t_{2b})^{(m-1 ) } } , \nonumber\end{aligned}\ ] ] we bring eq .
( [ lnp4 ] ) into the following form : @xmath115^n \nonumber \\ \cdot \ , \big\ { [ 1+r(z-1)]^n   + n \sum_{l=1}^{[n/2 ] } \frac{(-1)^lc^{n - l-1}_{l-1}}{l(2\cosh\delta\theta_0)^{2l } } \nonumber \\ \cdot \ , \big [ 1+r(z-1 ) \big]^{n-2l } \big\ } \big|^{z = e^{i\xi}}_{z=1 } .
\label{lnp5 } \end{aligned}\ ] ] equation ( [ lnp5 ] ) quantifies our previous conclusions about properties of the mzi charge transfer statistics ( [ lnps1 ] ) in the low - voltage regime . in the limit of strongly different contacts , @xmath116
, one finds that @xmath117 and @xmath118 , so that the charge transfer statistics ( [ lnp5 ] ) approaches that of one point contact @xcite characterized by @xmath104 ( i.e. , the strongest electron tunneling amplitude @xmath105 ) with corrections in @xmath119 also quantized in the electron charge units .
the @xmath120th term in the expansion of this statistics in powers of bias voltage @xmath99 corresponds in this case to tunneling of exactly @xmath120 electrons .
by contrast , the @xmath120th order term of the general mzi transfer statistics ( [ lnp5 ] ) involves transfer of all numbers of electrons up to @xmath120 . in the lowest order in @xmath99 ,
the mzi statistics reduces to the poisson distribution , with the coefficient in front of @xmath121 in eq .
( [ lnp5 ] ) is equal to the average electron tunneling current .
one can check this starting from eq .
( [ e2 ] ) by direct perturbative calculation @xcite , if the energy scales @xmath46 are expressed through the electron tunneling amplitudes @xmath34 : @xmath122 ^{{1/(m- 1 ) } } .
\label{w1}\ ] ]      at large voltages , @xmath123 , the combination large - voltage expansions of both single - contact generating functions @xmath67 in eq .
( [ lnps1 ] ) brings the total mzi generating function into the following form : @xmath124 \ , . \;\;\ ; \label{lnp6}\end{aligned}\ ] ] in terms of the parameters introduced in eq .  ( [ r ] )
, this equation can be rewritten as : @xmath125^{{n\over m } } \nonumber \\ \cdot\ ,   \sum_\pm \big [ 1+r ( z-1 ) \pm \big ( [ 1+r(z-1)]^{2 } \nonumber \\
- \cosh^{-2}\delta\theta_0 \big)^{{1 \over 2 } } \big]^{{n\over m } } \big|^{z = e^{i\xi}}_{z=1 } . \;\;\ ; \label{lnp7}\end{aligned}\ ] ] we can see again that in the asymmetric limit , @xmath126 , when @xmath117 and @xmath118 , eq .
( [ lnp7 ] ) for the charge transfer statistics reduces to that of a single point contact @xcite . the dominant contact is now characterized by the larger quasiparticle tunneling amplitude @xmath127 ( i.e. , smaller electron tunneling amplitude @xmath107 ) and corresponding energy scale @xmath106 , related to @xmath127 as : @xmath128^{m/(m-1)}. \label{w2}\ ] ] the @xmath120th order term in the expansion ( [ lnp7 ] ) of the generating function corresponds to the transfer of the fractional charge @xmath129 by @xmath120 quasiparticles
one can see , however , that the mzi transfer statistics ( [ lnp7 ] ) does not contain in general the terms @xmath130 that would correspond directly to transfer of individual quasiparticles of charge @xmath131 . in particular , the @xmath132 term of the expansion that gives the leading large - voltage contribution to the statistics , can not be interpreted as a poisson process of tunneling of independent quasiparticles , in contrast to the leading low - voltage term that did represent poisson process of individual electron tunneling events .
the reason for this is the @xmath2-state dynamics of the effective flux through the interferometer associated with the quasiparticle tunneling , which introduces correlations in the tunneling process .
these correlations can be most easily understood in the description of the quasiparticle tunneling based on kinetic equation .
such an equation is discussed in the next section .
to derive kinetic equation that reproduces the large - voltage asymptotics of the generating function ( [ lnp6 ] ) , and therefore provides a simple physical picture of the dynamics of quasiparticle tunneling , we start by rewriting this asymptotics in terms of the quasiparticle tunneling amplitudes @xmath133 . from the last relation in eq .
( [ r ] ) , we have @xmath134 . \label{der5}\ ] ] using this equation to transform the leading , @xmath132 , term in eq .
( [ lnp6 ] ) , and replacing the energy scales @xmath46 by the quasiparticle amplitudes @xmath133 with the help of eq .
( [ w2 ] ) , we have : @xmath135   \ , , \label{lnp8}\ ] ] where @xmath136 .    kinetic equation describing the quasiparticle tunneling can be written down based on the following considerations .
the general picture of the quasiparticle dynamics in the mzi discussed in the introduction implies that the quasiparticles create statistical contribution to the effective flux through the interferometer . because of this statistical contribution
, the mzi can be found in @xmath2 separate states which differ by the effective flux , with each tunneling quasiparticle changing successively the state @xmath17 into @xmath137 modulo @xmath2 .
since the total rates of the quasiparticle tunneling in the mzi depend on the effective flux , such dynamics of flux introduces correlations into quasiparticle transitions , separating naturally the processes of successive quasiparticle tunneling events into the groups of @xmath2 transitions . as usual , to cast the kinetic equation governing this flux dynamics into the form appropriate for the calculation of the cumulant - generating function , we multiply transition probabilities by a factor @xmath138 that keeps track of the transferred charge .
then , we introduce the probabilities @xmath139 that at time @xmath140 the mzi is in the state @xmath17 and @xmath120 quasiparticles have been transferred through it . combining these probabilities into an @xmath2-dimensional vector @xmath141 , one can write the kinetic equation in the following matrix form : @xmath142 according to the qualitative picture of quasiparticle tunneling discussed above
, the transition matrix has a simple form , with the only non - vanishing elements are those on the main diagonal , @xmath143 , and those with @xmath144 : @xmath145 here the kronecker symbol @xmath146 is defined modulo @xmath2 .
the leading large - time asymtotics of the generating function for the probability distribution evolving according to the kinetic equation ( [ kineq ] ) is ( see , e.g. , @xcite ) : @xmath147 , where @xmath148 is the maximum eigenvalue of the transition matrix ( [ m ] ) .
the structure of this matrix shows directly that the characteristic equation @xmath149 has the form : @xmath150 the maximum eigenvalue @xmath148 is the solution of this equation which goes to zero when @xmath151 , since all other eigenvalues of the matrix @xmath152 are negative .    before trying to establish the general relation between the generating function obtained from this equation and the generating function ( [ lnp8 ] )
, we consider the simple case @xmath20 . in this case , eq .  ( [ lnp8 ] ) for the large - voltage asymptotics of the generating function ( [ lnp6 ] ) , can be simplified further .
first , we have from eq .  ( [ r ] ) : @xmath153^{1/2 } .\ ] ] using this relation , eq .  ( [ r ] ) , and eq .  ( [ der5 ] ) with @xmath20
, we transform eq .
( [ lnp8 ] ) into : @xmath154^{1/2}\!\!\ ! -\sum_j w^2_j\big\ } .
\label{lnp9}\end{aligned}\ ] ] on the other hand , in the kinetic - equation approach , eq .
( [ lambda1 ] ) be solved readily for @xmath20 giving the following expression for @xmath148 : @xmath155^{1/2 } -\sum_j \gamma_j \big)/2 \ , .
\label{lambda2}\ ] ] this equation describes the generating function for the statistics of any transfer process consisting of the two steps with the rates @xmath156 , e.g. , incoherent charge transfer through a resonant level @xcite .
comparison of eqs .
( [ lambda2 ] ) and ( [ lnp9 ] ) shows that the generating function ( [ lambda2 ] ) obtained from the kinetic equation reproduces the large - voltage asymptotics ( [ lnp9 ] ) of the generating function ( [ lnp6 ] ) , if the tunneling rates are taken as @xmath157 equation ( [ emp1 ] ) for the tunneling rates agrees with the physical picture of quasiparticle tunneling discussed above . statistical contribution to the effective flux though
the mzi means that tunneling of each quasiparticle changes the phase between the interferometer contacts by @xmath158 in agreement with the quasiparticle anyonic exchange statistics .
equation ( [ emp1 ] ) also shows that the quasiparticles see a phase shift @xmath159 in addition to the phase @xmath160 induced by the external magnetic flux .
the origin of this phase shift is discussed below .    following this logic
, we look for solution of eq .
( [ lambda1 ] ) in the case of arbitrary @xmath2 taking the tunneling rates @xmath161 as : @xmath162 with some unknown @xmath163 . with these tunneling rates , eq .
( [ lambda1 ] ) reads @xmath164 where @xmath165 . to further transform eq .
( [ lambda3 ] ) , we use the basic identity : @xmath166 and the two identities that follow directly from it : @xmath167 = \cosh\delta\theta   \label{iden } \\ - \ , ( -)^m \cos(m\phi ) \ , .\nonumber\end{aligned}\ ] ] the first one is obtained essentially by taking @xmath168 in eq .
( [ identity ] ) , while the second one follows from the first if the sum of the cosines is transformed into their product . by direct comparison of the second identity in ( [ iden ] ) with eq .
( [ lambda3 ] ) we see that @xmath148 as defined by @xmath169 in eq .
( [ lnp8 ] ) indeed solves eq .
( [ lambda1 ] ) if @xmath170 making use of eq .
( [ identity ] ) one more time , we calculate the product on the right - hand - side of eq .
( [ rr ] ) : @xmath171 one can see that with this expression for the product , eq .
( [ rr ] ) precisely coincides with the definition of @xmath172 by eq .
( [ r ] ) . indeed , combining all three relations in eq .
( [ r ] ) one can express @xmath172 as @xmath173 replacing @xmath46 in eq .
( [ theta3 ] ) with the amplitudes @xmath133 through eq .
( [ w2 ] ) , we see that eq .
( [ theta3 ] ) precisely coincides with eqs .
( [ rr ] ) and ( [ productgammas ] ) , if @xmath174 , i.e. , if the phase @xmath163 is chosen to satisfy the condition @xmath175 the deviation of the phase @xmath163 from @xmath176 in this equation is important only for even @xmath2 . in this case
, it produces the shift @xmath177 of the interference phase in the quasiparticle tunneling rates from the value induced by the external magnetic field .
this shift coincides with the phase acquired by one of the two klein factors of the mzi quasiparticle tunneling action of the mzi @xcite in the flux - diagonal representation due to the @xmath2-power periodicity condition for the klein factors , which corresponds physically to the requirement of the proper exchange statistics between electrons and quasiparticles .
this phase shift ensures , actually , that there is no shift in the interference pattern of the tunnel current in the interferometer ( see the final results below ) between the regimes of electron and quasiparticle tunneling .
so far , we have established the interpretation of the low- and high - voltage asymptotic behavior of the charge transfer statistics in terms of , respectively , tunneling of individual electrons and quasiparticles . in this section ,
we calculate the charge transfer cumulants in these two limits , with the emphasis on the quasiparticle limit which exhibits the non - trivial behavior of the cumulants .
we also will use the generating function found in sec .
[ lnp ] for arbitrary voltages to calculate the full voltage dependence of the cumulants , and to study the crossover between the two asymptotic regimes of electron and quasiparticle tunneling .
the cumulants of the charge @xmath178 transferred through the interferometer during a large time interval @xmath140 can be found from the cumulant - generating function ( [ lnp2 ] ) by the standard relation : @xmath179 where the current @xmath79 is given by eq .
( [ ic ] ) .
the first cumulant gives the average tunneling current : @xmath180 . at arbitrary voltages ,
the average current was calculated before in @xcite .
its large-@xmath99 asymptotics can also be obtained directly from eqs .
( [ lnp7 ] ) and ( [ w2 ] ) : @xmath181 where the interference factor @xmath62 ( [ b1 ] ) can be expressed in term of the quasiparticle amplitudes @xmath182 as @xmath183 the second equality in eq .
( [ i1 ] ) provides direct interpretation of the asymptotics of the average current in terms of the quasiparticles transitions @xcite .
it can be proven formally by making use of the identity @xmath184 \
, , \label{sum1}\ ] ] which can be obtained by directly differentiating individual rates @xmath161 in this equation . on the other hand , differentiating the product of all @xmath161s together , as given by eq .
( [ productgammas ] ) , and using eq .
( [ fi ] ) , we obtained the second equality in ( [ i1 ] ) .
this equality agrees naturally with the simple solution of the quasiparticle kinetic equation , which gives the average tunneling rate as inverse of the average tunneling times in the different states of the interferometer .
the second cumulant @xmath185 defines the spectral density of the current fluctuations at zero frequency @xmath186 , which at zero temperature reflects the shot noise associated with the charge transfer processes . from eq .
( [ ic ] ) , one can write the first derivative of the current as @xmath187 \nonumber \\ -(\partial_u\delta\theta(u))^2 \sum_\pm \partial_{\bar \theta } i_{1/m}(\bar \theta\mp \delta\theta(u),v)\ , .
\label{dic } \end{aligned}\ ] ] substituting this formula into eq .
( [ cum ] ) and calculating the derivatives of @xmath188 ( [ theta2 ] ) for @xmath189 we obtain expression for the spectral density of current : @xmath190 it is convenient to characterize the short noise represented by this spectral density through the fano factor @xmath191 defined as @xmath192 .
in the case of mzi , the fano factor reflects both the charge and statistics of the tunneling excitations and illustrates the transition between the electron and quasiparticle regimes . in the low - voltage limit , @xmath193 as a result of the regular poisson process of electron tunneling . to find the fano factor in the quasiparticle , large - voltage , limit , we start with eq .
( [ i2 ] ) which gives the following general expression for @xmath191 : @xmath194^{-1 }   \big\ } \ , . \;\;\ ; \label{ff } \end{aligned}\ ] ] using eq .
( [ r ] ) , the fact that in the large - voltage limit only one quasiparticle tunneling term @xmath195 can be kept in the current @xmath88 , cf .
( [ i1 ] ) , and that with the parametrization of the energy scales @xmath46 with @xmath196 introduced above , @xmath197 , we get from eq .
( [ ff ] ) : @xmath198 the fano factor ( [ f ] ) corresponds to the dynamics of quasiparticle tunneling as described by the kinetic equation ( [ kineq ] ) .
this can be seen by following the steps similar to that taken above for the average current . applying the differential operator from eq .
( [ sum1 ] ) to individual terms in the sum of the inverse tunneling rates @xmath161 , one obtains directly the following identity : @xmath199 on the other hand , replacing the sum of inverse @xmath161s in this equation with the corresponding expression from eq .
( [ i1 ] ) : @xmath200 \ , , \ ] ] and performing differentiation , we see that the large - voltage asymptotics ( [ f ] ) of the fano factor can be written in terms of the tunneling rates @xmath161 as @xmath201 this result agrees with the calculation @xcite based directly on the kinetic equation . because of the complex nature of the quasiparticle tunneling dynamics characterized by @xmath2 different tunneling rates @xmath161 , @xmath191 is not equal simply to the quasiparticle charge @xmath131 but varies as a function of parameters , e.g. the interference phase @xmath18 , between @xmath131 and @xmath202 .
( 3.,2.6 ) ( -0.2,-0.15)=3.4 in    at arbitrary bias voltage @xmath99 , the fano factor @xmath191 should be plotted numerically .
figure 2 shows @xmath191 in the case @xmath203 which corresponds , e.g. , to tunneling between the two @xmath204 edges .
the curves are shown for different degrees of asymmetry between the two contacts of the interferometer and two values of interference phase , maximum constructive interference , @xmath205 , and complete destructive interference , @xmath206 .
the range of variation of @xmath191 with the interference phase @xmath18 decreases with increasing junction asymmetry . in general ,
the curves show the transition between electron tunneling with @xmath193 at small voltages @xmath99 to quasiparticle tunneling at large voltages . in the quasiparticle regime
, @xmath191 is still can be significantly different from @xmath207 because of the non - trivial three - state flux dynamics of the mzi .
in particular , for identical junctions , the three total quasiparticle tunneling rates ( [ emp2 ] ) satisfy the relation : @xmath208 . taking into account eq .
( [ fi ] ) for the non - statistical contribution to the interference phase , one can see that under the condition of destructive interference , @xmath209 , the tunneling rate in one of the flux states of the interferometer , @xmath210 , is much smaller than the rates in the two other states .
this mean that on the relevant large time scale set by the smallest rate , the three quasiparticles transition that transfer interferometer from state @xmath210 back to itself happen almost simultaneously , so that the three quasiparticle charges @xmath207 are effectively transferred together , restoring @xmath191 back to 1 .    finally , we study the third charge transfer cumulant that characterizes the asymmetry around average of the transferred charge distribution , and has been measured experimentally for electron tunneling in metallic tunnel junctions  see , e.g. , @xcite , and in quantum point contacts @xcite . as for the other cumulants ,
the large - time asympotic of the third cumulant is linear in time , and it can be characterized by the coefficient @xmath211 . to calculate this coefficient , we first find the second derivative of the tunnel current from eq .
( [ dic ] ) : @xmath212 ^ 2 \big)i - \nonumber \\
\sum_\pm \big [ 3\partial_u\delta\theta   \partial^2_u\delta\theta \partial_{\bar \theta } \mp ( \partial_u\delta\theta)^3 \partial^2_{\bar \theta } \big ] i_{1/m}(\bar \theta\mp \delta\theta , v ) .
\;\;\ ; \label{ddic } \end{aligned}\ ] ] derivatives of @xmath213 here can be found from eq .
( [ theta2 ] ) .
in particular , the coefficient in front of @xmath214 in eq .  ( [ ddic ] ) can be expressed as : @xmath215 ^ 2=1- 3\coth \delta \theta(u)\nonumber \\ \times \partial_u \delta\theta(u ) + [ 3\coth^2 \delta \theta(u)-1 ] ( \partial_u\delta\theta(u))^2\ , .
\label{coeff } \end{aligned}\ ] ] substitution of eqs .
( [ ddic ] ) , ( [ coeff ] ) , and ( [ theta2 ] ) into eq .
( [ cum ] ) gives us the coefficient @xmath216 of the third cumulant : @xmath217 i \nonumber \\   - \sum_{j=1,2 } \big[3 b^2 ( 1-b\coth \delta \theta_0 ) \partial_\theta i_{1/m } ( \theta_j , v ) \;\;\ ; \label{c3 - 2 } \\
+ ( -)^jb^3\partial^2_\theta i_{1/m}(\theta_j , v ) \big ] .
\nonumber   \end{aligned}\ ] ] the ratio @xmath218 has also been suggested @xcite as a possible alternative to the fano factor to characterize the charge of the tunneling particles . indeed , in a poisson process , @xmath219 is equal to the fano factor multiplied by the tunneling charge
. therefore , in the case of mzi , @xmath219 ( [ c3 - 2 ] ) reduces to 1 in the low - voltage limit , as a result of the regular poisson electron tunneling . in the quasiparticle , large - voltage , limit , repeating the calculation similar to that leading to eq .
( [ ff ] ) , we can relate both factors as follows : @xmath220
for general @xmath2 , the results for the cumulants of the charge transfer statistics in the mzi discussed above can not be presented in a finite analytical form for arbitrary bias voltages .
the situation is simpler for @xmath20 , when the kink - quasiparticles of the `` bulk '' sine - gordon model that provide the basis for the bethe - ansatz solution of the mzi transport are the regular fermions ( though carrying charge @xmath221 ) , and their distribution @xmath49 in eq .
( [ lnp1 ] ) is the fermi - dirac step - function . in practice
, the @xmath20 regime should take place in the mzi formed by the two edges of different quantum hall liquids , with filling factors @xmath204 and @xmath222 .
individual tunneling contacts of this type have been realized experimentally @xcite .
using the fermi - dirac property of the distribution function @xmath49 for @xmath20 , one can obtain the cumulant generating function @xmath223 either by integration in eq .
( [ lnp1 ] ) , or equivalently , by direct substitution into eq .
( [ lnp3 ] ) of the known single - contact tunnel conductance @xmath224 , which can be expressed simply as @xmath225/2 $ ] .
the resulting generating function is : @xmath226 \nonumber \\ + ( 1/2 ) \ln[1+y^2_j(u ) ] \big ) \big|^{u = i\xi}_{u=0 } \ , , \label{m2 } \end{aligned}\ ] ] where @xmath227 , and @xmath228 .
this generating function can be combined with eq .
( [ cum ] ) to calculate the cumulants of the transferred change distribution for @xmath20 using the same steps as in the previous section
. alternatively , one can substitute eq .
( [ lnp1 ] ) into ( [ cum ] ) and perform the integration directly . below we briefly discuss the first three cumulants obtained in this way . the first cumulant gives the average tunneling current in the mzi as @xcite : @xmath229 where @xmath230 , and @xmath231 are the characteristic quasiparticle tunneling rates in separate contacts .
using the fact that @xmath232 for @xmath20 , one can see explicitly that the current agrees at large voltages with eq .
( [ i1 ] ) that follows from the quasiparticle kinetic equation with the tunneling rates ( [ emp1 ] ) : @xmath233 it is interesting to note that this agreement relies strongly on the shift of the interference phase ( [ fi ] ) from the externally - induced phase @xmath18 .
( 3.,2.6 ) ( -0.2,-0.15)=3.4 in    the second cumulant gives the following expression for the fano factor at arbitrary voltages , including the transition region between electron and quasiparticle tunneling : @xmath234 \big\ } \ , .
\label{fm2}\end{aligned}\ ] ] this equation is plotted in fig .
( [ fonehalf ] ) and describes the transition from @xmath193 for electron tunneling at small voltages to @xmath235 for quasiparticle tunneling at large voltages .
one can see that the quasiparticle charge @xmath236 manifests itself most clearly for @xmath205 , when the total quasiparticle tunneling rates ( [ emp1 ] ) coincide , @xmath237 , regardless of the relation between the individual rates @xmath238 .
similarly to the case @xmath203 illustrated in fig .  2 , the fano factor reduces to electron value 1 even in the quasiparticle regime , if @xmath239 and @xmath206 . in this case ,
one of the total quasipartical tunneling rates @xmath240 is much smaller than the other , so on the relevant large time scale set by the smaller rate the quasiparticles effectively tunnel together , restoring @xmath191 to 1 .
( 3.,2.6 ) ( -0.2,-0.15)=3.4 in    calculating the third cumulant , we find the following expression for the `` alternative '' fano factor : @xmath241 ^ 2 } \big\ } \nonumber \ , , \end{aligned}\ ] ] where @xmath242 is defined by @xmath243 , @xmath244 .
equation ( [ last ] ) is plotted in fig .
[ f3onehalf ] , which shows @xmath219 as a function of the interference phase @xmath18 at several voltages .
voltage dependence of @xmath219 is qualitatively very similar to that of the fano factor shown in fig .
( [ fonehalf ] ) : it approaches 1 at small voltages , in agreement with the underlying poisson tunneling process of electrons . at large voltages , eq .
( [ last ] ) reduces to the following form @xmath245 this expression can be understood in terms of the same two - state tunneling dynamics of the quasiparticles that was discussed above . also similarly to the fano factor
, the voltage dependence of @xmath219 is non - monotonic , with a minimum between the regimes of electron and quasiparticle tunneling .
the main qualitative difference between the noise - related fano factor and its third - cumulant alternative is that the minimum of @xmath219 can become negative for some values of parameters ( regime not shown in fig .
[ f3onehalf ] ) .
starting from the exact solution of the tunneling model of symmetric mach - zender interferometer in the fqhe regime , we have calculated the statistics of the charge transfer between interferometer edges .
the obtained statistics shows the transition from electron tunneling at low voltages to tunneling of anyonic quasiparticles of the fractional charge @xmath5 and statistical angle @xmath19 at large voltages .
deep in the electron tunneling regime , the dynamics of charge transfer is represented by the standard poisson process .
dynamics of quasiparticle tunneling is more complicated and reflects the existence of @xmath2 effective flux states of the interferometer .
the interference phase between the quasiparticle tunneling amplitudes in two contacts of the interferometer contains a contribution from the quasiparticle exchange statistics , making the quasiparticle tunneling rates in different interferometer states different . in general , the transition from electron to quasiparticle tunneling is reflected in the fano factor @xmath191 or its third - cumulant alternative @xmath219 , which both reach minima in the transition region .
however , in the regime close to complete destructing interference ( interferometer phase @xmath206 and equal tunneling strength in the two contacts ) , both @xmath191 and @xmath219 have electron value 1 for all voltages .


the atlas and the cms collaborations have reported observations of a new particle with a mass of about 126 gev , which is considered to be the standard model ( sm ) higgs boson  @xcite . if the particle is indeed the higgs boson , the sm inevitably involves the hierarchy problem .
this unnaturalness indicates that there lies physics beyond the sm .
there is another indication for physics beyond the sm .
the precise measurement of the muon anomalous magnetic moment ( the muon @xmath1 )  @xcite has shown discrepancy from the sm prediction  @xcite . with dedicated efforts to determine hadronic contributions precisely , the latest result is @xmath2 where the hadronic vacuum polarization is quoted from ref .
@xcite , and the hadronic light - by - light contribution is from ref .
the muon @xmath1 anomaly indicates physics beyond the sm at more than @xmath0 level .
moreover , the difference is as large as the sm electroweak contribution , @xmath3  @xcite .
if new physics is responsible for the discrepancy , its contribution is naively estimated as @xmath4 , where @xmath5 is a coupling constant of new particles to the muon , and @xmath6 is a typical scale of their masses .
thus , new physics around the tev scale is required to involve strong couplings with the muon in order to solve the muon @xmath1 anomaly .
the supersymmetry ( susy ) is a promising candidate for the tev - scale new physics .
the minimal susy standard model ( mssm ) not only solves the hierarchy problem but also yields sizable contributions to the muon @xmath1  @xcite .
the susy contributions to the muon @xmath1 are naturally enhanced by @xmath7 . as we will see in sec .
[ sec : muon_g-2 ] , the muon @xmath1 anomaly is solved if the superparticles ( smuons , neutralinos , and charginos ) are around @xmath8 for @xmath9 .
superparticles are searched for at the lhc .
since none of them has been discovered , colored superparticles are considered to be heavier than @xmath10  @xcite . moreover
, the higgs boson mass of @xmath11 indicates the scalar tops to be as heavy as @xmath12  @xcite and/or to have a large trilinear coupling to the up - type higgs boson  @xcite .
these two results naively contradict with the indication of the muon @xmath1 anomaly that superparticles have a mass of @xmath8 .
in fact , some of the representative susy - breaking mediation mechanisms such as the minimal supergravity and the gauge mediation can not satisfy the above three phenomenological requirements simultaneously , ( i ) providing the higgs boson mass of @xmath11 , ( ii ) avoiding the constraints from the direct searches at the lhc , and ( iii ) solving the muon @xmath1 anomaly ( see ref .
@xcite for a study on the minimal supergravity ) . from the viewpoint of model building , the muon @xmath1 is reconciled with the higgs boson mass of @xmath11 by extending the mssm so that extra contributions to the higgs potential appear  @xcite .
the above inconsistency implies a split mass spectrum of the superparticles as @xmath13 where the squarks are favored to be relatively heavy to explain the higgs boson mass and to satisfy the lhc bounds , while the sleptons , neutralinos and charginos are to be light to explain the muon @xmath1 discrepancy .
this scenario forces us to change strategy of the susy search .
the searches based on productions of colored superparticles , which have been the standard , are not always promising , because the particles can be away from the lhc reach without diminishing a virtue of the susy as a solution to the muon @xmath1 anomaly . in this letter
, we investigate the lhc signatures assuming the mass hierarchy in eq .
( [ eq : rough_spectrum ] ) , with particular attention to the muon @xmath1 .
we examine the direct productions of the superparticles which are relevant for the muon @xmath1 .
the atlas and cms collaborations respectively reported results of their searches for this production channel  @xcite .
we will apply the result by the atlas to the above scenario , and will see that these searches are particularly important in the parameter regions where the muon @xmath1 is explained .    in sec .
[ sec : muon_g-2 ] , the susy contributions to the muon @xmath1 is reviewed , and the parameter regions in which the discrepancy is explained are clarified .
the detailed mass spectra and parameter spaces are introduced in sec .
[ sec : mass ] , and relevant susy searches at the lhc are summarized in sec .
[ sec : lhc ] . in sec .
[ sec : result ] , current lhc bounds on the muon @xmath1 parameter spaces are shown , and future prospects are discussed in sec .  [ sec : future ] . the last section is devoted to the conclusion .
the susy contributions to the muon @xmath1 are dominated by the chargino  sneutrino and the neutralino  smuon loop diagrams . at the leading order of @xmath14 and @xmath15 , where @xmath16 represents susy - breaking masses and the higgsino mass @xmath17 , they are evaluated as  @xcite @xmath18 where @xmath19 is the muon mass , while @xmath20 and @xmath21 are the fine structure constants of the sm u(1)@xmath22 and the su(2)@xmath23 gauge symmetries , respectively .
the loop functions are defined as and @xmath24 , are reduced from the functions , @xmath25 and @xmath26 , in ref .  @xcite . ] @xmath27\ , , \\ & f_n(x , y)= xy \left [ \frac{-3+x+y+xy}{(x-1)^2(y-1)^2 } + \frac{2x\log x}{(x - y)(x-1)^3 } -\frac{2y\log y}{(x - y)(y-1)^3 } \right]\,.\end{aligned}\ ] ] they satisfy @xmath28 and are monochromatically increasing for @xmath29 and @xmath30 . in the limit of degenerate masses ,
they satisfy @xmath31 and @xmath32 .
the arguments in the left - hand side of eqs .
([eq : bmulr ] ) show the superparticles which propagate in each loop diagram .
if one of them decouples , the corresponding susy contribution is suppressed .
( [ eq : whsnu ] ) comes from the chargino
sneutrino diagrams , and eqs .
( [ eq : whmul])([eq : bmulr ] ) are the neutralino  smuon contributions .
we can easily see that the parameters @xmath33 are relevant , where @xmath34 are the bino and wino masses , and @xmath35 are masses of the left- and right - handed smuons .
note that the sneutrino mass @xmath36 is related to @xmath37 by the su(2)@xmath23 symmetry as @xmath38 .
numerically , the susy contributions are evaluated as @xmath39    the susy contributions to the muon @xmath1 are enhanced when @xmath15 is large and @xmath16 is small . for @xmath40 and @xmath41 , they become @xmath42 , which can explain eq .  .
since the experimental value is larger than the sm prediction , @xmath43 is favored in most of the parameter region . throughout this letter
, the gaugino masses are taken positive .
let us investigate the susy contributions in detail to specify the parameter regions which should be searched for at the lhc .
we will study two representative cases when the muon @xmath1 discrepancy is explained .
other possibilities will be discussed in sec .
[ sec : future ] .
the first case is the one when the chargino  sneutrino contribution ( [ eq : whlnu_n ] ) dominates the susy contributions . as can be seen in eqs .
( [ eq : whlnu_n])([eq : blr_n ] ) , this is valid when the relevant mass parameters of eq .
are nearly degenerate .
in fact , in wide susy models which explain the anomaly with non - decoupling higgsinos , the chargino  sneutrino contribution ( [ eq : whlnu_n ] ) dominates .
as the higgsino mass @xmath17 increases , the chargino  sneutrino contribution decreases , while the neutralino
smuon contribution becomes relevant .
this is because the pure - bino contribution eq .
is enhanced under the presence of a large @xmath17 resulting in a large left
right mixing in the smuon mass matrix .
the other contributions in eqs .
( [ eq : whsnu])([eq : bhmur ] ) are suppressed , where the higgsino propagates in the diagrams . in this situation , which we will take as the second case , the muon @xmath1 discrepancy is explained when both @xmath44 and @xmath45 are @xmath8 as will be seen in the next section .    in this letter , we will study the susy searches at the lhc for the parameter regions of these two cases .
the soft susy - breaking parameters are set at the susy scale , and we will not specify the models to realize them .
it is interesting to mention that such parameter regions are realized in wide susy models which solve the muon @xmath1 anomaly  ( see e.g.  @xcite ) , in which much ingenuity is exercised in enhancing the higgs boson mass . as we take model independent approach , the following lhc analysis will be independent of the mechanisms which realize the higgs boson mass of @xmath11 .
let us specify the mass spectrum to study lhc phenomena .
we consider two representative cases in light of the muon @xmath1 : ( i ) the case when the chargino  sneutrino contribution to the muon @xmath1 dominates , and ( ii ) that with a large @xmath17 parameter where the neutralino  smuon diagram is relevant .
the former contribution is controlled by @xmath46 , and the latter is by @xmath47 , according to eqs .   and , respectively . in particular , the @xmath17 parameter is favored to be relatively small in the former case , while it is fairly large in the latter .
the leading susy contributions to the muon @xmath1 are proportional to @xmath15 , while the lhc sensitivities in the following sections do not depend much on it . mainly through the left - right mixing of the stau mass matrix .
@xcite [ footnote : stau ] ] it is fixed as @xmath48 throughout the lhc analysis . as @xmath15 increases , although the susy contributions to the muon @xmath1 are enhanced , the bottom and tau yukawa coupling constants are likely to blow up below the gut scale . if @xmath15 is lowered , the superparticle masses must be reduced to solve the muon @xmath1 anomaly .
then , the lhc searches become easier .
the soft masses of the left- and right - handed selectrons are respectively set to be the same as those of the smuons , which are denoted by @xmath49 and @xmath50 .
the results in the following sections are almost independent of the selectron masses .
the third - generation sleptons ( staus and tau - sneutrino ) are irrelevant for the susy contributions to the muon @xmath1 .
they are taken to be decoupled in the following analysis to simplify the lhc studies .
scenarios with light staus will be discussed in sec .
[ sec : future ] .    for the gaugino masses ,
the following two scenarios are analyzed : ( a ) they satisfy an approximate gut relation , @xmath51 , which is realized in a wide variety of susy models , and ( b ) the gluino is decoupled , while the bino and wino masses satisfy the same relation as ( a ) , @xmath52 .
the susy contributions to the muon @xmath1 are sensitive to @xmath53 and @xmath54 , while @xmath55 barely affects the muon @xmath1 and the higgs boson mass , but the lhc study depends on it . in the case
( a ) , superparticles can be searched for by analyzing events with hard - jets plus missing transverse energy , which come from the gluino pair - production , while electroweak productions of superparticles are relevant in the both cases ( a ) and ( b ) .
we will study both of ( a ) and ( b ) separately for the cases , ( i ) and ( ii ) .
the bino and wino masses are set to satisfy @xmath56 throughout the analysis , and general cases will be discussed in sec .
[ sec : future ] .
it is assumed that the squarks are decoupled from the lhc sensitivities in the analysis .
we set their soft masses as @xmath57 .
they satisfy the current bounds on the superparticle masses obtained at the lhc  @xcite , and also , heavy stops are preferred to realize the higgs boson mass of @xmath11  @xcite .
the following results do not change as long as the masses are large enough .
note that lighter squarks just tighten the lhc constraints , because their productions contribute to the susy events .
the remaining susy parameters are the cp - odd higgs mass , @xmath58 , and scalar trilinear couplings ( @xmath59-terms ) .
they are relevant neither for the muon @xmath1 nor the susy searches at the lhc . in the analysis , they are fixed as @xmath60 and @xmath61 , as a reference .      in the above parameter space , the lightest superparticle ( lsp )
is either the bino - like neutralino , the neutral higgsino , or the sneutrino .
the lsp is a candidate of the dark matter if the r - parity is conserved . in the last set of , @xmath66
is set to be slightly larger than @xmath67 in order to avoid a charged lsp in the whole parameter region , since the charged dark matter is excluded  @xcite .
the sneutrino lsp is also strongly disfavored by the direct searches  @xcite . in the analysis
, we will focus on the cosmologically favored regions .
note that other superparticles such as the gravitino can be lighter than them , and such cases will be discussed in sec .  [
sec : future ] .
in order to solve the muon @xmath1 anomaly , the bino and wino masses are as small as @xmath8 . then , if we respect the approximate gut relation , the gluino has a mass less than a few tev , which is within the reach of the lhc .
the gluinos are produced at the @xmath68 collisions , and decay into neutralinos or charginos with hard jets , @xmath69 similarly to the minimal supergravity models , searches for hard jets plus missing transverse energy are one of the most promising channels .
currently , the atlas collaboration yields the most stringent bound in this category  @xcite with analyzing their data corresponding to @xmath70 obtained at @xmath71 .
we interpret this result in this letter , which will be referred to as `` j - search . ''
the electroweak gaugino productions ( e.g. , @xmath72 ) are also searched for at the lhc .
when sleptons are relatively light , the gauginos can produce hard leptons as @xmath73 as the muon @xmath1 prefers such lighter sleptons , the multi - lepton signatures are another promising channel .
the atlas collaboration recently reported searches for events with three leptons plus missing transverse energy in the data of @xmath74 at @xmath71  @xcite .
we interpret the result , which will be referred to as `` l - search . ''
the multi - lepton signatures are particularly important , because they are yielded by the superparticles relevant for the muon @xmath1 , i.e. , sleptons , charginos , and neutralinos , and appear even in the decoupled gluino scenario ( b ) .
those two atlas results are interpreted in the analysis with monte carlo simulation . on each model point , the mass spectra generated with softsusy3.4  @xcite are passed to susy - hit1.3  @xcite to calculate decay tables of the superparticles . here , the soft susy - breaking parameters are set at the susy scale , @xmath75 .
the kinematical distributions of susy events are simulated by pythia6.4  @xcite with the atlas mc09 tune  @xcite .
the parton distribution functions ( pdfs ) are obtained from the cteq6l1 set  @xcite . as the squarks are assumed to be decoupled ,
superparticles are produced via the gluino pair - production ( @xmath76 ) and the electroweak productions of neutralinos , charginos and sleptons ( @xmath77 ) .
for the gluino productions , the cross sections are normalized by nlo @xmath78-factors obtained with prospino2  @xcite , where the cteq6l1 and cteq6.6 m pdfs  @xcite are used . for the electroweak channels ,
the normalization factor is set as @xmath79 , which is a typical value in the parameter space studied in the next section .
in addition , the susy contribution to the muon @xmath1 is evaluated by feynhiggs  @xcite , which includes two - loop contributions .
detector simulation is employed with delphes2.0  @xcite .
we assume trigger efficiencies to be @xmath80 for simplicity .
as for calorimeter configuration and resolutions , delphes original parameter card for the atlas detector is used .    following the atlas analysis ,
jet reconstruction is done with the anti-@xmath81 jet clustering algorithm  @xcite with the distance parameter @xmath82 .
fastjet  @xcite is utilized .
we also simulate @xmath83-tagging , used as @xmath83-veto in the l - search , with simplified efficiency and fake rates estimated from refs .
the missing transverse energy is read from delphes output without modification .
we take lepton detection efficiencies into account , which are important particularly for the l - search . for electrons ,
the atlas collaboration utilizes three identification criteria .
we consult ref .
@xcite , and set the efficiency of the `` medium '' criterion as @xmath84 for @xmath85 and @xmath86 for @xmath87 , and that of the `` tight '' criterion as @xmath88 for @xmath85 and @xmath89 for @xmath87 . @xmath90
dependence of the efficiency is ignored . for muons
, we interpret the results in ref .
@xcite , and set the efficiency ( for @xmath91 ) as @xmath92 for @xmath93 and @xmath94 for @xmath95 .
note that we further impose the lepton isolation criteria in the l - search , and the overlap removal in the j- and l - searches , in the same way as the atlas analysis .
the signal regions ( srs ) are defined to be the same as those in the original atlas analyses .
the details are summarized in the following subsections . throughout the analysis
, we ignore criteria on the primary vertex , origins of lepton tracks , and jet qualities , which are designed to reduce noise backgrounds .
the @xmath96 method is used to derive exclusions for each model point .
the numbers of the signal events in srs are compared to the corresponding upper bounds obtained in the atlas reports .
the analysis procedures are validated by comparing the simulations with the atlas results .
we checked that calculated effective mass distributions of the signal events agree with those of the atlas in the j - search , and those of the missing transverse energy are reproduced for the l - search .
exclusion plots are also simulated by using the upper bounds on the signal events obtained in the atlas analysis , and become consistent with those in the atlas papers .
.definition of srs of the j - search . among the original 12 srs  @xcite ,
the srs which are relevant for our exclusions are shown .
the effective mass @xmath97 ( @xmath98 ) is defined as a scalar sum of @xmath99 and @xmath100 s of the leading @xmath101-jets ( all jets with @xmath102 ) . in the last two rows of c and e ,
the three values mean different srs , called as ` tight ' , ` medium ' , and ` loose ' , respectively .
[ cols="^,^,^,^ " , ]     the l - search is designed to search for electroweak productions of charginos and neutralinos via the s - channel exchange of a virtual sm gauge boson ( e.g. , @xmath103 ) .
the event selections are summarized in table  [ tab:2012 - 154 ] .
events are required to have exactly three hard leptons .
trigger efficiencies are not taken into account for simplicity .
note that @xmath83-jet veto and @xmath104-cut are necessary to suppress events of sm electroweak backgrounds .
the orange ( yellow ) band shows the region where the susy contributions explain the muon @xmath1 discrepancy at the @xmath105 ( @xmath106 ) level . the dark gray regions in ( a ) and ( c ) are excluded by lep searches for the neutralinos and charginos .
the regions left to the blue dotted lines are excluded by the l - search . assuming the approximate gut relation for the gaugino masses , the regions left to the red lines are excluded by the j - search .
these exclusions are at @xmath107 cl , and the theoretical uncertainty of @xmath108 is included in the hatched regions .
the lsp is the lightest neutralino in the regions above the black thick lines , while the sneutrino is lightest below them .
sleptons become lighter than neutralinos below the black dashed or dotted lines ( see the text for details ) .
, title="fig:",scaledwidth=99.0% ] + ( a ) @xmath109 , @xmath110     explanations .
the orange ( yellow ) band shows the region where the susy contributions explain the muon @xmath1 discrepancy at the @xmath105 ( @xmath106 ) level .
the dark gray regions in ( a ) and ( c ) are excluded by lep searches for the neutralinos and charginos .
the regions left to the blue dotted lines are excluded by the l - search . assuming the approximate gut relation for the gaugino masses , the regions left to the red lines are excluded by the j - search .
these exclusions are at @xmath107 cl , and the theoretical uncertainty of @xmath108 is included in the hatched regions .
the lsp is the lightest neutralino in the regions above the black thick lines , while the sneutrino is lightest below them .
sleptons become lighter than neutralinos below the black dashed or dotted lines ( see the text for details ) . ,
title="fig:",scaledwidth=99.0% ] + ( b ) @xmath111 , @xmath110     +     explanations .
the orange ( yellow ) band shows the region where the susy contributions explain the muon @xmath1 discrepancy at the @xmath105 ( @xmath106 ) level . the dark gray regions in ( a ) and ( c ) are excluded by lep searches for the neutralinos and charginos .
the regions left to the blue dotted lines are excluded by the l - search . assuming the approximate gut relation for the gaugino masses , the regions left to the red lines are excluded by the j - search .
these exclusions are at @xmath107 cl , and the theoretical uncertainty of @xmath108 is included in the hatched regions .
the lsp is the lightest neutralino in the regions above the black thick lines , while the sneutrino is lightest below them .
sleptons become lighter than neutralinos below the black dashed or dotted lines ( see the text for details ) .
, title="fig:",scaledwidth=99.0% ] + ( c ) @xmath112 , @xmath110     explanations .
the orange ( yellow ) band shows the region where the susy contributions explain the muon @xmath1 discrepancy at the @xmath105 ( @xmath106 ) level . the dark gray regions in ( a ) and ( c ) are excluded by lep searches for the neutralinos and charginos .
the regions left to the blue dotted lines are excluded by the l - search .
assuming the approximate gut relation for the gaugino masses , the regions left to the red lines are excluded by the j - search .
these exclusions are at @xmath107 cl , and the theoretical uncertainty of @xmath108 is included in the hatched regions .
the lsp is the lightest neutralino in the regions above the black thick lines , while the sneutrino is lightest below them .
sleptons become lighter than neutralinos below the black dashed or dotted lines ( see the text for details ) . , title="fig:",scaledwidth=99.0% ] + ( d ) @xmath113 , @xmath114    the lhc constraints and the muon @xmath1 are shown in fig .
[ fig : result ] . in the orange ( yellow ) regions ,
the muon @xmath1 discrepancy is explained by the susy contributions at the @xmath105 ( @xmath106 ) level .
the chargino
sneutrino contribution dominates in the figures ( a ) , ( b ) and ( c ) , where the parameters are set as @xmath115 and @xmath116 , respectively .
it is found that the contribution is suppressed as the wino mass increases . on the other hand , the neutralino
smuon contribution dominates in the figure ( d ) with @xmath117 .
the muon @xmath1 is explained even with large @xmath54 , while the susy contribution decreases as the left- and right - handed slepton masses increase .
the lhc exclusion limits are displayed in fig .
[ fig : result ] with the red - solid and blue - dotted lines at the 95% confidence level ( cl ) . including 30% errors of pdf and scale uncertainties ,
the limits distribute in the red / blue hatched regions . assuming the approximate gut relation , the regions left to the red lines
are excluded by the j - search , and those to the blue lines are constrained by the l - search .
if the gluino is decoupled , the red lines disappear and only the l - search provides a constraint , because the j - search targets the gluino - pair productions .
the bounds from the l - search with the approximate gut relation are almost identical to those in the gluino decoupling limit .
it is found that the susy searches at the lhc start to constrain the muon @xmath1 ( orange or yellow ) regions .
the muon @xmath1 regions are also constrained by cosmology .
noting that the bino mass is set to be a half of the wino mass , the sneutrino becomes the lsp below the black solid lines in fig .
[ fig : result ] . due to radiative corrections , which are sizable since the soft susy - breaking parameters are set at @xmath118 .
] such parameter regions are disfavored by the dark matter searches  @xcite . on the other hand ,
the lsp in the regions above the lines is the lightest neutralino .
they are constrained by the j- and l - searches at the lhc .
the lhc constraints are understood as follows .
the j - search targets the gluino - pair productions . among the 12 srs defined in ref .
@xcite , those with higher multiplicity of jets are relevant for the exclusion ( see table .
[ tab:2012 - 109 ] ) .
some of the jets originate in the gluino decay into charginos and neutralinos .
others come from decays of the sm bosons that are decay products of the charginos and neutralinos .
also , the processes are often associated with initial / final state radiations . assuming the approximate gut relation , the gluino mass is related to the wino mass , and the red lines are almost independent of the other model parameters in fig .
[ fig : result ] .
the j - search limit @xmath119 means @xmath120 , leading to @xmath121 , where @xmath122 is the gluino pole mass .
this is almost the same as the constraint in the minimal supergravity models with heavy squarks  @xcite .
if the ratio @xmath123 is varied , the exclusion lines shift correspondingly .    the l - search targets the electroweak superparticle productions .
in particular , the wino productions are significant in the parameter regions in fig .
[ fig : result ] .
the search works only when hard leptons are produced at the wino decay .
in the regions below the black dashed lines , where sleptons are lighter than the wino - like chargino / neutralino , the wino can decay into an on - shell slepton with a lepton , and the slepton produces another lepton when it decays .
then , the sr1b provides a tight constraint . above the lines
, the wino mainly decays into the sm bosons as long as the channels are open . in such regions ,
no limit is obtained by the l - search .
the search becomes effective again if the mass splitting between the bino and the wino is small , where the wino decays into three - body final states , including leptons .
since @xmath53 is set to be @xmath124 , the limit is determined by the wino mass , as can be seen from the vertical blue lines in fig .
[ fig : result ] . here , the sr1a gives a limit . on the other hand
, the search becomes less effective if sleptons are degenerate with the lightest neutralino , because the lepton at the slepton decay becomes soft .
the charged slepton is degenerate with the lightest neutralino on the black dotted lines .
note that the regions below the black solid lines are disfavored by cosmology as discussed above .
it is emphasized that the l - search constraints are determined by the superparticles which are relevant for the muon @xmath1 .    in figs .
[ fig : result ] ( a)(c ) , the muon @xmath1 regions are partly excluded by the j- and/or l - searches . below the black dashed lines ,
tight constraints are obtained by the latter . above them ,
the limits are determined by the wino mass , i.e. , the gluino mass for the j - search and the splitting between @xmath53 and @xmath54 for the l - search . comparing the figures to each other
, it is found that the muon @xmath1 regions with larger @xmath17 are easily detected at the lhc . as @xmath17 increases , @xmath54 is required to be smaller to explain the muon @xmath1 .
if the sleptons are lighter than the wino , a wide region of the @xmath105 parameter space of the muon @xmath1 is already excluded in ( a ) and ( b ) . in contrast , almost all the muon @xmath1 regions are allowed by the lhc in ( c ) .
in particular , the lepton productions are relatively suppressed for @xmath125 even below the black dashed line , because the winos can decay into the higgsinos , and the higgsinos rarely produce leptons due to tiny yukawa couplings .
the sensitivity is slightly improved below the black dash - dotted line in ( c ) , where the sleptons are lighter than the second lightest neutralino . since @xmath126
, the second lightest neutralino does not decay into the sm bosons , but can produce the sleptons . here
, the sr1a provides a limit from the higgsino direct productions .    in figs .
[ fig : result ] ( a)(c ) , @xmath50 is set to be @xmath127 .
since the chargino  sneutrino contribution is insensitive to @xmath50 , the muon @xmath1 regions do not change so much even if @xmath66 is varied .
the lhc constraints do not depend much on @xmath50 either .
even when it is as low as @xmath49 , the right - handed selectron and smuon are rarely produced in the cascade decays of the wino and higgsino , because they are singlet under the sm su(2)@xmath23 symmetry and have tiny yukawa interactions .
if the mass is too small , the right - handed slepton becomes the lsp , and thus , excluded by the cosmology .
the lhc constraints in fig .
[ fig : result ] ( d ) are similar to those in ( b ) , where the wino composes the second lightest neutralino . above the black dashed line ,
since the sleptons are limited to be light by the muon @xmath1 , a large part of the @xmath105 region is already excluded by the j - search , while the l - search bound is determined similarly to ( a)(c ) . on the other hand ,
the muon @xmath1 is explained even with large @xmath54 .
such a parameter region is challenging to be searched for at the lhc .
( see also the discussion in the next section . )    in fig .
[ fig : result ] ( d ) , @xmath17 is set to be @xmath128 . as @xmath17 increases
, the muon @xmath1 is enhanced , and thus , larger @xmath53 is allowed . however , the upper limit on the slepton masses is not relaxed so much , as noticed from eq .  .
moreover , if @xmath17 becomes too large , the electroweak vacuum can be destabilized on the smuon
higgs plane . by naively scaling a bound on @xmath17 from the stau - higgs stability condition  @xcite
, the limit is estimated as @xmath129 for @xmath130 and @xmath131 .
in contrast , the lhc constraints do not change , because the higgsinos are decoupled . on the other hand , if @xmath17 is set to be smaller , the upper bound on the wino mass gets tightened , and the lhc detection becomes easier .
the lhc will be upgraded and run at @xmath132 or @xmath133 with the target luminosity of @xmath134 .
further upgrades are proposed to realize @xmath135 or @xmath136 .
there are several studies on future sensitivities of the susy searches . in this section ,
we discuss future searches for the susy parameter regions of the muon @xmath1 .
multi - jet signature : :    multi - jet searches are important for the approximate gut relation on    the gaugino masses . according to the future sensitivity in    ref .
@xcite based on the minimal supergravity models , the gluino    lighter than @xmath137 and    @xmath138 can be discovered at the    @xmath133 lhc with    @xmath139 and    @xmath140 , respectively .
these gluino masses    correspond to @xmath141 and    @xmath142 under the approximate gut relation ,    which means a large parameter regions in fig .
[ fig : result ] can be    covered .
since the sleptons are expected to be light particularly in a    large wino mass region , signatures of multi - jets plus leptons may    provide a better sensitivity than the expectations based on the    minimal supergravity .
multi - lepton signature : :    multi - lepton searches works very well for the two cases : when sleptons    are lighter than charginos and neutralinos , and when the wino mainly    decays into three bodies .
the former case is typically realized in the    regions below the black dashed lines with    @xmath143 as we saw in figs .
[ fig : result ] ( a ) ,    ( b ) , and ( d ) , where the multi - leptons are provided from the cascade    decays of winos . the latter case corresponds to the left - most regions    in ( a)(d ) , where the three - body decay provides the leptons .
it is    important to study future sensitivities of this channel .
@xmath144 signature from gaugino / higgsino decays : :    if charginos and neutralinos are lighter than sleptons and can decay    into the sm bosons , multi - lepton searches are not promising , but    rather , we should explore searches for the sm bosons from the gauginos    and/or higgsinos . for @xmath144 signature from gaugino pair    productions , a study within the minimal supergravity framework in    ref .
@xcite shows that the wino masses of    @xmath145 and @xmath146 can    be covered at the @xmath133 lhc with    @xmath139 and    @xmath140 , respectively , in one lepton plus    multi - bottoms signature . with this kind of searches , most of the muon    @xmath1 regions above the black dashed lines in    fig .
[ fig : result ] are expected to be examined at the lhc , where the    branching ratios of the sm boson channels are similar to those in the    minimal supergravity models .
+    the @xmath144 signature is also expected from higgsino decays    as long as higgsino is heavier than , but not degenerate to , the lsp .
they are produced by the @xmath147-channel exchange of a virtual    electroweak gauge boson , and their production cross sections are    comparable to those of the winos .
thus , the above sensitivity on the    wino mass can be naively applied to the higgsino mass . @xmath148 or @xmath149 signature : :    it may be possible to discover charginos and neutralinos via the    @xmath150 or @xmath151 channels with a large missing    transverse energy . since the bosons can be reconstructed by visible    decay products or have sizable branching ratios of the bottom quark    productions , the events could be discriminated from the sm    backgrounds .
we need further studies to explore such possibilities .
same - sign leptons : :    very recently , the authors in ref .
@xcite studied the events with a    same - sign lepton pair and no other leptons , obtained from same - sign    @xmath152-bosons .
they found that the signature is useful to    search for the electroweak superparticle events . in our setup
, this    signature appears more frequently because the sleptons are lighter ,    and thus , may cover the muon @xmath1 regions .
di - lepton channel search : :    even with combining the above techniques , the bottom right corners of    figs .
[ fig : result ] ( c ) and ( d ) are still challenging .
future    multi - lepton searches at the lhc may not be sensitive enough to cover    the regions completely since @xmath54 is very large . here
,    since sleptons are required to be light , di - lepton channels via the    pair productions of sleptons might be promising .
the search at the lhc    is challenging because it suffers from huge sm backgrounds and their    production cross sections are much smaller than those of the charginos    and the neutralinos in this parameter region .
the ilc is suitable to    search for the di - lepton signatures .
so far , we considered two representative scenarios to explain the muon @xmath1 discrepancy : with dominance of the chargino  sneutrino contribution , and of the pure - bino  smuon diagram .
we have also set that the stau is heavy , and the lsp is within the mssm .
let us touch on other possibilities .
we have imposed @xmath153 , but this is not mandatory .
as long as @xmath53 is smaller than , but not degenerate to , @xmath54 , and the winos are produced at the lhc , the above conclusions and discussions hold valid . on the other hand ,
if @xmath54 is too large with @xmath53 fixed , the muon @xmath1 regions disappear in figs .
[ fig : result ] ( a)(c ) , whereas those in fig .
[ fig : result ] ( d ) remain .
searches for di - lepton signatures are necessary as discussed above . for the case
when @xmath154 , the muon @xmath1 regions in figs .
[ fig : result ] ( a)(c ) does not depend much on @xmath53 , while the lhc signatures are changed .
signatures from the cascade decays of binos are less promising than those of winos due to a smaller cross section of the bino production .
rather , signatures from the productions of higgsinos and sleptons are hopeful . moreover
, a disappearing track at the lhc is expected to be observed .
the track comes from the charged wino , which is almost degenerate with the neutral wino and decays into a neutral wino with a soft pion  @xcite .
so far , we investigated the parameter regions where the susy contribution to the muon @xmath1 is dominated either by the chargino  sneutrino contribution or the pure - bino  smuon contribution .
let us mention other possibilities .
first of all , the contribution of eq .   can not take dominance , since it is always buried in that of eq .  .
next , if the left - handed smuon is decoupled , only eq
.   contributes to the muon @xmath1 , where the sign of @xmath155 is favored to be negative .
lastly , eq .   becomes dominant if the wino and the right - handed smuon are decoupled . however , the contributions of eqs .   and are quantitatively small as shown in eqs .   and , respectively .
thus , the superparticles must be very light , and the scenarios are expected to be examined by searching for the sm boson or di - lepton channels .
the staus and the tau sneutrino were supposed to be decoupled from the lhc phenomenology .
if their masses are comparable to those of the selectron and smuon , physics becomes much involved .
attention must be paid to the vacuum stability conditions on the stau - higgs plane especially when @xmath17 is large , i.e. , in fig .
[ fig : result ] ( d )  @xcite .
the decay branching ratio of the higgs boson into the di - photon may be affected in such parameter spaces ( see the footnote [ footnote : stau ] ) . as for the direct searches for the superparticles ,
since a stau tends to be light due to a large left - right mixing in the stau mass matrix , taus are likely to be produced in the decay chains of the superparticles .
then the tau reconstruction becomes a central subject , which is one of the challenging topics of the lhc . because of the rich phenomenology ,
these scenarios are worth investigated for future .
finally let us relax our restriction of the lsp .
we focused on the neutralino lsp to avoid the cosmological constraints .
however , they can be circumvented if the lsp is a superparticle which is not in the mssm , such as the gravitino and the axino . in those cases ,
charged particles or sneutrino can be the lsp among the mssm superparticles ( mssm - lsp ) , where the lhc signatures change correspondingly .
if a sneutrino is the mssm - lsp , hard leptons are less produced because the left - handed charged sleptons are degenerate to the sneutrinos .
although sneutrino decays into the lsp are not visible , the pair - production of the sleptons could be found by the mono - jet or mono - photon searches .
in addition , this case might be checked by reconstructing soft particles that are produced by the quasi - degenerate charged slepton when it decays into the sneutrino .
this search could be done at the ilc . on the other hand ,
when the mssm - lsp is a charged slepton , it leaves a track in the detectors if it is long - lived , or generates a hard lepton when it decays .
finally , if the lightest neutralino is the mssm - lsp and decays into the lsp in the detectors , the signature can be distinguished from the backgrounds by searching for events with multiple ( displaced ) photons .
all of these possibilities and future sensitivities on them will be studied elsewhere .
the susy not only solves the hierarchy problem between the planck and electroweak scales , but also explain the 34@xmath156 deviation of the muon @xmath1 .
the latter suggests that superparticles have masses of @xmath8 , whereas the higgs boson mass of @xmath11 and the lhc results on the direct susy searches indicate that they exist in @xmath12 . in this letter
, we examined the parameter regions where the muon @xmath1 is explained with squarks decoupled .
the muon @xmath1 anomaly is solved when electroweak superparticles are light and @xmath9 .
they are searched for by the lhc j- and l - searches .
we found that the wino is constrained to be heavier than @xmath157 by the j - search under the approximate gut relation on the gaugino masses , and than @xmath158 by the l - search irrespectively of the gluino mass .
in addition , if sleptons are produced by decays of charginos and neutralinos , the l - search provides a severe constraint , and the bound can be as tight as @xmath159 .
consequently , the muon @xmath1 parameter space starts to be constrained by the lhc independently of details of the susy models .
we also discussed future prospects on the collider searches for the superparticles which are relevant for the muon @xmath1 .
a wide parameter region is expected to be examined at the 13@xmath133 lhc . in particular ,
the sm boson channels of the electroweak superparticle decays such as the @xmath160 productions are promising for the study .
however , the searches are challenging when the higgsinos are degenerate with the lightest neutralino , or if both of the wino and the higgsino are heavy .
we need further studies to explore the whole parameter space where the muon @xmath1 anomaly is solved .
the study in this letter is motivated by the muon @xmath1 anomaly .
the @xmath161 deviation may originate in the uncertainties .
currently they are dominated by those from the measurement and from theoretical calculation of the hadronic contributions .
the former will be improved by the experiments at fermilab  @xcite and j - parc  @xcite .
the reduction of the uncertainty in the hadronic vacuum polarization contribution requires more experimental data for the cross section of @xmath162 , e.g. , at vepp-2000  @xcite .
the evaluations of the hadronic light - by - light contribution have been improving ( see ref .
@xcite for a recent review ) , and lattice calculations are in development  @xcite .
if the anomaly is confirmed in future , the susy is one of the most attractive candidates for the solution .
this letter is a first step towards the lhc searches for the scenario , and we hope it is useful for further studies .
this work was supported by jsps kakenhi grant no .
23740172 ( m.e . ) , no .
21740164 ( k.h . ) , no .
22244021 ( k.h . ) and no .
228132 ( s.i . ) . the work of t.y .
was supported by an advanced leading graduate course for photon science grant .
collaboration , http://dx.doi.org/10.1016/j.physletb.2012.08.020[phys .
( 2012 ) 129 ] .
collaboration , http://dx.doi.org/10.1016/j.physletb.2012.08.021[phys .
( 2012 ) 3061 ] . * muon g-2 * collaboration , http://dx.doi.org/10.1103/physrevd.73.072003[phys
.   rev .   * d73 * ( 2006 ) 072003 ] ; b.
l.  roberts , http://dx.doi.org/10.1088/1674-1137/34/6/021[chin .
* c34 * ( 2010 ) 741744 ] . k.  hagiwara , a.  d.  martin , d.  nomura , and t.  teubner , http://dx.doi.org/10.1016/j.physletb.2007.04.012[phys .   lett .
( 2007 ) 173179 ] ; t.  teubner , k.  hagiwara , r.  liao , a.  d.  martin , and d.  nomura , http://dx.doi.org/10.1088/1674-1137/34/6/019[chin .
* c34 * ( 2010 ) 728734 ] ; k.  hagiwara , r.  liao , a.  d.  martin , d.  nomura , and t.  teubner , http://dx.doi.org/10.1088/0954-3899/38/8/085003[j .
g * g38 * ( 2011 ) 085003 ] .
m.  davier , _ et al . _ , http://dx.doi.org/10.1140/epjc/s10052-009-1219-4[eur .
j. ( 2010 ) 127136 ] ; m.  davier , a.  hoecker , b.  malaescu , c.  z.  yuan , and z.  zhang , http://dx.doi.org/10.1140/epjc/s10052-010-1246-1[eur .
j. ( 2010 ) 19 ] ; m.  davier , a.  hoecker , b.  malaescu , and z.  zhang , http://dx.doi.org/10.1140/epjc/s10052-010-1515-z[eur .
j. ( 2011 ) 1515 ] [ http://arxiv.org/abs/1010.4180[arxiv:1010.4180 ] ] .
erratum http://dx.doi.org/10.1140/epjc/s10052-012-1874-8[ibid . *
c72 * ( 2012 ) 1874 ] .
j.  prades , e.  de  rafael , and a.  vainshtein , http://arxiv.org/abs/0901.0306[arxiv:0901.0306 ] ( 2009 ) .
a.  czarnecki , w.  j.  marciano , and a.  vainshtein , http://dx.doi.org/10.1103/physrevd.67.073006[phys .
* d67 * ( 2003 ) 073006 ] [ http://arxiv.org/abs/hep-ph/0212229[hep-ph/0212229 ] ] .
erratum http://dx.doi.org/10.1103/physrevd.73.119901[ibid .
* d73 * ( 2006 ) 119901 ] .
j.  l.  lopez , d.  v.  nanopoulos , and x.  wang , http://dx.doi.org/10.1103/physrevd.49.366[phys .
* d49 * ( 1994 ) 366372 ] .
u.  chattopadhyay and p.  nath , http://dx.doi.org/10.1103/physrevd.53.1648[phys .
* d53 * ( 1996 ) 16481657 ] . t.  moroi , http://dx.doi.org/10.1103/physrevd.53.6565 ,    10.1103/physrevd.56.4424[phys .
* d53 * ( 1996 ) 65656575 ] .
collaboration , http://cdsweb.cern.ch/record/1472710[atlasconf2012109 ] ( 2012 ) .    *
cms * collaboration , http://dx.doi.org/10.1007/jhep10(2012)018[jhep * 1210 * ( 2012 ) 018 ] ; * cms * collaboration , http://dx.doi.org/10.1103/physrevlett.109.171803[phys .
( 2012 ) 171803 ] .
y.  okada , m.  yamaguchi , and t.  yanagida , http://dx.doi.org/10.1143/ptp.85.1[prog .
* 85 * ( 1991 ) 16 ] ; j.  r.  ellis , g.  ridolfi , and f.  zwirner , http://dx.doi.org/10.1016/0370-2693(91)90863-l[phys .
* b257 * ( 1991 ) 8391 ] ; h.  e.  haber and r.  hempfling , http://dx.doi.org/10.1103/physrevlett.66.1815[phys .   rev .
( 1991 ) 18151818 ] ; j.  r.  ellis , g.  ridolfi , and f.  zwirner , http://dx.doi.org/10.1016/0370-2693(91)90626-2[phys .
* b262 * ( 1991 ) 477484 ] . y.  okada , m.  yamaguchi , and t.  yanagida , http://dx.doi.org/10.1016/0370-2693(91)90642-4[phys .   lett .
( 1991 ) 5458 ] .
m.  endo , k.  hamaguchi , s.  iwamoto , k.  nakayama , and n.  yokozaki , http://dx.doi.org/10.1103/physrevd.85.095006[phys .
* d85 * ( 2012 ) 095006 ] .
m.  endo , k.  hamaguchi , s.  iwamoto , and n.  yokozaki , http://dx.doi.org/10.1103/physrevd.84.075017[phys .
* d84 * ( 2011 ) 075017 ] ; m.  endo , k.  hamaguchi , s.  iwamoto , and n.  yokozaki , http://dx.doi.org/10.1103/physrevd.85.095012[phys .
rev .   * d85 * ( 2012 ) 095012 ] ; m.  endo , k.  hamaguchi , s.  iwamoto , and n.  yokozaki , http://dx.doi.org/10.1007/jhep06(2012)060[jhep
* 1206 * ( 2012 ) 060 ] ; m.  endo , k.  hamaguchi , k.  ishikawa , s.  iwamoto , and n.  yokozaki , http://dx.doi.org/10.1007/jhep01(2013)181[jhep * 1301 * ( 2013 ) 181 ] . j.  l.  evans , m.  ibe , and t.  t.  yanagida , http://dx.doi.org/10.1016/j.physletb.2011.10.031[phys .   lett .
( 2011 ) 342348 ] ; j.  l.  evans , m.  ibe , s.  shirai , and t.  t.  yanagida , http://dx.doi.org/10.1103/physrevd.85.095004[phys .
* d85 * ( 2012 ) 095004 ] .
cho , k.  hagiwara , y.  matsumoto , and d.  nomura , http://dx.doi.org/10.1007/jhep11(2011)068[jhep * 1111 * ( 2011 ) 068 ] .
collaboration , http://cdsweb.cern.ch/record/1493493[atlasconf2012154 ] ( 2012 ) .
m.  carena , i.  low , and c.  e.  m.  wagner , http://dx.doi.org/10.1007/jhep08(2012)060[jhep * 1208 * ( 2012 ) 060 ] .
t.  kitahara , http://dx.doi.org/10.1007/jhep11(2012)021[jhep * 1211 * ( 2012 ) 021 ] .
m.  carena , s.  gori , i.  low , n.  r.  shah , and c.  e.  m.  wagner , http://arxiv.org/abs/1211.6136[arxiv:1211.6136 ] ( 2012 ) .
t.  kitahara and t.  yoshinaga , http://arxiv.org/abs/1303.0461[arxiv:1303.0461 ] ( 2013 ) .
a.  kudo and m.  yamaguchi , http://dx.doi.org/10.1016/s0370-2693(01)00938-8[phys .
( 2001 ) 151155 ] .
t.  falk , k.  a.  olive , and m.  srednicki , http://dx.doi.org/10.1016/0370-2693(94)90639-4[phys .
( 1994 ) 248251 ] .
b.  c.  allanach , http://dx.doi.org/10.1016/s0010-4655(01)00460-x[comput .
commun .
* 143 * ( 2002 ) 305331 ] .
a.  djouadi , m.  m.  muhlleitner , and m.  spira , acta phys .
* b38 * ( 2007 ) 635644 .
t.  sjostrand , s.  mrenna , and p.  z.  skands , http://dx.doi.org/10.1088/1126-6708/2006/05/026[jhep * 05 * ( 2006 ) 026 ] .
collaboration , http://cds.cern.ch/record/1247375/[atlphyspub2010002 ] ( 2010 ) .
w.  beenakker , r.  hopker , m.  spira , and p.  m.  zerwas , http://dx.doi.org/10.1016/s0550-3213(97)00084-9[nucl .
( 1997 ) 51103 ] .
s.  heinemeyer , w.  hollik , and g.  weiglein , http://dx.doi.org/10.1016/s0010-4655(99)00364-1[comput .
commun .
* 124 * ( 2000 ) 7689 ] ; s.  heinemeyer , w.  hollik , and g.  weiglein , http://dx.doi.org/10.1007/s100529900006[eur .
.   j.   * c9 * ( 1999 ) 343366 ] ; g.  degrassi , s.  heinemeyer , w.  hollik , p.  slavich , and g.  weiglein , http://dx.doi.org/10.1140/epjc/s2003-01152-2[eur .   phys .   j. ( 2003 ) 133143 ] ; m.  frank , _ et al . _ , http://dx.doi.org/10.1088/1126-6708/2007/02/047[jhep * 0702 * ( 2007 ) 047 ] .
s.  ovyn , x.  rouby , and v.  lemaitre , http://arxiv.org/abs/0903.2225[arxiv:0903.2225 ] ( 2009 ) .
m.  cacciari , g.  p.  salam , and g.  soyez , http://dx.doi.org/10.1088/1126-6708/2008/04/063[jhep * 0804 * ( 2008 ) 063 ] .
m.  cacciari , g.  p.  salam , and g.  soyez , http://dx.doi.org/10.1140/epjc/s10052-012-1896-2[eur .
j. ( 2012 ) 1896 ] ; m.  cacciari and g.  p.  salam , http://dx.doi.org/10.1016/j.physletb.2006.08.037[phys .
( 2006 ) 5761 ] .
collaboration , https://atlas.web.cern.ch/atlas/groups/physics/confnotes/atlas-conf-2012-040/[atlasconf2012040 ] ( 2012 ) .
collaboration , http://dx.doi.org/10.1140/epjc/s10052-012-1909-1[eur .
j. ( 2012 ) 1909 ] .
collaboration , https://atlas.web.cern.ch/atlas/groups/physics/confnotes/atlas-conf-2011-063/[atlasconf2011063 ] ( 2011 ) .    j.  hisano and s.
sugiyama , http://dx.doi.org/10.1016/j.physletb.2010.12.013[phys .
( 2011 ) 9296 ] .
h.  baer , v.  barger , a.  lessa , and x.  tata , http://dx.doi.org/10.1103/physrevd.86.117701[phys .
* d86 * ( 2012 ) 117701 ] .
h.  baer , _ et al . _
, http://arxiv.org/abs/1302.5816[arxiv:1302.5816 ] ( 2013 ) .
collaboration , http://dx.doi.org/10.1007/jhep01(2013)131[jhep * 1301 * ( 2013 ) 131 ] .
collaboration , http://dx.doi.org/10.1016/j.nuclphysbps.2011.06.038[nucl .
* 218 * ( 2011 ) 237241 ] .
collaboration , http://dx.doi.org/10.1088/1742-6596/295/1/012032[j .
( 2011 ) 012032 ] .
y.  m.  shatunov , _ et al .
_ in _ proceedings of the 7th european particle accelerator conference _
, j.  l.  laclare , _ et al .
_ , eds . , conf .
* c0006262 * ( 2000 ) 439441 .
e.  de  rafael , http://arxiv.org/abs/1210.4705[arxiv:1210.4705 ] ( 2012 ) .
t.  blum , http://dx.doi.org/10.1103/physrevlett.91.052001[phys .  rev .
( 2003 ) 052001 ] ; * qcdsf * collaboration , http://dx.doi.org/10.1016/j.nuclphysb.2004.03.026[nucl .  phys .
( 2004 ) 135164 ] ; c.  aubin and t.  blum , http://dx.doi.org/10.1103/physrevd.75.114502[phys .
* d75 * ( 2007 ) 114502 ] ; x.  feng , k.  jansen , m.  petschlies , and d.  b.  renner , http://dx.doi.org/10.1103/physrevlett.107.081802[phys .
( 2011 ) 081802 ] ; p.  boyle , l.  del  debbio , e.  kerrane , and j.  zanotti , http://dx.doi.org/10.1103/physrevd.85.074504[phys .
* d85 * ( 2012 ) 074504 ] ; m.  della  morte , b.  jager , a.  juttner , and h.  wittig , http://dx.doi.org/10.1007/jhep03(2012)055[jhep * 1203 * ( 2012 ) 055 ] .