High energy accelerators may probe into dark matter and the seesaw neutrino mass scales if they are not much heavier than ~O(TeV). In the absence of supersymmetry, we extend a class of SO(10) models to predict well known cold dark matter candidates while achieving precision unification with experimentally testable proton lifetime. The most important prediction is a new radiative seesaw formula of Ma type accessible to accelerator tests while the essential small value of its quartic coupling also emerges naturally. This dominates over the high-scale seesaw contributions making a major impact on neutrino physics and dark matter, opening up high prospects as a theory of fermion masses.

We show how gauge coupling unification is successfully implemented through non-supersymmetric grand unified theory, $SO(10)\times G_f (~G_f=S_4, SO(3)_f, SU(3)_f)$, using low-scale flavor symmetric model of the type $SU(2)_L\times U(1)_Y$ $ \times SU(3)_C \times S_4$ recently proposed by Hagedorn, Lindner, and Mohapatra, while assigning matter-parity discrete symmetry for the dark matter stability. For gauge coupling unification in the single-step breaking case, we show that a color-octet fermion and a hyperchargeless weak-triplet fermionic dark matter are the missing particles needed to complete its MSSM-equivalent degrees of freedom. When these are included the model automatically predicts the nonsupersymmetric grand unification with a scale identical to the minimal supersymmetric standard model/grand unified theory scale. We also find a two-step breaking model with Pati-Salam intermediate symmetry where the dark matter and a low-mass color-octet scalar or the fermion are signaled by grand unification. The proton-lifetime predictions are found to be accessible to ongoing or planned searches in a number of models. We discuss grand unified origin of the light fermionic triplet dark matter, the color-octet fermion, and their phenomenology.

It is well known that in single step breaking of R-parity conserving SUSY SO(10) that needs the Higgs representations $126+\bar 126$ the GUT-gauge coupling violates the perturbative constraint at mass scales few times larger than the GUT scale. Therefore, if the SO(10) gauge coupling is to remain perturbative up the Planck scale(=2x10^{18} GeV), the scale $M_U$ of GUT symmetry breaking is to be bounded from below. The bound depends upon specific Higgs representations used for SO(10) symmetry breaking but, as we find, can not be lower than $1.5x10^{17}$ GeV. In order to obtain such high unification scale we propose a two-step SO(10) breaking through $SU(2)_LXSU(2)_RXU(1)_{B-L}XSU(3)_C(g_{2L} \neq g_{2R})$ intermediate gauge symmetry. We estimate potential threshold and gravitational corrections to the running of gauge couplings and show that they can make the picture of perturbative GUT- gauge coupling running consistent at least up to the Planck scale. We also show that when $SO(10) \to G_{2213}$ by $210 + 54$, gravitational corrections alone with negligible threshold effects may guarantee such perturbative gauge coupling. The lifetime of the proton is found to increase by nearly 6 orders over the current experimental limit for $p \to e^+\pi^0$. For the proton decay mediated by dim.5 operator a wide range of lifetimes is possible extending from the current experimental limit up to values 2-3 orders longer.

Including contributions of scale-dependent vacuum expectation values, we derive new analytic formulas and obtain substantially different numerical predictions for the running masses of quarks and charged-leptons at higher scales in the SM, 2HDM and MSSM. These formulas exhibit significantly different behaviours with respect to their dependence on gauge and Yukawa couplings than those derived earlier. At one-loop level the masses of the first two generations are found to be independent of Yukawa couplings of the third generation in all the three effective theories in the small mixing limit. Analytic formulas are also obtained for running $\tan\beta(\mu)$ in 2HDM and MSSM. Other numerical analyses include study of the third generation masses at high scales as functions of low-energy values of $\tan\beta$ and SUSY scale $M_S=M_Z-10^4$ GeV.