CPsuperH2.3: an Updated Tool for Phenomenology in the MSSM with Explicit CP Violation

We describe the Fortran code CPsuperH2.3, which incorporates the following updates compared with its predecessor CPsuperH2.0. It implements improved calculations of the Higgs-boson masses and mixing including stau contributions and finite threshold effects on the tau-lepton Yukawa coupling. It incorporates the LEP limits on the processes e^+ e^- to H_i Z, H_i H_j and the CMS limits on H_i to tau^+ tau^- obtained from 4.6/fb of data at a centre-of-mass energy of 7 TeV. It also includes the decay mode H_i to Z gamma and the Schiff-moment contributions to the electric dipole moments of Mercury and Radium225, with several calculational options for the case of Mercury. These additions make CPsuperH2.3 a suitable tool for analyzing possible CP-violating effects in the MSSM in the era of the LHC and a new generation of EDM experiments


Introduction
Supersymmetry is one of the most attractive possible scenarios for physics beyond the Standard Model in the TeV energy range, and the minimal supersymmetric extension of the Standard Model (MSSM) is the simplest framework that incorporates this scenario.The MSSM allows for many possible CP-violating phases, notably those in the SU(3), SU (2) and U(1) gaugino masses M 3,2,1 and in the trilinear soft supersymmetry-breaking parameters associated with the third-generation Yukawa couplings A t,b,τ .
CPsuperH [1,2] is an evolving software tool that incorporates the effects of these CPviolating phases into the calculation of Higgs boson masses, couplings and other properties.The previous version, CPsuperH2.0[1], also incorporated a number of B-physics observables, including the branching ratios for B s → µ + µ − , B d → τ + τ − , B u → τ ν, B → X s γ and the CP-violating asymmetry in the latter decay, A CP , as well as the supersymmetric contributions to the B 0 s,d − B0 s,d mass differences.Also included in CPsuperH2.0were calculations of two-loop supersymmetric Higgs contributions to the electric dipole moments (EDMs) of Thallium, the electron and muon † .
After the recent discovery at the LHC [4] of a new particle resembling a Standard-Model-like Higgs boson, Higgs physics and searches for associated CP-violating effects are entering a new era, in which the LHC experiments are now also probing directly the possible existences of heavier supersymmetric Higgs bosons.In parallel, the LHCb experiment is taking experiments on B-physics to a new level of precision, and a new round of experiments is set to improve significantly sensitivities to a wider range of EDMs, including also Mercury and Radium.This juncture in the exploration of physics at the TeV scale, both direct and indirect, is the appropriate moment to document the capabilities of a new update of CPsuperH.
The followings are the main new features incorporated in CPsuperH2.3.As described in Section 2, it features an improved treatment of CP-violating effects on Higgs-boson masses and mixing, with stau contributions included as in Ref. [5].We also take consistently into account finite radiative effects on the tau-lepton Yukawa coupling.As described in Section 3, CPsuperH2.3also features an implementation of the LEP limits on the processes e + e − → H i Z, H i H j [6], as well as the limits on H i → τ τ obtained by CMS with 4.6 fb −1 of LHC data at a centre-of-mass energy of 7 TeV [7], where H i (with i = 1, 2, 3) denotes collectively the three neutral Higgs bosons H 1,2,3 in the CP-violating MSSM.In the same section, we also present theoretical predictions for H i → Zγ, as these decay modes can become detectable as the integrated LHC luminosity increases.Section 4 outlines how CPsuperH2.3incorporates calculations of the EDMs of Mercury and 225 Ra that include estimates of the contributions due to Schiff moments [8].Each section includes figures that illustrate some typical results obtained using CPsuperH2.3.As explained in Section 5, an SLHA2 interface was created to facilitate the comparison and linkage of output data between CPsuperH2.3and other public codes, in accordance with the SUSY Les Houches Accord [10,11].Finally, Section 6 summarizes our updates to the code.

Stau Contributions to Higgs masses and Mixing
The scalar tau contributions to the masses and mixing of the neutral Higgs bosons have been included as in Ref. [5], similarly to the sbottom corrections with More precisely, in the (φ 1 , φ 2 , a) T basis, the scalar tau contributions to the entries of the neutral Higgs mass squared matrix from the one-loop effective potential are given by where ĝ2 ≡ (g 2 + g ′2 )/4 and the CP-violating quantity is The tan β-enhanced terms containing µ and A τ are included in R τ and R ′ τ , given by and the loop function Since the stau contributions to the mass splitting M 2 H ± − M 2 A are small, we have neglected them.For the same reason, we set Note that in the limit of small mass splitting one has 2 )/2.So, in this limit, the leading contributions proportional to tan 4 β are given by which are the same as those presented in [12][13][14].Note that here we take into account the τ 1,2 mass splitting.
In Fig. 1, we demonstrate the impact of the corrections from the light scalar taus for large µ = 1 TeV for several values of tan β.When m τ 1,2 ∼ 1 TeV, the stau corrections are negligible and the lightest Higgs boson mass increases by the small amount of ∼ 0.05 GeV even when tan β = 60.However, these corrections are larger for lighter staus, and the lightest Higgs boson mass may decrease by as much as ∼ 4 GeV, if tan β = 60 and m τ 1 ∼ 100 GeV.For comparison, as shown in Table 4, the masses and mixing matrix of the neutral Higgs bosons without including the stau effects are stored in the array RAUX H: We have fixed M H ± = 1 TeV.The thin short line around m τ 1 = 1 TeV for each tan β value shows the lightest Higgs boson mass before the inclusion of the stau effects.

LEP Limits
We have implemented the LEP limits from the processes e + e − → H i Z and e + e − → H i H j using Table 14 (a) and Table 17 (a) of Ref. [6], respectively, assuming that the Higgs decay patterns are not drastically different from those in the Standard Model and in the MHMAX scenario with tan β = 10.We have required that each of the Higgs couplings to the gauge boson(s) normalized to the Standard Model values, g 2 H i V V and g 2 H i H j Z , should be smaller than the corresponding values in the Tables.The result is saved in -RAUX H(430)=ILEP : 0 (Excluded) or 1 (Allowed).
For the more general cases, we refer to more refined tools such as HiggsBounds [15].

LHC Limits
We have also incorporated the recent CMS limit from the search for the MSSM neutral Higgs bosons decaying into tau pairs based on 4.6 fb −1 [7], and the result is saved in -RAUX H(440)= ILHC7 H→τ τ 4.6 : 0 (Excluded) or 1 (Allowed) .
In our implementation of the CMS limit, assuming the same K factors as in the SM, we approximate the production cross sections of the neutral Higgs bosons at the LHC as follows: where P = ggH i , bbH i , V V H i specifies each production process.The SM production cross sections are calculated by using HIGLU [16], BBH@NNLO [17] and HAWK [18]  § .In leading order, the process-dependent ratios are given by: where S g SM (M H i ) = f =b,t F sf (τ if ).These factors are saved in the array CAUX H(221-223) as shown in Table 6.
where the extra factor in the second term accounts for the phase-space difference between scalar and pseudo-scalar decays.In the program, we drop the extra factor since the ratio is involved in production process.
In Figs.2,3,4, we show the three process-dependent ratios as functions of the corresponding Higgs masses, assuming the MHMAX scenario.For completeness, also shown in Fig. 5 are the ratios of the Higgs-boson couplings to two photons defined by  (14), as functions of their masses for several values of tan β.In the CP-conserving limit, Φ ≡ h, H, A with h(H) and A being the lighter (heavier) CPeven and CP-odd Higgs bosons, respectively.We assume the MHMAX scenario: where , which are saved in the array CAUX H(231-233) as shown in Table 6.
We observe that the branching ratio B(H 1 → γγ) may be enhanced if the lighter stau has a low mass, particularly for large tan β as shown in Fig. 6.This enhancement is consistent with the signals observed [4], which may be larger than in the Standard Model.Furthermore, in Fig. 7 we show the correlation of B(H 1 → γγ) with B(H 1 → Zγ) in the MSSM, after the branching ratios are normalized to their SM predictions.Calculational details of B(H i → Zγ) are described in Appendix A. Once the production cross sections of the neutral MSSM Higgs bosons have been calculated, we may require the sum to be smaller than the observed limit for a given value M H i of a specified neutral Higgs boson H i .Here we are adding all Higgs production cross sections of φ = H i to the one of the H i boson when their mass difference is smaller than δM.We take δM = 0.21 * 130/2 ∼ 13 GeV, corresponding to the tau-pair mass resolution of ∼ 21 % at a Higgs boson mass of 130 GeV [7].Such a simple treatment is fairly accurate, as long as δM ≫ Γ H i , namely in the absence of strongly overlapping Higgs resonances [9].Under these considerations, Fig. 8 shows the LEP-and LHC-excluded regions in the M A -tan β plane for the MHMAX scenario.We find our results reasonably consistent with those presented in Fig. 4 of [7] when M A ≥ 90 GeV.
Finally, for the value given by M H SM =SMPARA H(20), the decay widths and branching ratios of the SM Higgs have been calculated, and the results are stored in RAUX H(600-700),  2 but for the couplings of the Higgs bosons to two vector bosons (16).
see Table 4.For comparison with other works, we display in Fig. 9 the decay widths into b quarks and photons (upper) and the branching ratios into W and Z bosons (lower) as functions of the SM Higgs-boson mass.We note that decay widths into off-shell vector bosons using the four-body phase-space have been implemented in CPsuperH2.3.

Mercury and Radium Electric Dipole Moments
CPsuperH2.0 included two-loop Higgs-mediated contributions to the Thallium and electron EDMs [1], and CPsuperH2.3extends these results to include calculations of the Mercury and 225 Ra EDMs that incorporate Schiff-moment contributions.In the case of the Mercury EDM, these are parameterized as follows [19]: where the ḡ(0), (1) πN N are the CP-odd πNN couplings, and we refer to Ref. [8] for further details.We note that d I Hg is basically the same as that calculated in CPsuperH2.0.In the case of the Radium 225 Ra, we estimate the EDM through [8]  -RAUX H(420) = d Ra as shown in Table 3.

SLHA2 interface
CPsuperH2.3 also provides an output in accordance with the SUSY Les Houches Accords 1 (SLHA1) [10] and 2 (SLHA2) [11].By taking * * • IFLAG H(30)=1 in the run file, the output file cpsuper2.3slha2.out is generated.The output file of the current version includes the following blocks: • MODSEL: In the block MODSEL, CP violation with completely general CP phases has been selected.
• SMINPUTS and VCKMIN: pole masses of top-quark, electron and muon, and m d,u,s (2 GeV) MS are given in the SMINPUTS block.In the VCKMIN block, the four parameters for the CKM mixing matrix λ, A, ρ, and η are given.
• EXTPAR and IMEXTPAR: The blocks for non-minimal parameters have been filled according to Section 2 of the SLHA2 writeup [11]  † † .For the input scale M input , we are taking Q tb , the scale of the heaviest third-generation squark [20], which is stored in RAUX H(13)= Q 2 tb .The imaginary parts of the gaugino masses, the trilinear couplings, and the µ parameters are listed in the block IMEXTPAR.• HCOUPLINGS, IMHCOUPLINGS: In these blocks, the seven Higgs-self couplings λ 1−7 calculated according to Ref. [21] are given from CAUX H(201−207).
• THRESHOLD: In this block, the measures of the threshold corrections to the Yukawa couplings h t,b are given from CAUX H(211,212).The quantity ∆ b is also available from CAUX H(10).
• DECAY: The total decay widths and the non-vanishing branching ratios of the neutral and charged Higgs bosons are stored in the block from the arrays GAMBRN(101,1,IH) and GAMBRN(IM,3,IH) with IH= 1−3 and GAMBRC(51,1) and GAMBRC(IM,3), respectively.The decay width and branching ratios of the top quark are also stored from the arrays RAUX H(50-53).
• FOBS: In this block, we show the branching ratios • FDIPOLE: Here, we list the EDMs of Thallium, neutron, Mercury, Deuteron, Radium, and muon, as well as the anomalous magnetic moment of muon, (g µ − 2).In the cases of the neutron and Mercury EDMs, we present 3 and 4 evaluations, respectively, for the estimation of the theoretical uncertainties.
• HiggsBoundsInputHiggsCouplingsBosons, HiggsBoundsInputHiggsCouplingsFermions: In these blocks, the Higgs couplings normalized to the corresponding SM couplings are listed for the interface to the program HiggsBounds [15].

Summary
Encouraged by the recent observation of a new particle resembling a SM-like Higgs boson at the CERN Large Hadron Collider [4], in CPsuperH2.3we have performed a number of updates to the CPsuperH code.In detail, we have improved the computation of CPviolating effects on Higgs-boson masses and mixing, by including stau contributions [5] and finite radiative effects on the tau-lepton Yukawa coupling.These effects play an important role in the decays H 1,2,3 → γγ, for large values of tan β [14,22].We have also implemented the LEP limits on the processes e + e − → H i Z, H i H j , including the bounds on H i → τ τ obtained by CMS with 4.6 fb −1 of LHC data at a centre-of-mass energy of 7 TeV.Finally, we have included the decay modes of the neutral Higgs bosons H 1,2,3 into a Z boson and a photon.These decay modes are expected to become observable as more luminosity is accumulated at the LHC.
We have also incorporated in CPsuperH2.3calculations of the EDMs of Mercury and 225 Ra, including estimates of the contributions due to Schiff moments.We have presented a number of numerical examples and figures for each of our updates, in order to illustrate some typical results obtained by CPsuperH2.3.To enhance the synergy and compatibility of CPsuperH2.3with other codes, we have created an SLHA2 interface in accordance with the SUSY Les Houches Accords.For comparison with other works, we include evaluations of the principal decay rates and branching ratios of a SM Higgs boson.In the Appendix, we exhibit a number of output Tables, in order to highlight the updates with respect to the older version of CPsuperH.
If indeed the new particle discovered by ATLAS and CMS turns out to be a Higgs boson, the joint task of theory and experiment will be to establish whether it is compatible with the SM or exhibits deviations characteristic of some extension of the SM such as the MSSM.The possibility of CP violation should be considered in studying the latter possibility, and CPsuperH2.3 is a suitable updated tool for this task.

Appendix
Most importantly, in addition to the SLHA2 interface, the array for the SM parameters SMPARA H has been extended to include the SM Higgs mass,M H SM , see Table 1.Because of these improvements, from CPsuperH2.2 to CPsuperH2.3, the following changes to the files run and cpsuperh2.fare needed: • In the run file: the entries SMPARA H (20) and IFLAG H(30) added • In the cpsuperh2.ffile: -REAL*8 SMPARA H (19),SSPARA H(38) =⇒ REAL*8 SMPARA H (20),SSPARA H(38) Furthermore, in the updated version, the contents of the array GAMBRN(IM,IWB=2,IH) for the neutral Higgs-boson decays now include the corresponding SM branching ratios.To reproduce Figs. 6 and 7, for example, one simply needs to use the following ratios of parameter arrays: .

A H i → Zγ
For the calculation of B(H i → Zγ), we closely follow Ref. [27].
The amplitude for the decay process H i → Z(k 1 , ǫ 1 ) γ(k 2 , ǫ 2 ) can be written as where k 1,2 are the momenta of the Z boson and the photon (we note that 2k The decay width is given by The scalar and pseudoscalar form factors are given by with , The loop functions are ‡ ‡ We follow the conventions and notations of CPsuperH [2] for the Higgs couplings to the SM and supersymmetric particles, and the relevant Z-boson interactions are given by the following Lagrangian terms: where g Z = e/(s W c W ) and , we refer to [27].
• Z-chargino-chargino [28] L where For completeness, we recall that the Z-boson couplings to the quarks and leptons are given by

Figure 2 :
Figure 2: The couplings of the MSSM Higgs bosons Φ to two gluons normalized to the Standard Model values, R H i gg(14), as functions of their masses for several values of tan β.In the CP-conserving limit, Φ ≡ h, H, A with h(H) and A being the lighter (heavier) CPeven and CP-odd Higgs bosons, respectively.We assume the MHMAX scenario:m Q3 = m Ũ3 = m D3 = m L3 = m Ẽ3 = M SUSY =1 TeV, µ = 200 GeV , M 1 = 100 GeV , M 2 = 200 GeV , M 3 = 800 GeV, and A t = √ 6 M SUSY + µ/ tan β with A b = A τ = A t .

Figure 3 :
Figure 3: The same as in Fig. 2 but for the couplings of the Higgs bosons to b quarks (15).

Figure 4 :
Figure4: The same as in Fig.2but for the couplings of the Higgs bosons to two vector bosons(16).

Figure 5 :
Figure 5: The same as in Fig. 2 but for the couplings of the Higgs bosons to two photons.

Fig. 10 Figure 6 :
Fig.10shows the four estimates of the Mercury EDM and the Radium EDM as functions of the parameter ρ that parameterizes the hierarchy between the first two and third generations in the CPX scenario.The new EDMs are stored in the array RAUX H:

Figure 7 :
Figure 7: Correlations of the MSSM-to-SM branching ratio of the lightest Higgs boson H 1 decaying into two photons with the respective one for the decay H 1 into a Z boson and a photon, for discrete choices of tan β.The scenario and lines are the same as in Fig. 1.
2 and α = ( t L , t R ), and similarly for sbottoms and staus.

Figure 9 :
Figure 9: The upper frames are for the decay widths of the SM Higgs boson into b quarks (upper left) and photons (upper right) as functions of its mass.The lower frames are for the branching ratios into W (lower left) and Z (lower right) bosons.In the lower frames, the solid (dotted) lines are obtained using 4 (3)-body phase space below the mass thresholds.The vertical lines correspond to M H SM = 125.5 GeV.

Fermilab
is operated by Fermi Research Alliance, LLC under Contract No. DE-AC02-07CH11359 with the U.S. Department of Energy.Work at ANL is supported in part by the U.S. Department of Energy under Contract No. DE-AC02-06CH11357.
eff|B 0 s SUSY given in CAUX H(150,151) to the SM contributions.

Table 1 :
The contents of the extended SMPARA H(IP).

Table 2 :
The contents of the extended SSPARA H(IP).

Table 7 :
The updated contents of the array GAMBRN(IM,IWB=1,IH) for the decays of the neutral Higgs bosons H IH .The entry GAMBRN(IM,IWB=1,IH) is for the decay width of the decay mode IM in GeV.And the entries GAMBRN(IM,IWB=2,IH) and GAMBRN(IM,IWB=3,IH) are the corresponding SM and the full SUSY branching ratios, respectively.