the surplus of matter over anti - matter within the universe can only be explained with a thorough understanding of cp violation .
the cp phase in the quark mixing matrix of the standard model , which has been confirmed by b - meson experiments  @xcite , is not sufficient to understand the baryon asymmetry of the universe  @xcite .
however , the minimal supersymmetric standard model ( mssm )  @xcite provides new physical phases that are manifestly cp - sensitive . after absorbing non - physical phases , we chose the complex parameters to be the higgsino mass parameter @xmath2 , the u(1 ) , and su(3 ) gaugino mass parameters @xmath1 , and @xmath4 , and the trilinear scalar coupling parameters @xmath5 of the third generation sfermions @xmath6 .
the corresponding phases violate cp and are generally constrained by experimental bounds on electric dipole moments ( edms )  @xcite . however , these restrictions are strongly model dependent  @xcite , such that additional measurements outside the low energy edm sector are required .
+ many cp observables have been proposed and studied in order to measure cp violation .
total cross sections  @xcite , masses  @xcite , and branching ratios  @xcite , are cp - even quantities .
for a direct evidence of cp violation , however , cp - odd ( t - odd ) observables are required .
examples are rate asymmetries of either branching ratios  @xcite , cross sections  @xcite , or angular distributions  @xcite .
since these rate asymmetries require the presence of absorptive phases , they are typically small , of the order of @xmath7 , if they are not resonantly enhanced  @xcite .
larger cp - odd observables which already appear at tree - level are desirable .
these are t - odd triple products of momenta and/or spins , from which cp - odd asymmetries can be constructed .
such triple product asymmetries are highly cp - sensitive , and have been intensively studied both at lepton and hadron colliders  @xcite .
+ third generation sfermions have a rich phenomenology at high energy colliders like the lhc  @xcite or ilc  @xcite due to a sizable mixing of left and right states .
in addition , the cp phases of the trilinear coupling parameters @xmath5 are rather unconstrained by the edms  @xcite
. the phases of @xmath8 and @xmath9 have been studied in stop  @xcite and sbottom  @xcite decays , respectively . since these are decays of a scalar particle , the spin - spin correlations have to be taken into account .
the triple product asymmetries can then be up to @xmath10 , for sizable squark mixing . similarly for probing the cp - violating phase of @xmath11 in the stau vertex , @xmath12-@xmath13-@xmath14 , it is essential to include the tau spin .
only then is there a sensitivity to the phase of @xmath11  @xcite .
if the spin of the tau is summed over , this crucial information is lost .
triple product asymmetries including the tau polarization have been studied in neutralino decays @xmath15  @xcite , and also in chargino decays @xmath16  @xcite .
it was shown that the normal tau polarization itself is cp - sensitive , and that the asymmetries are large and of the order of @xmath17 to @xmath3 .
+ we are thus motivated to study cp violation , including the tau polarization , in the two - body decay of a stau @xmath18 followed by the subsequent chain of two - body decays    [ eq : decaychain ] @xmath19    see fig .
[ feyn ] for a schematic picture of the entire stau decay .
this process is kinematically open for a mass hierarchy @xmath20 where the staus are heavier than the smuons and selectrons .
we thus work in mssm scenarios with heavier stau soft susy breaking parameters @xmath21 we show that the normal tau polarization , with respect to the plane spanned by the @xmath14 and @xmath22 momentum , is a triple product asymmetry which is sensitive to the phases of @xmath0 , @xmath1 , and @xmath2 in the stau - neutralino sector . for nearly degenerate stau masses ,
@xmath23 a strong stau mixing is obtained which results in tau polarization asymmetries of more than @xmath3 .
this should be measurable at colliders .
since the stau is a scalar particle , its particular production does not contribute to cp - sensitive spin - spin correlations , and can thus be considered separately .
this allows a collider - independent study , where we only discuss the boost dependence of the cp asymmetries .
+    the paper is organized as follows . in section  [ sec : formalism ] we review stau mixing and the stau - neutralino lagrangian with complex couplings .
we calculate the amplitude squared for the entire stau decay in the spin - density matrix formalism  @xcite .
we construct the cp asymmetry from the normal tau polarization , and discuss its mssm parameter dependence , as well its boost dependence for colliders like the ilc and lhc . in section  [ sec : results ] , we numerically study the phase and parameter dependence of the asymmetry , and the stau and neutralino branching ratios .
we comment on the impact of the @xmath24 decay in scenarios with nearly degenerate stau masses .
we summarize and conclude in section  [ sec : summary ] .
the appendices contain the definitions of momenta and spin vectors , the analytical expressions for the stau decay amplitudes in the spin - density matrix formalism , and formulae for the stau decay widths .
in the complex mssm , the stau mixing matrix in the @xmath25-basis is  @xcite @xmath26    cp violation is parameterized by the physical phase @xmath27 , \\ \lambda_{{\tilde{\tau } } } & = & a_\tau - \mu^ * \cot{\beta},\end{aligned}\ ] ] with the complex trilinear scalar coupling parameter @xmath0 , the complex higgisino mass parameter @xmath2 , and @xmath28 , the ratio of the vacuum expectation values of the two neutral higgs fields .
the left and right stau masses are @xmath29 with the real soft susy breaking parameters @xmath30 , the electroweak mixing angle @xmath31 , and the masses of the @xmath32 boson @xmath33 , and of the tau lepton , @xmath34 . + in the mass basis ,
the stop eigenstates are  @xcite @xmath35 with the diagonalization matrix @xmath36 and the stau mixing angle @xmath37    the stau mass eigenvalues are @xmath38 .
\end{aligned}\ ] ]      the relevant lagrangian terms for the stau decay @xmath39 are  @xcite @xmath40 with @xmath41 , and the weak coupling constant @xmath42 , @xmath43 . the couplings are defined as  @xcite @xmath44 the stau diagonalization matrix @xmath45 is given in eq .
( [ eq : decaychain ] ) , and @xmath46     h_{\tau i}^r \end{array}\right ) , \qquad    { \cal b}_i^\tau \equiv   \left(\begin{array}{ccc }      h_{\tau i}^l\\[2 mm ]      f_{\tau i}^r\end{array}\right ) .
\label{eq : ab}\ ] ] in the photino , zino , higgsino basis ( @xmath47 ) , we have @xmath48 , \qquad \label{eq : flt}\\[2 mm ]   f_{\tau i}^r & = & \sqrt{2 } \sin \theta_w                 \left(\tan\theta_w n_{i2}^*-n_{i1}^*\right ) , \label{eq : frt}\\[2 mm ]   h_{\tau i}^l & = & ( h_{\tau i}^r)^{\ast } =               -y_\tau ( n_{i3}^\ast\cos\beta+n_{i4}^\ast\sin\beta ) , \label{eq : hrlt}\\[2 mm ]        y_\tau   & = & \frac{m_\tau}{\sqrt{2}\,m_w \cos\beta } ,   \label{eq : yt}\end{aligned}\ ] ] with @xmath49 the mass of the @xmath50 boson , and @xmath51 the complex , unitary @xmath52 matrix that diagonalizes the neutralino mass matrix  @xcite @xmath53 the interaction lagrangian relevant for the neutralino decay @xmath54 , for @xmath55 is  @xcite @xmath56 with the couplings @xmath57 given in eqs .   and  .
the unnormalized , @xmath58 , hermitian , @xmath14 spin density matrix for stau decay , eqs .   and  , reads @xmath59 with the amplitude @xmath60 , and the lorentz invariant phase space element @xmath61 , for details see appendix  [ sec : phasespace ] .
the @xmath14 helicities are denoted by @xmath62 and @xmath63 . in the spin density matrix formalism
@xcite , the amplitude squared is given by @xmath64 & & \sum_{\lambda_i \lambda_i^\prime } \rho_d(\tilde\tau)^{\lambda_\tau\lambda_\tau^\prime } _
{ \lambda_i\lambda_i^\prime } \;\rho_{d_1}({\tilde{\chi}^0_i})^{\lambda_i^\prime\lambda_i } \;d_2(\tilde\ell ) , \label{rhoampli}\end{aligned}\ ] ] with the neutralino helicities @xmath65 .
the amplitude squared decomposes into the remnants of the propagators @xmath66 with mass @xmath67 , and width @xmath68 of particle @xmath69 or @xmath70 , and the unnormalized spin density matrices for stau decay @xmath71 , and neutralino decay @xmath72 .
the decay matrix of the spinless slepton is a factor since the polarizations of the final lepton and lsp are not accessible .
the corresponding amplitude is denoted by @xmath73 .
defining a set of spin basis vectors @xmath74 for the tau , see eqs .   in appendix
[ sec : momenta ] , and @xmath75 for the neutralino @xcite , the spin density matrices can be expanded in terms of the pauli matrices  @xmath76 @xmath77 \sigma_{\rm d}^b\,\delta^{\lambda_\tau\lambda_\tau^\prime }                     ( \sigma^b)_{\lambda_i\lambda_i^\prime }   + \sigma_{\rm d}^{ab}\,(\sigma^a)^{\lambda_\tau\lambda_\tau^\prime }                         ( \sigma^b)_{\lambda_i\lambda_i^\prime } , \label{rhod}\\[2 mm ] \rho_{{\rm d}_1}({\tilde\chi}_i^0)^{\lambda_i^\prime\lambda_i}= { \rm d}_1\,\delta^{\lambda_i^\prime\lambda_i } + \sigma_{{\rm d}_1}^b\,(\sigma^b)^{\lambda_i^\prime\lambda_i } , \phantom{spacespace } \label{rhod1}\end{aligned}\ ] ] with an implicit sum over @xmath78 , respectively .
the real expansion coefficients @xmath79 , @xmath80 , @xmath81 , @xmath82 , @xmath83 and @xmath84 contain the physical information of the process .
@xmath79 denotes the unpolarized part of the amplitude for stau decay @xmath85 , @xmath86 denotes the unpolarized part for neutralino decay @xmath87 , respectively .
@xmath88 gives the tau polarization , @xmath82 , and @xmath83 describe the contributions from the neutralino polarization , and @xmath89 is the spin - spin correlation term , which contains the cp - sensitive parts .
we give the expansion coefficients explicitly in appendix  [ sec : coefficients ] .
+ inserting the density matrices , eqs .   and , into eq .
, we get for the amplitude squared @xmath90{\rm                             d}_2,\quad   \label{eq : nosum } \end{aligned}\ ] ] with an implicit sum over @xmath78 .
the amplitude squared @xmath91 is now decomposed into an unpolarized part ( first summand ) , and into the part for the tau polarization ( second summand ) , in eq .  .
by using the completeness relations for the neutralino spin vectors , eq .
( [ eq : completeness2 ] ) , the products in eq .
( [ eq : nosum ] ) can be written , followed by @xmath92 .
the signs in parentheses in eqs .
( [ eq : product1 ] ) and ( [ eq : product2 ] ) hold for the charge conjugated stau decay @xmath93 ; @xmath92 . in order to obtain the terms for the decay @xmath94 , however , followed by the neutralino decay into a positively charged slepton , @xmath95 , one has to reverse the signs of eqs .
( [ eq : product1 ] ) and ( [ eq : product2 ] ) .
this is due to the sign change of @xmath96 , see eqs .
( [ eq : sigmabd1 ] ) . in appendix
[ sec : coefficients ] , we also give the terms for the neutralino decay into a left slepton , @xmath97 .
note that the term proportional to @xmath34 in eq .
is negligible at high particle energies @xmath98 .
] , @xmath99 & &     \left [        m_{\tilde\chi_i^0}^2(p_\tau\cdot p_{\ell_1 } )        -(p_{\tilde\chi_i^0}\cdot p_\tau)(p_{\ell_1}\cdot p_{\tilde\chi_i^0 } )    \right ] , \quad \label{eq : product1}\\ [ 2 mm ] \sigma^{ab}_{\rm d}~\sigma^{b}_{{\rm d}_1 } & = &        \,^{\;\,+}_{(- ) }         \frac{g^4}{2}\left(|a^{{\tilde{\tau}}}_{mi}|^2 + |b^{{\tilde{\tau}}}_{mi}|^2\right )      |f^{r}_{\ell i}|^2 m_\tau    \times\nonumber\\[2 mm ] & & \left [      ( s^a_\tau \cdot p_{\tilde\chi_i^0})(p_{\tilde\chi_i^0}\cdot p_{\ell_1 } )       -m_{\tilde\chi_i^0}^2   ( s^a_\tau \cdot p_{\ell_1 } ) \right ] \nonumber \\ [ 2 mm ] & &   \,^{\;\,+}_{(- ) } g^4 { \mathfrak{re}}\ { a^{{\tilde{\tau}}}_{mi } ( b^{{\tilde{\tau}}}_{mi})^\ast \ }    |f^{r}_{\ell i}|^2   m_{\tilde\chi_i^0 }    \times\nonumber\\ & & \left [      ( p_\tau \cdot p_{\tilde\chi_i^0})(s^a_\tau \cdot p_{\ell_1 } )       - ( p_\tau \cdot p_{\ell_1 } )   ( s^a_\tau \cdot p_{\tilde\chi_i^0 } ) \right ] \nonumber \\ [ 2 mm ] & &   \,^{\;\,+}_{(- ) } g^4     |f^{r}_{\ell i}|^2   m_{\tilde\chi_i^0 }   \times\nonumber\\[2 mm ] & & { \mathfrak{im}}\ { a^{{\tilde{\tau}}}_{mi } ( b^{{\tilde{\tau}}}_{mi})^\ast \ } [ p_{{\tilde{\tau}}},~ p_{\ell_1},~   p_{\tau } , ~s^a_\tau ] .
\label{eq : product2}\end{aligned}\ ] ] the spin - spin correlation term @xmath100 , eq .  , explicitly depends on the imaginary part @xmath101 of the stau - tau - neutralino couplings , eq .  .
thus this term is manifestly cp - sensitive , i.e. , it depends on the phases @xmath102 , @xmath103 , @xmath104 of the stau - tau - neutralino sector .
the imaginary part is multiplied by the totally anti - symmetric ( epsilon ) product , @xmath105 \equiv     \epsilon_{\mu\nu\rho\sigma } \
,         p_{{\tilde{\tau}}}^\mu \ ,        p_{\ell_1}^\nu \ ,        p_{\tau}^\rho \ ,         s^{a , \sigma}_{\tau},\label{eq : epsilon}\ ] ] with the convention @xmath106 .
since each of the spatial components of the four - momenta @xmath107 , or the spin vectors @xmath108 , changes sign under a time transformation , @xmath109 , the epsilon product @xmath110 is t - odd . in the stau rest frame , @xmath111 , the epsilon product reduces to the t - odd triple product @xmath112 @xmath113 = m_{{\tilde{\tau } } } \ ; \ , ( \mathbf{p}_{\ell_1 } \times \mathbf{p}_{\tau } )   \cdot    \mathbf{s}^a_{\tau } \equiv          m_{{\tilde{\tau } } } \ , { \mathcal t}^a .
\label{eq : tripleterm}\ ] ] the task in the next section is to define an observable that projects out from the amplitude squared the part proportional to @xmath110 ( or @xmath112 ) , in order to probe the cp - sensitive coupling combination @xmath114 .
the @xmath14 polarization is given by the expectation value of the pauli matrices @xmath76@xmath115
@xcite @xmath116 with the @xmath14 spin density matrix @xmath117 , as given in eq .  .
in our convention for the polarization vector @xmath118 @xmath119 , the components @xmath120 and @xmath121 are the transverse and longitudinal polarizations in the plane spanned by @xmath122 and @xmath123 , respectively , and @xmath124 is the polarization normal to that plane .
see our definition of the tau spin basis vectors @xmath74 in appendix  [ sec : momenta ] .
+ the normal @xmath14 polarization is equivalently defined as @xmath125 with the number of events @xmath51 with the @xmath14 spin up @xmath126 or down @xmath127 , with respect to the quantization axis @xmath128 , see eq .  .
the normal @xmath14 polarization can thus also be regarded as an asymmetry @xmath129 of the triple product @xmath130 where @xmath131@xmath132 is the direction of the @xmath14 spin vector for each event .
the triple product @xmath133 is included in the spin - spin correlation term @xmath100 , eq .
_ eq .  , and the asymmetry thus probes the term which contains the cp - sensitive coupling combination @xmath114 .
+ since under naive time reversal , @xmath134 , the triple product @xmath133 changes sign , the tau polarization @xmath124 , eq .
, is t - odd . due to cpt invariance  @xcite
, @xmath124 would thus be cp - odd at tree level . in general
, @xmath124 also has contributions from absorptive phases , e.g. from intermediate @xmath135-state resonances or final - state interactions , which do not signal cp violation .
although such absorptive contributions are a higher order effect , and thus expected to be small , they can be eliminated in the true cp asymmetry  @xcite @xmath136 where @xmath137 is the normal tau polarization for the charged conjugated process @xmath93 . for our analysis at tree level ,
where no absorptive phases are present , we find @xmath138 , see the sign change in eqs .   and , and thus @xmath139 .
we study @xmath140 in the following , which is , however , equivalent to @xmath124 at tree level .
+ inserting now the explicit form of the density matrix @xmath117 , eq .  , into eq .  ,
together with eq .
, we obtain the cp asymmetry @xmath141 where we have used the narrow width approximation for the propagators in the phase space element @xmath61 , see eq .  .
note that in the denominator of @xmath140 , eq .  , the spin correlation terms vanish , @xmath142 , see eq .
, when integrated over phase space
. in the numerator only the spin - spin correlation term @xmath143 for @xmath144 contributes , which contains the t - odd epsilon product @xmath110 , see eq .  .
+      to qualitatively understand the dependence of the asymmetry @xmath140 , eq .  , on the mssm parameters , we study in some detail its dependence on the @xmath145-@xmath14-@xmath146 couplings , @xmath147 and @xmath148 , see eq .  .
from the explicit form of the decay terms @xmath149 eq .  , and @xmath79 , @xmath150 , eqs .
, , respectively , we find that the asymmetry @xmath151                           \,{{\rm d}\mathscr{l}\!\textsl{ips } } }                  { ( p_{{\tilde\chi}_i^0}\cdot p_{\tau } )                     ( p_{{\tilde\chi}_i^0}\cdot p_{\ell_1})\,\int{{\rm d}\mathscr{l}\!\textsl{ips } } } , \label{finalmathcala}\end{aligned}\ ] ] with @xmath152 , and @xmath153 , is proportional to the decay coupling factor @xmath154 with @xmath155 $ ] .
we thus expect maximal asymmetries for equal moduli of left and right couplings , @xmath156 , which have a phase difference of about @xmath157 , where the coupling factor can be maximal @xmath158 , see eq .  .
+ to study the dependence of @xmath159 on the cp phase @xmath160 of the stau sector , and the stau mixing angle @xmath161 , we expand the imaginary part of the product of @xmath145-@xmath14-@xmath146 couplings @xmath162   \big\},\end{aligned}\ ] ] in terms of the stau mixing matrix @xmath163 , the gauge couplings @xmath164 and the higgs couplings @xmath165 . in particular , for a cp - conserving neutralino sector , @xmath166 , we have @xmath167 , \nonumber\\ \label{eq : expectedeta}\end{aligned}\ ] ] for @xmath168 , and the sign in parentheses holds for @xmath169 .
thus we expect a maximal @xmath159 and thus maximal asymmetries for maximal stau mixing and @xmath170 decay typically have similar magnitude but opposite sign , and thus might cancel .
see the discussion at the end of the numerics in section  [ sub : mlmr ] .
] , @xmath171 , and a maximal cp phase in the stau mixing matrix , @xmath172 .
note that , in particular , the dependence of @xmath160 on @xmath102 is strong for @xmath173
. we will study numerically the phase and parameter dependence on @xmath140 and @xmath159 further in section  [ sec : results ] .
[ fig : lhcboosts ]    the triple product asymmetry @xmath140 , eq .
( [ p2 ] ) , is not lorentz invariant but depends on the boost of the decaying stau , @xmath174 in fig .  2
, we show the boost dependence of the asymmetry @xmath175 , normalized by @xmath176 .
the susy parameters are given in table  [ tab : benchmark.scenario ] , and we have chosen three sets of different @xmath12 soft - breaking parameters @xmath177  gev ( solid , red ) ; @xmath178  gev ( dashed , green ) ; and @xmath179  gev ( dotted , blue ) .
the corresponding stau masses are @xmath180 ; @xmath181 ; @xmath182  gev , respectively .
the corresponding asymmetries in the stau rest frame are @xmath183 ; @xmath184 , @xmath185 .
note that we have chosen nearly degenerate stau masses which lead to an enhanced stau mixing and thus to maximal asymmetries ; see also the discussion in section  [ sec : results ] .
for the stau masses @xmath186  gev , the staus can be produced at the ilc with @xmath187  gev , and have a fixed boost of @xmath188 .
the corresponding asymmetry is then reduced to @xmath189 if the stau rest frame can not be reconstructed .
typical ilc cross section for these masses are of the order of some @xmath190  fb  @xcite .
if the staus are produced at the lhc , they will have a distinct boost distribution depending on their mass , which typically peaks at high values @xmath191 for stau masses of the order of a few
@xmath192  gev up to a @xmath193  tev , see e.g. refs  @xcite .
then the normal tau polarization in the laboratory frame is obtained by folding the boost dependent polarization @xmath140 with the normalized stau boost distribution  @xcite , @xmath194 with the production cross section @xmath195 .
the typical reduction of the normal tau polarization @xmath196 is of the order of two thirds of the asymmetry compared to that in the stau rest frame @xmath197 .
however , it has been recently shown ( for similar asymmetries in stop decays at the lhc ) , that the rest frame can be partly reconstructed event by event using on - shell mass conditions , see refs .
the lhc cross section for stau pair production , @xmath198 , also sensitively depends on the stau masses , e.g. , for our benchmark scenario in table  [ tab : benchmark.scenario ] , we find cross sections up to @xmath199  fb at @xmath200  tev  @xcite .    .
benchmark scenario .
the mass parameters @xmath201 , @xmath202 , @xmath0 , @xmath203 , @xmath204 @xmath205 , and @xmath206 are given in gev . [ cols="^,<,^,<,^,<,^,<,^,<,^,<,^",options="header " , ]     [ tab : benchmark.masses ]
we quantitatively study the tau polarization asymmetry , and the branching ratios for the two - body decay chain @xmath207 for @xmath208 .
the asymmetry probes the mssm phases @xmath209 , @xmath104 and @xmath102 , of the neutralino and stau sector .
we center our numerical discussion around a general mssm benchmark scenario , see table  [ tab : benchmark.scenario ] .
we choose heavier soft breaking parameters in the stau sector than in the @xmath210 sector , to enable the mass hierarchy @xmath211 further we choose almost degenerate staus which enhances their mixing , leading to maximal asymmetries .
we choose a large value of the trilinear scalar coupling parameter , @xmath173 is restricted bz the vacuum stability contidtion as @xmath212  @xcite .
] , to enhance the impact of @xmath102 in the stau sector .
finally , to reduce the number of mssm parameters , we use the ( gut inspired ) relation @xmath213  @xcite for the gaugino mass parameters .
the resulting masses of the staus , neutralinos and charginos are summarized in table  [ tab : benchmark.masses ] .      for the benchmark scenario given in table  [ tab : benchmark.scenario ] ,
we study the phase dependence of the asymmetry @xmath140 in the stau rest frame . in fig .
[ fig : phimu_m1_asym ] , we show the dependence on the cp phases in the neutralino sector , @xmath103 and @xmath104 . in fig .
[ fig : phimu_atau_asym ] , we show the dependence on the phases in the stau secotor @xmath102 and @xmath104 .
the asymmetry strongly depends on @xmath214 , which we expect for @xmath215 as in our benchmark scenario , see table  [ tab : benchmark.scenario ] . in particular for @xmath216 in fig .
[ fig : phimu_atau_asym ] , the asymmetry follows the approximation formula eq .  , and attains its maximal values at @xmath217 .      in fig .
[ fig : a_tanbeta_asym ] , we show the @xmath218 and @xmath219 dependence of the asymmetry @xmath140 in the stau rest frame .
we can observe that the asymmetry obtains its maximum , @xmath220 , where also the coupling factor is maximal , @xmath221 , see fig .
[ fig : a_tanbeta_eta_pi4 ] .
as discussed in subsection  [ sub : analytics ] , the imaginary part of the product of the stau couplings @xmath222 is maximal for a maximal cp phase @xmath223 in the stau sector , which we show in fig .  [
fig : a_tanbeta_phistau_pi4 ] . note that the location of the maximum of @xmath140 is not at maximal stau mixing , @xmath224 , since @xmath225 starts to decrease for increasing @xmath0 and @xmath219 .    to study the stau mixing , we show the @xmath226@xmath227 dependence of the asymmetry @xmath140 in fig .  [
fig : ml_mr_asym ] . in the entire @xmath226@xmath227 plane ,
the cp phase in the stau sector is almost maximal , @xmath228 .
however , the asymmetry obtains its maxima in the small corridor @xmath229 , where the stau mixing is maximal , @xmath230 .
we show the @xmath202@xmath201 dependence of the asymmetry @xmath140 in fig .
[ fig : mu_m2_asym ] .
the maxima of @xmath140 are obtained where the coupling factor @xmath159 is also maximal , see eq .
( [ etafactor ] ) .    in fig .
[ fig : mu_m2_br_stau ] , we show the corresponding stau branching ratio , @xmath231 , which can be as large as @xmath232%
. other competing channels can reach @xmath233 , and @xmath234 .
the stau decay into the chargino @xmath235 is always open since typically the second lightest neutralino and the lightest chargino are almost degenerate , @xmath236 .
the neutralino branching ratio @xmath237 , summed over @xmath238 , is shown in fig .  [
fig : mu_m2_br_chi ] , which reaches up to @xmath192% .
the other important competing decay channels are @xmath239 , and @xmath240 , which open around @xmath241  gev and @xmath242  gev , respectively , for @xmath243  gev .
note that in our benchmark scenario , see table  [ tab : benchmark.scenario ] , we have @xmath244 .      as we discussed in section  [ sub : atautanb ] , we find large asymmetries for nearly degenerate staus , where we naturally obtain a maximal stau mixing .
however , then typically the asymmetries for @xmath246 and @xmath170 decay are similar in magnitude , but opposite in sign .
for example in our benchmark scenario we find @xmath247 for @xmath246 decay , but @xmath248 for the decay of @xmath170 .
if the production and decay process of @xmath246 can not be experimentally disentangled from that of @xmath170 properly , the two asymmetries might cancel .
we show their sum in fig .
[ fig : ml - mr_asymm_k1-plus - k2 ] in the @xmath226@xmath227 plane . in fig .
[ fig : ml - mr_mstau2-mstau1 ] , we show the corresponding stau mass splitting .
note that also the stau branching ratios are similar in size ; for example in our benchmark scenario we have @xmath249 , and @xmath250 . for the @xmath226@xmath227 plane shown in fig .
[ fig : ml_mr _ ] , the decay branching ratio @xmath251 is at least @xmath252 , and that of @xmath253 is larger by roughly a factor of @xmath254 to @xmath255 .
we have analyzed the normal tau polarization and the corresponding cp asymmetry in the two - body decay chain of a stau @xmath256 the cp - sensitive parts appear only in the spin - spin correlations , which can be probed by the subsequent neutralino decay @xmath257 for @xmath208 .
the t - odd tau polarization normal to the plane spanned by the @xmath14 and @xmath22 momenta , can then be used to define a cp - odd tau polarization asymmetry .
it is based on a triple product , which probes the cp phases of the trilinear scalar coupling parameter @xmath0 , the higgsino mass parameter @xmath2 , and the u(1 ) gaugino mass parameter @xmath1 .
+ we have analyzed the analytical and numerical dependence of the asymmetry on these parameters in detail .
in particular , for nearly degenerate staus where the stau mixing is strong , the asymmetry obtains its maxima and can be larger than @xmath3 .
the normal tau polarization can thus be considered as an ideal cp observable to probe the cp phases in the stau and neutralino sector of the mssm .
+ since the cp - sensitive parts appear only in the subsequent stau decay products the stau production process can be separated .
thus both , ilc , and lhc collider studies are possible . concerning the kinematical dependence , the asymmetry is not lorentz invariant , since it is based on a triple product . at the lhc ,
staus are produced with a distinct boost distribution .
evaluated in the laboratory frame , the resulting tau polarization asymmetries get typically reduced by a factor of two thirds , compared to the stau rest frame . + we want to stress that a thorough experimental analysis , addressing background processes , detector properties , and event rate reconstruction efficiencies , will be needed in order to explore the measurability of cp phases in the stau sector at the lhc or ilc .
we hope that our work motivates such a study .
we thank m.  drees and f.  von  der  pahlen for enlightening discussions and helpful comments .
this work has been supported by micinn project fpa.2006 - 05294 .
am was supported by the _
konrad adenauer stiftung _ ,
bcgs , and a fellowship of bonn university .
hd was supported by the hemholtz alliance `` physics at the terascale '' and bmbf `` verbundprojekt hep - theorie '' under the contract 0509pde .
sk was supported by bcgs .
ok acknowledges support from cpan .
for the stau decay @xmath258 , we choose the coordinate frame in the laboratory ( lab ) system , such that the momentum of decaying @xmath259 points in the @xmath260-direction .
@xmath261    p^{\mu}_{\tau } & = &   e_\tau                       ( 1 , \sin\theta_{\tau } , 0 , \cos\theta_{\tau } ) ,   \ ] ] with the decay angle @xmath262 , and @xmath263 in the limit @xmath264 .
the momenta of the leptons from the subsequent neutralino decay @xmath265 ; @xmath266  ( [ eq : staudec ] ) , can be parameterized by @xmath267    p^{\mu}_{\ell_2 } & = & e_{\ell_2 }                ( 1 , \sin\theta_2\cos\phi_2 , \sin\theta_2\sin\phi_2,\cos\theta_2),\end{aligned}\ ] ] with the energies @xmath268 and the decay angles @xmath269 , @xmath270 , that is , @xmath271 \cos\theta_{d_2 } & = &      \frac{({\mathbf p}_{{\tilde{\tau}}}-{\mathbf p}_{\tau}-{\mathbf p}_{\ell_1 }     ) \cdot\hat{\mathbf p}_{\ell_2 } }          { | { \mathbf p}_{{\tilde{\tau}}}-{\mathbf p}_{\tau}-{\mathbf p}_{\ell_1 } | } , \end{aligned}\ ] ] with the unit momentum vector @xmath272 .
we define the tau spin vectors by @xmath273     s_{\tau}^{3,\mu}&=&\frac{1}{m_{\tau } }                          \left(|{\bf{p}}_{\tau}| ,                          \frac{e_{\tau}}{|{\bf{p}}_{\tau}| }                           { \bf{p}}_{\tau}\right ) .
\label{eq : spintau}\end{aligned}\ ] ] the spin vectors @xmath274 @xmath275 for the tau , and @xmath276 @xmath277 for the neutralino @xmath278 , fulfil completeness relations @xmath279 \sum_{b}s_{\tilde\chi_i^0}^{b,\,\mu } s_{\tilde\chi_i^0}^{b,\,\nu }      & = & -g^{\mu\nu }        + \frac{p_{\tilde\chi_i^0}^\mu p_{\tilde\chi_i^0}^\nu}{m_{\tilde\chi_i^0}^2 } , \label{eq : completeness2}\end{aligned}\ ] ] and they form orthonormal sets @xmath280            s_{{\tilde{\chi}^0_i}}^b \cdot s_{{\tilde{\chi}^0_i}}^c & = & -\delta^{bc } , \qquad             s_{{\tilde{\chi}^0_i}}^b \cdot \hat{p}_{{\tilde{\chi}^0_i } } = 0,\end{aligned}\ ] ] with @xmath281 .
note that the asymmetry @xmath140 , eq .
( [ p2 ] ) , does not depend on the explicit form of the neutralino spin vectors , since they are summed in the amplitude squared , see eq .
( [ eq : product1 ] ) , using the completeness relation .
the lorentz invariant phase - space element for the stau decay chain , see eqs
.   - , can be decomposed into two - body phase - space elements  @xcite @xmath282 \times   { \rm d}s_{{\tilde{\chi}^0_i}}\ , { { \rm d}\mathscr{l}\!\textsl{ips}}(s_{{\tilde{\chi}^0_i}};\,p_{\ell_1},p_{\tilde\ell})\ , { \rm d}s_{\tilde\ell}\ , { { \rm d}\mathscr{l}\!\textsl{ips}}(s_{\tilde\ell};\,p_{\ell_2},p_{{\tilde\chi}_1 ^ 0}).\quad\end{aligned}\ ] ] the different contributions are @xmath283 with @xmath284 and @xmath285 .
the coefficients of the stau decay matrix , eq .  , are @xmath286     & & -g^2 { \mathfrak{re}}\{a_{mi}^{{\tilde{\tau } } } \,(b_{mi}^{{\tilde{\tau}}})^\ast\ } m_{{\tilde{\chi}^0_i}}m_\tau , \label{eq : d}\\ \sigma_{\rm d}^a & = &           \,^{\;\,-}_{(+ ) }          \frac{g^2}{2}\left(|a^{{\tilde{\tau}}}_{mi}|^2-|b^{{\tilde{\tau}}}_{mi}|^2\right )          m_\tau(p_{\tilde\chi_i^0}\cdot s^a_{\tau } ) , \label{eq : sigmaad}\\ \sigma_{\rm d}^b & = &     \,^{\;\,-}_{(+ ) }       \frac{g^2}{2 } \left(|a^{{\tilde{\tau}}}_{mi}|^2-|b^{{\tilde{\tau}}}_{mi}|^2\right )        m_{{\tilde{\chi}^0_i } } ( p_\tau\cdot s_{{\tilde{\chi}^0_i}}^b ) , \label{eq : sigmabd }
\\ \sigma_{\rm d}^{ab } & = &    \frac{g^2}{2}\left(|a^{{\tilde{\tau}}}_{mi}|^2 + |b^{{\tilde{\tau}}}_{mi}|^2\right )           ( s^a_{\tau}\cdot s^b_{{\tilde\chi_i^0}})m_\tau m_{\tilde\chi_i^0 }          \nonumber \\ [ 2 mm ] & &           + g^2 { \mathfrak{re}}\ { a^{{\tilde{\tau}}}_{mi } ( b^{{\tilde{\tau}}}_{mi})^\ast \}\times   \nonumber \\ [ 2 mm ] & & \left [   ( s^a_{\tau}\cdot    p_{\tilde\chi_i^0 } ) ( s^b_{\tilde\chi_i^0}\cdot   p_\tau )          -(s^a_{\tau}\cdot s^b_{\tilde\chi_i^0 } ) ( p_{\tilde\chi_i^0 }   \cdot   p_\tau )   \right ]          \nonumber \\ [ 2 mm ] & &          -g^2 { \mathfrak{im}}\ { a^{{\tilde{\tau}}}_{mi } ( b^{{\tilde{\tau}}}_{mi})^\ast \ }           [ s^{a}_{\tau},~ p_\tau,~ s^{b}_{\tilde\chi_i^0},~p_{\tilde\chi_i^0 } ] .
\label{eq : sigmaabd }   \ ] ] the formulas are given for the decay of a negatively charged stau , @xmath287 .
the signs in parentheses hold for the charge conjugated decay @xmath93 .
+ note that the terms proportional to @xmath34 in eqs .
, , and , are negligible at high particle energies @xmath98 , in particular @xmath88 can be neglected .
+ the coefficients of the @xmath288 decay matrix , eq .  , are  @xcite @xmath289 \sigma_{{\rm d}_1}^b&= & \,^{\;\,+}_{(- ) } g^2|f_{\ell i}^r|^2 m_{{\tilde{\chi}^0_i}}(s_{{\tilde{\chi}^0_i}}^b\cdot p_{\ell_1 } ) , \label{eq : sigmabd1}\end{aligned}\ ] ] and the selectron decay factor is @xmath290 the signs in parentheses hold for the charge conjugated processes , that is @xmath95 in eq .
( [ eq : sigmabd1 ] ) . + for the decay into a left slepton @xmath291 , eqs .
( [ eq : d1 ] ) , ( [ eq : sigmabd1 ] ) , and ( [ eq : d2 ] ) read  @xcite @xmath292 \sigma^b_{{\rm d}_1 } & = & \,^{\;\,-}_{(+ ) } g^2 |f^{l}_{\ell i}|^2   m_{\tilde\chi_i^0 } ( s^b_{\tilde\chi_i^0 } \cdot p_{\ell_1 } ) , \label{eq : sigmabd1sell } \\ [ 2 mm ]          { \rm d}_2&= & g^2 |f^{l}_{\ell 1}|^2           ( m_{\tilde\ell}^2-m_{\tilde\chi_1 ^ 0}^2 ) , \label{eq : d2sell}\end{aligned}\ ] ] respectively .
the expressions for eqs .
( [ eq : product1 ] ) and ( [ eq : product2 ] ) have to be changed accordingly . the sign in parenthesis in eq .
( [ eq : sigmabd1sell ] ) holds for the charge conjugated process @xmath293 .
the partial decay width for the decay @xmath294 in the stau rest frame is  @xcite @xmath295 with the decay function @xmath296 given in eqs .
( [ eq : d ] ) , and the approximation @xmath297 . for the decay
@xmath298 the width is  @xcite @xmath299 with the stau - chargino - neutrino coupling  @xcite @xmath300 and the stau diagonalization matrix @xmath301 , eq .
( [ eq : rstau ] ) , the yukawa coupling @xmath302 , eq .
( [ eq : yt ] ) , and the matrix @xmath303 , that diagonalizes the chargino matrix  @xcite , @xmath304 the stau decay width for the entire decay chain , eqs
.   - , is then given by @xmath305 & = & \frac{1}{2m_{\tilde\tau } } \int\,|\mathcal{m}|^2\ , { { \rm d}\mathscr{l}\!\textsl{ips}}(s_{{\tilde{\tau}}};\,p_\tau , p_{\ell_1},p_{\ell_2},p_{\tilde\chi_1 ^ 0 } ) \\[2 mm ] & = & \gamma(\tilde\tau)\times { \rm br}(\tilde\tau\to\tau \tilde\chi_i^0 )                    \times { \rm br}(\tilde\chi_i^0\to \ell_1\tilde\ell ) \nonumber\\ & &                    \times { \rm br}(\tilde\ell\to \ell_2\tilde\chi_1 ^ 0 ) ,    \\\end{aligned}\ ] ] with the phase - space element @xmath306 , as given in the appendix  [ sec : momenta ] , the amplitude squared latexmath:[\[\begin{aligned }    obtained from eqs .   by summing the tau helicities
@xmath62 , @xmath308 . the neutralino branching ratios are given , for example , in ref .
@xcite , and we assume @xmath309 .
we use the narrow width approximation for the propagators @xmath310 which is justified for @xmath311 , which holds in our case with @xmath312 .
note , however , that in principle the naive @xmath313-expectation of the error can easily receive large off - shell corrections of an order of magnitude , and more , in particular at threshold , or due to interferences with other resonant , or non - resonant processes  @xcite .
+    0 belle collab .
, a.  abashian _ et al .
_ , phys .
lett . * 86 * , 2509 ( 2001 ) , hep - ex/0102018 ; + babar collab . , b.  aubert _ et al .
_ , phys .
lett . * 86 * , 2515 ( 2001 ) , hep - ex/0102030 .
a.  d.  sakharov , zh .
pisma * 5 * , 32 ( 1967 ) , jetp  lett . * 91b * , 24 ( 1967 ) . h.  e.  haber and g.  l.  kane , phys .
* 117 * ( 1985 ) 75 ; + h.  p.  nilles , phys .
* 110 * ( 1984 ) 1 ; + m.  drees , r.  godbole and p.  roy , `` theory and phenomenology of sparticles : an account of four - dimensional n=1 supersymmetry in high energy physics , '' _ hackensack , usa : world scientific ( 2004)_.      t.  ibrahim and p.  nath , phys .
d * 57 * ( 1998 ) 478 [ erratum - ibid .
d * 58 * ( 1998  errat , d60,079903.1999   errat , d60,119901.1999 ) 019901 ] [ arxiv : hep - ph/9708456 ] ; phys .
b * 418 * ( 1998 ) 98 [ arxiv : hep - ph/9707409 ] ; [ erratum - ibid .
d * 60 * ( 1999 ) 099902 ] [ arxiv : hep - ph/9807501 ] ; phys .
d * 61 * ( 2000 ) 093004 [ arxiv : hep - ph/9910553 ] ; + m.  brhlik , g.  j.  good and g.  l.  kane , phys .
d * 59 * ( 1999 ) 115004 [ arxiv : hep - ph/9810457 ] ; + s.  yaser ayazi and y.  farzan , phys .
d * 74 * , 055008 ( 2006 ) [ arxiv : hep - ph/0605272 ] .
see , e.g. , a.  bartl , t.  gajdosik , w.  porod , p.  stockinger and h.  stremnitzer , phys .
d * 60 * ( 1999 ) 073003 [ arxiv : hep - ph/9903402 ] ; + a.  bartl , t.  gajdosik , e.  lunghi , a.  masiero , w.  porod , h.  stremnitzer and o.  vives , phys .
d * 64 * ( 2001 ) 076009 [ arxiv : hep - ph/0103324 ] ; + l.  mercolli and c.  smith , nucl .
b * 817 * , 1 ( 2009 ) [ arxiv:0902.1949 [ hep - ph ] ] ; + v.  d.  barger , t.  falk , t.  han , j.  jiang , t.  li and t.  plehn , phys .
d * 64 * ( 2001 ) 056007 [ arxiv : hep - ph/0101106 ] .
a.  bartl , w.  majerotto , w.  porod and d.  wyler , phys .
d * 68 * ( 2003 ) 053005 [ arxiv : hep - ph/0306050 ] ; + k.  a.  olive , m.  pospelov , a.  ritz and y.  santoso , phys .
d * 72 * ( 2005 ) 075001 [ arxiv : hep - ph/0506106 ] ; + j.  r.  ellis , j.  s.  lee and a.  pilaftsis , jhep * 0810 * ( 2008 ) 049 [ arxiv:0808.1819 [ hep - ph ] ] ; + s.  a.  abel , a.  dedes and h.  k.  dreiner , jhep * 0005 * ( 2000 ) 013 [ arxiv : hep - ph/9912429 ] .
s.  y.  choi , a.  djouadi , m.  guchait , j.  kalinowski , h.  s.  song and p.  m.  zerwas , eur .
j.   c * 14 * ( 2000 ) 535 [ arxiv : hep - ph/0002033 ] ; + j.  l.  kneur and g.  moultaka , phys .
d * 61 * ( 2000 ) 095003 , arxiv : hep - ph/9907360 .
a.  pilaftsis _
et al . _ ,
b * 435 * ( 1998 ) 88 , arxiv : hep - ph/9805373 ; + d.  a.  demir _ et al .
_ , phys .
d * 60 * ( 1999 ) 055006 , arxiv : hep - ph/9807336 ; + a.  pilaftsis and c.  e.  m.  wagner , nucl .
b * 553 * ( 1999 ) 3 , arxiv : hep - ph/9902371 ; + s.  y.  choi , m.  drees , and j.  s.  lee , phys .
b * 481 * ( 2000 ) 57 , arxiv : hep - ph/0002287 ; + j.  l.  kneur and g.  moultaka , phys .
d * 59 * ( 1999 ) 015005 , arxiv : hep - ph/9807336 .
a.  bartl , s.  hesselbach , k.  hidaka , t.  kernreiter and w.  porod , phys .
d * 70 * ( 2004 ) 035003 [ arxiv : hep - ph/0311338 ] ; + k.  rolbiecki , j.  tattersall and g.  moortgat - pick , arxiv:0909.3196 [ hep - ph ] .
h.  eberl , t.  gajdosik , w.  majerotto and b.  schrausser , phys .
b * 618 * ( 2005 ) 171 [ arxiv : hep - ph/0502112 ] ; + e.  chistova , h.  eberl , w.  majerotto , and s.  kraml , nucl .
b * 639 * ( 2002 ) 263 , arxiv : hep - ph/0205227 ; + e.  chistova , h.  eberl , w.  majerotto , and s.  kraml , jhep * 12 * ( 2002 ) 021 , arxiv : hep - ph/0211063 ; + m.  frank and i.  turan , phys .  rev .
d * 76 * ( 2007 ) 016001 , arxiv : hep - ph/0703184 ; + m.  frank and i.  turan , phys .
d * 76 * ( 2007 ) 076008 [ arxiv:0708.0026 [ hep - ph ] ] .
e.  christova , h.  eberl , e.  ginina and w.  majerotto , phys .
d * 79 * , 096005 ( 2009 ) [ arxiv:0812.4392 [ hep - ph ] ] .
e.  christova , h.  eberl , e.  ginina and w.  majerotto , jhep * 0702 * ( 2007 ) 075 [ arxiv : hep - ph/0612088 ] .
a.  pilaftsis , nucl .
b * 504 * , 61 ( 1997 ) [ arxiv : hep - ph/9702393 ] ; + s.  y.  choi , j.  kalinowski , y.  liao and p.  m.  zerwas , eur .
j.   c * 40 * , 555 ( 2005 ) [ arxiv : hep - ph/0407347 ] ; + h.  k.  dreiner , o.  kittel and f.  von der pahlen , jhep * 0801 * , 017 ( 2008 ) [ arxiv:0711.2253 [ hep - ph ] ] ; + o.  kittel and f.  von der pahlen , jhep * 0808 * , 030 ( 2008 ) [ arxiv:0806.4534 [ hep - ph ] ] ; + m.  nagashima , k.  kiers , a.  szynkman , d.  london , j.  hanchey and k.  little , phys .
d * 80 * , 095012 ( 2009 ) [ arxiv:0907.1063 [ hep - ph ] ] . for studies with neutralino
3-body decays at the ilc , see + y.  kizukuri and n.  oshimo , phys .
b * 249 * ( 1990 ) 449 ; + s.  y.  choi , h.  s.  song and w.  y.  song , phys .
d * 61 * ( 2000 ) 075004 [ arxiv : hep - ph/9907474 ] ; + for further studies with neutralino 2-body and 3-body decays at the ilc , see + a.  bartl , h.  fraas , o.  kittel and w.  majerotto , phys .
d * 69 * ( 2004 ) 035007 [ arxiv : hep - ph/0308141 ] ; + a.  bartl , h.  fraas , o.  kittel and w.  majerotto , eur .
j.   c * 36 * ( 2004 ) 233 [ arxiv : hep - ph/0402016 ] ; + j.  a.  aguilar - saavedra , nucl .
b * 697 * ( 2004 ) 207 [ arxiv : hep - ph/0404104 ] ; + a.  bartl , h.  fraas , s.  hesselbach , k.  hohenwarter - sodek and g.  a.  moortgat - pick , jhep * 0408 * ( 2004 ) 038 [ arxiv : hep - ph/0406190 ] ; + s.  y.  choi , b.  c.  chung , j.  kalinowski , y.  g.  kim and k.  rolbiecki , eur .
j.   c * 46 * ( 2006 ) 511 [ arxiv : hep - ph/0504122 ] ; + for studies with chargino 2-body decays at the ilc , see + a.  bartl , h.  fraas , o.  kittel and w.  majerotto , phys .
b * 598 * ( 2004 ) 76 [ arxiv : hep - ph/0406309 ] ; + o.  kittel , a.  bartl , h.  fraas and w.  majerotto , phys .
d * 70 * ( 2004 ) 115005 [ arxiv : hep - ph/0410054 ] ; + for studies with chargino 3-body decays at the ilc , see + y.  kizukuri and n.  oshimo , arxiv : hep - ph/9310224 ; + a.  bartl , h.  fraas , s.  hesselbach , k.  hohenwarter - sodek , t.  kernreiter and g.  moortgat - pick , eur .
j.   c * 51 * ( 2007 ) 149 [ arxiv : hep - ph/0608065 ] . for recent reviews
see , for example , + g.  moortgat - pick , k.  rolbiecki , j.  tattersall and p.  wienemann , arxiv:0910.1371 [ hep - ph ] ; + o.  kittel , arxiv:0904.3241 [ hep - ph ] .
s.  abdullin _ et al . _
[ cms collaboration ] , j.  phys .
g * 28 * ( 2002 ) 469 [ arxiv : hep - ph/9806366 ] ; + atlas collab . , _ atlas detector and physics performance . technical design report .
vol . 2 _ , cern - lhcc-99 - 15 ; + g.  weiglein _ et al . _
[ lhc / lc study group ] , arxiv : hep - ph/0410364 . j.  brau _ et al .
_ [ ilc collaboration ] , arxiv:0712.1950 [ physics.acc-ph ] ; + j.  a.  aguilar - saavedra _ et al . _ [ ecfa / desy lc physics working group ] , arxiv : hep - ph/0106315 ; + t.  abe _ et al . _ [ american linear collider working group ] , arxiv : hep - ex/0106055 ; + k.  abe _ et al . _ [ acfa linear collider working group ] , arxiv : hep - ph/0109166 ; + j.  a.  aguilar - saavedra _ et al .
_ , eur .
j.  c * 46 * ( 2006 ) 43 [ arxiv : hep - ph/0511344 ] .
y.  k.  semertzidis , nucl .
* 131 * ( 2004 ) 244 [ arxiv : hep - ex/0401016 ] ; + j.  r.  ellis , s.  ferrara and d.  v.  nanopoulos , phys .
b * 114 * ( 1982 ) 231 ; + w.  buchmuller and d.  wyler , phys .
b * 121 * ( 1983 ) 321 ; + f.  del aguila , m.  b.  gavela , j.  a.  grifols and a.  mendez , phys .
b * 126 * ( 1983 ) 71 [ erratum - ibid .
b * 129 * ( 1983 ) 473 ] ; + d.  v.  nanopoulos and m.  srednicki , phys .
b * 128 * ( 1983 ) 61 ; + m.  dugan , b.  grinstein and l.  j.  hall , nucl .
b * 255 * ( 1985 ) 413 ; + c.  s.  huang and w.  liao , phys .
d * 62 * ( 2000 ) 016008 [ arxiv : hep - ph/0001174 ] .
s.  y.  choi , m.  drees and b.  gaissmaier , phys .
d * 70 * ( 2004 ) 014010 [ arxiv : hep - ph/0403054 ] .
f.  deppisch and o.  kittel , jhep * 0909 * , 110 ( 2009 ) [ erratum - ibid .
* 1003 * , 091 ( 2010 ) ] [ arxiv:0905.3088 [ hep - ph ] ] .
g.  moortgat - pick , k.  rolbiecki , j.  tattersall and p.  wienemann , jhep * 1001 * , 004 ( 2010 ) [ arxiv:0908.2631 [ hep - ph ] ] ; + g.  moortgat - pick , k.  rolbiecki and j.  tattersall , arxiv:1008.2206 [ hep - ph ] .
a.  bartl , e.  christova , k.  hohenwarter - sodek and t.  kernreiter , phys .
d * 70 * ( 2004 ) 095007 [ arxiv : hep - ph/0409060 ] .
j.  ellis , f.  moortgat , g.  moortgat - pick , j.  m.  smillie and j.  tattersall , eur .
j.   c * 60 * , 633 ( 2009 ) [ arxiv:0809.1607 [ hep - ph ] ] .
a.  bartl , e.  christova , k.  hohenwarter - sodek and t.  kernreiter , jhep * 0611 * ( 2006 ) 076 [ arxiv : hep - ph/0610234 ] .
f.  f.  deppisch and o.  kittel , jhep * 1006 * , 067 ( 2010 ) [ arxiv:1003.5186 [ hep - ph ] ] .
a.  bartl , h.  fraas , t.  kernreiter and o.  kittel , eur .
j.   c * 33 * ( 2004 ) 433 [ arxiv : hep - ph/0306304 ] .
a.  bartl , t.  kernreiter and o.  kittel , phys .
b * 578 * ( 2004 ) 341 [ arxiv : hep - ph/0309340 ] ; + o.  kittel , arxiv : hep - ph/0311169 .
h.  dreiner , o.  kittel and a.  marold , arxiv:1001.4714 [ hep - ph ] .
h.  e.  haber , proceedings of the 21st slac summer institute on particle physics , eds .
l.  deprocel , ch .
dunwoodie , stanford 1993 , 231 .
[ arxiv : hep - ph/9405376 ] .
a.  bartl , k.  hidaka , t.  kernreiter and w.  porod , phys .
b * 538 * ( 2002 ) 137 phys .
d * 66 * ( 2002 ) 115009 [ arxiv : hep - ph/0207186 ] .
o.  kittel , [ arxiv : hep - ph/0504183 ] .
f.  m.  renard , `` basics of electron positron collisions , '' _ dreux , france : editions frontieres ( 1981)_. g.  luders , kong .
sel .  mat .
fys .  med .
* 28n5 * ( 1954 ) 1 ; + w.  pauli , `` niels bohr and the development of physics , '' _ new york , usa : mc graw - hill ( 1955 ) _ ; + r.  jost , helv .
acta * 30 * ( 1957 ) 409 ; + r.  jost , helv .
acta * 36 * ( 1963 ) 77 ; + r.  f.  streater and a.  s.  wightman , `` pct , spin and statistics , and all that , '' _ redwood city , usa : addison - wesley ( 1989 ) ( advanced book classics ) .
_ j.  alwall _ et al .
_ , jhep * 0709 * , 028 ( 2007 ) [ arxiv:0706.2334 [ hep - ph ] ] . j.  m.  frere , d.  r.  t.  jones , and s.  rabi , nucl
.  phys .
* b222 * ( 1983 ) 11 ; m.  claudson , l.  j.  hall , and i.  hinchliffe , nucl .
* b228 * ( 1983 ) 501 ; c.  kounnas , a.  b.  lahanas , d.  v.  nanopoulos , and m.  quiros , nucl .
* b236 * ( 1984 ) 438 ; j.  f.  gunion , h.  e.  haber , and m.  sher , nucl .
phys .  * b306 * ( 1988 ) 1 .
e.  byckling , k.  kajantie , `` particle kinematics , '' _ london , england : john  wiley & sons ( 1973 ) _ ; + g.  costa _ et al .
_ , tepp `` kinematics and symmetries '' , bd.1 , ed . m.  nicoli , _ paris , france : institut national de physique nuclaire et de physique des particules ( 1979)_. k.  hagiwara _ et al .
_ , phys .
d * 73 * ( 2006 ) 055005 [ arxiv : hep - ph/0512260 ] ; + d.  berdine , n.  kauer and d.  rainwater , phys .
* 99 * ( 2007 ) 111601 [ arxiv : hep - ph/0703058 ; + n.  kauer , phys .
b * 649 * ( 2007 ) 413 [ arxiv : hep - ph/0703077 ] ; jhep * 0804 * ( 2008 ) 055 [ arxiv:0708.1161 [ hep - ph ] ] ; + c.  f.  uhlemann and n.  kauer , nucl .
b * 814 * ( 2009 ) 195 [ arxiv:0807.4112 [ hep - ph ] ] ; + m.  a.  gigg and p.  richardson , arxiv:0805.3037 [ hep - ph ] .