======================================== SARAH 4.6.0: Changelog ======================================== --------------------------------------------------------------- 1) GENERATING TEMPLATES FOR PARTICLES.m and PARAMETERS.m --------------------------------------------------------------- After the implementation of a new model, one can now run WriteTemplatesParFiles; This will create templates for the particles.m and parameters.m file... 1) ... including all new particles and parameters 2) ... auto-generating LaTeX names, Output names and FeynArts numbers 3) ... extending/generating PDGs Note, some physical information has to be adjusted by the user: 1) Electric charge of new particles 2) Goldstone corresponding to massive vectors The files are written to the model directory. --------------------------------------------------------------- 2) AUXILIARY YUKAWAS IN SPHENO.m --------------------------------------------------------------- For models in which the SM Yukawa couplings and/or THDM VEVs don't exist, but other Yukawa couplings/VEVs are present which are functions of them, new options for SPheno.m exist: - DEFINITION[MoreEWvevs]: Can be used to relate new VEVs to the ew VEV - DEFINITION[UseNonStandardYukwas]: Can be used to turn on the new functionality to relate other Yukwas than the common Yu, Yd, Ye to SM parameters - DEFINITION[NonStandardYukawas]: names of the new Yukawas - DEFINITION[NonStandardVEVs1L]: Can be used to turn on the shift in the new VEVs at the one-loop level - DEFINITION[NonStandardVEVs1Lrelations]: defines the relations to shift the VEVs at one-loop - AuxiliaryParametersSPheno: contains the initialisation of parameters needed for the above definitions EXAMPLE: Left-Right SUSY Model with two Higgs bi-doublets: a) The VEVs of the two generations of bi-doublets are defined by DEFINITION[MoreEWvevs] = { "vPhiu(1) = (TanBeta*TanBetaU*sqrt(vev2))/(Sqrt(1 + TanBeta**2)*Sqrt(1 + TanBetaU**2))", "vPhid(1) = (TanBetaD*sqrt(vev2))/(Sqrt(1 + TanBeta**2)*Sqrt(1 + TanBetaD**2))", "vPhiu(2) = (TanBeta*sqrt(vev2))/(Sqrt(1 + TanBeta**2)*Sqrt(1 + TanBetaU**2))", "vPhid(2) = sqrt(vev2)/(Sqrt(1 + TanBeta**2)*Sqrt(1 + TanBetaD**2))" }; b) The LR symmetric Yukawas are calculated using DEFINITION[UseNonStandardYukwas]=True; DEFINITION[NonStandardYukawas] = {YQ,YL}; DEFINITION[NonStandardYukawasRelations] = { "YQ(1:3,1:3,1) = -(vu * vPhid(2) * Transpose(Y_u) - vd * vPhiu(2) * Transpose(Y_d))/(vPhid(2) * vPhiu(1) - vPhid(1) * vPhiu(2))", "YQ(1:3,1:3,2) = (vu * vPhid(1) * Transpose(Y_u) - vd * vPhiu(1) * Transpose(Y_d))/(vPhid(2) * vPhiu(1) - vPhid(1) * vPhiu(2))", "YL(1:3,1:3,1) = ( vd * vPhiu(2) * Transpose(Y_l))/(vPhid(2) * vPhiu(1) - vPhid(1) * vPhiu(2))", "YL(1:3,1:3,2) = (- vd * vPhiu(1) * Transpose(Y_l))/(vPhid(2) * vPhiu(1) - vPhid(1) * vPhiu(2))" } c) The shift in the VEVs at the one-loop level are set via DEFINITION[NonStandardVEVs1L]=True; DEFINITION[NonStandardVEVs1Lrelations] = { "tanbQ=sqrt((vPhiu(1)**2+vPhiu(2)**2)/(vPhid(1)**2+vPhid(2)**2))", "tanbUQ = vPhiu(1)/vPhiu(2)", "tanbDQ = vPhid(1)/vPhid(2)", "vev2 = (4*((mZ2+dmz2)*(-4*(mZ2+dmz2) + (gBL**2 + gR**2)*vR**2)))/(-4*(g2**2 + gR**2)*mZ2 + (gBL**2*gR**2 + g2**2*(gBL**2 + gR**2))*vR**2)", "vPhiu(1) = (TanBQ*TanBUQ*sqrt(vev2))/(Sqrt(1 + TanBQ**2)*Sqrt(1 + TanBUQ**2))", "vPhid(1) = (TanBDQ*sqrt(vev2))/(Sqrt(1 + TanBQ**2)*Sqrt(1 + TanBDQ**2))", "vPhiu(2) = (TanBQ*sqrt(vev2))/(Sqrt(1 + TanBQ**2)*Sqrt(1 + TanBUQ**2))", "vPhid(2) = sqrt(vev2)/(Sqrt(1 + TanBQ**2)*Sqrt(1 + TanBDQ**2))" }; d) The new parameters are define via AuxiliaryParametersSPheno={ "Real(dp) :: tanbQ, tanBUQ, tanBDQ" }; --------------------------------------------------------------- 3) CHANGES IN MODEL FILES --------------------------------------------------------------- - The SPheno output of all files contains an example point which is used for the template of the Les Houches input - New model for THDM with CPV - UMSSM: has been extended by right-handed neutrinos superfields - THDM models: have been extended by Z_2 breaking term m12 H1.conj[H2] - B-L-SM/B-L-SM-IS: fixed particles.m to include width of new Higgs states - Changes in SPheno.m of several models (TMSSM, TNMSSM, SSDM, DiracGauginos, MRSSM) --------------------------------------------------------------- 4) BUGFIXES/IMPROVEMENTS/OTHER CHANGES --------------------------------------------------------------- Bugfixes: - Missing factor 1/2 in kovariant derivative for scalars defined as real from the very beginning in non-SUSY models (affects only TSM as public models) [Thank you to A. Vicente] - Fixed error messages from Mathematica in check of charge conservation of non-SUSY models - Fixed streaming errors appearing with Mathematica 10.0.2 - Fixed R-symmetry violating decays in SPheno (affects only MRSSM as public models) [Thank you to W. ] - Fixed problems with tadpoles at one-loop in SPheno in models with CPV [Thank you to A. Voigt] - Fixed problem in SPheno that renormalisation scale of non-SUSY models is sometimes put to 1 GeV: new default scale is 160 GeV; for any other scale, use MODSEL[12] - Fixed problem of UFO output for models with real scalars which appear in one generation only - Fixed problem of polarisation with WHIZARD Improvements: - Some improvements in LaTeX output - Improved handling of phases for fermion which don't mix - For the impatient users: new short introduction 'SARAH in a nutshell' as pdf included Other changes: - SPheno: 3-body decays with vector propagators adjusted to new conventions of SPheno-3.3.8