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.     

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     

micrOMEGAs_3: A program for calculating dark matter observables

micrOMEGAs  is a code to compute dark matter observables in generic extensions of the standard model. This new version of micrOMEGAs  is a major update which includes a generalization of the Boltzmann equations to accommodate models with asymmetric dark matter or with semi-annihilation and a first approach to a generalization of the thermodynamics of the Universe in the relic density computation. Furthermore a switch to include virtual vector bosons in the final states in the annihilation cross sections or relic density computations is added. Effective operators to describe loop-induced couplings of Higgses to two-photons or two-gluons are introduced and reduced couplings of the Higgs are provided allowing for a direct comparison with recent LHC results. A module that computes the signature of DM captured in celestial bodies in neutrino telescopes is also provided. Moreover the direct detection module has been improved as concerns the implementation of the strange “content” of the nucleon. New extensions of the standard model are included in the distribution.     

SuperIso v2.3: A program for calculating flavor physics observables in supersymmetry

We describe SuperIso v2.3 which is a public program for evaluation of flavor physics observables in the minimal supersymmetric extension of the Standard Model (MSSM). SuperIso v2.3, in addition to the isospin asymmetry of B→K^∗γ, which was the main purpose of the first version, incorporates new flavor observables such as the branching ratio of Bs→μ⁺μ⁻, the branching ratio of B→τντ, the branching ratio of B→Dτντ and the branching ratio of K→μνμ. The calculation of the branching ratio of B→Xsγ is also improved in this version, as it now includes NNLO Standard Model contributions in addition to partial NLO supersymmetric contributions. The program also computes the muon anomalous magnetic moment (g−2). Four sample models are included in the package, namely mSUGRA, NUHM, AMSB and GMSB. 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, or provided by the user. The calculation of the observables is detailed in the Appendices, where a suggestion for the allowed intervals for each observable is also provided.

Program summary

Program title: SuperIsoCatalogue identifier: AEAN_v2_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEAN_v2_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.: 5977No. of bytes in distributed program, including test data, etc.: 39 375Distribution format: tar.gzProgramming language: C (C99 Standard compliant)Computer: 32- or 64-bit PC, MacOperating system: Linux, MacOSRAM: less than 1 MbClassification: 11.6Catalogue identifier of previous version: AEAN_v1_0Journal reference of previous version: Comput. Phys. Comm. 178 (2008) 745External routines: ISASUGRA/ISAJET and/or SOFTSUSYDoes 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, in order to derive constraints on the supersymmetric parameter space.Solution method: SuperIso 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 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.Reasons for new version: This new version incorporates the calculation of several additional observables, and the inclusive branching ratio of b→sγ is now computed at NNLO accuracy for the Standard Model. The implemented routines are therefore extensively modified.Summary of revisions:

•

Compatibility with the SLHA2 input file format

•

Implementation of the calculation of the muon anomalous magnetic moment

•

Implementation of observables related to leptonic and semi-leptonic B meson decays

•

Implementation of observables related to K meson decays

•

Improvement of the calculations of the branching ratio of b→sγ (now at NNLO accuracy) and the isospin asymmetry of B→K∗γ

•

Update of parameters to their latest values

Unusual features: The code is very modular, and new routines for calculating new observables can be easily added.Running time: less than 1 sec     

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.     

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.     

EPW: A program for calculating the electron–phonon coupling using maximally localized Wannier functions

EPW (Electron–Phonon coupling using Wannier functions) is a program written in Fortran90 for calculating the electron–phonon coupling in periodic systems using density-functional perturbation theory and maximally localized Wannier functions. EPW can calculate electron–phonon interaction self-energies, electron–phonon spectral functions, and total as well as mode-resolved electron–phonon coupling strengths. The calculation of the electron–phonon coupling requires a very accurate sampling of electron–phonon scattering processes throughout the Brillouin zone, hence reliable calculations can be prohibitively time-consuming. EPW combines the Kohn–Sham electronic eigenstates and the vibrational eigenmodes provided by the Quantum ESPRESSO package (see Giannozzi et al., 2009 [1]) with the maximally localized Wannier functions provided by the wannier90 package (see Mostofi et al., 2008 [2]) in order to generate electron–phonon matrix elements on arbitrarily dense Brillouin zone grids using a generalized Fourier interpolation. This feature of EPW leads to fast and accurate calculations of the electron–phonon coupling, and enables the study of the electron–phonon coupling in large and complex systems.

Program summary

Program title: EPWCatalogue identifier: AEHA_v1_0Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEHA_v1_0.htmlProgram obtainable from: CPC Program Library, Queen's University, Belfast, N. IrelandLicensing provisions: GNU Public LicenseNo. of lines in distributed program, including test data, etc.: 304 443No. of bytes in distributed program, including test data, etc.: 1 487 466Distribution format: tar.gzProgramming language: Fortran 90Computer: Any architecture with a Fortran 90 compilerOperating system: Any environment with a Fortran 90 compilerHas the code been vectorized or parallelized?: Yes, optimized for 1 to 64 processorsRAM: Heavily system dependent, as small as a few MBSupplementary material: A copy of the “EPW/examples” directory containing the phonon binary files can be downloadedClassification: 7External routines: MPI, Quantum-ESPRESSO package [1], BLAS, LAPACK, FFTW. (The necessary Blas, Lapack and FFTW routines are included in the Quantum-ESPRESSO package [1].)Nature of problem: The calculation of the electron–phonon coupling from first-principles requires a very accurate sampling of electron–phonon scattering processes throughout the Brillouin zone; hence reliable calculations can be prohibitively timeconsuming.Solution method: EPW makes use of a real-space formulation and combines the Kohn–Sham electronic eigenstates and the vibrational eigenmodes provided by the Quantum-ESPRESSO package with the maximally localized Wannier functions provided by the wannier90 package in order to generate electron–phonon matrix elements on arbitrarily dense Brillouin zone grids using a generalized Fourier interpolation.Running time: Single processor examples typically take 5–10 minutes.References:

• [1] 

  P. Giannozzi, et al., J. Phys. Condens. Matter 21 (2009), 395502, http://www.quantum-espresso.org/.

     

eett6f v. 1.0: A program for top quark pair production and decay into 6 fermions at linear colliders

The first version of a computer program eett6f for calculating cross sections of e⁺e⁻→6 fermions processes relevant for a -pair production and decay at centre of mass energies typical for linear colliders is presented. eett6f v. 1.0 allows for calculating both the total and differential cross sections at tree level of the Standard Model (SM). The program can be used as the Monte Carlo generator of unweighted events as well.     

OptaDOS: A tool for obtaining density of states,core-level and optical spectra from electronic structure codes

We present OptaDOS, a program for calculating core-electron and low-loss electron energy loss spectra (EELS) and optical spectra along with total-, projected- and joint-density of electronic states (DOS) from single-particle eigenenergies and dipole transition coefficients. Energy-loss spectroscopy is an important tool for probing bonding within a material. Interpreting these spectra can be aided by first principles calculations. The spectra are generated from the eigenenergies through integration over the Brillouin zone. An important feature of this code is that this integration is performed using a choice of adaptive or linear extrapolation broadening methods which we show produces higher accuracy spectra than standard fixed-width Gaussian broadening. OptaDOS  may be straightforwardly interfaced to any electronic structure code. OptaDOS  is freely available under the GNU General Public licence from http://www.optados.org.     

A precise and robust quartz sensor based on tuning fork technology for (SF6)-gas density control

A novel sensor technology for measuring and monitoring gas density is described. The gas density sensor is self-calibrating and specifically designed to operate as a sulfur hexafluoride (SF₆) monitor and control device in gas insulated high-voltage switchgears (GIS), but can also be applied for density measurement of any kind of gas. It comprises a pair of tuning forks oscillating at their resonance frequency. One oscillator is exposed to the gas to be monitored, the other one is used for comparison and temperature compensation. Exposure to gas leads to a shift in the resonance frequency proportional to the gas density. A density standard based on a combined weight and mass measurement has been performed with a precision better than 0.02%, giving proof of exact sensor calibration. The gas density sensor works with a precision better than 0.5% over a range of 50 kg/m³ SF₆ gas density. Sensor response time is less than 50 ms. The real gas equation of Beattie and Bridgman fits the gauged values with a mean error of 0.125%.     

Comment on “A direct approach with computerized symbolic computation for finding a series of traveling waves to nonlinear equations” by Engui Fan and H.H. Dai: [Comput. Phys. Commun. 153 (2003) 17]

Fan and Dai [Comput. Phys. Commun. 153 (2003) 17] have found a series of traveling wave solutions for nonlinear equations by applying a direct approach with computerized symbolic computations. They have claimed that the proposed method, in comparison with most existing symbolic computation methods such as a tanh method and Jacobi function method, not only give new and more general solutions, but also provides a guideline to classify various types of the solution according to some parameters. We show that the claims by Fan and Dai are wrong since some of the solutions do not satisfy the differential equation that they have adopted for the algebraic method.