HiggsBounds 2.0.0: Confronting neutral and charged Higgs sector predictions with exclusion bounds from LEP and the Tevatron

HiggsBounds 2.0.0 is a computer code which tests both neutral and charged Higgs sectors of arbitrary models against the current exclusion bounds from the Higgs searches at LEP and the Tevatron. As input, it requires a selection of model predictions, such as Higgs masses, branching ratios, effective couplings and total decay widths. HiggsBounds 2.0.0 then uses the expected and observed topological cross section limits from the Higgs searches to determine whether a given parameter scenario of a model is excluded at the 95% C.L. by those searches. Version 2.0.0 represents a significant extension of the code since its first release (1.0.0). It includes now 28/53 LEP/Tevatron Higgs search analyses, compared to the 11/22 in the first release, of which many of the ones from the Tevatron are replaced by updates. As a major extension, the code allows now the predictions for (singly) charged Higgs bosons to be confronted with LEP and Tevatron searches. Furthermore, the newly included analyses contain LEP searches for neutral Higgs bosons (H) decaying invisibly or into (non-flavour tagged) hadrons as well as decay-mode independent searches for neutral Higgs bosons, LEP searches via the production modes τ⁺τ⁻H and , and Tevatron searches via . Also, all Tevatron results presented at the ICHEP?10 are included in version 2.0.0. As physics applications of HiggsBounds 2.0.0 we study the allowed Higgs mass range for model scenarios with invisible Higgs decays and we obtain exclusion results for the scalar sector of the Randall–Sundrum model using up-to-date LEP and Tevatron direct search results.

Program summary

Program title: HiggsBoundsCatalogue identifier: AEFF_v2_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEFF_v2_0.htmlProgram obtainable from: CPC Program Library, Queen?s University, Belfast, N. IrelandLicensing provisions: GNU General Public Licence version 3No. of lines in distributed program, including test data, etc.: 74 005No. of bytes in distributed program, including test data, etc.: 1 730 996Distribution format: tar.gzProgramming language: Fortran 77, Fortran 90 (two code versions are offered).Classification: 11.1.Catalogue identifier of previous version: AEFF_v1_0Journal reference of previous version: Comput. Phys. Comm. 181 (2010) 138External routines: HiggsBounds requires no external routines/libraries. Some sample programs in the distribution require the programs FeynHiggs 2.7.1 or CPsuperH2.2 to be installed.Does the new version supersede the previous version?: YesNature of problem: Determine whether a parameter point of a given model is excluded or allowed by LEP and Tevatron neutral and charged Higgs boson search results.Solution method: The most sensitive channel from LEP and Tevatron searches is determined and subsequently applied to test this parameter point. The test requires as input, model predictions for the Higgs boson masses, branching ratios and ratios of production cross sections with respect to reference values.Reasons for new version: This version extends the functionality of the previous version.Summary of revisions: List of included Higgs searches has been expanded, e.g. inclusion of (singly) charged Higgs boson searches. The input required from the user has been extended accordingly.Restrictions: Assumes that the narrow width approximation is applicable in the model under consideration and that the model does not predict a significant change to the signature of the background processes or the kinematical distributions of the signal cross sections.Running time: About 0.01 seconds (or less) for one parameter point using one processor of an Intel Core 2 Quad Q6600 CPU at 2.40 GHz for sample model scenarios with three Higgs bosons. It depends on the complexity of the Higgs sector (e.g. the number of Higgs bosons and the number of open decay channels) and on the code version.     

SPICE: Simulation Package for Including Flavor in Collider Events

We describe SPICE: Simulation Package for Including Flavor in Collider Events. SPICE takes as input two ingredients: a standard flavor-conserving supersymmetric spectrum and a set of flavor-violating slepton mass parameters, both of which are specified at some high “mediation” scale. SPICE then combines these two ingredients to form a flavor-violating model, determines the resulting low-energy spectrum and branching ratios, and outputs HERWIG and SUSY Les Houches files, which may be used to generate collider events. The flavor-conserving model may be any of the standard supersymmetric models, including minimal supergravity, minimal gauge-mediated supersymmetry breaking, and anomaly-mediated supersymmetry breaking supplemented by a universal scalar mass. The flavor-violating contributions may be specified in a number of ways, from specifying charges of fields under horizontal symmetries to completely specifying all flavor-violating parameters. SPICE is fully documented and publicly available, and is intended to be a user-friendly aid in the study of flavor at the Large Hadron Collider and other future colliders.

Program summary

Program title: SPICECatalogue identifier: AEFL_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEFL_v1_0.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 8153No. of bytes in distributed program, including test data, etc.: 67 291Distribution format: tar.gzProgramming language: C++Computer: Personal computerOperating system: Tested on Scientific Linux 4.xClassification: 11.1External routines: SOFTSUSY [1,2] and SUSYHIT [3]Nature of problem: Simulation programs are required to compare theoretical models in particle physics with present and future data at particle colliders. SPICE determines the masses and decay branching ratios of supersymmetric particles in theories with lepton flavor violation. The inputs are the parameters of any of several standard flavor-conserving supersymmetric models, supplemented by flavor-violating parameters determined, for example, by horizontal flavor symmetries. The output are files that may be used for detailed simulation of supersymmetric events at particle colliders.Solution method: Simpson's rule integrator, basic algebraic computation.Additional comments: SPICE interfaces with SOFTSUSY and SUSYHIT to produce the low energy sparticle spectrum. Flavor mixing for sleptons and sneutrinos is fully implemented; flavor mixing for squarks is not included.Running time: <1 minute. Running time is dominated by calculating the possible and relevant three-body flavor-violating decays of sleptons, which is usually 10-15 seconds per slepton.References:

[1]

B.C. Allanach, Comput. Phys. Commun. 143 (2002) 305, arXiv:hep-ph/0104145.

[2]

B.C. Allanach, M.A. Bernhardt, arXiv:0903.1805 [hep-ph].

[3]

A. Djouadi, M.M. Muhlleitner, M. Spira, Acta Phys. Pol. B 38 (2007) 635, arXiv:hep-ph/0609292.

     

FLY. A parallel tree N-body code for cosmological simulations

FLY is a parallel treecode which makes heavy use of the one-sided communication paradigm to handle the management of the tree structure. In its public version the code implements the equations for cosmological evolution, and can be run for different cosmological models.This reference guide describes the actual implementation of the algorithms of the public version of FLY, and suggests how to modify them to implement other types of equations (for instance, the Newtonian ones).

Program summary

Title of program:FLYCatalogue identifier: ADSCProgram summary URL:http://cpc.cs.qub.ac.uk/summaries/ADSCProgram obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandComputer for which the program is designed and others on which it has been tested: Cray T3E, Sgi Origin 3000, IBM SPOperating systems or monitors under which the program has been tested: Unicos 2.0.5.40, Irix 6.5.14, Aix 4.3.3Programming language used: Fortran 90, CMemory required to execute with typical data: about 100 Mwords with 2 million-particlesNumber of bits in a word: 32Number of processors used: parallel program. The user can select the number of processors ?1Has the code been vectorized or parallelized?: parallelizedNumber of bytes in distributed program, including test data, etc.: 4 615 604Distribution format: tar gzip fileKeywords: Parallel tree N-body code for cosmological simulationsNature of physical problem:FLY is a parallel collisionless N-body code for the calculation of the gravitational force.Method of solution: It is based on the hierarchical oct-tree domain decomposition introduced by Barnes and Hut (1986).Restrictions on the complexity of the program: The program uses the leapfrog integrator schema, but could be changed by the user.Typical running time: 50 seconds for each time-step, running a 2-million-particles simulation on an Sgi Origin 3800 system with 8 processors having 512 Mbytes RAM for each processor.Unusual features of the program:FLY uses the one-side communications libraries: the SHMEM library on the Cray T3E system and Sgi Origin system, and the LAPI library on IBM SP system     

Vscape V1.1.0. An interactive tool for metastable vacua

Vscape is an interactive tool for studying the one-loop effective potential of an ungauged supersymmetric model of chiral multiplets. The program allows the user to define a supersymmetric model by specifying the superpotential. The F-terms and the scalar and fermionic mass matrices are calculated symbolically. The program then allows you to search numerically for (meta)stable minima of the one-loop effective potential. Additional commands enable you to further study specific minima, by, e.g., computing the mass spectrum for those vacua. Vscape combines the flexibility of symbolic software, with the speed of a numerical package.

Program summary

Program title:Vscape 1.1.1Catalogue identifier: ADZW_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADZW_v1_0.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 80 507No. of bytes in distributed program, including test data, etc.: 6 708 938Distribution format: tar.gzProgramming language: C++Computer: Pentium 4 PC Computers: need (GNU) C++ compiler, Linux standard GNU installation (./configure; make; make install). A precompiled Windows XP version is included in the distribution packageOperating system: Linux, Windows XP using cygwinRAM: 10 MBWord size: 32 bitsClassification: 11.6External routines: GSL (http://www.gnu.org/software/gsl/), CLN (http://www.ginac.de/CLN/), GiNaC (http://directory.fsf.org/GiNaC.html)Nature of problem:Vscape is an interactive tool for studying the one-loop effective potential of an ungauged supersymmetric model of chiral multiplets. The program allows the user to define a supersymmetric model by specifying the superpotential. The F-terms and the scalar and fermionic mass matrices are calculated symbolically. The program then allows you to search numerically for (meta)stable minima of the one-loop effective potential. Additional commands enable you to further study specific minima, by, e.g., computing the mass spectrum for those vacua. Vscape combines the flexibility of symbolic software with the speed of a numerical package.Solution method: Coleman-Weinberg potential is computed using numerical matrix diagonalization. Minima of the one-loop effective potential are found using the Nelder and Mead simplex algorithm. The one-loop effective potential can be studied using numerical differentiation. Symbolic users interface implemented using flex and bison.Restrictions:N=1 supersymmetric chiral models onlyUnusual features: GiNaC (+CLN), GSL, ReadLib (not essential)Running time: Interactive users interface. Most commands execute in a few ms. Computationally intensive commands execute in order of minutes, depending on the complexity of the user defined model.     

SuperIso v3.0, flavor physics observables calculations: extension to NMSSM

SuperIso v3.0 is a public program for evaluation of flavor physics observables in the minimal supersymmetric extension of the Standard Model (MSSM) and the next to minimal supersymmetric extension of the Standard Model (NMSSM). SuperIso v3.0 incorporates many flavor observables such as the inclusive branching ratio of B→Xsγ, the isospin asymmetry of B→K^∗γ, the branching ratio of Bs→μ⁺μ⁻, the branching ratio of B→τντ, the branching ratio of B→Dτντ, the branching ratio of K→μνμ and the branching ratios of Ds→τντ and Ds→μνμ. The calculation of the branching ratio of B→Xsγ includes NNLO Standard Model contributions. The program also computes the muon anomalous magnetic moment (g−2). Seven sample models are included in the package, namely mSUGRA, NUHM, AMSB and GMSB for the MSSM, and CNMSSM, NGMSB and NNUHM for the NMSSM. SuperIso uses a SUSY Les Houches Accord file (SLHA1 or SLHA2) as input, which can be either generated automatically by the program via a call to external spectrum calculators (SOFTSUSY, ISAJET or NMSSMTools), or provided by the user.

New version program summary

Program title:SuperIso v3.0Catalogue identifier: AEAN_v3_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEAN_v3_0.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: GNU General Public LicenceNo. of lines in distributed program, including test data, etc.: 6869No. of bytes in distributed program, including test data, etc.: 42 627Distribution format: tar.gzProgramming language: C (C99 Standard compliant)Computer: 32- or 64-bit PC, MacOperating system: Linux, MacOSRAM: less than 1 MBClassification: 11.6External routines: ISASUGRA/ISAJET, SOFTSUSY and/or NMSSMToolsDoes the new version supersede the previous version?: YesNature of problem: Calculation of flavor physics observables as well as the muon anomalous magnetic moment in the Minimal Supersymmetric Standard Model with minimal flavor violation and in the Next to Minimal Supersymmetric Standard Model, in order to derive constraints on the supersymmetric parameter spaces.Solution method:SuperIso uses a SUSY Les Houches Accord (SLHA1 or SLHA2) file, which can be either generated automatically via a call to SOFTSUSY, ISAJET or NMSSMTools, or provided by the user. This file contains the masses, mixings and couplings of the supersymmetric particles. SuperIso then computes the most constraining flavor physics observables and the muon (g−2). SuperIso is able to perform the calculations in different supersymmetry breaking scenarios, such as mSUGRA, NUHM, AMSB and GMSB, as well as constrained NMSSM scenarios such as CNMSSM, NNUHM and NGMSB.Reasons for new version:SuperIso has been extended to the next to minimal supersymmetric extension of the Standard Model (NMSSM). The implemented routines are therefore extensively modified.Summary of revisions:

•

Improvement of the SLHA2 reader.

•

Replacement of “float” variables by “double”.

•

Implementation of an interface with NMSSMTools.

•

Extension of the calculation of flavor observables as well as the muon anomalous magnetic moment to NMSSM.

•

Addition of three different NMSSM scenarios: CNMSSM, NGMSB and NNUHM.

•

Three sample main programs have been added: cnmssm.c, ngmsb.c and nnuhm.c. Additional instructions to use them are given when running them without arguments.

Unusual features: The code is very flexible, and new observables can be added easily.Running time: Less than 1 sec     

Including R-parity violation in the numerical computation of the spectrum of the minimal supersymmetric standard model: SOFTSUSY3.0

Current publicly available computer programs calculate the spectrum and couplings of the minimal supersymmetric standard model under the assumption of R-parity conservation. Here, we describe an extension to the SOFTSUSY program which includes R-parity violating effects. The user provides a theoretical boundary condition upon the high-scale supersymmetry breaking R-parity violating couplings. Successful radiative electroweak symmetry breaking, electroweak and CKM matrix data are used as weak-scale boundary conditions. The renormalisation group equations are solved numerically between the weak scale and a high energy scale using a nested iterative algorithm. This paper serves as a manual to the R-parity violating mode of the program, detailing the approximations and conventions used.

Program summary

Program title:SOFTSUSY v3.0Catalogue identifier: ADPM_v2_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADPM_v2_0.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 75 927No. of bytes in distributed program, including test data, etc.: 570 916Distribution format: tar.gzProgramming language: C++, FortranComputer: Personal computerOperating system: Tested on Linux 4.xWord size: 32 bitsClassification: 11.6Catalogue identifier of previous version: ADPM_v1_0Journal reference of previous version: Comput. Phys. Comm. 143 (2002) 305Does the new version supersede the previous version?: YesNature of problem: Calculating supersymmetric particle spectrum and mixing parameters in the R-parity violating minimal supersymmetric standard model. The solution to the renormalisation group equations must be consistent with a high-scale boundary condition on supersymmetry breaking parameters and Rp parameters, as well as a weak-scale boundary condition on gauge couplings, Yukawa couplings and the Higgs potential parameters.Solution method: Nested iterative algorithmReasons for new version: This is an extension to the SOFTSUSY program which includes R-parity violating effects. The user provides a theoretical boundary condition upon the high-scale supersymmetry breaking R-parity violating couplings. Successful radiative electroweak symmetry breaking, electroweak and CKM matrix data are used as weak-scale boundary conditions. The renormalisation group equations are solved numerically between the weak scale and a high energy scale using a nested iterative algorithm. The paper serves as a manual to the R-parity violating mode of the program, detailing the approximations and conventions used.Restrictions:SOFTSUSY3.0 will provide a solution only in the perturbative regime and it assumes that all couplings of the MSSM are real (i.e. CP-conserving). The iterative SOFTSUSY algorithm will not converge if parameters are too close to a boundary of successful electroweak symmetry breaking, but a warning flag will alert the user to this fact.Running time: A few seconds per parameter point.     

SuSpect: A Fortran code for the Supersymmetric and Higgs particle spectrum in the MSSM

We present the Fortran code SuSpect version 2.3, which calculates the Supersymmetric and Higgs particle spectrum in the Minimal Supersymmetric Standard Model (MSSM). The calculation can be performed in constrained models with universal boundary conditions at high scales such as the gravity (mSUGRA), anomaly (AMSB) or gauge (GMSB) mediated supersymmetry breaking models, but also in the non-universal MSSM case with R-parity and CP conservation. Care has been taken to treat important features such as the renormalization group evolution of parameters between low and high energy scales, the consistent implementation of radiative electroweak symmetry breaking and the calculation of the physical masses of the Higgs bosons and supersymmetric particles taking into account the dominant radiative corrections. Some checks of important theoretical and experimental features, such as the absence of non-desired minima, large fine-tuning in the electroweak symmetry breaking condition, as well as agreement with precision measurements can be performed. The program is simple to use, self-contained and can easily be linked to other codes; it is rather fast and flexible, thus allowing scans of the parameter space with several possible options and choices for model assumptions and approximations.

Program summary

Title of program:SuSpectCatalogue identifier:ADYR_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADYR_v1_0Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandLicensing provisions:noneProgramming language used:FORTRAN 77Computer:Unix machines, PCNo. of lines in distributed program, including test data, etc.:21 821No. of bytes in distributed program, including test data, etc.:249 657Distribution format:tar.gzOperating system:Unix (or Linux)RAM:approximately 2500 KbytesNumber of processors used:1 processorNature of problem:SuSpect calculates the supersymmetric and Higgs particle spectrum (masses and some other relevant parameters) in the unconstrained Minimal Supersymmetric Standard Model (MSSM), as well as in constrained models (cMSSMs) such as the minimal Supergravity (mSUGRA), the gauge mediated (GMSB) and anomaly mediated (AMSB) Supersymmetry breaking scenarii. The following features and ingredients are included: renormalization group evolution between low and high energy scales, consistent implementation of radiative electroweak symmetry breaking, calculation of the physical particle masses with radiative corrections at the one- and two-loop level.Solution method:The main methods used in the code are: (1) an (adaptative fourth-order) Runge-Kutta type algorithm (following a standard algorithm described in “Numerical Recipes”), used to solve numerically a set of coupled differential equations resulting from the renormalization group equations at the two-loop level of the perturbative expansions; (2) diagonalizations of mass matrices; (3) some mathematical (Spence, etc) functions resulting from the evaluation of one and two-loop integrals using the Feynman graphs techniques for radiative corrections to the particle masses; (4) finally, some fixed-point iterative algorithms to solve non-linear equations for some of the relevant output parameters.Restrictions:(1) The code is limited at the moment to real input parameters. (2) It also does not include flavor non-diagonal terms which are possible in the most general soft supersymmetry breaking Lagrangian. (3) There are some (mild) limitations on the possible range of values of input parameter, i.e. not any arbitrary values of some input parameters are allowed: these limitations are essentially based on physical rather than algorithmic issues, and warning flags and other protections are installed to avoid as much as possible execution failure if unappropriate input values are used.Running time:between 1 and 3 seconds depending on options, with a 1 GHz processor.     

SPheno, a program for calculating supersymmetric spectra, SUSY particle decays and SUSY particle production at ee colliders

SPheno is a program that accurately calculates the supersymmetric particle spectrum within a high scale theory, such as minimal supergravity, gauge mediated supersymmetry breaking, anomaly mediated supersymmetry breaking, or string effective field theories. An interface exists for an easy implementation of other models. The program solves the renormalization group equations numerically to two-loop order with user-specified boundary conditions. The complete one-loop formulas for the masses are used which are supplemented by two-loop contributions in case of the neutral Higgs bosons and the μ parameter. The obtained masses and mixing matrices are used to calculate decay widths and branching ratios of supersymmetric particles as well as of Higgs bosons, b→sγ, Δρ and (g−2)_μ. Moreover, the production cross sections of all supersymmetric particle as well as Higgs bosons at e⁺e⁻ colliders can be calculated including initial state radiation and longitudinal polarization of the incoming electrons/positrons. The program is structured such that it can easily be extend to include non-minimal models and/or complex parameters.     

An updated version of wannier90: A tool for obtaining maximally-localised Wannier functions

wannier90  is a program for calculating maximally-localised Wannier functions (MLWFs) from a set of Bloch energy bands that may or may not be attached to or mixed with other bands. The formalism works by minimising the total spread of the MLWFs in real space. This is done in the space of unitary matrices that describe rotations of the Bloch bands at each k-point. As a result, wannier90  is independent of the basis set used in the underlying calculation to obtain the Bloch states. Therefore, it may be interfaced straightforwardly to any electronic structure code. The locality of MLWFs can be exploited to compute band-structure, density of states and Fermi surfaces at modest computational cost. Furthermore, wannier90  is able to output MLWFs for visualisation and other post-processing purposes. Wannier functions are already used in a wide variety of applications. These include analysis of chemical bonding in real space; calculation of dielectric properties via the modern theory of polarisation; and as an accurate and minimal basis set in the construction of model Hamiltonians for large-scale systems, in linear-scaling quantum Monte Carlo calculations, and for efficient computation of material properties, such as the anomalous Hall coefficient. We present here an updated version of wannier90, wannier90  2.0, including minor bug fixes and parallel (MPI) execution for band-structure interpolation and the calculation of properties such as density of states, Berry curvature and orbital magnetisation. wannier90  is freely available under the GNU General Public License from http://www.wannier.org/.     

micrOMEGAs 2.0: A program to calculate the relic density of dark matter in a generic model

micrOMEGAs 2.0 is a code which calculates the relic density of a stable massive particle in an arbitrary model. The underlying assumption is that there is a conservation law like R-parity in supersymmetry which guarantees the stability of the lightest odd particle. The new physics model must be incorporated in the notation of CalcHEP, a package for the automatic generation of squared matrix elements. Once this is done, all annihilation and coannihilation channels are included automatically in any model. Cross-sections at v=0, relevant for indirect detection of dark matter, are also computed automatically. The package includes three sample models: the minimal supersymmetric standard model (MSSM), the MSSM with complex phases and the NMSSM. Extension to other models, including non-supersymmetric models, is described.

Program summary

Title of program:micrOMEGAs2.0Catalogue identifier:ADQR_v2_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADQR_v2_0Program obtainable from: CPC Program Library, Queen's University of Belfast, N. IrelandComputers for which the program is designed and others on which it has been tested:PC, Alpha, Mac, SunOperating systems under which the program has been tested:UNIX (Linux, OSF1, SunOS, Darwin, Cygwin)Programming language used:C and FortranMemory required to execute with typical data:17 MB depending on the number of processes requiredNo. of processors used:1Has the code been vectorized or parallelized:noNo. of lines in distributed program, including test data, etc.:91 778No. of bytes in distributed program, including test data, etc.:1 306 726Distribution format:tar.gzExternal routines/libraries used:noCatalogue identifier of previous version:ADQR_v1_3Journal reference of previous version:Comput. Phys. Comm. 174 (2006) 577Does the new version supersede the previous version:yesNature of physical problem:Calculation of the relic density of the lightest stable particle in a generic new model of particle physics.Method of solution: In numerically solving the evolution equation for the density of dark matter, relativistic formulae for the thermal average are used. All tree-level processes for annihilation and coannihilation of new particles in the model are included. The cross-sections for all processes are calculated exactly with CalcHEP after definition of a model file. Higher-order QCD corrections to Higgs couplings to quark pairs are included.Reasons for the new version:There are many models of new physics that propose a candidate for dark matter besides the much studied minimal supersymmetric standard model. This new version not only incorporates extensions of the MSSM, such as the MSSM with complex phases, or the NMSSM which contains an extra singlet superfield but also gives the possibility for the user to incorporate easily a new model. For this the user only needs to redefine appropriately a new model file.Summary of revisions:

•

Possibility to include in the package any particle physics model with a discrete symmetry that guarantees the stability of the cold dark matter candidate (LOP) and to compute the relic density of CDM.

•

Compute automatically the cross-sections for annihilation of the LOP at small velocities into SM final states and provide the energy spectra for final states.

•

For the MSSM with input parameters defined at the GUT scale, the interface with any of the spectrum calculator codes reads an input file in the SUSY Les Houches Accord format (SLHA).

•

Implementation of the MSSM with complex parameters (CPV-MSSM) with an interface to CPsuperH to calculate the spectrum.

•

Routine to calculate the electric dipole moment of the electron in the CPV-MSSM.

•

In the NMSSM, new interface compatible with NMHDECAY2.1.

Typical running time:0.2 secUnusual features of the program:Depending on the parameters of the model, the program generates additional new code, compiles it and loads it dynamically.     

Automatic calculation of supersymmetric renormalization group equations and loop corrections

SARAH is a Mathematica package for studying supersymmetric models. It calculates for a given model the masses, tadpole equations and all vertices at tree-level. This information can be used by SARAH to write model files for CalcHep/CompHep or FeynArts/FormCalc. In addition, the second version of SARAH can derive the renormalization group equations for the gauge couplings, parameters of the superpotential and soft-breaking parameters at one- and two-loop level. Furthermore, it calculates the one-loop self-energies and the one-loop corrections to the tadpoles. SARAH can handle all N=1 SUSY models whose gauge sector is a direct product of SU(N) and U(1) gauge groups. The particle content of the model can be an arbitrary number of chiral superfields transforming as any irreducible representation with respect to the gauge groups. To implement a new model, the user has just to define the gauge sector, the particle, the superpotential and the field rotations to mass eigenstates.

Program summary

Program title: SARAHCatalogue identifier: AEIB_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEIB_v1_0.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.htmlNo. of lines in distributed program, including test data, etc.: 97 577No. of bytes in distributed program, including test data, etc.: 2 009 769Distribution format: tar.gzProgramming language: MathematicaComputer: All systems that Mathematica is available forOperating system: All systems that Mathematica is available forClassification: 11.1, 11.6Nature of problem: A supersymmetric model is usually characterized by the particle content, the gauge sector and the superpotential. It is a time consuming process to obtain all necessary information for phenomenological studies from these basic ingredients.Solution method: SARAH calculates the complete Lagrangian for a given model whose gauge sector can be any direct product of SU(N) gauge groups. The chiral superfields can transform as any, irreducible representation with respect to these gauge groups and it is possible to handle an arbitrary number of symmetry breakings or particle rotations. Also the gauge fixing terms can be specified. Using this information, SARAH derives the mass matrices and Feynman rules at tree-level and generates model files for CalcHep/CompHep and FeynArts/FormCalc. In addition, it can calculate the renormalization group equations at one- and two-loop level and the one-loop corrections to the one- and two-point functions.Unusual features: SARAH just needs the superpotential and gauge sector as input and not the complete Lagrangian. Therefore, the complete implementation of new models is done in some minutes.Running time: Measured CPU time for the evaluation of the MSSM on an Intel Q8200 with 2.33 GHz. Calculating the complete Lagrangian: 12 seconds. Calculating all vertices: 75 seconds. Calculating the one- and two-loop RGEs: 50 seconds. Calculating the one-loop corrections: 7 seconds. Writing a FeynArts file: 1 second. Writing a CalcHep/CompHep file: 6 seconds. Writing the LaTeX output: 1 second.     

2HDMC - two-Higgs-doublet model calculator

We describe the public C++ code 2HDMC which can be used to perform calculations in a general, CP-conserving, two-Higgs-doublet model (2HDM). The program features simple conversion between different parametrizations of the 2HDM potential, a flexible Yukawa sector specification with choices of different Z₂-symmetries or more general couplings, a decay library including all two-body - and some three-body - decay modes for the Higgs bosons, and the possibility to calculate observables of interest for constraining the 2HDM parameter space, as well as theoretical constraints from positivity and unitarity. The latest version of the 2HDMC code and full documentation is available from: http://www.isv.uu.se/thep/MC/2HDMC.

Program summary

Program title:2HDMCCatalogue identifier: AEFI_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEFI_v1_0.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: GNU GPLNo. of lines in distributed program, including test data, etc.: 12 032No. of bytes in distributed program, including test data, etc.: 90 699Distribution format: tar.gzProgramming language: C++Computer: Any computer running LinuxOperating system: LinuxRAM: 5 MbClassification: 11.1External routines: GNU Scientific Library (http://www.gnu.org/software/gsl/)Nature of problem: Determining properties of the potential, calculation of mass spectrum, couplings, decay widths, oblique parameters, muon g−2, and collider constraints in a general two-Higgs-doublet model.Solution method: From arbitrary potential and Yukawa sector, tree-level relations are used to determine Higgs masses and couplings. Decay widths are calculated at leading order, including FCNC decays when applicable. Decays to off-shell vector bosons are obtained by numerical integration. Observables are computed (analytically or numerically) as function of the input parameters.Restrictions: CP-violation is not treated.Running time: Less than 0.1 s on a standard PC     

Standard SANC modules

In this note we summarize the status of the standard SANC modules (in the EW and QCD sectors of the Neutral Current branch — version 1.21 and the Charged Current branch — version 1.20). All versions of the codes are accessible from the SANC project servers at Dubna http://sanc.jinr.ru and CERN http://pcphsanc.cern.ch.     

Cobra: A package for co-breaking analysis

Cobra, an econometric software package designed for CO-BReaking Analysis, is introduced. Cobra is programmed in Ox, the object-oriented statistical system, and consists of three modules: Firstly, CobraDgp is derived from the Database class and enables the user to generate a multivariate time series that is subject to multiple breaks in its intercept as well as linear trend. Secondly, the Cobra module is derived from the Modelbase class and implements two algorithms for the estimation of co-breaking relationships. Taking advantage of the capabilities provided by Modelbase, Cobra may be loaded into OxPack and used with the GUI GiveWin as front end. Finally, CobraSim is derived from the Simulation class and wrapped around CobraDgp and Cobra to allow a straightforward implementation of Monte Carlo experiments. A brief introduction of the concept of co-breaking is given before Cobra is illustrated by means of an empirical example as well as a simple Monte Carlo. Excerpts of Ox code as well as screenshots are provided.     

CPMC-Lab: A Matlab package for Constrained Path Monte Carlo calculations

We describe CPMC-Lab, a Matlab program for the constrained-path and phaseless auxiliary-field Monte Carlo methods. These methods have allowed applications ranging from the study of strongly correlated models, such as the Hubbard model, to ab initio calculations in molecules and solids. The present package implements the full ground-state constrained-path Monte Carlo (CPMC) method in Matlab with a graphical interface, using the Hubbard model as an example. The package can perform calculations in finite supercells in any dimensions, under periodic or twist boundary conditions. Importance sampling and all other algorithmic details of a total energy calculation are included and illustrated. This open-source tool allows users to experiment with various model and run parameters and visualize the results. It provides a direct and interactive environment to learn the method and study the code with minimal overhead for setup. Furthermore, the package can be easily generalized for auxiliary-field quantum Monte Carlo (AFQMC) calculations in many other models for correlated electron systems, and can serve as a template for developing a production code for AFQMC total energy calculations in real materials. Several illustrative studies are carried out in one- and two-dimensional lattices on total energy, kinetic energy, potential energy, and charge- and spin-gaps.     

GRACE/SUSY: Automatic generation of tree amplitudes in the minimal supersymmetric standard model

GRACE/SUSY is a program package for generating the tree-level amplitude and evaluating the corresponding cross section of processes of the minimal supersymmetric extension of the standard model (MSSM). The Higgs potential adopted in the system, however, is assumed to have a more general form indicated by the two-Higgs-doublet model. This system is an extension of GRACE for the standard model (SM) of the electroweak and strong interactions. For a given MSSM process the Feynman graphs and amplitudes at tree-level are automatically created. The Monte Carlo phase space integration by means of BASES gives the total and differential cross sections. When combined with SPRING, an event generator, the program package provides us with the simulation of the SUSY particle productions.