with the discovery of the standard model ( sm ) higgs - like scalar at the lhc @xcite , the sm particle content seems complete . in particular , the mass and couplings of the neutral higgs boson seem to disfavor an additional chiral generation of quarks and leptons @xcite .
however , additional vector - like fermions , in which an sm generation is paired with another one of opposite chirality , and with identical couplings , are less constrained , as there is no quadratic contribution to their masses .
these states appear as natural extensions of the sm particle content , in theories with warped or universal extra dimensions , as kaluza - klein excitations of the bulk field @xcite , in non - minimal supersymmetric extensions of the sm @xcite , in composite higgs models @xcite , in little higgs model @xcite and in gauged flavor groups @xcite .
vector fermions have identical left- and right - handed couplings and can have masses which are not related to their couplings to the higgs bosons @xcite .
depending on the dominant decay mode , the limits on new vector - like fermions range from @xmath2 gev @xcite , rendering them observable at the lhc .
vector quarks can modify both the production and decays of the higgs boson at the lhc , while vector leptons do not carry @xmath3 charge and can only modify the decay patterns of the higgs .
study of the latter would be a sensitive probe for new physics .
the lepton states contribute to self - energy diagrams for electroweak gauge boson masses and precision observables , and consistent limits on their masses and mixings have been obtained @xcite .
vector leptons have been studied in the context of the sm @xcite , but less so for models beyond the sm , where they can also significantly alter the phenomenology of the model . in the sm ,
introducing heavy fermions provides a contribution of the same magnitude and sign to the one of the top quark and interferes destructively with the dominant @xmath4-contribution , reducing the @xmath5 rate with respect to its sm value .
recent studies indicate that cancellations between scalar and fermionic contributions allow a wide range of yukawa and mass mixings among vector states @xcite .
an investigation of vector leptons in the two higgs doublet model @xcite showed that the presence of additional higgs bosons alleviates electroweak precision constraints .
introducing vector leptons into supersymmetry @xcite can improve vacuum stability and enhance the di - photon rate by as much as 50% , while keeping new particle masses above 100 gev and preserving vacuum stability conditions .
in the present work , we investigate vector leptons in the context of the higgs triplet model ( htm ) .
we do not deal with lhc phenomenology ( pair production and decay ) of the extra leptons , which has been discussed extensively in the literature @xcite , choosing instead to focus on signature features of this model .
we have previously shown that in the higgs triplet model , enhancement of the @xmath6 rate is possible only for the case where the doublet and triplet neutral higgs fields mix considerably @xcite .
we extend our analysis to include additional vector - like leptons in the model and investigate how these affect the higgs di - photon decay rate , with or without significant mixing in the neutral higgs sector .
originally , both cms and atlas experiments at lhc observed an enhancement of the higgs di - photon rate , while the di - boson rates ( @xmath7 ) have been roughly consistent with sm expectations . at present
cms observes @xmath8 times the sm rate , while atlas observes @xmath9 times the sm rate @xcite .
given these numbers , it is possible that either the sm value will be proven correct , or a modest enhancement will persist .
a further test of the sm is the correlation of the decay @xmath1 with one for @xmath10 .
we also include the prediction for the vector lepton effect on branching ratio of @xmath1 and comment on the relationship with the di - photon decay .
we have an additional reason to investigate the effects of vector leptons in the higgs triplet model .
the model includes doubly - charged higgs bosons , predicted by most models to be light .
being pair - produced , the doubly - charged higgs bosons are assumed to decay into a pair of leptons with the same electric charge , through majorana - type interactions @xcite . assuming a branching fraction of 100% decays into leptons , _
i.e. _ , neglecting possible decays into @xmath4-boson pairs , the doubly - charged higgs mass has been constrained to be larger than about 450 gev , or more , depending on the decay channel . however ,
if the vector leptons are light enough , which they can be , the doubly - charged higgs bosons can decay into them and thus evade the present collider bounds on their masses .
we investigate this possibility in the second part of this work .
our work is organized as follows .
we introduce the higgs triplet model with vector leptons in section [ sec : model ] . in section [ sec : neutral ] we analyze the effect of the vector leptons on the decays of the neutral higgs bosons , and discuss the constraints on the parameter space coming from requiring agreement with present lhc data .
we include both loop - dominated decays , @xmath0 in [ subsec : gamgam ] , and @xmath11 in [ subsec : zgam ] . in section [ sec :
doublycharged ] we analyze the effect of the vector leptons on the production and decay mechanisms of the doubly - charged higgs at the lhc .
we summarize our findings and conclude in section [ sec : conclusion ] .
the electroweak gauge symmetry is broken by the vevs of the neutral components of the doublet and triplet higgs fields , @xmath41 where @xmath42 and @xmath43 are the doublet higgs field and the triplet higgs field , with @xmath44 ( 246 gev)@xmath45 .
higgs masses and coupling constants in the presence of non - trivial mixing in the neutral sector have been obtained previously @xcite .
one can invoke new symmetries to restrict the interaction of the vector leptons with each other , or with the ordinary leptons , or disallow the presence of bare mass terms in the lagrangian .
for instance ,    1 .   if an additional @xmath46 symmetry under which the primed , double primed and the ordinary leptons have different charges , this would forbid explicit masses @xmath47 , @xmath48 , @xmath49 and @xmath50 in the lagrangian
. vector leptons would get masses only through couplings to the higgs doublet fields @xcite .
2 .   if one imposes a symmetry under which all the new @xmath51 singlet fields are odd , while the new @xmath51 doublets are even , this forces all yukawa couplings involving new leptons to vanish , @xmath52 , and the masses arise only from explicit terms in the lagrangian @xcite .
3 .   finally one can impose a new parity symmetry which disallows mixing between the ordinary leptons and the new vector lepton fields , under which all the new fields are odd , while the ordinary leptons are even @xcite , i.e. , such that @xmath53 ; alternatively one might choose these couplings to be very small .    in this analysis the focus will be on higgs decays .
we investigate the model subjected to symmetry conditions 1 and/or 2 ; and we neglect mixing between the ordinary and the new vector leptons . when allowed , stringent constraints exist on the masses and couplings with ordinary leptons .
new vector leptons are ruled out when they mediate flavor - changing neutral currents processes , generate sm neutrino masses or contribute to neutrinoless double beta decay .
recent studies of models which allow mixing between the ordinary leptons and the new ones exist @xcite , but there restrictions from lepton - flavor violating decays force the new leptons to be very heavy @xmath54 tev , or reduces the branching ratio for @xmath55 and of @xmath0 decay to 30 - 40% of the sm prediction , neither desirable features for our purpose here . in the higgs triplet model ,
distinguishing signals would be provided by lighter vector leptons . since imposing no mixing between ordinary and new leptons
allows new lepton masses to be as light as @xmath56 gev , perhaps without a reduction in the higgs di - photon branching ratio , we investigate this scenario here .    in the charged sector ,
the @xmath57 mass matrix @xmath58 is defined as @xcite @xmath59 with @xmath60 and @xmath61 , with @xmath62 the vev of the neutral component of the higgs doublet .
the mass matrix can be diagonalized as follows : @xmath63 the mass eigenvalues are @xmath64 ,   \label{eq : eigenvalues}\ ] ] with @xmath65 by convention , @xmath66 . for simplicity
we assume lepton yukawa couplings @xmath67 and @xmath68 are real so that the transformations that diagonalize the mass matrix are real orthogonal matrices : @xmath69 the angles @xmath70 are given by @xmath71 the eigenstates of the vector lepton mixing matrix enter in the evaluation of @xmath0 and @xmath1 in the next section .
the presence of the vector leptons affects the loop - dominated decays of the neutral higgs , @xmath0 and @xmath1 , and possible relationships between them . in the higgs triplet model , singly- and
doubly - charged bosons also enter in the loops , creating a different dynamic than in the sm .
we analyze these decays in turn , and look for possible correlations between them .
recently , the triplet higgs model has received renewed interest recently because of attempts to reconcile the excess of events in @xmath0 observed at the lhc over those predicted by the sm .
such an enhancement hints at the presence of additional particles , singlets under @xmath3 , but charged under @xmath72 which affect only the loop - dominated decay branching ratio , while leaving the production cross section and tree level decays largely unchanged .
vector leptons are prime candidates for such particles , so we study their contribution to the higgs decay branching . the decay width @xmath0 is @xmath73_{htm } & = & \frac{g_f\alpha^2 m_{h}^3 } { 128\sqrt{2}\pi^3 } \bigg| \sum_{f } n^f_c q_f^2 g_{h ff }   a_{1/2 } ( \tau^h_f ) + g_{h ww } a_1 ( \tau^h_w ) + \tilde{g}_{h h^\pm\,h^\mp } a_0(\tau^h_{h^{\pm } } ) \nonumber \\ & + &    4 \tilde{g}_{h h^{\pm\pm}h^{\mp\mp } } a_0(\tau^h_{h^{\pm\pm } } )    + \frac{\mu_{e_1}g_{hff}}{m_{e_1 } } a_{1/2}(\tau^h_{e_1 } ) + \frac{\mu_{e_2}g_{hff}}{m_{e_2 } } a_{1/2}(\tau^h_{e_2 } )        \bigg|^2 \ , , \label{partial_width_htm}\end{aligned}\ ] ] with @xmath74 given in eq . ( [ eq : eigenvalues ] ) , and where the loop functions for spin @xmath75 , spin @xmath76 and spin @xmath77 are given by : @xmath78\ , \tau^{-2 } \ , , \label{eq : ascalar}\\ a_{1/2}(\tau)&=&-\tau^{-1}\left[1+\left(1-\tau^{-1}\right)f\left(\tau^{-1}\right)\right ] , \label{eq : afermion}\\ a_1(\tau)&= & 1+\frac32 \tau^{-1}+4 \tau^{-1}\left(1-\frac12 \tau^{-1}\right)f\left(\tau^{-1}\right ) , \label{eq : avector } \end{aligned}\ ] ] and the function @xmath79 is given by : @xmath80 ^ 2 \hspace{0.5 cm } & \tau>1 \ , , \end{array } \right .
\label{eq : ftau } \end{aligned}\ ] ] with @xmath81 , and @xmath82 the mass of the particle running in the loop @xcite . in eq .
( [ eq : thm - h2gaga ] ) the first contribution is from the top quark , the second from the @xmath4 boson , the third from the singly - charged higgs boson , the fourth from the doubly - charged higgs boson , and the last two from the vector leptons .
we used the following expressions for the couplings of the higgs bosons with charged vector leptons : @xmath83 the couplings of @xmath84 to the vector bosons and fermions are as follows : @xmath85 with @xmath86 , and the scalar trilinear couplings are parametrized as follows @xmath87 \ , , \nonumber \\ & & \tilde{g}_{h h^+h^-}= \frac{m_w}{2 g m_{h^{\pm}}^2 }   \bigg\{\left[4v_\delta(\lambda_2 + \lambda_3 ) \cos^2{\beta_{\pm}}+2v_\delta\lambda_4\sin^2{\beta_{\pm}}- \sqrt{2}\lambda_5v_\phi   \cos{\beta_{\pm}}\sin{\beta_{\pm}}\right]\sin\alpha   \nonumber \\   & & + \left[4\lambda_1\,v_\phi   \sin{\beta_{\pm}}^2+{(2\lambda_{4}+\lambda_{5 } ) } v_\phi   \cos^2{\beta_{\pm}}+ ( 4\mu-\sqrt{2}\lambda_5v_\delta)\cos{\beta_{\pm}}\sin{\beta_{\pm}}\right]\cos\alpha\bigg\ }   \ , .
\label{eq : redgcallittlehhp}\end{aligned}\ ] ] since the new leptons do not affect the higgs production channels , the effect on the di - photon search channel at the lhc is expressed by the ratio @xmath88_{htm}}{[\sigma ( gg \to   \phi ) \times br(\phi \to \gamma \gamma ) ] _ { sm } } \nonumber \\ & = & \frac{[\sigma ( gg \to   h ) \times \gamma(h \to \gamma \gamma)]_{htm}}{[\sigma ( gg \to\phi ) \times \gamma(\phi \to \gamma \gamma)]_{sm } } \times \frac{[\gamma ( \phi)]_{sm}}{[\gamma(h)]_{htm } } , \end{aligned}\ ] ] where @xmath42 is the sm neutral higgs boson and where the ratio of the cross sections by gluon fusion is @xmath89 here @xmath90 is the mixing angle in the cp - even neutral sector : @xmath91 with @xmath92 in @xcite we investigated the higgs boson decay branching ratio into @xmath93 with respect to the sm considering the lightest higgs boson is the @xmath94 signal excess observed at the lep at 98 gev , while the heavier higgs boson is the boson observed at the lhc at 125 gev , in a higgs triplet model without vector leptons , and found that this is the only scenario which allows for enhancement of the @xmath0 branching fraction .
we thus set the values 125 gev and 98 gev for the @xmath84 and @xmath95 masses respectively , and adjust the parameters @xmath96 accordingly .
vector lepton masses and mixing parameters depend on @xmath47 and @xmath48 , the explicit mass parameters in the lagrangian ; and @xmath97 , the vector leptons yukawa parameters . in the limit of vanishing dirac mass terms @xmath47 and @xmath48 , the pre - factors @xmath98 in ( 11 ) go to one .
it then follows that there is destructive interference between the dominant @xmath4- boson contribution and the charged leptons loops @xcite . in this limit , despite possible enhancement from singly- and doubly - charged higgs bosons in the loop , we find a large suppression of the di - photon rate .
we present the plots for the relative signal strength of @xmath99 , defined in eq .
( [ eq : rgamgam ] ) as a function of @xmath100 ( or equivalently @xmath101 ) , for various values of doubly charged higgs bosons mass , in the left - side panel of fig .
[ fig : hgg ] , for @xmath102 . clearly , for this case ( no mixing ) the decay of the @xmath84 is suppressed significantly with respect to the value in sm over the whole region of the parameter space in @xmath103 .
-0.3 in    @xmath104{fig1leftpanel.pdf } & \hspace*{-1.5 cm }      \includegraphics[width=3.1in , height=3.1in]{fig1middlepanel.pdf } & \hspace*{-2.1 cm }      \includegraphics[width=3.1in , height=3.1in]{fig1rightpanel.pdf }          \end{array}$ ]    allowing mixing in the neutral higgs sector changes the relative contributions of the charged higgs to the di - photon decay .
we show decay rates for @xmath0 as a function of @xmath105 for different values of doubly - charged higgs boson mass considering @xmath106 gev ( and 200 gev ) in fig . [
fig : hgg ] middle ( and right ) panels respectively .
considerations for relative branching ratios are affected by the fact that the total width of higgs boson in the htm for @xmath107 is not the same as in the sm .
the relative widths factor is @xmath108_{htm}}{[\gamma(\phi)]_{sm}}=\frac{[\gamma ( h \to \sum\limits_f f { \bar f})+ \gamma ( h \to ww^ * ) + \gamma ( h \to zz^*)+ \gamma ( h \to \nu \nu)]_{htm}}{[\gamma ( \phi \to \sum\limits_f f { \bar f})+ \gamma ( \phi \to ww^ * ) + \gamma ( \phi \to zz^*)]_{sm } } .
\end{aligned}\ ] ] the plots in fig .
[ fig : hgg ] correspond to symmetry condition 1 in section [ sec : model ] , that is , @xmath109 .
however , if mixing with sm leptons is forbidden , but the vector leptons are still allowed to mix with each other , the pre - factors @xmath110 in eq . ( [ eq : thm - h2gaga ] ) are not one , and can modify the higgs di - photon decay . in the next plots we investigate the effect of non - zero mass parameters @xmath47 and @xmath48 , for fixed values of the yukawa couplings . in fig .
[ fig : hemlcont1 ] we present the contour plots of constant @xmath111 for @xmath112 in the plane of the explicit mass terms @xmath113 and @xmath48 , for various values of doubly - charged higgs bosons mass and @xmath114 .
the contours are labeled by the value of @xmath111 .
the vector lepton masses are restricted to values for which @xcite @xmath115 , where @xmath116 are given in eq . (
[ eq : eigenvalues ] ) .
the plots indicate that it is difficult to obtain any significant enhancement of the ratio @xmath111 for @xmath117 , and this does not depend on the chosen values for the doubly - charged higgs mass ; while for @xmath107 , enhancements are possible for various values of @xmath118 .
@xmath119{fig2upperleftpanel.pdf } & \hspace*{-0.4 cm }      \includegraphics[width=2.3in , height=2.2in]{fig2uppermiddlepanel.pdf } &    \hspace*{-0.4 cm }      \includegraphics[width=2.3 in , height=2.2in]{fig2upperrightpanel.pdf}\\ \hspace*{-0.9 cm }      \includegraphics[width=2.3in , height=2.2in]{fig2lowerleftpanel.pdf }      & \hspace*{-0.4 cm } \includegraphics[width=2.3in , height=2.2in]{fig2lowermiddlepanel.pdf } & \hspace*{-0.4 cm }      \includegraphics[width=2.3in , height=2.2in]{fig2lowerrightpanel.pdf }          \end{array}$ ]    in fig .
[ fig : hemlcont2 ] we investigate the dependence of @xmath111 on the yukawa couplings and vector lepton masses .
we show contour plots for fixed @xmath99 in a @xmath120 plane , with @xmath101 and @xmath121 , for various values of @xmath105 and doubly - charged higgs boson masses .
enhancements are possible here for all values of @xmath105 , but while for @xmath117 the decay @xmath0 is enhanced for large vector lepton masses and yukawa couplings , for @xmath107 we observe enhancements for light vector lepton masses and small yukawa couplings .
@xmath119{fig3upperleftpanel.pdf } & \hspace*{-0.4 cm }      \includegraphics[width=2.3in , height=2.2in]{fig3uppermiddlepanel.pdf } &    \hspace*{-0.4 cm }      \includegraphics[width=2.3 in , height=2.2in]{fig3upperrightpanel.pdf}\\ \hspace*{-0.9 cm }      \includegraphics[width=2.3in , height=2.2in]{fig3lowerleftpanel.pdf }      & \hspace*{-0.4 cm } \includegraphics[width=2.3in , height=2.2in]{fig3lowermiddlepanel.pdf } & \hspace*{-0.4 cm }      \includegraphics[width=2.3in , height=2.2in]{fig3lowerrightpanel.pdf }          \end{array}$ ]    if we wish to study the light vector leptons parameter space where @xmath0 is enhanced , @xmath107 is preferred .
the enhancement is affected by mixing in the vector lepton sector , the various values for doubly - charged higgs bosons mass and values of @xmath105 .
as the plots cover only a limited range of the parameter space , in the tables below we give the ranges for the values of the ratio @xmath122 for the various scenarios . in table
[ tab : y1/6_me_ml ] , we fix the value of the yukawa coupling to @xmath123 , allow the vector lepton masses to vary in the ( 100 - 500 ) gev range , and show the values for @xmath99 for different @xmath105 and doubly - charged higgs masses .
we note that the relative branching ratios are very sensitive to both doubly charged higgs mass values and values of @xmath105 .
enhancements in the branching ratio of @xmath0 are possible for light @xmath124 gev , and are much more pronounced at large @xmath105 . note that for @xmath117 , the result is independent of @xmath118 , in agreement with the results obtained in @xcite .
the reason is the following . in eq .
( [ eq : redgcallittlehhp ] ) , for @xmath117 , the coupling between neutral and doubly - charged higgs @xmath125 = \frac{m_w } { g m_{h^{\pm \pm}}^2}\bigg [    2\frac { m^2_{h^{\pm \pm}}}{v_\phi^2 } v_\phi \bigg]=2 \frac{m_w } { g v_\phi } , \end{aligned}\ ] ] where we used the expression for @xmath126 from @xcite , is independent of @xmath118 .    in table
[ tab : y1/6_h_e ] we allow , in addition to mass variations , variations in the yukawa coupling @xmath127 .
this means allowing both explicit ( dirac ) masses and additional contributions by electroweak symmetry breaking , @xmath128 .
the dependence on the yukawa coupling @xmath67 is much weaker than on @xmath105 or on @xmath118 .
however , one can see from the tables that modest enhancements of the ratio @xmath99 are possible for @xmath117 for large vector leptons yukawa couplings , unlike in the case of the triplet model without vector leptons .
this would then be a clear distinguishing feature : enhancements of the decay @xmath0 in the absence of mixing in the neutral sector .
the absence of mixing would manifest itself in observing tree - level decays ( @xmath129 and @xmath130 ) identical to those in the sm .
there seems to be a minimum value of @xmath99 for @xmath131 , where the contribution from the doubly - charged higgs bosons is important for small doubly - charged masses and counters the contribution from the vector leptons .
this is a suppression of the branching ratio for @xmath0 due to the fact that the vector lepton contribution interacts destructively with the dominant @xmath132 contribution . as a general feature ,
@xmath99 increases when we lower the doubly charged higgs mass and increase @xmath105 .
this rules out part of the parameter space .
for instance , for @xmath133 gev , the mixing can not be larger than @xmath134 , and for @xmath135 gev , mixings larger than @xmath136 are ruled out . if the value of @xmath137 is increased to @xmath138 gev , only modest enhancements are possible , and only for @xmath117 , for vector lepton explicit masses in the @xmath139 gev range and @xmath123 . increasing the vector leptons yukawa coupling increases the overall ratio @xmath99 .      in most models ,
the @xmath140 and @xmath141 partial decay widths are correlated or anti - correlated , though usually the enhancement / suppression in the @xmath142 channel is much smaller compared to that in the @xmath143 channel .
however , as in models with new loop - contributions to @xmath144 , sensitivity to both is expected , we study the correlation between the two here , in the presence of vector leptons .
investigation of the branching ratio of @xmath1 is also further justified by the recent results from cms and atlas @xcite , which indicate branching fractions consistent with the sm expectation at 1@xmath145 in the higgs boson @xmath84 mass region at 95% c.l . .
the decay width for @xmath1 is given by @xcite : @xmath146_{htm } & = & \frac{\alpha g_f^2   m_w^2 m_{h}^3 } { 64\pi^4 } \left ( 1-\frac{m_z^2}{m_h^2 } \right ) ^3   \bigg|\frac{1}{c_w }   \sum_{f } 2 n^f_c q_f ( i_3^f-2q_f s_w^2 ) g_{h ff }   a^h_{1/2 } ( \tau^h_f , \tau^z_f )   \nonumber \\   & + & \frac{(i_3^e-2q_e s_w^2)(2q_e)}{c_w } \bigg[\frac{\mu_{e_1}g_{hff}}{m_{e_1 } } a_{1/2}(\tau^h_{e_1 } , \tau^z_{e_1 } ) + \frac{\mu_{e_2}g_{hff}}{m_{e_2 } } a^h_{1/2}(\tau^h_{e_2 } , \tau^z_{e_2 } ) \bigg ]   \nonumber \\ & + & c_w g_{h ww } a^h_1 ( \tau^h_w , \tau^z_w)- 2s_w \tilde{g}_{h h^\pm\,h^\mp}g_{zh^\pm h^\mp } a^h_0(\tau^h_{h^{\pm } } , \tau^z_{h^\pm } ) \nonumber \\ & - &   4s_w \tilde{g}_{h h^{\pm\pm}h^{\mp\mp } } { g}_{z h^{\pm\pm}h^{\mp\mp } } a^h_0(\tau^h_{h^{\pm\pm } } , \tau^z_{h^{\pm \pm } } )         \bigg|^2 \ , , \end{aligned}\ ] ] where @xmath147 @xmath148 ( with @xmath149 ) , and the loop - factors are given by @xmath150i_1(\tau^h,\tau^z).\nonumber\end{aligned}\ ] ] the functions @xmath151 and @xmath152 are given by @xmath153\nonumber \\
& + & \frac{\tau^{h\,2}\tau^z}{\left(\tau^h-\tau^z\right)^2 } \left[g\left(\tau^{h\,-1}\right)-g\left(\tau^{z\,-1}\right)\right],\nonumber\\ i_2(\tau^h,\tau^z ) & = & -\frac{\tau^h\tau^z}{2(\tau^h-\tau^z)}\left[f\left(\tau^{h\,-1}\right)-f\left(\tau^{z\,-1}\right)\right],\end{aligned}\ ] ] where the function @xmath79 is defined in eq .
( [ eq : ftau ] ) , and the function @xmath154 is defined as @xmath155 , \hspace{0.5 cm } & ( \tau\geq 1 ) \ , .
\end{array } \right .
\label{eq : gtau } \end{aligned}\ ] ] in eq .
( [ eq : hzg ] ) we list , in order , the ordinary leptons , vector leptons , @xmath4 boson , singly - charged higgs , and doubly - charged higgs contribution . the scalar couplings @xmath156 and @xmath157 are given in eq .
( [ littleh1tt ] ) , and the scalar trilinear couplings @xmath158 and @xmath159 are given in eq .
( [ eq : redgcallittlehhp ] ) .
the remaining couplings in eq .
( [ eq : hzg ] ) are given by @xmath160 we proceed to perform a similar analysis as in sec .
[ subsec : gamgam ] .
we show first the variation of the branching ratio @xmath1 with the mass @xmath161 , for various values of the doubly - charged higgs masses , for the case of no mixing in the neutral sector , @xmath162 ( shown in fig .
[ fig : hgzme ] , left - hand panel ) , and as a function of the mixing angle @xmath105 for @xmath163 gev , and @xmath164 gev , in the middle and right panels of fig .
[ fig : hgzme ] , respectively .
we have chosen the same parameter values as in fig .
[ fig : hgg ] , for comparison .
it is clear that the branching ratio into @xmath165 is fairly independent of both @xmath166 and @xmath118 , and always just below the sm expectations .
note that the severe suppression seen in @xmath5 for @xmath117 ( fig .
[ fig : hgg ] , left side panel ) does not occur here , and the results of the left side of fig .
[ fig : hgzme ] are consistent with the data at lhc .
but the variation with the mixing angle @xmath90 is pronounced , and the branching ratio can reach almost twice its sm value for @xmath167 . however , correlated with our predictions from sec .
[ subsec : gamgam ] and lhc measurements for @xmath99 , the parameter space corresponding to an enhanced @xmath1 , for both @xmath168 gev and 200 gev , for doubly - charged higgs mass @xmath169 gev is ruled out .
for all other values considered , the value for @xmath170 is close to , or below the sm expectations .
this is general prediction of the model .
-0.3 in    @xmath171{fig4leftpanel.pdf } & \hspace*{-1.7 cm }      \includegraphics[width=3.1in , height=3.0in]{fig4middlepanel.pdf } &    \hspace*{-1.9 cm }      \includegraphics[width=3.1 in , height=3.0in]{fig4rightpanel.pdf }          \end{array}$ ]    for a large range of parameter space , the decay @xmath1 can be suppressed significantly with respect to the sm .
we plot decay rates for @xmath172 as a function of @xmath105 for different values of doubly charged higgs boson mass , considering @xmath173 gev and 200 gev in fig .
[ fig : hgzme ] , middle and right panel respectively .
again , considerations for relative branching ratios are affected by the fact that the total width of higgs boson in the htm is not the same as in the sm .
the widths are the same as those in the sm for @xmath84 for @xmath117 , while for @xmath107 we take into account the relative widths factors , eq.([eq : width ] ) .    in tables
[ tab : yz1/6_me_ml ] and [ tab : yz1/6_h_e ] we present the explicit ranges of @xmath170 for varying @xmath174 and for a range of @xmath67 parameters .
we choose a fixed value for @xmath123 in table [ tab : yz1/6_me_ml ] , as the preferred choice from other analyses @xcite , and for comparison with table [ tab : y1/6_me_ml ] .
comparison of tables [ tab : y1/6_me_ml ] and [ tab : yz1/6_me_ml ] shows that the decay @xmath1 is far more stable against variations in masses and values for @xmath105 than @xmath0 , making it a less sensitive indicator for the presence of vector lepton states .    in table
[ tab : yz1/6_h_e ] we also allow variations in the yukawa coupling @xmath127 . as before this amounts to allowing both explicit and contributions by electroweak symmetry breaking , @xmath128 , to vector lepton masses .
comparison of the tables [ tab : yz1/6_me_ml ] and [ tab : yz1/6_h_e ] indicates that the results are not very sensitive to variations in the yukawa coupling @xmath67 , or the vector lepton mass parameters @xmath174 .
however , the relative branching ratios are very sensitive to values of @xmath105 .
while the branching ratio into @xmath165 is almost always suppressed with respect to its sm value , there is a small region of the parameter space , light @xmath175 and @xmath176 where enhancement is possible ; but as discussed before , this region is ruled out by constraints from @xmath0 measurement , table [ tab : y1/6_h_e ] .
note that for @xmath117 the branching ratio is as before , independent of the mass of @xmath175 and about the same as in the sm .
the discovery of the doubly - charged higgs bosons would be one of the most striking signals of physics beyond sm , and a clear signature for the higgs triplet model .
the decay modes of @xmath175 depend on the value of the vev of the neutral triplet higgs component , @xmath177 . when @xmath178 mev , the dominant decay mode of @xmath175 is into lepton pairs .
if @xmath179 mev , the main decay modes of @xmath175 are into @xmath180 pairs are kinematically allowed . ] , and into @xmath181 , if kinematically allowed .
searches for @xmath175 were performed at large electron positron collider ( lep ) @xcite , the hadron electron ring accelerator ( hera ) @xcite and the tevatron @xcite .
the most up - to - date bounds have been more recently derived by atlas and cms collaborations at the lhc . assuming a drell - yan - like pair production , these collaborations have looked for long - lived doubly - charged states , and after analyzing 5 fb@xmath182 of lhc collisions at a center - of - mass energy @xmath183 tev and 18.8 fb@xmath182 of collisions at @xmath184 = 8 tev , they constrained the masses to lie above 685 gev @xcite .
the assumption is that the doubly - charged higgs bosons decay 100% into a pair of leptons with the same electric charge through majorana - type interactions @xcite , thus neglecting possible decays into @xmath4-boson pairs , shown to alter the pattern of @xmath175 branching fractions @xcite . in this work , we allow both decays into @xmath185 bosons , and also include the effects of decays into vector leptons , which , if light enough , would modify the decays of the doubly - charged higgs bosons further .
we take @xmath186 gev throughout our considerations is small enough to satisfy electroweak precision conditions , but large enough to allow decay into gauge boson and charged higgs @xcite . ] .
the main production mode for @xmath175 is the pair production @xmath187 and the associated production @xmath188 .
the production cross sections for both the vector boson fusion @xmath189 , and for weak boson associated production @xmath190 , are proportional to @xmath191 and much less significant for @xmath192 .
at hadron colliders the partonic cross section for the leading order ( lo ) production cross section for doubly charged higgs boson pair is @xmath193,\ ] ] where we defined @xmath194 with @xmath195 the third component of the isospin for particle @xmath196 , @xmath197 the square of the partonic center of mass energy , @xmath198 , and @xmath90 the qed coupling constant evaluated at scale @xmath199 .
the hadronic cross section is obtained by convolution with the partonic density functions of the proton @xmath200 where @xmath201 is the parton luminosity and @xmath202 ( @xmath203 is the total energy squared at the lhc ) .
the cross section for pair - production , including nlo corrections , has been evaluated in @xcite .
depending on mass parameters in the model , the doubly - charged higgs boson can decay into lepton pairs , including vector leptons , @xmath132 pairs , or @xmath204 states . in the higgs triplet model ,
the decay rate for @xmath175 into leptons is @xmath205 ^ 2 , \end{aligned}\ ] ] where @xmath82 is the mass of the @xmath206th lepton ( @xmath207 or @xmath208 ) and @xmath209 for @xmath210 , ( @xmath211 ) .
similarly the decay rate of @xmath175 into fourth generation vector leptons is , if kinematically allowed @xmath212 \frac{m_{h^{\pm \pm}}}{4 \pi}\left ( 1-\frac{m_{e_i}^2}{m_{h^{\pm \pm}}^2}-\frac{m_{e_j}^2}{m_{h^{\pm \pm}}^2 } \right ) \left [ \lambda \left(\frac{m_{e_i}^2}{m_{h^{\pm \pm}}^2 } , \frac{m_{e_j}^2}{m_{h^{\pm \pm}}^2 } \right ) \right]^2 , \nonumber \\\end{aligned}\ ] ] with @xmath213 the mass eigenvalue from eq .
( [ eq : eigenvalues ] ) .
in addition , we use the decay rates of @xmath175 into @xmath214 and @xmath215 : @xmath216^{3/2 } , \end{aligned}\ ] ] where @xmath217 is the mixing angle in the singly - charged higgs sector and @xmath218    we investigate the branching ratios of @xmath175 in two distinct parameter regions :    * condition 1 : when @xmath219 is not kinematically allowed .
then @xmath175 decays into leptons and @xmath132 pairs only : @xmath220 with the total width for condition 1 @xmath221 * condition 2 : when @xmath219 is kinematically allowed . then @xmath175 is able to decay into charged higgs and gauge bosons as well : @xmath222 with the total decay width for condition 2 @xmath223    we investigate the decay patterns of @xmath175 and present plots of the production cross section times the branching fractions under various conditions . to cover a wide range of parameter space , we distinguish two cases for each condition , depending on the vector lepton masses .
we set the yukawa coupling of the vector leptons with the doublet higgs bosons to be @xmath224 for both cases .
* case a corresponds to very light vector leptons : @xmath225 gev . for this case
we obtain for the eigenvalues , @xmath226 gev , @xmath227 gev , the latter of which is close to the allowed minimum .
* case b corresponds to intermediate mass vector leptons @xmath228 gev , @xmath229 gev .
for this case we obtain for the eigenvalues , @xmath230 gev , @xmath231 gev .    in fig . [
fig : condition1 ] we show the corresponding graphs for condition 1 ( when the decay @xmath219 is not kinematically allowed ) , case a on the top row and case b on the bottom row .
we plot @xmath232 with @xmath233 as specified in the attached panels , as functions of the doubly - charged masses . on the left side of the figure the yukawa couplings of vector lepton with the triplet higgs boson are @xmath234 and on the right side of the figure , @xmath235 .
we set @xmath236 throughout .    for small vector lepton coupling @xmath237 ,
the decay into @xmath238 is dominant at low @xmath175 masses ( this is because we assumed the couplings @xmath32 with ordinary leptons to be all equal ; should we have chosen them diagonal , the branching ratio into @xmath239 would dominate ) . at high @xmath240 masses the branching ratio into @xmath214 dominates and can reach 40% , branching ratios into ordinary leptons reaching at most 40% , and those into vector leptons remaining at or below 20% level .
if the triplet yukawa coupling to vector leptons is allowed to increase to @xmath241 , the decay into the two lightest vector leptons @xmath242 becomes dominant , and can reach 80 - 90% when kinematically allowed ( in the @xmath243 gev region for @xmath227 gev ) and overwhelms the other decay modes , which are now below 5% .
the only difference between case a and case b in this figure are threshold effects .
for case b , @xmath244 gev for decay into pairs of @xmath242 states , as @xmath231 gev ; otherwise the branching ratios are the same .
-0.6 in    @xmath245{fig5topleftcondition1case1yukawa001.pdf } & \hspace*{-0.8 cm }      \includegraphics[width=3.5in , height=3.5in]{fig5toprightcondition1case1yukawa01.pdf}\\      \hspace*{-1.2 cm }      \includegraphics[width=3.5 in , height=3.5in]{fig5bottomleftcondition1case2yukawa001.pdf } & \hspace*{-0.8 cm }      \includegraphics[width=3.5in , height=3.5in]{fig5bottomrightcondition1case2yukawa01.pdf }          \end{array}$ ]    in fig .
[ fig : condition2 ] we plot the same quantities for condition 2 ( when @xmath219 is kinematically allowed ) , also case a on the top row and case b on the bottom row , @xmath235 on the left hand side and @xmath241 on the right hand side .
the decay pattern is very different here , and it is dominated by @xmath219 . for small vector lepton coupling @xmath246
the branching ratios for decays into vector and ordinary lepton pairs , and @xmath132 pairs are very small , and reach at most 1% . increasing the vector lepton coupling to @xmath241 ,
the decay into @xmath204 still dominates throughout the parameter space where it is kinematically allowed and can reach a branching fraction of over 90% , while the decay into vector leptons can have branching ratios of up to 25% .
again , the decay rates into @xmath132 boson pairs and ordinary leptons are below 1% , and the only difference between case a and case b are , as in fig .
[ fig : condition1 ] , threshold effects .
-0.6 in    @xmath247{fig6topleftcondition2case1yukawa001.pdf } & \hspace*{-0.6 cm } \includegraphics[width=3.5in , height=3.5in]{fig6toprightcondition2case1yukawa01.pdf}\\ \hspace*{-1.2 cm }      \includegraphics[width=3.5 in , height=3.5in]{fig6bottomleftcondition2case2yukawa001.pdf } & \hspace*{-0.6 cm }      \includegraphics[width=3.5in , height=3.5in]{fig6bottomrightcondition2case2yukawa01.pdf }          \end{array}$ ]
despite no new signals of physics beyond the sm at the lhc , the sm can not be the complete theory of particle interactions .
an extension of the sm by additional vector - like leptons is not ruled out experimentally , and has been shown to provide a dark matter candidate . in models beyond the sm , the vector leptons can alter not only the phenomenology of the higgs , but also of other additional particle representations predicted by the models .
we provide an example within the higgs triplet model , where previously we showed that , in the absence of triplet - doublet higgs mixing in the neutral sector ( @xmath117 ) , there is no enhancement of the rate of decay of @xmath0 in this model with respect to the sm expectation .
introducing vector leptons does not affect any of the tree - level decays of the neutral higgs boson observed at the lhc , nor its production decay .
however , loop decays into electroweak particles , such as @xmath1 and @xmath0 would be affected .
we show that , for the no - mixing scenario ( @xmath117 ) the decays rates into @xmath93 and @xmath165 do not depend on the doubly charged higgs mass , and thus without the additional vector leptons , these decays would be unchanged from the sm expectations . with vector leptons , modest enhancements or suppressions
are possible , more so for @xmath0 , where for large yukawa couplings , the rate of decays could even double . under the same circumstances ,
the decay width for @xmath1 remains practically unchanged from its sm value .
the model thus presents a mechanism for enhancing one loop - decays and not the other , which seems to be consistent with the lhc data ( so far ) .    if @xmath107 , the effect of the doubly - charged higgs bosons is felt for both @xmath248 and @xmath249 , most spectacularly so for very light @xmath169 gev , which is ruled out for @xmath250 .
parameter space regions where light doubly - charged higgs masses ( 200 - 250 gev ) , _ and _ significant mixing in the neutral sector coexist , are disfavored . in general
, there are many parameter combinations for which the decay @xmath248 is enhanced , but few for an enhanced @xmath251 , and these regions are ruled out by the branching ratio for @xmath0 . however , if the decay @xmath248 is ( modestly ) enhanced , while @xmath1 is the same as the sm to 1@xmath145 , small mixing angles and light doubly charged higgs bosons @xmath252 gev are preferred .
the fact that there are no regions of the parameter space consistent with present measurements of the branching ratio for @xmath0 and an enhanced rate for @xmath1 is a feature of this model , valid over the whole explored range of the parameter space .
other than that , there are no definite correlations or anti - correlations between there two loop - dominated decays .    the intermediate mass doubly - charged higgs boson can decay into light vector leptons , which would alter its decay profile significantly .
we explored this possibility and found that , if the singly charged higgs mass is such that the decay @xmath219 is not kinematically accessible , dominant branching ratios into vector leptons , if kinematically accessible , are expected for triplet yukawa couplings @xmath253 , whereas for small @xmath254 branching ratios into ordinary leptons , vector leptons and @xmath132 pairs are comparable for @xmath255 gev . if and where the decay @xmath219 is kinematically accessible , its corresponding branching ratio is the largest , while the branching fraction into vector leptons could reach 20 - 25% for @xmath253 . under both these circumstances ,
the decay patterns of the doubly charged higgs bosons are significantly altered , raising the hope that they can be found at masses around 300 - 500 gev .
the work of s. b. and m.f . is supported in part by nserc under grant number sap105354 .
g.  aad _ et al . _ [ atlas collaboration ] , phys .
b * 716 * ( 2012 ) 1 ; s.  chatrchyan _ et al .
_ [ cms collaboration ] , phys .
b * 716 * ( 2012 ) 30 .
t.  aaltonen _ et al . _
[ cdf collaboration ] , phys .
lett .   * 104 * , 091801 ( 2010 ) ; a.  lister [ cdf collaboration ] , arxiv:0810.3349 [ hep - ex ] .
k.  agashe , g.  perez and a.  soni , phys .
d * 71 * , 016002 ( 2005 ) ; k.  agashe , g.  perez and a.  soni , phys .
d * 75 * , 015002 ( 2007 ) ; a.  carmona and j.  santiago , jhep * 1201 * , 100 ( 2012 ) ; g.  -y .
huang , k.  kong and s.  c.  park , jhep * 1206 * , 099 ( 2012 ) ; c.  biggio , f.  feruglio , i.  masina and m.  perez - victoria , nucl .
b * 677 * , 451 ( 2004 ) ; d.  e.  kaplan and t.  m.  p.  tait , jhep * 0006 * ( 2000 ) 020 ; h.  -c .
cheng , phys .
d * 60 * , 075015 ( 1999 ) .
t.  moroi and y.  okada , mod .
phys .  lett .  a * 7 * ( 1992 ) 187 ; t.  moroi and y.  okada , phys .
b * 295 * ( 1992 ) 73 ; k.  s.  babu , i.  gogoladze , m.  u.  rehman and q.  shafi , phys .
d * 78 * ( 2008 ) 055017 ; s.  p.
martin , phys .
d * 81 * ( 2010 ) 035004 ; p.
w.  graham , a.  ismail , s.  rajendran and p.  saraswat , phys .
d * 81 * ( 2010 ) 055016 ; s.  p.
martin , phys .
d * 82 * ( 2010 ) 055019 ; j.  kang , p.  langacker and b.  d.  nelson , phys .
d * 77 * ( 2008 ) 035003 .
b.  a.  dobrescu and c.  t.  hill , phys .
* 81 * ( 1998 ) 2634 ; r.  s.  chivukula , b.  a.  dobrescu , h.  georgi and c.  t.  hill , phys .  rev .
d * 59 * ( 1999 ) 075003 ; h.  collins , a.  k.  grant and h.  georgi , phys .  rev .
d * 61 * ( 2000 ) 055002 ; h.  -j .
he , c.  t.  hill and t.  m.  p.  tait , phys .
d * 65 * ( 2002 ) 055006 ; c.  t.  hill and e.  h.  simmons , phys .
* 381 * ( 2003 ) 235 [ erratum - ibid .
* 390 * ( 2004 ) 553 ] ; r.  contino , l.  da rold and a.  pomarol , phys .
d * 75 * ( 2007 ) 055014 ; phys .
d * 79 * ( 2009 ) 075003 ; k.  kong , m.  mccaskey and g.  w.  wilson , jhep * 1204 * ( 2012 ) 079 ; a.  carmona , m.  chala and j.  santiago , jhep * 1207 * ( 2012 ) 049 ; m.  gillioz , r.  grober , c.  grojean ,
m.  muhlleitner and e.  salvioni , jhep * 1210 * ( 2012 ) 004 . n.  arkani - hamed , a.  g.  cohen , e.  katz and a.  e.  nelson , jhep * 0207 * ( 2002 ) 034 ; t.  han , h.  e.  logan , b.  mcelrath and l.  -t .
wang , phys .
d * 67 * ( 2003 ) 095004 ; m.  perelstein , m.  e.  peskin and a.  pierce , phys .
d * 69 * ( 2004 ) 075002 ; m.  schmaltz and d.  tucker - smith , ann .
nucl .  part .
* 55 * ( 2005 ) 229 ; m.  s.  carena , j.  hubisz , m.  perelstein and p.  verdier , phys .  rev .
d * 75 * ( 2007 ) 091701 : s.  matsumoto , t.  moroi and k.  tobe , phys .
d * 78 * ( 2008 ) 055018 .
a.  davidson and k.  c.  wali , phys .
* 60 * ( 1988 ) 1813 ; k.  s.  babu and r.  n.  mohapatra , phys .
d * 41 * ( 1990 ) 1286 ; b.  grinstein , m.  redi and g.  villadoro , jhep * 1011 * ( 2010 ) 067 ; d.  guadagnoli , r.  n.  mohapatra and i.  sung , jhep * 1104 * ( 2011 ) 093 .
a.  azatov , m.  toharia and l.  zhu , phys .
d * 80 * , 035016 ( 2009 ) .
atlas collaboration atlas - conf-2013 - 019 ; s. chatrchyan et al .
[ cms collaboration ] , phys .
b * 718 * , 348 ( 2012 ) .
f.  del aguila , j.  de blas and m.  perez - victoria , phys .
d * 78 * , 013010 ( 2008 ) .
a.  joglekar , p.  schwaller and c.  e.  m.  wagner , jhep * 1212 * , 064 ( 2012 ) .
j.  kearney , a.  pierce and n.  weiner , phys .
d * 86 * , 113005 ( 2012 ) ; h.  davoudiasl , h.  -s .
lee and w.  j.  marciano , phys .
d * 86 * , 095009 ( 2012 ) ; k.  j.  bae , t.  h.  jung and h.  d.  kim , phys .  rev .
d * 87 * , 015014 ( 2013 ) ; m.  b.  voloshin , phys .
d * 86 * , 093016 ( 2012 ) ; d.  mckeen , m.  pospelov and a.  ritz , phys .
d * 86 * , 113004 ( 2012 ) ; h.  m.  lee , m.  park and w.  -i .  park ,
jhep * 1212 * , 037 ( 2012 ) ; c.  arina , r.  n.  mohapatra and n.  sahu , arxiv:1211.0435 [ hep - ph ] ; b.  batell , s.  jung and h.  m.  lee , jhep * 1301 * , 135 ( 2013 ) ; h.  davoudiasl , i.  lewis and e.  ponton , arxiv:1211.3449 [ hep - ph ] .
; j.  fan and m.  reece , arxiv:1301.2597 [ hep - ph ] ; a.  carmona and f.  goertz , arxiv:1301.5856 [ hep - ph ] ; c.  cheung , s.  d.  mcdermott and k.  m.  zurek , arxiv:1302.0314 [ hep - ph ] ; w.  -z .
feng and p.  nath , arxiv:1303.0289 [ hep - ph ] ; c.  englert and m.  mccullough , arxiv:1303.1526 [ hep - ph ]
. k.  ishiwata and m.  b.  wise , phys .
d * 84 * , 055025 ( 2011 ) .
s.  k.  garg and c.  s.  kim , arxiv:1305.4712 [ hep - ph ] .
a.  joglekar , p.  schwaller and c.  e.  m.  wagner , arxiv:1303.2969 [ hep - ph ] . g.  f.  giudice and o.  lebedev , phys .
b * 665 * , 79 ( 2008 ) ; e.  del nobile , r.  franceschini , d.  pappadopulo and a.  strumia , nucl .
b * 826 * , 217 ( 2010 ) ; s.  p.  martin , phys .  rev .
d * 81 * , 035004 ( 2010 ) .    f.  arbabifar , s.  bahrami and m.  frank , phys .
d * 87 * , 015020 ( 2013 ) .
[ cms collaboration ] , cms - pas - hig-13 - 001 and cms - pas - hig-13 - 005 ; g. aad et al .
[ atlas collaboration ] , [ arxiv:1307.1427 [ hep - ex ] ] and atlas - conf 2013 - 029 .
a.  g. akeroyd , mayumi aoki , and hiroaki sugiyama .
d * 77 * , 075010 ( 2008 ) ; a.  g. akeroyd , mayumi aoki , and hiroaki sugiyama .
d * 79 * , 113010 ( 2009 ) ; takeshi fukuyama , hiroaki sugiyama , and koji tsumura .
jhep * 03 * , 044 ( 2010 ) ; s.  t. petcov , h.  sugiyama , and y.  takanishi .
d * 80 * , 015005 ( 2009 ) ; takeshi fukuyama , hiroaki sugiyama , and koji tsumura .
d * 82 * , 036004 ( 2010 ) ; a.g .
akeroyd , cheng - wei chiang , and naveen gaur , jhep * 1011 * , 005 ( 2010 ) ; a.  arhrib , r.  benbrik , m.  chabab , g.  moultaka and l.  rahili , jhep * 1204 * , 136 ( 2012 ) ; a.  arhrib , r.  benbrik , m.  chabab , g.  moultaka and l.  rahili , [ arxiv:1202.6621 [ hep - ph ] ] ; a.  arhrib , r.  benbrik , m.  chabab , g.  moultaka , m.  c.  peyranere , l.  rahili and j.  ramadan , phys .
d * 84 * , 095005 ( 2011 ) ; e.  j.  chun , h.  m.  lee and p.  sharma , jhep * 1211 * , 106 ( 2012 ) ; a. melfo , m. nemevek , f. nesti , g. senjanovi and y. zhang , phys . rev .
d 85 , 055018 ( 2012 ) ; m.  aoki , s.  kanemura and k.  yagyu , phys .
d * 85 * , 055007 ( 2012 ) ; s.  kanemura and k.  yagyu , phys .  rev .
d * 85 * , 115009 ( 2012 ) ; m.  aoki , s.  kanemura , m.  kikuchi and k.  yagyu , [ arxiv:1211.6029 [ hep - ph ] ] .
k.  ishiwata and m.  b.  wise , arxiv:1307.1112 [ hep - ph ] .
r.  dermisek and a.  raval , arxiv:1305.3522 [ hep - ph ]
. s.  chatrchyan _ et al . _
[ cms collaboration ] , arxiv:1307.5515 [ hep - ex ] ; [ atlas collaboration ] , atlas - conf-2013 - 009 .    c.  -s .
chen , c.  -q .
geng , d.  huang and l.  -h .
tsai , phys .
d87 , * 075019 * ( 2013 ) ; c.  -s .
chen , c.  -q .
geng , d.  huang and l.  -h .
tsai , phys .
b * 723 * , 156 ( 2013 ) ; p.  s.  bhupal dev , d.  k.  ghosh , n.  okada and i.  saha , jhep * 1303 * , 150 ( 2013 ) [ erratum - ibid .
* 1305 * , 049 ( 2013 ) ] .
g.  abbiendi _ et al . _
[ opal collaboration ] , phys .
b * 526 * ( 2002 ) 221 ; g.  abbiendi _ et al .
_ [ opal collaboration ] , phys .
b * 577 * ( 2003 ) 93 ; p.  achard _ et al .
_ [ l3 collaboration ] , phys .
b * 576 * ( 2003 ) 18 ; j.  abdallah _ et al .
_ [ delphi collaboration ] , phys .
b * 552 * ( 2003 ) 127 .
j.  sztuk - dambietz [ h1 and zeus collaboration ] , j.  phys .
* 110 * ( 2008 ) 072042 ; a.  aktas _ et al . _
[ h1 collaboration ] , phys .
b * 638 * ( 2006 ) 432 ; v.  m.  abazov _ et al . _
[ d0 collaboration ] , phys .
* 101 * ( 2008 ) 071803 ; v.  m.  abazov _ et al .
_ [ d0 collaboration ] , phys .  rev .
* 108 * ( 2012 ) 021801 ; d.  e.  acosta _ et al . _
[ cdf collaboration ] , phys .
* 93 * ( 2004 ) 221802 ; t.  aaltonen _ et al . _
[ the cdf collaboration ] , phys .
* 101 * ( 2008 ) 121801 . s.  chatrchyan
_ [ cms collaboration ] , arxiv:1305.0491 [ hep - ex ] ; g.  aad _ et al . _ [ atlas collaboration ] , phys .
b 722 ( 2013 ) 305 ; s.  chatrchyan _ et al .
_ [ cms collaboration ] , eur .
j.  c * 72 * ( 2012 ) 2189 ; g.  aad _ et al .
_ [ atlas collaboration ] , eur .
j.  c * 72 * ( 2012 ) 2244 .
c.  -w .
chiang , t.  nomura and k.  tsumura , phys .
d * 85 * , 095023 ( 2012 ) ; s.  kanemura , k.  yagyu and h.  yokoya , arxiv:1305.2383 [ hep - ph ] . m.  muhlleitner and m.  spira , phys .  rev .
d * 68 * , 117701 ( 2003 ) .


edge localised modes ( elms ) are bursty , quasi - periodic expulsions of energy and particles from the plasma edge in the high - confinement state of toroidal magnetised plasmas @xcite .
elms originate from the steep pressure gradient pedestal in the outer closed field line region of a tokamak .
a substantial release of magnetic and potential energy , caused by a rapid rise of the rate of momentum transfer through the flux surfaces , leads to enhanced transport and loss of heat and particles into the open scrape - off layer ( sol ) field line region .
the filamented peak heat fluxes during elms on the bounding divertor plates present serious restrictions on the performance of future fusion experiments like iter @xcite and necessitate the development of suitable control techniques @xcite .
elms generically have been observed in all divertor tokamaks since the initial discovery of the high - confinement `` h - mode '' state in the early 1980s @xcite , while characteristics like frequency and intensity are found to vary widely depending on experimental conditions @xcite .
stellarator experiments also are reported to show similar bursty quasi - periodic edge activity in the presence of edge transport barriers @xcite .
indication that equilibrium drifts are part of the process has been found by experimental diagnosis of an asymmetry of the energy and particle fluxes on each of the divertor plates and evidence that this is sensitive to the directions of the toroidal current and magnetic field @xcite .
phenomenologically , edge localised mode events in tokamaks have been compared to solar flare eruptions @xcite , and , in the cataclysmic variability of observed radiation emanating from the plasma , they may seem to bear also some remote resemblance to the outbursts of pressure - driven dwarf novae @xcite , although the specific instability mechanisms behind these phenomena are clearly of a substantially different nature .
on the other hand , the nature of elms has also not yet been completely clarified , in the sense that there is no first - principles based predictive theory available .
a number of characteristically different variations of elms have been observed in experiments ( e.g. type i , ii or iii elms ) which may actually be caused by different instability mechanisms @xcite . ``
type iii '' elm observations show a remarkable similarity to global bursts found in computations of drift - wave turbulence , where the turbulence generated flows and fluxes are closely tied to the self - consistent equilibrium evolution in the plasma edge @xcite .
`` type i '' elms in the experiment specifically show stronger magnetic activity during the burst and are usually associated with the onset of a magneto - hydrodynamic ( mhd ) ideal ballooning mode ( ibm ) instability when a threshold pressure gradient is reached in the steep edge pedestal of an h - mode plasma @xcite .
the possible role of an additional current driven instability ( `` peeling - ballooning mode '' ) for `` type i '' elms has also received recent interest @xcite .
the ideal ballooning instability with nonlinear phases or aftermath has provided a paradigm for recent and current study of the elm phenomenon theoretically . explosive mhd instability and associated critical phenomena
have been advanced analytically @xcite .
early braginskii fluid simulations were given in support of this @xcite .
large , nonlinear mhd codes have been studying this in more detail , resolving several numerical problems and paying specific attention to the details of the magnetic geometry , with mixed results on the phenomenology of the instability beyond its linear phase @xcite .
whether or not the ideal ballooning instability or its peeling / ballooning variant is actually responsible for elms occuring in experiments ( here noting the absence of a demonstrably well - resolved l - to - h mode transition in computations from first principles ) , the nonlinear phase of this `` ideal mhd blowout '' or `` ibm blowout '' phenomenon is of physical interest .
the ability of the instability to destroy the edge pedestal layer is undisputed , and the term `` blowout '' is germane .
typical physical parameters lead to a situation wherein the ideal mhd interchange growth rate is comparable to the parallel alfvn transit time , one and the same with the basic ideal ballooning criterion .
this is fast with respect to mhd but not to microturbulence .
provided they can simultaneously treat global mhd ( i.e. , the spectrum covers both global and ion gyroradius scales ) , edge - turbulence computations can also address the phenomenon and are even in a position to treat parts of it which fall outside the paradigm of the mhd model .
here we present nonlinear gyrofluid computations of ibm blowout events localised to the edge / sol region of a tokamak plasma , resolved to below the ion gyroradius @xmath0 . as in the large mhd models
the self - consistent evolution of the equilibrium is not only included but is also an integral part of the dependent variables .
an electromagnetic 6-moment gyro - fluid model ( `` gemr '' ) for both electrons and ions is used , including an energy - conserving treatment of finite ion larmor radius effects and higher - moment terms which occur in the toroidal drift @xcite .
the energetic consistency is a key feature since the reaction of parts of the system such as flows and currents which have low energy content can play a central role in the indirect nature of nonlinear dynamics by serving as transfer channels .
although the gyrofluid model has a different formulation of nonlinear polarisation ( the generalisation of `` vorticity '' to a two - fluid setting ) than a fluid or mhd model , full correspondence in the regime of validity of the latter models has been shown @xcite . in this sense
the model is a superset of global reduced mhd and electromagnetic microturbulence .
the global dynamics is in the reduced mhd regime @xcite due to low absolute edge beta values ( @xmath1 ) but the small scales depend on treating the ion gyroradius to arbitrary order .
the model does not assume instabilities to occur at any particular scale , but when they occur the resulting spread of the spectrum is found to reach @xmath0 within a few eddy turnover times of the onset of nonlinearity in every case computed .
the work reported herein addressed the following points .
the relationship of the mhd dynamics to the ideal threshold is much discussed but in the gyrofluid model several mechanisms of microinstability are also present .
the threshold issue is obscured by the existence of ion temperature gradient ( itg ) drift instabilities at all beta values , with a smooth transition between them .
nevertheless , varying character in the saturation phases is observed at different values .
we also address the issue of possible `` rho - star '' dependence ( here the local rho - star is given by @xmath2 where @xmath3 is @xmath0 evaluated at @xmath4 and @xmath5 is the scale length of the @xmath6 profile , both evaluated at the mid - pedestal location in minor radius ) which would indicate a role for diamagnetic drifts .
the gradient is sufficiently steep that @xmath7 so that none of these effects can be ordered small .
we examine the necessary resolution to obtain converged cases , which gives the spectral range actively involved in nonlinear saturation .
the main instability at experimentally relevant beta is the ideal ballooning instability but the main saturation process is energy exchange with broadband electromagnetic drift wave / itg turbulence which the instability itself generates .
due mostly to this , the results turn out to be outside of both the magneto - hydrodynamic ( mhd ) and collisional braginskii paradigms , on which most previous approaches have been based . on the other hand ,
the gemr model still uses delta - f equations , so the actual profile phenomenology in the sol region is not well represented , and as in other work with gemr we concentrate on the properties of the nonlinear dynamics @xcite .
following sections give the details of the gemr model , a discussion of the physical difference between itg and mhd dynamics and their relative roles , the methods used to pre - set the axisymmetric equilibrium state so that the computations are done in the absence of axisymmetric oscillations whose decay times can be slower than the ibm rise time , the computations of the blowout phase itself , and the possibility and limitations of quantitative comparison to experiment .
the gyrofluid model used herein is given in ref .
it is based on the original derivations given in refs.@xcite , with corrections given as motivated by free energy conservation .
an alternative derivation using the conservation laws of the underlying gyrokinetic model is given in ref.@xcite .
the equations are normalised to spatial scale @xmath8 and time scale @xmath9 where @xmath8 is the minor radius and @xmath10 is the reference sound speed given by @xmath11 , that is , in terms of the electron temperature ( in energy units ) and ion mass .
the main parameters are the drift parameter , electron dynamical beta , and electron collisionality , respectively given by @xmath12 where @xmath13 and @xmath14 are the speed of light and fundamental charge , @xmath15 and @xmath16 are the pressure and braginskii collision time @xcite of the electrons , and @xmath17 is the drift scale . note that if @xmath4 then @xmath18 . the reference for @xmath19 is the deuterium mass @xmath20
. for the electrons the physical value of the mass @xmath21 is always used . with global geometry
there is no single magnetic field line connection length , but one with a profile with @xmath22 a function of the minor radius coordinate , where @xmath23 is a constant giving the reference value of the major radius .
each species @xmath24 has its own set of gyrofluid moment variable equations , one each for density @xmath25 , parallel velocity @xmath26 , parallel and perpendicular temperature @xmath27 and @xmath28 , and parallel and perpendicular energy components of the parallel heat flux @xmath29 and @xmath30 .
they are coupled to the electrostatic and parallel magnetic potentials @xmath31 and @xmath32 through self consistent field equations .
the tilde denotes a dependent variable to distinguish from constant parameters .
the correspondence of these equations to the braginskii fluid model is given in ref .
each species is characterised by a background charge density , temperature / charge ratio , and mass / charge ratio , given by the normalised parameters @xmath33 for electrons @xmath34 and @xmath35 . for the main ions
@xmath36 while @xmath37 is kept as a parameter . for the trace ions @xmath38 is always taken .
for hot trace ions @xmath39 is large while @xmath40 is still moderate .
the species gyroradius @xmath41 is given by @xmath42 , which is always small for electrons , close to unity for main ions , and can be moderate or large for trace ions .
for this study a simplified geometry is taken with @xmath43 except in curvature terms , so magnetic pumping and shaping effects @xcite are neglected . since dynamics at the electron gyroradius scale are neglected the gyroaveraging effects on @xmath32 are not treated .
the moment equations are given by @xmath44 @xmath45 @xmath46 @xmath47 @xmath48 @xmath49 the linearised pressures are given by @xmath50 collisional dissipation is controlled by a collision parameter @xmath51 analogous to @xmath52 and a set of numerical constants for each species , @xmath53 which are the thermal force , thermal conductivity , and viscosity coefficients involved in parallel dissipation , with values given by the braginskii model @xcite , though that regime is never reached in core turbulence . for electrons
these coefficients are @xmath54 and @xmath55 and @xmath56 , and for ions they are @xmath57 and @xmath58 and @xmath59 , respectively .
the resistive dissipation is given by @xmath60 with parallel current given by @xmath61 noting that it is the electron heat fluxes that enter @xmath62 and as with @xmath63 the @xmath62 term enters the parallel motion in the same way for every species .
the anisotropy dissipation is given by @xmath64 the thermal conduction components are given by @xmath65 @xmath66 @xmath67 the landau damping effects are modeled by @xmath68 with thermal velocity @xmath69 given by @xmath70 .
the gyroaveraging is done through pad approximants to operators which would act in fourier space , @xmath71 and then the gyroaveraged potentials are @xmath72 noting that they are species dependent . the main exb advection and
the flr correction pieces are given by @xmath73 - \nu_\perp\dpp^4 f + \nu_\parallel\ddpl f\ ] ] @xmath74\ ] ] where the @xmath75 $ ] symbols denote the nonlinear bracket defined below , and @xmath76 denote the artificial dissipation coefficients .
the parallel derivative is given by @xmath77\ ] ] the curvature operator is given by @xmath78\ ] ]      the species are coupled through two equations which give the self consistent response of the field potentials to the evolution of the gyrofluid moment variables . the electrostatic potential is governed by quasineutrality , which sets the total space charge density to zero .
the density for each species is given by a part due to the gyrocenters and another due to polarisation , which is what sets @xmath79 .
this polarisation equation is given by @xmath80 where @xmath81 gives the gyroscreening of the main ions .
in this version of the model the electron and trace ion contributions to gyroscreening are neglected , due to small @xmath82 and zero @xmath83 , respectively . on the right hand side
the electrons and main ions enter with oppositely signed @xmath83 values and again here the trace ions do not enter .
this establishes them as a trace species .    in a similar way the trace ions
are left out of the induction equation due to the zero @xmath83 , so that @xmath84 determines @xmath32 noting the way that the normalisation scales it with respect to @xmath85 .
more detail , including energy conservation and the relationship of these equations to it , is given in refs .
@xcite .      in these expressions
the operators @xmath86 and @xmath87 and the nonlinear brackets are determined by the representation of the geometry that is used .    the magnetic geometry is a simplified representation of an axisymmetric magnetic field using field aligned hamada coordinates .
the field representation is @xmath88 where @xmath89 is the physical toroidal angle about the symmetry axis , @xmath90 is the magnetic flux function , and @xmath91 is given by the constraint of mhd equilibrium .
without loss of generality one may define magnetic flux coordinates @xmath92 with the following properties , @xmath93 where @xmath94 is the volume enclosed by the surface with flux @xmath90 , then @xmath95 and @xmath96 are two further flux functions found , and the subscript denotes the partial derivative with respect to @xmath97 .
the ratio @xmath98 is another flux function which gives the pitch of the magnetic field .
the angles are cyclic on @xmath99 $ ] which determines the coordinate jacobian to be unity .
hence the metric tensor @xmath100 satisfies @xmath101 .
the function @xmath102 is found by the constraint that @xmath103 is cyclic on the unit interval @xmath104 where @xmath105 is a simple cyclic coordinate on the flux surface in the poloidal plane ( e.g. , path length on the curve , or the physical angle position about the magnetic axis with respect to any branch cut ) .
then @xmath106 is found @xmath107 where the angle brackets give the flux surface average defined as @xmath108 now , @xmath109 is also defined and the toroidal angle coordinate is given by @xmath110 which completes the prescription .    the field aligning is a one - to - one and onto coordinate transform from @xmath92 to @xmath111 given by @xmath112 where @xmath8 is a reference minor radius and @xmath113 is an arbitrary function of @xmath114 which is chosen to make the off diagonal perpendicular metric element @xmath115 vanish at a particular location .
this is called the shifted metric procedure @xcite .
the point is that at any position in @xmath116 where perpendicular drifts or @xmath86 is evaluated , the coordinate elements are rectangular , and the information contained in magnetic shear enters only in derivatives in the parallel coordinate .
the magnetic field components now satisfy @xmath117 at all locations in space for any choice of @xmath113 .
hence although the perpendicular coordinate elements are rectangular only for one particular location in @xmath116 the field aligning is exactly satisfied everywhere .
the coordinate metric elements @xmath118 are found in the @xmath92 representation .
then the grid locations @xmath119 are all given their own members of the family of these field aligned coordinates via @xmath120 then we have @xmath121 at @xmath119 .
the drift tensor elements are given by @xmath122 that is , @xmath123 is the only independent , nonvanishing component .
we have @xmath124 which defines the bracket @xmath125 $ ] for any scalar field quantities @xmath79 and @xmath126 .
with @xmath79 the electrostatic potential this gives the exb advection term , for electrons for which flr corrections are neglected by taking @xmath127 . for ions
this is generalised by the flr corrections as given above .    for local equations as those used here (
nonlinearities kept only in advection and magnetic flutter ) the geometry model must be modified to retain free energy conservation by the equations .
in particular , any function of @xmath114 such as a radially varying normalised parameter ( such as the temperature / charge ration @xmath39 ) multiplying any of the curvature terms will cause them to fail to conserve free energy .
the compressibility is already split away from advection and placed into the curvature terms , so the brackets must retain the properties of incompressible advection .
this means that the model for @xmath128 which is used must satisfy @xmath129 and to avoid any confusion the curvature terms are written with @xmath130 , which can be given any spatial dependence .
correspondingly , the mhd version of the continuity equation ( neglecting diamagnetic fluxes ) has already been manipulated according to @xmath131 @xmath132 and then terms such as @xmath133 are neglected due to the ordering .
the result is the same as the local continuity equation @xmath134 where @xmath135 is a normalising constant , if we identify the dependent variable @xmath136 with @xmath137 .
it is important to note that with the self consistent background being evolved this also includes the profile : @xmath138 .
the simplified geometry is now defined by its operators .
the laplacian and gyroaverging operators neglect @xmath87 .
the form in the ampere s law ( which does not involve the gyroradius ) is given by @xmath139 the form in the gyroaverging operations and polarisation is given by @xmath140 the brackets are defined by @xmath141=q[f , g]_{xy}+[f , g]_{xs}-q_{,x}(s - s_k)[f , g]_{ys}\ ] ] for any two scalar fields @xmath142 with each bracket piece defined by @xmath141_{ij } = \lp{\partial f\over\partial x^i}{\partial g\over\partial x^j }          - { \partial g\over\partial x^i}{\partial f\over\partial x^j}\rp\ ] ] with the third piece not used since these are always evaluated at @xmath119 .
nonlinear computations on ideal or peeling ballooning mode elm scenarios in tokamaks have previously been based on magnetohydrodynamic ( mhd ) models @xcite , and the respective codes had largely been tailored towards mode structures expected by linear analysis .
restriction to single - fluid mhd equations including two - fluid correction terms allows computations of large to mesoscale dynamics in realistically shaped 2d tokamak equilibria including x - point geometry crossing the separatrix , due to the lack of a need to use field - aligned coordinates .
recent mhd simulations are able to reproduce spatio - temporal mode structures , divertor footprints and other characteristics with good agreement to experimental diagnostics , in particular for the growth and immediate blowout phase of the finger - like instability @xcite .
however , as soon as the dynamics becomes nonlinear , the spectrum broadens to include scales normally associated with drift wave dynamics .
this turbulence is basically of a drift - alfvn type with strong edge itg character ( see ref .
@xcite for the signatures of the various mode structure types ) , made more powerful than otherwise by the energetic access to the long - wavelength mhd component .
the itg character results from the steeper logarithmic gradients in both temperature profiles , and the lack of involvement of ion temperature fluctuations in the parallel nonadiabatic response of the electrons . although edge turbulence of this type ( cf .
@xcite ) is strongly suppressed in the h - mode phase and is not initially involved in the instability phase , the experience of edge turbulence becomes relevant as these physical components become involved in all the nonlinear phases of the blowout .
the saturation and aftermath of the blowout , which in an experiment carries a large part of the actual transport losses that lead to a degradation of the pedestal , should be expected to involve physics not contained in the mhd model .
ion temperature gradient ( itg ) driven modes and ideal ballooning modes are both essentially caused by the gradient and curvature driven interchange instability and show similar character in their initial linear growth phases .
instability is achieved above their respective critical gradients , determined by the ratio @xmath143 between density gradient to ion temperature gradient scale lengths for itg , and by the ideal mhd ballooning parameter @xmath144 for ibm , where @xmath145 is the magnetic shear parameter and @xmath146 is a critical threshold which depends on beta , collisionality , and toroidicity . the principal difference between these modes is the relative role of the parallel alfvn responses which tend towards enforcement of an adiabatic response in the electrons ( parallel force balance , with the electron pressure gradient ) .
this is nonexistent in an mhd model and subdominant generally for an mhd instabililty .
but it constrains the electrons if the @xmath147 parameter is below threshold , still allowing the itg instability because the adiabatic response does not involve the ions .
the itg instability exists at all @xmath147 values , and furthermore for the edge situation the itg instability and mode structure transitions smoothly to a drift wave one in a nonlinear setting for @xmath148 .
hence there is actually no finite threshold in values of @xmath149 or @xmath147 for the nonlinear edge situation .
therefore , the existence of a threshold in the experiment is the same as the mechanism which maintains the h - mode , which is not yet well understood .
it follows that a gyrofluid model computation of an elm scenario with enough resolution to allow for ion - gyroradius based dynamics has to face the lack of a simple linear threshold in the parameter space .
the transition from initial ( micro-)instabilities to generic edge turbulence was studied in detail in ref .
as the most unstable linear modes crystallise out of an initial random bath of small - amplitude perturbations , the linear growth rate rises and becomes steady .
the maximum value of the instantaneous growth rate of total fluctuation free energy @xmath150 given by @xmath151 may be taken as the maximal linear growth rate .
the curve of @xmath152 then falls very sharply to zero ( over about @xmath153 ) as saturation occurs .
there is some structural adjustment over the next few @xmath154 as the spectrum fills out , and then the turbulence is fully developed .
but over the adjustment phase the value of @xmath155 is well below its previous maximum .
we will use the same diagnostic herein , except that the total ion exb heat flux @xmath156 averaged over the computational domain is used instead of @xmath150 because most of @xmath150 is represented by the self consistent profiles while @xmath156 is entirely due to fluctuations .
visualisation of the tokamak elm burst : the edge region of a circular poloidal tokamak cross - section is shown on the left with computational perpendicular @xmath116-domains ( @xmath157 2 , 5 , 8 , 11 , 14 of @xmath158 total ) mapped onto the circle during the maximum growth phase ( @xmath159 in units of @xmath9 ) of an ideal ballooning mode . on the right , the outer midplane area ( @xmath160 with size @xmath161 and @xmath162 ) is magnified and the radial domain stretched by a factor of 2 for better visibility at times @xmath159 ( max . flux ) and @xmath163 ( turbulent aftermath).,width=188 ]    the ibm elm blowout scenario is similar to this , initially , except that the instability is not a microinstability .
nevertheless , the scale differs by less than an order of magnitude : the toroidal mode numbers for linear itg instabilities are in the range of @xmath164 while the main ideal ballooning mode is near mode number @xmath165 , on the entire flux surface , for these typical @xmath166 values ( recall @xmath167 for @xmath168 , and the ratio @xmath169 is not far from unity in experimental cases ) .
the itg mode numbers are determined by the ion gyroradius and the dimensionless scale ratio @xmath170 , while the mhd values are determined by the width of the pedestal .
since @xmath5 is not larger than about @xmath171 in an h - mode pedestal , these scales are not disparate .
due to the closeness of the native scales , the ibm instability very quickly transfers energy to the ion gyroradius scale , which is only about two or three steps away in a turbulence cascade which acts at a factor of two at each step .
this brings the itg and drift alfvn physics into play , and in any nonlinear stages the latter involves stable shear alfvn component with substantial energy content .
the moment of saturation is defined as the time at which the ion exb heat flux averaged over the computational domain is maximum ; we refer to this as the `` peak - flux '' time .
the results shown below indicate that coincides with the establishment of the fully turbulent regime .
representation of this phase requires resolution of all scales and inclusion of the appropriate drift wave physics in the model .
how much resolution is actually needed depends on the strength of the blowout , which is situation dependent .
we therefore include a resolution scan on the toroidal mode number spectrum .
the result is that convergence is reached in the temporal phases just after peak flux only if the spectrum reaches down to @xmath172 .
the self - consistent equilibrium coupling in our gyrofluid model only allows treatment of shifted circular @xmath173 geometry ( the difference to mhd codes is that the high - resolution flr - gyrofluid dynamics necessitates field - aligned coordinates ) .
effects of flux surface shaping are postponed to later studies .
our focus here is on basic physics issues regarding the ability of the model to capture mhd phenomena ( see also ref .
@xcite concerning global alfvn oscillations with self consistent profiles ) , resolution , turbulent character , scalings and tests which are also accessible by @xmath173 geometry .
in particular , the theory of explosive instability has been formulated for circular geometry @xcite and is therefore accessible by the present model .
first - principles based local drift wave edge turbulence simulations are not able to obtain a realistic h - mode edge state with the known experimental characteristics : correct density and temperature pedestal profiles shapes or strength of flow shear are not obtained by self - consistent evolution by specifying core sources only , nor has a threshold transition character been found in any verified edge turbulence simulation @xcite . therefore some kind of `` modelling '' has to take place when the ibm instability ( as an h - mode phenomenon ) and its subsequent nonlinear evolution is simulated with a nonlinear gyrofluid turbulence code : although the realistic development of an edge transport barrier ( and thus a full elm cycle ) can not be directly obtained , one still may prescribe the h - mode pedestal profile before the onset of an elm , known from experimental data , as an initial state for the simulation . as a base case for the prescribed pedestal profiles
the well diagnosed edge characteristics of asdex upgrade h - mode shot # 17151 is used here @xcite .
the local parameters , taken as mid pedestal values , correspond to electron and ion temperatures @xmath174  ev , @xmath175  ev , densities @xmath176  m@xmath177 , magnetic field strength @xmath178  t , major torus radius @xmath179  m , aspect ratio @xmath180 , perpendicular temperature gradient length @xmath181  cm , density gradient length @xmath182  cm .
the profile of the safety factor @xmath183 is parabolic yielding local values at the last closed flux surface ( lcfs , @xmath184 ) of @xmath185 and @xmath186 .
the radial domain of the simulations covers a range of @xmath187 on either side of the lcfs .
the nominal pressure values are @xmath188 and @xmath189 .
all times are given in units of @xmath9 .
the computational grid is given in terms of the @xmath111 coordinate domains .
the spacing is equidistant . for each grid point in the third coordinate @xmath119 the @xmath190-coordinate is @xmath191 as defined in eqs.([eqcoords],[eqalphak ] ) so that for each grid point @xmath192 the coordinate system @xmath193 is the one which has @xmath194 at @xmath119 . for the nominal case @xmath162
perpendicular and @xmath195 parallel mesh points are used .
the radial domain with @xmath196 spans the plasma edge region between the h - mode pedestal top , with plasma core parameters as inner boundary values , and the outer bounded scrape - off layer region ( @xmath197 ) .
this represents a spatial range from the global scale to smaller than the ion gyroradius scale , with a ratio @xmath198 between ion gyroradius @xmath199 and minor torus radius @xmath8 . in the sol region ( @xmath200 )
the parallel boundary condition is replaced by a debye sheath model whose treatment is given in ref .
@xcite .
radial tokamak edge pedestal profiles of the plasma pressure before and after an elm blowout event across the separatrix into the scrape - off layer ( sol ) .
the area shaded in grey shows the present radial numerical simulation domain.,width=283 ]
the initial conditions are thus prescribed and are based on experimentally diagnosed radial temperature and density pedestal profiles @xmath201 and @xmath202 for each species ( electrons and ions ) .
a consistent electrostatic potential @xmath203 is derived by numerically solving the neoclassical balances ( parallel dynamics , toroidal drifts , and collisional dissipation , but not transport ) in a pre - processing step with a modified ( zonally frozen ) gemr setup , resulting in a time - steady 2d dissipative solution . the @xmath173 mhd equilibrium in turn
is set internally in gemr by computing the shafranov shift , and in the present setup would not be consistently described by prescribed in terms of external or coupled shaped equilibrium solvers @xcite .
the parallel and perpendicular electron and ion temperatures , @xmath204 , @xmath205 , @xmath206 and @xmath207 , are directly adopted and fixed from experimentally derived values by filtering the zonal component out of the total time derivative @xmath208 .
their values are zonally frozen by setting @xmath209 where @xmath210 is the zonal ( parallel and perpendicular ) average of all right - hand - side terms @xmath211 in the temperature evolution ( eqs .
5,6 ) including numerical dissipation .    on the other hand , gyrocenter densities @xmath212 and @xmath213 have to be set to obey relaxation relations that allow the vorticity to freely evolve into equilibrium .
this is achieved by freezing the zonal component of the sum @xmath214 , where @xmath215 , as the density part of the pressure during the equilibration phase , but allowing the difference ( i.e. , vorticity ) to evolve freely .
this ensures that the contribution of densities to the total pressure is zonally frozen through the relation @xmath216 while the densities relax regarding to @xmath217 and @xmath218    the numerical solution of the equilibration phase , starting directly from realistically steep pedestal profiles @xmath219 , into steady state is delayed by long , weakly damped global geodesic alfvn oscillations .
convergence is expedited by ramping up all of the gradients gradually from zero to prescribed value over the first @xmath220 of the run by @xmath221    this pre - processing equilibration phase is run until convergence with reduced perpendicular resolution ( @xmath222)=(@xmath223 ) , and without the exb and magnetic flutter nonlinearities , which allows establishment of the 2d structure in a smooth manner .
then , the resolution is increased to the nominal values ( @xmath224 ) , and a random turbulent bath with relative amplitude @xmath225 is added to the background pedestal profiles inside the closed flux surface region .
this procedure reduces transient alfvnic and geodesic acoustic ringing and prepares a reproducible initial state .
however , the following sudden release of the nonlinearities also leads to transient oscillations . depending on parameters , these may still be present at the onset of the instabilities , and can obscure a clean view on the nonlinear growth rates , which will be relevant for the discussion below on diagnosing linear or explosive instability .    a schematic sketch of the equilibrium pedestal profiles , representing an idealised asdex upgrade h - mode scenario @xcite , is shown in fig .
[ f : profiles ] together with the final relaxed state after the elm blowout phase ( which is discussed in the next section ) .
when this initialised pedestal pressure profile is ideal ballooning unstable , the ibm instability in gemr simulations is observed to be linearly growing in the pedestal region and at the onset of nonlinearity further overshoots and saturates , representing the beginning of the turbulent blowout phase during which a substantial fraction of the pedestal energy is thrown onto the sol .
subsequent evolution involves turbulence in both the pedestal and sol regions as the original pedestal energy is dissipated .
previous nonlinear approaches on elm ideal ballooning mode burst computations have treated only the initial growth and nonlinear phases , focused on low - wavelength modes and resolved only the mhd - relevant scales , excluding treatment or discussion of the ion gyroradius scales .
this however precludes the development of fully developed microturbulence which is caused by the onset of nonlinearity ( robust transfer of free energy to smaller scales ) .
ultimately , this microturbulence decides both the transition to nonlinearity and the eventual saturation ; that is , most of the post - peak transport curve .
this affects the mhd scenario of nonlinear explosive growth which most previous work advances .
however , we find that the exclusion of the ion gyroradius scales produces an under - resolved situation .
we measure the quantitative character of the growth in both linear and nonlinear phases with a growth curve @xmath226 defined in terms of the heat flux @xmath156 instead of the fluctuation free energy as explained above .
the heat flux is computed as a zonal ( flux surface ) average , @xmath227 @xmath228 and then the time trace @xmath229 , normalised in terms of @xmath230 at nominal parameters , is computed as a volume average over the central half of the radial domain .    the ibm instability , which follows the random seeding of the pre - processed pedestal profile in our computations , is very violent , growing for the nominal case at a rate @xmath231 , just below the ideal interchange rate .
on the other hand the growth curve @xmath152 appears qualitatively the same as in basic turbulence cases . at all time points in the nonlinear phase @xmath152
is well below its previous maximum . just after initial saturation
there is some nonlinear evolution in which @xmath152 crosses zero before settling down into long - term statistical saturation . at late times
the initial blowout no longer imprints the results : with a fixed source one merely finds bursty turbulence thereafter .
hence , there is no evidence for explosive instability in this nominal case , which , if present , would be visible in a finite time singularity in the fluctuating free energy @xmath232 rising with a power @xmath233 .    the burst and the resulting decaying turbulence act to transport the pedestal plasma across the lcfs , where energy and particles are lost within short times by parallel boundary outflow to the scrape - off layer ( sol ) limiter . without maintenance of a heat and particle source at the core boundary
, this leads to decay of the initial gradient with a rapid onset at the initial transition to nonlinearity at about @xmath234 .
this time point also corresponds to the peak - flux time .
total energy @xmath235 ( top figure ) and ion heat flux @xmath229 ( bottom figure ) during the elm blowout : without source maintenance ( dashed curve ) ; with the source term switched off at the time of maximum instability ( here : at @xmath236 ) ; and with the source maintained for the whole simulation time.,title="fig:",width=283 ]   total energy @xmath235 ( top figure ) and ion heat flux @xmath229 ( bottom figure ) during the elm blowout : without source maintenance ( dashed curve ) ; with the source term switched off at the time of maximum instability ( here : at @xmath236 ) ; and with the source maintained for the whole simulation time.,title="fig:",width=283 ]    however , the ion temperature profile already starts decaying during the low - amplitude linear growth phase , due to neoclassical ( finite larmor radius ) transport ( for detail on how this works in a fluid model see ref.@xcite ) .
this situation is shown in fig .
[ f : energy ] ( left frame , dashed line ) , where the total pedestal energy is reduced by a quarter of its initial content before the peak flux phase
. therefore we test sensitivity against sources .
maintenance of the ion temperature profile by a fixed ion heat source localised to the inner boundary ( following the method for driven cases in ref .
@xcite ) ensures nearly unchanged gradients during the linear growth phase . at the end of the simulation time ( @xmath237 ) after the elm burst
, the energy content is reduced to around 10% of the initial content when the source is absent or is switched off at the peak flux time ( @xmath234 ) , and to around 20% when the source is maintained for the whole simulation time .
the influence of either switching off the source or maintaining the zonal profiles on the nonlinear ion heat flux @xmath229 is shown in fig .
[ f : energy ] ( right ) : the linear growth phase and the peak flux is not significantly changed by the source .
as expected , the heat flux saturates for late times with maintained core boundary inflow and is otherwise decaying .     consistency check for different perpendicular resolutions by varying the grid point number @xmath238 between 64 and 512 .
convergence is seen for @xmath239 and higher.,width=283 ]    consistency checks on the perpendicular resolution are made by varying @xmath238 between 64 , 128 , 256 and 512 .
the evolution of @xmath229 for these cases is shown in fig .
[ f : nyres ] .
the @xmath240 case is clearly under - resolved and overestimates the total energy decay by a factor 2 . for @xmath241
the initial part of the post - peak decay phase disagrees with the @xmath162 case .
convergence is found for @xmath239 and higher .
the nonlinear saturation phase thus requires resolving the ion gyroradius scale .
this is the clearest indication that turbulence for toroidal mode numbers beyond about @xmath242 extending down to the @xmath243 range is involved in the blowout saturation process .
perpendicular mode number @xmath244 spectra for @xmath159 ( top ) and @xmath163 ( bottom ) around the peak flux phase .
the vorticity spectrum ( bold lines ) is already flattening down to the ion gyroradius scale for @xmath163 only a few @xmath245 after the maximum linear growth phase.,title="fig:",width=283 ]   perpendicular mode number @xmath244 spectra for @xmath159 ( top ) and @xmath163 ( bottom ) around the peak flux phase . the vorticity spectrum ( bold lines )
is already flattening down to the ion gyroradius scale for @xmath163 only a few @xmath245 after the maximum linear growth phase.,title="fig:",width=283 ]    the transition from linear instability to turbulence is studied during the peak flux phase around @xmath234 .
[ f : kyspecs ] shows perpendicular mode number @xmath244 spectra of the squared amplitude of various fluctuating plasma quantities ( density @xmath246 , ion temperature @xmath247 , electrostatic potential @xmath79 and vorticity @xmath248 ) for @xmath159 and @xmath163 .
this time difference corresponds to @xmath249 in physical units and around @xmath250 in local drift units , which is only slightly faster than the overshoot and saturation times known from edge microinstability cases @xcite .
initially , the ion temperature gradient ( itg ) driven microinstability and the ideal ballooning mode ( ibm ) compete in growth out of the random low - amplitude bath .
the itg mode grows strongest near the separatrix due to radially local steepening by parallel sol diffusion , which in our simulations may be overestimated by using the standard fluid bohm outflow boundary conditions .
for @xmath251 the linear ibm is clearly dominant near a toroidal mode number of 9 - 10 for the nominal aug parameters , consistent with experimental observations @xcite . around @xmath252
rapid formation of a turbulent cascade range in the spectra is observed and the vorticity spectrum is already spread out to the ion gyroradius scale ( @xmath253 ) .
this is a manifestation of the role of self generated drift wave turbulence in the saturation process .
because of the way both the diffusive mixing and vorticity scattering nonlinearities enter the drift wave physics @xcite the spectrum is held together as a unit , all scales down to @xmath0 are involved in the saturation phase of the overall elm blowout transport .
the involvement of @xmath0 at unit order , in turn , necessitates a nonlinear gyrofluid model ( and at some later date , a nonlinear gyrokinetic model ) .
the characteristics of the blowout at various plasma beta is studied for @xmath254 , with the absolute @xmath255 a factor of @xmath256 higher due to the @xmath257 value . a competition between unstable growth of itg high - n modes and ibm low - n modes
is observed : for low @xmath255 an early onset of itg turbulence degenerates the profile and influences the strength of the ibm . for large @xmath255
the nonlinear saturation is acting faster and stronger .
the two effects lead to a rather nondistinctive result concerning the peak elm flux as shown in fig .  [
f : beta ] , with a continuous transition between itg and ibm triggered onset of the turbulent aftermath .
competing growth of itg , microtearing and ibm instabilities determine onset times and intensities of the blowout time traces @xmath229 when scaling with the local plasma beta : for @xmath258 only an itg micro - instability is growing around @xmath259 and saturating on a low h - mode transport level . for @xmath260
the ibm instability dominates and leads to a clean elm signal around @xmath261 . in the nominal case ( @xmath262 ) itg and ibm mhd growth
compete initially with similar growth rates but the ideal ballooning mode takes over around @xmath263 and further determines the blowout . for @xmath264
a microtearing instability is saturating on a gyro - bohm transport level around @xmath261 and is transformed into more violent itg mhd turbulence after deterioration of the initial electron temperature gradient ( @xmath265 ) .
, width=283 ]    the presence or absence of background current gradient terms , which in this model are set by replacing @xmath266 everywhere the electron parallel velocity @xmath267 appears in the equations , with @xmath268 given by the @xmath106 profile , was found to have no discernible effect on the result .
indeed self - consistent inclusion of @xmath269 effects with possible role in the l - to - h transition was among the motivations of the gemr model in the first place @xcite .
however , the magnitude of even an impressive pedestal current peak is only one to three times the nominal saturation current given by @xmath270 , while values closer to @xmath271 are required to enter the energetics effectively .
since @xmath272 for the nominal case , current gradient effects are very weak .
the main conclusion of this study is that the qualitative nature of the saturation and aftermath of the initial ibm blowout is the same as for generic edge turbulence given a small - amplitude start .
the transition from linear mode structure and energetics to turbulence found for these blowout cases is the same as in ref .
only the nature of the linear mode itself differs .
the blowout saturates upon its own self generated drift - alfvn turbulence , with a strong ion temperature component given the gradients .
the vorticity spectrum reaches quickly to the ion gyroradius ( @xmath0 ) scale , requiring the gyrofluid model and explaining the numerical difficulties seen with braginskii models  it can be argued that the cases given in ref .
@xcite crash on entry to the nonlinear stage .
convergence in the aftermath requires resolving at least @xmath0 .
unfortunately , due to the lack of a self consistent h - mode state in a well resolved computation , no threshold is found . at lower beta values
one simply finds generic edge turbulence driven by the temperature gradients .
it is not clear that this scenario really describes an actual elm , although the energetic growth and decay curves are not unrealistic .
however , with a well resolved transition to nonlinearity in both the energetic peak and aftermath phases , our studies find no evidence of an explosive mhd phenomenology .
indeed it can be argued that nonlinear mhd processes are pre - empted by the efficient transfer to smaller scales through the two - fluid drift wave physics more commonly associated with microturbulence .
this work was supported by the austrian science fund fwf under contract y398 , by a junior research grant ( `` nachwuchsfrderung '' ) from university of innsbruck , and by the european communities under the contract of associations between euratom and the austrian academy of sciences , carried out within the framework of the european fusion development agreement . the views and opinions
herein do not necessarily reflect those of the european commission .
h. strauss , l. sugiyama l. , c.s .
chang , et al .
21st iaea fusion energy conf .
( chengdu , china , 2006 ) ( vienna : iaea ) cd - rom file th / p8 - 6 and http://www-naweb.iaea.org/napc/physics/fec/fec2006/html/index.htm