It has recently been pointed out that viable sparticle mass spectra can be generated in Yukawa unified $\mathrm{SO}\left(10\right)\mathrm{}$ supersymmetric grand unified models consistent with radiative breaking of electroweak symmetry. Model solutions are obtained only if $\tan \beta \sim 50,$ μ<0 and positive D-term contributions to scalar masses from $\mathrm{SO}\left(10\right)\mathrm{}$ gauge symmetry breaking are used. In this paper, we attempt to systematize the parameter space regions where solutions are obtained. We go on to calculate the relic density of neutralinos as a function of parameter space. No regions of the parameter space explored were actually cosmologically excluded, and very reasonable relic densities were found in much of parameter space. Direct neutralino detection rates could exceed 1 event/kg/day ${\mathrm{for}\mathrm{}\mathrm{a}\mathrm{}}^{73}\mathrm{Ge}$ detector, for low values of GUT scale gaugino mass $m_{1/2}.$ We also calculate the branching fraction for $\overrightarrow{b}s\gamma$ decays, and find that it is beyond the 95% C.L. experimental limits in much, but not all, of the parameter space regions explored. For the Fermilab Tevatron collider, significant regions of parameter space can be explored via $b\overline{b}A$ and $b\overline{b}H$ searches. There also exist some limited regions of parameter space where a trilepton signal can be seen at TeV33. Finally, there exist significant regions of parameter space where direct detection of bottom squark pair production can be made, especially for large negative values of the GUT parameter $A_0.$

• Received 3 May 2000

DOI:https://doi.org/10.1103/PhysRevD.63.015007