We explore the parameter space of the minimal supersymmetric extension of the Standard Model (MSSM), allowing the soft supersymmetry-breaking masses of the Higgs multiplets, m_{1,2}, to be non-universal (NUHM). Compared with the constrained MSSM (CMSSM) in which m_{1,2} are required to be equal to the soft supersymmetry-breaking masses m_0 of the squark and slepton masses, the Higgs mixing parameter mu and the pseudoscalar Higgs mass m_A, which are calculated in the CMSSM, are free in the NUHM model. We incorporate accelerator and dark matter constraints in determining allowed regions of the (mu, m_A), (mu, M_2) and (m_{1/2}, m_0) planes for selected choices of the other NUHM parameters. In the examples studied, we find that the LSP mass cannot be reduced far below its limit in the CMSSM, whereas m_A may be as small as allowed by LEP for large tan \beta. We present in Appendices details of the calculations of neutralino-slepton, chargino-slepton and neutralino-sneutrino coannihilation needed in our exploration of the NUHM.

We extend previous combinations of LEP and cosmological relic density constraints on the parameter space of the constrained MSSM, with universal input supersymmetry-breaking parameters, to large tan beta. We take account of the possibility that the lightest Higgs boson might weigh about 115 GeV, but also retain the possibility that it might be heavier. We include the most recent implementation of the b to s gamma constraint at large tan beta. We refine previous relic density calculations at large tan beta by combining a careful treatment of the direct-channel Higgs poles in annihilation of pairs of neutralinos chi with a complete treatment of chi - stau coannihilation, and discuss carefully uncertainties associated with the mass of the b quark. We find that coannihilation and pole annihilations allow the CMSSM to yield an acceptable relic density at large tan beta, but it is consistent with all the constraints only if m_chi > 140 (180) GeV for mu > 0 (mu < 0) for our default choices m_b(m_b) = 4.25 GeV, m_t = 175 GeV, and A_0 = 0.

Assuming that the lightest supersymmetric particle (LSP) is the lightest neutralino, we present a detailed exploration of neutralino-stau coannihilation channels, including analytical expressions and numerical results. We also include neutralino coannihilations with the selectron and smuon. We evaluate the implications of coannihilations for the cosmological relic density of the LSP, which is assumed to be stable, in the constrained minimal supersymmetric extension of the Standard Model (CMSSM), in which the soft supersymmetry-breaking parameters are universal at the supergravity GUT scale. We evaluate the changes due to coannihilations in the region of the MSSM parameter space that is consistent with the cosmological upper limit on the relic LSP density. In particular, we find that the upper limit on the mass of the neutralino is increased from about 200 GeV to about 600 GeV in the CMSSM, and estimate a qualitatively similar increase for gauginos in the general MSSM.

We analyze the accelerator constraints on the parameter space of the Minimal Supersymmetric extension of the Standard Model, comparing those now available from LEP II and anticipating the likely sensitivity of Tevatron Run II. The most important limits are those from searches for charginos, neutralinos and Higgs bosons at LEP, and searches for stop squarks, charginos and neutralinos at the Tevatron Collider. We also incorporate the constraints derived from b --> s + gamma decay, and discuss the relevance of charge- and colour-breaking minima in the effective potential. We combine and compare the different constraints on the Higgs-mixing parameter mu, the gaugino-mass parameter m_{1/2} and the scalar-mass parameter m0, incorporating radiative corrections to the physical particle masses. We focus on the resulting limitations on supersymmetric dark matter, assumed to be the lightest neutralino, incorporating coannihilation effects in the calculation of the relic abundance. We find that m_chi > 51 GeV and tan(beta) > 2.2 if all soft supersymmetry-breaking scalar masses are universal, including those of the Higgs bosons, and that these limits weaken to m_chi > 46 GeV and tan(beta) > 1.9 if non-universal scalar masses are allowed. Light neutralino dark matter cannot be primarily Higgsino in composition.

New CP violating phases in the MSSM can affect both the abundance and detection of neutralino dark matter. We discuss the effect of including cosmological constraints in the limits on new sources of CP violation in the MSSM and the effects of new CP violating parameters on dark matter densities and detection.

I discuss recent developments in the study of cosmological limits on the Minimal Supersymmetric Standard Model (MSSM). In particular, I focus on the effect of neutralino-stau coannihilation on the relic abundance of neutralinos, and I give examples where the cosmologically derived limits on the supersymmetric parameters are relaxed, and one example (CP violating phases) where they are not.

We compute the four-Fermi neutralino-quark interaction Lagrangian including contributions from the CP violating phases in the MSSM. We find that neutralino-nucleus scattering cross-sections relevant for direct detection experiments show a strong dependence on the value of the CP-violating phase associated with the $\mu$ parameter, $\theta_\mu$. In some cases, for a broad range of non-zero $\theta_\mu$, there are cancellations in the cross-sections which reduce the cross-section by more than an order of magnitude. In other cases, there may be enhancements as one varies $\theta_\mu$.

We consider the effects of neutralino-stau coannihilations on the cosmological relic density of the lightest supersymmetric particle (LSP) $\chi$ in the minimal supersymmetric extension of the Standard Model (MSSM), particularly in the constrained MSSM in which universal supergravity inputs at the GUT scale are assumed. For much of the parameter space in these models, $\chi$ is approximately a U(1) gaugino ${\tilde B}$, and constraints on the cosmological relic density $\Omega_{\tilde B} h^2$ yield an upper bound on the bino mass. We show that in regions of parameter space for which the cosmological bound is nearly saturated, coannihilations of the $\tilde B$ with the stau, the next lightest sparticle, are important and may reduce significantly the bino relic density. Including also bino coannihilations with the selectrons and smuons, we find that the upper limit on $m_{\chi}$ is increased from about 200 GeV to about 600 GeV in the constrained MSSM, with a similar new upper limit expected in the MSSM.

We analyze the LEP constraints from searches for charginos $\chi^{\pm}$ and neutralinos $\chi_i$, taking into account radiative corrections to the relations between their masses and the underlying Higgs-mixing and gaugino-mass parameters $\mu, m_{1/2}$ and the trilinear mass parameter $A_t$. Whilst radiative corrections do not alter the excluded domain in $m_{\chi^{\pm}}$ as a function of $m_{\chi^{\pm}} - m_{\chi}$, its mapping into the $\mu, m_{1/2}$ plane is altered. We update our previous lower limits on the mass of gaugino dark matter and on tan$\beta$, the ratio of Higgs vacuum expectation values, in the light of the latest LEP data and these radiative corrections. We also discuss the viability of Higgsino dark matter, incorporating co-annihilation effects into the calculation of the Higgsino relic abundance. We find that Higgsino dark matter is viable for only a very limited range of $\mu$ and $m_{1/2}$, which will be explored completely by upcoming LEP runs.

We examine the current state of neutralino dark matter and consider how the LEP constraints on the Minimal Supersymmetric Standard Model parameters are squeezing the available dark matter regions. We also show how cosmological constraints augment bounds coming from collider searches to further constrain the MSSM parameter space.

A significant lower limit on the mass of the lightest neutralino \chi can be obtained by combining the results from sparticle searches at LEP at centre-of-mass energies up to 172 GeV with cosmological considerations, if it is assumed that the \chi is stable. Exclusion domains from slepton searches close m_\chi \sim 0 loopholes that were left open by previous lower-energy LEP searches for charginos and neutralinos, leading to the lower limit $m_\chi \ga 17$ GeV. The constraints on supersymmetric parameter space are strengthened significantly if LEP constraints on supersymmetric Higgs bosons are taken into account, and further if the relic neutralino density is required to fall within the range favoured by astrophysics and cosmology. These bounds are considerably strengthened if universality at the GUT scale is assumed for soft supersymmetry-breaking scalar masses, including those of the Higgs bosons. In this case, the Higgs searches play a dramatic role, and we find that m_\chi \ga 40 GeV. Furthermore, we find that if tan\beta \la 1.7 for \mu<0, or tan\beta \la 1.4 for \mu >0, the cosmological relic density is too large for all values of m_\chi.

We discuss the lower limit on the mass of the neutralino $\chi$ that can be obtained by combining data from $e^+e^-$ annihilation at LEP and elsewhere with astrophysical and theoretical considerations. Loopholes in the purely experimental analysis of ALEPH data from the Z peak and LEP 1.5, which appear when $\mu<0$ for certain values of the sneutrino mass $m_{\tilde\nu}$ and the ratio $\tan\beta$ of supersymmetric Higgs vacuum expectation values, may be largely or totally excluded by data from lower-energy $e^+e^-$ data, the hypothesis that most of the cosmological dark matter consists of $\chi$ particles, and the assumption that electroweak symmetry breaking is triggered by radiative corrections due to a heavy top quark. The combination of these inputs imposes $m_{\chi} \ge 21.4~\gev$, if soft supersymmetry-breaking masses are assumed to be universal at the grand-unification scale.

In the standard inflationary scenario with inflaton potential $V(\Phi)=M^4-{1\over4}\lambda\Phi^4$, the resulting density perturbations $\delta\rho/\rho$ are proportional to $\lambda^{1/2}$. Upper bounds on $\delta\rho/\rho$ require $\lambda < 10^{-13}$. Ratra has shown that an alternative treatment of reheating results in $\delta\rho/\rho \propto \lambda^{-1}$, so that an upper bound on $\delta\rho/\rho$ does not put an obvious upper bound on $\lambda$. We verify that $\delta\rho/\rho \propto \lambda^{-1}$ is indeed a possibility, but show that $\lambda < 10^{-13}$ is still required.

Inflationary models predict a definite, model independent, angular dependence for the three-point correlation function of $\Delta T/T$ at large angles (greater than $\sim 1^\circ$) which we calculate. The overall amplitude is model dependent and generically unobservably small, but may be large in some specific models. We compare our results with other models of nongaussian fluctuations.