text:
# DESY 02–166 ZU-TH 12/02 The Reconstruction of Supersymmetric Theories at High Energy Scales

## 1 Introduction

Extending the Standard Model to a supersymmetric theory  is an attractive step which has provided the qualitative understanding of a diverse set of phenomena in particle physics. Supersymmetry stabilizes the gap between the Grand Unification scale / Planck scale and the electroweak scale . It allows the unification of the three gauge couplings at a scale $`M_U210^{16}`$ GeV in a straight forward way . Radiative electroweak symmetry breaking relates to the high value of the top mass . Moreover, the cold dark matter component in the universe can be identified with the lightest supersymmetric particle . Above all, local supersymmetry, requiring the existence of massless spin 2 fields, provides a rationale for gravity .

Supersymmetry is not an exact symmetry in Nature. Unravelling the breaking mechanism is therefore one of the central issues with this new concept. A variety of mechanisms have been proposed, based on rather different physical ideas. Among these schemes are supergravity theories  which have provided the framework for many phenomenological analyses. The suppression of flavour–changing neutral reactions is achieved in an automatic form within gauge mediated supersymmetry breaking . Supersymmetry is broken in these scenarios in a hidden sector at high and intermediate scales, respectively, and the breaking is mediated by gravity or gauge interactions to the visible sector. The breaking, however, may not be communicated by direct action from the hidden to the visible sector. This is realized in anomaly mediated supersymmetry breaking models  in which supersymmetric particle masses are a consequence of the superconformal anomaly. In gaugino meditated supersymmetry breaking , supersymmetry is broken on a 3-brane separated from the 3-brane of the visible sector, and the breaking is communicated by gauge and Higgs superfields propagating through the 5-dimensional bulk. While in many models of supersymmetry breaking the gaugino masses are assumed to be universal at the unification scale, superstring motivated models, in which the breaking is moduli dominated, as opposed to dilaton dominated scenarios, give rise to non-universal boundary conditions at the high scale for the gauginos as well as the sfermion mass parameters . They can be exploited to determine the parameters of the string effective field theories.

In this report we elaborate on earlier investigations of Ref.  in which elements of gravity and gauge mediated supersymmetry breaking have been considered in realistic experimental environments of the proton collider LHC  and prospective TeV $`e^+e^{}`$ linear colliders . We extend these investigations in several directions in the present report.

In supergravity inspired models we adopt a scenario close to the Snowmass Point SPS#1 . In a second step, the previous analysis, based on the Minimal Supersymmetric Standard Model, is extended to a left–right supersymmetric $`SO(10)`$ model . The $`SO(10)`$ symmetry is assumed to be realized at a scale between the standard $`SU(5)`$ scale $`M_U210^{16}`$, derived from the unification of the gauge couplings, and the Planck scale $`M_P10^{19}`$ GeV. The right–handed neutrinos are assumed heavy, with masses at intermediate scales between $`O(10^{10})`$ GeV and $`O(10^{15})`$ GeV, so that the observed light neutrino masses are generated by the see-saw mechanism in a natural way . A rough estimate of the intermediate scale follows from the evolution of the mass parameters to the low experimental scale if universality holds at the Grand Unification scale.

In the gauge mediated symmetry breaking scenario, the fundamental scale is expected to be in the range from O(10 TeV) to O($`10^6`$ TeV). We present an update and an extension of the earlier analysis. In particular, the effective supersymmetry breaking scales, the messenger and supersymmetric mass scales, can be reconstructed at the point where the masses of the sparticles carrying the same quantum numbers become identical, the characteristic regularity of gauge mediated supersymmetry breaking.

The anomaly mediated as well as the gaugino mediated SUSY breaking are technically equivalent to the mSUGRA case and will therefore not be treated explicitly again.

Among the most exciting schemes rank superstring induced scenarios (see e.g  and references therein). In this report a string effective field theory, based on orbifold compactification of the heterotic string, will be analyzed in which, though dominated by the vacuum expectation values of the dilaton field, supersymmetry breaking is also affected by the moduli fields. Such a mechanism gives rise to gaugino mass parameters with small but noticeable departure from universality, and non-universal sfermion mass parameters. From these mass parameters the fundamental parameters of the string effective field theory, such as the vacuum expectation values of the dilaton and the moduli fields, the moduli/dilaton mixing angle as well as the modular weights can be derived. In this way high–precision experiments can provide access to elements which are directly induced by superstrings .

Extrapolations over many orders of magnitude from the electroweak scale to scales near the Planck scale require high–precision measurements at the laboratory scale . Such extrapolations can be performed in practice as demonstrated in the analysis of the electroweak and strong couplings at LEP and elsewhere . The unification of these couplings provides the most compelling argument, derived from experiment, in support of supersymmetry. An initial set of precision data on supersymmetric particles is expected from LHC experiments if favourable cascade decays can be exploited to measure mass differences very precisely . A globally comprehensive high–precision analysis can only be performed at lepton colliders . They are expected to be realized in a first phase up to an energy of about 1 TeV, and in a subsequent second phase up to about 5 TeV. $`e^+e^{}`$ linear collider designs for the first phase are being worked out for JLC, NLC and TESLA while the second phase may be realized in the CLIC technology. TESLA, in particular, can be operated at very high luminosity. A large number of threshold scans can therefore be performed which allows model independent high–precision measurements of the masses of supersymmetric particles. Chargino, neutralino and slepton masses are expected to be measured with accuracies at the per–mille level. Very heavy squarks and gluinos, on the other hand, may be analyzed in detail at CLIC after their discovery and first analysis at LHC. However, the accuracy is presumably reduced to the per–cent level as a consequence of the decreasing production cross sections, the non–zero widths of the heavy particles and the increasing energy smearing due to beam–strahlung.

Starting with observed numbers at the electroweak scale, the bottom–up approach exhausts all experimental information to the maximal extent possible in the empirical reconstruction of the underlying supersymmetric theory at the high scale. Finally, the parameters of the fundamental high-scale theory will become accessible in this way. This exploration of GUT and Planck scale physics by combining high precision with high energy in experiments at hadron and lepton colliders, is complemented by only a very small number of other methods, notabene proton decay, likely neutrino physics, textures of mass matrices, and cosmology. In all these individual approaches only scarce information on the underlying physical theoretical structures at the GUT/Planck scale can be extracted. Any of these methods should therefore be exploited in the maximal form in order to shed light on the boundaries of the physics area where gravity may affect properties and interactions of particles observed in the laboratory at the electroweak scale. In this way consequences of incorporating the fourth of the fundamental forces into the particle system could become accessible at laboratory experiments.

## 2 Gravity Mediated SUSY Breaking

### 2.1 Minimal Supergravity – mSUGRA

Supersymmetry cannot be broken spontaneously in our eigen–world without risking conflict with experimental results. The Ferrara–Girardello–Palumbo mass sum rule  requires supersymmetric particle masses below the corresponding Standard Model particle masses in this case — in obvious disagreement with observations. The elegant concept of spontaneous symmetry breaking, by non–perturbative gluino condensation for instance, can be realized, however, in a hidden sector which interacts with our eigen–world only by gravity. Gravitational interactions generate the soft supersymmetry breaking terms near the Grand Unification scale / Planck scale. Not compulsory but suggestive, the soft terms may be universal, i.e. the gaugino mass parameters and the scalar mass parameters111Universality may be broken GUT-scale threshold corrections, see e.g. Ref.. The bottom-up approach should enable us to explore this domain in a systematic way to the maximum extent possible.. Being flavour blind, the suppression of flavour–changing neutral processes can be realized in a natural way. Moreover, for a heavy top mass $`m_t174`$ GeV the breaking of the electroweak symmetry $`SU(2)_L\times U(1)_YU(1)_{EM}`$ can be generated radiatively. While at the universality scale all scalar masses squared are positive, the Higgs mass parameter $`M_{H_2}^2`$ turns negative at a scale of about $`10^4`$ TeV. This induces spontaneous electroweak symmetry breaking at the electroweak scale where the sum $`M_{H_2}^2+|\mu |^2`$ becomes negative, leaving however the strong and electromagnetic gauge symmetries $`SU(3)_C`$ and $`U(1)_{EM}`$ unbroken.

The minimal supergravity scenario mSUGRA is characterized by the universal parameters

| gaugino mass parameter | : | $`M_{1/2}`$ |
| --- | --- | --- |
| scalar mass parameter | : | $`M_0`$ |
| trilinear coupling | : | $`A_0`$ |

complemented by the phase of $`\mu `$, the modulus $`|\mu |`$ determined by radiative symmetry breaking, and the mixing angle $`\mathrm{tan}\beta `$ in the Higgs sector. The mass parameters $`M_{1/2}`$, $`M_0`$ and the trilinear coupling $`A_0`$ are defined to be universal at the Grand Unification scale $`M_U`$. The unified gauge coupling is denoted by $`\alpha _U`$ at $`M_U`$, for the sake of simplicity taken real. These parameters at the GUT scale are related to the low energy parameters by the supersymmetric renormalization group equations  which to leading order generate the evolution for

| gauge couplings | : | $`\alpha _i=Z_i\alpha _U`$ | (1) |
| --- | --- | --- | --- |
| gaugino mass parameters | : | $`M_i=Z_iM_{1/2}`$ | (2) |
| scalar mass parameters | : | $`M_{\stackrel{~}{j}}^2=M_0^2+c_jM_{1/2}^2+_{\beta =1}^2c_{j\beta }^{}\mathrm{\Delta }M_\beta ^2`$ | (3) |
| trilinear couplings | : | $`A_k=d_kA_0+d_k^{}M_{1/2}`$ | (4) |

The index $`i`$ runs over the gauge groups $`i=SU(3)`$, $`SU(2)`$, $`U(1)`$. To leading order, the gauge couplings, and the gaugino and scalar mass parameters of soft–supersymmetry breaking depend on the $`Z`$ transporters

$`Z_i=[1+b_i{\displaystyle \frac{\alpha _U}{4\pi }}\mathrm{log}\left({\displaystyle \frac{M_U}{M_Z}}\right)^2]^1`$ (5)
with $`b[SU_3,SU_2,U_1]=3,\mathrm{\hspace{0.17em}1},\mathrm{\hspace{0.17em}33}/5`$; the scalar mass parameters depend also on the Yukawa couplings $`h_t`$, $`h_b`$, $`h_\tau `$ of the top quark, bottom quark and $`\tau `$ lepton. The coefficients $`c_j`$ \[$`j=L_l,E_l,Q_l,U_l,D_l,H_{1,2}`$; $`l=1,2,3`$\] for the slepton and squark doublets/singlets of generation $`l`$, and for the two Higgs doublets are linear combinations of the evolution coefficients $`Z_i`$; the coefficients $`c_{j\beta }^{}`$ are of order unity. The shifts $`\mathrm{\Delta }M_\beta ^2`$ are nearly zero for the first two families of sfermions but they can be rather large for the third family and for the Higgs mass parameters, depending on the coefficients $`Z_i`$, the universal parameters $`M_0^2`$, $`M_{1/2}`$ and $`A_0`$, and on the Yukawa couplings $`h_t`$, $`h_b`$, $`h_\tau `$. The coefficients $`d_k`$ of the trilinear couplings $`A_k`$ \[$`k=t,b,\tau `$\] depend on the corresponding Yukawa couplings and they are approximately unity for the first two generations while being O($`10^1`$) and smaller if the Yukawa couplings are large; the coefficients $`d_k^{}`$, depending on gauge and Yukawa couplings, are of order unity. Beyond the approximate solutions shown explicitly, the evolution equations have been solved numerically in the present analysis to two–loop order  and threshold effects have been incorporated at the low scale .

These parameters enter the mass–matrices for the various particles. In the case of charginos $`\stackrel{~}{\chi }_m^+`$ \[$`m=1,2`$\] the $`2\times 2`$ mass matrix reads as

$`M_{\stackrel{~}{\chi }^+}=\left(\begin{array}{cc}M_2& \sqrt{2}m_W\mathrm{cos}\beta \\ \sqrt{2}m_W\mathrm{sin}\beta & \mu \end{array}\right)`$ (8)
while the mass matrix for neutralinos $`\stackrel{~}{\chi }_n^0`$ \[$`n=1,..,4`$\] is a $`4\times 4`$ matrix,

$`M_{\stackrel{~}{\chi }^0}=\left(\begin{array}{cccc}M_1& 0& m_Z\mathrm{cos}\beta \mathrm{sin}\theta _W& m_Z\mathrm{sin}\beta \mathrm{sin}\theta _W\\ 0& M_2& m_Z\mathrm{cos}\beta \mathrm{cos}\theta _W& m_Z\mathrm{sin}\beta \mathrm{cos}\theta _W\\ m_Z\mathrm{cos}\beta \mathrm{sin}\theta _W& m_Z\mathrm{cos}\beta \mathrm{cos}\theta _W& 0& \mu \\ m_Z\mathrm{sin}\beta \mathrm{sin}\theta _W& m_Z\mathrm{sin}\beta \mathrm{cos}\theta _W& \mu & 0\end{array}\right)`$ (13)
Exploiting all the information available from a linear collider, both mass matrices can be reconstructed even the in case of complex parameters . For large values, $`\mathrm{tan}\beta `$ needs supplementary analyses in the Higgs sector .

Assuming that the sfermion generations mix only weakly the mass matrices of the third generation sfermions can be written as:

$`M_{\stackrel{~}{f}}^2=\left(\begin{array}{cc}m_{\stackrel{~}{f}_L}^2& a_fm_f\\ a_fm_f& m_{\stackrel{~}{f}_R}^2\end{array}\right)`$ (16)
with

$`m_{\stackrel{~}{f}_L}^2=M_{\stackrel{~}{F}_L}^2+(T_f^3e_f\mathrm{sin}^2\theta _W)\mathrm{cos}2\beta m_Z^2+m_f^2,`$ (17)
$`m_{\stackrel{~}{f}_R}^2=M_{\stackrel{~}{F}_R}^2+e_f\mathrm{sin}^2\theta _W\mathrm{cos}2\beta m_Z^2+m_f^2,`$ (18)
$`a_tA_t\mu \mathrm{cot}\beta ,a_bA_b\mu \mathrm{tan}\beta ,a_\tau A_\tau \mu \mathrm{tan}\beta ,`$ (19)
where $`e_f`$ and $`T_f^3`$ are the electric charge and the third component of the weak isospin of the sfermion $`\stackrel{~}{f}`$; $`M_{\stackrel{~}{F}_L}=M_{\stackrel{~}{Q}}`$ for $`\stackrel{~}{f}_L=\stackrel{~}{t}_L,\stackrel{~}{b}_L`$, $`M_{\stackrel{~}{F}}=M_{\stackrel{~}{L}}`$ for $`\stackrel{~}{f}_L=\stackrel{~}{\tau }_L,\stackrel{~}{\nu }_\tau `$; $`M_{\stackrel{~}{F}_R}=M_{\stackrel{~}{U}},M_{\stackrel{~}{D}},M_{\stackrel{~}{E}}`$ for $`\stackrel{~}{f}_R=\stackrel{~}{t}_R,\stackrel{~}{b}_R,\stackrel{~}{\tau }_R`$, respectively; $`m_f`$ is the mass of the corresponding fermion. Also in this case it has been shown that the mass matrix can be reconstructed . The mass matrices for the first two generation sfermions have the same structure. However, due to the small fermion masses the mixing between the $`L/R`$ sfermions can be neglected in general. In the fits for the parameters we have used the complete one–loop mass matrices as given in . For the Higgs bosons also the two–loop contributions  are included.

The mSUGRA point we have analyzed in detail, was chosen close to the Snowmass Point SPS#1 , except for the scalar mass parameter $`M_0`$ which was taken slightly larger for merely illustrative purpose: $`M_{1/2}=250`$ GeV, $`M_0=200`$ GeV, $`A_0=100`$ GeV, $`\mathrm{tan}\beta =10`$ and $`sign(\mu )=+`$. The initial “experimental” values, have been generated by evolving the universal parameters down to the electroweak scale according to standard procedures .

The parameters chosen are compatible with the present results of low–energy experiments which they affect by virtual contributions, and they are also compatible with dark–matter estimates : BR($`bs\gamma )=2.710^4`$, $`\mathrm{\Delta }[g2]_\mu =1710^{10}`$, $`\mathrm{\Delta }\rho =3810^5`$ and $`\mathrm{\Omega }h^2=0.4`$. We have used the formulas given in  for the computation of $`bs\gamma `$, those given in  for $`\mathrm{\Delta }\rho `$ and those given in  for $`\mathrm{\Delta }[g2]_\mu `$. $`\mathrm{\Omega }h^2`$ has been calculated using the program of Ref..

The five basic parameters define the experimental observables, including supersymmetric particle masses and production cross sections. They are endowed with errors as expected for threshold scans as well as measurements in the continuum at $`e^+e^{}`$ linear colliders (LC). Major parts of the LC analysis can be performed for energies below 1 TeV, some of the squarks require energies above 1 TeV. Estimates are based on integrated LC luminosities of 1 ab-1. The errors given in Ref. are scaled in proportion to the masses of the spectrum. Typical examples are shown in Table 1. The LC errors on the squark masses, see e.g. Ref., are set to an average value of 10 GeV \[similar errors may also be obtained if the precisely measured mass differences at the LHC are combined with high-precision measurements of the low-lying states at the LC\]; varying this error within a factor two does not change the conclusions significantly since the measurement of the cross sections provides the maximal sensitivity in this sector. For the cross-sections we use purely statistical errors, while assuming a (conservative) reconstruction efficiency of 20%. In addition the mass errors on the lightest gauginos were inflated with respect to earlier analyses to be conservative in advance of detailed experimental analyses of models with higher values of $`\mathrm{tan}\beta `$. \[Parameter combinations from the fits to the spectrum and the cross sections which lead to charge and/or color breaking minima , are not accepted.\] These observables are interpreted as the experimental input values for the evolution of the mass parameters in the bottom-up approach to the Grand Unification scale.

#### 2.1.1 Gauge Couplings

The presumably strongest support, though indirect, for supersymmetry is related to the tremendous success of this theory in predicting the unification of the gauge couplings . The precision, being at the per–cent level, is surprisingly high after extrapolations over fourteen orders of magnitude in the energy from the electroweak scale to the unification scale $`M_U`$. Conversely, the electroweak mixing angle has been predicted in this approach at the per–mille level. The evolution of the gauge couplings from low energy to the GUT scale $`M_U`$ is carried out in the two–loop accuracy. The gauge couplings $`g_1`$, $`g_2`$, $`g_3`$ and the Yukawa couplings are calculated in the $`\overline{DR}`$ scheme by adopting the shifts given in . These parameters are evolved to $`M_U`$ using 2–loop RGEs . At 2-loop order the gauge couplings do not meet exactly , the differences owing to threshold effects at the unification scale $`M_U`$ which leave us with an ambiguity in the definition of $`M_U`$. In this report we define $`M_U`$ as the scale, ad libitum, where $`g_1=g_2`$ in the RGE evolution. The non–zero difference $`g_1g_3`$ at this scale is then attributed to threshold effects of particles with masses of order $`M_U`$. The quantitative evolution implies important constraints on the particle content at $`M_U`$ -.

Based on the set of low–energy gauge and Yukawa parameters $`\{\alpha (m_Z)`$, $`\mathrm{sin}^2\theta _W`$, $`\alpha _s(m_Z)`$, $`Y_t(m_Z)`$, $`Y_b(m_Z)`$, $`Y_\tau (m_Z)\}`$ the evolution of the inverse couplings $`\alpha _i^1`$ $`[i=U(1)`$, $`SU(2)`$, $`SU(3)]`$ is depicted in Fig. 1a. The evolution is performed for the mSUGRA reference point defined above. Unlike earlier analyses, the low–energy thresholds of supersymmetric particles can be calculated in this framework exactly without reference to effective SUSY scales. The broken error ellipse in Fig. 1b, derived for \[$`M_U,\alpha _U`$\] by requiring $`g_1=g_2`$, corresponds to the present experimental accuracy of the gauge couplings : $`\mathrm{\Delta }\{\alpha ^1(m_Z)`$, $`\mathrm{sin}^2\theta _W`$, $`\alpha _s(m_Z)\}`$ $`=\{0.03,1.710^4,310^3\}`$. The full ellipse demonstrates the improvement for the absolute errors $`\{10^3,10^5,10^3\}`$ after operating GigaZ . The expected accuracies in $`M_U`$ and $`\alpha _U`$ are summarized in the values given in Table 2. The difference between the unification point in the ellipse and the value of $`g_3`$ is accounted for by contributions from high scale physics, colour-triplet Higgs fields, for example. Thus, for a typical set of SUSY parameters, the evolution of the gauge couplings from low to high scales leads to a precision of 1.5 per–cent of the Grand Unification picture.

#### 2.1.2 Gaugino and Scalar Mass Parameters

The results for the evolution of the mass parameters to the GUT scale $`M_U`$ are shown in Fig. 2. Fig. 2(a) presents the evolution of the gaugino parameters $`M_i^1`$ which clearly is under excellent control, as is the extrapolation of the slepton mass parameter squared of the first (and second) and the third generation in Fig. 2(c) and (d), respectively. The accuracy deteriorates for the squark mass parameters and for the Higgs mass parameter $`M_{H_2}^2`$. The origin of the differences between the errors for slepton, squark and Higgs mass parameters can be traced back to the numerical size of the coefficients in Eqs. (2.1). Typical examples using the formulas presented in Appendix B evaluated at $`Q=500`$ GeV read as follows:

$`M_{\stackrel{~}{L}_1}^2`$ $``$ $`M_0^2+0.47M_{1/2}^2`$ (20)
$`M_{\stackrel{~}{Q}_1}^2`$ $``$ $`M_0^2+5.0M_{1/2}^2`$ (21)
$`M_{\stackrel{~}{H}_2}^2`$ $``$ $`0.03M_0^21.34M_{1/2}^2+1.5A_0M_{1/2}+0.6A_0^2`$ (22)
$`|\mu |^2`$ $``$ $`0.03M_0^2+1.17M_{1/2}^22.0A_0M_{1/2}0.9A_0^2`$ (23)
While the coefficients for sleptons are of order unity, the coefficient $`c_j`$ for squarks grows very large, $`c_j5.0`$, so that small errors in $`M_{1/2}^2`$ are magnified by nearly an order of magnitude in the solution for $`M_0`$. By close inspection of Eq.(2.1) for the Higgs mass parameter it turns out that the formally leading $`M_0^2`$ part is nearly cancelled by the $`M_0^2`$ part of $`c_{j,\beta }^{}\mathrm{\Delta }M_\beta ^2`$. Inverting Eq.(2.1) for $`M_0^2`$ therefore gives rise to large errors in the Higgs case. A representative set of mass values and the associated errors, after evolution from the electroweak scale to $`M_U`$, is presented in Table 3. The corresponding error ellipses for the unification of the gaugino masses are shown in Fig. 2(b).

Extracting the trilinear parameters $`A_k`$ is difficult and more refined analyses based on sfermion cross sections and Higgs and/or sfermion decays are necessary to determine these parameters accurately. $`A_t`$ can be obtained from the mixing angle of the stop sector by measuring the stop production cross section in $`e^+e^{}`$ annihilation with different electron and/or positron polarizations . In the cases $`A_b`$ and $`A_\tau `$ the situation is more difficult, because these parameters influence the mixing angle in the sbottom and stau sector only weakly as evident from Eq.(19). In these cases the $`\stackrel{~}{b}`$ and $`\stackrel{~}{\tau }`$ couplings to the Higgs bosons must be measured, because these couplings include terms directly proportional to $`A_k\mathrm{tan}\beta `$. For instance, by analysing the decays $`\stackrel{~}{\tau }_2A^0\stackrel{~}{\tau }_1`$, $`h^0\stackrel{~}{\tau }_1`$ and $`H^0\stackrel{~}{\tau }_1`$, $`A_\tau `$ can be extracted within 10% . If these modes are kinematically forbidden, the couplings can either be measured in the decays of the heavier Higgs bosons, as $`H^0,A^0\stackrel{~}{\tau }_1\stackrel{~}{\tau }_1`$, or by means of the cross sections for processes such as $`e^+e^{}\stackrel{~}{\tau }_1\stackrel{~}{\tau }_1h^0`$. Similar procedures are expected to apply for $`A_b`$. In certain areas of SUSY parameter space, the trilinear couplings can also be extracted from measurements of the degree of fermion polarization  in sfermion decays $`\stackrel{~}{t},\stackrel{~}{b}`$ and $`\stackrel{~}{\tau }`$.

The unified value $`A_0`$ of the $`A_t`$ coupling, the best measured coupling among the $`A_k`$ parameters, is shielded by the pseudo–fixed point behaviour of $`A_t`$  since $`d_t0.2`$ is small compared to $`d_t^{}2`$. The impact of the other trilinear couplings on physical observables is weak so that large experimental errors are expected. As a result, the universal character of the fundamental parameter $`A_0`$ cannot be determined as precisely as the other parameters at the GUT scale. Although the trilinear couplings $`A_b`$ and $`A_\tau `$ have only little impact on physical observables, they do strongly influence the running of the third generation sfermion mass parameters as well as the Higgs mass parameters. The error propagation is stabilized if $`A_\tau `$ and $`A_b`$ can be measured in the way outlined above. \[Otherwise the errors would increase by an order of magnitude.\] The detailed analysis in this report has been based on the auxiliary assumption that $`A_b`$ and $`A_\tau `$ are within 1$`\sigma `$ of $`A_t=A_0`$ at $`M_U`$; this assumption is conservative if the envisaged experimental analyses of $`A_\tau `$ and $`A_b`$ can be performed at the electroweak scale in practice.

Even though the auxiliary assumption seems conservative, given the size of the error on $`A_0`$ determined from $`A_t`$, dedicated phenomenological and experimental analyses of the third family must be developed, as indicated above, to improve the measurement of the associated parameters, in particular in view of the evolution of the Higgs mass parameter which induces electroweak symmetry breaking.

Inspecting Figs. 2(c) and (d) leads to the conclusion that a blind top-down approach eventually may generate an incomplete picture. Global fits based on mSUGRA without allowing for deviations from universality, are dominated by $`M_{1,2}`$ and the slepton mass parameters due to the pseudo-fixed point behaviour of the squark mass parameters. Therefore, the structure of the theory in the squark sector is not scrutinized stringently at the unification scale in the top-down approach let alone the Higgs sector. By contrast, the bottom-up approach demonstrates very clearly the extent to which the theory can be tested at the high scale quantitatively. The quality of the test is apparent from Table 3, in which the evolved gaugino values should reproduce the universal mass $`M_{1/2}=250`$ GeV and all the scalars the mass $`M_0=200`$ GeV. They are compared with the global mSUGRA fit of the universal parameters in Table 4.

### 2.2 Left–Right Supergravity

It is generally accepted that neutrinos are massive particles, though at a very low scale. Supersymmetric scenarios like MSSM and mSUGRA must therefore be extended to incorporate the right–handed neutrino degrees of freedom. Since the complexity grows strongly with the rising number of parameters, it is useful, in a first attempt, to analyze the system in characteristic scenarios based on compelling physical assumptions. In particular, we will assume that the small neutrino masses are generated by the seesaw mechanism . Moreover, we will assume hierarchies for the heavy neutrino masses as well as the neutrino Yukawa couplings similar to the up–type particles in the quark sector; such a scheme, suggested by $`SO(10)`$ GUT, is compatible with the data collected in low–energy neutrino experiments .

This scenario can be embedded in a Grand Unified $`SO(10)`$ theory with the following breaking pattern of the symmetries. The $`SO(10)`$ symmetry is realized between the Planck scale $`M_P`$ and a scale $`M_{SO(10)}`$ at which the symmetry breaks to $`SU(5)`$. The scale $`M_{SO(10)}`$ is assumed above the scale $`M_U`$ where $`SU(5)`$ breaks to the symmetry group $`SU(3)_C\times SU(2)_L\times U(1)_Y`$ of the Standard Model. At the scale $`M_U`$ the gauge couplings split and the effective theory is the MSSM plus right–handed neutrinos with masses of order $`10^9`$ to $`10^{15}`$ GeV. Below this mass scale the right–handed neutrinos freeze out and the MSSM is effectively realized in its standard form. The relevant SUSY parameters are summarized in Table 5. It is less obvious that $`M_U`$ associated with the $`SU(5)`$ symmetry is the scale where the gaugino and scalar mass parameters are universal. The supporting argument for this point is derived empirically from the unification of the gauge couplings. Nevertheless, the subsequent analysis will be based on this hypothesis which, of course, is a clear target for confirmation or rejection in the bottom–up approach we investigate222Potential sources of deviations from this picture can easily be illustrated by assuming $`M_{SO(10)}`$ as the scale proper of universality: The Yukawa interactions contribute differently to the running of the $`M_{\stackrel{~}{10}}^2`$, $`M_{\stackrel{~}{5}}^2`$, $`M_{\stackrel{~}{\nu }_R}^2`$; the same holds true for the $`A`$ parameters . Moreover, different $`D`$–term contributions to the scalar masses are in general generated by the breaking mechanism from $`SO(10)`$ to $`SU(5)`$ ..

In this left–right supergravity point, called LR–SUGRA for short333 Other left–right scenarios will be presented in a forthcoming publication., we have probed the same SUSY parameters as above, complemented by the same universal parameters in the right–handed sneutrino sector. The sneutrinos $`\stackrel{~}{\nu }_L`$ and $`\stackrel{~}{\nu }_R`$ mix by the (large) Yukawa interactions in the $`\widehat{\nu }_R`$ sector of the superpotential to form the mass eigenstates $`\stackrel{~}{\nu }_1`$ and $`\stackrel{~}{\nu }_2`$. Also in this sector an effective seesaw mechanism is induced by the large $`\nu _R`$ mass, as can be most easily seen by considering the one generation case:

$`m^2=\left(\begin{array}{cc}M_{\stackrel{~}{L}}^2+\frac{1}{2}m_Z^2\mathrm{cos}2\beta & \frac{1}{\sqrt{2}}Y_\nu (A_\nu v_2\mu v_1)\\ \frac{1}{\sqrt{2}}Y_\nu (A_\nu v_2\mu v_1)& M_{\stackrel{~}{\nu }_R}^2+M_{\nu _R}^2\end{array}\right)`$ (26)
In this mass matrix $`M_{\nu _R}`$ is the \[GUT-scale\] mass of the right–handed neutrino, $`M_{\stackrel{~}{\nu }_R}`$ the scalar \[TeV-scale\] mass parameter of the right sneutrino, $`Y_\nu `$ and $`A_\nu `$ the neutrino Yukawa coupling and the neutrino trilinear coupling, respectively. $`v_1`$ and $`v_2`$ are the vacuum expectation values of the Higgs field with isospin $`1/2`$ and isospin $`1/2`$, respectively. The approximate eigenvalues of the sneutrino mass matrix read

$`m_{\stackrel{~}{\nu }_1}^2`$ $``$ $`M_{\stackrel{~}{L}}^2+{\displaystyle \frac{1}{2}}m_Z^2\mathrm{cos}2\beta Y_\nu ^2{\displaystyle \frac{(A_\nu v_2\mu v_1)^2}{2M_{\nu _R}^2}}`$ (27)
$`m_{\stackrel{~}{\nu }_2}^2`$ $``$ $`M_{\nu _R}^2+M_{\stackrel{~}{\nu }_R}^2+Y_\nu ^2{\displaystyle \frac{(A_\nu v_2\mu v_1)^2}{2M_{\nu _R}^2}}`$ (28)
Therefore, to a very good approximation, $`\stackrel{~}{\nu }_1`$ coincides with $`\stackrel{~}{\nu }_L`$ and $`\stackrel{~}{\nu }_2`$ with $`\stackrel{~}{\nu }_R`$. The heavy right–handed neutrino masses are calculated by identifying the Yukawa couplings with the up-type quark couplings in the quark sector at the GUT scale (largely equivalent to the $`SO(10)`$ scale in this regard) and by identifying the light neutrino masses with the neutrino mass differences in the large mixing angle for the solar neutrino problem solution: $`m_{\nu _{L1}}=10^5`$ eV, $`m_{\nu _{L2}}=310^3`$ eV, $`m_{\nu _{L3}}=610^2`$ eV; $`M_{\nu _{R_1}}=310^9`$ GeV, $`M_{\nu _{R_2}}=1.410^{11}`$ GeV, $`M_{\nu _{R_3}}=1.710^{14}`$ GeV.

The impact on the evolution of the mass parameters is rather simple. In the analysis of the first two generations the Yukawa interactions involving the heavy neutrinos and the R-sneutrinos are so small that their effect is not noticeable in practice. The evolution of the gaugino and scalar mass parameters is not affected by the left–right extension of the system in the present form as is evident from Figs. 3(a) and (c). This is only different for the third generation and for $`M_{H_2}^2`$ owing to the enhanced Yukawa coupling in this case as shown in Figs. 3(b) and (d). The sensitivity to the intermediate $`\nu _R`$ scales is rather weak because neutrino Yukawa couplings affect the evolution of the sfermion mass parameters only mildly. Since the $`\nu _R`$ of the third generation is unfrozen only beyond $`Q=M_{\nu _R}`$ the impact of the LR extension becomes visible in the evolution only at very high scales. In Fig. 3(b) we display the evolution of $`M_{\stackrel{~}{E}_3}^2`$, $`M_{\stackrel{~}{L}_3}^2`$ and $`M_{H_2}^2`$ for illustrative purposes. The full lines include the effects of the right–handed neutrino, which are to be compared with the dashed lines where the $`\nu _R`$ effects are removed. The scalar mass parameter $`M_{\stackrel{~}{E}_3}^2`$ appears unaffected by the right–handed sector, while $`M_{\stackrel{~}{L}_3}^2`$ and $`M_{H_2}^2`$ clearly are. Only the picture including $`\nu _R`$, $`\stackrel{~}{\nu _R}`$ is compatible with the unification assumption. The kinks in the evolution of $`M_{\stackrel{~}{L}_3}^2`$ and $`M_{\stackrel{~}{H_2}}^2`$ can be traced back to the fact that around $`10^{14}`$ GeV the third generation (s)neutrinos become quantum mechanically effective, given a large enough neutrino Yukawa coupling to influence the evolution of these mass parameters. A much better understanding of the third generation family must be achieved to draw quantitative conclusions beyond the rough estimates of the $`\nu _R`$ scales sketched in the present analysis.

## 3 Gauge Mediated Supersymmetry Breaking

Motivated by the observed suppression of flavour changing neutral transitions, supersymmetry breaking mediated by gauge interactions from a secluded sector to the visible eigen–world, offers an automatic solution to this problem . The scalar and the F components of a Standard–Model singlet superfield $`S`$ acquire vacuum expectation values $`S`$ and $`F_S`$ through interactions with fields in the secluded sector, thus breaking supersymmetry.444 A solution of the doublet–triplet splitting problem can be found in GMSB by introducing two different $`S`$ fields. The masses of supersymmetric particles are less constrained in this approach than in the one-scale model, and they depend on the values of the two $`\mathrm{\Lambda }_i`$ parameters. In particular, the approximate equality of the gaugino masses at the GUT scale is lifted, see  for details. Vector-like messenger fields $`M`$, carrying non–zero $`SU(3)\times SU(2)\times U(1)`$ charges and coupling to $`S`$, transport the supersymmetry breaking to the eigen–world. The system is characterized by the mass $`M_M`$ $`S`$ of the messenger fields and the mass scale $`\mathrm{\Lambda }=F_S/S`$ setting the seize of the gaugino and scalar masses. $`M_M`$ is expected to be in the range of 10 to $`10^6`$ TeV and $`\mathrm{\Lambda }`$ has to be smaller than $`M_M`$. The gaugino masses

$$M_i(M_M)=(N_5+3N_{10})g\left(\mathrm{\Lambda }/M_M\right)\alpha _i(M_M)\mathrm{\Lambda }$$
(29)
are generated by loops of scalar and fermionic messenger component fields; $`N_i`$ is the multiplicity of messengers in the $`5+\overline{5}`$ and $`10+\overline{10}`$ vector-like multiplets, and

$`g(x)`$ $`=`$ $`{\displaystyle \frac{1+x}{x^2}}\mathrm{log}(1+x)+(xx)`$ (30)
is the messenger–scale threshold function  which approaches unity for $`\mathrm{\Lambda }M_M`$. Masses of the scalar fields in the visible sector are generated by 2-loop effects of gauge/gaugino and messenger fields:

$$M_{\stackrel{~}{j}}^2(M_M)=2(N_5+3N_{10})f\left(\mathrm{\Lambda }/M_M\right)\underset{i=1}{\overset{3}{}}k_iC_i^j\alpha _i^2(M_M)\mathrm{\Lambda }^2$$
(31)
with $`k_i=1,1,3/5`$ for $`SU(3)`$, $`SU(2)`$, and $`U(1)`$, respectively; the coefficients $`C_i^j`$ are the quadratic Casimir invariants, being 4/3, 3/4, and $`Y^2/4`$ for the fundamental representations $`\stackrel{~}{j}`$ in the groups $`i=SU(3),SU(2)`$ and $`U(1)`$, with $`Y=2(QI_3)`$ denoting the usual hypercharge; also the threshold function

$`f(x)`$ $`=`$ $`{\displaystyle \frac{1+x}{x^2}}\left[\mathrm{log}(1+x)2\mathrm{L}\mathrm{i}_2\left({\displaystyle \frac{x}{1+x}}\right)+{\displaystyle \frac{1}{2}}\mathrm{Li}_2\left({\displaystyle \frac{2x}{1+x}}\right)\right]+(xx)`$ (32)
approaches unity for $`\mathrm{\Lambda }M_M`$. As evident from Eq.(31) scalar particles with identical Standard–Model charges squared have equal masses at the messenger scale $`M_M`$. In the minimal version of GMSB, the $`A`$ parameters are generated at 3-loop level and they are practically zero at $`M_M`$.

We have investigated this scheme for the point $`\mathrm{\Lambda }=100`$ TeV, $`M_M=200`$ TeV, $`N_5=1`$, $`N_{10}=0`$, $`\mathrm{tan}\beta =15`$ and $`\mu >0`$ corresponding to the Snowmass Point SPS#8. We find for the low energy data: BR($`bs\gamma )=3.710^4`$, $`\mathrm{\Delta }[g2]_\mu =1510^{10}`$, $`\mathrm{\Delta }\rho =6410^5`$. The evolution555The same formulae as in Appendix B.1 apply for the GMSB boundary conditions at the electroweak scale by replacing $`M_U`$ by $`M_M`$, the GMSB scale. of the gaugino and sfermion mass parameters of the first and third generation as well as the Higgs mass parameters, including 2-loop $`\beta `$–functions, are presented in Fig. 4. Owing to the influence of the $`A`$–parameters in the 2-loop RGEs for the gaugino mass parameters, the gaugino mass parameters do not meet at the same point as the gauge couplings in this scheme. It is obvious from the figure that the GMSB scenario cannot be confused with the universal supergravity scenario666A comparison of the mass characteristics at the low scale between mSUGRA and GMSB in a top-down approach is presented in Ref... \[Specific experimental signatures generated in the decays of the next to lightest supersymmetric particle, the neutralino $`\stackrel{~}{\chi }_1^0`$ or the stau $`\stackrel{~}{\tau }_1`$, to gravitinos which are very light in GMSB, provide a complementary experimental discriminant, see \].

The bands of the slepton $`L`$–doublet mass parameter $`M_{\stackrel{~}{L}}^2`$ and the Higgs parameter $`M_{H_2}^2`$, which carry the same moduli of standard–model charges, cross at the scale $`M_M`$. The crossing, which is indicated by an arrow in Fig. 4(c), is a necessary condition (in the minimal form) for the GMSB scenario to be realized.

The two scales $`\mathrm{\Lambda }`$ and $`M_M`$ can be extracted from the spectrum of the gaugino and scalar particles. Combining the two species allows to determine the multiplicity coefficient $`(N_5+3N_{10})`$ in addition. The dependence of the spectra on $`\mathrm{\Lambda }`$ is, trivially, very strong. The messenger scale $`M_M`$ can be determined only from the running of the masses between the messenger scale and the electroweak scale ; despite being governed by weakly varying logarithms, the accuracy in determining $`M_M`$ is surprisingly good. For the point analyzed in the example above, the following accuracy for the mass parameters and the messenger multiplicity has been found:

$`\mathrm{\Lambda }`$ $`=`$ $`(1.01\pm 0.03)10^2\mathrm{TeV}`$ (33)
$`M_M`$ $`=`$ $`(1.92\pm 0.24)10^2\mathrm{TeV}`$ (34)
$`N_5+3N_{10}`$ $`=`$ $`0.978\pm 0.056`$ (35)
The correlation between $`\mathrm{\Lambda }`$ and $`M_M`$ is shown in Fig. 4(b).

While the gaugino masses in GMSB evolve nearly in the same way as in mSUGRA, the running of the scalar masses is quite different in both theories. Moreover, at the messenger scale the ratios of scalar masses squared in the simplest version of GMSB are determined solely by group factors and gauge couplings, being independent of the specific GMSB characteristics, i. e. messenger multiplicities and $`\mathrm{\Lambda }`$ mass scale:

$$M_{\stackrel{~}{j}}^2(M_M)/M_{\stackrel{~}{j}^{}}^2(M_M)=\underset{i=1}{\overset{3}{}}k_iC_i^j\alpha _i^2(M_M)/\underset{i=1}{\overset{3}{}}k_iC_i^j^{}\alpha _i^2(M_M).$$
(36)
The predictions for these ratios are listed in Table 6. The ratios in GMSB are distinctly different from the ratios in mSUGRA, taken at the scale where the upper boundary of the $`2\sigma `$ band for $`H_2^2/L_1^2`$ approaches unity from below. \[Ideally all ratios approach unity only at the grand unification scale $`M_U`$ in mSUGRA.\] The distinct differences between GMSB and mSUGRA are clearly visible in Figs. 5(a) vs. (b).

## 4 String Induced Supersymmetry Breaking

In the previously analyzed SUGRA and GMSB models the structure of the supersymmetry breaking mechanisms sui generis and the fields involved in the hidden sectors are shielded from the eigen–world. Four–dimensional strings naturally give rise to a minimal set of fields for inducing supersymmetry breaking; they play the rôle of the fields in the hidden sectors: the dilaton $`S`$ and the moduli $`T_m`$ chiral superfields which are generically present in large classes of 4–dimensional heterotic string theories777For other scenarios see Ref... The vacuum expectation values of $`S`$ and $`T_m`$, generated by genuinely non–perturbative effects, determine the soft supersymmetry breaking parameters. In this approach, grand unification at the standard scale can be reconciled with the higher string scale by moduli dependent string loop corrections.

In the following we assume that all moduli fields get the same vacuum expectation values and that they couple in the same way to matter fields. Therefore, we omit the index $`m`$ and take only one moduli field $`T`$. The properties of the supersymmetric theories are quite different for dilaton and moduli dominated scenarios. This can be quantified by introducing a mixing angle $`\theta `$, characterizing the $`\stackrel{~}{S}`$ and $`\stackrel{~}{T}`$ wave functions of the Goldstino, which is associated with the breaking of supersymmetry and which is absorbed to generate the mass of the gravitino:

$`\stackrel{~}{G}`$ $`=`$ $`\mathrm{sin}\theta \stackrel{~}{S}+\mathrm{cos}\theta \stackrel{~}{T}`$ (37)
The mass scale is set by the second parameter of the theory, the gravitino mass $`m_{3/2}`$.

A dilaton dominated scenario, i.e. $`\mathrm{sin}\theta 1`$, leads to universal boundary conditions of the soft supersymmetry breaking parameters. Universality is broken888For other mechanisms of breaking universality see e.g. Ref.. only slightly by small loop effects. On the other hand, in moduli field dominated scenarios, $`\mathrm{cos}\theta 1`$, the gaugino mass parameters are universal to lowest order \[broken only in higher orders\], but universality is not realized for the scalar mass parameters in general. The breaking is quantified by modular weights $`n_j`$ characterizing the couplings between the matter and the moduli fields in orbifold compactifications. Within one generation significant differences between left and right field components and between sleptons and squarks can occur; since these patterns are modified only by small loop effects between different generations, flavour–changing neutral effects remain suppressed.

In leading order, and next–to–leading order parameterized by the quantities $`\mathrm{\Delta }M`$, the masses  are given by the following expressions for the gaugino sector,

$`M_i`$ $`=`$ $`g_i^2m_{3/2}s\sqrt{3}\mathrm{sin}\theta +\mathrm{\Delta }M_i`$ (38)
$`\mathrm{\Delta }M_i`$ $`=`$ $`g_i^2m_{3/2}\{b_i+s\sqrt{3}\mathrm{sin}\theta g_s^2(C_i{\displaystyle \underset{j}{}}C_i^j)`$ (39)
$`+2t\mathrm{cos}\theta G_2(t)[\delta _{\mathrm{GS}}+b_i2{\displaystyle \underset{j}{}}C_i^j(1+n_j)]\}/16\pi ^2`$
and for the scalar sector,

$`M_{\stackrel{~}{j}}^2`$ $`=`$ $`m_{3/2}^2\left(1+n_j\mathrm{cos}^2\theta \right)+\mathrm{\Delta }M_{\stackrel{~}{j}}^2`$ (40)
$`\mathrm{\Delta }M_{\stackrel{~}{j}}^2`$ $`=`$ $`m_{3/2}^2\{\gamma _j+2t\mathrm{cos}\theta G_2(t){\displaystyle \underset{km}{}}\gamma _j^{km}(n_j+n_k+n_m+3)`$ (41)
$`+2\sqrt{3}s\mathrm{sin}\theta [{\displaystyle \underset{i}{}}\gamma _j^ig_i^2{\displaystyle \frac{1}{2s}}{\displaystyle \underset{km}{}}\gamma _j^{km}]\}`$
while the $`A`$ parameters read

$`A_{jkm}`$ $`=`$ $`m_{3/2}\left[2t\mathrm{cos}\theta (n_j+n_k+n_m+3)G_2(t){\displaystyle \frac{\mathrm{sin}\theta }{\sqrt{3}}}\right]+\mathrm{\Delta }A_{jkm}`$ (42)
$`\mathrm{\Delta }A_{jkm}`$ $`=`$ $`m_{3/2}(\gamma _j+\gamma _k+\gamma _m)`$ (43)
The mass $`m_{3/2}`$ is the gravitino mass introduced earlier. The gravitino mass can be expressed in terms of the Kähler potential $`K`$ and the superpotential $`W`$, which include the (non-perturbative) solutions of all the fields at the string scale: $`m_{3/2}=\mathrm{exp}(K/2)\overline{W}`$. $`s=S`$ is the vacuum expectation values of the dilaton field. $`t=T/m_{3/2}`$ is the vacuum expectation value of the moduli field(s), and $`G_2(t)=2\zeta (t)+1/2t`$ is the non-holomorphic Eisenstein function with $`\zeta `$ denoting the Riemann zeta function. $`\delta _{GS}`$ is the parameter of the Green-Schwarz counterterm. $`\gamma _j`$ are the anomalous dimensions of the matter fields, the $`\gamma _j^i`$ and $`\gamma _j^{km}`$ are their gauge and Yukawa parts, respectively. $`C_i`$, $`C_i^j`$ are the quadratic Casimir operators for the gauge group $`G_i`$, respectively, in the adjoint representation and in the matter representation.

In the case of the gaugino mass parameters the next–to–leading order effects induce a splitting proportional to the beta functions $`b_i`$ which is large enough to be “measured” at future collider experiments as demonstrated in Fig. 6(a) and (b). In the case of the scalar mass parameters the next–to–leading order contributions generate small departures from non-universality between the generations even if the corresponding modular weight is generation independent. These departures are proportional to the Yukawa couplings squared so that the third generation, in particular the stop sector, is mainly affected.

Scenarios have been found in which the phenomenological unification of the three gauge couplings can be reconciled with a string mass scale which is an order of magnitude larger than the unification scale :

$`\alpha _i^1(M_U)`$ $`=`$ $`\alpha ^1(M_{String})+\mathrm{\Delta }\alpha _i^1`$ (44)
The corrections $`\mathrm{\Delta }\alpha _i^1`$ depend on the value of the moduli fields and the modular weights:

$`\mathrm{\Delta }\alpha _i^1={\displaystyle \frac{1}{4\pi }}(b_i^{}b_{GS})\mathrm{log}|\eta (t)|^4`$ (45)
where $`\eta (t)`$ is the Dedekind $`\eta `$ function and

$`b_3^{}`$ $`=`$ $`9+{\displaystyle \underset{i=1}{\overset{3}{}}}\left(2n_{Q_i}+n_{U_i}+n_{D_i}\right)`$ (46)
$`b_2^{}`$ $`=`$ $`15+{\displaystyle \underset{i=1}{\overset{3}{}}}\left(3n_{Q_i}+n_{L_i}\right)+n_{H_1}+n_{H_2}`$ (47)
$`b_1^{}`$ $`=`$ $`{\displaystyle \frac{99}{5}}+{\displaystyle \frac{1}{5}}{\displaystyle \underset{i=1}{\overset{3}{}}}\left(n_{Q_i}+8n_{U_i}+2n_{D_i}+3n_{L_i}+6n_{E_i}\right)+{\displaystyle \frac{3}{5}}\left(n_{H_1}+n_{H_2}\right)`$ (48)
as compared to the one–loop $`SU(3)\times SU(2)\times U(1)`$ $`\beta `$–functions $`(b_3,b_2,b_1)`$ $`=`$ $`33/5,1,3`$.

We have analyzed a mixed dilaton/moduli superstring scenario but with dominating dilaton field component, $`\mathrm{sin}^2\theta =0.9`$, and with different couplings of the moduli field to the (L,R) sleptons, the (L,R) squarks and to the Higgs fields, corresponding to O–I representation $`n_{L_i}=3`$, $`n_{E_i}=1`$, $`n_{H_1}=n_{H_2}=1`$, $`n_{Q_i}=0`$, $`n_{D_i}=1`$ and $`n_{U_i}=2`$, that is one out of several assignments that is adopted quite frequently in the literature. The gravitino mass is chosen to be 180 GeV in this analysis. We find for the low energy data: BR($`bs\gamma )=3.110^4`$, $`\mathrm{\Delta }[g2]_\mu =1410^{10}`$, $`\mathrm{\Delta }\rho =1310^5`$; and $`\mathrm{\Omega }h^2=0.25`$.

Given this set of superstring induced parameters, the evolution of the gaugino and scalar mass parameters is displayed in Fig. 6. The pattern of the trajectories is remarkably different from other scenarios. The breaking of universality in the gaugino sector, induced by string threshold corrections, is pronounced, cf. Table 7.

In fact, the differences can be exploited to determine superstring parameters as argued above. The number of observables in the set of gauge couplings $`g_a`$, gaugino masses $`M_a`$ and scalar masses $`M_{\stackrel{~}{j}}`$ exceeds the number of parameters in the superstring effective field theory: the gravitino mass $`m_{3/2}`$, the dilaton/moduli mixing angle $`\mathrm{sin}\theta `$, the ground–state value of the moduli field $`T`$ and the ground–state value of the dilaton field $`S`$. The latter is at tree–level directly related to the string coupling: $`1/g_s^2=S`$.

Based on the “experimental” input observables, the fundamental parameters of the string effective field theory can be reconstructed; the reconstructed values are compared with the ideal values in Table 8. The errors for the basic parameters $`\mathrm{sin}\theta `$, $`T/m_{3/2}`$ are displayed in Figs. 6(b).

Thus, high-precision measurements at high energy proton and $`e^+e^{}`$ linear colliders provide access to crucial derivative parameters in string theories.

## 5 Conclusions

In supersymmetric theories stable extrapolations can be performed from the electroweak scale to the Grand Unification scale close to the Planck scale. This feature has been compellingly demonstrated in the evolution of the three gauge couplings to the unification point in the minimal supersymmetric theory.

Such extrapolations are made possible by high-precision measurements of the low-energy parameters. The operation of the $`e^+e^{}`$ colliders LEP and SLC has been crucial in this context. In the near future an enormous extension of this area will be possible if measurements at LHC and prospective $`e^+e^{}`$ linear colliders are combined to draw, if realized in Nature, a comprehensive high-precision picture of the supersymmetric particles and their interactions. Based merely on measurements at low energies, the parameters of the theory can be evolved to high scales by means of renormalization group techniques.

Supersymmetric theories and their breaking mechanisms have the simplest structures and the greatest regularities at high scales. Extrapolations to high scales are therefore crucial to uncover the regularities. The bottom-up approach in the extrapolation of parameters measured at low scales to the high scales provides the most transparent picture. In this way the basis of the SUSY breaking mechanism can be explored and the crucial elements of the fundamental supersymmetric theory can be reconstructed. The method can thus be used to explore particle physics phenomena at a scale where, eventually, particle physics is linked to gravity.

Apart from other examples, we have focused on two interesting scenarios in this approach. The universality of gaugino and scalar mass parameters in minimal supergravity can be demonstrated very clearly if realized in the supersymmetric theory. Small deviations from universality, on the other hand, may be exploited to measure the fundamental parameters in superstring effective field theories, i.e. the strength of dilaton and moduli fields, their mixing and the modular weights. In this way, high-precision extrapolations of gauge and supersymmetric parameters can establish direct contact between superstring theory and experiment.

Many more refinements of the theoretical calculations and future experimental analyses will be necessary to expand the picture we have drawn in this first attempt. However, the prospect of exploring elements of the ultimate unification of the interactions provides a strong impetus to this direction.

## Acknowledgements

We are very grateful for enlightening discussions with J. Kalinowski, G.L. Kane, B.D. Nelson and H.P. Nilles. G.A.B. is grateful for the hospitality extended to him in long-term visits to CERN and DESY, for support from the British Council ARC programme and Intas grant 00-00679. W.P. is supported by the Erwin Schrödinger fellowship No. J2095 of the ‘Fonds zur Förderung der wissenschaftlichen Forschung’ of Austria FWF and partly by the Swiss ‘Nationalfonds’.

## Appendix A One–loop RGEs

In this first appendix we collect the one–loop renormalization group equations (RGEs) including right-handed neutrinos.

Using the notation for the gauge and Yukawa couplings

$`\alpha _i={\displaystyle \frac{g_i^2}{16\pi ^2}},i=1,2,3;Y_k={\displaystyle \frac{y_k^2}{16\pi ^2}},k=t,b,\tau ,\nu `$ (49)
the one-loop RG equations can be written as

$`\dot{\alpha }_i`$ $`=`$ $`b_i\alpha _i^2,`$ (50)
$`\dot{Y}_k`$ $`=`$ $`Y_k({\displaystyle \underset{i}{}}c_{ki}\alpha _i{\displaystyle \underset{l}{}}a_{kl}Y_l),`$ (51)
where the dot denotes the derivative with respect to $`t=\mathrm{log}M_{GUT}^2/Q^2`$, and

$`b_i`$ $`=`$ $`\{33/5,1,3\},`$ (52)
$`c_{ti}`$ $`=`$ $`\{13/15,3,16/3\},c_{bi}=\{7/15,3,16/3\},`$ (53)
$`c_{\tau i}`$ $`=`$ $`\{9/5,3,0\},c_{\nu i}=\{3/5,3,0\},`$ (54)
$`a_{tl}`$ $`=`$ $`\{6,1,0,1\},a_{bl}=\{1,6,1,0\},`$ (55)
$`a_{\tau l}`$ $`=`$ $`\{0,3,4,1\},a_{\nu l}=\{3,0,1,4\},`$ (56)
while the RGEs for the gaugino mass parameters and the $`A`$-parameters read

$`\dot{M}_i`$ $`=`$ $`b_i\alpha _iM_i,`$ (57)
$`\dot{A}_k`$ $`=`$ $`{\displaystyle \underset{i}{}}c_{ki}\alpha _iM_i{\displaystyle \underset{l}{}}a_{kl}A_l.`$ (58)
The RGEs for the soft SUSY breaking mass parameters of the third generation and the Higgs mass parameters are given by:

$`\dot{M}_{L_3}`$ $`=`$ $`2Y_\tau X_\tau 2Y_\nu X_\nu +{\displaystyle \frac{6}{5}}\alpha _1M_1^2+6\alpha _2M_2^2+{\displaystyle \frac{3}{5}}S,`$ (59)
$`\dot{M}_{\nu _{R3}}`$ $`=`$ $`4Y_\nu X_\nu ,`$ (60)
$`\dot{M}_{E_3}`$ $`=`$ $`4Y_\tau X_\tau +{\displaystyle \frac{24}{5}}\alpha _1M_1^2{\displaystyle \frac{6}{5}}S,`$ (61)
$`\dot{M}_{Q_3}`$ $`=`$ $`2Y_bX_b2Y_tX_t+{\displaystyle \frac{2}{15}}\alpha _1M_1^2+6\alpha _2M_2^2+{\displaystyle \frac{16}{3}}\alpha _3M_3^2{\displaystyle \frac{1}{5}}S,`$ (62)
$`\dot{M}_{U_3}`$ $`=`$ $`4Y_tX_t+{\displaystyle \frac{32}{15}}\alpha _1M_1^2+{\displaystyle \frac{16}{3}}\alpha _3M_3^2+{\displaystyle \frac{4}{5}}S,`$ (63)
$`\dot{M}_{D_3}`$ $`=`$ $`4Y_bX_b+{\displaystyle \frac{8}{15}}\alpha _1M_1^2+{\displaystyle \frac{16}{3}}\alpha _3M_3^2{\displaystyle \frac{2}{5}}S,`$ (64)
$`\dot{M}_{H_1}`$ $`=`$ $`6Y_bX_b2Y_\tau X_\tau +{\displaystyle \frac{6}{5}}\alpha _1M_1^2+6\alpha _2M_2^2+{\displaystyle \frac{3}{5}}S,`$ (65)
$`\dot{M}_{H_2}`$ $`=`$ $`6Y_tX_t2Y_\nu X_\nu +{\displaystyle \frac{6}{5}}\alpha _1M_1^2+6\alpha _2M_2^2{\displaystyle \frac{3}{5}}S,`$ (66)
with

$`X_t`$ $`=`$ $`M_{Q_3}^2+M_{U_3}^2+M_{H_2}^2+A_t^2,`$ (67)
$`X_b`$ $`=`$ $`M_{Q_3}^2+m_D^2+M_{H_1}^2+A_b^2,`$ (68)
$`X_\tau `$ $`=`$ $`M_{L_3}^2+M_{E_3}^2+M_{H_1}^2+A_\tau ^2,`$ (69)
$`X_\nu `$ $`=`$ $`M_{L_3}^2+M_{\nu _{R3}}^2+M_{H_2}^2+A_\nu ^2,`$ (70)
$`S`$ $`=`$ $`M_{H_2}^2M_{H_2}^2+{\displaystyle \underset{i=1}{\overset{3}{}}}\left(M_{Q_i}^2M_{L_i}^22M_{U_i}^2+M_{D_i}^2+M_{E_i}^2\right).`$ (71)
The evolution equations for the first two generations are obtained by replacing appropriately the corresponding parameters and Yukawa couplings.

## Appendix B Solutions of the one–loop RGEs

In the following subsections we present the analytical solutions to the 1-loop RGEs including Yukawa couplings using the procedure of Ref. . We also include the generic trace term $`S`$ (see Eq.(71)) in the solutions which had been neglected in Ref. . In this appendix we mark all quantities defined at the the GUT-scale $`M_{GUT}`$ with a subscript G.

### B.1 mSUGRA Boundary Conditions at the GUT-scale

The solutions for the case of the MSSM are summarized first for proper reference. The solution for the gauge couplings and Yukawa couplings are given by:

$`\alpha _i(t)`$ $`=`$ $`{\displaystyle \frac{\alpha _{i,G}}{1+b_i\alpha _{i,G}t}}`$ (72)
$`Y_k(t)`$ $`=`$ $`{\displaystyle \frac{Y_{k,G}u_k}{1+a_{kk}Y_{k,G}_0^tu_k}},`$ (73)
where the functions $`u_k`$ obey the integral system of equations

$`u_t`$ $`=`$ $`{\displaystyle \frac{E_t}{(1+6Y_{b,G}_0^tu_b)^{1/6}}},`$ (74)
$`u_b`$ $`=`$ $`{\displaystyle \frac{E_b}{(1+6Y_{t,G}_0^tu_t)^{1/6}(1+4Y_{\tau ,G}_0^tu_\tau )^{1/4}}},`$ (75)
$`u_\tau `$ $`=`$ $`{\displaystyle \frac{E_\tau }{(1+6Y_{b,G}_0^tu_b)^{1/2}}},`$ (76)
and the functions $`E_k`$ denote the products

$$E_k=\underset{i=1}{\overset{3}{}}(1+b_i\alpha _{i,G}t)^{c_{ki}/b_i}.$$
(77)
The system of integral equations can be solved iteratively and a discussion on the convergence can be found in Ref. .

The gaugino mass parameters and the $`A_k`$ parameters are given by

$`M_i(t)`$ $`=`$ $`{\displaystyle \frac{M_{i,G}}{1+b_i\alpha _{i,G}t}}={\displaystyle \frac{\alpha _i(t)}{\alpha _{i,G}}}M_{i,G},`$ (78)
$`A_k`$ $`=`$ $`e_k+{\displaystyle \frac{A_{k,G}/Y_{k,G}+a_{kk}u_ke_k}{1/Y_{k,G}+a_{kk}u_k}},`$ (79)
with the coefficients

$`e_t`$ $`=`$ $`\stackrel{~}{F}_t+{\displaystyle \frac{A_{b,G}u_bu_be_b}{1/Y_{b,G}+6u_b}},`$ (80)
$`e_b`$ $`=`$ $`\stackrel{~}{F}_b+{\displaystyle \frac{A_{t,G}u_tu_te_t}{1/Y_{t,G}+6u_t}}+{\displaystyle \frac{A_{\tau ,G}u_\tau u_\tau e_\tau }{1/Y_{\tau ,G}+4u_\tau }},`$ (81)
$`e_\tau `$ $`=`$ $`\stackrel{~}{F}_\tau +3{\displaystyle \frac{A_{b,G}u_bu_be_b}{1/Y_{b,G}+6u_b}},`$ (82)
$`\stackrel{~}{F}_k`$ $`=`$ $`t{\displaystyle \underset{i=1}{\overset{3}{}}}c_{ki}M_{i,G}\alpha _i(t).`$ (83)
The mass parameters of the first two generations ($`k=1,2`$) can be expressed as

$`M_{L_k}^2(t)`$ $`=`$ $`M_{L_k,G}^2+{\displaystyle \frac{3}{2}}f_2(t)+{\displaystyle \frac{3}{10}}f_1(t)+{\displaystyle \frac{3}{5}}S^{}(t)`$ (84)
$`M_{E_k}^2(t)`$ $`=`$ $`M_{E_k,G}^2+{\displaystyle \frac{6}{5}}f_1(t){\displaystyle \frac{6}{5}}S^{}(t)`$ (85)
$`M_{Q_k}^2(t)`$ $`=`$ $`M_{Q_k,G}^2+{\displaystyle \frac{8}{3}}f_3(t)+{\displaystyle \frac{3}{2}}f_2(t)+{\displaystyle \frac{1}{30}}f_1(t){\displaystyle \frac{1}{5}}S^{}(t)`$ (86)
$`M_{U_k}^2(t)`$ $`=`$ $`M_{U_k,G}^2+{\displaystyle \frac{8}{3}}f_3(t)+{\displaystyle \frac{8}{15}}f_1(t)+{\displaystyle \frac{4}{5}}S^{}(t)`$ (87)
$`M_{D_k}^2(t)`$ $`=`$ $`M_{D_k,G}^2+{\displaystyle \frac{8}{3}}f_3(t)+{\displaystyle \frac{2}{15}}f_1(t){\displaystyle \frac{2}{5}}S^{}(t)`$ (88)
with

$`f_i(t)`$ $`=`$ $`{\displaystyle \frac{M_{i,G}^2}{b_i}}\left(1{\displaystyle \frac{1}{(1+\alpha _{i,G}b_it)^2}}\right)`$ (89)
$`S^{}(t)`$ $`=`$ $`{\displaystyle \frac{1}{2b_1}}\left(S(t)S(M_U)\right)`$ (90)
$`S(t)`$ $`=`$ $`S(M_U)(1+\beta _1t)^2`$ (91)
$`S(M_U)`$ $`=`$ $`M_{H_1,G}^2M_{H_2,G}^2+{\displaystyle \underset{i=1}{\overset{3}{}}}\left(M_{Q_i,G}^2M_{L_i,G}^22M_{U_i,G}^2+M_{D_i,G}^2+M_{E_i,G}^2\right),`$ (92)
in agreement with Ref. . The mass parameters for the third generation and the Higgs mass parameters are involved owing to the Yukawa couplings:

$`M_{L_3}^2`$ $`=`$ $`M_{L_3,G}^2+{\displaystyle \frac{80f_3+123f_2103/5f_1}{122}}{\displaystyle \frac{3}{5}}S^{}(t)+{\displaystyle \frac{3\mathrm{\Delta }\mathrm{\Sigma }_t18\mathrm{\Delta }\mathrm{\Sigma }_b+35\mathrm{\Delta }\mathrm{\Sigma }_\tau }{122}},`$ (93)
$`M_{E_3}^2`$ $`=`$ $`M_{E_3,G}^2+{\displaystyle \frac{80f_360f_2+16f_1}{61}}+{\displaystyle \frac{6}{5}}S^{}(t)+{\displaystyle \frac{3\mathrm{\Delta }\mathrm{\Sigma }_t18\mathrm{\Delta }\mathrm{\Sigma }_b+35\mathrm{\Delta }\mathrm{\Sigma }_\tau }{61}},`$ (94)
$`M_{Q_3}^2`$ $`=`$ $`M_{Q_3,G}^2+{\displaystyle \frac{128f_3+87f_211f_1}{122}}+{\displaystyle \frac{1}{5}}S^{}(t)+{\displaystyle \frac{17\mathrm{\Delta }\mathrm{\Sigma }_t+20\mathrm{\Delta }\mathrm{\Sigma }_b5\mathrm{\Delta }\mathrm{\Sigma }_\tau }{122}},`$ (95)
$`M_{U_3}^2`$ $`=`$ $`M_{U_3,G}^2+{\displaystyle \frac{72f_354f_2+72/5f_1}{61}}{\displaystyle \frac{4}{5}}S^{}(t)+{\displaystyle \frac{21\mathrm{\Delta }\mathrm{\Sigma }_t4\mathrm{\Delta }\mathrm{\Sigma }_b+\mathrm{\Delta }\mathrm{\Sigma }_\tau }{61}},`$ (96)
$`M_{D_3}^2`$ $`=`$ $`M_{D_3,G}^2+{\displaystyle \frac{56f_342f_2+56/5f_1}{61}}+{\displaystyle \frac{2}{5}}S^{}(t)+{\displaystyle \frac{4\mathrm{\Delta }\mathrm{\Sigma }_t+24\mathrm{\Delta }\mathrm{\Sigma }_b6\mathrm{\Delta }\mathrm{\Sigma }_\tau }{61}},`$ (97)
$`M_{H_1}^2`$ $`=`$ $`M_{H_1,G}^2+{\displaystyle \frac{240f_33f_257/5f_1}{122}}{\displaystyle \frac{3}{5}}S^{}(t)+{\displaystyle \frac{9\mathrm{\Delta }\mathrm{\Sigma }_t+54\mathrm{\Delta }\mathrm{\Sigma }_b+17\mathrm{\Delta }\mathrm{\Sigma }_\tau }{122}},`$ (98)
$`M_{H_2}^2`$ $`=`$ $`M_{H_2,G}^2+{\displaystyle \frac{272f_3+21f_289/5f_1}{122}}+{\displaystyle \frac{3}{5}}S^{}(t)+{\displaystyle \frac{63\mathrm{\Delta }\mathrm{\Sigma }_t12\mathrm{\Delta }\mathrm{\Sigma }_b+3\mathrm{\Delta }\mathrm{\Sigma }_\tau }{122}},`$ (99)
with

$`\mathrm{\Delta }\mathrm{\Sigma }_k`$ $`=`$ $`\mathrm{\Sigma }_k(t)\mathrm{\Sigma }_{k,G},`$ (100)
$`\mathrm{\Sigma }_t`$ $`=`$ $`M_{Q_3}^2+M_{U_3}^2+M_{H_2}^2,`$ (101)
$`\mathrm{\Sigma }_b`$ $`=`$ $`M_{Q_3}^2+M_{D_3}^2+M_{H1}^2,`$ (102)
$`\mathrm{\Sigma }_\tau `$ $`=`$ $`M_{L_3}^2+M_{E_3}^2+M_{H1}^2.`$ (103)
The explicit solution for $`\mathrm{\Sigma }_k`$ read as:

$`\mathrm{\Sigma }_k`$ $`=`$ $`\xi _k+A_k^2+2e_kA_k{\displaystyle \frac{A_{k,G}^2/Y_{k,G}\mathrm{\Sigma }_{k,G}/Y_{k,G}+a_{kk}u_k\xi _k}{1/Y_{k,G}+a_{kk}u_k}},`$ (104)
with

$`\xi _t`$ $`=`$ $`\stackrel{~}{E}_t+2\stackrel{~}{F}_t{\displaystyle \frac{A_{b,G}u_bu_be_b}{1/Y_{b,G}+6u_b}}+7\left({\displaystyle \frac{A_{b,G}u_bu_be_b}{1/Y_{b,G}+6u_b}}\right)^2`$ (105)
$``$ $`{\displaystyle \frac{(\mathrm{\Sigma }_{b,G}+A_{b,G}^2)u_b2A_{b,G}u_be_b+u_b\xi _b}{1/Y_{b,G}+6u_b}},`$
$`\xi _b`$ $`=`$ $`\stackrel{~}{E}_b+2\stackrel{~}{F}_b\left[{\displaystyle \frac{A_{t,G}u_tu_te_t}{1/Y_{t,G}+6u_t}}+{\displaystyle \frac{A_{\tau ,G}u_\tau u_\tau e_\tau }{1/Y_{\tau ,G}+4u_\tau }}\right]`$ (106)
$`+`$ $`7\left({\displaystyle \frac{A_{t,G}u_tu_te_t}{1/Y_{t,G}+6u_t}}\right)^2+5\left({\displaystyle \frac{A_{\tau ,G}u_\tau u_\tau e_\tau }{1/Y_{\tau ,G}+4u_\tau }}\right)^2`$
$`+`$ $`2\left({\displaystyle \frac{A_{t,G}u_tu_te_t}{1/Y_{t,G}+6u_t}}\right)\left({\displaystyle \frac{A_{\tau ,G}u_\tau u_\tau e_\tau }{1/Y_{\tau ,G}+4u_\tau }}\right){\displaystyle \frac{(\mathrm{\Sigma }_{t,G}+A_{t,G}^2)u_t2A_{t,G}u_te_t+u_t\xi _t}{1/Y_{t,G}+6u_t}}`$
$``$ $`{\displaystyle \frac{(\mathrm{\Sigma }_{\tau ,G}+A_{\tau ,G}^2)u_\tau 2A_{\tau ,G}u_\tau e_\tau +u_\tau \xi _\tau }{1/Y_{\tau ,G}+4u_\tau }},`$
$`\xi _\tau `$ $`=`$ $`\stackrel{~}{E}_\tau +6\stackrel{~}{F}_\tau {\displaystyle \frac{A_{b,G}u_bu_be_b}{1/Y_{b,G}+6u_b}}+27\left({\displaystyle \frac{A_{b,G}u_bu_be_b}{1/Y_{b,G}+6u_b}}\right)^2`$ (107)
$``$ $`3{\displaystyle \frac{(\mathrm{\Sigma }_{b,G}+A_{b,G}^2)u_b2A_{b,G}u_be_b+u_b\xi _b}{1/Y_{b,G}+6u_b}}`$
$`\stackrel{~}{E}_k`$ $`=`$ $`t^2\left({\displaystyle \underset{i=1}{\overset{3}{}}}c_{ki}\alpha _iM_{i,G}\right)^2+2t{\displaystyle \underset{i=1}{\overset{3}{}}}c_{ki}\alpha _iM_{i,G}^2t^2{\displaystyle \underset{i=1}{\overset{3}{}}}c_{ki}b_i\alpha _i^2M_{i,G}^2.`$ (108)
Finally we express $`t_Z=\mathrm{log}(M_{GUT}^2/m_Z^2)`$ and $`\alpha _G`$ in terms of observables at the electroweak scale, using Eq.(72), by

$`t_Z={\displaystyle \frac{4\pi }{(b_1b_2)\alpha (m_Z)}}\left({\displaystyle \frac{3\mathrm{cos}^2\vartheta _W}{5}}\mathrm{sin}^2\vartheta _W\right)`$ (109)
and similarly for the gauge coupling at the GUT scale:

$`\alpha _G={\displaystyle \frac{5\alpha (m_Z)}{3}}{\displaystyle \frac{b_1b_2}{\frac{5}{3}b_1\mathrm{sin}^2\mathrm{\Theta }_Wb_2\mathrm{cos}^2\mathrm{\Theta }_W}}`$ (110)
### B.2 Universal SUGRA Boundary Conditions at the GUT-scale including Right-handed Neutrinos

Those formulae are noted in this subsection which are changed compared to the previous section in the range between $`M_U`$ and $`M_{\nu _R}`$. Below $`M_{\nu _R}`$, these quantities have the same form as given above. In addition we note also the equations related to the right-handed neutrinos,

$`u_t`$ $`=`$ $`{\displaystyle \frac{E_t}{(1+6Y_{b,G}_0^tu_b)^{1/6}(1+4Y_{\nu ,G}_0^tu_\nu )^{1/4}}},`$ (111)
$`u_\tau `$ $`=`$ $`{\displaystyle \frac{E_\tau }{(1+6Y_{b,G}_0^tu_b)^{1/2}(1+4Y_{\nu ,G}_0^tu_\nu )^{1/4}}},`$ (112)
$`u_\nu `$ $`=`$ $`{\displaystyle \frac{E_\nu }{(1+6Y_{t,G}_0^tu_t)^{1/2}(1+4Y_{\tau ,G}_0^tu_\tau )^{1/4}}},`$ (113)
$`e_t`$ $`=`$ $`\stackrel{~}{F}_t+{\displaystyle \frac{A_{b,G}u_bu_be_b}{1/Y_{b,G}+6u_b}}+{\displaystyle \frac{A_{\nu ,G}u_\nu u_\nu e_\nu }{1/Y_{\nu ,G}+4u_\nu }},`$ (114)
$`e_\tau `$ $`=`$ $`\stackrel{~}{F}_\tau +3{\displaystyle \frac{A_{b,G}u_bu_be_b}{1/Y_{b,G}+6u_b}}+{\displaystyle \frac{A_{\nu ,G}u_\nu u_\nu e_\nu }{1/Y_{\nu ,G}+4u_\nu }},`$ (115)
$`e_\nu `$ $`=`$ $`\stackrel{~}{F}_\nu +3{\displaystyle \frac{A_{t,G}u_tu_te_t}{1/Y_{t,G}+6u_t}}+{\displaystyle \frac{A_{\tau ,G}u_\tau u_\tau e_\tau }{1/Y_{\tau ,G}+4u_\tau }}.`$ (116)
$`\xi _t`$ $`=`$ $`\stackrel{~}{E}_t+2\stackrel{~}{F}_t\left({\displaystyle \frac{A_{b,G}u_bu_be_b}{1/Y_{b,G}+6u_b}}+{\displaystyle \frac{A_{\nu ,G}u_\nu u_\nu e_\nu }{1/Y_{\nu ,G}+4u_\nu }}\right)+7\left({\displaystyle \frac{A_{b,G}u_bu_be_b}{1/Y_{b,G}+6u_b}}\right)^2`$ (117)
$`+`$ $`5\left({\displaystyle \frac{A_{\nu ,G}u_\nu u_\nu e_\nu }{1/Y_{\nu ,G}+4u_\nu }}\right)^2+2\left({\displaystyle \frac{A_{b,G}u_bu_be_b}{1/Y_{b,G}+6u_b}}\right)\left({\displaystyle \frac{A_{\nu ,G}u_\nu u_\nu e_\nu }{1/Y_{\nu ,G}+4u_\nu }}\right)`$
$``$ $`{\displaystyle \frac{(\mathrm{\Sigma }_{b,G}+A_{b,G}^2)u_b2A_{b,G}u_be_b+u_b\xi _b}{1/Y_{b,G}+6u_b}},`$
$``$ $`{\displaystyle \frac{(\mathrm{\Sigma }_{\nu ,G}+A_{\nu ,G}^2)u_\nu 2A_{\nu ,G}u_\nu e_\nu +u_\nu \xi _\nu }{1/Y_{\nu ,G}+4u_\nu }},`$
$`\xi _\tau `$ $`=`$ $`\stackrel{~}{E}_\tau +2\stackrel{~}{F}_\tau \left(3{\displaystyle \frac{A_{b,G}u_bu_be_b}{1/Y_{b,G}+6u_b}}+{\displaystyle \frac{A_{\nu ,G}u_\nu u_\nu e_\nu }{1/Y_{\nu ,G}+4u_\nu }}\right)`$ (118)
$`+`$ $`27\left({\displaystyle \frac{A_{b,G}u_bu_be_b}{1/Y_{b,G}+6u_b}}\right)^2+5\left({\displaystyle \frac{A_{\nu ,G}u_\nu u_\nu e_\nu }{1/Y_{\nu ,G}+4u_\nu }}\right)^2`$
$`+`$ $`6\left({\displaystyle \frac{A_{b,G}u_bu_be_b}{1/Y_{b,G}+6u_b}}\right)\left({\displaystyle \frac{A_{\nu ,G}u_\nu u_\nu e_\nu }{1/Y_{\nu ,G}+4u_\nu }}\right)`$
$``$ $`3{\displaystyle \frac{(\mathrm{\Sigma }_{b,G}+A_{b,G}^2)u_b2A_{b,G}u_be_b+u_b\xi _b}{1/Y_{b,G}+6u_b}}`$
$``$ $`{\displaystyle \frac{(\mathrm{\Sigma }_{\nu ,G}+A_{\nu ,G}^2)u_\nu 2A_{\nu ,G}u_\nu e_\nu +u_\nu \xi _\nu }{1/Y_{\nu ,G}+4u_\nu }},`$
$`\xi _\nu `$ $`=`$ $`\stackrel{~}{E}_\nu +2\stackrel{~}{F}_\nu \left(3{\displaystyle \frac{A_{t,G}u_tu_te_t}{1/Y_{t,G}+6u_t}}+{\displaystyle \frac{A_{\tau ,G}u_\tau u_\tau e_\tau }{1/Y_{\tau ,G}+4u_\tau }}\right)`$ (119)
$`+`$ $`27\left({\displaystyle \frac{A_{t,G}u_tu_te_t}{1/Y_{t,G}+6u_t}}\right)^2+5\left({\displaystyle \frac{A_{\tau ,G}u_\tau u_\tau e_\tau }{1/Y_{\tau ,G}+4u_\tau }}\right)^2`$
$`+`$ $`6\left({\displaystyle \frac{A_{t,G}u_tu_te_t}{1/Y_{t,G}+6u_t}}\right)\left({\displaystyle \frac{A_{\tau ,G}u_\tau u_\tau e_\tau }{1/Y_{\tau ,G}+4u_\tau }}\right)`$
$``$ $`3{\displaystyle \frac{(\mathrm{\Sigma }_{t,G}+A_{t,G}^2)u_t2A_{t,G}u_te_t+u_t\xi _t}{1/Y_{t,G}+6u_t}}`$
$``$ $`{\displaystyle \frac{(\mathrm{\Sigma }_{\tau ,G}+A_{\tau ,G}^2)u_\tau 2A_{\tau ,G}u_\tau e_\tau +u_\tau \xi _\tau }{1/Y_{\tau ,G}+4u_\tau }}.`$
$`\mathrm{\Sigma }_\nu `$ $`=`$ $`M_{L_3}^2+M_{\nu _{3,R}}^2+M_{H_2}^2.`$ (120)