Ten weeks of daily imaging of the giant elliptical galaxy M87 with the Hubble Space Telescope (HST) has yielded 41 nova light curves of unprecedented quality for extragalactic cataclysmic variables. We have recently used these light curves to demonstrate that the observational scatter in the so-called Maximum-Magnitude Rate of Decline (MMRD) relation for classical novae is so large as to render the nova-MMRD useless as a standard candle. Here we demonstrate that a modified Buscombe - de Vaucouleurs hypothesis, namely that novae with decline times t2 > 10 days converge to nearly the same absolute magnitude about two weeks after maximum light in a giant elliptical galaxy, is supported by our M87 nova data. For 13 novae with daily-sampled light curves, well determined times of maximum light in both the F606W and F814W filters, and decline times $t2 > 10 days we find that M87 novae display M(606W,15) = -6.37 +/- 0.46 and M(814W,15) = -6.11 +/- 0.43. If very fast novae with decline times t2 < 10 days are excluded, the distances to novae in elliptical galaxies with stellar binary populations similar to those of M87 should be determinable with 1 sigma accuracies of +/-20% with the above calibrations.

The extensive grid of numerical simulations of nova eruptions of Yaron et al.(2005) first predicted that some classical novae might deviate significantly from the Maximum Magnitude - Rate of Decline (MMRD) relation, which purports to characterise novae as standard candles. Kasliwal et al. (2011) have announced the observational detection of an apparently new class of faint, fast classical novae in the Andromeda galaxy. These objects deviate strongly from the MMRD relationship, exactly as predicted by Yaron et al. (2005). Shara et al. (2016) recently reported the first detections of faint, fast novae in M87. These previously overlooked objects are as common in the giant elliptical galaxy M87 as they are in the giant spiral M31; they comprise about 40% of all classical nova eruptions and greatly increase the observational scatter in the MMRD relation. We use the extensive grid of nova simulations of Yaron et al. (2005) to identify the underlying causes of the existence of faint, fast novae. These are systems which have accreted, and can thus eject, only very low mass envelopes, of order 10^-7 - 10^-8 Msun, on massive white dwarfs. Such binaries include, but are not limited to, the recurrent novae. These same models predict the existence of ultrafast novae which display decline times t2 as short as five hours. We outline a strategy for their future detection.

The Hubble Space Telescope has imaged the central part of M87 over a 10 week span, leading to the discovery of 32 classical novae and nine fainter, likely very slow and/or symbiotic novae. In this first in a series of papers we present the M87 nova finder charts, and the light and color curves of the novae. We demonstrate that the rise and decline times, and the colors of M87 novae are uncorrelated with each other and with position in the galaxy. The spatial distribution of the M87 novae follows the light of the galaxy, suggesting that novae accreted by M87 during cannibalistic episodes are well-mixed. Conservatively using only the 32 brightest classical novae we derive a nova rate for M87: $363_{-45}^{+33}$ novae/yr. We also derive the luminosity-specific classical nova rate for this galaxy, which is $7.88_{-2.6}^{+2.3} /yr/ 10^{10}L_\odot,_{K}$. Both rates are 3-4 times higher higher than those reported for M87 in the past, and similarly higher than those reported for all other galaxies. We suggest that most previous ground-based surveys for novae in external galaxies, including M87, miss most faint, fast novae, and almost all slow novae near the centers of galaxies.

Gravitational lenses on galaxy scales are plausibly modelled as having ellipsoidal symmetry and a universal dark matter density profile, with a Sersic profile to describe the distribution of baryonic matter. Predicting all lensing effects requires knowledge of the total lens potential: in this work we give analytic forms for that of the above hybrid model. Emphasising that complex lens potentials can be constructed from simpler components in linear combination, we provide a recipe for attaining elliptical symmetry in either projected mass or lens potential. We also provide analytic formulae for the lens potentials of Sersic profiles for integer and half-integer index. We then present formulae describing the gravitational lensing effects due to smoothly-truncated universal density profiles in cold dark matter model. For our isolated haloes the density profile falls off as radius to the minus fifth or seventh power beyond the tidal radius, functional forms that allow all orders of lens potential derivatives to be calculated analytically, while ensuring a non-divergent total mass. We show how the observables predicted by this profile differ from that of the original infinite-mass NFW profile. Expressions for the gravitational flexion are highlighted. We show how decreasing the tidal radius allows stripped haloes to be modelled, providing a framework for a fuller investigation of dark matter substructure in galaxies and clusters. Finally we remark on the need for finite mass halo profiles when doing cosmological ray-tracing simulations, and the need for readily-calculable higher order derivatives of the lens potential when studying catastrophes in strong lenses.

This is a report on the findings of the dark matter science working group for the white paper on the status and future of TeV gamma-ray astronomy. The white paper was commissioned by the American Physical Society, and the full white paper can be found on astro-ph (arXiv:0810.0444). This detailed section discusses the prospects for dark matter detection with future gamma-ray experiments, and the complementarity of gamma-ray measurements with other indirect, direct or accelerator-based searches. We conclude that any comprehensive search for dark matter should include gamma-ray observations, both to identify the dark matter particle (through the charac- teristics of the gamma-ray spectrum) and to measure the distribution of dark matter in galactic halos.

We present an analysis of a 50 orbit HST ACS observation of the M87 globular cluster system. We use the extraordinary depth of this dataset to test whether the colors and magnitudes show evidence for a mass-metallicity relation in globular cluster populations. We find only a weak or absent relation between the colors and magnitudes of the metal poor subpopulation of globular clusters. The weakness or absence of a color-magnitude relation is established over a wide range in luminosity from $M_V=-11$ to $M_V=-6$, encompassing most of the M87 globular clusters. The constancy of the colors of the metal-poor subpopulation seen in our 50 orbit observation is in contrast to suggestions from single orbit ACS data that the metal-poor globular clusters in M87 and several other galaxies show a "blue tilt." The formal best fit for the mass-metallicity relation for the metal-poor subpopulation in our much deeper data is $Z\propto M^{0.08\pm0.05}$. Our analysis of these data also shows a possible small "red tilt" in the metal-rich globular cluster subpopulation. While either of these small tilts may be real, they may also illustrate the limit to which mass-metallicity relations can be determined, even in such extraordinarily deep data. We specifically test for a wide range of systematic effects and find that while small tilts cannot be confirmed or rejected, the data place a strong upper limit to any tilt of $|0.20|\pm0.05$. This upper limit is much smaller than some earlier claims from single orbit data, and strongly limits self-enrichment within globular clusters. This mass-metallicity relation for globular clusters is also shallower than the relation for galaxies, suggesting that the formation mechanisms for these two types of objects are different.

We investigate the sensitivity of the Gamma-ray Large Area Space Telescope (GLAST) to indirectly detect weakly interacting massive particles (WIMPs) through the $\gamma$-ray signal that their pair annihilation produces. WIMPs are among the favorite candidates to explain the compelling evidence that about 80% of the mass in the Universe is non-baryonic dark matter (DM). They are serendipitously motivated by various extensions of the standard model of particle physics such as Supersymmetry and Universal Extra Dimensions (UED). With its unprecedented sensitivity and its very large energy range (20 MeV to more than 300 GeV) the main instrument on board the GLAST satellite, the Large Area Telescope (LAT), will open a new window of discovery. As our estimates show, the LAT will be able to detect an indirect DM signature for a large class of WIMP models given a cuspy profile for the DM distribution. Using the current state of the art Monte Carlo and event reconstruction software developed within the LAT collaboration, we present preliminary sensitivity studies for several possible sources inside and outside the Galaxy. We also discuss the potential of the LAT to detect UED via the electron/positron channel. Diffuse background modeling and other background issues that will be important in setting limits or seeing a signal are presented.

In this review, we discuss the role of the various experimental programs taking part in the broader effort to identify the particle nature of dark matter. In particular, we focus on electroweak scale dark matter particles and discuss a wide range of search strategies being carried out and developed to detect them. These efforts include direct detection experiments, which attempt to observe the elastic scattering of dark matter particles with nuclei, indirect detection experiments, which search for photons, antimatter and neutrinos produced as a result of dark matter annihilations, and collider searches for new TeV-scale physics. Each of these techniques could potentially provide a different and complementary set of information related to the mass, interactions and distribution of dark matter. Ultimately, it is hoped that these many different tools will be used together to conclusively identify the particle or particles that constitute the dark matter of our universe.

This Conference is the fifth of a series of Workshops on High Energy Gamma- ray Experiments, following the Conferences held in Perugia 2003, Bari 2004, Cividale del Friuli 2005, Elba Island 2006. This year the focus was on the use of gamma-ray to study the Dark Matter component of the Universe, the origin and propagation of Cosmic Rays, Extra Large Spatial Dimensions and Tests of Lorentz Invariance.

We present here the prospects for the GLAST Large Area Telescope (LAT) detection of the signature of the lightest Kaluza-Klein particle (LKP). It decays by direct annihilation into electron-positron pairs that may be detectable in the high energy electron flux. We discuss the LAT capability for detecting the high energy (20 GeV - ~1 TeV) cosmic ray electron flux and we analyze the LAT sensitivity to detect LKP-produced electrons for various particle masses. We include an analysis of the diffusive propagation of the electrons in the galaxy.

This paper has been withdrawn.

This paper has been withdrawn.

A study of the current status of WIMP dark matter searches has been made in the context of scientific and technical planning for a Deep Underground Science and Engineering Laboratory (DUSEL) in the U.S. The table of contents follows: 1. Overview 2. WIMP Dark Matter: Cosmology, Astrophysics, and Particle Physics 3. Direct Detection of WIMPs 4. Indirect Detection of WIMPs 5. Dark Matter Candidates and New Physics in the Laboratory 6. Synergies with Other Sub-Fields 7. Direct Detection Experiments: Status and Future Prospects 8. Infrastructure 9. International Context 10. Summary and Outlook 11. Acknowledgments

The nature of the cosmic dark matter is unknown. The most compelling hypothesis is that dark matter consists of weakly interacting massive particles (WIMPs) in the 100 GeV mass range. Such particles would annihilate in the galactic halo, producing high-energy gamma rays which might be detectable in gamma ray telescopes such as the GLAST satellite. We investigate the ability of GLAST to distinguish between WIMP annihilation sources and astrophysical sources. Focusing on the galactic satellite halos predicted by the cold dark matter model, we find that the WIMP gamma-ray spectrum is nearly unique; separation of the brightest WIMP sources from known source classes can be done in a convincing way by including spectral and spatial information. Candidate WIMP sources can be further studied with Imaging Atmospheric Cerenkov Telescopes. Finally, Large Hadron Collider data might have a crucial impact on the study of galactic dark matter.

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.

We present the luminosity function to very faint magnitudes for the globular clusters in M87, based on a 30 orbit \textit{Hubble Space Telescope (HST)} WFPC2 imaging program. The very deep images and corresponding improved false source rejection allow us to probe the mass function further beyond the turnover than has been done before. We compare our luminosity function to those that have been observed in the past, and confirm the similarity of the turnover luminosity between M87 and the Milky Way. We also find with high statistical significance that the M87 luminosity function is broader than that of the Milky Way. We discuss how determining the mass function of the cluster system to low masses can constrain theoretical models of the dynamical evolution of globular cluster systems. Our mass function is consistent with the dependence of mass loss on the initial cluster mass given by classical evaporation, and somewhat inconsistent with newer proposals that have a shallower mass dependence. In addition, the rate of mass loss is consistent with standard evaporation models, and not with the much higher rates proposed by some recent studies of very young cluster systems. We also find that the mass-size relation has very little slope, indicating that there is almost no increase in the size of a cluster with increasing mass.

We present results of a microlensing survey toward the Andromeda Galaxy (M31) carried out during four observing seasons at the Isaac Newton Telescope (INT). This survey is part of the larger microlensing survey toward M31 performed by the Microlensing Exploration of the Galaxy and Andromeda (MEGA) collaboration. Using a fully automated search algorithm, we indentify 14 candidate microlensing events, three of which are reported here for the first time. Observations obtained at the Mayall telescope are combined with the INT data to produce composite lightcurves for these candidates. The results from the survey are compared with theoretical predictions for the number and distribution of events. These predictions are based on a Monte Carlo calculation of the detection efficiency and disk-bulge-halo models for M31. The models provide the full phase-space distribution functions for the lens and source populations and are motivated by dynamical and observational considerations. They include differential extinction and span a wide range of parameter space characterised primarily by the mass-to-light ratios for the disk and bulge. For most models, the observed event rate is consistent with the rate predicted for self-lensing -- a MACHO halo fraction of 30% or higher can be ruled at the 95% confidence level. The event distribution does show a large near-far asymmetry hinting at a halo contribution to the microlensing signal. Two candidate events are located at particularly large projected radii on the far side of the disk. These events are difficult to explain by self lensing and only somewhat easier to explain by MACHO lensing. A possibility is that one of these is due to a lens in a giant stellar stream.

We investigate $HST$/ACS and WFPC2 images at the positions of five candidate microlensing events from a large survey of variability in M31 (MEGA). Three closely match unresolved sources, and two produce only flux upper limits. All are confined to regions of the color-magnitude diagram where stellar variability is unlikely to be easily confused with microlensing. Red variable stars cannot explain these events (although background supernova are possible for two). If these lenses arise in M31's halo, they are due to masses $0.15 < m / M_\odot < 0.49$ (95% certainty, for a $\delta$-function mass distribution), brown dwarfs for disk lenses, and stellar masses for bulge lenses.

It is now widely accepted that most of mass--energy in the universe is unobserved except by its gravitational effects. Baryons make only about 4% of the total, with "dark matter" making up about 23% and the "dark energy" responsible for the accelerated expansion of the universe making up the remainder. We focus on the dark matter, which is the primary constituent of galaxies. We outline the observed properties of this material, enumerating some candidates covering 90 orders of magnitude in mass. Finally, we argue that the weak scale (100 GeV) is relevant to new physics, including the dark matter problem.

Kaluza-Klein dark matter particles can annihilate efficiently into electron-positron pairs, providing a discrete feature (a sharp edge) in the cosmic $e^+ e^-$ spectrum at an energy equal to the particle's mass (typically several hundred GeV to one TeV). Although this feature is probably beyond the reach of satellite or balloon-based cosmic ray experiments (those that distinguish the charge and mass of the primary particle), gamma ray telescopes may provide an alternative detection method. Designed to observe very high-energy gamma-rays, ACTs also observe the diffuse flux of electron-induced electromagnetic showers. The GLAST satellite, designed for gamma ray astronomy, will also observe any high energy showers (several hundred GeV and above) in its calorimeter. We show that high-significance detections of an electron-positron feature from Kaluza-Klein dark matter annihilations are possible with GLAST, and also with ACTs such as HESS, VERITAS or MAGIC.

The Andromeda Galaxy (M31) is the closest large galaxy to the Milky Way, thus it is an important laboratory for studying massive dark objects in galactic halos (MACHOs) by gravitational microlensing. Such studies strongly complement the studies of the Milky Way halo using the the Large and Small Magellanic Clouds. We consider the possibilities for microlensing surveys of M31 using the next generation of wide field imaging telescopes with fields of view in the square degree range. We consider proposals for such imagers both on the ground and in space. For concreteness, we specialize to the SNAP proposal for a space telescope and the LSST proposal for a ground based telescope. We find that a modest space-based survey of 50 visits of one hour each is considerably better than current ground based surveys covering 5 years. Crucially, systematic effects can be considerably better controlled with a space telescope because of both the infrared sensitivity and the angular resolution. To be competitive, 8 meter class wide-field ground based imagers must take exposures of several hundred seconds with several day cadence.

Our fundamental lack of understanding of the acceleration of the Universe suggests that we consider a kinematic description. The simplest formalism involves the third derivative of the scale factor through a jerk parameter. A new approach is presented for describing the results of astronomical observations in terms of the contemporary jerk parameter and this is related to the equation of state approach. Simple perturbative expansions about $\Lambda$CDM are given.

We explore the full parameter space of Minimal Supergravity (mSUGRA), allowing all four continuous parameters (the scalar mass m_0, the gaugino mass m_1/2, the trilinear coupling A_0, and the ratio of Higgs vacuum expectation values tan beta) to vary freely. We apply current accelerator constraints on sparticle and Higgs masses, and on the b -> s gamma branching ratio, and discuss the impact of the constraints on g_mu-2. To study dark matter, we apply the WMAP constraint on the cold dark matter density. We develop Markov Chain Monte Carlo (MCMC) techniques to explore the parameter regions consistent with WMAP, finding them to be considerably superior to previously used methods for exploring supersymmetric parameter spaces. Finally, we study the reach of current and future direct detection experiments in light of the WMAP constraint.

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.

Annihilating dark matter particles produce roughly as much power in electrons and positrons as in gamma ray photons. The charged particles lose essentially all of their energy to inverse Compton and synchrotron processes in the galactic environment. We discuss the diffuse signature of dark matter annihilations in satellites of the Milky Way (which may be optically dark with few or no stars), providing a tail of emission trailing the satellite in its orbit. Inverse Compton processes provide X-rays and gamma rays, and synchrotron emission at radio wavelengths might be seen. We discuss the possibility of detecting these signals with current and future observations, in particular EGRET and GLAST for the gamma rays.

We report results from intensive monitoring of two fields on either side of M31,emphasizing microlensing involving stars and masses in M31. These results stem from the three-year VATT/Columbia survey of variability on 3d to 2m timescales. Observations were conducted intensively from 1997-1999, with baselines 1995-now,at the Vatican Advanced Technology Telescope, MDM 1.3-meter, and Isaac Newton telescopes, totaling about 200 epochs. The two fields cover 560 square arcmin total, along the minor axis on either side of M31. Candidate microlensing events are subject to tests to distinguish microlenses from variable stars. Fourprobable microlensing events total, compared to carefully computed event rate and efficiency models, indicate a marginally significant microlensing activity above expectations for the stellar lenses alone in M31 (and the Galaxy) acting as lenses. Maximum likelihood analyses of the event distribution in timescale and across the face of M31 indicates a microlensing dark matter halo fraction consistent with that seen in our Galaxy towards the LMC (Alcock et al. 2000). For a nearly singular isothermal sphere model, we find a microlensing halo mass fraction f=0.29 +0.30/-0.13 of total dark matter, and a poorly constrained lensing component mass 0.02 to 1.5 solar mass, 1 sigma limits). This study is the prototype for a larger one approaching completion; between them there is significant evidence for asymmetry in the distribution of microlensing events across the face of M31, and a large population of halo microlensing objects. (Abridged)

Only one certain classical nova eruption has ever been detected inside a globular cluster - nova 1860 A.D. (T Sco) in M80. During a survey of M87 we have detected an erupting star coincident (to within 0.08 pixels) with a globular cluster of that giant elliptical galaxy. We are able to discount variables in the foreground or background of M87. The light curve and color of the erupting star match those expected for a nova at the distance of M87. The chance superposition of an M87 field nova on the globular cluster is very unlikely but cannot be completely ruled out.Our detection hints at a globular cluster nova frequency $f \sim .004$ novae/cluster/year, much higher than previous observations have suggested.

The position of the giant elliptical galaxy M87 at the center of the Virgo Cluster means that the inferred column density of dark matter associated with both the cluster halo and the galaxy halo is quite large. This system is thus an important laboratory for studying massive dark objects in elliptical galaxies and galaxy clusters by gravitational microlensing, strongly complementing the studies of spiral galaxy halos performed in the Local Group. We have performed a microlensing survey of M87 with the WFPC2 instrument on the Hubble Space Telescope. Over a period of thirty days, with images taken once daily, we discover seven variable sources. Four are variable stars of some sort, two are consistent with classical novae, and one exhibits an excellent microlensing lightcurve, though with a very blue color implying the somewhat disfavored possibility of a horizontal branch source being lensed. Based on sensitivity calculations from artificial stars and from artificial lightcurves, we estimate the expected microlensing rate. We find that the detection of one event is consistent with a dark halo with a 20% contribution of microlensing objects for both M87 and the Virgo Cluster, similar to the value found from observations in the Local Group. Further work is required to test the hypothesized microlensing component to the cluster.

The microlensing interpretation of the optical afterglow of GRB 000301C seems naively surprising, since a simple estimate of the stellar microlensing rate gives less than one in four hundred for a flat Omega_Lambda=0.7 cosmology, whereas one event was seen in about thirty afterglows. Considering baryonic MACHOs making up half of the baryons in the universe, the microlensing probability per burst can be roughly 5% for a GRB at redshift z=2. We explore two effects that may enhance the probability of observing microlensed gamma-ray burst afterglows: binary lenses and double magnification bias. We find that the consideration of binary lenses can increase the rate only at the ~15% level. On the other hand, because gamma-ray bursts for which afterglow observations exist are typically selected based on fluxes at widely separated wavebands which are not necessarily well correlated (e.g. localization in X-ray, afterglow in optical/infrared), magnification bias can operate at an enhanced level compared to the usual single-bias case. We find that existing estimates of the slope of the luminosity function of gamma-ray bursts, while as yet quite uncertain, point to enhancement factors of more than three above the simple estimates of the microlensing rate. We find that the probability to observe at least one microlensing event in the sample of 27 measured afterglows can be 3-4% for stellar lenses, or as much as 25 Omega_lens for baryonic MACHOs. We note that the probability to observe at least one event over the available sample of afterglows is significant only if a large fraction of the baryons in the universe are condensed in stellar-mass objects. (ABRIDGED)

We present the first M31 candidate microlensing events from the Microlensing Exploration of the Galaxy and Andromeda (MEGA) survey. MEGA uses several telescopes to detect microlensing towards the nearby Andromeda galaxy, M31, in order to establish whether massive compact objects are a significant contribution to the mass budget of the dark halo of M31. The results presented here are based on observations with the Isaac Newton Telescope on La Palma, during the 1999/00 and 2000/01 observing seasons. In this data set, 14 variable sources consistent with microlensing have been detected, 12 of which are new and 2 have been reported previously by the POINT-AGAPE group. A preliminary analysis of the spatial and timescale distributions of the candidate events support their microlensing nature. We compare the spatial distributions of the candidate events and of long-period variable stars, assuming the chances of finding a long-period variable and a microlensing event are comparable. The spatial distribution of our candidate microlensing events is more far/near side asymmetric than expected from the detected long-period variable distribution. The current analysis is preliminary and the asymmetry not highly significant, but the spatial distribution of candidate microlenses is suggestive of the presence of a microlensing halo.

In the nuclear and particle astrophysics session of CIPANP 2003 we heard talks on a number of topics, focused for the most part into four broad areas. Here we outline the discussions of the standard cosmological model, dark matter searches, cosmic rays, and neutrino astrophysics. The robustness of theoretical and experimental programs in all of these areas is very encouraging, and we expect to have many questions answered, and new ones asked, in time for CIPANP 2006.

Weakly interacting massive particles (WIMPs) are among the favored candidates for cold dark matter in the universe. The phenomenology of supersymmetric WIMPs has been quite developed during recent years. However, there are other possibilities which have not been discussed as much. One example is a right-handed massive neutrino, which has recently been proposed in the context of a version of the Zee model for massive neutrinos. This TeV-scale, leptonic WIMP (or LIMP, for short) may at first sight appear to be essentially undetectable. However, we point out that the radiatively induced annihilation rate into leptons and photons is bound to be substantial, and provides a conspicuous gamma-ray signature for annihilations in the galactic halo. This gives a window of opportunity for Air Cerenkov Telescopes with ability to observe the galactic center, such as the HESS and CANGAROO arrays, and also for the GLAST space telescope. In addition, the contribution to the positron cosmic ray flux is in principle detectable, but this would require very strong local density enhancements in the dark matter halo distribution.

Using a new instrument, the HEAT collaboration has confirmed the excess of cosmic ray positrons that they first detected in 1994. We explore the possibility that this excess is due to the annihilation of neutralino dark matter in the galactic halo. We confirm that neutralino annihilation can produce enough positrons to make up the measured excess only if there is an additional enhancement to the signal. We quantify the `boost factor' that is required in the signal for various models in the Minimal Supersymmetric Standard Model parameter space, and find that a boost factor >30 provides good fits to the HEAT data. Such an enhancement in the signal could arise if we live in a clumpy halo.

The question of the nature of 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. We here describe DarkSUSY, an advanced numerical FORTRAN package for supersymmetric dark matter calculations. With DarkSUSY one can: (i) compute masses and compositions of various supersymmetric particles; (ii) compute the relic density of the lightest neutralino, using accurate methods which include the effects of resonances, pair production thresholds and coannihilations; (iii) check accelerator bounds to identify allowed supersymmetric models; and (iv) obtain neutralino detection rates for a variety of detection methods, including direct detection and indirect detection through antiprotons, gamma-rays and positrons from the Galactic halo or neutrinos from the center of the Earth or the Sun.

We provide a set of microlensing event rate maps for M31, the Andromeda Galaxy. Rates for M31 microlensing were calculated on the basis of a four component model of the lens and source populations: disk and bulge sources lensed by bulge, disk, M31 halo and Galactic halo lenses. We confirm the high rate gradient along the minor axis of M31 due to a dark halo of lenses. Furthermore, we compute the timescale distributions of events, for both Einstein times and full-width at half-maximum times. We explore how the rate contours and the timescale distributions can be used to measure the shape and extent of the microlensing halo. With one year of twice--weekly sampling, or three observing seasons, a halo MACHO fraction as small as 5% can be detected with modest ground based telescopes.

The new measurement of the anomalous magnetic moment of the muon by the Brookhaven AGS experiment 821 again shows a discrepancy with the Standard Model value. We investigate the consequences of these new data for neutralino dark matter, updating and extending our previous work [E. A. Baltz and P. Gondolo, Phys. Rev. Lett. 86, 5004 (2001)]. The measurement excludes the Standard Model value at 2.6sigma confidence. Taking the discrepancy as a sign of supersymmetry, we find that the lightest superpartner must be relatively light and it must have a relatively high elastic scattering cross section with nucleons, which brings it almost within reach of proposed direct dark matter searches. The SUSY signal from neutrino telescopes correlates fairly well with the elastic scattering cross section. The rate of cosmic ray antideuterons tends to be large in the allowed models, but the constraint has little effect on the rate of gamma ray lines. We stress that being more conservative may eliminate the discrepancy, but it does not eliminate the possibility of high astrophysical detection rates.

XENON is a novel liquid xenon experiment concept for a sensitive dark matter search using a 1-tonne active target, distributed in an array of ten independent time projection chambers. The design relies on the simultaneous detection of ionization and scintillation signals in liquid xenon, with the goal of extracting as much information as possible on an event-by-event basis, while maintaining most of the target active. XENON is expected to have effective and redundant background identification and discrimination power, higher than 99.5%, and to achieve a very low threshold, on the order of 4 keV visible recoil energy. Based on this expectation and the 1-tonne mass of active xenon, we project a sensitivity of 0.0001 events/kg/day, after 3 yr operation in an appropriate underground location. The XENON experiment has been recently proposed to the National Science Foundation (NSF) for an initial development phase leading to the development of the 100 kg unit module.

Using a new instrument, the HEAT collaboration has confirmed the excess of cosmic ray positrons that they first detected in 1994. We explore the possibility that this excess is due to the annihilation of neutralino dark matter in the galactic halo. We confirm that neutralino annihilation can produce enough positrons to make up the measured excess only if there is an additional enhancement to the signal. We quantify the `boost factor' that is required in the signal for various models in the Minimal Supersymmetric Standard Model parameter space, and study the dependence on various parameters. We find models with a boost factor greater than 30. Such an enhancement in the signal could arise if we live in a clumpy halo. We discuss what part of supersymmetric parameter space is favored (in that it gives the largest positron signal), and the consequences for other direct and indirect searches of supersymmetric dark matter.

We propose a novel antiparticle detector. The gaseous antiparticle spectrometer (GAPS) effects particle identification through the characteristic X-rays emitted by antiparticles when they form exotic atoms in gases. GAPS obtains particularly high grasp (effective area--solid angle product) at lower particle energies, where conventional schemes are most limited in their utility. The concept is simple and lightweight, so it can be readily employed on balloon and space-based missions. An extremely powerful potential application of GAPS is a space-based search for the neutralino through the detection of a neutralino annihilation by-product -- the antideuteron. Paradoxically, this space-based search for the neutralino is capable of achieving comparable sensitivity to as yet unrealized third generation, underground dark matter experiments. And GAPS can obtain this performance in a very modest satellite experiment. GAPS can also provide superior performance in searches for primary antiprotons produced via neutralino annihilation, black hole evaporation and in probing subdominant contributions to the antiproton flux at low energies. In a deep space mission GAPS will obtain higher sensitivity for a given weight and power than BGO calorimeters.

Gravitinos with a mass in the keV range are an interesting candidate for warm dark matter. Recent measurements of the matter density of the universe and of cosmic structures at the dwarf galaxy scale rule out the simplest gravitino models with thermal freeze-out. We construct a model where the decay of the messenger particles that transmit the supersymmetry breaking to the observable sector generates the required entropy to dilute the gravitino relic density by the required factor of a few to come in line with observations. The model is natural, and requires only that the coupling of the messenger sector to the standard model be set so that the decay happens at the appropriate time.

Pixel microlensing, i.e. gravitational microlensing of unresolved stars, can be used to explore distant stellar systems, and as a bonus may be able to detect extragalactic planets. In these studies, binary-lens events with multiple high-magnification peaks are crucial. Considering only those events which exhibit caustic crossings, we estimate the fraction of binary events in several example pixel microlensing surveys and compare them to the fraction of binary events in a classical survey with resolved stars. We find a considerable enhancement of the relative rate of binary events in pixel microlensing surveys, relative to surveys with resolved sources. We calculate the rate distribution of binary events with respect to the time between caustic crossings. We consider possible surveys of M31 with ground-based telescopes and of M87 with HST and NGST. For the latter, a pixel microlensing survey taking one image a day may observe of order one dozen binary events per month.

The anomalous magnetic moment of the muon has recently been measured to be in conflict with the Standard Model prediction with an excess of 2.6 sigma. Taking this result as a measurement of the supersymmetric contribution, we find that at 95% confidence level it imposes an upper bound of about 500 GeV on the neutralino mass and forbids higgsino dark matter. More interestingly, it predicts an accessible lower bound on the direct detection rate, and it strongly favors models detectable by neutrino telescopes. Cosmic ray antideuterons may also be an interesting probe of such models.

The anomalous magnetic moment of the muon has recently been measured to be in conflict with the Standard Model prediction with an excess of 2.6 sigma. Taking the excess at face value as a measurement of the supersymmetric contribution, we find that at 95% confidence level it imposes an upper bound of 500 GeV on the neutralino mass and forbids higgsinos as being the bulk of cold dark matter. Other implications for the astrophysical detection of neutralinos include: an accessible minimum direct detection rate, lower bounds on the indirect detection rate of neutrinos from the Sun and the Earth, and a suppression of the intensity of gamma-ray lines from neutralino annihilations in the galactic halo.

We argue that gravitational microlensing is a feasible technique for measuring the mass function of brown dwarf stars in distant galaxies. Microlensing surveys of the bulge of M31, and of M87 in the Virgo cluster, may provide enough events to differentiate the behaviour of the mass function of lenses below the hydrogen burning limit (though we find that M87 is a more favourable target). Such objects may provide a significant supply of baryonic dark matter, an interesting possibility for the study of galactic dynamics. Furthermore, these systems have different metallicities than the solar neighbourhood, which may affect the mass function. These considerations are relevant in the context of star formation studies.

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. We here describe DarkSUSY, an advanced numerical FORTRAN package for supersymmetric dark matter calculations which we release for public use. With the help of this package, the masses and compositions of various supersymmetric particles can be computed, for given input parameters of the minimal supersymmetric extension of the Standard Model (MSSM). For the lightest neutralino, the relic density is computed, using accurate methods which include the effects of resonances, pair production thresholds and coannihilations. Accelerator bounds are checked to identify viable dark matter candidates. Finally, detection rates are computed for a variety of detection methods, such as direct detection and indirect detection through antiprotons, gamma-rays and positrons from the Galactic halo or neutrinos from the center of the Earth or the Sun.

Are there other planetary systems in our Universe? Indirect evidence has been found for planets orbiting other stars in our galaxy: the gravity of orbiting planets makes the star wobble, and the resulting periodic Doppler shifts have been detected for about a dozen stars. But are there planets in other galaxies, millions of light years away? Here we suggest a method to search for extragalactic planetary systems: gravitational microlensing of unresolved stars. This technique may allow us to discover planets in elliptical galaxies, known to be far older than our own Milky Way, with broad implications for life in the Universe.

We consider neutralino annihilation in dense extragalactic systems known to be dominated by dark matter, in particular M87 and several local dwarf spheroidal galaxies. These annihilations can produce energetic gamma rays which may be visible to atmospheric Cerenkov telescopes. We explore the supersymmetric parameter space, and compute the expected flux of gamma--rays coming from these objects. It is shown that some parts of the parameter space lead to a signal observable with the next generation of Cerenkov telescopes, provided the supersymmetric dark matter has a clumpy structure, as may be expected in a hierarchical scenario for structure formation.

It has been shown that a flux--weighted full width at half maximum timescale of a microlensing event can be used in an unbiased estimator of the optical depth. For the first time, this allows a physical parameter to be easily estimated from pixel microlensing data. We derive analytic expressions for the observed rate of pixel lensing events as a function of the full width at half maximum timescale. This contrasts work in the literature which express rates in terms of an ``event duration'' or Einstein time, which require knowledge of the magnification, which is difficult to determine in a pixel event. The full width at half maximum is the most directly measured timescale. We apply these results to possible pixel lensing surveys, using HST for M87, and CFHT for M31. We predict M87 microlensing rates for the HST Advanced Camera and for NGST, and demonstrate that one will be able to probe the stellar IMF. Next, we describe a new method by which a crude measurement of the magnification can be made in the regime of magnifications A~10-100. This in turn gives a crude measurement of the Einstein time. This program requires good photometry and sampling in the low magnification tails of an event, but is feasible with today's technology.

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.

We constrain the hadronic R-parity violating couplings in extensions to the minimal supersymmetric standard model. These interactions violate baryon and lepton number, and allow the lightest superpartner (LSP) to decay into standard model particles. The observed flux of cosmic ray antiprotons places a strong bound on the lifetime of the LSP in models where the lifetime is longer than the age of the universe. We exclude 10^-18<lambda''<10^-15 and 2x10^-18<lambda'<2x10^-15 except in the case of a top quark, where we can only exclude 4x10^-19<lambda'<4x10^-16.

Emission lines in hydrogen can be used to measure the approximate redshift of the reionization of the universe. This is an important measurement given the lack of a convincing theoretical prediction of this epoch. There is a rapid change in absorption at the reionization event, which can leave sharp features in the spectrum of recombining hydrogen. These features in the near infrared background permit us to measure the redshift at which reionization occurs. We have calculated these signals based on numerical simulations of the Lambda-CDM cosmogony.

Mica can store (for >1 Gy) etchable tracks caused by atoms recoiling from WIMPs. Ancient mica is a directional detector despite the complex motions it makes with respect to the WIMP "wind". We can exploit the properties of directionality and long integration time to probe for structure in the dark matter halo of our galaxy. We compute a sample of possible signals in mica for a plausible model of halo structure.