# Multi-TeV Scalars are Natural in Minimal Supergravity

IASSNS-HEP-99-78

FERMILAB-PUB-99/226-T

For a top quark mass fixed to its measured value, we find natural regions of minimal supergravity parameter space where all squarks, sleptons, and heavy Higgs scalars have masses far above 1 TeV and are possibly beyond the reach of the Large Hadron Collider at CERN. This result is simply understood in terms of “focus point” renormalization group behavior and holds in any supergravity theory with a universal scalar mass that is large relative to other supersymmetry breaking parameters. We highlight the importance of the choice of fundamental parameters for this conclusion and for naturalness discussions in general.

The standard model with a fundamental Higgs boson suffers from a large and unexplained hierarchy between the weak and Planck scales . Because supersymmetric theories are free of quadratic divergences, however, this hierarchy is stabilized in supersymmetric extensions of the standard model when the scale of superpartner masses is roughly of order the weak scale $`M_{\mathrm{weak}}`$ . The promise of providing a natural solution to the gauge hierarchy problem is the primary phenomenological motivation for supersymmetry.

Because the requirement of naturalness places upper bounds on superpartner masses, this criterion has important experimental implications. In a model-independent analysis, naturalness constraints are weak for some superpartners, e.g., the squarks and sleptons of the first two generations . However, in widely studied scenarios where the scalar masses are unified at some high scale, such as minimal supergravity, it is commonly assumed that squark and slepton masses must all be $`\stackrel{<}{}1\mathrm{TeV}`$. This bound places all scalar superpartners within the reach of present and near future colliders, and is a source of optimism in the search for supersymmetry at the high energy and high precision frontiers. We show here, however, that this assumption is invalid, and in fact, it is precisely in supergravity theories with a universal scalar mass that all squark and slepton masses may naturally be far above 1 TeV.

Supersymmetric theories are considered natural if the weak scale is not unusually sensitive to small variations in the fundamental parameters. Although the criterion of naturalness is inherently subjective, its importance for supersymmetry has motivated several groups to provide quantitative definitions of naturalness . In this analysis, we adopt the following prescription:

(1) We consider the minimal supergravity framework with its 4+1 input parameters

$$\left\{P_{\mathrm{input}}\right\}=\{m_0,M_{1/2},A_0,\mathrm{tan}\beta ,\mathrm{sign}(\mu )\},$$
(1)   
where $`m_0`$, $`M_{1/2}`$, and $`A_0`$ are the universal scalar mass, gaugino mass, and trilinear coupling, respectively, $`\mathrm{tan}\beta =H_u^0/H_d^0`$ is the ratio of Higgs expectation values, and $`\mu `$ is the Higgsino mass parameter. The first three parameters are at the grand unified theory (GUT) scale $`M_{\mathrm{GUT}}2\times 10^{16}\mathrm{GeV}`$, i.e., the scale where the U(1)<sub>Y</sub> and SU(2) coupling constants meet.

(2) The naturalness of each point $`𝒫\{P_{\mathrm{input}}\}`$ is then calculated by first determining all the parameters of the theory (Yukawa couplings, soft supersymmetry breaking masses, etc.), consistent with low energy constraints. Renormalization group (RG) equations are used to relate high and low energy boundary conditions. In particular, at the weak scale, proper electroweak symmetry breaking requires<sup>1</sup><sup>1</sup>1The tree-level conditions are displayed here for clarity of presentation. In all numerical results presented below, we use the full one-loop Higgs potential , minimized at the scale $`m_0/2`$, approximately where one-loop corrections are smallest, as well as two-loop RG equations , including all low-energy thresholds .

$`{\displaystyle \frac{1}{2}}m_Z^2`$ $`=`$ $`{\displaystyle \frac{m_{H_d}^2m_{H_u}^2\mathrm{tan}^2\beta }{\mathrm{tan}^2\beta 1}}\mu ^2`$ (2)   
$``$ $`f(m_{H_d}^2,m_{H_u}^2,\mathrm{tan}\beta )\mu ^2,`$   
$`2B\mu `$ $`=`$ $`\mathrm{sin}2\beta (m_{H_d}^2+m_{H_u}^2+2\mu ^2),`$ (3)   
where $`m_{H_u}^2`$ and $`m_{H_d}^2`$ are the soft scalar Higgs masses, and $`B\mu `$ is the bilinear scalar Higgs coupling.

(3) We choose to consider the following set of (GUT scale) parameters to be free, independent, and fundamental:

$$\{a_i\}=\{m_0,M_{1/2},A_0,B_0,\mu _0\}.$$
(4)   
(4) All observables, including the $`Z`$ boson mass, are then reinterpreted as functions of the fundamental parameters $`a_i`$, and the sensitivity of the weak scale to small fractional variations in these parameters is measured by the sensitivity coefficients 

$$c_i\left|\frac{\mathrm{ln}m_Z^2}{\mathrm{ln}a_i}\right|.$$
(5)   
(5) Finally, we form the fine-tuning parameter

$$c=\mathrm{max}\{c_i\},$$
(6)   
which is taken as a measure of the naturalness of point $`𝒫`$, with large $`c`$ corresponding to large fine-tuning.

As is clear from the description above, several subjective choices have been made, as they must be in any definition of naturalness. The choice of minimal supergravity in step (1), and particularly the assumption of a universal scalar mass, plays a crucial role. Deviations from this assumption will be considered below.

The choice of fundamental parameters in step (3) is also important and varies throughout the literature. An appealingly simple choice (see, e.g., Ref. ) is $`\{a_i\}=\{\mu \}`$, where $`\mu `$ is to be evaluated at the weak scale. This is equivalent to using $`\mu ^2`$ as a fine-tuning measure, since Eqs. (2) and (5) imply $`c_\mu =4\mu ^2/m_Z^2`$. While generally adequate, this definition is insensitive to large fine-tunings in the function $`f`$ of Eq. (2), as we will see below; such fine-tunings are accounted for in the more sophisticated choice of Eq. (4).

The top quark Yukawa $`Y_t`$ (sometimes along with other standard model parameters, such as the strong coupling) is included among the fundamental parameters in some studies  and not in others . This choice typically attracts little comment, and attitudes toward it are at best ambivalent . This ambiguity reflects, perhaps, a diversity of prejudices concerning the fundamental theory of flavor. It is important to note, however, that unlike the parameters of Eq. (4), $`Y_t`$ is not expected to be related to supersymmetry breaking and is, in some sense, now measured, as it is strongly correlated with the top quark mass $`m_t`$. For these reasons, we find it reasonable to assume that in some more fundamental theory, $`Y_t`$ is fixed to its measured value in a flavor sector separate from the supersymmetry breaking sector, and we therefore do not include it among the $`a_i`$. This choice is critical for our conclusions, as will be discussed below.

In step (5), various other choices are also possible. For example, the $`c_i`$ may be combined linearly or in quadrature; we follow the most popular convention. In other prescriptions, the $`c_i`$ are combined after first dividing them by some suitably defined average $`\overline{c}_i`$ to remove artificial appearances of fine-tuning . We have not done this, but note that such a normalization procedure typically reduces the fine-tuning measure and would only strengthen our conclusions.

Given the prescription for measuring naturalness described above, we may now present our results. In Fig. 1, contours of constant $`c`$, along with squark mass contours, are presented for $`\mathrm{tan}\beta =10`$. Moving from low to high $`m_0`$, the contours are determined successively by $`c_{\mu _0}`$, $`c_{M_{1/2}}`$ and $`c_{m_0}`$. The naturalness requirement $`c<25`$ ($`c<50`$) allows regions of parameter space with $`m_02\mathrm{TeV}(2.4\mathrm{TeV})`$. More importantly, regions with $`m_0\stackrel{>}{}2\mathrm{TeV}`$, where all squarks and sleptons have masses well above 1 TeV, are as natural as the region with $`(m_0,M_{1/2})\stackrel{<}{}(1000\mathrm{GeV},400\mathrm{GeV})`$, where squark masses are below 1 TeV.

The naturalness of multi-TeV $`m_0`$, though perhaps surprising, may be simply understood as a consequence of a “focus point” in the RG behavior of $`m_{H_u}^2`$ , which renders its value at $`M_{\mathrm{weak}}`$ highly insensitive to its value in the ultraviolet. Note that for moderate and large $`\mathrm{tan}\beta `$, Eq. (2) implies that $`m_Z^2`$ is insensitive to $`m_{H_d}^2`$ and is determined primarily by $`m_{H_u}^2`$.

Consider any set of minimal supergravity input parameters. These generate a particular set of RG trajectories, $`m_i^2|_\mathrm{p}(t),M_i|_\mathrm{p}(t),A_i|_\mathrm{p}(t),\mathrm{}`$, where $`t\mathrm{ln}(Q/M_{\mathrm{GUT}})`$ and $`Q`$ is the renormalization scale. Now consider another set of boundary conditions that differs from the first by shifts in the scalar masses. The new scalar masses $`m_i^2=m_i^2|_\mathrm{p}+\delta m_i^2`$ satisfy the RG equations

$$\frac{d}{dt}m_i^2\frac{1}{16\pi ^2}\left[g^2M_{1/2}^2+Y^2A^2+\underset{j}{}Y^2m_j^2\right]$$
(7)   
at one-loop, where positive numerical coefficients have been omitted, and the sum is over all chiral fields $`\varphi _j`$ interacting with $`\varphi _i`$ through the Yukawa coupling $`Y`$. However, because the $`m_i^2|_\mathrm{p}`$ are already a particular solution to these RG equations, the deviations $`\delta m_i^2`$ obey the homogeneous equations

$$\frac{d}{dt}\delta m_i^2\frac{1}{16\pi ^2}\underset{j}{}Y^2\delta m_j^2.$$
(8)   
Such equations are easily solved. Assume for the moment that the only large Yukawa coupling is $`Y_t`$, i.e., $`\mathrm{tan}\beta `$ is not extremely large. Then $`\delta m_{H_u}^2`$ is determined from

$$\frac{d}{dt}\left[\begin{array}{c}\delta m_{H_u}^2\\ \delta m_{U_3}^2\\ \delta m_{Q_3}^2\end{array}\right]=\frac{Y_t^2}{8\pi ^2}\left[\begin{array}{ccc}3& 3& 3\\ 2& 2& 2\\ 1& 1& 1\end{array}\right]\left[\begin{array}{c}\delta m_{H_u}^2\\ \delta m_{U_3}^2\\ \delta m_{Q_3}^2\end{array}\right],$$
(9)   
where $`Q_3`$ and $`U_3`$ denote the third generation squark SU(2) doublet and up-type singlet representations, respectively. The solution corresponding to the universal initial condition $`\delta m_0^2(1,1,1)^T`$ is

$$[\begin{array}{c}\delta m_{H_u}^2\\ \delta m_{U_3}^2\\ \delta m_{Q_3}^2\end{array}]=\frac{\delta m_0^2}{2}\left\{[\begin{array}{c}3\\ 2\\ 1\end{array}]\mathrm{exp}\left[_0^t\frac{6Y_t^2}{8\pi ^2}𝑑t^{}\right][\begin{array}{c}\hfill 1\\ \hfill 0\\ \hfill 1\end{array}]\right\}.$$
(10)   
For $`t`$ and $`Y_t`$ such that $`\mathrm{exp}\left[\frac{6}{8\pi ^2}_0^tY_t^2𝑑t^{}\right]=1/3`$, $`\delta m_{H_u}^2=0`$, i.e., $`m_{H_u}^2`$ is independent of $`\delta m_0^2`$.

The RG evolution of $`m_{H_u}^2`$ in minimal supergravity is shown for several values of $`m_0`$ in Fig. 2. As expected, the RG curves exhibit a focus (not a fixed) point, where $`m_{H_u}^2`$ is independent of its ultraviolet value. Remarkably, however, for the physical top mass of $`m_t175\mathrm{GeV}`$, the focus point is very near the weak scale. Thus, the weak scale value of $`m_{H_u}^2`$ and, with it, the fine-tuning parameter $`c`$ are highly insensitive to $`m_0`$. If the particular solution is natural (say, with all input parameters near the weak scale), the new solution, even with very large $`m_0`$, is also natural.

We have also checked numerically that the focusing effect persists even for very large values of $`\mathrm{tan}\beta `$. Indeed, in the limit $`Y_t=Y_bY_\tau `$, Eq. (8) can be similarly solved analytically, and one finds that focusing occurs for $`\mathrm{exp}\left[\frac{7}{8\pi ^2}_0^tY_t^2𝑑t^{}\right]=2/9`$. For the experimentally preferred range of top masses, the focus point is again tantalizingly close to $`M_{\mathrm{weak}}`$ .

The naturalness of multi-TeV $`m_0`$ has important implications for collider searches. Although $`m_{H_u}^2`$ is focused to the weak scale, all other soft masses remain of order $`m_0`$. From Eqs. (8) and (10), we find that for $`m_0M_{1/2},A_0`$, the physical masses of squarks, sleptons, and heavy Higgs scalars are well-approximated by

$`\stackrel{~}{t}_R:\sqrt{1/3}m_0`$ $`\mathrm{All}\mathrm{other}\stackrel{~}{\mathrm{q}},\stackrel{~}{\mathrm{}}:\mathrm{m}_0`$   
$`\stackrel{~}{t}_L,\stackrel{~}{b}_L:\sqrt{2/3}m_0`$ $`H^\pm ,A,H^0:m_0.`$ (11)   
Exact values of $`m_{\stackrel{~}{u}_L}`$ are presented in Fig. 1. All squarks, sleptons, and heavy Higgs scalars may therefore have masses $`\stackrel{>}{}12\mathrm{TeV}`$, and may be beyond the reach of the Large Hadron Collider (LHC) and proposed linear colliders. The discovery of such heavy scalars then requires some even more energetic facility, such as the envisioned muon or very large hadron colliders.

As may be seen from Fig. 1, however, fine-tuning constraints do not allow multi-TeV $`M_{1/2}`$. A similar conclusion applies to $`\mu `$, as may be seen in Fig. 3. We therefore expect all gauginos and Higgsinos to be within the kinematic reach of the LHC. Note that some regions of low $`\mu `$ are unnatural. In these regions, large cancellations in the function $`f`$ of Eq. (2) occur, and the simple definition $`c\mu ^2`$ is inadequate.

In addition to the gauginos and Higgsinos, the lightest Higgs boson is, of course, still required to be light. Contours of lightest Higgs mass $`m_h`$ are also presented in Fig. 3. Very heavy top and bottom squarks increase $`m_h`$ through radiative corrections: for low $`M_{1/2}`$, $`m_h`$ increases by roughly $`6`$ GeV as $`m_0`$ increases from 500 GeV to 2 TeV. However, in the multi-TeV $`m_0`$ scenario, naturalness requires $`A_0M_{\mathrm{weak}}`$ (see below), and so left-right squark mixing is suppressed. The upper bound on $`m_h`$ in Fig. 3 is thus approximately $`120\mathrm{GeV}`$, well below limits achieved for TeV squarks with maximal left-right mixing, and within the 3-5$`\sigma `$ discovery range of Higgs searches at the Tevatron with luminosity $`1030\mathrm{fb}^1`$ .

The focus point analysis presented above (for small $`Y_b`$) relied heavily on the universality of the $`H_u`$, $`U_3`$ and $`Q_3`$ soft masses. It is not hard to show, however, that GUT scale boundary conditions of the form $`(m_{H_u}^2,m_{U_3}^2,m_{Q_3}^2)=(1,1x,1+x)`$, for any $`x`$, also exhibit the focus point behavior. With respect to the other supersymmetry breaking parameters, the focus point is fairly robust. The mechanism is independent of all other scalar masses. Also, in the analysis above, any natural particular solution would do. Arbitrary and non-universal gaugino masses and trilinear couplings of order $`M_{\mathrm{weak}}`$ are therefore allowed. (Similarly, deviations in $`m_{H_u}^2`$, $`m_{U_3}^2`$, and $`m_{Q_3}^2`$ of order $`M_{\mathrm{weak}}^2`$ do not destabilize the focus point.) Note, however, that multi-TeV gaugino masses and $`A`$ parameters are not allowed. The required hierarchy between the scalar masses and the gaugino mass, $`A`$, and $`\mu `$ parameters may result from an approximate U(1)<sub>R+PQ</sub> symmetry or from the absence of singlet $`F`$ terms . $`B\mu `$ may also be suppressed by such a symmetry, and so leads to an experimentally viable scenario with naturally large $`\mathrm{tan}\beta m_{H_d}^2/(B\mu )`$, which is typically difficult to realize .

Although the focus point mechanism depends on a relation between $`m_t`$ and $`\mathrm{ln}(\frac{M_{\mathrm{GUT}}}{M_{\mathrm{weak}}})`$, it is not extraordinarily sensitive to these values. The focus point is still near the weak scale if $`m_t`$ is varied within its experimental uncertainty of 5 GeV, and, in fact, natural regions with multi-TeV $`m_0`$ are also possible if the high scale is raised to $`10^{18}\mathrm{GeV}`$ .

We stress, however, that if $`Y_t`$ is included among the free and fundamental parameters, multi-TeV $`m_0`$ would be considered unnatural. For example, for $`\mathrm{tan}\beta =10`$ and $`A_0=0`$, $`c_{Y_t}<25`$ (50) corresponds to $`m_0\stackrel{<}{}500\mathrm{GeV}`$ (800 GeV) . We have presented above our rationale for not including $`Y_t`$ among the $`a_i`$, although a definitive resolution of this issue most likely requires an understanding of the fundamental theory of flavor.

In conclusion, for moderate and large $`\mathrm{tan}\beta `$, multi-TeV scalars are natural in minimal supergravity. In view of this result, the discovery of squarks, sleptons, and heavy Higgs scalars may be extremely challenging even at the LHC. In addition, it is not surprising that these scalars have so far escaped detection, as present bounds are far from excluding most of the natural parameter space. Finally, it is tempting to speculate that what appears to be an accidental conspiracy between $`m_t`$ and the ratio of high to weak scales may find some fundamental explanation. If gauginos and Higgsinos are discovered, but all supersymmetric scalars escape detection at the LHC, the preservation of the naturalness motivation for supersymmetry, as currently understood, will require either an explanation of large cancellations between supersymmetry breaking soft masses at the weak scale, or the above scenario with a top mass fixed to be near 175 GeV. The latter possibility is, in our view, far more compelling and is supported by experimental data.

Acknowledgments — We are grateful to K. Agashe, M. Drees, and L. Hall for correspondence and conversations, and to the Aspen Center for Physics for hospitality. This work was supported in part by DOE under contracts DE–FG02–90ER40542 and DE–AC02–76CH03000, by the NSF under grant PHY–9513835, through the generosity of Frank and Peggy Taplin (JLF), and by a Marvin L. Goldberger Membership (TM).