Abstract: In this work, we explore the flavor-changing decays H-i -> bs in a general supersymmetric scenario. In these models the flavor-changing decays arise at loop level, but-because they originate from a dimension-four operator-they do not decouple and may provide a first sign of new physics for heavy masses beyond the reach of colliders. In the framework of the minimal supersymmetric extension of the Standard Model, we find that the largest branching ratio of the lightest Higgs (H-1) is O(10(-6)) after imposing present experimental constraints, while heavy Higgs states may still present branching ratios O(10(-3)). In a more general supersymmetric scenario, where additional Higgs states may modify the Higgs mixings, the branching ratio BR(H-1 -> bs) can reach values O(10(-4)), while heavy Higgses still remain at O(10(-3)). Although these values are clearly out of reach for the LHC, a full study in a linear collider environment could be worth pursuing.

Abstract: We introduce FlexibleSUSY, a Mathematica and C++ package, which generates a fast, precise C++ spectrum generator for any SUSY model specified by the user. The generated code is designed with both speed and modularity in mind, making it easy to adapt and extend with new features. The model is specified by supplying the superpotential, gauge structure and particle content in a SARAH model file; specific boundary conditions e.g. at the GUT, weak or intermediate scales are defined in a separate FlexibleSUSY model file. From these model files, FlexibleSUSY generates C++ code for self-energies, tadpole corrections, renormalization group equations (RGEs) and electroweak symmetry breaking (EWSB) conditions and combines them with numerical routines for solving the RGEs and EWSB conditions simultaneously. The resulting spectrum generator is then able to solve for the spectrum of the model, including loop-corrected pole masses, consistent with user specified boundary conditions. The modular structure of the generated code allows for individual components to be replaced with an alternative if available. FlexibleSUSY has been carefully designed to grow as alternative solvers and calculators are added. Predefined models include the MSSM, NMSSM, E6SSM, USSM, R-symmetric models and models with right-handed neutrinos. Program Summary Program title: FlexibleSUSY Catalogue identifier: AEVIv10 Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEVIv10.html obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: GNU General Public License, version 3 No. of lines in distributed program, including test data, etc.: 129406 No. of bytes in distributed program, including test data, etc.: 854831 Distribution format: tar.gz Programming language: C++, Wolfram/Mathematica, FORTRAN, Bourne shell. Computer: Personal computer. Operating system: Tested on Linux 3.x, Mac OS X. Classification: 11.1, 11.6, 6.5. External routines: SARAH 4.0.4, Boost library, Eigen, LAPACK Nature of problem: Determining the mass spectrum and mixings for any supersymmetric model. The generated code must find simultaneous solutions to constraints which are specified at two or more different renormalization scales, which are connected by renormalization group equations forming a large set of coupled first-order differential equations. Solution method: Nested iterative algorithm and numerical minimization of the Higgs potential. Restrictions: The couplings must remain perturbative at all scales between the highest and the lowest boundary condition. FlexibleSUSY assumes that all couplings of the model are real (i.e. CP-conserving). Due to the modular nature of the generated code, adaption and extension to overcome restrictions in scope is quite straightforward. Running time: 0.06-0.2 seconds per parameter point.