Constraints on the CP-Violating MSSM

We discuss the prospects for observing CP violation in the MSSM with six CP-violating phases, using a geometric approach to maximise CP-violating observables subject to the experimental upper bounds on electric dipole moments. We consider constraints from Higgs physics, flavour physics, the dark matter relic density and spin-independent scattering cross section with matter.


Introduction
In the minimal supersymmetric extension of the Standard Model (MSSM), there are possibilities of having new sources of CP violation beyond the Cabibbo Kobayashi Maskawa (CKM) phase. However, experimental constraints in particular from the measurements of electric dipole moments (EDMs, see Tab. 1) limit strongly the values of the additional phases. Here we consider the maximally CP-violating, minimally flavour-violating (MCPMFV) model that contains 6 new CP-violating phases: 3 phases Φ 1,2,3 in the masses of the U(1), SU(2) and SU(3) gauginos, and 3 phases Φ t,b,τ in the trilinear soft supersymmetry breaking couplings A t,b,τ of the third-generation stop, sbottom and stau sfermions, respectively. Brute force to sample randomly the values of the phases imposing all the experimental constraints (from Higgs, flavour physics, dark matter and EDMs) appears to be inefficient and one would need an optimised procedure such as a geometric approach [1,2]. We study signatures of CP violation considering the lightest neutral Higgs boson to be $ Based on the talk by A.A. at the Sixth Workshop on Theory, Phenomenology and Experiments in Flavour Physics, Capri, June 2016.
Email address: alexandre.arbey@ens-lyon.fr (A. Arbey  the observed one and the lightest supersymmetric particle to be the lightest neutralino, and study in particular the CP asymmetry in b → sγ, B s meson mixing ∆M B s , and CP-violating couplings of the heavier neutral Higgs bosons in several representative MSSM scenarios. Extended discussions can be found in [2].

Methodology
We consider here two different MSSM scenarios: the CPV-CMSSM and the CPV-pMSSM, which are extensions of the standard CMSSM and pMSSM scenarios incorporating six additional phases Φ 1,2,3,t,b,τ to account arXiv:1611.05061v1 [hep-ph] 15 Nov 2016 for CP violation. We perform flat random scans on the standard parameters, similarly to [7][8][9]. The phases however are very severely constrained by the electric dipole moment (EDM) measurements given in Tab. 1, so that the use of an optimised geometric technique becomes necessary. The muon EDM has not been used in our analysis since its experimental upper bound [10] only provides very weak constraints that are not competitive with the constraints from the other EDMs.
Let us consider the four EDMs, E a,b,c,d , of Tab. 1 in the small phase approximation, with where Φ ≡ Φ α = Φ 1,2,3,t,b,τ and E i ≡ ∂E i /∂Φ, and an additional CP-violating observable O in the the small phase approximation, such that The optimal direction that maximises the observable O and is orthogonal to the EDM vectors E a,b,c,d α is given by with an unknown normalisation factor.
For each choice of the CP-conserving parameters, we fix the phases to 0 • or ±180 • , and compute the optimal direction using Eq. (3). Sets of phases are chosen randomly along this direction, then moved by 20 • along the favoured direction. Subsequently, the favoured direction is recomputed at this new position. We iterate this procedure up to 100 • .
The SUSY mass spectra and couplings, and EDM constraints are computed with CPsuperH [11][12][13]. The thorium monoxide EDM is computed following [14]. Flavour observables are calculated with SuperIso [15,16] and CPsuperH. The dark matter relic density is computed with SuperIso Relic [17] and micrOMEGAs [18], and the later is also used to compute dark matter direct detection observables. Finally, HiggsBounds [19] is used to impose the Higgs constraints. In contrast to the CP-conserving MSSM scenarios, here the couplings and mixing matrices can be complex, and the three Higgs bosons can also be mixed and have scalar and pseudoscalar components, leading to three states h 1,2,3 , named in the order of increasing masses.

Constraints in the CMSSM
We consider the CMSSM as a benchmark scenario to study the effects of the CP phases. For this purpose, we focus on the CP-conserving CMSSM best-fit point found in a global analysis [20]:  and vary the phases subsequently. To test the efficiency of the geometric approach, we vary the phases randomly, then employ the geometric approach and finally apply the constraints. 600.000 points are generated for each case. Imposing the constraints removes about 85% of the points in the purely random approach, while only 50% of the points are removed with the geometric approach, showing the efficiency of the latter. Fig. 1 shows the distribution of the phases obtained by imposing the EDM constraints. Clearly, the phases Φ 1 , Φ 2 , Φ 3 are strongly constrained. The EDM constraints also severely restrict the possible Higgs masses, as shown in Fig. 2. In addition, the EDM constraints impose the Higgs signal strengths to be very close to 1, as demonstrated in Fig. 3.
More generically, when varying all the CMSSM parameters, the EDM measurements strongly limit the CPviolating CMSSM to be very close to the CP-conserving CMSSM (see [2] for more details). This conclusion however does not hold in an unconstrained scenario such as the phenomenological MSSM (pMSSM).

Constraints in the pMSSM
To study the pMSSM, we first vary the 19 CPconserving pMSSM parameters. We then use the geometric approach to choose points with phases between −180 • and +180 • . 40 million such points have been generated with gluino and squark masses compatible with the LHC limits. We then impose the h 1 mass to lie between 121 and 129 GeV and the lightest neutralino to be the lightest supersymmetric particle so that it constitutes a dark matter candidate. This reduces the number of points to about one million. After imposing the EDM constraints, we are left with 15.000 points. The distribution of the CP phases is shown in Fig. 4.
As can be seen, the Φ 2 phase is particularly constrained, and the other phases can take any values. In Fig. 5, the phases of the points compatible with the EDM constraints are drawn in 2 dimensional plots. They reveal that the EDM constraints impose correlations between the phases related to the gluino and thirdgeneration squarks, with prefered directions that are revealed by the geometric approach.
In the following, we study the influence of CP viola-  tion on the dark matter, Higgs and flavour sectors.

Dark matter sector
We first compute the neutralino relic density and compare the results to the cosmological cold dark matter density Ω c h 2 ∼ 0.11 [21]. The results are presented in the upper panel of Fig. 6. We notice that the points with CP conservation and the ones with CP violation are similarly spread, such that the CP phases do not add much content to the picture. The points above the cosmological measurements have mostly bino-like neutralino, and the two strips below the dark matter line correspond to wino and higgsino-like neutralinos. The points above the line can be considered as excluded, since they lead to an overdensity of dark matter, provided the pre-Big Bang nucleosynthesis period is dominated by radiation [22,23], so that most of the bino-like neutralinos are excluded. CP violation allows for more such points to sur- vive, without changing the picture. For the points below the line the neutralino relic density could only partially account for the whole dark matter density, and they can therefore be considered as still allowed.
We also consider dark matter direct detection and the latest upper limits provided by the PandaX [24] and LUX [25] experiments. The results are shown in the lower panel of Fig. 6. Similarly to the case of the relic density, we see that while CP violation allows for some more spread of the points in comparison to the CP-conserving case, it does not add much to the picture. This shows that dark matter observables are unlikely to disentangle CP-violating scenarios form the CP-conserving ones.

Higgs sector
First, we assume the lightest Higgs boson to be the one discovered at the LHC by imposing its mass to lie between 121 and 129 GeV, so as to allow for theoretical uncertainties in its calculation, and we study the modifications to the signal strengths due to CP violation. The signal strength in the decay channel h 1 → XX is defined as the ratio of the cross section production of the h 1 times its decay branching ratio to XX over the Standard Model values. The LHC results for the measured decay modes show that the signal strengths are compatible with 1 [26]. Our results are summarised in Fig. 7. As can be seen, the set of points with CP violation are not distinct from the one with CP conservation, so that the signal strength measurement will not allow to probe the difference between the two cases.
The Higgs sector however provides other observables which can allow to probe CP violation. The Higgs bosons can be admixtures of scalar and pseudoscalar components. On the one side, the measurements of the h 1 properties showed that it is mostly a scalar, and the pseudoscalar component could only be tiny, and adding the EDM constraints limits the pseudoscalar component to be negligible. On the other side, the two heavier Higgses can be mixtures of scalar and pseudoscalar components. We denote g h if f S ,P the scalar and pseudoscalar couplings of h i tof f and define The angles φ h i τ and φ h i t characterise the scalar and pseudoscalar couplings to the τ leptons and top quarks, respectively. An angle close to 0 means a mostly scalar coupling, while π/2 means a pure pseudoscalar coupling. In Fig. 8, we show the distributions of the angles of the h 2 after imposing the EDM constraints. Similar plots can be obtained for h 3 [2]. These results demonstrate that strongly-mixed states are still allowed for the heavy Higgses, and that a measurement of the ττ and tt spin correlations in h 2,3 decays may reveal the CPviolating nature of the MSSM.

Flavour sector
CP violation is intimately related to flavour physics, in particular because of the CP-violating phase present in the CKM matrix. Flavour-changing neutral currents constitute important observables to probe for new physics and CP violation. In particular, the inclusive decay B → X s γ measured by B-factories provides valuable constraints. In Fig. 9, we present the constraints imposed by the measurements of the branching ratio and CP asymmetry of this decay on our pMSSM points passing the EDM constraints. For comparison, we show the current [27] and prospective [28] limits. As can be seen, the current limits are superseded by the EDM constraints, but the perspectives to discover CP violation in the MSSM at Belle-II are promising.
Another important flavour observable is the B s meson mixing ∆M B s . As can be seen in Fig. 10, the CPviolating MSSM contributions to ∆M NP B s are in general below the present upper limit, which is dominated by theoretical uncertainties. If these could be reduced, ∆M B s could also provide an interesting and complementary constraint on the phases, enabling them to be determined experimentally, in principle.

Conclusion
We have explored the effects of CP violation in the MSSM using an iterative geometrical approach within  Bs for the pMSSM points. The gray curve corresponds to the points before the EDM constraints, and the black curve after the EDM constraints. The red dashed line shows the current limit and the green dashed line the prospective limit after a factor 10 improvement in the form factor determination. the maximally CP-violating, minimal flavour-violating framework with 6 CP-violating phases, applying constraints from the Higgs mass and signal strengths, flavour physics, dark matter relic density and scattering cross section with matter.
In the CMSSM scenario, we found relatively little scope for large deviations from the CP-conserving case, e.g., in the masses of the Higgs bosons and the spinindependent dark matter scattering cross section. Moreover, we found that only very small values of A CP would be possible in this case, and the new physics contribution to B s meson mixing would not be observable.
In the pMSSM scenario, the A CP could be as large as ∼ 3%, within the reach of the Belle-II experiment. We found in this scenario that ∆M NP B s could be large enough to be observable with a prospective reduction in the theoretical uncertainty in the Standard Model calculation of B s mixing.
The CP-violating phases in the h 1 τ + τ − and h 1t t cou-plings are in general small in the MSSM. However, the phases in the h 2,3 τ + τ − and h 2,3t t couplings can be quite large, and may present interesting prospects for future experiments.
In general, we showed that the EDM constraints do not force all the six CP-violating phases to be small, and in some variants of the MSSM there could be observable signatures of CP violation beyond the Standard Model, such as A CP in the B → X s γ decay and ∆M NP B s .