Showing papers by "Klaus Dolag published in 2013"

Planck 2013 results. XVI. Cosmological parameters

Abstract: We present the first results based on Planck measurements of the CMB temperature and lensing-potential power spectra. The Planck spectra at high multipoles are extremely well described by the standard spatially-flat six-parameter LCDM cosmology. In this model Planck data determine the cosmological parameters to high precision. We find a low value of the Hubble constant, H0=67.3+/-1.2 km/s/Mpc and a high value of the matter density parameter, Omega_m=0.315+/-0.017 (+/-1 sigma errors) in excellent agreement with constraints from baryon acoustic oscillation (BAO) surveys. Including curvature, we find that the Universe is consistent with spatial flatness to percent-level precision using Planck CMB data alone. We present results from an analysis of extensions to the standard cosmology, using astrophysical data sets in addition to Planck and high-resolution CMB data. None of these models are favoured significantly over standard LCDM. The deviation of the scalar spectral index from unity is insensitive to the addition of tensor modes and to changes in the matter content of the Universe. We find a 95% upper limit of r<0.11 on the tensor-to-scalar ratio. There is no evidence for additional neutrino-like relativistic particles. Using BAO and CMB data, we find N_eff=3.30+/-0.27 for the effective number of relativistic degrees of freedom, and an upper limit of 0.23 eV for the summed neutrino mass. Our results are in excellent agreement with big bang nucleosynthesis and the standard value of N_eff=3.046. We find no evidence for dynamical dark energy. Despite the success of the standard LCDM model, this cosmology does not provide a good fit to the CMB power spectrum at low multipoles, as noted previously by the WMAP team. While not of decisive significance, this is an anomaly in an otherwise self-consistent analysis of the Planck temperature data.

The Hot and Energetic Universe: A White Paper presenting the science theme motivating the Athena+ mission

Abstract: This White Paper, submitted to the recent ESA call for science themes to define its future large missions, advocates the need for a transformational leap in our understanding of two key questions in astrophysics: 1) How does ordinary matter assemble into the large scale structures that we see today? 2) How do black holes grow and shape the Universe? Hot gas in clusters, groups and the intergalactic medium dominates the baryonic content of the local Universe. To understand the astrophysical processes responsible for the formation and assembly of these large structures, it is necessary to measure their physical properties and evolution. This requires spatially resolved X-ray spectroscopy with a factor 10 increase in both telescope throughput and spatial resolving power compared to currently planned facilities. Feedback from supermassive black holes is an essential ingredient in this process and in most galaxy evolution models, but it is not well understood. X-ray observations can uniquely reveal the mechanisms launching winds close to black holes and determine the coupling of the energy and matter flows on larger scales. Due to the effects of feedback, a complete understanding of galaxy evolution requires knowledge of the obscured growth of supermassive black holes through cosmic time, out to the redshifts where the first galaxies form. X-ray emission is the most reliable way to reveal accreting black holes, but deep survey speed must improve by a factor ~100 over current facilities to perform a full census into the early Universe. The Advanced Telescope for High Energy Astrophysics (Athena+) mission provides the necessary performance (e.g. angular resolution, spectral resolution, survey grasp) to address these questions and revolutionize our understanding of the Hot and Energetic Universe. These capabilities will also provide a powerful observatory to be used in all areas of astrophysics.

Planck 2013 results. XXIX. Planck catalogue of Sunyaev-Zeldovich sources

Abstract: We describe the all-sky Planck catalogue of clusters and cluster candidates derived from Sunyaev--Zeldovich (SZ) effect detections using the first 15.5 months of Planck satellite observations. The catalogue contains 1227 entries, making it over six times the size of the Planck Early SZ (ESZ) sample and the largest SZ-selected catalogue to date. It contains 861 confirmed clusters, of which 178 have been confirmed as clusters, mostly through follow-up observations, and a further 683 are previously-known clusters. The remaining 366 have the status of cluster candidates, and we divide them into three classes according to the quality of evidence that they are likely to be true clusters. The Planck SZ catalogue is the deepest all-sky cluster catalogue, with redshifts up to about one, and spans the broadest cluster mass range from (0.1 to 1.6) 10^{15}Msun. Confirmation of cluster candidates through comparison with existing surveys or cluster catalogues is extensively described, as is the statistical characterization of the catalogue in terms of completeness and statistical reliability. The outputs of the validation process are provided as additional information. This gives, in particular, an ensemble of 813 cluster redshifts, and for all these Planck clusters we also include a mass estimated from a newly-proposed SZ-mass proxy. A refined measure of the SZ Compton parameter for the clusters with X-ray counter-parts is provided, as is an X-ray flux for all the Planck clusters not previously detected in X-ray surveys.

Planck 2013 results. I. Overview of products and scientific results

Abstract: The ESA's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14 May 2009 and has been scanning the microwave and submillimetre sky continuously since 12 August 2009. This paper gives an overview of the mission and its performance, the processing, analysis, and characteristics of the data, the scientific results, and the science data products and papers in the release. The science products include maps of the CMB and diffuse extragalactic foregrounds, a catalogue of compact Galactic and extragalactic sources, and a list of sources detected through the SZ effect. The likelihood code used to assess cosmological models against the Planck data and a lensing likelihood are described. Scientific results include robust support for the standard six-parameter LCDM model of cosmology and improved measurements of its parameters, including a highly significant deviation from scale invariance of the primordial power spectrum. The Planck values for these parameters and others derived from them are significantly different from those previously determined. Several large-scale anomalies in the temperature distribution of the CMB, first detected by WMAP, are confirmed with higher confidence. Planck sets new limits on the number and mass of neutrinos, and has measured gravitational lensing of CMB anisotropies at greater than 25 sigma. Planck finds no evidence for non-Gaussianity in the CMB. Planck's results agree well with results from the measurements of baryon acoustic oscillations. Planck finds a lower Hubble constant than found in some more local measures. Some tension is also present between the amplitude of matter fluctuations derived from CMB data and that derived from SZ data. The Planck and WMAP power spectra are offset from each other by an average level of about 2% around the first acoustic peak.

Planck intermediate results: V. Pressure profiles of galaxy clusters from the Sunyaev-Zeldovich effect

Abstract: Taking advantage of the all-sky coverage and broad frequency range of the Planck satellite, we study the Sunyaev-Zeldovich (SZ) and pressure profiles of 62 nearby massive clusters detected at high significance in the 14-month nominal survey. Careful reconstruction of the SZ signal indicates that most clusters are individually detected at least out to R500. By stacking the radial profiles, we have statistically detected the radial SZ signal out to 3R500, i.e., at a density contrast of about 50-100, though the dispersion about the mean profile dominates the statistical errors across the whole radial range. Our measurement is fully consistent with previous Planck results on integrated SZ fluxes, further strengthening the agreement between SZ and X-ray measurements inside R500. Correcting for the effects of the Planck beam, we have calculated the corresponding pressure profiles. This new constraint from SZ measurements is consistent with the X-ray constraints from xmm in the region in which the profiles overlap (i.e., [0.1-1] R500), and is in fairly good agreement with theoretical predictions within the expected dispersion. At larger radii the average pressure profile is shallower than the predictions. Combining the SZ and X-ray observed profiles into a joint fit to a generalised pressure profile gives best-fit parameters [P0, c500, gamma, alpha, beta] = [6.41, 1.81, 0.31, 1.33, 4.13]. Using a reasonable hypothesis for the gas temperature in the cluster outskirts we reconstruct from our stacked pressure profile the gas mass fraction profile out to 3R500. Within the temperature driven uncertainties, our Planck constraints are compatible with the cosmic baryon fraction and expected gas fraction in halos.

The pre-launch Planck Sky Model: a model of sky emission at submillimetre to centimetre wavelengths

Abstract: We present the Planck Sky Model (PSM), a parametric model for the generation of all-sky, few arcminute resolution maps of sky emission at submillimetre to centimetre wavelengths, in both intensity and polarisation. Several options are implemented to model the cosmic microwave background, Galactic diffuse emission (synchrotron, free-free, thermal and spinning dust, CO lines), Galactic H-II regions, extragalactic radio sources, dusty galaxies, and thermal and kinetic Sunyaev-Zeldovich signals from clusters of galaxies. Each component is simulated by means of educated interpolations/extrapolations of data sets available at the time of the launch of the Planck mission, complemented by state-of-the-art models of the emission. Distinctive features of the simulations are: spatially varying spectral properties of synchrotron and dust; different spectral parameters for each point source; modeling of the clustering properties of extragalactic sources and of the power spectrum of fluctuations in the cosmic infrared background. The PSM enables the production of random realizations of the sky emission, constrained to match observational data within their uncertainties, and is implemented in a software package that is regularly updated with incoming information from observations. The model is expected to serve as a useful tool for optimizing planned microwave and sub-millimetre surveys and to test data processing and analysis pipelines. It is, in particular, used for the development and validation of data analysis pipelines within the planck collaboration. A version of the software that can be used for simulating the observations for a variety of experiments is made available on a dedicated website.

Structure finding in cosmological simulations: the state of affairs

Abstract: The ever increasing size and complexity of data coming from simulations of cosmic structure formation demand equally sophisticated tools for their analysis. During the past decade, the art of object finding in these simulations has hence developed into an important discipline itself. A multitude of codes based upon a huge variety of methods and techniques have been spawned yet the question remained as to whether or not they will provide the same (physical) information about the structures of interest. Here we summarize and extent previous work of the 'halo finder comparison project': we investigate in detail the (possible) origin of any deviations across finders. To this extent, we decipher and discuss differences in halo-finding methods, clearly separating them from the disparity in definitions of halo properties. We observe that different codes not only find different numbers of objects leading to a scatter of up to 20 per cent in the halo mass and V-max function, but also that the particulars of those objects that are identified by all finders differ. The strength of the variation, however, depends on the property studied, e.g. the scatter in position, bulk velocity, mass and the peak value of the rotation curve is practically below a few per cent, whereas derived quantities such as spin and shape show larger deviations. Our study indicates that the prime contribution to differences in halo properties across codes stems from the distinct particle collection methods and - to a minor extent - the particular aspects of how the procedure for removing unbound particles is implemented. We close with a discussion of the relevance and implications of the scatter across different codes for other fields such as semi-analytical galaxy formation models, gravitational lensing and observables in general.

Revisiting scaling relations for giant radio halos in galaxy clusters

Abstract: Many galaxy clusters host megaparsec-scale radio halos, generated by ultrarelativistic electrons in the magnetized intracluster medium. Correlations between the synchrotron power of radio halos and the thermal properties of the hosting clusters were established in the last decade, including the connection between the presence of a halo and cluster mergers. The X-ray luminosity and redshift-limited Extended GMRT Radio Halo Survey provides a rich and unique dataset for statistical studies of the halos. We uniformly analyze the radio and X-ray data for the GMRT cluster sample, and use the new Planck Sunyaev-Zel'dovich (SZ) catalog to revisit the correlations between the power of radio halos and the thermal properties of galaxy clusters. We find that the radio power at 1.4 GHz scales with the cluster X-ray (0.1-2.4 keV) luminosity computed within R(sub 500) as P(sub 1.4) approx. L(2.1+/-0.2) - 500). Our bigger and more homogenous sample confirms that the X-ray luminous (L(sub 500) > 5 × 10(exp 44) erg/s)) clusters branch into two populations-radio halos lie on the correlation, while clusters without radio halos have their radio upper limits well below that correlation. This bimodality remains if we excise cool cores from the X-ray luminosities. We also find that P(sub 1.4) scales with the cluster integrated SZ signal within R(sub 500), measured by Planck, as P(sub 1.4) approx. Y(2.05+/-0.28) - 500), in line with previous findings. However, contrary to previous studies that were limited by incompleteness and small sample size, we find that "SZ-luminous" Y(sub 500) > 6×10(exp −5) Mpc(exp 2) clusters show a bimodal behavior for the presence of radio halos, similar to that in the radio-X-ray diagram. Bimodality of both correlations can be traced to clusters dynamics, with radio halos found exclusively in merging clusters. These results confirm the key role of mergers for the origin of giant radio halos, suggesting that they trigger the relativistic particle acceleration.

Planck 2013 results. XXI. All-sky Compton parameter power spectrum and high-order statistics

Abstract: We have constructed the first all-sky map of the thermal Sunyaev-Zeldovich (tSZ) effect by applying specifically tailored component separation algorithms to the 100 to 857 GHz frequency channel maps from the Planck survey. These maps show an obvious galaxy cluster tSZ signal that is well matched with blindly detected clusters in the Planck SZ catalogue. To characterize the signal in the tSZ map we have computed its angular power spectrum. At large angular scales ($\ell 500$) the clustered Cosmic Infrared Background (CIB) and residual point sources are the major contaminants. These foregrounds are carefully modelled and subtracted. We measure the tSZ power spectrum in angular scales, $0.17^{\circ} \lesssim \theta \lesssim 3.0^{\circ}$, that were previously unexplored. The measured tSZ power spectrum is consistent with that expected from the Planck catalogue of SZ sources, with additional clear evidence of signal from unresolved clusters and, potentially, diffuse warm baryons. We use the tSZ power spectrum to obtain the following cosmological constraints: $\sigma_8(\Omega_{\mathrm{m}}/0.28)^{3.2/8.1}=0.784 \pm 0.016 (68% C.L.). Marginalized band-powers of the Planck tSZ power spectrum and the best-fit model are given. The non-Gaussianity of the Compton parameter map is further characterized by computing its 1D probability distribution function and its bispectrum. These are used to place additional independent constraints on $\sigma_{8}$.

Planck 2013 results. XIX. The integrated Sachs-Wolfe effect

Abstract: Based on CMB maps from the 2013 Planck Mission data release, this paper presents the detection of the ISW effect, i.e., the correlation between the CMB and large-scale evolving gravitational potentials. The significance of detection ranges from 2 to 4 sigma, depending on which method is used. We investigate three separate approaches, which cover essentially all previous studies, as well as breaking new ground. (i) Correlation of the CMB with the Planck reconstructed gravitational lensing potential (for the first time). This detection is made using the lensing-induced bispectrum; the correlation between lensing and the ISW effect has a significance close to 2.5 sigma. (ii) Cross-correlation with tracers of LSS, yielding around 3 sigma significance, based on a combination of radio (NVSS) and optical (SDSS) data. (iii) Aperture photometry on stacked CMB fields at the locations of known large-scale structures, which yields a 4 sigma signal when using a previously explored catalogue, but shows strong discrepancies in amplitude and scale compared to expectations. More recent catalogues give more moderate results, ranging from negligible to 2.5 sigma at most, but with a more consistent scale and amplitude, the latter being still slightly above what is expected from numerical simulations within LCMD. Where they can be compared, these measurements are compatible with previous work using data from WMAP, which had already mapped these scales to the limits of cosmic variance. Planck's broader frequency coverage confirms that the signal is achromatic, bolstering the case for ISW detection. As a final step we use tracers of large-scale structure to filter the CMB data, presenting maps of the ISW temperature perturbation. These results provide complementary and independent evidence for the existence of a dark energy component that governs the current accelerated expansion of the Universe.

Baryon Census in Hydrodynamical Simulations of Galaxy Clusters

Abstract: We carry out an analysis of a set of cosmological SPH hydrodynamical simulations of galaxy clusters and groups aimed at studying the total baryon budget in clusters, and how this budget is shared between the hot diffuse component and the stellar component. Using the TreePM+SPHGADGET-3 code, we carried out one set of non‐radiative simulations, and two sets of simulations including radiative cooling, star form ation and feedback from supernovae (SN), one of which also accounting for the effect of feedback from active galactic nuclei (AGN). The analysis is carried out with the twofold aim of studying the implication of stellar and hot gas content on the relative role played by SN and AGN feedback, and to calibrate the cluster baryon fraction and its evolution as a cosmological tool. With respect to previous similar analysis, the simulations used in this study provid e us with a sufficient statistics of massive objects and including an efficient AGN feedback. We fi nd that both radiative simulation sets predict a trend of stellar mass fraction with clu ster mass that tends to be weaker than the observed one. However this tension depends on the particular set of observational data considered. Including the effect of AGN feedback alleviates this tension on the stellar mass and predicts values of the hot gas mass fraction and total baryon fraction to be in closer agreement with observational results. We further compute the ratio between the cluster baryon content and the cosmic baryon fraction, Yb, as a function of cluster-centric radius and redshift. At R500 we find for massive clusters with M500 > 2×10 14 h −1 M⊙ that Yb is nearly independent of the physical processes included and characterized by a negligible redshift evolution: Yb,500 = 0.85 ± 0.03 with the error accounting for the intrinsic r.m.s. scatter w ithin the set of simulated clusters. At smaller radii, R2500, the typical value of Yb slightly decreases, by an amount that depends on the physics included in the simulations, while its scatter increases by about a factor of two. These results have interesting implications for the cosmological applications of the baryon fraction in clusters.

Toward Unbiased Galaxy Cluster Masses from Line-of-sight Velocity Dispersions

Abstract: We study the use of red-sequence-selected galaxy spectroscopy for unbiased estimation of galaxy cluster masses by using a publicly available simulated galaxy catalog. We explore the impact of selection using galaxy color, projected separation from the cluster center, galaxy luminosity, and spectroscopic redshift. We identify and characterize each of the following sources of bias and scatter in velocity dispersion at fixed mass: the intrinsic properties of halos in the form of halo triaxiality, sampling noise, the presence of multiple kinematic populations within the cluster, and the effect of interlopers. We show that even in red-sequence and spectroscopically selected galaxy samples, the interloper fraction is significant, and that the variations in the interloper population from cluster to cluster provide the dominant contribution to the velocity dispersion scatter at fixed mass. We present measurements of the total scatter in dispersion at fixed mass as a function of the number of redshifts. Results indicate that improvements in scatter are modest beyond samples of ~30 redshifts per cluster. Our results show that while cluster velocity dispersions extracted from a few dozen red-sequence-selected galaxies do not provide precise masses on a single cluster basis, an ensemble of cluster velocity dispersions can be combined to produce a precise calibration of a cluster survey-mass-observable relation. Currently, disagreements in the literature on simulated subhalo velocity dispersion-mass relations place a systematic floor on velocity dispersion mass calibration at the 5% level in dispersion.

Rise and fall of radio haloes in simulated merging galaxy clusters

Abstract: We present the first high resolution MHD simulation of cosmic-ray electron reacceleration by turbulence in cluster mergers. We use an idealised model for cluster mergers, combined with a numerical model for the injection, cooling and reacceleration of cosmic-ray electrons, to investigate the evolution of cluster scale radio emission in these objects. In line with theoretical expectations, we for the first time, show in a simulation that reacceleration of CRe has the potential to reproduce key observables of radio halos. In particular, we show that clusters evolve being radio loud or radio quiet, depending on their evolutionary stage during the merger. We thus recover the observed transient nature of radio halos. In the simulation the diffuse emission traces the complex interplay between spatial distribution of turbulence injected by the halo infall and the spatial distribution of the seed electrons to reaccelerate. During the formation and evolution of the halo the synchrotron emission spectra show the observed variety: from power-laws with spectral index of 1 to 1.3 to curved and ultra-steep spectra with index > 1:5.

Planck intermediate results - X. Physics of the hot gas in the Coma cluster

Abstract: We present an analysis of Planck satellite data on the Coma cluster observed via the Sunyaev-Zeldovich effect. Thanks to its great sensitivity, Planck is able, for the first time, to detect SZ emission up to r ≈ 3 × R500. We test previously proposed spherically symmetric models for the pressure distribution in clusters against the azimuthally averaged data. In particular, we find that the Arnaud et al. (2010, AA by r = 2 × R500 it underestimates the observed y profile by a factor of ≃2. This may indicate that at these larger radii either: i) the cluster SZ emission is contaminated by unresolved SZ sources along the line of sight; or ii) the pressure profile of Coma is higher at r > R500 than the mean pressure profile predicted by the simulations used to constrain the models. The Planck image shows significant local steepening of the y profile in two regions about half a degree to the west and to the south-east of the cluster centre. These features are consistent with the presence of shock fronts at these radii, and indeed the western feature was previously noticed in the ROSAT PSPC mosaic as well as in the radio. Using Plancky profiles extracted from corresponding sectors we find pressure jumps of 4.9-0.2+0.4 and 5.0-0.1+1.3 in the west and south-east, respectively. Assuming Rankine-Hugoniot pressure jump conditions, we deduce that the shock waves should propagate with Mach number Mw = 2.03-0.04+0.09 and Mse = 2.05-0.02+0.25 in the west and south-east, respectively. Finally, we find that the y and radio-synchrotron signals are quasi-linearly correlated on Mpc scales, with small intrinsic scatter. This implies either that the energy density of cosmic-ray electrons is relatively constant throughout the cluster, or that the magnetic fields fall off much more slowly with radius than previously thought.

Planck intermediate results - VIII. Filaments between interacting clusters

Abstract: Context. About half of the baryons of the Universe are expected to be in the form of filaments of hot and low-density intergalactic medium. Most of these baryons remain undetected even by the most advanced X-ray observatories, which are limited in sensitivity to the diffuse low-density medium. Aims. The Planck satellite has provided hundreds of detections of the hot gas in clusters of galaxies via the thermal Sunyaev-Zel'dovich (tSZ) effect and is an ideal instrument for studying extended low-density media through the tSZ effect. In this paper we use the Planck data to search for signatures of a fraction of these missing baryons between pairs of galaxy clusters. Methods. Cluster pairs are good candidates for searching for the hotter and denser phase of the intergalactic medium (which is more easily observed through the SZ effect). Using an X-ray catalogue of clusters and the Planck data, we selected physical pairs of clusters as candidates. Using the Planck data, we constructed a local map of the tSZ effect centred on each pair of galaxy clusters. ROSAT data were used to construct X-ray maps of these pairs. After modelling and subtracting the tSZ effect and X-ray emission for each cluster in the pair, we studied the residuals on both the SZ and X-ray maps. Results. For the merging cluster pair A399-A401 we observe a significant tSZ effect signal in the intercluster region beyond the virial radii of the clusters. A joint X-ray SZ analysis allows us to constrain the temperature and density of this intercluster medium. We obtain a temperature of kT = 7.1 +/- 0.9 keV (consistent with previous estimates) and a baryon density of (3.7 +/- 0.2) x 10(-4) cm(-3). Conclusions. The Planck satellite mission has provided the first SZ detection of the hot and diffuse intercluster gas.

Planck intermediate results. XIII. Constraints on peculiar velocities

Abstract: Using \Planck\ data combined with the Meta Catalogue of X-ray detected Clusters of galaxies (MCXC), we address the study of peculiar motions by searching for evidence of the kinetic Sunyaev-Zeldovich effect (kSZ). By implementing various filters designed to extract the kSZ generated at the positions of the clusters, we obtain consistent constraints on the radial peculiar velocity average, root mean square (rms), and local bulk flow amplitude at different depths. For the whole cluster sample of average redshift 0.18, the measured average radial peculiar velocity with respect to the cosmic microwave background (CMB) radiation at that redshift, i.e., the kSZ monopole, amounts to $72 \pm 60$ km s$^{-1}$. This constitutes less than 1% of the relative Hubble velocity of the cluster sample with respect to our local CMB frame. While the linear $\Lambda$CDM prediction for the typical cluster radial velocity rms at $z=0.15$ is close to 230km s$^{-1}$, the upper limit imposed by \Planck\ data on the cluster subsample corresponds to 800 km s$^{-1}$ at 95% confidence level, i.e., about three times higher. \Planck\ data also set strong constraints on the local bulk flow in volumes centred on the Local Group. There is no detection of bulk flow as measured in any comoving sphere extending to the maximum redshift covered by the cluster sample. A blind search for bulk flows in this sample has an upper limit of 254 km s$^{-1}$ (95% confidence level) dominated by CMB confusion and instrumental noise, indicating that the Universe is largely homogeneous on Gpc scales. In this context, in conjunction with supernova observations, \Planck\ is able to rule out a large class of inhomogeneous void models as alternatives to dark energy or modified gravity. The \Planck\ constraints on peculiar velocities and bulk flows are thus consistent with the $\Lambda$CDM scenario.

The dark halo—spheroid conspiracy and the origin of elliptical galaxies

Abstract: Dynamical modeling and strong-lensing data indicate that the total density profiles of early-type galaxies are close to isothermal, i.e., {rho}{sub tot}{proportional_to}r {sup {gamma}} with {gamma} Almost-Equal-To -2. To understand the origin of this universal slope we study a set of simulated spheroids formed in isolated binary mergers as well as the formation within the cosmological framework. The total stellar plus dark matter density profiles can always be described by a power law with an index of {gamma} Almost-Equal-To -2.1 with a tendency toward steeper slopes for more compact, lower-mass ellipticals. In the binary mergers the amount of gas involved in the merger determines the precise steepness of the slope. This agrees with results from the cosmological simulations where ellipticals with steeper slopes have a higher fraction of stars formed in situ. Each gas-poor merger event evolves the slope toward {gamma} {approx} -2, once this slope is reached further merger events do not change it anymore. All our ellipticals have flat intrinsic combined stellar and dark matter velocity dispersion profiles. We conclude that flat velocity dispersion profiles and total density distributions with a slope of {gamma} {approx} -2 for the combined system of stars and dark matter act as a naturalmore » attractor. The variety of complex formation histories as present in cosmological simulations, including major as well as minor merger events, is essential to generate the full range of observed density slopes seen for present-day elliptical galaxies.« less

Sussing merger trees: the Merger Trees Comparison Project

Abstract: Merger trees follow the growth and merger of dark-matter haloes over cosmic history. As well as giving important insights into the growth of cosmic structure in their own right, they provide an essential backbone to semi-analytic models of galaxy formation. This paper is the first in a series to arise from the Sussing Merger Trees Workshop in which 10 different tree-building algorithms were applied to the same set of halo catalogues and their results compared. Although many of these codes were similar in nature, all algorithms produced distinct results. Our main conclusions are that a useful merger-tree code should possess the following features: (i) the use of particle IDs to match haloes between snapshots; (ii) the ability to skip at least one, and preferably more, snapshots in order to recover subhaloes that are temporarily lost during merging; (iii) the ability to cope with (and ideally smooth out) large, temporary fluctuations in halo mass. Finally, to enable different groups to communicate effectively, we defined a common terminology that we used when discussing merger trees and we encourage others to adopt the same language. We also specified a minimal output format to record the results.

Large-scale inhomogeneities of the intracluster medium: improving mass estimates using the observed azimuthal scatter

Abstract: Achieving a robust determination of the gas density profile i n cluster outskirts is a crucial point in order to measure their baryonic content and to use them as cosmological probes. The difficulty in obtaining this measurement lies not only in the low surface brightness of the ICM, but also in the inhomogeneities of the gas associated to clumps, asymmetries and accretion patterns. Using a set of hydrodynamical simulations of 62 galaxy clusters and groups we study this kind of inhomogeneities, focusing on the ones on the large scale that, unlike clumps, are the most difficult to identify. To this purpose we introduce the concept of residual clumpiness, CR, that quantifies the large-scale inhomogeneity of the ICM. A fter showing that this quantity can be robustly defined for relaxed systems, we characterize how it varies with radius, mass and dynamical state of the halo. Most importantly, we observe that it introduces an overestimate in the determination of the density profile f rom the X-ray emission, which translates into a systematic overestimate of 6 (12) per cent in the measurement of Mgas at R200 for our relaxed (perturbed) cluster sample. At the same time, the increase of CR with radius introduces also a �2 per cent systematic underestimate in the measurement of the hydrostatic-equilibrium mass (Mhe), which adds to the previous one generating a systematic �8.5 per cent overestimate in fgas in our relaxed sample. Since the residual clumpiness of the ICM is not directly observable, we study its correlation with the azimuthal scatter in the X-ray surface brightness of the halo, a quantity that is well -constrained by current measurements, and in the y-parameter profiles that is at reach of the forthcoming SZ exp eriments. We find that their correlation is highly significant ( rS = 0.6 0.7), allowing to define the azimuthal scatter measured in the X-ray surface brightness profile and in the y-parameter as robust proxies of CR. After providing a function that connects the two quantitie s, we obtain that correcting the observed gas density profiles using the a zimuthal scatter eliminates the bias in the measurement of Mgas for relaxed objects, which becomes 0 ± 2 per cent up to 2R200, and reduces it by a factor of 3 for perturbed ones. This method allows also to eliminate the systematics on the measurements of Mhe and fgas, although a significant halo to halo scatter remains.

Investigating the velocity structure and X-ray observable properties of simulated galaxy clusters with PHOX

Abstract: Non-thermal motions in the intra-cluster medium (ICM) are believed to play a non-negligible role in the pressure support to the total gravitating mass of galaxy clusters. Future X-ray missions, such as ASTRO-H and ATHENA, will eventually allow us to directly detect the signature of these motions from high-resolution spectra of the ICM. In this paper, we present a study on a set of clusters extracted from a cosmological hydrodynamical simulation, devoted to explore the role of non-thermal velocity amplitude in characterising the cluster state and the relation between observed X-ray properties. In order to reach this goal, we apply the X-ray virtual telescope PHOX to generate synthetic observations of the simulated clusters with both Chandra and ATHENA, the latter used as an example for the performance of very high-resolution X-ray telescopes. From Chandra spectra we extract global properties, e.g. luminosity and temperature, and from ATHENA spectra we estimate the gas velocity dispersion along the line of sight from the broadening of heavy-ion emission lines (e.g. Fe). We further extend the analysis to the relation between non-thermal velocity dispersion of the gas and the L_X-T scaling law for the simulated clusters. Interestingly, we find a clear dependence of slope and scatter on the selection criterion for the clusters, based on the level of significance of non-thermal motions. Namely, the scatter in the relation is significantly reduced by the exclusion of the clusters, for which we estimate the highest turbulent velocities. Such velocity diagnostics appears therefore as a promising independent way to identify disturbed clusters, in addition to the commonly used morphological inspection.

Planck intermediate results - III. The relation between galaxy cluster mass and Sunyaev-Zeldovich signal

Abstract: We examine the relation between the galaxy cluster mass M and Sunyaev-Zeldovich ({SZ)} effect signal {DA2} Y500 for a sample of 19 objects for which weak lensing ({WL)} mass measurements obtained from Subaru Telescope data are available in the literature. Hydrostatic X-ray masses are derived from {XMM-Newton} archive data, and the {SZ} effect signal is measured from Planck all-sky survey data. We find an {MWL} - {DA2} Y500 relation that is consistent in slope and normalisation with previous determinations using weak lensing masses; however, there is a normalisation offset with respect to previous measures based on hydrostatic X-ray mass-proxy relations. We verify that our {SZ} effect measurements are in excellent agreement with previous determinations from Planck data. For the present sample, the hydrostatic X-ray masses at R500 are on average {\textasciitilde} 20 percent larger than the corresponding weak lensing masses, which is contrary to expectations. We show that the mass discrepancy is driven by a difference in mass concentration as measured by the two methods and, for the present sample, that the mass discrepancy and difference in mass concentration are especially large for disturbed systems. The mass discrepancy is also linked to the offset in centres used by the X-ray and weak lensing analyses, which again is most important in disturbed systems. We outline several approaches that are needed to help achieve convergence in cluster mass measurement with X-ray and weak lensing observations. Appendices are available in electronic form at http://www.aanda.org

Galactic winds in cosmological simulations of the circumgalactic medium

Abstract: We explore new observationally-constrained sub-resolution models of galactic outflows and investigate their impact on the circumgalactic medium (CGM) in the redshift range z = 2 4. We perform cosmological hydrodynamic simulations, including star formation, chemical enrichment, and four cases of SNe-driven outflows: no wind (NW), an energy-driven constant velocity wind (CW), a radially varying wind (RVWa) where the outflow velocity has a positive correlation with galactocentric distance (r), and a RVW with additional dependence on halo mass (RVWb). Overall, we find that the outflows expel metal-enriched gas away from galaxies, significantly quench the star formation, reduce the central galactic metallicity and enrich the CGM. At z = 2, the radial profiles of gas properties around galaxy centers are most sensitive to the choice of the wind model for halo masses in the range (10 9 10 11 )M⊙. We infer that outflows in the RVWb model are least effective, with results similar to the NW case, except that the CGM is enriched more. Moreover, we find that the models CW and RVWa are similar, both showing the impact of effective winds, with the following notable differences. RVWa causes a greater suppression of star formation rate at z 6 5, and has a higher fraction of low-density (δ < 10), warm-hot (10 4 10 6 K) gas than in CW. Outflows in CW produce a higher and earlier enrichment of some IGM phases than in RVWa. By visual inspection, we note that the RVWa model shows galactic disks more pronounced than all the other wind models. We predict that some observational diagnostics are more promising to distinguish between different outflow driving mechanisms in galaxies: ZC of the CGM gas at r � (30 300)h −1 kpc comoving, and CIV fraction of the inner gas at r < (4 5)h −1 kpc comoving.

The Hot and Energetic Universe: The missing baryons and the warm-hot intergalactic medium

Abstract: The backbone of the large-scale structure of the Universe is determined by processes on a cosmological scale and by the gravitational interaction of the dominant dark matter. However, the mobile baryon population shapes the appearance of these structures. Theory predicts that most of the baryons reside in vast unvirialized filamentary structures that connect galaxy groups and clusters, but the observational evidence is currently lacking. Because the majority of the baryons are supposed to exist in a large-scale, hot and dilute gaseous phase, X-rays provide the ideal tool to progress our understanding. Observations with the Athena+ X-ray Integral Field Unit will reveal the location, chemical composition, physical state and dynamics of the active population of baryons.

Virial scaling of galaxies in clusters: bright to faint is cool to hot

Abstract: By combining galaxy tracers from high-resolution N-body and hydrodynamical simulations, we present a consistent picture of the behaviour of galaxy velocities in massive clusters. In haloes above ~ 10^14 Msun, the brightest satellite galaxies are slightly cooler compared to the dark matter, while fainter satellites are hotter. Within the virial radius of a cluster, the mean velocity dispersion based on the 100 brightest galaxies is a factor of 1.065 +/- 0.005 (stat) +/- 0.027 (sys) higher than that of the dark matter (corresponding to a ~10-15 per cent bias in the dynamical mass estimate) while that based on only the five brightest galaxies is 0.868 +/- 0.039 (stat) +/- 0.035 (sys). These trends are approximately independent of redshift. The velocity structure is sensitive to the modelling of galaxies in clusters, indicative of the complex interplay of tidal stripping, dynamical friction, and merging. Velocity dispersions derived from instantaneous subhalo properties are larger than those employing either peak subhalo properties or hydrodynamical galaxy tracers. The latter two methods are consistent, implying that stacked spectroscopic analysis of cluster samples should, after correction for projection, show a trend towards slightly higher velocities when fainter galaxies are included, with an unbiased measure of dark matter velocity dispersion coming from approximately 30 galaxies per cluster. We show evidence that the velocity distribution function of bright galaxies near the cluster centre has a low-velocity tail due to strong dynamical friction.

Numerical simulations of the possible origin of the two sub-parsec scale and counterrotating stellar disks around SgrA*

Abstract: We present a high-resolution simulation of an idealized model to explain the origin of the two young, counterrotating, sub-parsec scale stellar disks around the supermassive black hole SgrA* at the center of the Milky Way In our model, the collision of a single molecular cloud with a circumnuclear gas disk (similar to the one observed presently) leads to multiple streams of gas flowing toward the black hole and creating accretion disks with angular momentum depending on the ratio of cloud and circumnuclear disk material The infalling gas creates two inclined, counterrotating sub-parsec scale accretion disks around the supermassive black hole with the first disk forming roughly 1 Myr earlier, allowing it to fragment into stars and get dispersed before the second counterrotating disk forms Fragmentation of the second disk would lead to the two inclined, counterrotating stellar disks which are observed at the Galactic center A similar event might be happening again right now at the Milky Way Galactic center Our model predicts that the collision event generates spiral-like filaments of gas, feeding the Galactic center prior to disk formation with a geometry and inflow pattern that is in agreement with the structure of the so-called mini spiral that has been detected in the Galactic center

Hydrodynamical simulations of galaxy clusters in dark energy cosmologies – II. c–M relation

Abstract: We use hydrodynamical simulations of different dark energy cosmologies to investigate the concentration-mass (c M ) relation in galaxy clusters. In particular, we consider a r eference �CDM model, two quintessence models with inverse power-law potentials (RP and SUGRA), and two extended quintessence models, one with positive and one with negative coupling with gravity (EQp and EQn respectively). All the models are normalized in order to match CMB data from WMAP3. We fit both the dark matter only and the total m ass profile with a NFW profile, and recover the concentration of each halo from the fi t using different definition. We consider both the complete catalog of clusters and groups and subsamples of objects at different level of relaxation. We find that the definition itself of the concentration can lead to differences up to 20% in its value and that these differences are smaller when more relaxed ˜

Strong magnetic fields and large rotation measures in protogalaxies from supernova seeding

Abstract: We present a model for the seeding and evolution of magnetic fi elds in protogalaxies. Super- nova (SN) explosions during the assembly of a protogalaxy self-consistently provide magnetic seed fields, which are subsequently amplified by compression , shear flows and random mo- tions. Our model explains the origin of strong magnetic field s ofG amplitude within the first starforming protogalactic structures shortly after the fir st stars have formed. We implement the model into the MHD version of the cosmological N-body / SPH simulation code GADGET and we couple the magnetic seeding directly to the underlying multi-phase description of star formation. We perform simulations of Milky Way-like galactic halo for- mation using a standardCDM cosmology and analyse the strength and distribution of the subsequent evolving magnetic field. Within starforming regions and given typical dimensions and magnetic field strengths in canonical SN remnants, we inject a dipole-shape magnetic fie ld at a rate of �10 9 G Gyr 1 . Subsequently, the magnetic field strength increases expone ntially on timescales of a few ten million years within the innermost regions of the halo. Furt hermore, turbulent diffusion, shocks and gas motions transport the magnetic field towards t he halo outskirts. At redshift z�0, the entire galactic halo is magnetized and the field amplit ude is of the order of a fewG in the center of the halo and �10 9 G at the virial radius. Additionally, we analyse the intrinsic rotation measure (R M) of the forming galactic halo over redshift. The mean halo intrinsic RM peaks between redshifts z�4 and z�2 and reaches abso- lute values around 1000 rad m 2 . While the halo virializes towards redshift z�0, the intrinsic RM values decline to a mean value below 10 rad m 2 . At high redshifts, the distribution of individual starforming and thus magnetized regions is widespread. This leads to a widespread distribution of large intrinsic RM values. In our model for the evolution of galactic magnetic fields, th e seed magnetic field amplitude and distribution is no longer a free parameter, but determin ed self-consistently by the star formation process occuring during the formation of cosmic structures. Thus, it provides a solution to the seed field problem.

Planck intermediate results - IV. The XMM-Newton validation programme for new Planck galaxy clusters

Abstract: We present the final results from the XMM-Newton validation follow-up of new Planck galaxy cluster candidates. We observed 15 new candidates, detected with signal-to-noise ratios between 4.0 and 6.1 in the 15.5-month nominal Planck survey. The candidates were selected using ancillary data flags derived from the ROSAT All Sky Survey (RASS) and Digitized Sky Survey all-sky maps, with the aim of pushing into the low SZ flux, high-z regime and testing RASS flags as indicators of candidate reliability. Fourteen new clusters were detected by XMM-Newton, ten single clusters and two double systems. Redshifts from X-ray spectroscopy lie in the range 0.2 to 0.9, with six clusters at z > 0.5. Estimated masses (M500) range from 2.5 × 10 14 to 8 × 10 14 M� . We discuss our results in the context of the full XMM-Newton validation programme, in which 51 new clusters have been detected. This includes four double and two triple systems, some of which are chance projections on the sky of clusters at different redshifts. We find that association with a source from the RASS-Bright Source Catalogue is a robust indicator of the reliability of a candidate, whereas association with a source from the RASS-Faint Source Catalogue does not guarantee that the SZ candidate is a bona fide cluster. Nevertheless, most Planck clusters appear in RASS maps, with a significance greater than 2σ being a good indication that the candidate is a real cluster. Candidate validation from association with SDSS galaxy overdensity at z > 0.5 is also discussed. The full sample gives a Planck sensitivity threshold of Y500 ∼ 4 × 10 −4 arcmin 2 , with indication for Malmquist bias in the YX–Y500 relation below this threshold. The corresponding mass ][ �

Numerical simulations of the possible origin of the two sub-parsec scale and counter-rotating stellar disks around SgrA*

Abstract: We present a high resolution simulation of an idealized model to explain the origin of the two young, counter-rotating, sub-parsec scale stellar disks around the supermassive black hole SgrA* at the Center of the Milky Way. In our model, the collision of a single molecular cloud with a circum-nuclear gas disk (similar to the one observed presently) leads to multiple streams of gas flowing towards the black hole and creating accretion disks with angular momentum depending on the ratio of cloud and circum-nuclear disk material. The infalling gas creates two inclined, counter-rotating sub-parsec scale accretion disks around the supermassive black hole with the first disk forming roughly 1 Myr earlier, allowing it to fragment into stars and get dispersed before the second, counter-rotating disk forms. Fragmentation of the second disk would lead to the two inclined, counter-rotating stellar disks which are observed at the Galactic Center. A similar event might be happening again right now at the Milky Way Galactic Center. Our model predicts that the collision event generates spiral-like filaments of gas, feeding the Galactic Center prior to disk formation with a geometry and inflow pattern that is in agreement with the structure of the so called mini-spiral that has been detected in the Galactic Center.

Galaxies going MAD: the galaxy-finder comparison project

Abstract: With the ever-increasing size and complexity of fully self-consistent simulations of galaxy formation within the framework of the cosmic web, the demands upon object finders for these simulations have simultaneously grown. To this extent we initiated the Halo-Finder Comparison Project that gathered together all the experts in the field and has so far led to two comparison papers, one for dark matter field haloes, and one for dark matter subhaloes. However, as state-of-the-art simulation codes are perfectly capable of not only following the formation and evolution of dark matter but also accounting for baryonic physics, i.e. gas hydrodynamics, star formation, stellar feedback, etc., object finders should also be capable of taking these additional physical processes into consideration. Here we report - for the first time - on a comparison of codes as applied to the Constrained Local UniversE Simulation (CLUES) of the formation of the Local Group which incorporates much of the physics relevant for galaxy formation. We compare both the properties of the three main galaxies in the simulation (representing the Milky Way, Andromeda and M33) and their satellite populations for a variety of halo finders ranging from phase space to velocity space to spherical overdensity based codes, including also a mere baryonic object finder. We obtain agreement amongst codes comparable to (if not better than) our previous comparisons - at least for the total, dark and stellar components of the objects. However, the diffuse gas content of the haloes shows great disparity, especially for low-mass satellite galaxies. This is primarily due to differences in the treatment of the thermal energy during the unbinding procedure. We acknowledge that the handling of gas in halo finders is something that needs to be dealt with carefully, and the precise treatment may depend sensitively upon the scientific problem being studied.