We study the radio signals generated when an axion star enters into the magnetosphere of a neutron star. As the axion star moves through the resonant region where the plasma-induced photon mass becomes equal to the axion mass, the axions can efficiently convert into photons, giving rise to an intense, transient radio signal. The energy released is determined by the axion star mass and conversion probability. Similarly, the peak frequency of the emitted radio signal is fixed by the axion mass, while cosmological redshift and Doppler shift could give rise to a wide range of frequencies. In particular, we show that a dense axion star with a mass $\sim 10^{-13}M_{\odot}$ composed of $\sim 10\: \mu$eV axions can account for most of the mysterious fast radio bursts in a wide frequency range.

We consider a gauge-singlet complex scalar field $\Phi$ with a global $U(1)$ symmetry that is spontaneously broken at some high energy scale $f_a$. As a result, the angular part of the $\Phi$-field becomes an axion-like particle (ALP). We show that if the $\Phi$-field has a non-zero coupling $\kappa$ to the Standard Model Higgs boson, there exists a certain region in the $\left(f_a, \kappa\right)$ parameter space where the global $U(1)$ symmetry-breaking induces a strongly first order phase transition, thereby producing stochastic gravitational waves that are potentially observable in current and future gravitational-wave detectors. In particular, we find that future gravitational-wave experiments such as TianQin, BBO and Cosmic Explorer could probe a broad range of the energy scale $10^3 \, {\rm GeV} \lesssim f_a \lesssim 10^{8} \, {\rm GeV}$, independent of the ALP mass. Since all the ALP couplings to the Standard Model particles are proportional to inverse powers of the energy scale $f_a$ (up to model-dependent ${\cal O}(1)$ coefficients), the gravitational-wave detection prospects are largely complementary to the current laboratory, astrophysical and cosmological probes of the ALP scenarios.

We study the generation of magnetic field seeds during a first-order electroweak phase transition, by numerically evolving the classical equations of motion of the bosonic electroweak theory on the lattice. The onset of the transition is implemented by the random nucleation of bubbles with an arbitrarily oriented Higgs field in the broken phase. We find that about 10% of the latent heat is converted into magnetic energy, with most of the magnetic fields being generated in the last stage of the phase transition when the Higgs oscillates around the true vacuum. The energy spectrum of the magnetic field has a peak that shifts towards larger length scales as the phase transition unfolds. By the end of our runs the peak wavelength is of the order of the bubble percolation scale, or about a third of our lattice size.

We propose a mechanism to generate Primordial Black Holes (PBHs) which is independent of cosmological inflation and occurs slightly below the QCD phase transition. Our setup relies on the collapse of long-lived string-domain wall networks and is naturally realized in QCD axion models with domain wall number $N_{DW}>1$ and Peccei-Quinn symmetry broken after inflation. In our framework, dark matter is mostly composed of axions in the meV mass range along with a small fraction, $\Omega_{\text{PBH}} \gtrsim 10^{-6} \Omega_{\text{CDM}} $ of heavy $M \sim 10^4-10^7 M_\odot$ PBHs. The latter could play a role in alleviating some of the shortcomings of the $\Lambda$CDM model on sub-galactic scales. The scenario has distinct signatures in ongoing axion searches as well as gravitational wave observatories.

TeV emission from blazars can be used to probe the intergalactic magnetic fields and measure their intensity, coherence length, and helicity. Intergalactic magnetic fields deflect the electron-positron pairs produced by very-high-energy gamma-rays from blazars, resulting in broadened beams of cascade gamma-rays developing along the projected direction of the blazar jet. We present an analysis that uses for the first time the jet orientation of 12 high-synchrotron-peaked (HSP) BL~Lac-type blazars from VLBA radio observations to search for signatures of pair haloes in the Fermi-LAT data. Our search improves the sensitivity of previous studies by taking the asymmetry of the pair haloes into account, increasing the signal to noise and reducing the possibility of systematics in determining the point spread function. Although there is no significant detection, a hint for an offset halo with a global significance of 2$\sigma$ is found in the $30-300$ GeV energy range, corresponding to an intergalactic magnetic field with $\sim 10^{-15}$ Gauss, consistent with the inferred field from a prior study using stacked HSP BL Lac objects.

The growth of a massive black hole will steepen the cold dark matter density at the center of a galaxy into a dense spike, enhancing the prospects for indirect detection. We study the impact of black hole spin on the density profile using the exact Kerr geometry of the black whole in a fully relativistic adiabatic growth framework. We find that, despite the transfer of angular momentum from the hole to the halo, rotation increases significantly the dark matter density close to the black hole. The gravitational effects are still dominated by the black hole within its influence radius, but the larger dark matter annihilation fluxes might be relevant for indirect detection estimates.

We investigate if the topology of pure gauge fields in the electroweak vacuum can play a role in classical dynamics at the electroweak phase transition. Our numerical analysis shows that magnetic fields are produced if the initial vacuum has non-trivial Chern-Simons number, and the fields are helical if the Chern-Simons number changes during the phase transition.

This white paper summarizes the workshop "U.S. Cosmic Visions: New Ideas in Dark Matter" held at University of Maryland on March 23-25, 2017.

We propose a method that allows to place an upper limit on the dark matter elastic scattering cross section with nucleons which is independent of the velocity distribution. Our approach combines null results from direct detection experiments with indirect searches at neutrino telescopes, and goes beyond previous attempts to remove astrophysical uncertainties in that it directly constrains the particle physics properties of the dark matter. The resulting halo-independent upper limits on the scattering cross section of dark matter are remarkably strong and reach $\sigma_{\text{SI}}^{p} \lesssim 10^{-43} \, (10^{-42})~\text{cm}^2$ and $\sigma_{\text{SD}}^{p} \lesssim 10^{-37} \, (3\times 10^{-37 })~\text{cm}^2$, for dark matter particles of $m_{\text{DM}}\sim 1~\text{TeV}$ annihilating into $W^+W^-$ ($b\bar b$), assuming $\rho_\text{loc}=0.3\text{ GeV}/\text{cm}^3$.

We report on the results of a search for $\gamma$-ray pair halos with a stacking analysis of low-redshift blazars using data from the Fermi Large Area Telescope. For this analysis we used a number of a-priori selection criteria, including the spatial and spectral properties of the Fermi sources. The angular distribution of $\sim$ 1GeV photons around 24 stacked isolated high-synchrotron-peaked BL Lacs with redshift $z<0.5$ shows an excess over that of point-like sources. A statistical analysis yields a Bayes factor of $\mathrm{log}_{10}B_{10}>2$, providing evidence in favor of extended emission against the point-source hypothesis, consistent with expectations for pair halos produced in the IGMF with strength $B_{\mathrm{IGMF}}\sim 10^{-17}-10^{-15}\mathrm{G}$.

The VERITAS array of Cherenkov telescopes has carried out a deep observational program on the nearby dwarf spheroidal galaxy Segue 1. We report on the results of nearly 48 hours of good quality selected data, taken between January 2010 and May 2011. No significant $\gamma$-ray emission is detected at the nominal position of Segue 1, and upper limits on the integrated flux are derived. According to recent studies, Segue 1 is the most dark matter-dominated dwarf spheroidal galaxy currently known. We derive stringent bounds on various annihilating and decaying dark matter particle models. The upper limits on the velocity-weighted annihilation cross-section are $\mathrm{<\sigma v >^{95% CL} \lesssim 10^{-23} cm^{3} s^{-1}}$, improving our limits from previous observations of dwarf spheroidal galaxies by at least a factor of two for dark matter particle masses $\mathrm{m_{\chi}\gtrsim 300 GeV}$. The lower limits on the decay lifetime are at the level of $\mathrm{\tau^{95% CL} \gtrsim 10^{24} s}$. Finally, we address the interpretation of the cosmic ray lepton anomalies measured by ATIC and PAMELA in terms of dark matter annihilation, and show that the VERITAS observations of Segue 1 disfavor such a scenario.

Non-vanishing parity-odd correlators of gamma ray arrival directions observed by Fermi-LAT indicate the existence of a helical intergalactic magnetic field (Tashiro et al.2013). We successfully test this hypothesis using more stringent cuts of the data, Monte Carlo simulations with Fermi-LAT time exposure information, separate analyses for the northern and southern galactic hemispheres, and confirm predictions made in Tashiro & Vachaspati (2014). With some further technical assumptions, we show how to reconstruct the magnetic helicity spectrum from the parity-odd correlators.

The existence of a cosmological magnetic field could be revealed by the effects of non-trivial helicity on large scales. We evaluate a CP odd statistic, $Q$, using gamma ray data obtained from Fermi satellite observations at high galactic latitudes to search for such a signature. Observed values of $Q$ are found to be non-zero; the probability of a similar signal in Monte Carlo simulations is $\sim 0.2\%$. Contamination from the Milky Way does not seem to be responsible for the signal since it is present even for data at very high galactic latitudes. Assuming that the signal is indeed due to a helical cosmological magnetic field, our results indicate left-handed magnetic helicity and field strength $\sim 10^{-14}~{\rm G}$ on $\sim 10~{\rm Mpc}$ scales.

The next generation of large scale WIMP direct detection experiments have the potential to go beyond the discovery phase and reveal detailed information about both the particle physics and astrophysics of dark matter. We report here on early results arising from the development of a detailed numerical code modeling the proposed DARWIN detector, involving both liquid argon and xenon targets. We incorporate realistic detector physics, particle physics and astrophysical uncertainties and demonstrate to what extent two targets with similar sensitivities can remove various degeneracies and allow a determination of dark matter cross sections and masses while also probing rough aspects of the dark matter phase space distribution. We find that, even assuming dominance of spin-independent scattering, multi-ton scale experiments still have degeneracies that depend sensitively on the dark matter mass, and on the possibility of isospin violation and inelasticity in interactions. We find that these experiments are best able to discriminate dark matter properties for dark matter masses less than around 200 GeV. In addition, and somewhat surprisingly, the use of two targets gives only a small improvement (aside from the advantage of different systematics associated with any claimed signal) in the ability to pin down dark matter parameters when compared with one target of larger exposure.

Report of the CF6 Working Group at Snowmass 2013. Topics addressed include ultra-high energy cosmic rays, neutrinos, gamma rays, baryogenesis, and experiments probing the fundamental nature of spacetime.

As part of the Snowmass process, the Cosmic Frontier Indirect-Detection subgroup (CF2) has drawn on input from the Cosmic Frontier and the broader Particle Physics community to produce this document. The purposes of this report are to identify opportunities for dark matter science through indirect detection, to give an overview of the primary scientific drivers for indirect searches for dark matter, and to survey current and planned experiments that have, as a large part of their scientific program, the goal of searching for indirect (or astrophysical) signatures of dark matter. We primarily address existing experiments with a large U.S. role, or future experiments where a U.S. contribution is sought. We also address the limitations of this technique, and answer the tough questions relevant to this subgroup posed by the HEP community through the Snowmass process.

We study the indirect detection of dark matter when the local dark matter velocity distribution depends upon position, as expected for the Milky Way and its dwarf spheroidal satellites, and the annihilation cross-section is not purely s-wave. Using a phase-space distribution consistent with the dark matter density profile, we present estimates of cosmic and gamma-ray fluxes from dark matter annihilations. The expectations for the indirect detection of dark matter can differ significantly from the usual calculation that assumes that the velocity of the dark matter particles follows a Maxwell-Boltzmann distribution.

The cold dark matter at the center of a galaxy will be redistributed by the presence of a massive black hole. The redistribution may be determined using an approach pioneered by Gondolo and Silk: begin with a model distribution function for the dark matter, and ``grow'' the black hole adiabatically, holding the adiabatic invariants of the motion constant. Unlike the approach of Gondolo and Silk, which adopted Newtonian theory together with ad hoc correction factors to mimic general relativistic effects, we carry out the calculation fully relativistically, using the exact Schwarzschild geometry of the black hole. We find that the density of dark matter generically vanishes at r=2R_S, not 4R_S as found by Gondolo and Silk, where R_S is the Schwarzschild radius, and that the spike very close to the black hole reaches significantly higher densities. We apply the relativistic adiabatic growth framework to obtain the final dark matter density for both cored and cusped initial distributions. Besides the implications of these results for indirect detection estimates, we show that the gravitational effects of such a dark matter spike are significantly smaller than the relativistic effects of the black hole, including frame dragging and quadrupolar effects, for stars orbiting close to the black hole that might be candidates for testing the black hole no-hair theorems.

The late-time accelerated expansion of the universe could be caused by a scalar field that is screened on small scales, as in chameleon or symmetron scenarios. We present an analogy between thin shell configurations of such scalar fields and electrostatics, which allows calculation of the field profile for general extended bodies. Interestingly, the field demonstrates a `lightning rod' effect, where it becomes enhanced near the ends of a pointy or elongated object. Drawing from this correspondence, we show that non-spherical test bodies immersed in a background field will experience a net torque caused by the scalar field. This effect, with no counterpart in the gravitational case, can be potentially tested in future experiments.

We derive new bounds on hidden sector gauge bosons which could produce new energy loss mechanisms in supernovae, enlarging the excluded region in mass-coupling space by a significant factor compared to earlier estimates. Both considerations of trapping and possible decay of these particles need to be incorporated when determining such bounds, as does scattering on both neutrons and protons. For masses and couplings near the region which saturates current bounds, a significant background of such gauge bosons may also be produced due to the cumulative effects of all supernovae over cosmic history.

The surface detector array of the Pierre Auger Observatory consists of 1600 water-Cherenkov detectors, for the study of extensive air showers (EAS) generated by ultra-high-energy cosmic rays. We describe the trigger hierarchy, from the identification of candidate showers at the level of a single detector, amongst a large background (mainly random single cosmic ray muons), up to the selection of real events and the rejection of random coincidences. Such trigger makes the surface detector array fully efficient for the detection of EAS with energy above $3\times 10^{18}$ eV, for all zenith angles between 0$^\circ$ and 60$^\circ$, independently of the position of the impact point and of the mass of the primary particle. In these range of energies and angles, the exposure of the surface array can be determined purely on the basis of the geometrical acceptance.

Noticeable deviations from the prediction of the fiducial LCDM cosmology are found in the angular power spectrum of the CMB. Besides large-angle anomalies, the WMAP 1st year data revealed a dip in the power spectrum at l \sim 200, which seemed to disappear in the 3rd year and subsequent angular power spectra. Using the WMAP 1st, 3rd, and 5th year data as well as the 5 year coadded data, we study the intensity and spatial distribution of this feature in order to unveil its origin and its implications for the cosmological parameters. We show that in all WMAP data releases there is a substantial suppression of the first Doppler peak in a region near the north ecliptic pole.

We explore the dark matter sector in extensions of the Minimal Supersymmetric Standard Model (MSSM) that can provide a good fit to the PAMELA cosmic ray positron excess, while at the same time addressing the little hierarchy problem of the MSSM. Adding a singlet Higgs superfield, S, can account for the observed positron excess, as recently discussed in the literature, but we point out that it requires a fine-tuned choice for the parameters of the model. We find that including an additional singlet allows both a reduction of the weak-scale fine-tuning, and an interpretation of the cosmic ray observations in terms of dark matter annihilations in the galactic halo. Our setup contains a light axion, but does not require light CP-even scalars in the spectrum.

If the dark matter sector in the universe is composed by metastable particles, galaxies and galaxy clusters are expected to undergo significant secular evolution from high to low redshift. We show that the decay of dark matter, with a lifetime compatible with cosmological constraints, can be at the origin of the observed evolution of the Tully-Fisher relation of disk galaxies and alleviate the problem of the size-evolution of elliptical galaxies, while being consistent with the current observational constraints on the gas fraction of clusters of galaxies.

The Pierre Auger Observatory is a hybrid detector for ultra-high energy cosmic rays. It combines a surface array to measure secondary particles at ground level together with a fluorescence detector to measure the development of air showers in the atmosphere above the array. The fluorescence detector comprises 24 large telescopes specialized for measuring the nitrogen fluorescence caused by charged particles of cosmic ray air showers. In this paper we describe the components of the fluorescence detector including its optical system, the design of the camera, the electronics, and the systems for relative and absolute calibration. We also discuss the operation and the monitoring of the detector. Finally, we evaluate the detector performance and precision of shower reconstructions.

Studies of the cosmic ray energy spectrum at the highest energies with the Pierre Auger Observatory.

Technical reports on operations and features of the Pierre Auger Observatory, including ongoing and planned enhancements and the status of the future northern hemisphere portion of the Observatory. Contributions to the 31st International Cosmic Ray Conference, Lodz, Poland, July 2009.

Studies of the correlations of ultra-high energy cosmic ray directions with extra-Galactic objects, of general anisotropy, of photons and neutrinos, and of other astrophysical effects, with the Pierre Auger Observatory. Contributions to the 31st ICRC, Lodz, Poland, July 2009.

Reports on the atmospheric monitoring, calibration, and other operating systems of the Pierre Auger Observatory. Contributions to the 31st International Cosmic Ray Conference, Lodz, Poland, July 2009.

Studies of the composition of the highest energy cosmic rays with the Pierre Auger Observatory, including examination of hadronic physics effects on the structure of extensive air showers. Submissions to the 31st ICRC, Lodz, Poland (July 2009).

From direct observations of the longitudinal development of ultra-high energy air showers performed with the Pierre Auger Observatory, upper limits of 3.8%, 2.4%, 3.5% and 11.7% (at 95% c.l.) are obtained on the fraction of cosmic-ray photons above 2, 3, 5 and 10 EeV (1 EeV = 10^18 eV) respectively. These are the first experimental limits on ultra-high energy photons at energies below 10 EeV. The results complement previous constraints on top-down models from array data and they reduce systematic uncertainties in the interpretation of shower data in terms of primary flux, nuclear composition and proton-air cross-section.

In string gas cosmology, extra dimensions are stabilised by a gas of strings. In the matter-dominated era, competition between matter pushing the extra dimensions to expand and the string gas pulling them back can lead to oscillations of the extra dimensions and acceleration in the visible dimensions. We fit this model to supernova data, taking into account the Big Bang Nucleosynthesis constraint on the energy density of the string gas. The fit to the Union set of supernova data is acceptable, but the fit to the ESSENCE data is poor.

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.

Many models of baryogenesis rely on anomalous particle physics processes to give baryon number violation. By numerically evolving the electroweak equations on a lattice, we show that baryogenesis in these models creates helical cosmic magnetic fields. After a transitory period, electroweak dynamics is found to conserve the Chern-Simons number and the total electromagnetic helicity. We argue that baryogenesis could lead to magnetic fields of nano-Gauss strength today on astrophysical length scales. In addition to being astrophysically relevant, such helical magnetic fields can provide an independent probe of baryogenesis and CP violation in particle physics.

We discuss an inflationary scenario based on Lovelock terms. These higher order curvature terms can lead to inflation when there are more than three spatial dimensions. Inflation will end if the extra dimensions are stabilised, so that at most three dimensions are free to expand. This relates graceful exit to the number of large dimensions.

We explore the prospects for indirect detection of neutralino dark matter in supersymmetric models with an extended Higgs sector (NMSSM). We compute, for the first time, one-loop amplitudes for NMSSM neutralino pair annihilation into two photons and two gluons, and point out that extra diagrams (with respect to the MSSM), featuring a potentially light CP-odd Higgs boson exchange, can strongly enhance these radiative modes. Expected signals in neutrino telescopes due to the annihilation of relic neutralinos in the Sun and in the Earth are evaluated, as well as the prospects of detection of a neutralino annihilation signal in space-based gamma-ray, antiproton and positron search experiments, and at low-energy antideuteron searches. We find that in the low mass regime the signals from capture in the Earth are enhanced compared to the MSSM, and that NMSSM neutralinos have a remote possibility of affecting solar dynamics. Also, antimatter experiments are an excellent probe of galactic NMSSM dark matter. We also find enhanced two photon decay modes that make the possibility of the detection of a monochromatic gamma-ray line within the NMSSM more promising than in the MSSM.

Observations of the Milky Way by the SPI/INTEGRAL satellite have confirmed the presence of a strong 511 KeV gamma-ray line emission from the bulge, which require an intense source of positrons in the galactic center. These observations are hard to account for by conventional astrophysical scenarios, whereas other proposals, such as light DM, face stringent constraints from the diffuse gamma-ray background. Here we suggest that light superconducting strings could be the source of the observed 511 KeV emission. The associated particle physics, at the ~ 1 TeV scale, is within reach of planned accelerator experiments, while the distinguishing spatial distribution, proportional to the galactic magnetic field, could be mapped by SPI or by future, more sensitive, satellite missions.

A classical analysis suggests that an external magnetic field can cause trajectories of charge carriers on a superconducting domain wall or cosmic string to bend, thus expelling charge carriers with energy above the mass threshold into the bulk. We study this process by solving the Dirac equation for a fermion of mass $m_f$ and charge $e$, in the background of a domain wall and a magnetic field of strength $B$. We find that the modes of the charge carriers get shifted into the bulk, in agreement with classical expectations. However the dispersion relation for the zero modes changes dramatically -- instead of the usual linear dispersion relation, $\omega_k =k$, the new dispersion relation is well fit by $\omega \approx m_f tanh(k/k_*)$ where $k_*=m_f$ for a thin wall in the weak field limit, and $k_*=eBw$ for a thick wall of width $w$. This result shows that the energy of the charge carriers on the domain wall remains below the threshold for expulsion even in the presence of an external magnetic field. If charge carriers are expelled due to an additional perturbation, they are most likely to be ejected at the threshold energy $\sim m_f$.

We study the stability of extra dimensions in string gas cosmology at late times. Vacuum energy and, interestingly, baryons lead to decompactification after they become dynamically important. The string gas can stabilise the effect of baryons, but not that of vacuum energy. However, we find that the interplay of baryons and strings can lead to acceleration in the visible dimensions, without the need for vacuum energy.

We show that a tangle of light superconducting strings in the Milky Way could be the source of the observed 511 KeV emission from electron-positron annihilation in the Galactic bulge. The scenario predicts a flux that is in agreement with observations if the strings are at the ~ 1 TeV scale making the particle physics within reach of planned accelerator experiments. The emission is directly proportional to the galactic magnetic field and future observations should be able to differentiate the superconducting string scenario from other proposals.

Ground based Atmospheric Cerenkov Telescopes have recently unveiled a TeV gamma-ray signal from the direction of the Galactic Centre. We examine whether these gamma-rays, observed by the VERITAS, CANGAROO-II and HESS collaborations, may arise from annihilations of dark matter particles. Emission from nearby dwarf spheroidals, such as Sagittarius, could provide a test of this scenario.

Scalar mediated interactions among baryons extend well above the Compton wavelength, when they are embedded in a Bose-Einstein condensate composed of the mediating particles. Indeed, this non-trivial environment results in an infinite-ranged interaction. We show that if the Dark Matter of the Universe is composed of such a condensate, the imprints of an interaction between baryonic and Dark Matter could be manifest as anomalies in the peak structure of the Cosmic Microwave Background.

It has been recently proposed that the observed 511 keV emission from the Galactic bulge could be the product of very light (1-100 MeV) annihilating dark matter particles. Other possible explanations for this signal are associated with stellar objects, such as hypernovae. In order to distinguish between annihilating light dark matter scenario and more conventional astrophysical sources for the bulge emission, we here propose the study of dwarf spheroidals such as Sagittarius. These galaxies have typical luminosities of $10^5 L_\odot$ but mass-to-light ratios of $\sim 100$. As there are comparatively few stars, the prospects for 511 keV emission from standard astrophysical scenarios are minimal. The dwarf spheroidals do, however, contain copious amounts of dark matter. INTEGRAL/SPI has observed the Sagittarius region. Analysis of this data for 511 keV emission will provide a test of MeV dark matter which can distinguish between annihilating dark matter and more standard astrophysics. The observation of such a signal from Sagittarius should be a ``smoking gun'' for MeV dark matter.

We study cosmological perturbations in the case that present-day matter consists of a mixture of inflaton and curvaton decay products. We calculate how the curvaton perturbations are transferred to its decay products in the general case when it does not behave like dust. Taking into account that the decay products of the inflaton can also have perturbations results in an interesting mixture of correlated adiabatic and isocurvature perturbations. In particular, negative correlation can improve the fit to the CMB data by lowering the angular power in the Sachs-Wolfe plateau without changing the peak structure. We do an 11-parameter fit to the WMAP data. We find that the best-fit is not the 'concordance model', and that well-fitting models do not cluster around the best-fit, so that cosmological parameters cannot be reliably estimated. We also find that in our model the mean quadrupole (l=2) power is l(l+1) C_l/2pi = 1081 muK^2, much lower than in the pure adiabatic LCDM model, which gives 1262 muK^2.

We calculate the flux of high energy gamma-rays from annihilation of neutralino dark matter in the centre of the Milky Way and the three nearest dwarf spheroidals (Sagittarius, Draco and Canis Major), using realistic models of the dark matter distribution.

Two ground-based experiments have recently independently detected TeV $\gamma$-rays from the direction of the Galactic center. The observations made by the VERITAS and CANGAROO collaborations are unexpected, although not impossible to interpret in terms of astrophysical sources. Here we examine in detail whether the observed $\gamma$-rays may arise from the more exotic alternative of annihilations of dark matter particles clustered in the center of the Galaxy.

We reiterate that there is no evidence that BL Lacs are sources of ultrahigh energy cosmic rays.

We provide a ``Baedecker'' or travel guide to the directions on the sky where the dark matter annihilation signal may be expected. We calculate the flux of high energy gamma-rays from annihilation of neutralino dark matter in the centre of the Milky Way and the three nearest dwarf spheroidals (Sagittarius, Draco and Canis Major), using realistic models of the dark matter distribution. Other investigators have used cusped dark halo profiles (such as the Navarro-Frenk-White) to claim a significant signal. This ignores the substantial astrophysical evidence that the Milky Way is not dark-matter dominated in the inner regions. We show that the annihilation signal from the Galactic Centre falls by two orders of magnitude on substituting a cored dark matter density profile for a cusped one. The present and future generation of high energy gamma-ray detectors, whether atmospheric Cerenkov telescopes or space missions like GLAST, lack the sensitivity to detect any of the monochromatic gamma-ray annihilation lines. The continuum gamma-ray signal above 1 GeV and above 50 GeV may however be detectable either from the dwarf spheroidals or from the Milky Way itself. If the density profiles of the dwarf spheroidals are cusped, then the best prospects are for detecting Sagittarius and Canis Major. However, if the dwarf spheroidals have milder, cored profiles, then the annihilation signal is not detectable. For GLAST, an attractive strategy is to exploit the wide field of view and observe the Milky Way at medium latitudes, as suggested by Stoehr et al. This is reasonably robust against changes in the density profile.

The sky distribution of cosmic rays with energies above the 'GZK cutoff' holds important clues to their origin. The AGASA data, although consistent with isotropy, shows evidence for small-angle clustering, and it has been argued that such clusters are aligned with BL Lacertae objects, implicating these as sources. It has also been suggested that clusters can arise if the cosmic rays come from the decays of very massive relic particles in the Galactic halo, due to the expected clumping of cold dark matter. We examine these claims and show that both are in fact not justified.

Ultra-high energy cosmic rays (UHECRs) may originate from the decay of massive relic particles in the dark halo of the Galaxy, or they may be produced by supermassive black holes in the nuclei of nearby galaxies. The anisotropy in the arrival directions is studied in four dark halo models (cusped, isothermal, triaxial and tilted) and in four galaxy samples (galaxies intrinsically brighter than Centaurus A within 50 and 100 Mpc, and brighter than M32 within 50 and 100 Mpc). In decaying dark matter models, the amplitude of the anisotropy is controlled by the size of the Galactic halo, while the phase is controlled by the shape. In the northern hemisphere, the amplitude is about 0.5 for cusped haloes, but falls to roughly 0.3 for isothermal haloes. The phase points in the direction of the Galactic Centre, with deviations of up to 30 degrees possible for triaxial and tilted haloes. The effect of the halo of M31 is too weak to provide conclusive evidence for the decaying dark matter origin of UHECRs. In extragalactic models, samples of galaxies brighter than Centaurus A produce substantial anisotropies (roughly 1.8), much larger than the limits set by the available data. If all galaxies brighter than M32 contribute, then the anisotropy is more modest (roughly 0.5) and is directed toward mass concentrations in the supergalactic plane, like the Virgo cluster. Predictions are made for the southern hemisphere station of the Pierre Auger Observatory. If the UHECRs have a Galactic origin, then the phase points towards the Galactic Centre. If they have an extragalactic origin, then it points in the rough direction of the Fornax cluster. This provides an unambiguous discriminant and requires about 350-500 events at South Auger.