as the large hadron collider ( lhc ) continues to successfully probe the nature of electroweak symmetry breaking ( a recent breakthrough being the discovery of a higgs candidate  @xcite ) , a solution to the gauge hierarchy problem , if it exists , has evaded our efforts so far .
an elegant candidate for such a solution and one which addresses several other open issues as well , _
e.g. _ dark matter and gauge coupling unification , is supersymmetry ( susy ) .
because of its theoretical appeal , it is important to understand all the guises that susy may adopt in order to recognize it if it is produced at the lhc .    an important open issue in susy , which has strong ramifications for its lhc phenomenology , as well as its cosmology , is the gauge invariance of lepton and baryon number violating interactions . aside from introducing many new unknown parameters , these interactions also lead to rapid proton decay .
the most common solution to address this issue is the imposing of a discrete symmetry , r - parity , defined as @xmath9 ( see @xcite for a review ) .
this forbids all tree - level lepton and baryon number violating terms and also causes the lightest supersymmetric particle ( lsp ) to be stable and therefore neutral .
this lsp can then play the role of dark matter and would manifest itself at colliders as missing energy .
however , proton decay requires the elimination of only the lepton number or baryon number violating terms and from a theoretical perspective , r - parity can be _
ad hoc_. furthermore , an open - mindedness to possible signals at the lhc should push us to consider alternatives .
finally , when r - parity is violated the stringent constraints on superpartner masses , derived from negative collider searches for missing energy events at the tevatron and the lhc @xcite , can be relaxed .    a systematic study of all possible r - parity violating terms and their effects on phenomenology is an arduous task and furthermore
, one would like a mechanism for understanding why proton decay is significantly suppressed . a well - motivated solution to both of these issues
are models which can predict the fate of r - parity .
a natural framework for this endeavor is in the context of @xmath10 symmetries , see  @xcite for early examples can be also used to construct models where the rpv couplings , arising from effective operators , are intrinsically small @xcite . ] .
while some @xmath11 models predict r - parity conservation  @xcite , the most minimal ones ( in terms of particle content ) require r - parity violation ( rpv )  @xcite .
even some non - minimal models prefer rpv from considerations of the renormalization group evolution of the soft masses  @xcite .
these and other models of spontaneous r - parity violation , such as  @xcite , have the common feature that they can be described , in an effective field theory way , by bilinear r - parity violation ( brpv ) : the only r - parity violating terms are the mixings between the lepton doublets and up - type higgs doublet in the superpotential .
this makes brpv a powerful tool for studying possible signatures of spontaneous rpv .
furthermore , proton decay is highly suppressed and it is important to mention that a gravitino lsp can be a dark matter candidate in such models  @xcite .    once r - parity is broken the lsp is no longer stable and therefore astrophysical constraints on its nature do not apply @xcite .
accordingly , from a purely phenomenological point of view any superpartner can be the lsp , and studies of the different possibilities in brpv models ( and of sleptons and sneutrinos in general models  @xcite ) and their relation to neutrino masses and mixings have been carried out in e.g. @xcite .
in this paper we extend upon previous results by considering sneutrino lsp decays in brpv into `` non - conventional '' modes : pairs of heavy standard model ( sm ) states ( @xmath12 , @xmath13 , @xmath14 , @xmath15 ) , invisible modes ( @xmath16 ) and different - flavor charged lepton final states ( @xmath17 , @xmath18 ) .
such states , to our knowledge , had not received much attention before despite the fact that they can dominate sneutrino decays and can yield unique and unanticipated susy signals . specifically , neutrino masses force rpv to be small so its only effect is on the decay of the lsp . therefore ,
for a sneutrino lsp , every susy event will eventually decay into two sneutrinos which could then decay into one of the aforementioned states .
as we will show , this typically translates into into different - flavor charged leptons or heavy sm modes , depending on the sneutrino generation .
the main goal of this paper is to study the sneutrino decays into these `` non - conventional '' modes and to show that they can dominate over the traditionally considered @xmath19 final state .
we find that the latter dominate roughly when the brpv term , @xmath20 is smaller than the vacuum expectation value ( vev ) of the sneutrino , @xmath21 , for a given flavor of sneutrino , @xmath22 . in the case
when brpv is the only source of neutrino masses , this possibility can only hold for one sneutrino flavor , however , when this assumption is relaxed , it can hold true for all three generations .
therefore , if two or more generations of sneutrinos decay via rpv , it might be possible to rule out brpv as the sole generator of neutrino masses .    the rest of the paper is organized as follows . in section [ sec : generalities ]
we discuss the generalities of the brpv model , in particular those related with the neutral scalar sector . in section
[ sec : sneutrino - mix - analyr ] we derive formulas for brpv induced mixings . in section  [ sec : sneutrino - decays - analyticalr ] we write the relevant couplings for brpv sneutrino decays , give analytical formulas for the different partial decay widths , analyze the constraints on parameter space enforced by neutrino data and present our results . in section  [ sec : conclusions ] we summarize and present our conclusions .
in what follows we will briefly describe the main features of the bilinear r - parity breaking model , in particular those related with the neutral scalar sector .
we shall closely follow the notation used in @xcite and assume that all parameters are real . throughout the text matrices
will be denoted in bold - face .
in addition to the mssm r - parity conserving superpotential ( where we have suppressed relevant indices ) : @xmath23 the bilinear r - parity breaking model also contains the following terms : @xmath24 where @xmath25 is the @xmath26 completely antisymmetric tensor , @xmath27 runs over the sm fermion generations and @xmath20 is the r - parity and lepton number breaking bilinear parameter with units of mass .
consistency then requires a new set of soft susy breaking terms in the scalar potential , namely @xmath28 neglecting soft flavor mixing , the scalar potential relevant for neutral scalars is @xmath29 with @xmath30 and @xmath31 .
electroweak symmetry is broken once the higgs and slepton acquire a vev , @xmath32 and @xmath33 , with @xmath34  gev and @xmath35 .
the doublets are parameterized as @xmath36 with the neutral components given by @xmath37 here we introduce the notation @xmath38 to differentiate the cp - even sneutrinos from the cp - odd .    in the basis
@xmath39 the linear part of the neutral scalar potential can be written as @xmath40 where the @xmath41 s are the so - called tadpoles . at tree - level
these are ( see eqs .
( [ eq : full - ssb - scalar - potential ] ) and ( [ eq : vlinear ] ) ) @xmath42 where @xmath43 and @xmath44 .
the minimization of the potential , @xmath45 , requires the tadpoles to vanish .
thus , the vevs can be determined from the system of equations in ( [ eq : tree - lvel - tadpoles ] ) by imposing @xmath46 ( @xmath47 ) .
in particular , considering only leading order brpv terms , the sneutrino vevs can be written as @xmath48 where @xmath49 ( the tree level sneutrino mass ) and @xmath50 .
without conserved r - parity , there are no quantum numbers to distinguish the leptons and sleptons from the gauginos and higgsinos and higgs bosons respectively . therefore ,
in brpv several mixings between supersymmetric and non supersymmetric particles exist : ( @xmath22 ) neutralinos mix with neutrinos , ( @xmath51 ) charginos mix with charged leptons , ( @xmath52 ) higgs bosons mix with the sneutrinos and ( @xmath53 ) charged higss bosons mix with charged sleptons .
these mixings are important for calculating lsp decays , especially mixings of type ( @xmath22 ) as they allow to fix  via experimental neutrino data  the size of the brpv parameters .
since mixings of type ( @xmath53 ) are not of interest for sneutrino decays , in what follows we will only discuss analytical approximations for mixings of type ( @xmath22)-(@xmath52 ) .      in the basis
@xmath54 the neutral fermion mass matrix can be written as @xmath55 with @xmath56 where @xmath57 is the neutralino mass matrix : @xmath58 and @xmath59 and @xmath60 are the soft masses for the bino and wino respectively .
@xmath61 is the @xmath62 neutralino - neutrino mixing matrix given by @xmath63 in the weyl mass eigenstate basis , defined as . ]
@xmath64 the mass matrix becomes @xmath65 due to the smallness of the brpv parameters , at order @xmath20 , @xmath66 can be block diagonalized by decomposing the diagonalizing matrix @xmath67 as follows @xcite : @xmath68 the matrix @xmath69 block diagonalizes @xmath66 to the form @xmath70 , where @xmath71 is the tree - level light neutrino mass matrix
. the mixing parameters @xmath72 can thus be determined to be @xmath73 where @xmath74 ( with @xmath75 , @xmath76 being the weak mixing angle ) and @xmath77 finally the block diagonal mixing matrices @xmath78 and @xmath79 in ( [ eq : n - matrix ] ) diagonlize the neutralino and neutrino effective mass matrix , which reads @xmath80 since this matrix has two vanishing eigenvalues it can be diagonalized by only two rotation matrices , namely @xmath81 with @xmath82 for sneutrino decays the relevant part of the @xmath83 mixing turns out to be the @xmath84 block , that from eqs .
( [ eq : neutralino - neu - mixing - parameters ] ) and ( [ eq : leptonic - mix - mat ] ) can be written as @xcite @xmath85 where @xmath86 and @xmath87 with @xmath88 .
the coefficients @xmath89 are given by @xmath90 with @xmath91 .
taking into account eqs .
( [ eq : leptonic - mix - mat ] ) and ( [ eq : mixing - angles ] ) , explicitly @xmath92 ( @xmath93 ) and @xmath94 are given by @xmath95      in the bases @xmath96 the chargino charged lepton mass matrix is determined by the following lagrangian @xmath97 where @xmath98 and in the basis in which the charged lepton mass matrix is diagonal @xmath99 can be written as @xmath100 the block diagonal matrices correspond to the mssm chargino and charged lepton mass matrices whereas the off - diagonal mass matrices read @xmath101 defining the weyl mass eigenstates as . ]
@xmath102 the diagonal mass matrix @xmath103 is obtained through the biunitary transformation : @xmath104 approximate analytical expressions for the mixing matrices @xmath105 have been discussed in @xcite . here
we describe the method for completeness .
the off - diagonal block matrix @xmath106 , being proportional to the charged lepton yukawa couplings , can be neglected , and due to the smallness of the brpv parameters the mixing matrices @xmath105 can be written according to @xmath107 where , in first approximation in the brpv parameters , the matrices @xmath108 block - diagonalize the mass matrix , namely @xmath109 and the matrices @xmath110 diagonalize in turn the chargino mass matrix , with the rotation angles given by @xmath111 from equation ( [ eq : block - diag - mc ] ) the matrices @xmath112 are found to be @xmath113 which implies @xmath114 is suppressed with respect @xmath115 by a factor @xmath116 and thus can be neglected ( @xmath117 ) . explicitly @xmath115 can be written in terms of the brpv parameters and the coefficients entering in the chargino mass matrix : @xmath118      in the basis @xmath119 the mass matrix of the cp - even neutral scalars @xmath120 is determined by the following quadratic terms @xmath121    in what follows we will discuss approximate analytical formulas for the @xmath122 mixing . to our knowledge , this is the first time that they are explicitly given .
the entries of the mass matrix in ( [ eq : cp - even - mm ] ) involve the parameters @xmath123 , @xmath124 , the soft susy breaking coefficients @xmath125 and @xmath126 , and the r - parity breaking parameters @xmath20 , @xmath21 and @xmath127 .
we use the minimization conditions of the scalar potential ( @xmath128 ) to remove the parameters @xmath127 . in doing so
the matrix @xmath129 can be written in terms of the cp - odd neutral scalar mass @xmath130 , the sneutrino masses @xmath131 , @xmath132 and the brpv parameters @xmath20 and @xmath21 , namely @xmath133 where the following relations have been used @xmath134 note that phenomenologically consistency requires the inclusion of the one - loop correction in the @xmath135 entry , which in the limit of no stop mixing is @xmath136 see  @xcite for a more complete expression .    in the absence of brpv , @xmath129 corresponds to the neutral cp - even higgs mass matrix of the r - parity conserving mssm , which can be diagonalized via the rotation matrix @xmath137 the elements of the @xmath138 right - lower block sneutrino mass matrix in ( [ eq : cp - even - mm ] ) are given by @xmath139 finally , for the @xmath122 mixing @xmath140 block we have @xmath141 in the mass eigenstate basis defined as @xmath142 the lagrangian in ( [ eq : cp - even - mm ] ) becomes @xmath143 where @xmath144 is diagonal .
assuming real parameters @xmath145 can be parameterized as @xmath146 where the @xmath147 are @xmath148 rotation matrices .
if the @xmath149 mixing is small  as expected due to the smallness of the brpv parameters required by neutrino data  a perturbative diagonalization of the mass matrix in ( [ eq : cp - even - mm ] ) can be done . by neglecting the brpv parameters in @xmath129 and @xmath150
the rotation matrix reduces to @xmath151 when acting on @xmath152 the matrix @xmath153 diagonalizes the @xmath154 block @xmath129 according to @xmath155 ( where @xmath156 and @xmath157 are the light and heavy cp - even higgs bosons ) and modifies the brpv mixing matrix @xmath158 : @xmath159 for the mixing angle we have @xmath160 .
the matrices @xmath161 eliminate the first row entries in ( [ eq : brpv - mm - modified ] ) and up to @xmath162 leave the block - diagonal matrices @xmath163 and @xmath164 diagonal .
the rotation angles are given by @xmath165 the matrices @xmath166 instead eliminate the second row elements in ( [ eq : brpv - mm - modified ] ) leaving again , up to order @xmath167 , the block - diagonal matrices diagonal .
the rotation angles in this case read @xmath168 with these results at hand and neglecting terms of @xmath169 the full rotation matrix in ( [ eq : rot - mat - approx - s0 ] ) can be written as @xmath170 where @xmath171 and @xmath172 account for the sneutrino and cp - even neutral higgs components of the mass eigenstates , and read @xmath173
with the results of section  [ sec : sneutrino - mix - analyr ] we are now in a position to discuss approximate formulas for sneutrino decays . from now on we will consider only cp - even sneutrino decays , @xmath174 , and so will drop the superscript @xmath175 .
possible tree - level two - body sneutrino final states include fermionic modes @xmath176 , @xmath177 and @xmath178 ; electroweak gauge bosons modes @xmath179 and @xmath180 and higgs bosons modes @xmath181 , @xmath182 , @xmath183 and @xmath184 .    the goal of this paper and section is not a full study of the parameter space and all decays but rather to show that these heavy states can dominate the traditionally considered decay to @xmath19 and to identify the relevant parameter space for this dominance .
taking into account our approximate results for the chargino - right - handed lepton mixing , @xmath185 , mixing with right - handed leptons is zero ( @xmath186 , see eq .
( [ eq : xir - xil ] ) ) , the full lagrangian for sneutrino decays into charged leptons is reduced to  @xcite @xmath187 where the mixing matrices are given by eqs .
( [ eq : xil - explicitly ] ) and ( [ eq : full - rot - s0-final ] ) . due to the smallness of the brpv induced mixing we take @xmath188 .
note that while the first term in ( [ eq : charged - lepton - modes - lag ] ) necessarily leads to two same - flavor opposite - sign charged lepton final states the second term can lead to different - flavor signatures .
the corresponding feynman diagrams for these processes are depicted in figure  [ fig : sneutrino - decays](@xmath189,@xmath190 ) .
being proportional to sm charged lepton yukawa couplings these decays are dominated by final states involving @xmath191 s .
sneutrino decays are induced by sneutrino - higgs mixing while in figure  @xmath190 by chargino - charged lepton or neutralino - neutrino mixing , depending on whether the final states involves charged lepton or neutrinos .
the open circles with a cross inside indicate a brpv mixing insertion.__,width=377,height=132 ]    for sneutrino decays into neutrinos ( invisible decays @xmath192 ) , the relevant interactions are given by @xmath193 where the coupling reads @xmath194 where @xmath195 and @xmath196 the elements of the neutral fermion mixing matrix correspond to the entries of the matrix given in eq .
( [ eq : n - matrix ] ) .
these interactions induce invisible sneutrino decays as the ones shown in figure  [ fig : sneutrino - decays]@xmath190 .        in the up and down quark mass eigenstate basis
the sneutrino - quark - quark interactions are dictated by @xmath197 as for the charged lepton final states , these decays are controlled by sm quark yukawa couplings and thus are dominated by @xmath198 and @xmath15 , the last one if kinematically allowed .
the lagrangian for gauge boson final states is given by @xmath199 with @xmath200 and @xmath201 .
for higgs final states we write the lagrangian involving @xmath156 and @xmath157 and @xmath202 decays and so do not present the lagrangians that govern these interactions . ] : @xmath203 where the fully symmetric coupling @xmath204 is given by @xmath205\,,\ ] ] with @xmath206 and @xmath207 the interactions in ( [ eq : lagrangian - scalars ] ) for @xmath208 lead to decays of type @xmath189 in figure  [ fig : sneutrino - decays2 ] while for @xmath209 to those shown in figure  [ fig : sneutrino - decays2]@xmath190 .
fermionic final states are dominated by third generation quark and charged leptons . due to the structure of the higgs - sneutrino mixing , @xmath210 , @xmath211 and @xmath15 final states
are possible independently of the sneutrino flavor , whereas @xmath212 final states are only sizable for tau and mu sneutrinos . neglecting the final state masses , the partial decay widths for @xmath212 , @xmath210 and @xmath211 decays can be written as @xmath213 ^ 2\,,\\    \label{eq : bbar - ttbar - partial - widths2 }    \gamma(\tilde\nu_i\to b\bar b)&=    \frac{3 m_b^2\,g_f}{4\sqrt{2}\pi\,c_\beta^2 }    m_{\tilde\nu_i}(\boldsymbol{r^{s^0}})^2_{(i+2)1}\,,\end{aligned}\ ] ] where @xmath214 is the fermi constant . for invisible decay modes the partial decay width , summing over lepton flavors , can be written according to @xmath215 with @xmath216 given by ( [ invisible - coupling ] ) . for @xmath15 final states the phase space factors are relevant ,
accordingly the corresponding decay width reads @xmath217    for gauge boson final states the partial decay width is given by @xmath218 where @xmath219 , @xmath220 and the amplitude @xmath221 reads @xmath222 for @xmath223 final states the partial decay width is given by @xcite @xmath224 with @xmath225 and @xmath226 the integration variables @xmath227 lie in the ranges @xmath228 $ ] and @xmath229 $ ] .
the parameter @xmath230 ( @xmath231 being the total decay width of the gauge boson @xmath232 ) allows a smooth transition in the threshold region where the off - shell gauge boson becomes on - shell ( @xmath233 ) , the calculation of @xmath234 thus requires numerical integration over the variables @xmath227 in the transition region . outside that region ,
i.e. for @xmath235 it can be written according to @xcite @xmath236 where @xmath237 is a kinematic function given by @xmath238    finally for cp - even higgs bosons final states the partial decay width can be written as @xmath239 where @xmath240 , @xmath241 and the dimensionful coupling @xmath204 is given in ( [ eq : gijk - coupling ] ) . for the particular case @xmath242 ( @xmath14 and @xmath243 final states ) the width reduces to @xmath244      neutrino data plays an important role in brpv , since brpv parameters contribute to neutrino masses . therefore
, even if brpv is not solely responsible for neutrino masses , _
@xcite , one must not saturate neutrino masses via brpv
. therefore it would be worthwhile to discuss this correlation briefly here .
we proceed by assuming brpv as the sole source of neutrino masses and later consider relaxing this assumption .    at tree - level
, brpv allows for only one massive neutrino , as mentioned earlier .
one - loop contributions then leave only one generation massless . furthermore ,
for the approximations made in this paper , the tree - level mass is larger than the one - loop mass .
the upshot of all this is that only the so - called normal hierarchy of neutrino masses is allowed .
recent neutrino data  @xcite then specifies the following neutrino masses : @xmath245 in terms of the model parameters @xmath246 the experimental values of the atmospheric and reactor angles yield two additional constraints given by eq .
( [ eq : mixing - angles ] ) .
therefore , the atmospheric sector entirely fixes the @xmath247 parameters .
( @xmath27 ) parameters and sneutrino vevs as required for explaining neutrino data for the r - parity conserving parameters discussed in the text .
the horizontal ( vertical ) solid lines in the right - panel plot indicate the values where @xmath248 is dominated by @xmath249 ( @xmath250 ) , i.e. where @xmath251 ( @xmath252).__,title="fig:",width=283,height=245 ]   ( @xmath27 ) parameters and sneutrino vevs as required for explaining neutrino data for the r - parity conserving parameters discussed in the text .
the horizontal ( vertical ) solid lines in the right - panel plot indicate the values where @xmath248 is dominated by @xmath249 ( @xmath250 ) , i.e. where @xmath251 ( @xmath252).__,title="fig:",width=291,height=253 ]    constraints on the @xmath253 parameters arise from the solar sector . as long as the one - loop contribution to neutrino masses is smaller than the tree - level one ( an assumption we use throughout ) , @xmath254
is accurately determined by @xmath255 one - loop diagrams @xcite : @xmath256 where @xmath257 stands for sbottom mixing , @xmath258 are defined in eq .
( [ eq : epsilon - tilde - parameters ] ) and @xmath259 ( with @xmath260 is a scalar passarino - veltman function @xcite ) .
note that due to @xmath261 , @xmath262 moreover the solar mixing angle can be written as @xcite @xmath263 thus , eqs .
( [ eq : solar - mass - scale ] ) and ( [ eq : solar - mixing ] ) provide two constraints and determine , from eqs .
( [ eq : epsilon - tilde - parameters ] ) , @xmath264 ( @xmath265 or @xmath266 ) as a function of @xmath250 ( @xmath267 or @xmath268 ) . once the @xmath269 s and @xmath253 s are fixed the sneutrino vevs
are automatically fixed as well through @xmath270 see eq .
( [ eq : alignment - vector ] ) .
this in turn fixes all the relevant mixings for sneutrino decays ( @xmath271 , @xmath272 , @xmath273 ) once the r - parity conserving supersymmetric parameters are specified . with this knowledge in hand , we start to explore the consequences of the neutrino sector on the brpv parameters .
we note that our results have been verified using spheno @xcite ] .
figure  [ fig : epsilons - vevs ] shows typical values for @xmath264 and sneutrino vevs as a function of @xmath250 .
the plots were obtained by fixing @xmath274 , @xmath275 and @xmath276 where we have assumed a conservative lower bound on the sneutrino mass of 100  gev , see  @xcite for more details .
we will use this parameter point throughout the paper .
neutrino observable @xmath277 are varied in their 3@xmath278 experimental range @xcite .
we further assume that @xmath250 is positive . for these assumptions
we see from the left - hand side that @xmath279 and @xmath280 have a lower bound and from the right - handed side that this lower bound is larger than @xmath281 and @xmath282 respectively . on the other hand ,
since @xmath250 is undetermined , @xmath249 dominates it up to about @xmath283 . while allowing negative values for the parameters changes this picture , it leaves one important qualitative property the same : * only in one generation , @xmath284 , can the sneutrino vev , @xmath285 , be larger than the bilinear mixing term , @xmath286 , when brpv is the sole contributor to neutrino masses .
* this will have important consequences in the next section .    also of note in the left - panel plot
, @xmath250 obeys only an upper bound determined by the condition @xmath287 , which in terms of the bilinear r - parity violating parameters translates into @xmath288 .
in contrast the @xmath264 parameters , due to solar neutrino physics constraints , are forced to lie in a `` narrow '' range and are such that a region where @xmath289 exist .
consequently , while @xmath290 are mostly determined by @xmath264 , @xmath248 is controlled by @xmath249 in the region where @xmath291 , as demonstrated by the horizontal solid lines in figure  [ fig : epsilons - vevs ] ( right - panel ) which correspond to @xmath251 .    the bilinear r - parity breaking parameters selected as described above satisfy neutrino data , and thus lead to brpv models that can account for neutrino masses and mixings . however , it might be that these parameters are not sufficiently large to account for the neutrino mass scales . in that case
their contribution to the atmospheric and solar masses are still determined by eqs .
( [ eq : m - atmospheric ] ) and ( [ eq : solar - mass - scale ] ) but are such that , for example , @xmath292 and @xmath293 . figure [ fig : vevs - noatmospheric - nosolar ] shows the results for an illustrative case where the brpv model fits the atmospheric mass scale as well as the atmospheric and reactor angles in their @xmath294 experimental range , but the contribution to @xmath295 is subdominant . note that we plot only @xmath296 , as we have found that @xmath297 .
this result as well as @xmath298 in the region of small @xmath265 are due to the constraints arising from fitting @xmath299 and @xmath300
.     fits the experimental @xmath294 range but @xmath254 falls below the measured value ( see the text for more details).__,width=340,height=253 ]    in summary , if data is not explained by the brpv parameters all the @xmath253 s can be small mainly due to the absence of the solar data constraint , thus implying that , in this case , a region where the three sneutrino vevs are large ( @xmath301 ) exist .
we are finally ready to address the main aim of this paper : showing that sneutrino lsp decays into heavy sm pairs , invisible modes and different lepton flavor final states _ i.e. _ @xmath302 ; @xmath16 and @xmath303 ( @xmath304 ) can dominate over the traditionally considered @xmath211 mode , which was typically considered as the dominate mode  @xcite .
we begin with a very general study ignoring all neutrino constrains ( susy parameters can always be chosen in such a way so that brpv does not saturate the neutrino masses ) and exploring in the @xmath20-@xmath21 parameter space where the aforementioned modes dominate the @xmath19 decays of @xmath305 . without neutrino constraints , these decay properties are independent of generation .
we scan over the following values : @xmath306     and @xmath16 with the @xmath211 final state .
it has been assumed that brpv does not saturate neutrino masses .
blue points indicate where @xmath307 ( left - hand side plot ) or @xmath308 ( right - hand side plot ) while the red points the opposite .
the invisible mode involves a sum over final state lepton flavors , and so its importance is determined by the relative size between @xmath20 and any of the three sneutrino vevs .
the vertical axes in the right - hand side plot corresponds to the largest sneutrino vev obtained when randomly scanning the parameter space .
these plots therefore show an important general result : the @xmath309 and invisible modes have larger branching fractions than the @xmath211 mode only when the sneutrino vev is much larger than its corresponding bilinear @xmath253 parameter.__,title="fig:",width=283,height=245 ]   and @xmath16 with the @xmath211 final state .
it has been assumed that brpv does not saturate neutrino masses .
blue points indicate where @xmath307 ( left - hand side plot ) or @xmath308 ( right - hand side plot ) while the red points the opposite .
the invisible mode involves a sum over final state lepton flavors , and so its importance is determined by the relative size between @xmath20 and any of the three sneutrino vevs .
the vertical axes in the right - hand side plot corresponds to the largest sneutrino vev obtained when randomly scanning the parameter space .
these plots therefore show an important general result : the @xmath309 and invisible modes have larger branching fractions than the @xmath211 mode only when the sneutrino vev is much larger than its corresponding bilinear @xmath253 parameter.__,title="fig:",width=283,height=245 ]    and @xmath123 below 1000  gev but large enough so that the sneutrino is the lsp .
we also use @xmath310  gev as suggested by recent lhc results . for @xmath12 and @xmath16
modes the results are displayed in figure  [ fig : sneu - ww - bb ] where blue dots indicate the points at which these channels have partial decay widths much larger than the @xmath211 mode , while the red points show the opposite .
it is striking that the ratio of the partial widths are relatively independent of r - parity conserving susy parameters and depend only on the brpv parameters . from this figure
, one can conclude that the @xmath309 and the @xmath16 partial widths dominate the @xmath19 width only when the sneutrino vev is much larger than the bilinear mixing parameter and , as will be shown below , this condition holds for the other heavy sm final states .
couple this with the main result of the previous section : satisfying neutrino masses solely through brpv means that @xmath311 can only be satisfied in one generation , indicates that heavy final states as well as invisible modes can only have branching fractions larger than the @xmath211 mode for one generation of sneutrinos .
we now turn to the discussion of the different lepton flavor final states .
these modes ( @xmath312 ) , being induced by chargino charged lepton mixing , are mainly controlled by the bilinear mixing parameters for the flavor of the final lepton which differs from the sneutrino flavor ( see eq .
( [ eq : xil - explicitly ] ) and ( [ eq : bbar - ttbar - partial - widths ] ) ) _ i.e. _ @xmath313 not @xmath314 . the @xmath211 mode , on the other hand ,
is governed by @xmath314 ( same flavor as the sneutrino ) if @xmath315 , or @xmath316 , if @xmath317 .    therefore , when @xmath315 ( region where @xmath318 modes dominate @xmath319 ) , neutrino constraints imply all three @xmath253s are of the same order ( see fig .
[ fig : epsilons - vevs ] ) .
therefore the different lepton flavor final states become subdominant to the @xmath318 modes due to charged lepton yukawa suppression . in contrast , when @xmath317 ( @xmath319 modes dominate @xmath318 modes ) the different lepton flavor final states can have larger partial widths than the @xmath211 mode provided @xmath320 ^ 2 $ ] exceeds the corresponding yukawa suppression @xmath321 , which turns out to be possible only for second and third generation sneutrinos as shown in fig .
[ fig : epsilons - vevs ] , where it can be seen that @xmath320 ^ 2 $ ] can be at most @xmath322 ( note that @xmath323 , @xmath324 while @xmath325 ) . as in the case of heavy sm and invisible modes , these different lepton flavor final states
can dominant over the @xmath211 final state only in a single sneutrino generation .
we continue by considering the specific susy point discussed in the last section and examining first and third generation sneutrino decays into @xmath19 , @xmath326 , @xmath327 and different lepton flavor modes assuming that brpv is solely responsible for neutrino masses .
we have also calculated the partial decay widths for the @xmath14 and @xmath13 final states , finding that they have a similar dependence on @xmath328 and are of the same order as the @xmath12 mode ( the @xmath329 channel is suppressed due to the off - shell top for this sneutrino mass ) .
the results are displayed in figure  [ fig : width - snu1-snu3 ] versus @xmath328 , where it can be seen that @xmath330 exhibits the usual bahavior quoted in the literature , that is to say , a dominant @xmath211 mode followed by a slightly yukawa suppressed @xmath331 mode and strongly suppressed @xmath16 , @xmath319 , @xmath332 and @xmath333 final states .
the huge supression of the different lepton flavor mode , @xmath334 , in the small @xmath250 region is due to the combined effect of the electron yukawa suppression and the smallness of @xmath335 .
the @xmath332 final state being controlled by @xmath336 involves only the electron yukawa suppression , and so it has a larger partial decay width than the @xmath337 mode .
the situation for @xmath338 , however , is quite different . in the region of large @xmath250
it features a @xmath211 dominant mode and apart from the @xmath339 mode it behaves like @xmath330 . in the region of small @xmath250 , in contrast , the `` non - conventional '' modes @xmath340 ( @xmath341 ) , @xmath16 , @xmath12 , @xmath342 and @xmath343 have the largest partial widths by far , greatly surpassing the @xmath211 final state .
in particular , the different lepton flavor channels exceed this mode by more than two - orders of magnitude , thus drastically changing the sneutrino lsp phenomenology .
some words are in order as regards the size of the @xmath331 partial width in @xmath338 decays .
this decay mode is a unequivocal signal of brpv models , but its size is strongly dependent on the nature of the bilinear parameters . in r - parity violating models satisfying the conditions @xmath344 and @xmath345 at certain high energy scale , the bilinears can be rotated away .
they are , however , subsequently induced at the tev scale by rge running effects @xcite therefore inducing the @xmath331 mode , but with a quite suppressed partial width ( typically several orders of magnitude smaller than that of the @xmath211 channel ) . in a mssm based analysis , as the one performed here , where all the brpv parameters are treated as free parameters subject only to phenomenological constraints this is not necessarily the case , as demonstrated in the right - hand side plot in fig .
[ fig : width - snu1-snu3 ] .
( left ) and @xmath346 ( right ) to @xmath211 , @xmath319 , @xmath331 , @xmath327 , @xmath347 , @xmath348 and @xmath349 final states as a function of @xmath350 .
brpv parameters have been fixed by using neutrino oscillation data ( see the text for more details ) .
the remaining heavy sm modes , @xmath351 and @xmath343 , have the same behavior and are of the same order than the @xmath319 mode ( are not shown to avoid crowding the plots ) .
decays of @xmath352 are not shown since they are very similar to those of @xmath353.__,title="fig:",width=283,height=245 ]   ( left ) and @xmath346 ( right ) to @xmath211 , @xmath319 , @xmath331 , @xmath327 , @xmath347 , @xmath348 and @xmath349 final states as a function of @xmath350 .
brpv parameters have been fixed by using neutrino oscillation data ( see the text for more details ) .
the remaining heavy sm modes , @xmath351 and @xmath343 , have the same behavior and are of the same order than the @xmath319 mode ( are not shown to avoid crowding the plots ) .
decays of @xmath352 are not shown since they are very similar to those of @xmath353.__,title="fig:",width=283,height=253 ]    figure  [ fig : width - snu1-snu3 ] reinforces the arguments given in the previous paragraphs . reiterating : even if bilinear r - parity violation accounts for neutrino data , the @xmath211 channel is not necessarily the dominant decay mode .
we have checked that for third and second generation sneutrino lsps the different lepton flavor final states dominate ( @xmath341 , @xmath354 for @xmath338 and @xmath355 , @xmath356 for @xmath352 ) while first generation sneutrinos decays are driven by invisible and heavy sm modes ( @xmath351 and @xmath319 ) . as already pointed out and exemplified for @xmath338 in fig .
[ fig : width - snu1-snu3 ] , these decays can dominate only in a single sneutrino generation .
so , observing at least two different sneutrino flavors decaying to the aforementioned modes ( something that turns out to be viable if the sneutrino mass splittings do not allow the heavier sneutrinos to decay into the lsp ) will prove that brpv does not account for neutrino data , because this would require @xmath357 for those generations , which is not possible if neutrino physics is determined solely by brpv , see fig .
[ fig : vevs - noatmospheric - nosolar ] .
focusing further on those final states of @xmath358 decays ( excluding the @xmath343 mode since it behaves like the @xmath319 mode and is smaller by an order of 1/2 ) , we study @xmath358 decays as a function of @xmath359 in figure  [ fig : brxxoverbb - mnu ] allowing a @xmath358 heavy enough so that the @xmath329 channel is also open .
all the r - parity conserving parameters have been fixed according to the supersymmetric point used in the previous section .
the r - parity breaking parameters have been fixed by adjusting neutrino observables to their best fit point values @xcite and @xmath360 .
we have checked that the values are insensitive to changes in @xmath124 in the range [ 2,30 ] .
off - shell calculations of the higgs and top states have been neglected since the different lepton flavor final states dominant over the full sneutrino mass range .
both plots support our claim : even if brpv is the sole generator of neutrino masses a sneutrino lsp does not necessarily decay predominantly into @xmath211 final states .     and @xmath340 ) , invisible ( @xmath327 ) , @xmath319 , @xmath351 and @xmath15 final states normalized to @xmath361 as a function of @xmath362 .
right plot : decay branching ratios for dominant ( excluding @xmath343 ) and @xmath211 modes as a function of @xmath362 .
the branching ratios have been calculated in the small @xmath250 region ( @xmath363 ) and for bilinear r - parity breaking parameters fixed by neutrino parameters according to their best fit point values @xcite .
both plots show the main result of this paper : sneutrino lsp decays are not necessarily dominated by the @xmath211 mode.__,title="fig:",width=283,height=245 ]   and @xmath340 ) , invisible ( @xmath327 ) , @xmath319 , @xmath351 and @xmath15 final states normalized to @xmath361 as a function of @xmath362 .
right plot : decay branching ratios for dominant ( excluding @xmath343 ) and @xmath211 modes as a function of @xmath362 .
the branching ratios have been calculated in the small @xmath250 region ( @xmath363 ) and for bilinear r - parity breaking parameters fixed by neutrino parameters according to their best fit point values @xcite .
both plots show the main result of this paper : sneutrino lsp decays are not necessarily dominated by the @xmath211 mode.__,title="fig:",width=283,height=245 ]    one can extrapolate from figure 6 the behavior of the r - parity violating sneutrino decays for the first two generations , in the analogous region @xmath364 . in these cases ,
only the leptonic different flavor states shift and they do so by a factor of approximately @xmath365 for @xmath366 and @xmath367 for @xmath368 .
therefore , @xmath352 decays will still be dominated by different flavor leptonic final states , while @xmath369 will decay predominately into heavy sm states and neutrinos .
relaxing the assumption that bilinear r - parity violation is the sole generator of neutrino masses ( as in the second case discussed in section  [ sec : constraints - neu - physics ] ) will not drastically change the results in this section .
this is because the relative importance of the ratios @xmath370 is determined by @xmath371 .
however , a quantity that becomes relevant in this case is the sneutrino decay length @xmath372 , since the size of the bilinear r - parity breaking parameters can be small enough to cause the sneutrino to decay outside of the detector . neglecting the lorentz boost factor ,
we have found that as long as the r - parity violating parameters fit the neutrino data , in general @xmath372 is well below 1 mm ( inline with  @xcite ) .
we have also checked that if @xmath373 the decay length is generically below @xmath374 cm .
note that values of @xmath375 in the range @xmath376 $ ] might be in conflict with neutrino data when the contributions from the mechanism responsible for neutrino masses are taken into account ( for example a standard seesaw ) .
once the contributions to the atmospheric and solar mass scales fall below @xmath377 the decay length can eventually exceed 1 m almost independently of the sneutrino mass ( as expected due to the large partial decay widths the different lepton flavor modes have ) . figure
[ fig : decay - length ] shows the decay length for @xmath378 .
, where @xmath379 have been taken according to their best fit point values @xcite.__,width=340,height=253 ]    thus far , this section has focused on sneutrino decays . at this point ,
we briefly comment on production . if the sneutrino is indeed the lsp , all produced supersymmetric particles produced will cascade decay down to it .
the dominate production depends on the relative masses , with strong production  @xcite dominating unless the gluinos and squarks are too heavy , in which case relatively light sneutrinos can still be pair produced through electroweak processes .
regardless , the parameter region explored in this paper yields compelling signals since the dominate different flavor leptonic final states would have little background if the intermediate sneutrino can be reconstructed . meanwhile
, the @xmath211 mode corresponding to the parameter space typically considered in the literature would suffer from larger qcd background and would therefore be more difficult to detect . furthermore ,
if the brpv contributions to neutrino masses are subdominant to some other mechanism , sneutrino decays could involve large displaced vertices ( see fig .
[ fig : decay - length ] ) thus rendering its identification even more plausible .
brpv is a good effective theory for models of spontaneous r - parity violation , in which proton decay as well as the number of new parameters is under control compared to explicit r - parity violation . in the presence of brpv
, sneutrino lsp decays can in some cases resemble a heavy higgs .
however , as we have discussed , sneutrinos ( third and second generation ) can dominantly decay to different - flavor charged leptons final states or heavy sm modes ( first generation ) , thus completely altering the `` conventional '' sneutrino lsp phenomenology driven by the @xmath211 mode .
we have studied these decays into @xmath12 , @xmath343 , @xmath14 , @xmath15 , @xmath16 and @xmath303 ( @xmath304 ) and found that long as @xmath317 , they have larger branching ratios than the @xmath198 mode . as discussed in section  [ sec : constraints - neu - physics ] , for models where the brpv parameters are fixed by neutrino data , this is possible only for a single sneutrino flavor , due to the constraints from the solar sector .
therefore , if 2 or more sneutrino generations are degenerate enough to decay via brpv into these states one could rule out brpv as the sole source of neutrino masses .    for models where the bilinear r - parity breaking parameters do not contribute significantly to neutrino masses ,
the constraint @xmath380 for all flavors is viable .
accordingly , in these brpv models large sneutrino lsp branching fractions into either @xmath303 or @xmath12 , @xmath343 , @xmath14 and @xmath15 are possible . in this case , however , special attention has to be paid to the sneutrino decay length @xmath372 .
we calculated @xmath372 , finding that as long as the brpv contributions to neutrino masses are not below @xmath381 the decay length is generically below @xmath374 cm .
regardless of these considerations , sneutrino decays into different - flavor charged leptons or heavy sm final states ( @xmath303 or @xmath12 , @xmath343 , @xmath14 and @xmath15 ) are an interesting phenomenological possibility that has not received much attention before . because brpv only affects the decays of the lsp , assuming a sneutrino lsp means that susy events could cascade decay to either different flavor charged leptons or two pairs of heavy sm states , which constitutes a unique and novel susy signature .
we want to thank martin hirsch for his comments in the early stage of this project and werner porod for spheno support .
ss would also like to thank pavel fileviez perez for comments .
das is supported by a belgian fnrs postdoctoral fellowship .
dr has been partially supported by sostenibilidad - udea and colciencias through the grant number 111556934918 .
99 g.  aad _ et al . _ [ atlas collaboration ] , phys .
b * 716 * ( 2012 ) 1 [ arxiv:1207.7214 [ hep - ex ] ] .
s.  chatrchyan _ et al . _
[ cms collaboration ] , phys .
b * 716 * ( 2012 ) 30 [ arxiv:1207.7235 [ hep - ex ] ] .
r.  barbier , c.  berat , m.  besancon , m.  chemtob , a.  deandrea , e.  dudas , p.  fayet and s.  lavignac _ et al .
_ , phys .
rept .   * 420 * , 1 ( 2005 ) [ hep - ph/0406039 ] .
s.  chatrchyan _ et al . _
[ cms collaboration ] , arxiv:1206.3949 [ hep - ex ] .
g.  aad _ et al . _ [ atlas collaboration ] , phys .
b * 709 * , 137 ( 2012 ) [ arxiv:1110.6189 [ hep - ex ] ] .
l.  m.  carpenter , d.  e.  kaplan and e.  -j .
rhee , phys .
lett .   * 99 * , 211801 ( 2007 ) [ hep - ph/0607204 ] .
l.  m.  carpenter , d.  e.  kaplan and e.  j.  rhee , arxiv:0804.1581 [ hep - ph ] .
p.  w.  graham , d.  e.  kaplan , s.  rajendran and p.  saraswat , jhep * 1207 * , 149 ( 2012 ) [ arxiv:1204.6038 [ hep - ph ] ] . c.  s.  aulakh and r.  n.  mohapatra , phys .
b * 119 * ( 1982 ) 136 .
m.  j.  hayashi and a.  murayama , phys .
b * 153 * ( 1985 ) 251 .
r.  n.  mohapatra , phys .
* 56 * ( 1986 ) 561 .
j.  m.  mira , e.  nardi , d.  a.  restrepo and j.  w.  f.  valle , phys .
b * 492 * , 81 ( 2000 ) [ hep - ph/0007266 ] .
d.  aristizabal sierra , d.  restrepo and o.  zapata , phys .
d * 80 * , 055010 ( 2009 ) [ arxiv:0907.0682 [ hep - ph ] ] . c.  s.  aulakh , a.  melfo , a.  rasin and g.  senjanovic , phys .
b * 459 * ( 1999 ) 557 [ hep - ph/9902409 ] . c.  s.  aulakh , b.  bajc , a.  melfo , a.  rasin and g.  senjanovic , nucl .
b * 597 * ( 2001 ) 89 [ hep - ph/0004031 ] .
p.  fileviez perez and s.  spinner , phys .
b * 673 * ( 2009 ) 251 [ arxiv:0811.3424 [ hep - ph ] ] .
v.  barger , p.  fileviez perez and s.  spinner , phys .
* 102 * ( 2009 ) 181802 [ arxiv:0812.3661 [ hep - ph ] ] . m.  ambroso and b.  ovrut , jhep * 0910 * ( 2009 ) 011 [ arxiv:0904.4509 [ hep - th ] ] ; m.  ambroso and b.  a.  ovrut , int .
j.  mod .
phys .  a * 25 * ( 2010 ) 2631 [ arxiv:0910.1129 [ hep - th ] ] .
p.  fileviez perez and s.  spinner , phys .
d * 83 * ( 2011 ) 035004 [ arxiv:1005.4930 [ hep - ph ] ] . a.  masiero and j.  w.  f.  valle , phys .
b * 251 * , 273 ( 1990 ) ; w.  buchmuller , l.  covi , k.  hamaguchi , a.  ibarra and t.  yanagida , jhep * 0703 * , 037 ( 2007 ) [ hep - ph/0702184 [ hep - ph ] ] . f.  takayama and m.  yamaguchi , `` gravitino
dark matter without r - parity , '' phys .
b * 485 * , 388 ( 2000 ) [ hep - ph/0005214 ] .
s.  dimopoulos , d.  eichler , r.  esmailzadeh and g.  d.  starkman , phys .
d * 41 * , 2388 ( 1990 ) .
a.  bartl , m.  hirsch , t.  kernreiter , w.  porod and j.  w.  f.  valle , jhep * 0311 * , 005 ( 2003 ) [ hep - ph/0306071 ] .
d.  aristizabal sierra , m.  hirsch and w.  porod , jhep * 0509 * , 033 ( 2005 ) [ hep - ph/0409241 ] .
w.  porod , m.  hirsch , j.  romao and j.  w.  f.  valle , phys .
d * 63 * , 115004 ( 2001 ) [ hep - ph/0011248 ] .
m.  hirsch and w.  porod , phys .
d * 68 * , 115007 ( 2003 ) [ hep - ph/0307364 ] .
s.  roy , b.  mukhopadhyaya and , phys .
d * 55 * , 7020 ( 1997 ) [ hep - ph/9612447 ] .
b.  mukhopadhyaya , s.  roy , f.  vissani and , phys .
b * 443 * , 191 ( 1998 ) [ hep - ph/9808265 ] .
s.  bar - shalom , g.  eilam and b.  mele , phys .
d * 64 * , 095008 ( 2001 ) [ hep - ph/0106053 ] .
y.  grossman , h.  e.  haber and , phys .
d * 59 * , 093008 ( 1999 ) [ hep - ph/9810536 ] .
m.  hirsch , m.  a.  diaz , w.  porod , j.  c.  romao and j.  w.  f.  valle , phys .
d * 62 * , 113008 ( 2000 ) [ erratum - ibid .
d * 65 * , 119901 ( 2002 ) ] [ hep - ph/0004115 ] .
m.  a.  diaz , m.  hirsch , w.  porod , j.  c.  romao and j.  w.  f.  valle , phys .
d * 68 * , 013009 ( 2003 ) [ erratum - ibid .
d * 71 * , 059904 ( 2005 ) ] [ hep - ph/0302021 ] .
m.  nowakowski and a.  pilaftsis , nucl .
b * 461 * , 19 ( 1996 ) [ hep - ph/9508271 ] .
m.  hirsch and j.  w.  f.  valle , nucl .
b * 557 * , 60 ( 1999 ) [ hep - ph/9812463 ] .
m.  drees and m.  m.  nojiri , phys .
d * 45 * , 2482 ( 1992 ) .
f.  de campos , m.  a.  garcia - jareno , a.  s.  joshipura , j.  rosiek , j.  w.  f.  valle , nucl .
b * 451 * , 3 ( 1995 ) [ hep - ph/9502237 ] .
a.  djouadi , j.  kalinowski and p.  m.  zerwas , z.  phys .
c * 70 * , 435 ( 1996 ) [ hep - ph/9511342 ] .
a.  djouadi , phys .
rept .   * 457 * , 1 ( 2008 ) [ hep - ph/0503172 ] .
e.  ma , mod .
a * 17 * , 1259 ( 2002 ) [ hep - ph/0205025 ] .
d.  aristizabal sierra , m.  hirsch , j.  w.  f.  valle and a.  villanova del moral , phys .
d * 68 * , 033006 ( 2003 ) [ hep - ph/0304141 ] .
d.  k.  ghosh , g.  senjanovic and y.  zhang , phys .
b * 698 * , 420 ( 2011 ) [ arxiv:1010.3968 [ hep - ph ] ] .
v.  barger , p.  fileviez perez and s.  spinner , phys .
b * 696 * , 509 ( 2011 ) [ arxiv:1010.4023 [ hep - ph ] ] .
d.  v.  forero , m.  tortola and j.  w.  f.  valle , phys .
d * 86 * , 073012 ( 2012 ) [ arxiv:1205.4018 [ hep - ph ] ] .
m.  c.  gonzalez - garcia , m.  maltoni , j.  salvado and t.  schwetz , arxiv:1209.3023 [ hep - ph ] .
g.  l.  fogli , e.  lisi , a.  marrone , d.  montanino , a.  palazzo and a.  m.  rotunno , phys .
d * 86 * , 013012 ( 2012 ) [ arxiv:1205.5254 [ hep - ph ] ] .
g.  passarino and m.  j.  g.  veltman , nucl .
b * 160 * , 151 ( 1979 ) .
w.  porod , comput .
commun .   *
153 * , 275 ( 2003 ) [ hep - ph/0301101 ] ; w.  porod and f.  staub , comput .
commun .
* 183 * , 2458 ( 2012 ) [ arxiv:1104.1573 [ hep - ph ] ] .
j.  abdallah _ et al .
_ [ delphi collaboration ] , eur .
j.  c * 31 * , 421 ( 2003 ) [ hep - ex/0311019 ] .
b.  de carlos and p.  l.  white , phys .
d * 55 * , 4222 ( 1997 ) [ hep - ph/9609443 ] .
b.  de carlos and p.  l.  white , phys .  rev .
d * 54 * , 3427 ( 1996 ) [ hep - ph/9602381 ] .
e.  nardi , phys .
d * 55 * , 5772 ( 1997 ) [ hep - ph/9610540 ] .
m.  kramer , a.  kulesza , r.  van der leeuw , m.  mangano , s.  padhi , t.  plehn , x.  portell and , arxiv:1206.2892 [ hep - ph ] .