HiggsBounds: Confronting arbitrary Higgs sectors with exclusion bounds from LEP and the Tevatron

HiggsBounds is a computer code that tests theoretical predictions of models with arbitrary Higgs sectors against the exclusion bounds obtained from the Higgs searches at LEP and the Tevatron. The included experimental information comprises exclusion bounds at 95% C.L. on topological cross sections. In order to determine which search topology has the highest exclusion power, the program also includes, for each topology, information from the experiments on the expected exclusion bound, which would have been observed in case of a pure background distribution. Using the predictions of the desired model provided by the user as input, HiggsBounds determines the most sensitive channel and tests whether the considered parameter point is excluded at the 95% C.L. HiggsBounds is available as a Fortran 77 and Fortran 90 code. The code can be invoked as a command line version, a subroutine version and an online version. Examples of exclusion bounds obtained with HiggsBounds are discussed for the Standard Model, for a model with a fourth generation of quarks and leptons and for the Minimal Supersymmetric Standard Model with and without CP-violation. The experimental information on the exclusion bounds currently implemented in HiggsBounds will be updated as new results from the Higgs searches become available.

Program summary

Program title: HiggsBoundsCatalogue identifier: AEFF_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEFF_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.: 55 733No. of bytes in distributed program, including test data, etc.: 1 986 213Distribution format: tar.gzProgramming language: Fortran 77, Fortran 90 (two code versions are offered).Computer: HiggsBounds can be built with any compatible Fortran 77 or Fortran 90 compiler. The program has been tested on x86 CPUs running under Linux (Ubuntu 8.04) and with the following compilers: The Portland Group Inc. Fortran compilers (pgf77, pgf90), the GNU project Fortran compilers (g77, gfortran).Operating system: LinuxRAM: minimum of about 6000 kbytes (dependent on the code version)Classification: 11.1External routines: HiggsBounds requires no external routines/libraries. Some sample programs in the distribution require the programs FeynHiggs 2.6.x or CPsuperH2 to be installed (see “Subprograms used”).Subprograms used:

Cat Id	Title	Reference
ADKT_v2_0	FeynHiggsv2.6.5	CPC 180(2009)1426
ADSR_v2_0	CPsuperH2.0	CPC 180(2009)312

Full-size table

     

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     

FIESTA 2: Parallelizeable multiloop numerical calculations

The program FIESTA has been completely rewritten. Now it can be used not only as a tool to evaluate Feynman integrals numerically, but also to expand Feynman integrals automatically in limits of momenta and masses with the use of sector decompositions and Mellin–Barnes representations. Other important improvements to the code are complete parallelization (even to multiple computers), high-precision arithmetics (allowing to calculate integrals which were undoable before), new integrators, Speer sectors as a strategy, the possibility to evaluate more general parametric integrals.

Program summary

Program title:FIESTA 2Catalogue identifier: AECP_v2_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AECP_v2_0.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: GNU GPL version 2No. of lines in distributed program, including test data, etc.: 39 783No. of bytes in distributed program, including test data, etc.: 6 154 515Distribution format: tar.gzProgramming language: Wolfram Mathematica 6.0 (or higher) and CComputer: From a desktop PC to a supercomputerOperating system: Unix, Linux, Windows, Mac OS XHas the code been vectorised or parallelized?: Yes, the code has been parallelized for use on multi-kernel computers as well as clusters via Mathlink over the TCP/IP protocol. The program can work successfully with a single processor, however, it is ready to work in a parallel environment and the use of multi-kernel processor and multi-processor computers significantly speeds up the calculation; on clusters the calculation speed can be improved even further.RAM: Depends on the complexity of the problemClassification: 4.4, 4.12, 5, 6.5Catalogue identifier of previous version: AECP_v1_0Journal reference of previous version: Comput. Phys. Comm. 180 (2009) 735External routines: QLink [1], Cuba library [2], MPFR [3]Does the new version supersede the previous version?: YesNature of problem: The sector decomposition approach to evaluating Feynman integrals falls apart into the sector decomposition itself, where one has to minimize the number of sectors; the pole resolution and epsilon expansion; and the numerical integration of the resulting expression.Solution method: The sector decomposition is based on a new strategy as well as on classical strategies such as Speer sectors. The sector decomposition, pole resolution and epsilon-expansion are performed in Wolfram Mathematica 6.0 or, preferably, 7.0 (enabling parallelization) [4]. The data is stored on hard disk via a special program, QLink [1]. The expression for integration is passed to the C-part of the code, that parses the string and performs the integration by one of the algorithms in the Cuba library package [2]. This part of the evaluation is perfectly parallelized on multi-kernel computers.Reasons for new version:

• 1. 

  The first version of FIESTA had problems related to numerical instability, so for some classes of integrals it could not produce a result.

• 2. 

  The sector decomposition method can be applied not only for integral calculation.

Summary of revisions:

• 1. 

  New integrator library is used.

• 2. 

  New methods to deal with numerical instability (MPFR library).

• 3. 

  Parallelization in Mathematica.

• 4. 

  Parallelization on multiple computers via TCP-IP.

• 5. 

  New sector decomposition strategy (Speer sectors).

• 6. 

  Possibility of using FIESTA to for integral expansion.

• 7. 

  Possibility of using FIESTA to discover poles in d.

• 8. 

  New negative terms resolution strategies.

Restrictions: The complexity of the problem is mostly restricted by CPU time required to perform the evaluation of the integralRunning time: Depends on the complexity of the problemReferences:

• [1] 

  http://qlink08.sourceforge.net, open source.

• [2] 

  http://www.feynarts.de/cuba/, open source.

• [3] 

  http://www.mpfr.org/, open source.

• [4] 

  http://www.wolfram.com/products/mathematica/index.html.

     

SARAH 4: A tool for (not only SUSY) model builders

We present the new version of the Mathematica package SARAH which provides the same features for a non-supersymmetric model as previous versions for supersymmetric models. This includes an easy and straightforward definition of the model, the calculation of all vertices, mass matrices, tadpole equations, and self-energies. Also the two-loop renormalization group equations for a general gauge theory are now included and have been validated with the independent Python code PyR@TE. Model files for FeynArts, CalcHep/CompHep, WHIZARD and in the UFO format can be written, and source code for SPheno for the calculation of the mass spectrum, a set of precision observables, and the decay widths and branching ratios of all states can be generated. Furthermore, the new version includes routines to output model files for Vevacious for both, supersymmetric and non-supersymmetric, models. Global symmetries are also supported with this version and by linking Susyno the handling of Lie groups has been improved and extended.     

spinney: A Form library for helicity spinors

In this work, the library spinney is presented, which provides an implementation of helicity spinors and related algorithms for the symbolical manipulation program Form. The package is well suited for symbolic amplitude calculations both in traditional, Feynman diagram based approaches and unitarity-based techniques.

Program summary

Program title: spinneyCatalogue identifier: AEJQ_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEJQ_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.: 21 128No. of bytes in distributed program, including test data, etc.: 377 589Distribution format: tar.gzProgramming language: FormComputer: Any supporting the Form languageOperating system: Any supporting the Form languageClassification: 4.4, 5, 11.1Nature of problem: Implementation of the spinor-helicity formalismSolution method: Form implementationRunning time: From actual calculations of all six-point one-loop diagrams of the process bounds of 50 ms<t?71 s for the simplest and the most complicated diagram respectively have been derived on an Intel Xeon 3.20 GHz using Form 3.3.     

Synthesis of Moduli of Uniform Continuity by the Monotone Dialectica Interpretation in the Proof-system MinLog

We extract on the computer a number of moduli of uniform continuity for the first few elements of a sequence of closed terms t of Gödel's T of type . The generic solution may then be quickly inferred by the human. The automated synthesis of such moduli proceeds from a proof of the hereditarily extensional equality (≈) of t to itself, hence a proof in a weakly extensional variant of Berger-Buchholz-Schwichtenberg's system Z of . We use an implementation on the machine, in Schwichtenberg's MinLog proof-system, of a non-literal adaptation to Natural Deduction of Kohlenbach's monotone functional interpretation. This new version of the Monotone Dialectica produces terms in NbE-normal form by means of a recurrent partial NbE-normalization. Such partial evaluation is strictly necessary.     

C-library raft: Reconstruction algorithms for tomography. Applications to X-ray fluorescence tomography

There are many reconstruction algorithms for tomography, raft for short, and some of them are considered “classic” by researchers. The so-called raft library, provide a set of useful and basic tools, usually needed in many inverse problems that are related to medical imaging. The subroutines in raft are free software and written in C language; portable to any system with a working C compiler. This paper presents source codes written according to raft routines, applied to a new imaging modality called X-ray fluorescence tomography.

Program summary

Program title: raftCatalogue identifier: AEJY_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEJY_v1_0.htmlProgram obtainable from: CPC Program Library, Queen?s University, Belfast, N. IrelandLicensing provisions: GNU General Public Licence, version 2No. of lines in distributed program, including test data, etc.: 218 844No. of bytes in distributed program, including test data, etc.: 3 562 902Distribution format: tar.gzProgramming language: Standard C.Computer: Any with a standard C compilerOperating system: Linux and WindowsClassification: 2.4, 2.9, 3, 4.3, 4.7External routines:

•  

  raft:   autoconf 2.60 or later – http://www.gnu.org/software/autoconf/   GSL scientific library – http://www.gnu.org/software/gsl/   Confuse parser library – http://www.nongnu.org/confuse/

raft-fun: gengetopt – http://www.gnu.org/software/gengetopt/gengetopt.htmlNature of problem: Reconstruction algorithms for tomography, specially in X-ray fluorescence tomography.Solution method: As a library, raft covers the standard reconstruction algorithms like filtered backprojection, Novikov?s inversion, Hogan?s formula, among others. The input data set is represented by a complete sinogram covering a determined angular range. Users are allowed to set solid angle range for fluorescence emission at each algorithm.Running time: 1 second to 15 minutes, depending on the data size.     

CPsuperH2.0: An improved computational tool for Higgs phenomenology in the MSSM with explicit CP violation

We describe the Fortran code CPsuperH2.0, which contains several improvements and extensions of its predecessor CPsuperH. It implements improved calculations of the Higgs-boson pole masses, notably a full treatment of the 4×4 neutral Higgs propagator matrix including the Goldstone boson and a more complete treatment of threshold effects in self-energies and Yukawa couplings, improved treatments of two-body Higgs decays, some important three-body decays, and two-loop Higgs-mediated contributions to electric dipole moments. CPsuperH2.0 also implements an integrated treatment of several B-meson observables, including the branching ratios of Bs→μ⁺μ⁻, Bd→τ⁺τ⁻, Bu→τν, B→Xsγ and the latter's CP-violating asymmetry ACP, and the supersymmetric contributions to the mass differences. These additions make CPsuperH2.0 an attractive integrated tool for analyzing supersymmetric CP and flavour physics as well as searches for new physics at high-energy colliders such as the Tevatron, LHC and linear colliders.¹

Program summary

Program title: CPsuperH2.0Catalogue identifier: ADSR_v2_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/ADSR_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.: 13 290No. of bytes in distributed program, including test data, etc.: 89 540Distribution format: tar.gzProgramming language: Fortran 77Computer: PC running under Linux and computers in Unix environmentOperating system: LinuxRAM: 32 MbytesClassification: 11.1Catalogue identifier of the previous version: ADSR_v1_0Journal reference of the previous version: CPC 156 (2004) 283Does the new version supersede the previous version?: YesNature of problem: The calculations of mass spectrum, decay widths and branching ratios of the neutral and charged Higgs bosons in the Minimal Supersymmetric Standard Model with explicit CP violation have been improved. The program is based on recent renormalization-group-improved diagrammatic calculations that include dominant higher-order logarithmic and threshold corrections, b-quark Yukawa-coupling resummation effects and improved treatment of Higgs-boson pole-mass shifts. The couplings of the Higgs bosons to the Standard Model gauge bosons and fermions, to their supersymmetric partners and all the trilinear and quartic Higgs-boson self-couplings are also calculated. The new implementations include a full treatment of the 4×4(2×2) neutral (charged) Higgs propagator matrix together with the center-of-mass dependent Higgs-boson couplings to gluons and photons, two-loop Higgs-mediated contributions to electric dipole moments, and an integrated treatment of several B-meson observables.Solution method: One-dimensional numerical integration for several Higgs-decay modes, iterative treatment of the threshold corrections and Higgs-boson pole masses, and the numerical diagonalization of the neutralino mass matrix.Reasons for new version: Mainly to provide a coherent numerical framework which calculates consistently observables for both low- and high-energy experiments.Summary of revisions: Improved treatment of Higgs-boson masses and propagators. Improved treatment of Higgs-boson couplings and decays. Higgs-mediated two-loop electric dipole moments. B-meson observables.Running time: Less than 0.1 seconds.     

HidSecSOFTSUSY: Incorporating effects from hidden sectors in the numerical computation of the MSSM spectrum

SOFTSUSY is a software designed to solve the RG equations of the MSSM and compute its low energy spectrum. HidSecSOFTSUSY is an extension of the SOFTSUSY package which modifies the beta functions to include contributions from light dynamic fields in the hidden sector.

Program summary

Program title: HidSecSOFTSUSYCatalogue identifier: AEHP_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEHP_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.: 4167No. of bytes in distributed program, including test data, etc.: 141 411Distribution format: tar.gzProgramming language: C++, FortranComputer: Personal computerOperating system: Tested on GNU/LinuxWord size: 32 bitsClassification: 11.6External routines: Requires an installed version of SOFTSUSY (http://projects.hepforge.org/softsusy/)Nature of problem: Calculating supersymmetric particle spectrum and mixing parameters while incorporating dynamic modes from the hidden sector into the renormalization group equations. The solution to the equations must be consistent with a high-scale boundary condition on supersymmetry breaking parameters, as well as a weak-scale boundary condition on gauge couplings, Yukawa couplings and the Higgs potential parameters.Solution method: Nested iterative algorithm.Running time: A few seconds per parameter point.     

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     

yambo: An ab initio tool for excited state calculations

yambo is an ab initio code for calculating quasiparticle energies and optical properties of electronic systems within the framework of many-body perturbation theory and time-dependent density functional theory. Quasiparticle energies are calculated within the GW approximation for the self-energy. Optical properties are evaluated either by solving the Bethe-Salpeter equation or by using the adiabatic local density approximation. yambo is a plane-wave code that, although particularly suited for calculations of periodic bulk systems, has been applied to a large variety of physical systems. yambo relies on efficient numerical techniques devised to treat systems with reduced dimensionality, or with a large number of degrees of freedom. The code has a user-friendly command-line based interface, flexible I/O procedures and is interfaced to several publicly available density functional ground-state codes.

Program summary

Program title:yamboCatalogue identifier: AEDH_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEDH_v1_0.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: GNU General Public Licence v2.0No. of lines in distributed program, including test data, etc.: 149 265No. of bytes in distributed program, including test data, etc.: 2 848 169Distribution format: tar.gzProgramming language: Fortran 95, CComputer: any computer architecture, running any flavor of UNIXOperating system: GNU/Linux, AIX, Irix, OS/XHas the code been vectorised or parallelized?: YesRAM: 10-1000 MbytesClassification: 7.3, 4.4, 7.2External routines:

•

BLAS (http://www.netlib.org/blas/)

•

LAPACK (http://www.netlib.org/lapack/)

•

MPI (http://www-unix.mcs.anl.gov/mpi/) is optional.

•

BLACS (http://www.netlib.org/scalapack/) is optional.

•

SCALAPACK (http://www.netlib.org/scalapack/) is optional.

•

FFTW (http://www.fftw.org/) is optional.

•

netCDF (http://www.unidata.ucar.edu/software/netcdf/) is optional.

Nature of problem: Calculation of excited state properties (quasiparticles, excitons, plasmons) from first principles.Solution method: Many body perturbation theory (Dyson equation, Bethe Salpeter equation) and time-dependent density functional theory. Quasiparticle approximation. Plasmon-pole model for the dielectric screening. Plane wave basis set with norm conserving pseudopotentials.Unusual features: During execution, yambo supplies estimates of the elapsed and remaining time for completion of each runlevel. Very friendly shell-based user-interface.Additional comments:yambo was known as “SELF” prior to GPL release. It belongs to the suite of codes maintained and used by the European Theoretical Spectroscopy Facility (ETSF) [1].Running time: The typical yambo running time can range from a few minutes to some days depending on the chosen level of approximation, and on the property and physical system under study.References:[1] The European Theoretical Spectroscopy Facility, http://www.etsf.eu.     

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 Relic: A program for calculating relic density and flavor physics observables in Supersymmetry

We describe SuperIso Relic, a public program for evaluation of relic density and flavor physics observables in the minimal supersymmetric extension of the Standard Model (MSSM). SuperIso Relic is an extension of the SuperIso program which adds to the flavor observables of SuperIso the computation of all possible annihilation and coannihilation processes of the LSP which are required for the relic density calculation. All amplitudes have been generated at the tree level with FeynArts/FormCalc, and widths of the Higgs bosons are computed with FeynHiggs at the two-loop level. SuperIso Relic also provides the possibility to modify the assumptions of the cosmological model, and to study their consequences on the relic density.

Program summary

Program title: SuperIso RelicCatalogue identifier: AEGD_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEGD_v1_0.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: yesNo. of lines in distributed program, including test data, etc.: 2 274 720No. of bytes in distributed program, including test data, etc.: 6 735 649Distribution format: tar.gzProgramming language: C (C99 Standard compliant) and FortranComputer: 32- or 64-bit PC, MacOperating system: Linux, MacOSRAM: 100 MbClassification: 1.9, 11.6External routines: ISASUGRA/ISAJET and/or SOFTSUSY, FeynHiggsDoes the new version supersede the previous version?: No (AEAN_v2_0)Nature of problem: Calculation of the lightest supersymmetric particle relic density, as well as flavor physics observables, in order to derive constraints on the supersymmetric parameter space.Solution method: SuperIso Relic uses a SUSY Les Houches Accord file, which can be either generated automatically via a call to SOFTSUSY or ISAJET, or provided by the user. This file contains the masses and couplings of the supersymmetric particles. SuperIso Relic then computes the lightest supersymmetric particle relic density as well as the most constraining flavor physics observables. To do so, it calculates first the widths of the Higgs bosons with FeynHiggs, and then it evaluates the squared amplitudes of the diagrams needed for the relic density calculation. These thousands of diagrams have been previously generated with the FeynArts/FormCalc package. SuperIso Relic is able to perform the calculations in different supersymmetry breaking scenarios, such as mSUGRA, NUHM, AMSB and GMSB.Reasons for new version: This version incorporates the calculation of the relic density, which is often used to constrain Supersymmetry.Summary of revisions:

•

Addition of the relic density calculation

•

Replacement of "float" type by "double".

Unusual features: SuperIso Relic includes the possibility of altering the underlying cosmological model and testing the influence of the cosmological assumptions.Additional comments: This program is closely associated with the "SuperIso" program - CPC Program Library, Catalogue Id. AEAN.Running time:Compilation time: a couple of hours for the statically linked version, a few minutes for the dynamically linked version. Running time: about 1 second, or a few seconds if libraries need to be compiled on the fly.