A non-relativistic effective field theory (NREFT) offers a bottom-up framework to classify Dark Matter (DM)-nucleon interactions relevant for scattering at direct detection experiments by organizing the interactions in powers of the momentum transfer $\vec{q}$ and DM velocity $\vec{v}$. This approach generates a number of operators including P-odd and T-odd operators; these can only be generated from a relativistic theory with CP violating interactions. We consider the leading order P-odd, T- odd operators viz. $\mathcal{O}_{10}$, $\mathcal{O}_{11}$ and $\mathcal{O}_{12}$ and compare the constraints on these operators from leading direct detection searches and from the bound on the neutron EDM (nEDM). We perform our analysis using simplified models with charged mediators and compute the loop diagrams contributing to the nEDM. We find that constraints on the DM scattering cross section from the bound on the nEDM are several orders of magnitude stronger than the limits from direct searches, and even well below the neutrino floor for such NREFT operators, for the entire sub-GeV to TeV DM mass range. This indicates that these operators need not be considered when analyzing data from present or future direct dark matter detection experiments.

The self coupling $\lambda$ of the Higgs boson in the Standard Model may show critical behavior, i.e. the Higgs potential may have a point at an energy scale $\sim 10^{17-18}$ GeV where both the first and second derivatives (almost) vanish. In this case the Higgs boson can serve as inflaton even if its nonminimal coupling to the curvature scalar is only ${\cal O}(10)$, thereby alleviating concerns about the perturbative unitarity of the theory. We find that just before the Higgs as inflaton enters the flat region of the potential the usual slow--roll conditions are violated. This leads to "overshooting" behavior, which in turn strongly enhances scalar curvature perturbations because of the excitation of entropic (non--adiabatic) perturbations. For appropriate choice of the free parameters these large perturbations occur at length scales relevant for the formation of primordial black holes. Even if these perturbations are not quite large enough to trigger copious black hole formation, they source second order tensor perturbations, i.e. primordial gravitational waves; the corresponding energy density can be detected by the proposed space-based gravitational wave detectors DECIGO and BBO.

We investigate the production of neutralino dark matter in a cosmological scenario featuring an early matter dominated era ending at a relatively low reheating temperature. In such scenarios different production mechanisms of weakly interacting massive particles (WIMPs), besides the well--studied thermal production, can be important. This opens up new regions of parameter space where the lightest neutralino, as the best--known supersymmetric (SUSY) WIMP, obtains the required relic abundance. Many of these new sets of parameters are also compatible with current limits from colliders as well as direct and indirect WIMP searches. In particular, in standard cosmology bino--like neutralinos, which emerge naturally as lightest neutralino in many models, can have the desired relic density only in some finetuned regions of parameter space where the effective annihilation cross section is enhanced by co--annihilation or an $s-$channel pole. In contrast, if the energy density of the universe was dominated by long--lived PeV--scale particles (e.g. moduli or Polonyi fields), bino--like neutralinos can obtain the required relic density over wide regions of supersymmetric parameter space. We identify the interesting ranges of mass and decay properties of the heavy long--lived particles, carefully treating the evolution of the temperature of the thermal background.

We consider quantum corrections to co-annihilation processes of Weakly Interacting Massive Particles (WIMPs) due to the exchange of light bosons in the initial state ("Sommerfeld corrections"). We work at one-loop level, i.e. we assume that these corrections can be treated perturbatively. Co-annihilation is important if there is at least one additional new particle with mass close to the lightest WIMP, which is a Dark Matter candidate. In this case the exchange of a (relatively light) boson in the initial state can change the identity of the annihilating particles. The corrections we are interested in factorize, as in the case of WIMP self-annihilation treated previously, but they can mix different tree-level amplitudes. Moreover, even small mass splittings between the external particles and those in the loop can change the relevant loop functions significantly. We find exact analytical expressions for these functions, and illustrate the effects by considering the cases of wino- or higgsino-like neutralinos as examples.

A broad range of single field models of inflation are analyzed in light of all relevant recent cosmological data, checking whether they can lead to the formation of long--lived Primordial Black Holes (PBHs) as candidate for dark matter. To that end we calculate the spectral index of the power spectrum of primordial perturbations as well as its first and second derivatives. PBH formation is possible only if the spectral index $n_S(k_0)$ increases significantly at small scales. Since current data indicate that the first derivative $\alpha_S$ of the spectral index is negative at the pivot scale, PBH formation is only possible in the presence of a sizable and positive second derivative ("running of the running") $\beta_S$. Among the three small-field and five large-field models we analyze, only the "running-mass" model allows PBH formation, for a narrow range of parameters.

A broad range of single field models of inflation are analyzed in light of all relevant recent cosmological data, checking whether they can lead to the formation of long-lived Primordial Black Holes (PBHs). To that end we calculate the spectral index of the power spectrum of primordial perturbations as well as its first and second derivatives. PBH formation is possible only if the spectral index increases significantly at small scales, i.e. large wave number $k$. Since current data indicate that the first derivative $\alpha_S$ of the spectral index $n_S(k_0)$ is negative at the pivot scale $k_0$, PBH formation is only possible in the presence of a sizable and positive second derivative ("running of the running") $\beta_S$. Among the three small-field and five large-field models we analyze, only one small-field model, the "running mass" model, allows PBH formation, for a narrow range of parameters. We also note that none of the models we analyze can accord for a large and negative value of $\alpha_S$, which is weakly preferred by current data.

We investigate the relic abundance of asymmetric Dark Matter particles that were in thermal equilibrium in the early universe. The standard analytic calculation of the symmetric Dark Matter is generalized to the asymmetric case. We calculate the asymmetry required to explain the observed Dark Matter relic abundance as a function of the annihilation cross section. We show that introducing an asymmetry always reduces the indirect detection signal from WIMP annihilation, although it has a larger annihilation cross section than symmetric Dark Matter. This opens new possibilities for the construction of realistic models of MeV Dark Matter.

We revisit the question whether the running-mass inflation model allows the formation of Primordial Black Holes (PBHs) that are sufficiently long-lived to serve as candidates for Dark Matter. We incorporate recent cosmological data, including the WMAP 7-year results. Moreover, we include "the running of the running" of the spectral index of the power spectrum, as well as the renormalization group "running of the running" of the inflaton mass term. Our analysis indicates that formation of sufficiently heavy, and hence long-lived, PBHs still remains possible in this scenario. As a by-product, we show that the additional term in the inflaton potential still does not allow significant negative running of the spectral index.

We compute one-loop corrections to the annihilation of non--relativistic particles $\chi$ due to the exchange of a (gauge or Higgs) boson $\phi$ with mass $\mu$ in the initial state. In the limit $m_\chi \gg \mu$ this leads to the "Sommerfeld enhancement" of the annihilation cross section. However, here we are interested in the case $\mu \lsim m_\chi$, where the one--loop corrections are well--behaved, but can still be sizable. We find simple and accurate expressions for annihilation from both $S-$ and $P-$wave initial states; they differ from each other if $\mu \neq 0$. In order to apply our results to the calculation of the relic density of Weakly Interacting Massive Particles (WIMPs), we describe how to compute the thermal average of the corrected cross sections. We apply this formalism to scalar and Dirac fermion singlet WIMPs, and show that the corrections are always very small in the former case, but can be very large in the latter. Moreover, in the context of the Minimal Supersymmetric Standard Model, these corrections can decrease the relic density of neutralinos by more than 1%, if the lightest neutralino is a strongly mixed state.

Approximate analytical solutions of the Boltzmann equation for particles that are either extremely relativistic or non-relativistic when they decouple from the thermal bath are well established. However, no analytical formula for the relic density of particles that are semi-relativistic at decoupling is yet known. We propose a new ansatz for the thermal average of the annihilation cross sections for such particles, and find a semi-analytical treatment for calculating their relic densities. As examples, we consider Majorana- and Dirac-type neutrinos. We show that such semi-relativistic relics cannot be good cold Dark Matter candidates. However, late decays of meta-stable semi-relativistic relics might have released a large amount of entropy, thereby diluting the density of other, unwanted relics.

In this talk we present data analysis methods for reconstructing the mass and couplings of Weakly Interacting Massive Particles (WIMPs) by using directly future experimental data (i.e., measured recoil energies) from direct Dark Matter detection. These methods are independent of the model of Galactic halo as well as of WIMPs. The basic ideas of these methods and the feasibility and uncertainties of applying them to direct detection experiments with the next generation detectors will be discussed.

Weakly Interacting Massive Particles (WIMPs) are one of the leading candidates for Dark Matter. For understanding the properties of WIMPs and identifying them among new particles produced at colliders (hopefully in the near future), determinations of their mass and their couplings on nucleons from direct Dark Matter detection experiments are essential. Based on our method for determining the WIMP mass model-independently from experimental data, we present a way to also estimate the spin-independent (SI) WIMP-nucleon coupling by using measured recoil energies directly. This method isindependent of the as yet unknown velocity distribution of halo WIMPs. In spite of the uncertainty of the local WIMP density (of a factor of ~ 2), at least an upper limit on the SI WIMP-nucleon coupling could be given, once two (or more) experiments with different target nuclei obtain positive signals. In a background-free environment, for a WIMP mass of 100 GeV its SI coupling on nucleons could in principle be estimated with a statistical error of only ~ 15% with just 50 events from each experiment.

Weakly Interacting Massive Particles (WIMPs) are one of the leading candidates for Dark Matter. We developed a model-independent method for determining the WIMP mass by using data (i.e., measured recoil energies) of direct detection experiments. Our method is independent of the as yet unknown WIMP density near the Earth, of the form of the WIMP velocity distribution, as well as of the WIMP-nucleus cross section. It requires however positive signals from at least two detectors with different target nuclei. At the first phase of this work we found a systematic deviation of the reconstructed WIMP mass from the real one for heavy WIMPs. Now we improved this method so that this deviation can be strongly reduced for even very high WIMP mass. The statistical error of the reconstructed mass has also been reduced. In a background-free evironment, a WIMP mass of ~ 50 GeV could in principle be determined with an error of ~ 35% with only 2 times 50 events.

Weakly Interacting Massive Particles (WIMPs) are one of the leading candidates for Dark Matter. We develop a model-independent method for determining the mass $m_\chi$ of the WIMP by using data (i.e., measured recoil energies) of direct detection experiments. Our method is independent of the as yet unknown WIMP density near the Earth, of the form of the WIMP velocity distribution, as well as of the WIMP-nucleus cross section. However, it requires positive signals from at least two detectors with different target nuclei. In a background-free environment, $m_\chi \sim 50$ GeV could in principle be determined with an error of $\sim 35%$ with only $2 \times 50$ events; in practice upper and lower limits on the recoil energy of signal events, imposed to reduce backgrounds, can increase the error. The method also loses precision if $m_\chi$ significantly exceeds the mass of the heaviest target nucleus used.

We investigate the relic density n_\chi of non-relativistic long-lived or stable particles \chi in non-standard cosmological scenarios. We calculate the relic abundance starting from arbitrary initial temperatures of the radiation-dominated epoch, and derive the lower bound on the initial temperature T_0 \geq m_\chi/23, assuming that thermally produced \chi particles account for the dark matter energy density in the universe; this bound holds for all \chi annihilation cross sections. We also investigate cosmological scenarios with modified expansion rate. Even in this case an approximate formula similar to the standard one is capable of predicting the final relic abundance correctly. Choosing the \chi annihilation cross section such that the observed cold dark matter abundance is reproduced in standard cosmology, we constrain possible modifications of the expansion rate at T \sim m_\chi/20, well before Big Bang nucleosynthesis.

Weakly interacting massive particles (WIMPs) are one of the leading candidates for Dark Matter. So far we can use direct Dark Matter detection to estimate the mass of halo WIMPs only by fitting predicted recoil spectra to future experimental data. Here we develop a model-independent method for determining the WIMP mass by using experimental data directly. This method is independent of the as yet unknown WIMP density near the Earth as well as of the WIMP-nuclear cross section and can be used to extract information about WIMP mass with O(50) events.

Weakly interacting massive particles (WIMPs) are one of the leading candidates for dark matter. Currently, the most promising method to detect many different WIMP candidates is the direct detection of the recoil energy deposited in a low-background laboratory detector due to elastic WIMP-nucleus scattering. So far the usual procedure has been to predict the event rate of direct detection of WIMPs based on some model(s) of the galactic halo. The aim of our work is to invert this process. That is, we study what future direct detection experiment can teach us about the WIMP halo. As the first step we consider a time-averaged recoil spectrum, assuming that no directional information exists. We develop a method to construct the (time-averaged) one-dimensional velocity distribution function from this spectrum. Moments of this function, such as the mean velocity and velocity dispersion of WIMPs, can also be obtained directly from the recoil spectrum. The only input needed in addition to a measured recoil spectrum is the mass of the WIMP; no information about the scattering cross section or WIMP density is required.

We investigate the relic density n_\chi of non-relativistic long-lived or stable particles \chi in cosmological scenarios in which the temperature T is too low for \chi to achieve full chemical equilibrium. The case with a heavier particle decaying into \chi is also investigated. We derive approximate solutions for n_\chi(T) which accurately reproduce numerical results when full thermal equilibrium is not achieved. If full equilibrium is reached, our ansatz no longer reproduces the correct temperature dependence of the \chi number density. However, it does give the correct final relic density, to an accuracy of about 3% or better, for all cross sections and initial temperatures.

We provide an up-to-date analysis of the parameter space of the minimal supergravity model (mSUGRA). Novel features include the new central value of the top quark mass, an improved calculation of the masses of the supersymmetric particles and the neutral Higgs bosons, constraints from b to s \ell^+ \ell^- decays, and a careful treatment of the most important experimental and theoretical uncertainties. In addition to the by now traditional plots of the allowed region in the (m_0, m_{1/2}) plane, we show allowed regions in the planes spanned by pairs of {\em physical} sparticle or Higgs boson masses. Moreover, we search for the minimal allowed masses of new particles for various sets of constraints. We find that in many cases the direct experimental limits from collider and Dark Matter searches can be saturated even in this minimal model, and even after including the by now quite restrictive constraint on the Dark Matter relic density.

I summarize some recent work on supersymmetric neutralinos as candidates for cold Dark Matter in the Universe. This includes a new scan of mSUGRA parameter space, with special emphasis on neutralinos annihilating predominantly through exchange of the light CP--even Higgs boson, and on bounds on sparticle masses. Next, prospects of testing models with TeV higgsino--like Dark Matter at colliders are discussed. Finally, I briefly comment on extensions of the mSUGRA model, and on scenarios with non--standard cosmology.

The requirement that the lightest neutralino $\tilde\chi_1^0$ has the right thermal relic density to explain all Dark Matter in the universe strongly constrains the parameter space of supersymmetric models in general, and of the mSUGRA model in particular. Recently improved calculations of the mass of the light CP-even Higgs boson $h$ present in this model, and the increased central value of the mass of the top quark, have re--opened the possibility that $2 \mlsp \lsim m_h$. In this ``$h-$pole region'' the LSP annihilation cross section is enhanced by near-resonant $h$ exchange in the $s-$channel, reducing the relic density to acceptable values. We delineate the corresponding region of mSUGRA parameter space, and explore its phenomenology. In particular, we find strong upper bounds on the masses of the gluino, lighter chargino and LSP.

We investigate the possibility that the entropy producing decay of a right-handed sneutrino condensate can simultaneously be the source of the baryon asymmetry, of gravitino dark matter, and of cosmological density perturbations. For generic values of soft supersymmetry breaking terms in the visible sector of 1-10 TeV, condensate decay can yield the dark matter abundance for gravitinos in the mass range 1 MeV to 1 TeV, provided that the resulting reheat temperature is below $10^6$ GeV. The abundance of thermally produced gravitinos before and after sneutrino decay is then negligible. We consider different leptogenesis mechanisms to generate a sufficient asymmetry, and find that low-scale soft leptogenesis works most naturally at such temperatures. The condensate can easily generate sufficient density perturbations if its initial amplitude is $\sim {\cal O}(M_{\rm GUT})$, for a Hubble expansion rate during inflation $> 10^9$ GeV. Right-handed sneutrinos may therefore at the same time provide a source for baryogenesis, dark matter and the seed of structure formation.

The origin of Ultra--High Energy ($E \gsim 10^{20}$ eV) Cosmic Rays (UHECR) remains mysterious. I discuss ``top--down'' models, where UHECR originate from the decay of very massive, long--lived particles. I summarize the calculation of the spectrum of decay products, discuss possible problems with this scenario, and describe ways to test it by searching for very energetic neutrinos and neutralinos.

Decays of superheavy X particles with mass M_X ~ 10^12 - 10^16 GeV have been proposed as origin of the observed ultra high energy cosmic rays (UHECR). We describe in detail the physics involved in the different steps of the decay of such a particle. In particular, we give for the first time the complete set of splitting functions needed to model a parton shower in the minimal supersymmetric extension of the Standard Model (MSSM). We present our results in the form of fragmentation functions of any (s)particle of the MSSM into any final stable particle (proton, photon, electron, three types of neutrino, lightest superparticle LSP) at a virtuality Q = M_X, over a scaled energy range x = 2E/M_X in [10^{-13}, 1]. Extending the coverage to such small fractional energies is necessary since the energy region around 10^18 eV and below could be of considerable interest in testing this kind of model for generating UHECR. We explicitly demonstrate that our treatment conserves energy, and discuss the dependence of the final result on SUSY parameters. We also show that our results are essentially independent of the necessary extrapolation of the input fragmentation functions, which are known only for x >= 0.1, towards small x. Finally, we added a new treatment of the color coherence effects at very small x, using the analytic ``MLLA'' solution. Our computer code will soon be made available.

In this talk I discuss the interplay between collider physics and four topics of astro-particle physics: neutrino oscillations, electroweak baryogenesis, LSP Dark Matter, and ultra-high energy cosmic rays (UHECR). Some astrophysical scenarios can (only) be tested decisively at colliders. In other cases input from collider experiments is required to sharpen predictions for future astro--particle physics experiments, e.g. for the LSP detection rate or the UHECR spectrum in top--down models.

We discuss thermalization through perturbative inflaton decay at the end of inflation. We find that a thermal plasma should form well before all inflatons have decayed, unless all gauge symmetries are badly broken during that epoch. However, before they thermalize, the very energetic inflaton decay products can contribute to the production of massive stable particles, either through collisions with the thermal plasma, or through collisions with each other. If such reactions exist, the same massive particles can also be produced directly in inflaton decay, once higher--order processes are included. We show that these new, non--thermal production mechanisms often significantly strengthen constraints on the parameters of models containing massive stable particles; for example, stable charged particles with mass below the inflaton mass seem to be essentially excluded.

We point out that inflaton decays can be a copious source of stable or long--lived particles $\chi$ with mass exceeding the reheat temperature $T_R$. Once higher order processes are included, this statement is true for any $\chi$ particle with renormalizable (gauge or Yukawa) interactions. This contribution to the $\chi$ density often exceeds the contribution from thermal $\chi$ production, leading to significantly stronger constraints on model parameters than those resulting from thermal $\chi$ production alone. For example, we all but exclude models containing stable charged particles with mass less than half the mass of the inflaton.

Top-down models of cosmic rays produce more neutrinos than photons and more photons than protons. In these models, we reevaluate the fluxes of neutrinos associated with the highest energy cosmic rays in light of mounting evidence that they are protons and not gamma rays. While proton dominance at EeV energies can possibly be achieved by efficient absorption of the dominant high-energy photon flux on universal and galactic photon and magnetic background fields, we show that the associated neutrino flux is inevitably increased to a level where it should be within reach of operating experiments such as AMANDA II, RICE and AGASA. In future neutrino telescopes, tens to a hundred, rather than a few neutrinos per kilometer squared per year, may be detected above 1 PeV.

We present a new and complete numerical analysis of the decay of super-heavy X particles, assumed to be the origin of cosmic rays with energy beyond the GZK cut-off. The decay of X initiates a ``parton shower'', where we include all degrees of freedom contained in the Minimal Supersymmetric Standard Model (MSSM). Since at energies near $M_X$ all gauge couplings are of similar magnitude, we include all of them, as well as third generation Yukawa couplings. Technically the shower development is described through the DGLAP evolution of the relevant fragmentation functions (FFs). We also carefully treat the decay of the superparticles as well as heavy SM particles created in the shower. Nonperturbative physics is parameterized through the input values of the FFs, which we take from the literature. The final result is the the complete spectrum of all stable particles at the very end of the shower : protons, electrons, neutrinos, photons and neutralinos, for an energy range from $10^{-7} M_X$ to $M_X$. In particular, the flux of high-energy neutralinos is sizable ; it might serve as ``smoking gun'' signature for this kind of scenario.

We calculate the one-loop induced couplings of two gaugino-like neutralinos to the $Z$ and Higgs bosons in the Minimal Supersymmetric Standard Model. These couplings, which vanish at the tree level, can be generated through loops involving fermions and sfermions. We show that, while the neutralino contribution to the invisible $Z$ boson decay width remains small, the loop induced couplings to the lightest Higgs boson might be sufficiently large to yield a rate of invisible decays of this Higgs boson that should be detectable at future $e^+e^-$ colliders. We also study the implications of these couplings for direct searches of Dark Matter and show that they can modify appreciably the neutralino-nucleon elastic cross section for some parameter range.

We examine the viability of the Affleck-Dine mechanism for baryogenesis under radiatively induced running of soft breaking $({\rm mass})^2$ of the flat directions stemming from non-zero energy density of the inflaton during inflation. A major difference from analogous phenomenological studies is that the horizon radius provides a natural infrared cut-off to the quantum corrections in this case. We identify different scenarios which may arise and point out that the $H_{u}L$ flat direction remains the most promising flat direction, since it is largely independent of uncertainties about high scale physics and details of the inflationary model.

We perform a complete analysis of the supersymmetric particle spectrum in the Minimal Supergravity (mSUGRA) model where the soft SUSY breaking scalar masses, gaugino masses and trilinear couplings are unified at the GUT scale, so that the electroweak symmetry is broken radiatively. We show that the present constraints on the Higgs boson and superparticle masses from collider searches and precision measurements still allow for large regions of the mSUGRA parameter space where charginos, neutralinos, sleptons and top squarks as well as the heavier Higgs particles, are light enough to be produced at the next generation of $\ee$ linear colliders with center of mass energy around $\sqrt{s} \sim 800$ GeV, with sizeable cross sections. An important part of this parameter space remains even when we require that the density of the lightest neutralinos left over from the Big Bang, which we calculate using standard assumptions, falls in the range favored by current determinations of the Dark Matter density in the Universe. Already at a c.m. energy of 500 GeV, SUSY particles can be accessible in some parameter range, and if the energy is increased to $\sqrt{s} \simeq 1.2$ TeV, the $e^+e^-$ collider will have a reach for high precision studies of SUSY particles in a range that is comparable to the discovery range of the LHC.

We compare predictions for the spin-independent contribution to the neutralino-proton scattering cross section $\sig$ and for the anomalous magnetic dipole moment of the muon, $a_\mu = (g_\mu-2)/2$, in models with gravity-mediated supersymmetry breaking. We nearly always find a positive correlation between these two measurables, i.e. scenarios with larger $a_\mu$ also tend to have larger $\sig$, but the detailed prediction differs greatly between models. In particular, we find that for the popular mSUGRA scenario with universal soft breaking masses at the scale of Grand Unification, measurements of $a_\mu$ currently seem more promising. On the other hand, if scalar soft breaking masses at the GUT scale receive sizable contributions from SO(10) D-terms, one often finds scenarios with large $\sig$ but $a_\mu$ below the currently foreseen sensitivity. A string-inspired model with non-universal scalar spectrum at the GUT scale falls between these two cases.

We show that LHC experiments might well be able to determine all the parameters required for a prediction of the present density of thermal LSP relics from the Big Bang era. If the LSP is an almost pure bino we usually only need to determine its mass and the mass of the SU(2) singlet sleptons. This information can be obtained by reconstructing the cascade $\tilde{q}_L \to \tilde{\chi}_2^0 q \to \tilde{\ell}_R \ell q \to \tilde{\chi}_1^0 \ell^+ \ell^- q$. The only requirement is that $m_{\tilde{\ell}_R} < m_{\tilde{\chi}_2^0}$, which is true for most of the cosmologically interesting parameter space. If the LSP has a significant higgsino component, its predicted thermal relic density is smaller than for an equal--mass bino. We show that in this case squark decays also produce significant numbers of $\tilde{\chi}_4^0$ and $\tilde{\chi}_2^\pm$. Reconstructing the corresponding decay cascades then allows to determine the higgsino component of the LSP.

It has been suggested that a sterile neutrino \nu_s which mixes with standard neutrinos can form nonthermal ``cool'' Dark Matter if its mass and mixing angle fall in the ranges 0.1 keV \lsim m_s \lsim 10 keV and 10^{-10} \lsim \sin^2 \theta \lsim 10^{-4}, respectively. We point out that the required mixing makes these heavy neutrinos unstable. The dominant decay mode is into three light neutrinos, but the most stringent constraint comes from the non-observation of radiative decays into a single light neutrino and a photon. Moreover, we point out that the density of thermal relics of such \nu_{s}'s would be too high, unless the reheating temperature after inflation was below \sim 10 GeV.

We calculate the dominant loop corrections from both standard and supersymmetric QCD to the effective coupling of neutralinos to nucleons. The potentially largest corrections come from gluino-squark loop contributions to the Higgs boson couplings to quarks; these corrections also affect the leading spin-independent squark exchange contribution. Ordinary QCD corrections to the effective coupling of CP-even Higgs bosons to two gluons are also sizable. For large $\tanb$ values, i.e. in the region of parameter space probed by current and near-future direct Dark-Matter search experiments, the total corrections can exceed a factor of three.

I point out that the sterile neutrinos suggested as candidates of "cool" Dark Matter will decay through their mixing with light neutrinos. This leads to an upper bound of about 200 keV on the mass of the sterile neutrinos, but might facilitate their detection.

We present a complete calculation of the cross section for neutralino annihilation into the two-gluon final state. This channel can be quite important for the phenomenology of neutralino annihilation due to the well-known helicity suppression of neutralino annihilation into light quarks and leptons. In addition, we calculate the cross section for annihilation of neutralinos into a gluon and quark-antiquark pair, and discuss QCD corrections to the tree-level cross sections for neutralino annihilation into quarks. If the neutralino is lighter than the top quark, the effect of these results on high-energy neutrino signals from neutralino annihilation in the Sun and in the Earth can be significant, especially if the neutralino is primarily gaugino. On the other hand, our results should have little effect on calculations of the cosmological abundance of neutralinos. We also briefly discuss implications for cosmic-ray antiprotons from neutralino annihilation in the galactic halo.