Elsevier Science Home
Computer Physics Communications Program Library
Full text online from Science Direct
Programs in Physics & Physical Chemistry
CPC Home

[Licence| Download | New Version Template] aegv_v2_5.tar.gz(600 Kbytes)
Manuscript Title: SUSY FLAVOR v2.5: a computational tool for FCNC and CP-violating processes in the MSSM
Authors: J. Rosiek, P. Chankowski, A. Crivellin, A. Dedes, S. Jäger, P. Tanedo
Program title: SUSY FLAVOR v2.5
Catalogue identifier: AEGV_v2_5
Distribution format: tar.gz
Journal reference: Comput. Phys. Commun. 188(2015)208
Programming language: Fortran 77.
Computer: Any.
Operating system: Any, tested on Linux.
Keywords: Supersymmetry, K physics, B physics, rare decays, CP-violation.
PACS: 12.60.Jv, 13.20.He.
Classification: 11.6.

Does the new version supersede the previous version?: Yes

Nature of problem:
Predicting CP-violating observables, meson mixing parameters and branching ratios for a set of rare processes in the general R-parity conserving MSSM.

Solution method:
We use standard quantum theoretical methods to calculate Wilson coefficients in MSSM and at one loop including QCD corrections at higher orders when this is necessary and possible.

Reasons for new version:
The input/output routines have been rewritten to make them more flexible and compatible with the SLHA2 standard [1]. Calculations of the several processes implemented in earlier SUSY_FLAVOR versions have been corrected. New observables have been added. A number of bugs have been corrected.

Summary of revisions:
  1. Modified initialization routines. Currently the program should be able to read without modifications most of the SLHA2-compatible output files produced by other publicly available libraries calculating observables related to the MSSM phenomenology. In addition, new optional input block SFLAV_HADRON has been defined to facilitate modifications of the parameters related to the hadronic and QCD sector.
    The initialization sequence now goes through the following steps:

    • Before reading the file, all parameters are set to some initial values (which can be changed by editing the values given in the subroutine sflav_defaults in file sflav_io.f).

    • Subsequently, the user-defined data are read from the file with the default name susy_flavor.in. Data are grouped in Blocks following the SLHA2 specification or extensions described in [2]. Blocks are read in the following order: SOFTINP, SMINPUTS, VCKMIN, MINPAR (tan β only, other entries are ignored), EXTPAR, IMEXTPAR, MSL2IN, IMMSL2IN, MSE2IN, IMMSE2IN, TEIN, IMTEIN, TEINH, IMTEINH, MSQ2IN, IMMSQ2IN, MSU2IN, IMMSU2IN, MSD2IN, IMMSD2IN, TUIN, IMTUIN, TUINH, IMTUINH, TDIN, IMTDIN, TDINH, IMTDINH, SFLAV_HADRON.

    • The presence of any Block is optional - if some Block is absent, the program falls back to the default parameter values. At least flavor-diagonal SUSY mass parameters have to be defined, otherwise the vanishing default values cause the program to crash.

    • If a parameter is multiply defined in several Blocks, the value from Block read as latest in the list above overwrites (without warning!) the values from preceding Blocks.

    • Blocks do not need to be complete and to contain all the entries described in the SLHA2 specification - it is sufficient to define a minimal set of the parameters relevant for a given problem, others are filled with the default values.

    • The "non-holomorphic" LR mixing terms are not included in the SLHA2 specification and by default are set to 0. They can be initialized to the non-trivial values in the blocks TXINH and IMTXINH (X = E, D, U)

    • The new input block SFLAV_HADRON allows the modification of the hadronic- and QCD-related quantities used by SUSY_FLAVOR. The structure of this block and the default values of the hadronic parameters are shown at the end of the sample input file susy_flavor.in attached to the SUSY_FLAVOR distribution.

  2. New control variables have been added, allowing the seperate switching of contributions from various MSSM sectors on or off. They can be set by the following statement at the beginning of the driver program:
    call set_active_sector(ih,ic,in,ig)
    where the variables ih, ic, in, ig can take values 0 or 1 and they control, respectively, the inclusion in the total result of the diagrams with gauge and Higgs bosons, charginos, neutralinos and gluinos exchanged in the loops. Note that diagrams with Higgs and gauge bosons are always added together and currently cannot be disentangled, so setting ih=1, ic=in=ig=0 does not reproduce the SM result. By default, if no call to set_active_sector is made, all control variables are assumed to be equal to 1, so that all contributions are included.

  3. Added or modified processes:
    • The expressions used to calculate the neutron Electric Dipole Moment have been modified.

    • The branching ratios for the radiative decays of the heavy lepton into the lighter lepton and the photon, μ → eγ and τ → eγ, μγ, are now normalized to the total heavy lepton decay width (previously they were normalized to the decay width into leptonic channels).

    • The routines calculating branching ratios of B → τν and BDτν decays have been generalized to include more general structure of the effective Higgs boson-fermion couplings. In addition the routine calculating Br(BD* τν) has been added.

    • The routines for rare decays of the top quark to the CP-even Higgs boson and the light quarks, t → ch, uh, have been added, based on Ref. [3] (program can calculate also the decay rates of the top quark to the heavier CP-even Higgs boson H, assuming that such decays are allowed kinematically).

    • The routine calculating the approximate 2-loop estimate of the neutral CP-even Higgs mass mh has been added, based on Ref. [4]. Note that for the more precise calculations of this mass other publicly available SUSY codes should be used.

    • Default values of numerous quantities which are treated by SUSY_FLAVOR as the external parameters, mainly the values of hadronic parameters obtained from lattice calculations and results of experimental measurements, have been updated to accommodate the latest published results.

  4. SUSY_FLAVOR's output is now written to the file named susy_flavor.out. It has "SLHA-like" structure, i.e. it is divided into "data blocks", however these blocks are SUSY_FLAVOR specific and do not follow the common SLHA2 standards. The output file contains the following data blocks:
    • SFLAV_CONTROL: control variables and error code status.
    • SFLAV_MASS: MSSM mass spectrum after mass matrix diagonalization.
    • SFLAV_CHIRAL_YUKAWA: relative size of the resummed chiral corrections to the Yukawa couplings.
    • SFLAV_CHIRAL_CKM: relative size of the resummed chiral corrections to the CKM matrix elements.
    • SFLAV_DELTA_F0: ΔF = 0 observables: leptonic EDMs and g - 2 anomalies, neutron EDM.
    • SFLAV_DELTA_F1: ΔF = 1 observables: decay rates of ll'γ, K → πνν, B+ → τ+ν, BDτν, BD*τν, BXsγ, Bd,sli+lj-, t → uh, t → ch.
    • SFLAV_DELTA_F2: ΔF = 2 observables: εK, ΔmK, ΔmD, ΔmBd, ΔmBs.

    Blocks SFLAV_CHIRAL_YUKAWA and SFLAV_CHIRAL_CKM show the relative differences of bare and physical Yukawa couplings and CKM matrix elements after the resummation of chiral corrections. If they are large, ≥ O(1), the perturbation expansion is not reliable and the remaining program output may not be correct.

  5. The new integrated manual for SUSY_FLAVOR_v2.5 has been created, including the detailed description of the modifications listed above. It is attached to the SUSY_FLAVOR distribution. Regular code distribution updates and bug fixes (between the major revisions submitted to Computer Physics Communications) can be found on the program web page www.fuw.edu.pl/susy_flavor.

The results apply only to the case of MSSM with R-parity con- servation and without heavy right neutrino sector [5].

Additional comments:
This program has been catalogued as AEGV_v2_5 to conform to the manuscript version number. There are no programs catalogued as, AEGV_v2_1, AEGV_v2_2, AEGV_v2_3, AEGV_v2_4, in the CPC Program Library.

Running time:
For a single parameter set, under 1s on a personal computer.

[1] B. Allanach et al., Comput. Phys. [arXiv:0801.0045 [hep-ph]]. Commun. 180 (2009) 8
[2] J. Rosiek, P. Chankowski, A. Dedes, S. Jager and P. Tanedo, Comput. Phys. Commun. 181 (2010) 2180 [arXiv:1003.4260 [hep-ph]]; A. Crivellin, J. Rosiek et. al., Comput. Phys. Commun. 184 (2013) 1004, [arXiv:1203.5023" [hep-ph]].
[3] A. Dedes, M. Paraskevas, J. Rosiek, K. Suxho and K. Tamvakis, arXiv:1409.6546 [hep-ph].
[4] S. Heinemeyer, W. Hollik and G. Weiglein, Phys.Lett. B455 (1999) 179-191 [hep-ph/9903404]
[5] A. Dedes, H. Haber and J. Rosiek, [arXiv:0707.3718 [hep-ph]].