The kinematic footprints of five stellar streams in Andromeda's halo

(abridged) We present a spectroscopic analysis of five stellar streams (`A', `B', `Cr', `Cp' and `D') as well as the extended star cluster, EC4, which lies within streamC, all discovered in the halo of M31 from our CFHT/MegaCam survey. These spectroscopic results were initially serendipitous, making use of our existing observations from the DEep Imaging Multi-Object Spectrograph mounted on the Keck II telescope, and thereby emphasizing the ubiquity of tidal streams that account for ~70% of the M31 halo stars in the targeted fields. Subsequent spectroscopy was then procured in streamCr/p and streamD to trace the velocity gradient along the streams. For the cluster EC4, candidate member stars with average [Fe/H]~-1.4 (Fe/H_spec=-1.6), are found at v_{hel}=-285 km/s suggesting it could be related to streamCp. No similarly obvious cold kinematic candidate is found for streamD, although candidates are proposed in both of two spectroscopic pointings along the stream (both at -400 km/s). Spectroscopy near the edge of streamB suggests a likely kinematic detection, while a candidate kinematic detection of streamA is found (plausibly associated to M33 rather than M31). The low dispersion of the streams in kinematics, physical thickness, and metallicity makes it hard to reconcile with a scenario whereby these stream structures as an ensemble are related to the giant southern stream. We conclude that the M31 stellar halo is largely made up of multiple kinematically cold streams.Comment: 19 pages, 12 figures, accepted in MNRAS. High resolution version, with fig10 here: http://www.ast.cam.ac.uk/~schapman/streams.pd

The role of mergers in driving morphological transformation over cosmic time

Accepted for publication in MNRASUnderstanding the processes that trigger morphological transformation is central to understanding how and why the Universe transitions from being disc-dominated at early epochs to having the morphological mix that is observed today. We use Horizon-AGN, a cosmological hydrodynamical simulation, to perform a comprehensive study of the processes that drive morphological change in massive (M*/M ⊙ > 10 10) galaxies over cosmic time. We show that (1) essentially all the morphological evolution in galaxies that are spheroids at z = 0 is driven by mergers with mass ratios greater than 1: 10; (2) major mergers alone cannot produce today's spheroid population - minor mergers are responsible for a third of all morphological transformation over cosmic time and are its dominant driver after z ~ 1; (3) prograde mergers trigger milder morphological transformation than retrograde mergers - while both types of event produce similar morphological changes at z > 2, the average change due to retrograde mergers is around twice that due to their prograde counterparts at z ~ 0; (4) remnant morphology depends strongly on the gas fraction of a merger, with gas-rich mergers routinely re-growing discs; and (5) at a given stellar mass, discs do not exhibit drastically different merger histories from spheroids - disc survival in mergers is driven by acquisition of cold gas (via cosmological accretion and gas-rich interactions) and a preponderance of prograde mergers in their merger histories.Peer reviewedFinal Accepted Versio

Are Brightest Halo Galaxies Central Galaxies?

It is generally assumed that the central galaxy in a dark matter halo, that is, the galaxy with the lowest specific potential energy, is also the brightest halo galaxy (BHG), and that it resides at rest at the centre of the dark matter potential well. This central galaxy paradigm (CGP) is an essential assumption made in various fields of astronomical research. In this paper we test the validity of the CGP using a large galaxy group catalogue constructed from the Sloan Digital Sky Survey. For each group we compute two statistics, ${\cal R}$ and ${\cal S}$, which quantify the offsets of the line-of-sight velocities and projected positions of brightest group galaxies relative to the other group members. By comparing the cumulative distributions of $|{\cal R}|$ and $|{\cal S}|$ to those obtained from detailed mock group catalogues, we rule out the null-hypothesis that the CGP is correct. Rather, the data indicate that in a non-zero fraction $f_{\rm BNC}(M)$ of all haloes of mass $M$ the BHG is not the central galaxy, but instead, a satellite galaxy. In particular, we find that $f_{\rm BNC}$ increases from $\sim 0.25$ in low mass haloes (10^{12} h^{-1} {\rm M_{\odot}} \leq M \lsim 2 \times 10^{13} h^{-1}{\rm M_{\odot}}) to $\sim 0.4$ in massive haloes (M \gsim 5 \times 10^{13} h^{-1} {\rm M_{\odot}}). We show that these values of $f_{\rm BNC}$ are uncomfortably high compared to predictions from halo occupation statistics and from semi-analytical models of galaxy formation. We end by discussing various implications of a non-zero $f_{\rm BNC}(M)$, with an emphasis on the halo masses inferred from satellite kinematics.Comment: 17 pages, 11 figures. Axes labels of Fig. 3 fixe

Cosmological Galaxy Formation Simulations Using SPH

We present the McMaster Unbiased Galaxy Simulations (MUGS), the first 9 galaxies of an unbiased selection ranging in total mass from 5$\times10^{11}$ M$_\odot$ to 2$\times10^{12}$ M$_\odot$ simulated using n-body smoothed particle hydrodynamics (SPH) at high resolution. The simulations include a treatment of low temperature metal cooling, UV background radiation, star formation, and physically motivated stellar feedback. Mock images of the simulations show that the simulations lie within the observed range of relations such as that between color and magnitude and that between brightness and circular velocity (Tully-Fisher). The greatest discrepancy between the simulated galaxies and observed galaxies is the high concentration of material at the center of the galaxies as represented by the centrally peaked rotation curves and the high bulge-to-total ratios of the simulations determined both kinematically and photometrically. This central concentration represents the excess of low angular momentum material that long has plagued morphological studies of simulated galaxies and suggests that higher resolutions and a more accurate description of feedback will be required to simulate more realistic galaxies. Even with the excess central mass concentrations, the simulations suggest the important role merger history and halo spin play in the formation of disks.Comment: 16 pages, 16 figures, submitted to MNRAS, movies available at http://mugs.mcmaster.ca . Comments welcome

Cuspy No More: How Outflows Affect the Central Dark Matter and Baryon Distribution in Lambda CDM Galaxies

We examine the evolution of the inner dark matter (DM) and baryonic density profile of a new sample of simulated field galaxies using fully cosmological, Lambda CDM, high resolution SPH + N-Body simulations. These simulations include explicit H2 and metal cooling, star formation (SF) and supernovae (SNe) driven gas outflows. Starting at high redshift, rapid, repeated gas outflows following bursty SF transfer energy to the DM component and significantly flatten the originally `cuspy' central DM mass profile of galaxies with present day stellar masses in the 10^4.5 -- 10^9.8 Msolar range. At z=0, the central slope of the DM density profile of our galaxies (measured between 0.3 and 0.7 kpc from their centre) is well fitted by rhoDM propto r^alpha with alpha \simeq -0.5 + 0.35 log_10(Mstar/10^8Msolar) where Mstar is the stellar mass of the galaxy and 4 < log_10 Mstar < 9.4. These values imply DM profiles flatter than those obtained in DM--only simulations and in close agreement with those inferred in galaxies from the THINGS and LITTLE THINGS survey. Only in very small halos, where by z=0 star formation has converted less than ~ 0.03% of the original baryon abundance into stars, outflows do not flatten the original cuspy DM profile out to radii resolved by our simulations. The mass (DM and baryonic) measured within the inner 500 pc of each simulated galaxy remains nearly constant over four orders of magnitudes in stellar mass for Mstar 10^9 Msolar. This finding is consistent with estimates for faint Local Group dwarfs and field galaxies. These results address one of the outstanding problems faced by the CDM model, namely the strong discrepancy between the original predictions of cuspy DM profiles and the shallower central DM distribution observed in galaxies.Comment: MNRAS in press. Accepted version, a few references added. 12 pages. Animation at http://youtu.be/FbcgEovabDI?hd=

Constraints on the steady and pulsed very high energy gamma-ray emission from observations of PSR B1951+32/CTB 80 with the MAGIC Telescope

We report on very high energy gamma-observations with the MAGIC Telescope of the pulsar PSR B1951+32 and its associated nebula, CTB 80. Our data constrain the cutoff energy of the pulsar to be less than 32 GeV, assuming the pulsed gamma-ray emission to be exponentially cut off. The upper limit on the flux of pulsed gamma-ray emission above 75 GeV is 4.3*10^-11 photons cm^-2 sec^-1, and the upper limit on the flux of steady emission above 140 GeV is 1.5*10^-11 photons cm^-2 sec^-1. We discuss our results in the framework of recent model predictions and other studies.Comment: 7 pages, 7 figures, replaced with published versio

First bounds on the very high energy gamma-ray emission from Arp 220

Using the Major Atmospheric Gamma Imaging Cherenkov Telescope (MAGIC), we have observed the nearest ultra-luminous infrared galaxy Arp 220 for about 15 hours. No significant signal was detected within the dedicated amount of observation time. The first upper limits to the very high energy $\gamma$-ray flux of Arp 220 are herein reported and compared with theoretical expectations.Comment: Accepted for publication in Ap

Systematic search for VHE gamma-ray emission from X-ray bright high-frequency BL Lac objects

All but three (M87, BL Lac and 3C 279) extragalactic sources detected so far at very high energy (VHE) gamma-rays belong to the class of high-frequency peaked BL Lac (HBL) objects. This suggested to us a systematic scan of candidate sources with the MAGIC telescope, based on the compilation of X-ray blazars by Donato et al. (2001). The observations took place from December 2004 to March 2006 and cover sources on the northern sky visible under small zenith distances zd < 30 degrees at culmination. The sensitivity of the search was planned for detecting X-ray bright F(1 keV) > 2 uJy) sources emitting at least the same energy flux at 200 GeV as at 1 keV. In order to avoid strong gamma-ray attenuation close to the energy threshold, the redshift of the sources was constrained to values z<0.3. Of the fourteen sources observed, 1ES 1218+304 and 1ES 2344+514 have been detected in addition to the known bright TeV blazars Mrk 421 and Mrk 501. A marginal excess of 3.5 sigma from the position of 1ES 1011+496 was observed and has been confirmed as a source of VHE gamma-rays by a second MAGIC observation triggered by a high optical state (Albert et al. 2007). For the remaining sources, we present here the 99% confidence level upper limits on the integral flux above ~200 GeV. We characterize the sample of HBLs (including all HBLs detected at VHE so far) by looking for correlations between their multi-frequency spectral indices determined from simultaneous optical, archival X-ray, and radio luminosities, finding that the VHE emitting HBLs do not seem to constitute a unique subclass. The absorption corrected gamma-ray luminosities at 200 GeV of the HBLs are generally not higher than their X-ray luminosities at 1 keV.Comment: 15 pages, 7 figures, 5 tables, submitted to ApJ (revised version

Observation of Pulsed Gamma-rays Above 25 GeV from the Crab Pulsar with MAGIC

One fundamental question about pulsars concerns the mechanism of their pulsed electromagnetic emission. Measuring the high-end region of a pulsar's spectrum would shed light on this question. By developing a new electronic trigger, we lowered the threshold of the Major Atmospheric gamma-ray Imaging Cherenkov (MAGIC) telescope to 25 GeV. In this configuration, we detected pulsed gamma-rays from the Crab pulsar that were greater than 25 GeV, revealing a relatively high cutoff energy in the phase-averaged spectrum. This indicates that the emission occurs far out in the magnetosphere, hence excluding the polar-cap scenario as a possible explanation of our measurement. The high cutoff energy also challenges the slot-gap scenario.Comment: Slight modification of the analysis: Fitting a more general function to the combined data set of COMPTEL, EGRET and MAGIC. Final result and conclusion is unchange

Upper limit for gamma-ray emission above 140 GeV from the dwarf spheroidal galaxy Draco

The nearby dwarf spheroidal galaxy Draco with its high mass to light ratio is one of the most auspicious targets for indirect dark matter searches. Annihilation of hypothetical DM particles can result in high-energy gamma-rays, e.g. from neutralino annihilation in the supersymmetric framework. With the MAGIC telescope a search for a possible DM signal originating from Draco was performed during 2007. The analysis of the data results in a flux upper limit of 1.1x10^-11 photons cm^-2 sec^-1 for photon energies above 140 GeV, assuming a point like source. Furthermore, a comparison with predictions from supersymmetric models is given. While our results do not constrain the mSUGRA phase parameter space, a very high flux enhancement can be ruled out.Comment: Accepted for publication by Astrophysical Journa