Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: 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 ) .
Abstract:

Answer: we analyze the phenomenological implications of introducing vector - like leptons on the higgs sector in the higgs triplet model . we impose only a parity symmetry which disallows mixing between the new states and the ordinary leptons . if the vector leptons are allowed to be relatively light , they enhance or suppress the decay rates of loop - dominated neutral higgs bosons decays @xmath0 and @xmath1 , and affect their correlation . an important consequence is that , for light vector leptons , the decay patterns of the the doubly - charged higgs boson will be altered , modifying the restrictions on their masses .


Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: convection is a major process taking place in reg giant ( rg ) envelopes . in a previous letter ( @xcite , hereafter at97 ) , results of a study of this subject were reported briefly . in the present paper , more results of this study are presented .    throughout the paper , the reference one - dimensional ( 1d ) model of convection is the mixing length theory ( mlt , @xcite ) , according to the prescription of @xcite . laboratory experiments , observations of the outer layers of the solar convective zone , and theoretical understanding have demonstrated that convection is much more complicated than mlt assumes . analysis of experiments of convection in a box , in which a temperature gradient was established by using helium at about 5k and controlling the boundary conditions , have shown that convective flow is non - local with defined narrow streams going from the top down ( and vice versa ) , i.e. the existence of a downstream at a given depth at a given time depends upon the existence of a downstream in a higher point at previous time and not upon local parameters only . this phenomenon is known as hard turbulence ( see @xcite ) . observation of the outer layers of the sun have shown the existence of two dimensional cells of upstreams , in a variety of sizes and shapes , with the boundaries of these cells being the downstreams , i.e. the flow is not symmetric ( see @xcite ) .    during the years , many studies were made in order to overcome the simplifications of mlt . several non - local and non - instantaneous theories were proposed ( @xcite ) . a common feature in these theories is using turbulence theory techniques in which the hydrodynamical quantities at each point are constructed of two components : a constant value and a varying value . equations are formulated for each of the varying quantities which are then used to determine the convective velocity and flux . these equations include creation and diffusion terms and thus constitute a set of coupled partial difference equations . maybe because of this complexity , there is no wide usage of these theories . xiong ( 1998a , b ) have used one of these theories to compute linear pulsation times of rr lyrae and long period variables . @xcite are using a simplified 1d model of turbulent energy diffusion and convection with eight free parameters to describe cepheid and rr lyrae variables . @xcite used turbulence models to construct an alternative model for mlt which includes a full spectrum of eddy sizes ( see also @xcite ) . this model is still local and instantaneous , but is easy to include in astrophysical codes . it was found that this model yields higher values of convective flux than mlt at efficient convection regimes , and lower values of convective flux at inefficient convection regimes . this model was found to be preferable to mlt in many studies ( see @xcite ) . these studies included red giant evolution ( @xcite ) , heliosesmology ( @xcite ) and calculations of uvby colors of a f and g stars ( @xcite ) . the superiority of this model , for example in modeling red giant evolution , was the usage of a constant parameter for various masses , while when using mlt there was a need to change the mixing length parameter with mass . however , @xcite have shown that a varying calibration is needed for this model too . alternatives to the usual mixing length used in mlt ( which is a free parameter times the pressure scale height ) were suggested . these included using the density scale height ( @xcite ) and the distance from the boundary of the convection zone as in @xcite .    for a correct modeling of red giant variables , the lack of time dependence and the locality of mlt are major disadvantages . the convection velocities , as predicted by mlt , and the radial velocities associated with pulsation are of the same order of magnitude ( of @xmath4 ) . thus , interaction of the two phenomena is expected . ad hoc prescriptions were made in order to take into account time dependence in dynamical calculations with mlt ( for example , pulsation calculation of miras , see @xcite ) . in these calculations , the convective flux at a given time is determined by using mlt and the flux at a previous time assuming a typical convective time scale estimated by convection velocity . in other calculations , like evolution of massive stars , a non locality is presented to mlt in particular when treating penetration ( or overshoot ) of convection to stable neighboring zones . a common assumption is that this overshoot region is proportional to the mixing length at the edge of the convection region . these methods are obviously not rigorous and not intended to describe these processes precisely ( for example see @xcite . being mainly a hydrodynamical phenomenon , an appealing method for modeling convection is by using numerical multidimensional simulations . in this approach , no explicit model of convection is needed , and thus no mixing length parameter is assumed . time dependence and non locality are naturally present . however , convection is multidimensional and the flow has many length scales and time scales , that might be smaller by far than other typical scales of the problem . this is why such calculations are not suitable for all problems with convection . numerical simulations of the hard turbulence experiments mentioned above where made ( werne 1993,1994 , @xcite ) . these calculations reproduced the typical narrow streams of the experiments , and revealed that the origin of the flow was the thin boundary layers at the bottom and at the top of the experimental box . many studies were made in order to understand general features of convective flow ( @xcite ) . typically , these studies assumed simplified physical assumptions , such as ideal gas equation of state , constant opacity and constant gravitation . the systems were confined in a box with an initial unstable polytropic profile . boundary conditions were periodic on the sides of the box , no flow was assumed in and out of the bottom and top of the box , a constant temperature gradient at the bottom and a constant temperature at the top . some of the simulations assumed planar symmetry using 2d codes , while other used three - dimensional ( 3d ) codes .    summary of some the results of these studies is : there is an asymmetry in the flow with narrow rapid downstreams and wider upstreams and the 3d picture resembles the observations of the sun with cells of upstreams ; the flow in the wider upstreams is somewhat more turbulent ; the size of the cells increases with depth so downstreams merge . this picture is revealed when using high resolution , while low resolution simulations show laminar flow with big convection cells . the internal energy flux in the downstreams is relatively large , but a large kinetic energy flux with opposite sign almost cancels the net energy flux in the downstreams . near the outer boundary , convection might be supersonic . other `` simulation in a box '' studies were dealing with invoking of p - modes in the convection zone ( @xcite ) , producing g - modes below ( and above ) the convection region ( @xcite ) and penetration of convection to neighboring layers ( @xcite singh 1994,1995 ) . attempts have been made to simulate the outer layers of the solar convection region ( stein & nordlund 1989,1998 , @xcite ) . since the computation box in these simulations included only a fraction of the convection region , penetrative boundary conditions were applied at the bottom of the computation box . the outer layers which were included in the simulations and in which a superadiabatic gradient exist , are the most sensitive to the convection model . these simulations reproduced the granulation pattern of convective streams of solar convection , as well as the heat flux and spectrum ( @xcite ) . comparison of various measurements form the simulations ( such as maximal temperature gradient near the outer boundary ) enables calibration of the mlt parameter . @xcite and @xcite have published such calibration for various models of solar - type stars . they found that mlt parameter is 1.2 to 1.8 for starts with effective temperatures of 7100k to 4300k ( with dependency upon g - gravity and metallicity as well ) . they also investigated the @xcite model and found it needs varying calibration as well . simulations of convection in zz - ceti type white dwarfs were done by @xcite . they found that different calibration models yield different mlt parameter values : 1.5 according to one measurement , and 4 according to another .    two - dimensional ( 2d ) calculations of rr - lyra were performed by deupree ( 1975,1977 , 1985 ) ; the ( very ) low resolution simulations of the narrow convection region in these variables stars have shown that convection has a damping effect upon pulsation , and vice versa . in these stars , there is little difference between a radiative model and a convective model , due to the narrow convection region . thus , the relaxation time in the 2d simulations is relatively small , and might be calculated . this is not true for rg in which the convection region is wide , and that is why simulating convection in rg s envelopes is much more difficult . the above discussion , as well as the topic of this paper , concerns convective envelopes , in which the convection is not efficient enough to enforce an almost adiabatic profile throughout the convective region . multidimensional studies of other astrophysical cases have been published . this includes core convection ( @xcite ) , shell burning ( @xcite , bazan & arnett 1994,1998 ) , nova bursts ( @xcite ) , core collapse ( see @xcite and ref . therein ) , and accretion disks ( for example @xcite ) . in these studies , some of the interesting questions are : the amount of mixing and penetration due to convective streams , nuclear energy sources and sinks in the convection region , the coupling of convection and rotation . an interesting question in the context of multidimensional simulations of convection , is the validity of 2d calculations . while it is obvious that 3d simulations are preferable , it is not clear what are the quantitative differences between 2d and 3d results . in many 2d studies of convection , there are large convective cells ( eddies ) that occupy the whole convective region ( see for example @xcite ) , this phenomena is not seen in 3d simulations . however , it does not mean that there are no small size eddies in 2d simulations , when increasing spatial resolution large eddies tend to brake up to small eddies ( see @xcite ) , and even low resolution 2d simulations have small eddies ( see @xcite and discussion in @xcite ) . the qualitative features of convective flow , revealed by 3d simulations , are present in 2d simulations as well ( see @xcite and ref . therein ) . there are few quantitative measurements of the differences between 2d and 3d simulations : in a study of convection in proto - neutron star , @xcite have shown that the importance of small scale features , and the kinetic energies are different when comparing 3d to 2d simulations ; @xcite have simulated novae in 3d and compared the results to their 2d simulations ( @xcite ) , they found the difference in the spectrum of eddy sizes to have an influence on the amount of overshoot and mixing , thus yielding a different answer to the question of the existence of thermonuclear runaway ; in a study of convection in solar type stars , which is more comparable to this paper , @xcite have results of first 3d calculations , these results indicate little difference between the calibrated mixing length parameter from 3d simulations , relative to the calibrated mixing length parameter from 2d simulations . modeling the convection region of rg envelopes is necessary when studying phenomena related to these stars . for example , consider the observational relation of the pulsation time to the luminosities of mirae : there is ambiguity concerning the pulsation mode ( is it the fundamental mode or the first mode ? ) . the period in the first mode is about half of the period in the fundamental mode , however changing the mixing length parameter from 1 to 1.5 is enough to change the calculated period by a factor of 2 ( @xcite ) . the value of the mixing length parameter affects the calculated age of globular clusters ( see @xcite and ref . therein ) .    due to the importance of modeling convection in studying rg s , and the lack of previous hydrodynamical calculations of this phenomena , i focused on 2d simulation of rg s envelope convection . in 2 i describe the code used for the simulations . summary of the results published in at97 is presented in  3 and more results of the at97 simulations are summarized in  4 . a survey of the importance of the various parameters of the simulation techniques is presented in  5 and  6 . a short summary of this research is given in  7 . as was described in at97 , the hydrodynamical code vulcan ( @xcite ) was used for this study , after utilizing several modifications to the code . this version of the code solves the hydrodynamic equations and radiative diffusion equations by splitting each time step into hydrodynamic time step followed by a radiative diffusion time step . both steps are solved implicitly to have an accurate and stable solution . the code uses an arbitrary lagrangian eulerian ( ale ) scheme in which each lagrangian step is followed by a relaxation of the numerical grid . in this relaxation stage , the grid is modified to eliminate distortion and to allow calculation of flow with vorticity . a quasi 1d lagrangian relaxation was used , in which each radial row of cells kept its mass . this was done for each row of cells by calculating an average outer radius , so that the sum of the mass fluxes through it would be zero . using this method , one was able to follow the average radial expansion or shrinking of the envelope . the mass , momentum and energy fluxes in the relaxation stage were calculated by a second order donor scheme , in which the gradiants of these variables were used to estimate the fluxes at the boundaries of the cells ( according to @xcite ) . . like other numerical studies of convection ( @xcite ) , i implemented subscale viscosity and diffusion terms using the @xcite model . the inclusion of these terms in the equations solved by the code should compensate for the use of a finite number of cells which is much less than needed when simulating turbulent flow if one wishes to simulate the effects of all scales up to the scale of physical dissipation . in the simple smagorinsky model these terms are the usual viscosity terms with a coefficient which depends upon derivatives of the velocities and a typical length scale of the numerical grid . when there is a scalar field that is characteristic of the fluid , one should add a diffusion term of this scalar field with a coefficient which is proportional to the viscosity coefficient . in convective flow , this scalar is the entropy ( @xcite , ludwig 1994 , @xcite ) and thus entropy diffusion was added to the code . the two parameters in these large eddies simulation ( les ) terms were the smagorinsky coefficient ( which multiplies the viscosity coefficient ) and the prandtl number ( which is the ratio of the viscosity coefficient to the diffusion coefficient ) . both the viscosity and the entropy diffusion were added as explicit terms since they were relatively small . the inner boundary in the simulations was located deeper than the convection zone to minimize the effects of not including the entire envelope . however , in preliminary simulations , hydrodynamical waves were invoked below the convection zone , and there was a convective flux near the inner boundary . by comparing simulations in which the inner boundary was placed at different depths , it was found that this convective flux was not physical because it was always located near the inner boundary . realizing this was due to an improper boundary condition ( in which only the radial velocities were set to zero for points on the inner boundary ) , i changed the condition and set both components of the velocities at the inner boundary to zero and allowed only radial velocities at the lowest three rows of points . with this new boundary condition , the flux disappeared . the inner boundary condition for the radiative diffusion was of a constant entering flux equal to the luminosity of the static model ( i.e. equivalent to @xmath2 ) .    for the hydrodynamic step a constant ( usually zero ) outer pressure was used , without limiting the direction of the velocities at the outer boundary . however , at the mesh relaxation stage , an average radius was calculated for the outer boundary so that the mass flux through the boundary was zero . the remapping fluxes ( i.e. the fluxes at the mesh relaxation stage ) , through this boundary were calculated using the values of the bundary cells ( i.e. , using an acceptor scheme for the incoming flow ) . for the radiative diffusion , an outgoing flux was calculated using @xmath5 where t is the temperature of the boundary cell and @xmath6 is the stefan constant ( i.e. , assuming it to be at optical depth of about two thirds ) .    for the sides of the computational sector either reflective boundary conditions , or periodic boundary conditions were used . when using the reflective boundary conditions , all fluxes through these boundaries were set to be zero , and thus only downstreams or upstreams existed on the boundaries . this was different when using periodic boundary condition , for which each boundary cell has an artificial neighbor at the other side of the sector ( though , since i assumed cylindrical symmetry , there is no direct physical meaning to periodic boundary conditions ) .    in order to analyze the results of the simulations i used printouts of time averaged energy fluxes for each radial layer . the fluxes included radiative flux as was calculated by the radiative diffusion step , remapping fluxes of internal energy and kinetic energy calculated by the grid relaxation stage , and hydrodynamical fluxes that can be calculated using @xmath7 where @xmath8 is the total energy density , @xmath9 is the pressure , @xmath10 is the velocity vector and @xmath11 is the mass density . the convective flux is the sum of the remapping and the hydrodynamical fluxes . both the radiative and the remapping fluxes were calculated explicitly in the simulations and had only to be time averaged . however , the hydrodynamical fluxes were intrinsic and thus had to be estimated . i used two methods : a. explicitly using the above formula , and b. by energy balance of each radial layer in the hydrodynamical step ( i.e. the flux from one layer to the next in each time step is equal to the change in total energy of the layer during the time step minus the flux from the previous layer ) . both methods had inconsistencies : the first uses a formula which is not one of the formulae that control the simulation and the second is correct only if full conservation of energy exists in the simulation . additional discussion of this issue appears in  4 , where a flux plot is presented .    as a test case for the simulation techniques i reproduced the results of a `` convection in a box '' study ( published by @xcite ) . the configuration for these simulations was of two stable layers with an unstable layer in between . the velocity field and the averaged fluxes in the simulations were similar to the published existing results . in addition to the 2d hydrodynamic code , two 1d codes were used : a static code ( see @xcite for description ) , and a hydrodynamic code ( @xcite ) . the static code was used to integrate initial models of the envelope , and the dynamic code was used to follow the hydrodynamic evolution of the models for comparison with the 2d results ( as explained in the following sections ) . both codes use @xcite prescription of mlt to calculate convective flux . the equation of state and opacities used by these 1d codes are identical to those used by the 2d code . the equation of state is based on tabular form and includes sum of components for electrons , ions , molecules of hydrogen and helium and radiation ( see @xcite for details ) . opacities are calculated according to @xcite with water molecules contribution according to @xcite ( taken from @xcite ) . the study was focused on models of an envelope of a red giant with total mass of @xmath0 core mass of @xmath1 , luminosity of @xmath2 and metallicity of @xmath3 . the models were integrated using the mixing length theory with the mixing length proportional to the pressure scale height . five models were integrated using different values of the proportion coefficient ( -mixing length parameter ) @xmath12 : @xmath13 ( i.e. no convection ) , @xmath14 , @xmath15 , @xmath16 , and @xmath17 . these five models had different total energy as well as different energetic structure as an outcome of different energy transport efficiency , as was dictated by the value of the mixing length parameter . the aim of the study was to find which of these five models was more consistent with the results of 2d dynamical simulations . the final configuration reached by the end of 2d dynamical simulations of each of these five models should be the same , since they all correspond to the same physical envelope ( the same mass , luminosity and structure below the convection zone ) . the difference in the simulations is in the relaxation time needed for the different models to adjust their total energy and structure and to reach thermal equilibrium . due to the long thermal time scale of these envelopes , and to computational limitations , this relaxation phase could not be fully simulated , however , as was demonstrated in figure 1 , there was enough data in the beginning of the simulations to reach conclusions about the preferred model .    in this figure the luminosity as a function of time for each of the five simulations was presented . the most obvious result in this figure is that in the simulations that started from models with @xmath18 the envelope lost energy while for simulations started from models with @xmath19 the envelope gained energy . that was the main reason for the conclusion that the 1d model integrated with @xmath20 was the most consistent with the 2d simulations . a consistency check of this result was done by using 1d dynamic simulation with @xmath20 mlt starting from the above different initial models . it was found out that the 1d simulations mimic the 2d simulations , i.e. the luminosity ( and radius ) as a function of time for the different initial models looked similar . a more precise analysis of the results revealed that the final configuration of the _ dynamical _ simulations ( both one and 2d ) had slightly less total energy than the _ static _ 1d model of @xmath20 . a typical velocity field in the simulations is presented in figure 2 . one can identify narrow downstreams and wider upstreams in the flow . the length scale of the streams increased with depth ( as well as the pressure scale height ) . the existence of hydrodynamical waves below the convection zone is apparent .    as mentioned before , by comparison of the five simulations a conclusion could be made that the @xmath20 model was the most consistent with the final 2d configuration , however from figure 3 it is clear that this model was not in a full thermal equilibrium . in this figure , the ( time ) averaged fluxes in the 2d simulation of this model are presented , including the radiative flux and the total flux ( the sum of radiative and convective fluxes ) . for comparison , the 1d ( static ) radiative flux is presented , too . the radiative flux in the 2d simulation resembled the static model , though a modest widening of the convective zone is clear , and an increase in the radiative flux above the constant luminosity is apparent near the bottom of the convective zone ( see @xcite ) . the value of the total flux is indicative of the thermal equilibrium of the model . the total flux was slightly less than @xmath21 below the convective zone , and reached a maximum of about @xmath22 near the upper boundary and then it declined to the surface luminosity of about @xmath23 . while it is clear that this is a result of a yet non equilibrium configuration , one should nt forget that due to the lack of algebric energy conservation , this plot is influenced by the uncertainty in the value of the hydrodynamical flux ( as was described in  2 ) : the energy conservation in each time step was of the order of @xmath24 of the energy . this seems to be small , but when one divides it by the time step ( about 80 seconds ) , one gets @xmath25 . i think it is important to note that this has little effect upon the simulation itself , which i confirmed by adding artificial non - energy - conserving term to the _ 1d _ dynamic code and comparing the results of several simulations where i changed the value of this term . the different values of the luminosity in the simulations ( figure 1 ) were dominated by the values of of the outer radii ( compare this figure with figure 3 in at97 ) . in the initial configurations there was a relatively large difference in the effective temperatures : 3675k for @xmath26 , 4600k for @xmath27 and 5110k for @xmath28 . however , in the dynamical simulations , very quickly ( less than a month ) the difference shrank dramatically , and the effective temperature became : 4520k , 4600k , and 4650k ( for the same models as above ) with time fluctuations of about 100k . in other words , using the dynamical simulations , one gets almost immediately , an effective temperature which is very close to the effective temperature of the @xmath20 model . this gives more confidence in the above conclusion that the @xmath20 is the preferred model . the 1d consistency simulations featured the same phenomenon . the five ( 2d ) simulations used similar numerical parameters . the grid consisted of 36 angular cells in each row , with 53 to 72 rows ( for the simulations of models with @xmath19 , i included deeper regions of the envelope - down to eight @xmath29 - because of the deeper convection zone ) . the angular sector was @xmath30 near the equatorial plane . the fluxes in the relaxation stage were calculated using a second order donor scheme , and the energy flux consisted of both internal and kinetic energies . i used reflective boundary conditions on the sides of the computational domain . les viscosity was used , with smagorinsky coefficient of 0.5 and without entropy diffusion . i made many tests of the dynamical simulation of the model that started with the @xmath20 profile . i changed the numerical scheme , the number of cells and the angular width of the computational domain and the parameters of the les terms . the changes in numerical scheme included : calculating the remapping fluxes using full donor ( i.e. without taking the gradients of the hydrodynamical variables into account ) , using internal energy fluxes in the grid relaxation stage ( without adding the kinetic energy ) and the use of periodic boundary conditions at the sides . i found that the full donor simulation yielded a rapid expansion of the envelope - implying it is necessary to include gradients in the calculation of the fluxes . the use of the internal energy ( only ) for the remapping energy fluxes caused a small constant loss of kinetic energy and thus a modest shrinking of the envelope without changing the luminosity . the velocity field of the simulation with periodic boundary conditions is presented in figure 4 . the non radial velocities at the sides are expected , but no other significant change can be seen , compared to the velocity field of reflecting boundary conditions ( figure 2 ) , i.e. the size of the eddies and the changes in the velocities with depth look similar . the luminosity and outer radius as a function of time look the same as well , as can be seen in figure 5 . changes in the angular width or number of cells had small effects . increasing both the angular width to @xmath31 and the number of cells in each row to 54 yielded similar results as the nominal simulation , as presented in the velocity field - figure 6 , and in the luminosity and outer radius - figure 7 . when i decreased the angular width to @xmath32 i got a larger variability in the luminosity and outer radius ( see figure 8) and a small increase in the average luminosity , as well . a small increase in the average luminosity was seen in a simulation with larger number of cells ( = 54 ) in each row ( with the nominal angular width of @xmath30 ) , as presented in figure 9 . increasing the number of rows to 81 , had no significant effect . the variations in the les terms included changing the numerical coefficient and adding an entropy diffusion term . when i had zero or small ( smagorinsky coefficient of 0.01 ) les viscosity i got somewhat larger variability in the luminosity and outer radius , but without changes of the average values . however , including an entropy diffusion term ( with prandtl number of one ) increased the average luminosity by a small amount ( see figure 10 ) . when i had this term , a further increase in the luminosity was evident when the smagorinsky coefficient was increased to one a summary of these tests reveals that there was a tendency to increase the luminosity in two cases : when i increased the angular resolution , and when i added an entropy diffusion term in the les scheme . these two cases have one physical meaning : an increase in the efficiency of convection by adding smaller scale turbulent fluxes . since there was a very small luminosity increase in the above simulations , i performed yet another simulation where i increased the resolution by a ( relatively ) large factor . i started with a profile of the nominal simulation ( of @xmath33 cells ) and increased the number of cells almost by an order of magnitude ( @xmath34 cells ) , without the entropy diffusion term . this simulation , being much more cpu demanding , was followed for a small physical time but nevertheless a clear conclusion can be made , as can be seen in figure 11 . the luminosity increased quite significantly from about @xmath23 in the nominal simulation to about @xmath35 . an increase in luminosity means that the envelope loses more energy , and from the five simulation with different initial models , this means that the 1d model which is consistent with the 2d simulations should have a larger value of @xmath12 , i.e. @xmath19 . to test this conclusion i performed high resolution simulation starting with the initial model based upon @xmath28 . again , i started with a profile of the nominal simulation ( with @xmath36 cells ) and increased the resolution to @xmath37 cells . this yielded an increase in the luminosity from about @xmath38 to about @xmath39 , as expected . in order to check systematically the dependence of the luminosity in the spatial resolution and in the les terms , i performed the following simulations , which started with the same initial profile . the initial profile was a 2d profile with an already developed convective flow based on the initial 1d model with @xmath20 . from this profile i interpolated four equivalent profiles with different spatial resolution : @xmath40 cells , @xmath41 cells , @xmath42 cells and @xmath43 cells , thus having factor of about 64 between the number of cells in the finer grid and in the coarser grid . from the above discussion it is clear that the significant term in the les scheme is the entropy diffusion term , so i used the nominal 0.5 value of the smagorinsky coefficient , and changed the importance of the entropy diffusion term by changing the prandtl number , using the values : @xmath44 ( i.e. the last value corresponds to increasing the diffusion term by a factor of five compared to the first value ) . i used these values of the prandtl number for each of the four resolutions , and thus having sixteen simulations . the flow in four of these sixteen simulations is presented in figure 12 , for different resolutions . the flow is presented by contours of the vorticity and indeed increasing the resolution yielded more small scale eddies ( and higher values of vorticity ) . the luminosities of these simulations are summarized in figure 13 , and in table 1 i present the average luminosities . from these results one can conclude that the best estimate of the luminosity is @xmath45 , and that the better the resolution is , the less important the value of prandtl number is . in a perfect les scheme one should get the converged result for each resolution for given set of numerical parameters ( i.e. prandtl number in our case ) . however , the value of about @xmath46 was obtained in each resolution with a different value of prandtl number , this value seems to increase when decreasing the resolution .    for two of these 16 simulations i checked the importance of the les entropy diffusion term comparing to the les viscosity terms . i performed simulations where i changed both smagorinsky coefficient and prandtl number so that the amplitude of the diffusion term would be the same , while changing the amplitude of the viscosity terms . i performed a @xmath42 cell simulation with smagorinsky coefficient of 1.5 and prandtl number of 1 and got an average luminosity of @xmath47 ( equivalent to a prandtl number @xmath48 simulation ) , and a @xmath49 cell simulation with smagorinsky coefficient of 1 . and prandtl number of 1 where i got an average luminosity of @xmath46 ( equivalent to a prandtl number @xmath14 simulation ) . the summary of these results is that the 2d simulations of an initial model calculated with @xmath20 mlt , yields an energy - losing envelope with an average luminosity of @xmath46 . to find what value of @xmath12 can reproduce this i performed several 1d dynamic simulations . in each of these simulations , i started with the same @xmath20 profile , and assuming different values for the @xmath12 in the dynamical model i got different luminosities . for the following values of @xmath12 : 1.2 , 1.3 , 1.4 and 1.5 i got average luminosities of about : 230 , 240 , 250 and 270 ( @xmath50 ) . from this i conclude that the consistent value of @xmath12 is 1.4 . as was presented in at97 the convective envelope of a red giant ( of mass @xmath0 , core mass of @xmath1 , luminosity of @xmath2 , and metallicity @xmath3 ) was studied by comparing 2d simulations of five different initial models . while a full convergence of these simulations to a balanced final configuration was not possible , due to the large thermal time scale of this envelope , it was found that an initial model calculated with mlt using @xmath20 , was almost energetically balanced , and thus was the best estimate of the final configuration .      after a comprehensive survey of the importance of the various parameters used in the simulation technique , a sensitivity to the spatial resolution and to the les entropy diffusion term was found . this sensitivity was tested by using four resolutions and four values of the les term , and it was concluded that the preferred value of @xmath12 is 1.4 , for the envelope tested in this study . it was found that , at least for this code and for these kind of simulations , one needs to calibrate the les terms .    the conclusion for the value of @xmath12 , to be used in modeling the ged giant envelope studied in this research , being larger than one is consistent with the results of @xcite where they calculated evolutionary sequences and compared them to observation . this value of 1.4 is also consistent with _ extrapolation _ of the results of @xcite and @xcite where numerical simulations were applied to the outer convective zone of solar like stars . this conclusion is important for studies in which an accurate modeling of red giants is necessary , like obtaining pulsation times of miras ( @xcite ) or calculating isochrones of globular clusters(@xcite ) . an obvious limitation of this study is the use of 2d simulations ( and not 3d simulations ) . the spectrum of eddy sizes is probably different in 3d , but the influence of this upon the calibrated mixing length parameter is not known . however , @xcite have first estimation of the change in the calibrated mixing length parameter for their study of solar type convection , due to the difference between 2d and 3d simulations . they conclude that the calibrated mixing length parameter for the sun from 3d simulation is higher by 0.07 than from 2d simulation . i am grateful to yitzchak tuchman for his encouraging cooperation and to eli livne for his code . i thank ami glasner , yossi stein and zalman barkat for many fruitful discussions , and dave arnett and the referee for useful suggestions to the manuscript . abbett , w.p . , beaner , m. , davids , b. , georgobiani , d. , rathbun , p. , & stein r.f . 1997 , , 480 , 395 arnett , w.d . 1994 , , 427 , 932 asida , s.m . , & tuchman , y. 1997 , , 491 , l47 basu , s. & antia , h. m. 1994 , journal of astrophysics and astronomy , 15 , 143 bazan , g. , & arnett , w.d . 1994 , , 433 , l41 bazan , g. , & arnett , w.d . 1998 , , 496 , 316 bogdan , t.j . , cattaneo , f. , & malagoni , a. , 1993 , , 407 , 316 bohm - vitense e. , 1958 , z. astrophys . 46 , 108 buchler , j.r . , yecko , p. , koll@xmath51th , z. , & goupil , m .- j . 1999 , astro - ph/9905025 canuto , v.m . , 1996 , , 467 , 385 canuto , v.m . , 1997 , , 482 , 827 canuto , v.m . , & christensen - dalsgaard , j. , 1998 , annu fluid mech . , 30 , 167 canuto , v.m . , & dubovikov , m. , 1998 , , 493 , 834 canuto , v.m . , goldman , i. , & mazzitelli , i. , 1996 , , 473 , 550 canuto , v.m . , & mazzitelli , i. , 1991 , , 370 , 295 castaing , b. , gunaratne , g. , & heslot , f. , 1989 , j. fluid mech . , 204 , 1 cattaneo , f. , brummell , n.h . , toomre , j. , malagoni , a. , & hurlburt , n.e . 1991 , , 370 , 282 cattaneo , f. , hurlburt , n.e . , & toomre , j. 1990 , 349 , l63 chan , k.l . , & gigas , d. 1992 , , 389 , l87 chan , k.l . , & sofia , s. 1986 , , 307 , 222 chan , k.l . , sofia , s. , & wolff , c.l . 1982 , , 263 , 935 cox , a.n . , & stewart , j.n . 1970 , , 19 , 243 cox , j.p . , & giuli , r.t . 1968 , principles of stellar structure ( new york : gordon & breach ) demarque . p. , guenther , d.b . , & kim , y .- c . 1997 , , 474 , 790 deupree , r.g . 1975 , 201 , 183 deupree , r.g . 1977 , 211 , 509 deupree , r.g . 1985 , 296 , 160 deupree , r.g . 1998 , 499 , 340 freytag , b. , ludwig , h .- g . , & steffen , m. 1999 , in proceedings contribution to the workshop on stellar structure : `` theory and tests of convective energy transport '' granada , spain , 1998 , eds . a.  gimenez , e.  guinan , b.  montesinos freytag , b. & salaris , m. 1999 , , 513 , l49 gehmeyr , m. & winkler , k. -h . a. 1992 , , 253 , 92 glasner , s.a . , & livne , e. 1995 , , 445 , l149 glasner , s.a . , & livne , e. , & truran , j.w . 1997 , , 475 , 754 grossman , s. a. 1996 , , 279 , 305 grossman , s.a . , narayan , r. 1993 , , 89 , 361 grossman , s.a . , narayan , r. , & arnett , d. 1993 , , 407 , 284 helsot , f. , castaing , b. , & libchaber , a. 1987 , , 36 , 5870 hossain , m. , & mullan , d.j . 1991 , , 380 , 631 hurlburt , n. e. , toomre , j. , & massaguer , j. m. 1984 , , 282 , 557 hurlburt , n. e. , toomre , j. , & massaguer , j. m. 1986 , , 311 , 563 hurlburt , n. e. , toomre , j. , massaguer , j. m. , & zahn , j. -p . 1994 , , 421 , 245 kercek , a. , hillebrant , w. , & truran , j.w . 1998 , , 337 , 379 kercek , a. , hillebrant , w. , & truran , j.w . 1998 , , 345 , 831 kerr , r.m . 1996 , j. fluid mech . , 310 , 139 kim , y. -c . , fox , p. a. , sofia , s. , & demarque , p. 1995 , , 442 , 422 kim , y. -c . , fox , p. a. , demarque , p. , & sofia , s. 1996 , , 461 , 499 livne , e. 1993 , , 412 , 634 ludwig , h .- g . , freytag , b. , & steffen , m. 1999 , , 346 , 111 ludwig , h. -g . , jordan , s. & steffen , m. 1994 , , 284 , 105 maeder , a. , & conti , p.s . 1994 , , 32 , 227 malagoni , a. , cattaneo , f. , & brummell , n.h . 1990 , , 361 , l33 mezzacappa , a. , calder , a. c. , bruenn , s. w. , blondin , j. m. , guidry , m. w. , strayer , m. r. & umar , a. s. 1998 , , 495 , 911 m@xmath52ller , e. , & janka h .- 1997 , , 317 , 140 nordlund , a. , & stein , r.f . 1995 , stellar evolution : what can be done , 32@xmath53 liege int . nordlund , a. , stein , r.f . , & brandenburg , a. 1996 , bull astron . india , 24 paczinski , b. 1969 , acta astr . 19 , 1 porter , d.h . , & woodward , p.r . 1994 , , 93 , 309 rast , m.p . , & toomre , j. 1993 , , 419 , 224 schwarzchild , m. 1961 , , 134 , 1 singh , h. p. , & chan , k. l. 1993 , , 279 , 107 singh , h. p. , roxburgh , i. w. , & chan , k. l. 1994 , , 281 , l73 singh , h. p. , roxburgh , i. w. , & chan , k. l. 1995 , , 295 , 703 smagorinsky , j. 1963 , mon . weather rev . , 91(3 ) , 99 smalley , b. , & kupka , f. 1997 , , 328 , 349 spiegel , e.a . 1963 , , 138 , 216 spruit , h.c . , nordlund , a. , & title , a.m. 1990 , , 28 , 263 stein r.f . , & nordlund , a. 1989 , 342 , l95 stein r.f . , & nordlund , a. 1998 , 499 , 914 stellingwerf , r.f . 1982 , , 262 , 330 stone , j. m. , & balbus , s. a. 1996 , , 464 , 364 stothers , r. b. , & chin , c. -w . 1995 , , 440 , 297 stothers , r. b. , & chin , c. -w . 1997 , , 478 , l103 tligner , a. , belmonte , a. , & libchaber , a. 1993 , , 47 , r2253 tuchman , y. 1980 , `` dynamical phenomena in red giant envelopes '' , phd . thesis tuchman , y. , sack , n. , & barkat , z. 1978 , , 219 , 183 tuchman , y. , sack , n. , & barkat , z. 1979 , , 234 , 217 ulrich , r.k . 1970 , ap . space sci . , 7 , 71 ulrich , r.k . 1976 , , 207 , 564 van leer b. 1979 , j. comput . phys . , 32 , 101 werne , j. 1993 , , 48 , 1020 werne , j. 1994 , , 49 , 4072 wu , x .- z . , & libchaber , a. 1992 , , 45 , 842 xie , x. , & toomre , j. 1993 , , 405 , 747 xiong , d. r. , cheng , q. l. , & deng , l. 1997 , , 108 , 529 xiong , d. r. , deng , l. , & cheng , q. l. 1998 , , 499 , 355 xiong , d. r. , cheng , q. l. , & deng , l. 1998 , , 500 , 449 yaari , a. , & tuchman , y. 1998 , private communication zahn , j. -p . 1991 , , 252 , 179
Abstract:

Answer: two - dimensional hydrodynamical simulations were made to calibrate the mixing length parameter for modeling red giant s convective envelope . as was briefly reported in @xcite , a comparison of simulations starting with models integrated with different values of the mixing length parameter , has been made . in this paper more results are presented , including tests of the spatial resolution and large eddy simulation terms used by the numerical code . the consistent value of the mixing length parameter was found to be 1.4 , for a red giant of mass @xmath0 , core mass of @xmath1 , luminosity of @xmath2 , and metallicity @xmath3 .