Annihilation of positrons in the Galaxy

Observations of electron-positron annihilation radiation from the Galactic Center region with the SPI instrument aboard INTEGRAL are summarized. The measured width of the 511 keV line and inferred fraction of positrons annihilating through positronium formation are consistent with the annihilation taking place in the warm ISM phase, although combinations of several ISM phases are also allowed by the data. The spatial distribution of 511 keV emission suggests that positron sources are concentrated toward the Galactic bulge and avoid the Galactic disk.     

On the origin of the 511-keV emission in the Galactic Centre

Diffuse 511-keV line emission, from the annihilation of cold positrons, has been observed in the direction of the Galactic Centre for more than 30 yr. The latest high-resolution maps of this emission produced by the SPI instrument on INTEGRAL suggest at least one component of the emission is spatially coincident with the distribution of ∼70 luminous, low-mass X-ray binaries detected in the soft gamma-ray band. The X-ray band, however, is generally a more sensitive probe of X-ray binary populations. Recent X-ray surveys of the Galactic Centre have discovered a much larger population (>4000) of faint, hard X-ray point sources. We investigate the possibility that the positrons observed in the direction of the Galactic Centre originate in pair-dominated jets generated by this population of fainter accretion-powered X-ray binaries. We also consider briefly whether such sources could account for unexplained diffuse emission associated with the Galactic Centre in the microwave (the Wilkinson Microwave Anisotropy Probe 'haze') and at other wavelengths. Finally, we point out several unresolved problems in associating Galactic Centre 511-keV emission with the brightest X-ray binaries.     

Constraining dark matter properties with SPI

Using the high-resolution spectrometer SPI on board the International Gamma-Ray Astrophysics Laboratory ( INTEGRAL ), we search for a spectral line produced by a dark matter (DM) particle with a mass in the range  40 keV < M _DM < 14 MeV  , decaying in the DM halo of the Milky Way. To distinguish the DM decay line from numerous instrumental lines found in the SPI background spectrum, we study the dependence of the intensity of the line signal on the offset of the SPI pointing from the direction toward the Galactic Centre. After a critical analysis of the uncertainties of the DM density profile in the inner Galaxy, we find that the intensity of the DM decay line should decrease by at least a factor of 3 when the offset from the Galactic Centre increases from 0° to 180°. We find that such a pronounced variation of the line flux across the sky is not observed for any line, detected with a significance higher than 3σ in the SPI background spectrum. Possible DM decay origin is not ruled out only for the unidentified spectral lines, having low (∼3σ) significance or coinciding in position with the instrumental ones. In the energy interval from 20 keV to 7 MeV, we derive restrictions on the DM decay line flux, implied by the (non-)detection of the DM decay line. For a particular DM candidate, the sterile neutrino of mass M _DM, we derive a bound on the mixing angle.     

Annihilation radiation from a dark matter spike at the Galactic Centre

We study annihilation radiation of neutralinos in the Galactic Centre, assuming the existence of a 'spike' in the dark matter density profile, caused by adiabatic accretion on to the massive black hole lying at the Galactic Centre. Under this assumption we find that it is possible to reproduce the observed SgrA* emission at radio and gamma-ray frequencies in a consistent scenario with a magnetic field close to the equipartition strength and with values of γ (the density profile power-law index) around 0.1.     

Dark stars at the Galactic Centre – the main sequence

In regions of very high dark matter density such as the Galactic Centre, the capture and annihilation of WIMP dark matter by stars has the potential to significantly alter their evolution. We describe the dark stellar evolution code D ark S tars , and present a series of detailed grids of WIMP-influenced stellar models for main-sequence stars. We describe the changes in stellar structure and main-sequence evolution which occur as a function of the rate of energy injection by WIMPs, for masses of  0.3–2.0 M_⊙  and metallicities   Z = 0.0003–0.02  . We show what rates of energy injection can be obtained using realistic orbital parameters for stars at the Galactic Centre, including detailed consideration of the velocity and density profiles of dark matter. Capture and annihilation rates are strongly boosted when stars follow elliptical rather than circular orbits. If there is a spike of dark matter induced by the supermassive black hole at the Galactic Centre, single solar mass stars following orbits with periods as long as 50 yr and eccentricities as low as 0.9 could be significantly affected. Binary systems with similar periods about the Galactic Centre could be affected on even less eccentric orbits. The most striking observational effect of this scenario would be the existence of a binary consisting of a low-mass protostar and a higher mass evolved star. The observation of low-mass stars and/or binaries on such orbits would either provide a detection of WIMP dark matter, or place stringent limits on the combination of the WIMP mass, spin-dependent nuclear-scattering cross-section, halo density and velocity distribution near the Galactic Centre. In some cases, the derived limits on the WIMP mass and spin-dependent nuclear-scattering cross-section would be of comparable sensitivity to current direct-detection experiments.     

Cuspy dark matter haloes and the Galaxy

The microlensing optical depth to Baade's Window constrains the minimum total mass in baryonic matter within the Solar circle to be greater than ∼     , assuming the inner Galaxy is barred with viewing angle ∼20°. From the kinematics of solar neighbourhood stars, the local surface density of dark matter is ∼     . We construct cuspy haloes normalized to the local dark matter density and calculate the circular-speed curve of the halo in the inner Galaxy. This is added in quadrature to the rotation curve provided by the stellar and ISM discs, together with a bar sufficiently massive so that the baryonic matter in the inner Galaxy reproduces the microlensing optical depth. Such models violate the observational constraint provided by the tangent-velocity data in the inner Galaxy (typically at radii     . The high baryonic contribution required by the microlensing is consistent with implications from hydrodynamical modelling and the pattern speed of the Galactic bar. We conclude that the cuspy haloes favoured by the cold dark matter cosmology (and its variants) are inconsistent with the observational data on the Galaxy.     

On the prospect of inferring the halo structure and the masses of dark objects through parallax microlensing

We study the proposed use of parallax microlensing in the direction of the Large Magellanic Cloud (LMC) to separate the effects of the mass function of dark massive halo objects (MHOs or 'machos') on the one hand, and their spatial distribution and kinematics on the other. This disentanglement is supposed to allow a much better determination of the two than could be achieved entirely on the basis of the durations of events. We restrict our treatment to the same class of power-law spherical models for the halo of MHOs studied in a previous paper by Marković 38 Sommer-Larsen, and assume that one can eliminate microlensing events caused by massive objects outside the halo (e.g., the LMC halo). Whereas the duration-based error in the average MHO mass, μ¯ ≡  M ¯/M, exceeds (at N  = 100 events) μ¯ by a factor of 2 or more, parallax microlensing remarkably brings it down to 15–20 per cent of μ¯, regardless of the shape of the mass function. In addition, the slope α of the mass function, d n /dμ ∝ μ^α, can be inferred relatively accurately (σ_α < 0.4) for a broader range, −3 < α < 0. The improvement in the inference of the halo structure is also significant: the index γ of the density profile ( ρ ∼  R ^−γ) can be obtained with the error σ_γ < 0.4. While in a typical situation the errors for the parameters specifying the velocity dispersion profile are of about the same magnitude as the parameters themselves, virtually all the uncertainty is 'concentrated' in linear combinations of the parameters that may have little influence on the profile, thus allowing its reasonably accurate inference.     

Constraints on decaying dark matter from XMM–Newton observations of M31

We derive constraints on the parameters of the radiatively decaying dark matter (DM) particle, using the XMM–Newton EPIC spectra of the Andromeda galaxy (M31). Using the observations of the outer (5–13 arcmin) parts of M31, we improve the existing constraints. For the case of sterile neutrino DM, combining our constraints with the latest computation of abundances of sterile neutrinos in the Dodelson–Widrow (DW) scenario, we obtain the lower mass limit   m _s < 4 keV  , which is stronger than the previous one   m _s < 6 keV  , obtained recently by Asaka, Laine & Shaposhnikov. Comparing this limit with the most recent results on Lyman α forest analysis of Viel et al.  ( m _s > 5.6 keV  ), we argue that the scenario in which all the DM is produced via the DW mechanism is ruled out. We discuss, however, other production mechanisms and note that the sterile neutrino remains a viable candidate for DM, either warm or cold.     

Gamma rays from the Galactic Centre region: A review

During the last decades, increasingly precise astronomical observations of the Galactic Centre (GC) region at radio, infrared, and X-ray wavelengths laid the foundations to a detailed understanding of the high energy astroparticle physics of this most remarkable location in the Galaxy. Recently, observations of this region in high energy (HE, 10 MeV–100 GeV) and very high energy (VHE, > 100 GeV) γ-rays added important insights to the emerging picture of the Galactic nucleus as a most violent and active region where acceleration of particles to very high energies – possibly up to a PeV – and their transport can be studied in great detail. Moreover, the inner Galaxy is believed to host large concentrations of dark matter (DM), and is therefore one of the prime targets for the indirect search for γ-rays from annihilating or decaying dark matter particles. In this article, the current understanding of the γ-ray emission emanating from the GC is summarised and the results of recent DM searches in HE and VHE γ-rays are reviewed.     

Analytical approach to subhalo population in dark matter haloes

In the standard model of cosmic structure formation, dark matter haloes form by gravitational instability. The process is hierarchical: smaller systems collapse earlier, and later merge to form larger haloes. The galaxy clusters, hosted by the largest dark matter haloes, are at the top of this hierarchy and representing the largest as well as the last structures formed in the Universe, while the smaller and first haloes are those Earth-sized dark subhaloes that have been both predicted by theoretical considerations and found in numerical simulations, though there do not exist any observational hints of their existence. The probability that a halo of mass m at redshift z will be part of a larger halo of mass M at the present time can be described in the frame of the extended Press & Schecter theory making use of the progenitor (conditional) mass function. Using the progenitor mass function, we calculate analytically, at redshift zero, the distribution of subhaloes in mass, formation epoch and rarity of the peak of the density field at the formation epoch. That is done for a Milky Way size system, assuming both a spherical and an ellipsoidal collapse model. Our calculation assumes that small progenitors do not lose mass due to dynamical processes after entering the parent halo, and that they do not interact with other subhaloes. For a Λ cold dark matter power spectrum, we obtain a subhalo mass function  d n /d m   proportional to   m ^−α  with a model-independent  α∼ 2  . Assuming that the dark matter is a weakly interacting massive particle, the inferred distributions are used to test the feasibility of an indirect detection in the γ-ray energy band of such a population of subhaloes with a Gamma-ray Large Area Space Telescope like satellite.     

Diffuse Galactic γ-rays: Constraining Cosmic-Ray Origin and Propagation

We have developed a model which aims to reproduce observational data of many kinds related to cosmic-ray (CR) origin and propagation: direct measurements of nuclei, antiprotons, electrons and positrons,γ-rays, and synchrotron radiation. Our main results include evaluation of diffusion/convection and reacceleration models, estimates of the halo size, calculations of the interstellar positron and antiproton spectra, evaluation of alternative hypotheses of nucleon and electron interstellar spectra, and computation of the Galactic diffuseγ-ray emission. Recently our CR propagation code has been generalized to include fragmentation networks of arbitrary complexity. The code can now provide an alternative to leaky-box calculations for full isotopic abundance calculations and has the advantage of including the spatial dimension which is essential for radioactive nuclei. Preliminary predictions for sub-Fe/Fe, ¹⁰Be/⁹Be and ²⁶Al/²⁷Al are presented in anticipation of new experimental isotopic data. We show that combining information from classical CR studies with γ-ray and other data leads to tighter constraints on CR origin and propagation. This revised version was published online in July 2006 with corrections to the Cover Date.     

The Planetary Nebula Spectrograph elliptical galaxy survey: the dark matter in NGC 4494

We present new Planetary Nebula Spectrograph observations of the ordinary elliptical galaxy NGC 4494, resulting in positions and velocities of 255 planetary nebulae out to seven effective radii (25 kpc). We also present new wide-field surface photometry from MMT/Megacam, and long-slit stellar kinematics from VLT/FORS2. The spatial and kinematical distributions of the planetary nebulae agree with the field stars in the region of overlap. The mean rotation is relatively low, with a possible kinematic axis twist outside  1 R _e  . The velocity dispersion profile declines with radius, though not very steeply, down to  ∼70 km s⁻¹  at the last data point.
We have constructed spherical dynamical models of the system, including Jeans analyses with multi-component Λ cold dark matter (CDM) motivated galaxies as well as logarithmic potentials. These models include special attention to orbital anisotropy, which we constrain using fourth-order velocity moments. Given several different sets of modelling methods and assumptions, we find consistent results for the mass profile within the radial range constrained by the data. Some dark matter (DM) is required by the data; our best-fitting solution has a radially anisotropic stellar halo, a plausible stellar mass-to-light ratio and a DM halo with an unexpectedly low central density. We find that this result does not substantially change with a flattened axisymmetric model.
Taken together with other results for galaxy halo masses, we find suggestions for a puzzling pattern wherein most intermediate-luminosity galaxies have very low concentration haloes, while some high-mass ellipticals have very high concentrations. We discuss some possible implications of these results for DM and galaxy formation.     

Constraining the dark matter annihilation cross-section with Cherenkov telescope observations of dwarf galaxies

The presence of dark matter in the halo of our Galaxy could be revealed through indirect detection of its annihilation products. Dark matter annihilation is one possible interpretation of the recently measured excesses in positron and electron fluxes, provided that boost factors of the order of 10³ or more are taken into account. Such boost factors are actually achievable through the velocity-dependent Sommerfeld enhancement of the annihilation cross-section. Here, we study the expected γ-ray flux from two local dwarf galaxies for which Cherenkov telescope measurements are available, namely Draco and Sagittarius. We use recent stellar kinematical measurements to model the dark matter haloes of the dwarfs and the results of numerical simulations to model the presence of an associated population of subhaloes. We incorporate the Sommerfeld enhancement of the annihilation cross-section. We compare our predictions with the observations of Draco and Sagittarius performed by MAGIC and HESS, respectively, and derive exclusion limits on the effective annihilation cross-section. We also study the sensitivities of Fermi and of the future Cherenkov telescope array to cross-section enhancements. We find that the boost factor due to the Sommerfeld enhancement is already constrained by the MAGIC and HESS data, with enhancements greater than ∼10⁴ being excluded.     

On the interaction between cosmic rays and dark matter molecular clouds in the Milky Way – I. Basic considerations

We explore some basic observational consequences of assuming that the dark matter in the Milky Way consists mainly of molecular clouds, and that cosmic rays can penetrate these clouds. In a favoured model of the clouds, this penetration would have the following consequences, all of which agree with observation.
(i) Cosmic ray nuclei would be fragmented when they enter a cloud, giving them a lifetime in the Galaxy of ∼10¹⁵ s (for relativistic nuclei).
(ii) Pionic γ -rays emitted by the clouds, after proton–proton (pp) collisions, would have a diffuse flux in the Galactic plane comparable to the flux from known sources for photon energies ≳1 GeV .
(iii) The heat input into the clouds from cosmic rays would be re-radiated mainly in the far-infrared. The resulting radiation background agrees, in both intensity and spectrum in different directions, with a known excess in the far‐infrared background of the galaxy over emission by warm dust.     

Numerical action reconstruction of the dynamical history of dark matter haloes in N-body simulations

We test the ability of the numerical action method (NAM) to recover the individual orbit histories of mass tracers in an expanding universe, given the masses and redshift-space coordinates at the present epoch. The mass tracers are represented by dark matter (DM) haloes identified in a region of radius  26  h ⁻¹ Mpc  of a high-resolution N -body simulation of the standard Λ cold dark matter (CDM) cosmology. Since previous tests of NAM at this scale have traced the underlying distribution of DM particles rather than extended haloes, our study offers an assessment of the accuracy of NAM in a scenario which more closely approximates the complex dynamics of actual galaxy haloes. We show that NAM can recover the present-day halo distances with typical distance errors of less than 3 per cent and radial peculiar velocities with a dispersion of  ∼130 km s⁻¹  . The accuracy of individual orbit reconstructions was limited by the inability of NAM, in some instances, to correctly model the positions of haloes at early times solely on the basis of the redshifts, angular positions and masses of the haloes at the present epoch. Improvements in the quality of NAM reconstructions may be possible using the present-day three-dimensional halo velocities and distances to further constrain the dynamics. This velocity data is expected to become available for nearby galaxies in the coming generations of observations by Space Interferometry Mission ( SIM ) and Global Astrometric Interferometer for Astrophysics ( GAIA ).     

The signature of dark energy on the local Hubble flow

Using N -body simulations of flat, dark energy-dominated cosmologies, we show that galaxies around simulated binary systems resembling the Local Group (LG) have low peculiar velocities, in good agreement with observational data. We have compared results for LG-like systems selected from large, high-resolution simulations of three cosmologies: a ΛCDM model, a ΛWDM model with a 2-keV warm dark matter candidate, and a quintessence (QCDM) model with an equation-of-state parameter   w =−0.6  . The Hubble flow is significantly colder around LGs selected in a flat, Λ-dominated cosmology than around LGs in open or critical models, showing that a dark energy component manifests itself on the scales of nearby galaxies, cooling galaxy peculiar motions. Flows in the ΛWDM and QCDM models are marginally colder than in the ΛCDM one.
The results of our simulations have been compared to existing data and to a new data set of 28 nearby galaxies with robust distance measures (Cepheids and surface brightness fluctuations). The measured line-of-sight velocity dispersion is given by  σ_H= (88 ± 20  km s⁻¹) × ( R /7 Mpc)  . The best agreement with observations is found for LGs selected in the ΛCDM cosmology in environments with  −0.1 < δρ/ρ < 0.6  on scales of 7 Mpc, in agreement with existing observational estimates on the local matter density. These results provide new, independent evidence for the presence of dark energy on scales of a few megaparsecs, corroborating the evidence gathered from observations of distant objects and the early Universe.     

Searching for neutrinos from WIMP annihilations in the Galactic stellar disc

Weakly interacting massive particles (WIMPs) are a viable candidate for the relic abundance of dark matter (DM) produced in the early universe. So far, WIMPs have eluded direct detection through interactions with baryonic matter. Neutrino emission from accumulated WIMP annihilations in the solar core has been proposed as a signature of DM, but has not yet been detected. These null results may be due to small-scale DM density fluctuations in the halo with the density of our local region being lower than the average  (∼0.3 GeV cm⁻³)  . However, the accumulated neutrino signal from WIMP annihilations in the Galactic stellar disc would be insensitive to local density variations. Inside the disc, DM can be captured by stars causing an enhanced annihilation rate and therefore a potentially higher neutrino flux than what would be observed from elsewhere in the halo. We estimate a neutrino flux from the WIMP annihilations in the stellar disc to be enhanced by more than an order of magnitude compared to the neutrino fluxes from the halo. We offer a conservative estimate for this enhanced flux, based on the WIMP–nucleon cross-sections obtained from direct-detection experiments by assuming a density of  ∼0.3 GeV cm⁻³  for the local DM. We also compare the detectability of these fluxes with a signal of diffuse high-energy neutrinos produced in the Milky Way by the interaction of cosmic rays with the interstellar medium. These comparative signals should be observable by large neutrino detectors.     

Gamma-ray emission from individual classical novae

Classical novae are important producers of radioactive nuclei, such as ⁷Be, ¹³N, ¹⁸F, ²²Na and ²⁶Al. The disintegration of these nuclei produces positrons (except for ⁷Be) that through annihilation with electrons produce photons of energies 511 keV and below. Furthermore, ⁷Be and ²²Na decay producing photons with energies of 478 and 1275 keV, respectively, well in the γ-ray domain. Therefore, novae are potential sources of γ-ray emission. We have developed two codes in order to analyse carefully the γ-ray emission of individual classical novae: a hydrodynamical one, which follows both the accretion and the explosion stages, and a Monte Carlo one, able to treat both production and transfer of γ-ray photons. Both codes have been coupled in order to simulate realistic explosions. The properties of γ-ray spectra and γ-ray light curves (for the continuum and for the lines at 511, 478 and 1275 keV) have been analysed, with a special emphasis on the difference between carbon–oxygen and oxygen–neon novae. Predictions of detectability of individual novae by the future SPI spectrometer on board the INTEGRAL satellite are made. Concerning ²⁶Al, its decay produces photons of 1809 keV but this occurs on a time-scale much longer than the typical time interval between nova outbursts in the Galaxy, making it undetectable in individual novae. The accumulated emission of ²⁶Al from many Galactic novae has not been modelled in this paper.     

Inferring dark halo structure from observed scaling laws of late-type galaxies and LSBs

We re-examine the Fall & Efstathiou scenario for galaxy formation, including the dark halo gravitational reaction to the formation of the baryon disc, as well as continuous variations in the intrinsic halo density profile. The recently published rotation curves of low surface brightness (LSB) and dwarf galaxies together with previously known scaling relations provide sufficient information on the present-day structure of late-type disc galaxies to invert the problem. By requiring that the models reproduce all the observational restrictions we can fully constrain the initial conditions of galaxy formation, with a minimum of assumptions, in particular without the need to specify a cold dark matter (CDM) halo profile. This allows one to solve for all the initial conditions, in terms of the halo density profile, the baryon fraction and the total angular momentum. We find that a unique initial halo shape is sufficient to accurately reproduce the rotation curves of both LSB and normal late-type spiral galaxies. This unique halo profile differs substantially from that found in standard CDM models. A galactic baryon fraction of 0.065 is found. The initial value of the dimensionless angular momentum is seen to be the principal discriminator between the galaxy classes we examine. The present-day scalings between structural parameters are seen to originate in the initial conditions.     

The dark matter distribution in disc galaxies

We use high-quality optical rotation curves of nine low-luminosity disc galaxies to obtain the velocity profiles of the surrounding dark matter haloes. We find that they increase linearly with radius at least out to the edge of the stellar disc, implying that, over the entire stellar region, the density of the dark halo is about constant.
The properties of the mass structure of these haloes are similar to those found for a number of dwarf and low surface brightness galaxies, but provide a more substantial evidence of the discrepancy between the halo mass distribution predicted in the cold dark matter scenario and those actually detected around galaxies. We find that the density law proposed by Burkert reproduces the halo rotation curves, with halo central densities ( ρ ₀∼1–4×10⁻²⁴ g cm⁻³) and core radii ( r ₀∼5–15 kpc) scaling as ρ ₀∝ r ₀^−2/3.