The question of the nature of the dark matter in the Universe remains one of the most outstanding unsolved problems in basic science. One of the best motivated particle physics candidates is the lightest supersymmetric particle, assumed to be the lightest neutralino -a linear combination of the supersymmetric partners of the photon, the Z boson and neutral scalar Higgs particles. Here we describe DarkSUSY, a publicly-available advanced numerical package for neutralino dark matter calculations. In DarkSUSY one can compute the neutralino density in the Universe today using precision methods which include resonances, pair production thresholds and coannihilations. Masses and mixings of supersymmetric particles can be computed within DarkSUSY or with the help of external programs such as FeynHiggs, ISASUGRA and SUSPECT. Accelerator bounds can be checked to identify viable dark matter candidates. DarkSUSY also computes a large variety of astrophysical signals from neutralino dark matter, such as direct detection in low-background counting experiments and indirect detection through antiprotons, antideuterons, gamma-rays and positrons from the Galactic halo or high-energy neutrinos from the center of the Earth or of the Sun. Here we describe the physics behind the package. A detailed manual will be provided with the computer package.

Supersymmetric neutralinos are one of the most promising candidates for the dark matter in the Universe. If they exist, they should make up some fraction of the Milky Way halo. We investigate the fluxes of positrons expected at the Earth from neutralino annihilation in the halo. Positron propagation is treated in a diffusion model including energy loss. The positron source function includes contributions from both continuum and monochromatic positrons. We find that, for a "canonical" halo model and propagation parameters, the fluxes are generally too low to be visible. Given the large uncertainties in both propagation and halo structure, it is however possible to obtain observable fluxes. We also investigate the shapes of the positron spectra, including fits to a feature indicated by the results of the Heat experiment. 95.35.+d, 14.80.Ly, 98.70.Sa

If the cosmic dark matter consists of weakly-interacting massive particles, these particles should be produced in reactions at the next generation of high-energy accelerators. Measurements at these accelerators can then be used to determine the microscopic properties of the dark matter. From this, we can predict the cosmic density, the annihilation cross sections, and the cross sections relevant to direct detection. In this paper, we present studies in supersymmetry models with neutralino dark matter that give quantitative estimates of the accuracy that can be expected. We show that these are well matched to the requirements of anticipated astrophysical observations of dark matter. The capabilities of the proposed International Linear Collider (ILC) are expected to play a particularly important role in this study.