Light and compressed gluinos at the LHC via string theory

TL;DR

This paper explores how string theory models, specifically large volume scenarios with D3-branes, predict light gluinos that could evade current LHC detection, linking theoretical models to experimental data.

Contribution

It demonstrates the use of global fits of string theory parameters to low energy data to identify potential signatures of light gluinos at the LHC.

Findings

LVS models predict light, quasi-degenerate gluinos and neutralinos

Such gluinos may evade current LHC search strategies

Global fits connect string theory parameters to observable physics

Abstract

In this article, we show that making global fits of string theory model parameters to data is an interesting mechanism for probing, mapping and forecasting connections of the theory to real world physics. We considered a large volume scenario (LVS) with D3-branes matter fields and supersymmetry breaking. A global fit of the parameters to low energy data shows that the set of LVS models are associated with light gluinos which are quasi-degenerate with the neutralinos and charginos they can promptly decay into and thus possibly hidden to current LHC gluino search strategies.

Tables

Table 1. Table 1: Summary for the central values and errors for the electroweak precision, B-physics and cold dark matter relic density constraints.

Observable	Constraint	Observable	Constraint
$m_W$ [GeV]	$80.399\pm 0.023$ verzo	$A^l=A^e$	$0.1513\pm 0.0021$ :2005ema
${\Gamma }_Z$ [GeV]	$2.4952\pm 0.0023$ :2005ema	$A^b$	$0.923\pm 0.020$ :2005ema
${\sin }^2{\theta }_{eff}^{lep}$	$0.2324\pm 0.0012$ :2005ema	$A^c$	$0.670\pm 0.027$ :2005ema
$R_l^0$	$20.767\pm 0.025$ :2005ema	$Br\left(B_s\to {\mu }^{+}{\mu }^{-}\right)$	${3.2}_{-1.2}^{+1.5}\times {10}^{-9}$ Aaij:2012nna
$R_b^0$	$0.21629\pm 0.00066$ :2005ema	$\Delta M_{B_s}$	$17.77\pm 0.12$ ps-1 Abulencia:2006ze
$R_c^0$	$0.1721\pm 0.0030$ :2005ema	$R_{Br\left(B_u\to \tau \nu \right)}$	$1.49\pm 0.3091$ Aubert:2004kz
$A_{\text{FB}}^b$	$0.0992\pm 0.0016$ :2005ema	$\Delta M_{B_d}$	$0.507\pm 0.005$ ps-1Barberio:2008fa
$A_{\text{FB}}^c$	$0.0707\pm 0.0035$ :2005ema	${\Omega }_{CDM}{h}^2$	$0.11\pm 0.02$ 0803.0547
$m_h$ [GeV]	$125.6\pm 3.0$ 1207.7214 ; 1207.7235	$Br\left(B\to X_s\gamma \right)$	$\left(3.52\pm 0.25\right)\times {10}^{-4}$ Barberio:2007cr

Equations

$$V$$

$$D_{i}W$$

$$G^{i\bar{j}}$$

$$K=K_0+K_{1}X+{\bar{K}}_1\bar{X}+K_{2}X\bar{X},W=\rho X+W_0.$$

$$K=-2\log \left(\mathcal{V}-\hat{\xi }\right)+{\tilde{K}}_i\varphi \bar{\varphi }+{\tilde{Z}}_{i}X\bar{X}+{\tilde{H}}_i\varphi \bar{\varphi }X\bar{X}+....$$

$${\tilde{K}}_i=\frac{{\alpha }_0}{{\mathcal{V}}^{2/3}}\left(1-{\alpha }_1\frac{\xi s^{3/2}}{\mathcal{V}}\right),{\tilde{Z}}_i=\frac{{\beta }_0}{{\mathcal{V}}^{2/3}}\left(1-{\beta }_1\frac{\xi s^{3/2}}{\mathcal{V}}\right),{\tilde{H}}_i=\frac{{\gamma }_0}{{\mathcal{V}}^{4/3}}\left(1-{\gamma }_1\frac{\xi s^{3/2}}{\mathcal{V}}\right).$$

$$W=W_0+\rho X+A{e}^{-a_s{T}_s}$$

$$m_0^2=\frac{5}{4}\frac{s^{3/2}\xi }{\mathcal{V}}(3{\alpha }_1-1)m_{3/2}^2+\frac{9}{8}\frac{s^{3/2}\xi }{\mathcal{V}}\frac{1}{5a_s{\tau }_s}\left(1-\frac{3{\gamma }_0}{{\alpha }_0{\beta }_0}\right)m_{3/2}^2,$$

$$M_{1/2}=sign(W_0)\frac{3}{4}\frac{s^{3/2}\xi }{\mathcal{V}}\left[3-2{\omega }_s\right]m_{3/2},\text{ and }$$

$$A_0=-(1-y)M_{1/2},y=s{\partial }_s\log Y_{ijk}^{(0)}.$$

$$\underline{\theta }\equiv \{x,{\alpha }_1,y,{\omega }_s,\tan \beta ,{\log }_{10}{m}_{3/2},m_Z,m_t,m_b,{\alpha }_{em},{\alpha }_s\}.$$

$$\underline{O}$$

$$p(\underline{d}|\underline{\theta })=L(x)\prod _i\frac{\exp \left[-(O_i-{\mu }_i{)}^2/2{\sigma }_i^2\right]}{\sqrt{2\pi {\sigma }_i^2}}$$

$$L(x)=\{\begin{array}{ll}1/(y+\sqrt{\pi s^2/2})& \text{ if }x<y\\ \exp \left[-(x-y{)}^2/2s^2\right]/(y+\sqrt{\pi s^2/2})& \text{ if }x\ge y\end{array}$$

$$p(\underline{\theta }|\underline{d})=\frac{p(\underline{d}|\underline{\theta })\times p(\underline{\theta })}{p(\underline{d})}$$

$$\tau \approx 8{\left(\frac{\tilde{m}}{10^9\text{ GeV}}\right)}^4{\left(\frac{1\text{ TeV}}{m_{\tilde{g}}}\right)}^{5}s$$

Full text

1 Introduction

Flux compactification with stabilised moduli sets a connection of string theory models to real world physics hep-th/0610102 ; hep-th/0701050 . The parameters describing the various sets of string theory models lead to different physical properties. This is a feature associated to the so-called “landscape” of string theory vacua phenomenon. The KKLT hep-th/0301240 ; hep-th/0503216 and large volume scenarios hep-th/0502058 ; hep-th/0610129 are famous examples with each representing sets of string theory models. Within each of these, there are additional model variations including differences in gauge groups, particle representations and the cosmological constants for representing the physical world. Our target in this article is about relating the sets of models’ predictions to collider observables and, eventually, to collider phenomenology.

The vacuum that we inhabit should be just one from the set of landscape vacua possibilities. An interesting question goes as follows. Is it possible to have an explicit set of string theory based models which is in complete agreement with our real world? Here real world can be taken to mean the standard models of particle physics and cosmology, or even including their extensions such as with supersymmetry breaking. For the later case, if one assumes that supersymmetry exists in nature but broken at some energy scale accessible to colliders, then predictions for supersymmetry-breaking from flux compactification models can be used for assessing the string theory models’ parameters 111Throughout this article, we assume that it is possible to find a compactification with exactly the MSSM matter content. See Braun:2005nv ; Cicoli:2013mpa ; Cicoli:2012vw ; Romao:2015jrh for instance and develoments along this direction.. This is also true for cosmological observables that can be related to moduli fields. But in this article, the focus is on the particle physics aspect.

Sampling string theory based model parameters for making global fit to certain experimental data, representing the real world, could be the best way for finding explicit models in agreement with observations. For this purpose, here, the Bayesian techniques for fitting models to data AbdusSalam:2008uv ; AbdusSalam:2009qd were employed and applied to the string theory models. Based on the experimental results used, statistical weights will be assigned to points for mapping and making forecasts about the parameter space. This approach is new within string theory phenomenological research (see 1612.01569 for an overview).

For instance, consider a IIB string theory model with the dilaton and complex structure moduli stabilised via flux superpotential $W$ hep-th/9906070 ; Giddings:2001yu . Its vacuum expectation value (vev), $W_{0}$, affects the physical properties of the string vacuum such as the magnitude of the cosmological constant and scale of supersymmetry breaking. A scan over the landscape of flux compactifications can be achieved by varying the value of the flux superpotential $W_{0}$. In AbdusSalam:2007pm , by considering a particular Calabi-Yau manifold with two moduli fields and starting from the large volume scenario (LVS) hep-th/0502058 ; hep-th/0610129 limit, with $W_{0}\sim 1$, or the KKLT limit, with $W_{0}\ll 1$ hep-th/0301240 ; hep-th/0503216 a set of flux-dependent AdS and dS vacua (without uplift terms) were obtained. Some of the minima have supersymmetry spontaneously broken by the fluxes for matter fields on D7 branes. One can go further by computing the full-fledged supersymmetry spectrum via RG-running of the predicted supersymmetry-breaking terms to the electroweak scale such as in 0704.3403 . Given the supersymmetry spectrum, other observables such the Higgs boson mass, dark matter relic density, electroweak precision and B-physics observables can be computed. A non-agreement between the predicted observables and their corresponding experimentally measured values or limits can be used for marking the point $W_{0}$ as not possible for representing our world. This way, statistical methods for exploration and fitting the landscape parameters can be used for mapping and physics forecasting.

In this article, we consider a particular LVS model Aparicio:2015psl with matter fields on D3-brane and supersymmetry broken via a chiral superfield $X$ with nilpotent constraint ($X^{2}=0$) Kallosh:2014wsa ; Bergshoeff:2015jxa ; Kallosh:2015nia . In section 2 we briefly introduce the LVS and the supersymmetry-breaking terms derived in Aparicio:2015psl for setting the context of the article and defining an LVS subset of the minimal supersymmetric standard model (MSSM), with R-parity conserved, called LMSSM-6 as the phenomenological frame for our analysis. The Bayesian method for fitting the LMSSM-6 to low-energy physics is presented in section 3. The results of the global fit indicate an abundance of light gluinos, neutralinos and charginos together with ${\cal O}(10)TeV$ squarks. The electroweak-inos are within one to few top-quark heavy. The mass difference between the gluinos, charginos or heavy neutralinos and the lightest supersymmetric particle (LSP) turns out to be small thereby making a lower sub-TeV quasi-degenerate spectrum. As such the gluinos and electroweak-inos are possibly surviving current LHC limits. Detailed reinterpretation or recasting analyses should ultimately lead to a robust conclusion about the status of such spectrum at the LHC. In section 4 we round up the the article with conclusions and outlook.

2 Large volume scenario (LVS) SUSY-breaking, LMSSM-6

The shape and volume of string theory internal space represent massless scalar fields from a $4$-dimensional theory point of view. These moduli fields must be stabilised since no such extra scalar fields are observed in nature. There are many moduli stabilisation work in the literature. Here we concentrate on a LVS scenario where all the moduli fields are fixed at an exponentially large internal volume. The LVS supergravity effective theory couples to the particle content of the MSSM. In the effective theory, the Kähler and superpotentials are generated from the superstring theory. This forms the hidden sector which then couples via gravity mediation to the MSSM sector. 4d $N=1$ supergravity is specified up to two derivatives by the Kähler potential $K$, superpotential $W$ and gauge kinetic function $f_{a}$. With these, the scalar potential is given by

[TABLE]

Here $i,j$ run over the two Kähler moduli $T_{s}$ and $T_{b}$ whose real parts $\tau_{b}$ and $\tau_{s}$ determine the internal space volume, $\mathcal{V}=\tau_{b}^{3/2}-\tau_{s}^{3/2}$. At the minimum of the potential the moduli fields acquire vevs and provide non vanishing auxiliary field that spontaneously breaks supersymmetry and generate the breaking terms in the visible sector. Other dynamics such as the presence of D3-brane within the superstring construction can also lead to the breaking of supersymmetry.

The LVS supersymmetry-breaking terms induced by the presence of nilpotent superfield $X$ Aparicio:2015psl is presented as follows. The Kähler potential and superpotentials for $X$ can be taken as

[TABLE]

where $K_{0},K_{1},K_{2},\rho,W_{0}$ are coefficients and can be functions of other low-energy fields. Assuming the dilaton and complex structure moduli have been fixed and integrated out at high scale, the LVS Kähler potential becomes

[TABLE]

Here $\hat{\xi}=\frac{s^{3/2}}{\xi}{2}$, $\xi$ is a Calabi-Yau manifold related constant of order one Becker:2002nn . $s=1/g_{s}$ is the real-part of the axion-dilaton field. $\phi$ represents matter, $\tilde{K}_{i}$ and $\tilde{Z}_{i}$ are respectively the D3-branes matter and the nilpotent goldstino metrics. $\tilde{H}_{i}$ represents the nilpotent goldstino and matter fields quartic interaction. These are parameterised for the $\alpha^{{}^{\prime}}$-corrected potential as

[TABLE]

The LVS superpotential is given by

[TABLE]

where $A$, and $a_{s}$ are gaugino condensation parameters. The Kähler and superpotentials can then be used for computing the soft supersymmetry-breaking terms via a standard method as done in Aparicio:2015psl . The scalar $m_{0}$, gaugino $M_{1/2}$, and trilinear coupling $A_{0}$ soft supersymmetry-breaking terms for visible sector fields on D3-brane are

[TABLE]

[TABLE]

[TABLE]

Here $m_{3/2}=e^{K/2}|W|$ is the gravitino mass. $\omega_{s}\lesssim 1$ parametrise corrections such that $D_{s}W\sim\frac{e^{-a\,\tau}}{2s}\omega_{s}$ used in computing the soft supersymmetry-breaking terms.

Given a set of LMSSM-6 parameters, the supersymmetry-breaking terms Eq.(8) to Eq.(10) can be computed. These are then set as the boundary conditions for renormalisation group (RG) running from the supersymmetry-breaking scale, $m_{3/2}$, to the weak-scale. For doing this we chose a particular specialisation of gravity-mediated supersymmetry-breaking with a non-universal scalar mass terms for the MSSM Higgs doublets such that $m_{H_{1}}=m_{H_{2}}=0$ at symmetry-breaking scale. Other possibilities include varying the non-universal Higgs doublet mass terms and the minimal supergravity but are not considered here. The set of LVS parameters, $\{\,x,\,\alpha_{1},\,y,\,\omega_{s},\,\tan\beta,\,\log_{10}m_{3/2}\}$, in the soft supersymmetry-breaking terms above together the RG to the weak scale with MSSM sparticle content is referred to as the LVS MSSM with six parameters (LMSSM-6). $\tan\beta$ is the ratio of the Higgs doublet vevs. For making a global fit of the LMSSM-6 parameters to data we add 5 Standard Model nuisance parameters that are for some of the precision observables used. As such we will be exploring a total of 11 parameters

[TABLE]

The Standard Model nuisance parameters are the Z-boson mass $m_{Z}=91.1876\pm 0.0021$ GeV, the top quark mass $m_{t}=172.6\pm 1.4$ GeV, the bottom quark mass $m_{b}=4.2\pm 0.07$ GeV, the electromagnetic coupling constant $\alpha_{em}^{-1}=127.918\pm 0.018$, and the strong interactions coupling constant $\alpha_{s}=0.1172\pm 0.002$. These were all set to vary in a Gaussian manner with central values and deviations according to the experimental results. The LVS parameters are were allowed to vary as: $x\equiv\frac{3\gamma_{0}}{\alpha_{0}\beta_{0}}\in[0.01,100.0]$, $\alpha_{1}\in[0.01,100.0]$, $y\in[-1000,1000]$, $\omega_{s}\in[0.01,100.0]$, $\tan\beta\in[2,60]$, $log_{10}m_{3/2}\in[3,19]$. The other LVS parameters which were not varied but enter the soft terms are $\xi=1.0,W_{0}=\pm 1.0,s=25.0$, and ${\mathcal{V}}=1.0\times 10^{9}$ units. In the next section we describe the procedure for fitting the parameters to low energy data.

3 The LMSSM-6 fit to low energy physics

We consider a flat (Gaussian) prior probability density distribution, $p(\underline{\theta})$, for the LMSSM-6 (Standard Model) parameters in Eq.(11) for exploring the space using the MultiNest Feroz:2007kg ; Feroz:2008xx package which implements Nested Sampling algorithm Skilling . During the model exploration, at each parameter space point the following packages were used via the SLHA1 Skands:2003cj interface: the supersymmetry spectrum generator SOFTSUSY Allanach:2001kg for computing the sparticle masses, mixing angles and couplings, micrOMEGAs Belanger:2008sj for computing the neutralino cold dark matter relic density and the anomalous magnetic moment of the muon $\delta a_{\mu}$, SuperIso Mahmoudi:2007vz for predicting the branching ratios $BR(B_{s}\rightarrow\mu^{+}\mu^{-})$, $BR(B\rightarrow s\gamma)$ and the isospin asymmetry, $\Delta_{0-}$ in the decay $B\rightarrow K^{*}\gamma$, and susyPOPE Heinemeyer:2006px ; Heinemeyer:2007bw for computing precision observables that include the $W$-boson mass $m_{W}$, the effective leptonic mixing angle variable $\sin^{2}\theta^{lep}_{eff}$, the total $Z$-boson decay width, $\Gamma_{Z}$, and the other electroweak observables whose experimentally determined central values and associated errors are summarised in Table 1. The experimental central values ($\mu_{i}$) and errors ($\sigma_{i}$) for these make the set of the predictable observables, $\underline{O}$:

[TABLE]

The compatibility of the LMSSM-6 parameter-space points with the data is quantified by the likelihood, $p(\underline{d}|\underline{\theta})$. Assuming the observables are independent, the combined likelihood can calculated as

[TABLE]

where the index $i$ runs over the list of observables $\underline{O}$. Here, $x$ represents the predicted value of neutralino cold dark matter (CDM) relic density and

[TABLE]

where $y=0.11$ is the CDM relic density central value and $s=0.02$ the corresponding inflated (to allow for theoretical uncertainties) error.

The outcome of the Bayesian fit of the LMSSM-6 to data is the posterior probability density

[TABLE]

shows that the sfermions are order tens of TeV heavy while the gluino and electroweak-inos are light, in few hundreds of GeV range. The posterior probability distributions for the LMSSM-6 supersymmetry-breaking parameters, the Higgs boson mass for reference purpose, the electroweak-inos and the gluino-neutralino mass difference are shown in Figure 1.

The heavy squarks can potentially make the gluinos long-lived or to turn into R-hadrons before decaying. In Fig. 2, the 2-dimensional posterior distribution for gluino decay length estimate versus its mass is given. The gluino decay life-time for squark at order 10 TeV is estimated as Dawson:1983fw ; Hewett:2004nw

[TABLE]

It turns-out that the gluino and electroweak-inos are compressed relative to the LSP mass. For this reason, one cannot automatically rule out LMSSM-6 light gluinos.

The direct search for gluinos at the LHC puts the lower limit on gluino mass well into the TeV towards multi-TeV region (see for instance refs 1507.05525 ; SUS-16-014 ; SUS-16-015 .) The limits come from search channels with multijets plus zero leptons and can probe gluinos decaying into gluons and neutralinos. The same is the case for the gluino decay to two quarks and a neutralino; and for the topologies with initial state radiations. However it is not clear what the status will be for the light gluinos which are quasi-degenerate with the neutralinos or charginos they can decay to. Gluino production cross-section at the LHC is much higher than squarks’ of the same mass. So in principle as long as the gluino decays to gluon and neutralino or to a quark pair and neutralino then then, for instance, the CMS limits from searches with multi-jet final states SUS-16-014 ; SUS-16-015 should apply for the LMSSM-6 gluinos. But, that and the application of other relevant LHC limits necessarily require rigorous and dedicated reinterpretation analysis of the results.

In 1701.07664 various experimental search strategies for finding gluinos quasi-degenerate with neutralinos have been addressed. These include searches for displaced vertices and disappearing tracks. The analysis was however based on a fixed gluino mass at 1.5 TeV. It will be interesting to see what the outcome of a similar analysis will be for light gluinos such as for the LMSSM-6 shown in this article. We again emphasis that a decisive conclusion requires detailed reinterpretation of the experimental results within the LMSSM-6 context.

4 Conclusion and outlook

Models from string theory compactification with moduli stabilised by fluxes can predict supersymmetry-breaking at the TeV scale. The generation of supersymmetry spectrum allows the connection of the string theory or landscape parameters to experimental observables for constraining the the theory. In this article, we have applied Bayesian global fit technique for string theory inspired phenomenology. We introduce LMSSM-6 for representing the set of supersymmetry-breaking parameters predicted from a large volume scenario with matter fields on D3-branes. The field content were considered to be that of the minimal supersymmetric standard model with R-parity conserved and neutralino lightest supersymmetric particle as dark matter candidate. This way, the LMSSM-6 parameters which originate from the string theory setting can be constrained via the low energy properties of the associated sparticle spectrum. This feature was used for forecasting and mapping the string theory based parameter space.

The low energy constrains used are the Higgs boson mass, dark matter relic density, electroweak precision and B-physics experimentally measured observables. The global fit of the LMSSM-6 parameters to these showed that the bulk of the posterior distribution yielded heavy squarks at order 10 TeV. These come together with light (much less than or of the order of the top-quark mass) but promptly decaying gluinos. Such gluinos would have been already ruled out even before the LHC if not for being quasi-degenerate with neutralinos and charginos. A dedicated reinterpretation of the experimental results that probe gluino productions at the LHC such as in 1507.05525 ; SUS-16-014 ; SUS-16-015 is needed for determining the status of the gluinos. This is an interesting work to be done but is beyond the scope of this article. Should the LMSSM-6 gluinos be ruled out, then the corresponding class of large volume scenario models cannot represent our real world.

As an outlook, the following are interesting questions or work to be done based on the concepts addressed in this article.

• •

The methodology presented can be applied to other string theory phenomenology frames such as the famous KKLT hep-th/0301240 ; hep-th/0503216 scenarios in comparison to the large volume ones.

• •

The posterior sample from the global fit of parameters to data can be used in recasting collider results for establishing the status of the considered models (such as the LMSSM-6 presented here).

• •

There is an interesting complementarity between the supersymmetry spectrum such as from the LMSSM-6 global fit to data presented in this article, and similar spectra that could be constructed via the simplified models approach to supersymmetry phenomenology. With the global fit approach, the resulting spectra are guaranteed to be in agreement with the experimentally measured valued of the observables used for the fitting procedure. This is not necessarily the case for simplified model spectra. It is interesting to explore further and contrast this and similar characteristics.

Acknowledgment

Thanks to M.Pierini and F.Quevedo for discussions; and ICTP for hospitality during an early stage of the research presented here. An early stage of this research was supported by an ERC fund under the European Union’s Seventh Framework Programme (FP/2007-2013)/ERC Grant Agreement no. 279972 NPFlavour.