Source author record

K. Dolag

K. Dolag appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

104works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

104 published item(s)

preprint2025arXiv

Magnetic Fields in the Shapley Supercluster Core with POSSUM: Challenging Model Predictions

Faraday Rotation Measure (RM) Grids provide a sensitive means to trace magnetized plasma across a wide range of cosmic environments. We study the RM signal from the Shapley Supercluster Core (SSC), in order to constrain the magnetic field properties of the gas. The SSC region consists of two galaxy clusters A3558 and A3562, and two galaxy groups between them, at $z\simeq 0.048$. We combine RM Grid data with thermal Sunyaev-Zeldovich effect data, obtained from the POSSUM pilot survey, and Planck, respectively. To robustly determine the gas density, its magnetic field properties, and their correlation, we study the RM scatter in the SSC region and its behavior as a function of distance to the nearest cluster/group. We compare observational results with semi-analytic Gaussian random field models and more realistic cosmological MHD simulations. With a sky-density of 36 RMs/deg$^{2}$, we detect an excess RM scatter of $30.5\pm 4.6 \, \mathrm{rad/m^2}$ in the SSC region. Comparing with models, we find an average magnetic field strength of 1-3 $μ$G (in the groups and clusters). The RM scatter profile, derived from data ranging from 0.3-1.8 $r_{500}$ for all objects, is systematically flatter than expected compared to models, with $η<0.5$ being favored. Despite this discrepancy, we find that cosmological MHD simulations matched to the SSC structure most closely align with scenarios where the magnetic field is amplified by the turbulent velocity in the intercluster regions on scales $\lesssim 0.8\,r_{500}$. The dense RM grid and precision provided by POSSUM allows us to probe magnetized gas in the SSC clusters and groups on scales within and beyond their $r_{500}$. Flatter-than-expected RM scatter profiles reveal a significant challenge in reconciling observations with even the most realistic predictions from cosmological MHD simulations in the outskirts of interacting clusters.

preprint2024arXiv

Average X-ray properties of galaxy groups. From Milky Way-like halos to massive clusters

This study examines the average X-ray properties of massive halos at z< 0.2, covering the largest halo mass range to date, from Milky Way-like halos to massive clusters. The analysis is based on stacking in the eFEDS area of the GAMA galaxy group sample, validated with synthetic data that mimic observed eROSITA X-ray and GAMA optical data using Magneticum lightcones. Stacking was conducted in halo mass bins and tested for AGN and X-ray binary contamination, systematics in the halo mass proxy, and uncertainties in optical group centers. The study provides average X-ray surface brightness profiles in six mass bins, spanning Milky Way-like systems to poor clusters. The scatter in the X-ray luminosity-mass (LX-M) relation is attributed to gas concentration: low X-ray luminosity systems at fixed halo mass exhibit lower central gas concentrations than high-luminosity systems, consistent with Magneticum predictions. However, discrepancies in dark matter concentration arise, with Magneticum predicting undetected groups as older and more relaxed, while observations suggest the opposite. New LX-M relations are presented covering three decades of halo mass. These relations fit a single power law, aligning with previous studies. Magneticum matches observed gas distributions across all masses, whereas IllustrisTNG, EAGLE, Simba, and FLAMINGO exhibit significant discrepancies at various mass scales. Simulations calibrated on local galaxy properties accurately reproduce central galaxies but fail to capture gas properties. Conversely, simulations like Magneticum excel in gas predictions but produce overly massive central galaxies. Further exploration of gas and dark matter distributions and their effects on galaxy properties is critical to comprehending the role of gravitational forces and feedback in shaping large-scale structure and galaxy evolution.

preprint2024arXiv

The perils of stacking optically selected groups in eROSITA data. The Magneticum perspective

Hydrodynamical simulation predictions are often compared with observational data without fully accounting for systematics and biases specific to observational techniques. Using the magnetohydrodynamical simulation Magneticum, we generate mock eROSITA eRASS:4 data, combined with GAMA-like spectroscopic surveys and optically selected galaxy catalogs from the same light-cone, to analyze hot gas properties in galaxy groups via a stacking technique. This study aims to (i) incorporate observational systematics into predictions and (ii) evaluate the reliability of stacking techniques for determining average X-ray properties of galaxy groups. Our analysis provides X-ray emission predictions from Magneticum, including contributions from AGN, X-ray binaries (XRBs), and the Intra-Group Medium (IGM) as a function of halo mass, covering Milky Way (MW)-like groups to poor clusters. We find that AGN and XRBs dominate the X-ray surface brightness profiles of low-mass halos. The reliability of stacking techniques is tested by reproducing input X-ray surface brightness and electron density profiles, accounting for completeness and contamination of prior samples, miscentering of optical group centers, uncertainties in X-ray emissivity due to gas temperature and metallicity assumption, and systematics in halo mass proxies. The halo mass proxy emerges as the primary source of systematics, affecting X-ray surface brightness and scaling relations. We show that stacked X-ray luminosity-mass relations are flatter than input relations but consistent with observations. Additionally, the retrieved hot gas fraction-mass relation aligns well with observational data. These results highlight the need to account for systematic errors when comparing stacking techniques to other methods using different prior catalogs or predictions.

preprint2022arXiv

Mapping "out-of-the-box" the properties of the baryons in massive halos

We study the distributions of the baryons in massive halos ($M_{vir} > 10^{13} \ h^{-1}M_{\odot}$) in the $Magneticum$ suite of Smoothed Particle Hydrodynamical cosmological simulations, out to the unprecedented radial extent of $10 R_{500,\mathrm c}$. We confirm that, under the action of non-gravitational physical phenomena, the baryon mass fraction is lower in the inner regions ($<R_{500,\mathrm c}$) of increasingly less massive halos, and rises moving outwards, with values that spans from 51% (87%) in the regions around $R_{500,\mathrm c}$ to 95% (100%) at $10R_{500,\mathrm c}$ of the cosmological value in the systems with the lowest (highest; $M_{vir} \sim 5 \times 10^{14} \ h^{-1}M_{\odot}$) masses. The galaxy groups almost match the gas (and baryon) fraction measured in the most massive halos only at very large radii ($r>6 R_{500,\mathrm c}$), where the baryon depletion factor $Y_{\rm bar} = f_{\rm bar} / (Ω_{\rm b}/Ω_{\rm m})$ approaches the value of unity, expected for "closed-box" systems. We find that both the radial and mass dependency of the baryon, gas, and hot depletion factors are predictable and follow a simple functional form. The star mass fraction is higher in less massive systems, decreases systematically with increasing radii, and reaches a constant value of $Y_{\rm star} \approx 0.09$, where also the gas metallicity is constant, regardless of the host halo mass, as a result of the early ($z>2$) enrichment process.

preprint2021arXiv

Studying the merging cluster Abell 3266 with eROSITA

Abell 3266 is one of the X-ray brightest galaxy clusters in the sky and is a well-known merging system. Using the ability of the eROSITA telescope onboard SRG (Spectrum Röntgen Gamma) to observe a wide field with a single pointing, we analysed a new observation of the cluster out to a radius of R_200. The X-ray images highlight sub-structures present in the cluster, including the north-east--south-west merger seen in previous ASCA, Chandra, and XMM-Newton data, a merging group towards the north-west, and filamentary structures between the core and one or more groups towards the west. We compute spatially resolved spectroscopic maps of the thermodynamic properties of the cluster, including the metallicity. The merging subclusters are seen as low entropy material within the cluster. The filamentary structures could be the rims of a powerful outburst of an active galactic nucleus, or most likely material stripped from the western group(s) as they passed through the cluster core. Seen in two directions is a pressure jump at a radius of 1.1 Mpc, which is consistent with a shock with a Mach number of ~1.5-1.7. The eROSITA data confirm that the cluster is not a simple merging system, but it is made up of several subclusters which are merging or will shortly merge. We computed a hydrostatic mass from the eROSITA data, finding good agreement with a previous XMM-Newton result. With this pointing we detect several extended sources, where we find secure associations between z=0.36-1.0 for seven of them, that is background galaxy groups and clusters, highlighting the power of eROSITA to find such systems.

preprint2021arXiv

The eROSITA Final Equatorial-Depth Survey (eFEDS): Catalog of galaxy clusters and groups

The eROSITA Final Equatorial-Depth Survey has been carried out during the PV phase of the SRG/eROSITA telescope and completed in November 2019. This survey is designed to provide the first eROSITA-selected sample of galaxy clusters and groups and to test the predictions for the all-sky survey in the context of cosmological studies with clusters. In the 140 deg$^2$ area covered by eFEDS, 542 candidate clusters and groups are detected as extended X-ray sources, down to a flux of $\sim10^{-14} $erg/s/cm$^2$ in the soft band (0.5-2 keV) within 1'. In this work, we provide the catalog of candidate galaxy clusters and groups in eFEDS. We perform imaging and spectral analysis on the eFEDS clusters with eROSITA X-ray data, and study the properties of the sample. The clusters are distributed in the redshift range [0.01, 1.3], with the median redshift at 0.35. We obtain the ICM temperature measurement with $>2σ$ c.l. for $\sim$1/5 (102/542) of the sample. The average temperature of these clusters is $\sim$2 keV. Radial profiles of flux, luminosity, electron density, and gas mass are measured from the precise modeling of the imaging data. The selection function, the purity and completeness of the catalog are examined and discussed in detail. The contamination fraction is $\sim1/5$ in this sample, dominated by misidentified point sources. The X-ray Luminosity Function of the clusters agrees well with the results obtained from other recent X-ray surveys. We also find 19 supercluster candidates in eFEDS, most of which are located at redshifts between 0.1 and 0.5. The eFEDS cluster and group catalog provides a benchmark proof-of-concept for the eROSITA All-Sky Survey extended source detection and characterization. We confirm the excellent performance of eROSITA for cluster science and expect no significant deviations from our pre-launch expectations for the final All-Sky Survey.

preprint2020arXiv

The Abell 3391/95 galaxy cluster system: A 15 Mpc intergalactic medium emission filament, a warm gas bridge, infalling matter clumps, and (re-) accelerated plasma discovered by combining SRG/eROSITA data with ASKAP/EMU and DECam data

We used dedicated SRG/eROSITA X-ray, ASKAP/EMU radio, and DECam optical observations of a 15 sq.deg region around the interacting galaxy cluster system A3391/95 to study the warm-hot gas in cluster outskirts and filaments, the surrounding large-scale structure and its formation process. We relate the observations to expectations from cosmological hydrodynamic simulations from the Magneticum suite. We trace the irregular morphology of warm-hot gas of the main clusters from their centers out to well beyond their characteristic radii, $r_{200}$. Between the two main cluster systems, we observe an emission bridge; thanks to eROSITA's unique soft response and large field of view, we discover tantalizing hints for warm gas. Several matter clumps physically surrounding the system are detected. For the "Northern Clump," we provide evidence that it is falling towards A3391 from the hot gas morphology and radio lobe structure of its central AGN. Many of the extended sources in the field detected by eROSITA are known clusters or new clusters in the background, including a known SZ cluster at redshift z=1. We discover an emission filament north of the virial radius, $r_{100}$, of A3391 connecting to the Northern Clump and extending south of A3395 towards another galaxy cluster. The total projected length of this continuous warm-hot emission filament is 15 Mpc, running almost 4 degrees across the entire eROSITA observation. The DECam galaxy density map shows galaxy overdensities in the same regions. The new datasets provide impressive confirmation of the theoretically expected structure formation processes on the individual system level, including the surrounding warm-hot intergalactic medium distribution compared to the Magneticum simulation. Our spatially resolved findings show that baryons indeed reside in large-scale warm-hot gas filaments with a clumpy structure.

preprint2019arXiv

The Three Hundred Project: correcting for the hydrostatic-equilibrium mass bias in X-ray and SZ surveys

Accurate and precise measurements of masses of galaxy clusters are key to derive robust constraints on cosmological parameters. Rising evidence from observations, however, confirms that X-ray masses, obtained under the assumption of hydrostatic equilibrium, might be underestimated, as previously predicted by cosmological simulations. We analyse more than 300 simulated massive clusters, from `The Three Hundred Project', and investigate the connection between mass bias and several diagnostics extracted from synthetic X-ray images of these simulated clusters. We find that the azimuthal scatter measured in 12 sectors of the X-ray flux maps is a statistically significant indication of the presence of an intrinsic (i.e. 3D) clumpy gas distribution. We verify that a robust correction to the hydrostatic mass bias can be inferred when estimates of the gas inhomogeneity from X-ray maps (such as the azimuthal scatter or the gas ellipticity) are combined with the asymptotic external slope of the gas density or pressure profiles, which can be respectively derived from X-ray and millimetric (Sunyaev-Zeldovich effect) observations. We also obtain that mass measurements based on either gas density and temperature or gas density and pressure result in similar distributions of the mass bias. In both cases, we provide corrections that help reduce both the dispersion and skewness of the mass bias distribution. These are effective even when irregular clusters are included leading to interesting implications for the modelling and correction of hydrostatic mass bias in cosmological analyses of current and future X-ray and SZ cluster surveys.

preprint2016arXiv

Geometrical on-the-fly shock detection in SPH

We present an on-the-fly geometrical approach for shock detection and Mach number calculation in simulations employing smoothed particle hydrodynamics (SPH). We utilize pressure gradients to select shock candidates and define up- and downstream positions. We obtain hydrodynamical states in the up- and downstream regimes with a series of normal and inverted kernel weightings parallel and perpendicular to the shock normals. Our on-the-fly geometrical Mach detector incorporates well within the SPH formalism and has low computational cost. We implement our Mach detector into the simulation code GADGET and alongside many SPH improvements. We test our shock finder in a sequence of shock-tube tests with successively increasing Mach numbers exceeding by far the typical values inside galaxy clusters. For all shocks, we resolve the shocks well and the correct Mach numbers are assigned. An application to a strong magnetized shock-tube gives stable results in full magnetohydrodynamic set-ups. We simulate a merger of two idealized galaxy clusters and study the shock front. Shock structures within the merging clusters as well as the cluster shock are well-captured by our algorithm and assigned correct Mach numbers.

preprint2016arXiv

On the nature of hydrostatic equilibrium in galaxy clusters

In this paper we investigate the level of hydrostatic equilibrium (HE) in the intra-cluster medium of simulated galaxy clusters, extracted from state-of-the-art cosmological hydrodynamical simulations performed with the Smoothed-Particle-Hydrodynamic code GADGET-3. These simulations include several physical processes, among which stellar and AGN feedback, and have been performed with an improved version of the code that allows for a better description of hydrodynamical instabilities and gas mixing processes. Evaluating the radial balance between the gravitational and hydrodynamical forces, via the gas accelerations generated, we effectively examine the level of HE in every object of the sample, its dependence on the radial distance from the center and on the classification of the cluster in terms of either cool-coreness or dynamical state. We find an average deviation of 10-20% out to the virial radius, with no evident distinction between cool-core and non-cool-core clusters. Instead, we observe a clear separation between regular and disturbed systems, with a more significant deviation from HE for the disturbed objects. The investigation of the bias between the hydrostatic estimate and the total gravitating mass indicates that, on average, this traces very well the deviation from HE, even though individual cases show a more complex picture. Typically, in the radial ranges where mass bias and deviation from HE are substantially different, the gas is characterized by a significant amount of random motions (>~30 per cent), relative to thermal ones. As a general result, the HE-deviation and mass bias, at given interesting distance from the cluster center, are not very sensitive to the temperature inhomogeneities in the gas.

preprint2016arXiv

Planck 2015 results. XI. CMB power spectra, likelihoods, and robustness of parameters

This paper presents the Planck 2015 likelihoods, statistical descriptions of the 2-point correlations of CMB data, using the hybrid approach employed previously: pixel-based at $\ell<30$ and a Gaussian approximation to the distribution of spectra at higher $\ell$. The main improvements are the use of more and better processed data and of Planck polarization data, and more detailed foreground and instrumental models, allowing further checks and enhanced immunity to systematics. Progress in foreground modelling enables a larger sky fraction. Improvements in processing and instrumental models further reduce uncertainties. For temperature, we perform an analysis of end-to-end instrumental simulations fed into the data processing pipeline; this does not reveal biases from residual instrumental systematics. The $Λ$CDM cosmological model continues to offer a very good fit to Planck data. The slope of primordial scalar fluctuations, $n_s$, is confirmed smaller than unity at more than 5σ from Planck alone. We further validate robustness against specific extensions to the baseline cosmology. E.g., the effective number of neutrino species remains compatible with the canonical value of 3.046. This first detailed analysis of Planck polarization concentrates on E modes. At low $\ell$ we use temperature at all frequencies and a subset of polarization. The frequency range improves CMB-foreground separation. Within the baseline model this requires a reionization optical depth $τ=0.078\pm0.019$, significantly lower than without high-frequency data for explicit dust monitoring. At high $\ell$ we detect residual errors in E, typically O($μ$K$^2$); we recommend temperature alone as the high-$\ell$ baseline. Nevertheless, Planck high-$\ell$ polarization allows a separate determination of $Λ$CDM parameters consistent with those from temperature alone.

preprint2016arXiv

Planck 2015 results. XIII. Cosmological parameters

We present results based on full-mission Planck observations of temperature and polarization anisotropies of the CMB. These data are consistent with the six-parameter inflationary LCDM cosmology. From the Planck temperature and lensing data, for this cosmology we find a Hubble constant, H0= (67.8 +/- 0.9) km/s/Mpc, a matter density parameter Omega_m = 0.308 +/- 0.012 and a scalar spectral index with n_s = 0.968 +/- 0.006. (We quote 68% errors on measured parameters and 95% limits on other parameters.) Combined with Planck temperature and lensing data, Planck LFI polarization measurements lead to a reionization optical depth of tau = 0.066 +/- 0.016. Combining Planck with other astrophysical data we find N_ eff = 3.15 +/- 0.23 for the effective number of relativistic degrees of freedom and the sum of neutrino masses is constrained to < 0.23 eV. Spatial curvature is found to be |Omega_K| < 0.005. For LCDM we find a limit on the tensor-to-scalar ratio of r <0.11 consistent with the B-mode constraints from an analysis of BICEP2, Keck Array, and Planck (BKP) data. Adding the BKP data leads to a tighter constraint of r < 0.09. We find no evidence for isocurvature perturbations or cosmic defects. The equation of state of dark energy is constrained to w = -1.006 +/- 0.045. Standard big bang nucleosynthesis predictions for the Planck LCDM cosmology are in excellent agreement with observations. We investigate annihilating dark matter and deviations from standard recombination, finding no evidence for new physics. The Planck results for base LCDM are in agreement with BAO data and with the JLA SNe sample. However the amplitude of the fluctuations is found to be higher than inferred from rich cluster counts and weak gravitational lensing. Apart from these tensions, the base LCDM cosmology provides an excellent description of the Planck CMB observations and many other astrophysical data sets.

preprint2016arXiv

Planck 2015 results. XIX. Constraints on primordial magnetic fields

We compute and investigate four types of imprint of a stochastic background of primordial magnetic fields (PMFs) on the cosmic microwave background (CMB) anisotropies: the impact of PMFs on the CMB spectra; the effect on CMB polarization induced by Faraday rotation; the impact of PMFs on the ionization history; magnetically-induced non-Gaussianities; and the magnetically-induced breaking of statistical isotropy. Overall, Planck data constrain the amplitude of PMFs to less than a few nanogauss. In particular, individual limits coming from the analysis of the CMB angular power spectra, using the Planck likelihood, are $B_{1\,\mathrm{Mpc}}< 4.4$ nG (where $B_{1\,\mathrm{Mpc}}$ is the comoving field amplitude at a scale of 1 Mpc) at 95% confidence level, assuming zero helicity, and $B_{1\,\mathrm{Mpc}}< 5.6$ nG for a maximally helical field.For nearly scale-invariant PMFs we obtain $B_{1\,\mathrm{Mpc}}<2.0$ nG and $B_{1\,\mathrm{Mpc}}<0.9$ nG if the impact of PMFs on the ionization history of the Universe is included. From the analysis of magnetically-induced non-Gaussianity we obtain three different values, corresponding to three applied methods, all below 5 nG. The constraint from the magnetically-induced passive-tensor bispectrum is $B_{1\,\mathrm{Mpc}}< 2.8$ nG. A search for preferred directions in the magnetically-induced passive bispectrum yields $B_{1\,\mathrm{Mpc}}< 4.5$ nG, whereas the the compensated-scalar bispectrum gives $B_{1\,\mathrm{Mpc}}< 3$ nG. The analysis of the Faraday rotation of CMB polarization by PMFs uses the Planck power spectra in $EE$ and $BB$ at 70 GHz and gives $B_{1\,\mathrm{Mpc}}< 1380$ nG. In our final analysis, we consider the harmonic-space correlations produced by Alfvén waves, finding no significant evidence for the presence of these waves. Together, these results comprise a comprehensive set of constraints on possible PMFs with Planck data.

preprint2016arXiv

Planck intermediate results. XLII. Large-scale Galactic magnetic fields

Recent models for the large-scale Galactic magnetic fields in the literature have been largely constrained by synchrotron emission and Faraday rotation measures. We use three different but representative models to compare their predicted polarized synchrotron and dust emission with that measured by the Planck satellite. We first update these models to match the Planck synchrotron products using a common model for the cosmic-ray leptons. We discuss the impact on this analysis of the ongoing problems of component separation in the Planck microwave bands and of the uncertain cosmic-ray spectrum. In particular, the inferred degree of ordering in the magnetic fields is sensitive to these systematic uncertainties, and we further show the importance of considering the expected variations in the observables in addition to their mean morphology. We then compare the resulting simulated emission to the observed dust polarization and find that the dust predictions do not match the morphology in the Planck data but underpredict the dust polarization away from the plane. We modify one of the models to roughly match both observables at high latitudes by increasing the field ordering in the thin disc near the observer. Though this specific analysis is dependent on the component separation issues, we present the improved model as a proof of concept for how these studies can be advanced in future using complementary information from ongoing and planned observational projects.

preprint2015arXiv

An improved SPH scheme for cosmological simulations

We present an implementation of smoothed particle hydrodynamics (SPH) with improved accuracy for simulations of galaxies and the large-scale structure. In particular, we combine, implement, modify and test a vast majority of SPH improvement techniques in the latest instalment of the GADGET code. We use the Wendland kernel functions, a particle wake-up time-step limiting mechanism and a time-dependent scheme for artificial viscosity, which includes a high-order gradient computation and shear flow limiter. Additionally, we include a novel prescription for time-dependent artificial conduction, which corrects for gravitationally induced pressure gradients and largely improves the SPH performance in capturing the development of gas-dynamical instabilities. We extensively test our new implementation in a wide range of hydrodynamical standard tests including weak and strong shocks as well as shear flows, turbulent spectra, gas mixing, hydrostatic equilibria and self-gravitating gas clouds. We jointly employ all modifications; however, when necessary we study the performance of individual code modules. We approximate hydrodynamical states more accurately and with significantly less noise than standard SPH. Furthermore, the new implementation promotes the mixing of entropy between different fluid phases, also within cosmological simulations. Finally, we study the performance of the hydrodynamical solver in the context of radiative galaxy formation and non-radiative galaxy cluster formation. We find galactic disks to be colder, thinner and more extended and our results on galaxy clusters show entropy cores instead of steadily declining entropy profiles. In summary, we demonstrate that our improved SPH implementation overcomes most of the undesirable limitations of standard SPH, thus becoming the core of an efficient code for large cosmological simulations.

preprint2015arXiv

Connecting Angular Momentum and Galactic Dynamics: The complex Interplay between Spin, Mass, and Morphology

The evolution and distribution of the angular momentum of dark matter (DM) halos have been discussed in several studies over the past decades. In particular, the idea arose that angular momentum conservation should allow to infer the total angular momentum of the entire DM halo from measuring the angular momentum of the baryonic component, which is populating the center of the halo, especially for disk galaxies. To test this idea and to understand the connection between the angular momentum of the DM halo and its galaxy, we use the Magneticum simulations. We successfully produce populations of spheroidal and disk galaxies self-consistently. Thus, we are able to study the dependence of galactic properties on their morphology. We find that (1) the specific angular momentum of stars in disk and spheroidal galaxies as a function of their stellar mass compares well with observational results; (2) the specific angular momentum of the stars in disk galaxies is slightly smaller compared to the specific angular momentum of the cold gas, in good agreement with observations; (3) simulations including the baryonic component show a dichotomy in the specific stellar angular momentum distribution when splitting the galaxies according to their morphological type (this dichotomy can also be seen in the spin parameter, where disk galaxies populate halos with slightly larger spin compared to spheroidal galaxies); (4) disk galaxies preferentially populate halos in which the angular momentum vector of the DM component in the central part shows a better alignment to the angular momentum vector of the entire halo; and (5) the specific angular momentum of the cold gas in disk galaxies is approximately 40 percent smaller than the specific angular momentum of the total DM halo and shows a significant scatter.

preprint2015arXiv

Constraints on the distribution and energetics of fast radio bursts using cosmological hydrodynamic simulations

We present constraints on the origins of fast radio bursts (FRBs) using large cosmological simulations. We calculate contributions to FRB dispersion measures (DMs) from the Milky Way, from the local Universe, from cosmological large-scale structure, and from potential FRB host galaxies, and then compare these simulations to the DMs of observed FRBs. We find that the Milky Way contribution has previously been underestimated by a factor of ~2, and that the foreground-subtracted DMs are consistent with a cosmological origin, corresponding to a source population observable to a maximum redshift z~0.6-0.9. We consider models for the spatial distribution of FRBs in which they are randomly distributed in the Universe, track the star-formation rate of their host galaxies, track total stellar mass, or require a central supermassive black hole. Current data do not discriminate between these possibilities, but the predicted DM distributions for different models will differ considerably once we begin detecting FRBs at higher DMs and higher redshifts. We additionally consider the distribution of FRB fluences, and show that the observations are consistent with FRBs being standard candles, each burst producing the same radiated isotropic energy. The data imply a constant isotropic burst energy of ~7e40 erg if FRBs are embedded in host galaxies, or ~9e40 erg if FRBs are randomly distributed. These energies are 10-100 times larger than had previously been inferred. Within the constraints of the available small sample of data, our analysis favours FRB mechanisms for which the isotropic radiated energy has a narrow distribution in excess of 1e40 erg.

preprint2015arXiv

Cool Core Clusters from Cosmological Simulations

We present results obtained from a set of cosmological hydrodynamic simulations of galaxy clusters, aimed at comparing predictions with observational data on the diversity between cool-core (CC) and non-cool-core (NCC) clusters. Our simulations include the effects of stellar and AGN feedback and are based on an improved version of the smoothed particle hydrodynamics code GADGET-3, which ameliorates gas mixing and better captures gas-dynamical instabilities by including a suitable artificial thermal diffusion. In this Letter, we focus our analysis on the entropy profiles, the primary diagnostic we used to classify the degree of cool-coreness of clusters, and on the iron profiles. In keeping with observations, our simulated clusters display a variety of behaviors in entropy profiles: they range from steadily decreasing profiles at small radii, characteristic of cool-core systems, to nearly flat core isentropic profiles, characteristic of non-cool-core systems. Using observational criteria to distinguish between the two classes of objects, we find that they occur in similar proportions in both simulations and in observations. Furthermore, we also find that simulated cool-core clusters have profiles of iron abundance that are steeper than those of NCC clusters, which is also in agreement with observational results. We show that the capability of our simulations to generate a realistic cool-core structure in the cluster population is due to AGN feedback and artificial thermal diffusion: their combined action allows us to naturally distribute the energy extracted from super-massive black holes and to compensate for the radiative losses of low-entropy gas with short cooling time residing in the cluster core.

preprint2015arXiv

Planck 2013 results. XXIX. The Planck catalogue of Sunyaev-Zeldovich sources: Addendum

We update the all-sky Planck catalogue of 1227 clusters and cluster candidates (PSZ1) published in March 2013, derived from Sunyaev-Zeldovich (SZ) effect detections using the first 15.5 months of Planck satellite observations. Addendum. We deliver an updated version of the PSZ1 catalogue, reporting the further confirmation of 86 Planck-discovered clusters. In total, the PSZ1 now contains 947 confirmed clusters, of which 214 were confirmed as newly discovered clusters through follow-up observations undertaken by the Planck Collaboration. The updated PSZ1 contains redshifts for 913 systems, of which 736 (~80.6%) are spectroscopic, and associated mass estimates derived from the Y_z mass proxy. We also provide a new SZ quality flag, derived from a novel artificial neural network classification of the SZ signal, for the remaining 280 candidates. Based on this assessment, the purity of the updated PSZ1 catalogue is estimated to be 94%. In this release, we provide the full updated catalogue and an additional readme file with further information on the Planck SZ detections.

preprint2015arXiv

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

The European Space Agency's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14~May 2009 and scanned the microwave and submillimetre sky continuously between 12~August 2009 and 23~October 2013. In February~2015, ESA and the Planck Collaboration released the second set of cosmology products based on data from the entire Planck mission, including both temperature and polarization, along with a set of scientific and technical papers and a web-based explanatory supplement. This paper gives an overview of the main characteristics of the data and the data products in the release, as well as the associated cosmological and astrophysical science results and papers. The science products include maps of the cosmic microwave background (CMB), the thermal Sunyaev-Zeldovich effect, and diffuse foregrounds in temperature and polarization, catalogues of compact Galactic and extragalactic sources (including separate catalogues of Sunyaev-Zeldovich clusters and Galactic cold clumps), and extensive simulations of signals and noise used in assessing the performance of the analysis methods and assessment of uncertainties. The likelihood code used to assess cosmological models against the Planck data are described, as well as a CMB lensing likelihood. Scientific results include cosmological parameters deriving from CMB power spectra, gravitational lensing, and cluster counts, as well as constraints on inflation, non-Gaussianity, primordial magnetic fields, dark energy, and modified gravity.

preprint2015arXiv

Planck 2015 results. XII. Full Focal Plane simulations

We present the 8th Full Focal Plane simulation set (FFP8), deployed in support of the Planck 2015 results. FFP8 consists of 10 fiducial mission realizations reduced to 18144 maps, together with the most massive suite of Monte Carlo realizations of instrument noise and CMB ever generated, comprising $10^4$ mission realizations reduced to about $10^6$ maps. The resulting maps incorporate the dominant instrumental, scanning, and data analysis effects; remaining subdominant effects will be included in future updates. Generated at a cost of some 25 million CPU-hours spread across multiple high-performance-computing (HPC) platforms, FFP8 is used for the validation and verification of analysis algorithms, as well as their implementations, and for removing biases from and quantifying uncertainties in the results of analyses of the real data.

preprint2015arXiv

Planck 2015 results. XXII. A map of the thermal Sunyaev-Zeldovich effect

We have constructed all-sky y-maps of the thermal Sunyaev-Zeldovich (tSZ) effect by applying specifically tailored component separation algorithms to the 30 to 857 GHz frequency channel maps from the Planck satellite survey. These reconstructed y-maps are delivered as part of the Planck 2015 release. The y-maps are characterised in terms of noise properties and residual foreground contamination, mainly thermal dust emission at large angular scales and CIB and extragalactic point sources at small angular scales. Specific masks are defined to minimize foreground residuals and systematics. Using these masks we compute the y-map angular power spectrum and higher order statistics. From these we conclude that the y-map is dominated by tSZ signal in the multipole range, 20-600. We compare the measured tSZ power spectrum and higher order statistics to various physically motivated models and discuss the implications of our results in terms of cluster physics and cosmology.

preprint2015arXiv

Planck 2015 results. XXVII. The Second Planck Catalogue of Sunyaev-Zeldovich Sources

We present the all-sky Planck catalogue of Sunyaev-Zeldovich (SZ) sources detected from the 29 month full-mission data. The catalogue (PSZ2) is the largest SZ-selected sample of galaxy clusters yet produced and the deepest all-sky catalogue of galaxy clusters. It contains 1653 detections, of which 1203 are confirmed clusters with identified counterparts in external data-sets, and is the first SZ-selected cluster survey containing > $10^3$ confirmed clusters. We present a detailed analysis of the survey selection function in terms of its completeness and statistical reliability, placing a lower limit of 83% on the purity. Using simulations, we find that the Y5R500 estimates are robust to pressure-profile variation and beam systematics, but accurate conversion to Y500 requires. the use of prior information on the cluster extent. We describe the multi-wavelength search for counterparts in ancillary data, which makes use of radio, microwave, infra-red, optical and X-ray data-sets, and which places emphasis on the robustness of the counterpart match. We discuss the physical properties of the new sample and identify a population of low-redshift X-ray under- luminous clusters revealed by SZ selection. These objects appear in optical and SZ surveys with consistent properties for their mass, but are almost absent from ROSAT X-ray selected samples.

preprint2015arXiv

Planck intermediate results. XL. The Sunyaev-Zeldovich signal from the Virgo cluster

The Virgo cluster is the largest Sunyaev-Zeldovich (SZ) source in the sky, both in terms of angular size and total integrated flux. Planck's wide angular scale and frequency coverage, together with its high sensitivity, allow a detailed study of this large object through the SZ effect. Virgo is well resolved by Planck, showing an elongated structure, which correlates well with the morphology observed from X-rays, but extends beyond the observed X-ray signal. We find a good agreement between the SZ signal (or Compton paranmeter, y_c) observed by Planck and the expected signal inferred from X-ray observations and simple analytical models. Due to its proximity to us, the gas beyond the virial radius can be studied with unprecedented sensitivity by integrating the SZ signal over tens of square degrees. We study the signal in the outskirts of Virgo and compare it with analytical models and a constrained simulation of the environment of Virgo. Planck data suggest that significant amounts of low-density plasma surround Virgo out to twice the virial radius. We find the SZ signal in the outskirts of Virgo to be consistent with a simple model that extrapolates the inferred pressure at lower radii while assuming that the temperature stays in the keV range beyond the virial radius. The observed signal is also consistent with simulations and points to a shallow pressure profile in the outskirts of the cluster. This reservoir of gas at large radii can be linked with the hottest phase of the elusive warm/hot intergalactic medium. Taking the lack of symmetry of Virgo into account, we find that a prolate model is favoured by the combination of SZ and X-ray data, in agreement with predictions.

preprint2015arXiv

Planck intermediate results. XXXVII. Evidence of unbound gas from the kinetic Sunyaev-Zeldovich effect

By looking at the kinetic Sunyaev-Zeldovich effect (kSZ) in Planck nominal mission data, we present a significant detection of baryons participating in large-scale bulk flows around central galaxies (CGs) at redshift $z\approx 0.1$. We estimate the pairwise momentum of the kSZ temperature fluctuations at the positions of the CGC (Central Galaxy Catalogue) samples extracted from Sloan Digital Sky Survey (DR7) data. For the foreground-cleaned maps, we find $1.8$-$2.5σ$ detections of the kSZ signal, which are consistent with the kSZ evidence found in individual Planck raw frequency maps, although lower than found in the WMAP-9yr W band ($3.3σ$). We further reconstruct the peculiar velocity field from the CG density field, and compute for the first time the cross-correlation function between kSZ temperature fluctuations and estimates of CG radial peculiar velocities. This correlation function yields a $3.0$-$3.7$$σ$ detection of the peculiar motion of extended gas on Mpc scales, in flows correlated up to distances of 80-100 $h^{-1}$ Mpc. Both the pairwise momentum estimates and kSZ temperature-velocity field correlation find evidence for kSZ signatures out to apertures of 8 arcmin and beyond, corresponding to a physical radius of $> 1$ Mpc, more than twice the mean virial radius of halos. This is consistent with the predictions from hydro simulations that most of the baryons are outside the virialized halos. We fit a simple model, in which the temperature-velocity cross-correlation is proportional to the signal seen in a semi-analytic model built upon N-body simulations, and interpret the proportionality constant as an "effective" optical depth to Thomson scattering. We find $τ_T=(1.4\pm0.5)\times 10^{-4}$; the simplest interpretation of this measurement is that much of the gas is in a diffuse phase, which contributes little signal to X-ray or thermal SZ observations.

preprint2014arXiv

Cluster Radio Halos at the crossroads between astrophysics and cosmology in the SKA era

Giant Radio Halos (RH) are diffuse, Mpc-sized, synchrotron radio sources observed in a fraction of merging galaxy clusters. The current scenario for the origin of RHs assumes that turbulence generated during cluster mergers re-accelerates pre-existing fossil and/or secondary electrons in the intra-cluster-medium (ICM) to the energies necessary to produce the observed radio emission. Moreover, more relaxed clusters could host diffuse "off state" halos produced by secondary electrons. In this Chapter we use Monte Carlo simulations, that combine turbulent-acceleration physics and the generation of secondaries in the ICM, to calculate the occurrence of RHs in the Universe, their spectral properties and connection with properties of the hosting clusters. Predictions for SKA1 surveys are presented at low (100-300 MHz) and mid (1-2 GHz) frequencies assuming the expected sensitivities and spatial resolutions of SKA1. SKA1 will step into an unexplored territory allowing us to study the formation and evolution of RHs in a totally new range of cluster masses and redshift, allowing firm tests of the current theoretical hypothesis. In particular, the combination of SKA1-LOW and SUR will allow the discovery of ~1000 ultrasteep- spectrum halos and to detect for the very first time "off state" RHs. We expect that at least ~2500 giant RHs will be discovered by SKA1-LOW surveys up to z~0.6. Remarkably these surveys will be sensitive to RHs in a cluster mass range (down to ~10^14 solar masses) and redshifts (up to ~1) that are unexplored by current observations. SKA1 surveys will be highly competitive with present and future SZ-surveys in the detection of high-redshift massive objects.

preprint2014arXiv

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

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.

preprint2014arXiv

Planck 2013 results. XVI. Cosmological parameters

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.

preprint2014arXiv

Planck 2013 results. XX. Cosmology from Sunyaev-Zeldovich cluster counts

We present constraints on cosmological parameters using number counts as a function of redshift for a sub-sample of 189 galaxy clusters from the Planck SZ (PSZ) catalogue. The PSZ is selected through the signature of the Sunyaev--Zeldovich (SZ) effect, and the sub-sample used here has a signal-to-noise threshold of seven, with each object confirmed as a cluster and all but one with a redshift estimate. We discuss the completeness of the sample and our construction of a likelihood analysis. Using a relation between mass $M$ and SZ signal $Y$ calibrated to X-ray measurements, we derive constraints on the power spectrum amplitude $σ_8$ and matter density parameter $Ω_{\mathrm{m}}$ in a flat $Λ$CDM model. We test the robustness of our estimates and find that possible biases in the $Y$--$M$ relation and the halo mass function are larger than the statistical uncertainties from the cluster sample. Assuming the X-ray determined mass to be biased low relative to the true mass by between zero and 30%, motivated by comparison of the observed mass scaling relations to those from a set of numerical simulations, we find that $σ_8=0.75\pm 0.03$, $Ω_{\mathrm{m}}=0.29\pm 0.02$, and $σ_8(Ω_{\mathrm{m}}/0.27)^{0.3} = 0.764 \pm 0.025$. The value of $σ_8$ is degenerate with the mass bias; if the latter is fixed to a value of 20% we find $σ_8(Ω_{\mathrm{m}}/0.27)^{0.3}=0.78\pm 0.01$ and a tighter one-dimensional range $σ_8=0.77\pm 0.02$. We find that the larger values of $σ_8$ and $Ω_{\mathrm{m}}$ preferred by Planck's measurements of the primary CMB anisotropies can be accommodated by a mass bias of about 40%. Alternatively, consistency with the primary CMB constraints can be achieved by inclusion of processes that suppress power on small scales relative to the $Λ$CDM model, such as a component of massive neutrinos (abridged).

preprint2014arXiv

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

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 < 60$), the major foreground contaminant is the diffuse thermal dust emission. At small 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 θ\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: $σ_8(Ω_{\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 $σ_{8}$.

preprint2014arXiv

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

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.

preprint2014arXiv

Planck intermediate results. XIII. Constraints on peculiar velocities

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 $Λ$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 $Λ$CDM scenario.

preprint2014arXiv

Planck intermediate results. XXX. The angular power spectrum of polarized dust emission at intermediate and high Galactic latitudes

The polarized thermal emission from Galactic dust is the main foreground present in measurements of the polarization of the cosmic microwave background (CMB) at frequencies above 100GHz. We exploit the Planck HFI polarization data from 100 to 353GHz to measure the dust angular power spectra $C_\ell^{EE,BB}$ over the range $40<\ell<600$ well away from the Galactic plane. These will bring new insights into interstellar dust physics and a precise determination of the level of contamination for CMB polarization experiments. We show that statistical properties of the emission can be characterized over large fractions of the sky using $C_\ell$. For the dust, they are well described by power laws in $\ell$ with exponents $α^{EE,BB}=-2.42\pm0.02$. The amplitudes of the polarization $C_\ell$ vary with the average brightness in a way similar to the intensity ones. The dust polarization frequency dependence is consistent with modified blackbody emission with $β_d=1.59$ and $T_d=19.6$K. We find a systematic ratio between the amplitudes of the Galactic $B$- and $E$-modes of 0.5. We show that even in the faintest dust-emitting regions there are no "clean" windows where primordial CMB $B$-mode polarization could be measured without subtraction of dust emission. Finally, we investigate the level of dust polarization in the BICEP2 experiment field. Extrapolation of the Planck 353GHz data to 150GHz gives a dust power $\ell(\ell+1)C_\ell^{BB}/(2π)$ of $1.32\times10^{-2}μ$K$_{CMB}^2$ over the $40<\ell<120$ range; the statistical uncertainty is $\pm0.29$ and there is an additional uncertainty (+0.28,-0.24) from the extrapolation, both in the same units. This is the same magnitude as reported by BICEP2 over this $\ell$ range, which highlights the need for assessment of the polarized dust signal even in the cleanest windows of the sky.

preprint2013arXiv

A numerical study of the effects of primordial non-Gaussianities on weak lensing statistics

While usually cosmological initial conditions are assumed to be Gaussian, inflationary theories can predict a certain amount of primordial non-Gaussianity which can have an impact on the statistical properties of the lensing observables. In order to evaluate this effect, we build a large set of realistic maps of different lensing quantities starting from light-cones extracted from large dark-matter only N-body simulations with initial conditions corresponding to different levels of primordial local non-Gaussianity strength $f_{\rm NL}$. Considering various statistical quantities (PDF, power spectrum, shear in aperture, skewness and bispectrum) we find that the effect produced by the presence of primordial non-Gaussianity is relatively small, being of the order of few per cent for values of $|f_{\rm NL}|$ compatible with the present CMB constraints and reaching at most 10-15 per cent for the most extreme cases with $|f_{\rm NL}|=1000$. We also discuss the degeneracy of this effect with the uncertainties due to the power spectrum normalization $σ_8$ and matter density parameter $Ω_{\rm m}$, finding that an error in the determination of $σ_8$ ($Ω_{\rm m}$) of about 3 (10) per cent gives differences comparable with non-Gaussian models having $f_{\rm NL}=\pm 1000$. These results suggest that the possible presence of an amount of primordial non-Gaussianity corresponding to $|f_{\rm NL}|=100$ is not hampering a robust determination of the main cosmological parameters in present and future weak lensing surveys, while a positive detection of deviations from the Gaussian hypothesis is possible only breaking the degeneracy with other cosmological parameters and using data from deep surveys covering a large fraction of the sky.

preprint2013arXiv

Characterizing Diffused Stellar Light in simulated galaxy clusters

[Abridged] In this paper, we carry out a detailed analysis of the performance of two different methods to identify the diffuse stellar light in cosmological hydrodynamical simulations of galaxy clusters. One method is based on a dynamical analysis of the stellar component. The second method is closer to techniques commonly employed in observational studies. Both the dynamical method and the method based on the surface brightness limit criterion are applied to the same set of hydrodynamical simulations for a large sample about 80 galaxy clusters. We find significant differences between the ICL and DSC fractions computed with the two corresponding methods, which amounts to about a factor of two for the AGN simulations, and a factor of four for the CSF set. We also find that the inclusion of AGN feedback boosts the DSC and ICL fractions by a factor of 1.5-2, respectively, while leaving the BCG+ICL and BCG+DSC mass fraction almost unchanged. The sum of the BCG and DSC mass stellar mass fraction is found to decrease from ~80 per cent in galaxy groups to ~60 per cent in rich clusters, thus in excess of what found from observational analysis. We identify the average surface brightness limits that yields the ICL fraction from the SBL method close to the DSC fraction from the dynamical method. These surface brightness limits turn out to be brighter in the CSF than in the AGN simulations. This is consistent with the finding that AGN feedback makes BCGs to be less massive and with shallower density profiles than in the CSF simulations. The BCG stellar component, as identified by both methods, are slightly older and more metal-rich than the stars in the diffuse component.

preprint2013arXiv

Constraints on the CMB Temperature Evolution using Multi-Band Measurements of the Sunyaev Zel'dovich Effect with the South Pole Telescope

The adiabatic evolution of the temperature of the cosmic microwave background (CMB) is a key prediction of standard cosmology. We study deviations from the expected adiabatic evolution of the CMB temperature of the form $T(z) =T_0(1+z)^{1-α}$ using measurements of the spectrum of the Sunyaev Zel'dovich Effect with the South Pole Telescope (SPT). We present a method for using the ratio of the Sunyaev Zel'dovich signal measured at 95 and 150 GHz in the SPT data to constrain the temperature of the CMB. We demonstrate that this approach provides unbiased results using mock observations of clusters from a new set of hydrodynamical simulations. We apply this method to a sample of 158 SPT-selected clusters, spanning the redshift range $0.05 < z < 1.35$, and measure $α= 0.017^{+0.030}_{-0.028}$, consistent with the standard model prediction of $α=0$. In combination with other published results, we constrain $α= 0.011 \pm 0.016$, an improvement of $\sim 20\%$ over published constraints. This measurement also provides a strong constraint on the effective equation of state in models of decaying dark energy $w_\mathrm{eff} = -0.987^{+0.016}_{-0.017}$.

preprint2013arXiv

Is the Sunyaev-Zeldovich effect responsible for the observed steepening in the spectrum of the Coma radio halo ?

The spectrum of the radio halo in the Coma cluster is measured over almost two decades in frequency. The current radio data show a steepening of the spectrum at higher frequencies, which has implications for models of the radio halo origin. There is an on-going debate on the possibility that the observed steepening is not intrinsic to the emitted radiation, but is instead caused by the SZ effect. Recently, the Planck satellite measured the SZ signal and its spatial distribution in the Coma cluster allowing to test this hypothesis. Using the Planck results, we calculated the modification of the radio halo spectrum by the SZ effect in three different ways. With the first two methods we measured the SZ-decrement within the aperture radii used for flux measurements of the halo at the different frequencies. First we adopted the global compilation of data from Thierbach et al. and a reference aperture radius consistent with those used by the various authors. Second we used the available brightness profiles of the halo at different frequencies to derive the spectrum within two fixed apertures, and derived the SZ-decrement using these apertures. As a third method we used the quasi-linear correlation between the y and the radio-halo brightness at 330 MHz discovered by Planck to derive the modification of the radio spectrum by the SZ-decrement in a way that is almost independent of the adopted aperture radius. We found that the spectral modification induced by the SZ-decrement is 4-5 times smaller than that necessary to explain the observed steepening. Consequently a break or cut-off in the spectrum of the emitting electrons is necessary to explain current data. We also show that, if a steepening is absent from the emitted spectrum, future deep observations at 5 GHz with single dishes are expected to measure a halo flux in a 40 arcmin radius that would be 7-8 times higher than currently seen.

preprint2013arXiv

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

Using a set of hydrodynamical simulations of 62 galaxy clusters and groups we study the ICM 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, C_R, that quantifies the large-scale inhomogeneity of the ICM. After 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 from the X-ray emission, which translates into a systematic overestimate of 6 (12)% in the measurement of M_gas at R_200 for our relaxed (perturbed) cluster sample. At the same time, the increase of C_R with radius introduces also a ~2% systematic underestimate in the measurement of the hydrostatic-equilibrium mass (M_he), which adds to the previous one generating a systematic ~8.5% overestimate in f_gas 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 and in the y-parameter profiles. We find that their correlation is highly significant (r_S = 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 C_R. After providing a function that connects the two quantities, we obtain that correcting the observed gas density profiles using the azimuthal scatter eliminates the bias in the measurement of M_gas for relaxed objects, which becomes (0+/-2)% up to 2R_200, and reduces it by a factor of 3 for perturbed ones. This method allows also to eliminate the systematics on the measurements of M_he and f_gas, although a significant halo to halo scatter remains. (abridged)

preprint2013arXiv

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

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.

preprint2013arXiv

On the role of AGN feedback on the thermal and chemodynamical properties of the hot intra-cluster medium

We present an analysis of the properties of the ICM in an extended set of cosmological hydrodynamical simulations of galaxy clusters and groups performed with the TreePM+SPH GADGET-3 code. Besides a set of non-radiative simulations, we carried out two sets of simulations including radiative cooling, star formation, metal enrichment and feedback from supernovae, one of which also accounts for the effect of feedback from AGN resulting from gas accretion onto super-massive black holes. These simulations are analysed with the aim of studying the relative role played by SN and AGN feedback on the general properties of the diffuse hot baryons in galaxy clusters and groups: scaling relations, temperature, entropy and pressure radial profiles, and ICM chemical enrichment. We find that simulations including AGN feedback produce scaling relations that are in good agreement with X-ray observations at all mass scales. However, our simulations are not able to account for the observed diversity between CC and NCC clusters: unlike for observations, we find that temperature and entropy profiles of relaxed and unrelaxed clusters are quite similar and resemble more the observed behaviour of NCC clusters. As for the pattern of metal enrichment, we find that an enhanced level of iron abundance is produced by AGN feedback with respect to the case of purely SN feedback. As a result, while simulations including AGN produce values of iron abundance in groups in agreement with observations, they over-enrich the ICM in massive clusters. The efficiency of AGN feedback in displacing enriched gas from halos into the inter-galactic medium at high redshift also creates a widespread enrichment in the outskirts of clusters and produces profiles of iron abundance whose slope is in better agreement with observations.

preprint2013arXiv

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

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.

preprint2013arXiv

Revisiting scaling relations for giant radio halos in galaxy clusters

Many galaxy clusters host Megaparsec-scale radio halos, generated by ultrarelativistic electrons in the magnetized intracluster medium. Correlations between the 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 SZ catalog, to revisit the correlations between the power of halos and the thermal properties of galaxy clusters. We find that the radio power of halos at 1.4 GHz scales with the cluster X-ray (0.1--2.4 keV) luminosity computed within R_500 as P_1.4 L_500^2.0. Our bigger and more homogenous sample confirms that the X-ray luminous (L_500 > 5x10^44 erg/s) clusters branch into two populations --- radio halos lie on the correlation, while clusters with upper limits to radio-halo emission are well below that correlation. This bimodality remains if we excise cool cores from the X-ray luminosities. Correlating with Planck data, we find that P_1.4 scales with the cluster integrated SZ signal within R_500 as P_1.4 Y_500^2.1, in line with previous findings. However, contrary to previous studies that were limited by incompleteness and small sample size, we find that the "SZ-luminous" Y_500 > 6x10^-5 Mpc^2 clusters show a bimodal behavior 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.

preprint2013arXiv

The Hot and Energetic Universe: The evolution of galaxy groups and clusters

Major astrophysical questions related to the formation and evolution of structures, and more specifically of galaxy groups and clusters, will still be open in the coming decade and beyond: what is the interplay of galaxy, supermassive black hole, and intergalactic gas evolution in the most massive objects in the Universe - galaxy groups and clusters? What are the processes driving the evolution of chemical enrichment of the hot diffuse gas in large-scale structures? How and when did the first galaxy groups in the Universe, massive enough to bind more than 10^7 K gas, form? Focussing on the period when groups and clusters assembled (0.5<z<2.5), we show that, due to the continuum and line emission of this hot intergalactic gas at X-ray wavelengths, Athena+, combining high sensitivity with excellent spectral and spatial resolution, will deliver breakthrough observations in view of the aforementioned issues. Indeed, the physical and chemical properties of the hot intra-cluster gas, and their evolution across time, are a key to understand the co-evolution of galaxy and supermassive black hole within their environments.

preprint2012arXiv

A divergence-cleaning scheme for cosmological SPMHD simulations

In magnetohydrodynamics (MHD), the magnetic field is evolved by the induction equation and coupled to the gas dynamics by the Lorentz force. We perform numerical smoothed particle magnetohydrodynamics (Spmhd) simulations and study the influence of a numerical magnetic divergence. For instabilities arising from divergence B related errors, we find the hyperbolic/parabolic cleaning scheme suggested by Dedner et al. 2002 to give good results and prevent numerical artifacts from growing. Additionally, we demonstrate that certain current Spmhd implementations of magnetic field regularizations give rise to unphysical instabilities in long-time simulations. We also find this effect when employing Euler potentials (divergenceless by definition), which are not able to follow the winding-up process of magnetic field lines properly. Furthermore, we present cosmological simulations of galaxy cluster formation at extremely high resolution including the evolution of magnetic fields. We show synthetic Faraday rotation maps and derive structure functions to compare them with observations. Comparing all the simulations with and without divergence cleaning, we are able to confirm the results of previous simulations performed with the standard implementation of MHD in Spmhd at normal resolution. However, at extremely high resolution, a cleaning scheme is needed to prevent the growth of numerical errors at small scales.

preprint2012arXiv

Lensing and X-ray mass estimates of clusters (SIMULATION)

[Abridged] We present a comparison between weak-lensing (WL) and X-ray mass estimates of a sample of numerically simulated clusters. The sample consists on the 20 most massive objects at redshift z=0.25 and Mvir > 5 x 10^{14} Msun h^{-1}. They were found in a cosmological simulation of volume 1 h^{-3} Gpc^3, evolved in the framework of a WMAP-7 normalized cosmology. Each cluster has been resimulated at higher resolution and with more complex gas physics. We processed it thought Skylens and X-MAS to generate optical and X-ray mock observations along three orthogonal projections. The optical simulations include lensing effects on background sources. Standard observational tools and methods of analysis are used to recover the mass profiles of each cluster projection from the mock catalogues. Given the size of our sample, we could also investigate the dependence of the results on cluster morphology, environment, temperature inhomogeneity, and mass. We confirm previous results showing that WL masses obtained from the fit of the cluster tangential shear profiles with NFW functionals are biased low by ~ 5-10% with a large scatter (~10-25%). We show that scatter could be reduced by optimally selecting clusters either having regular morphology or living in substructure-poor environment. The X-ray masses are biased low by a large amount (~25-35%), evidencing the presence of both non-thermal sources of pressure in the ICM and temperature inhomogeneity, but they show a significantly lower scatter than weak-lensing-derived masses. The X-ray mass bias grows from the inner to the outer regions of the clusters. We find that both biases are weakly correlated with the third-order power ratio, while a stronger correlation exists with the centroid shift. Finally, the X-ray bias is strongly connected with temperature inhomogeneities.

preprint2012arXiv

Planck Early Results. II. The thermal performance of Planck

The performance of the Planck instruments in space is enabled by their low operating temperatures, 20K for LFI and 0.1K for HFI, achieved through a combination of passive radiative cooling and three active mechanical coolers. The scientific requirement for very broad frequency coverage led to two detector technologies with widely different temperature and cooling needs. Active coolers could satisfy these needs; a helium cryostat, as used by previous cryogenic space missions (IRAS, COBE, ISO, Spitzer, AKARI), could not. Radiative cooling is provided by three V-groove radiators and a large telescope baffle. The active coolers are a hydrogen sorption cooler (<20K), a 4He Joule-Thomson cooler (4.7K), and a 3He-4He dilution cooler (1.4K and 0.1K). The flight system was at ambient temperature at launch and cooled in space to operating conditions. The HFI bolometer plate reached 93mK on 3 July 2009, 50 days after launch. The solar panel always faces the Sun, shadowing the rest of Planck, andoperates at a mean temperature of 384K. At the other end of the spacecraft, the telescope baffle operates at 42.3K and the telescope primary mirror operates at 35.9K. The temperatures of key parts of the instruments are stabilized by both active and passive methods. Temperature fluctuations are driven by changes in the distance from the Sun, sorption cooler cycling and fluctuations in gas-liquid flow, and fluctuations in cosmic ray flux on the dilution and bolometer plates. These fluctuations do not compromise the science data.

preprint2012arXiv

Planck Intermediate Results II: Comparison of Sunyaev-Zeldovich measurements from Planck and from the Arcminute Microkelvin Imager for 11 galaxy clusters

A comparison is presented of Sunyaev-Zeldovich measurements for 11 galaxy clusters as obtained by Planck and by the ground-based interferometer, the Arcminute Microkelvin Imager. Assuming a universal spherically-symmetric Generalised Navarro, Frenk & White (GNFW) model for the cluster gas pressure profile, we jointly constrain the integrated Compton-Y parameter (Y_500) and the scale radius (theta_500) of each cluster. Our resulting constraints in the Y_500-theta_500 2D parameter space derived from the two instruments overlap significantly for eight of the clusters, although, overall, there is a tendency for AMI to find the Sunyaev-Zeldovich signal to be smaller in angular size and fainter than Planck. Significant discrepancies exist for the three remaining clusters in the sample, namely A1413, A1914, and the newly-discovered Planck cluster PLCKESZ G139.59+24.18. The robustness of the analysis of both the Planck and AMI data is demonstrated through the use of detailed simulations, which also discount confusion from residual point (radio) sources and from diffuse astrophysical foregrounds as possible explanations for the discrepancies found. For a subset of our cluster sample, we have investigated the dependence of our results on the assumed pressure profile by repeating the analysis adopting the best-fitting GNFW profile shape which best matches X-ray observations. Adopting the best-fitting profile shape from the X-ray data does not, in general, resolve the discrepancies found in this subset of five clusters. Though based on a small sample, our results suggest that the adopted GNFW model may not be sufficiently flexible to describe clusters universally.

preprint2012arXiv

Planck Intermediate Results. I. Further validation of new Planck clusters with XMM-Newton

We present further results from the ongoing XMM-Newton validation follow-up of Planck cluster candidates, detailing X-ray observations of eleven candidates detected at a signal-to-noise ratio of 4.5<S/N<5.3 in the same 10-month survey maps used in the construction of the Early SZ sample. The sample was selected in order to test internal SZ quality flags, and the pertinence of these flags is discussed in light of the validation results. Ten of the candidates are found to be bona fide clusters lying below the RASS flux limit. Redshift estimates are available for all confirmed systems via X-ray Fe-line spectroscopy. They lie in the redshift range 0.19<z<0.94, demonstrating Planck's capability to detect clusters up to high z. The X-ray properties of the new clusters appear to be similar to previous new detections by Planck at lower z and higher SZ flux: the majority are X-ray underluminous for their mass, estimated using Y_X as mass proxy, and many have a disturbed morphology. We find tentative indication for Malmquist bias in the Y_SZ-Y_X relation, with a turnover at Y_SZ \sim 4 e-4 arcmin^2. We present additional new optical redshift determinations with ENO and ESO telescopes of candidates previously confirmed with XMM-Newton. The X-ray and optical redshifts for a total of 20 clusters are found to be in excellent agreement. We also show that useful lower limits can be put on cluster redshifts using X-ray data only via the use of the Y_X vs. Y_SZ and X-ray flux F_X vs. Y_SZ relations.

preprint2012arXiv

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

We examine the relation between the galaxy cluster mass M and Sunyaev-Zeldovich (SZ) effect signal D_A^2 Y 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 M_WL-D_A^2 Y 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 R_500 are on average ~ 20 per cent larger than the corresponding weak lensing masses, at odds with 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, the mass discrepancy and difference in mass concentration is 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.

preprint2012arXiv

Planck Intermediate Results. IV. The XMM-Newton validation programme for new Planck galaxy clusters

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. 14 new clusters were detected by XMM, including 2 double systems. Redshifts lie in the range 0.2 to 0.9, with 6 clusters at z>0.5. Estimated M500 range from 2.5 10^14 to 8 10^14 Msun. We discuss our results in the context of the full XMM validation programme, in which 51 new clusters have been detected. This includes 4 double and 2 triple systems, some of which are chance projections on the sky of clusters at different z. We find that association with a RASS-BSC source is a robust indicator of the reliability of a candidate, whereas association with a FSC source 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 sigma being a good indication that the candidate is a real cluster. 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 level. The corresponding mass threshold depends on z. Systems with M500 > 5 10^14 Msun at z > 0.5 are easily detectable with Planck. The newly-detected clusters follow the YX-Y500 relation derived from X-ray selected samples. Compared to X-ray selected clusters, the new SZ clusters have a lower X-ray luminosity on average for their mass. There is no indication of departure from standard self-similar evolution in the X-ray versus SZ scaling properties. (abridged)

preprint2012arXiv

Planck Intermediate Results. V. Pressure profiles of galaxy clusters from the Sunyaev-Zeldovich effect

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 3 x R500, 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-Newton 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 slightly flatter than most predictions from numerical simulations. 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 3 x R500. Within the temperature driven uncertainties, our Planck constraints are compatible with the cosmic baryon fraction and expected gas fraction in halos.

preprint2012arXiv

Planck intermediate results. VI: The dynamical structure of PLCKG214.6+37.0, a Planck discovered triple system of galaxy clusters

The survey of galaxy clusters performed by Planck through the Sunyaev-Zeldovich effect has already discovered many interesting objects, thanks to the whole coverage of the sky. One of the SZ candidates detected in the early months of the mission near to the signal to noise threshold, PLCKG214.6+37.0, was later revealed by XMM-Newton to be a triple system of galaxy clusters. We have further investigated this puzzling system with a multi-wavelength approach and we present here the results from a deep XMM-Newton re-observation. The characterisation of the physical properties of the three components has allowed us to build a template model to extract the total SZ signal of this system with Planck data. We partly reconciled the discrepancy between the expected SZ signal from X-rays and the observed one, which are now consistent at less than 1.2 sigma. We measured the redshift of the three components with the iron lines in the X-ray spectrum, and confirmed that the three clumps are likely part of the same supercluster structure. The analysis of the dynamical state of the three components, as well as the absence of detectable excess X-ray emission, suggest that we are witnessing the formation of a massive cluster at an early phase of interaction.

preprint2012arXiv

Planck intermediate results. VIII. Filaments between interacting clusters

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. 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. 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 select physical pairs of clusters as candidates. Using the Planck data we construct a local map of the tSZ effect centered on each pair of galaxy clusters. ROSAT data is used to construct X-ray maps of these pairs. After having modelled and subtracted the tSZ effect and X-ray emission for each cluster in the pair we study the residuals on both the SZ and X-ray maps. 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)x10^-4, cm^-3. The Planck satellite mission has provided the first SZ detection of the hot and diffuse intercluster gas.

preprint2012arXiv

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

We present an analysis of Planck satellite data on the Coma Cluster observed via the Sunyaev-Zeldovich effect. Planck is able, for the first time, to detect SZ emission up to r ~ 3 X R_500. We test previously proposed models for the pressure distribution in clusters against the azimuthally averaged data. We find that the Arnaud et al. universal pressure profile does not fit Coma, and that their pressure profile for merging systems provides a good fit of the data only at r<R_500: by r=2XR_500 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>R_500 than the mean pressure profile predicted by the simulations. 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 Planck y profiles extracted from corresponding sectors we find pressure jumps of 4.5+0.4-0.2 and 5.0+1.3-0.1 in the west and southeast, respectively. Assuming Rankine-Hugoniot pressure jump conditions, we deduce that the shock waves should propagate with Mach number M_w=2.03+0.09-0.04 and M_se=2.05+0.25-0.02 in the West and Southeast, 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.

preprint2012arXiv

Resolving infall caustics in dark matter halos

We have found that the phase-space of a dark matter particles assembling a galactic halo in cosmological N-body simulations has well defined fine grained structure. Recently accreted particles form distinctive velocity streams with high density contrast. For fixed observer position these streams lead to peaks in velocity distribution. Overall structure is close to that emerging in the secondary infall model.

preprint2012arXiv

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

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.

preprint2012arXiv

Toward Unbiased Galaxy Cluster Masses from Line of Sight Velocity Dispersions

We study the use of red sequence selected galaxy spectroscopy for unbiased estimation of galaxy cluster masses. We use the publicly available galaxy catalog produced using the semi-analytic model of De Lucia & Blaizot (2007) on the Millenium Simulation (Springel et al. 2005). We explore the impacts on selection using galaxy color, projected separation from the cluster center, and galaxy luminosity. We study the relationship between cluster mass and velocity dispersion and identify and characterize the following sources of bias and scatter: halo triaxiality, dynamical friction of red luminous galaxies and interlopers. We show that due to halo triaxiality the intrinsic scatter of estimated line of sight dynamical mass is about three times larger (30-40%) than the one estimated using the 3D velocity dispersion (~12%) and a small bias (~1%) is induced. We find evidence of increasing scatter as a function of redshift and provide a fitting formula to account for it. We characterize the amount of bias and scatter introduced by dynamical friction when using subsamples of red-luminous galaxies to estimate the velocity dispersion. We study the presence of interlopers in spectroscopic samples and their effect on the estimated cluster dynamical mass. 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 15% level in mass. We show that the selection related uncertainties are small by comparison, providing hope that with further improvements this systematic floor can be reduced.

preprint2011arXiv

A Comparison of Cosmological Codes: Properties of Thermal Gas and Shock Waves in Large Scale Structures

[...] We present results for the statistics of thermal gas and the shock wave properties for a large volume simulated with three different cosmological numerical codes: the Eulerian total variations diminishing code TVD, the Eulerian piecewise parabolic method-based code ENZO, and the Lagrangian smoothed-particle hydrodynamics code GADGET. Starting from a shared set of initial conditions, we present convergence tests for a cosmological volume of side-length 100 Mpc/h, studying in detail the morphological and statistical properties of the thermal gas as a function of mass and spatial resolution in all codes. By applying shock finding methods to each code, we measure the statistics of shock waves and the related cosmic ray acceleration efficiencies, within the sample of simulations and for the results of the different approaches. We discuss the regimes of uncertainties and disagreement among codes, with a particular focus on the results at the scale of galaxy clusters. Even if the bulk of thermal and shock properties are reasonably in agreement among the three codes, yet some significant differences exist (especially between Eulerian methods and smoothed particle hydrodynamics). In particular, we report: a) differences of huge factors (10-100) in the values of average gas density, temperature, entropy, Mach number and shock thermal energy flux in the most rarefied regions of the simulations between grid and SPH methods; b) the hint of an entropy core inside clusters simulated in grid codes; c) significantly different phase diagrams of shocked cells in grid codes compared to SPH; d) sizable differences in the morphologies of accretion shocks between grid and SPH methods.

preprint2011arXiv

A non-ideal MHD Gadget: Simulating massive galaxy clusters

Magnetic fields in the intra-cluster medium of galaxy clusters have been studied in the past years through different methods. In the next years the up-coming generation of radio telescopes is going to provide new data that have the potential of setting constraints on the properties of magnetic fields in galaxy clusters. Here we present zoomed-in simulations for a set of massive galaxy clusters (M_v > 10^15 M_sun/h). This is an ideal sample to study the evolution of magnetic field during the process of structure formation in detail. Turbulent motions of the gas within the ICM will manifest themselves in a macroscopic magnetic resistivity eta_m, which has to be taken explicitly into account, especially at scales below the resolution limit. We have adapted the MHD GADGET code by Dolag & Stasyszyn (2009) to include the treatment of the magnetic resistivity and for the first time we have included non-ideal MHD equations to better follow the evolution of the magnetic field within galaxy clusters. We investigate which value of the magnetic resistivity eta_m is required to match the magnetic field profile derived from radio observations. We find that a value of eta_m ~ 6 10^27 cm^2/s is necessary to recover the shape of the magnetic field profile inferred from radio observations of the Coma cluster. This value agrees well with the expected level of turbulent motions within the ICM at our resolution limit. The magnetic field profiles can be fitted by a beta-model like profile (Cavaliere & Fusco-Femiano 1976), with small dispersion of the parameters. We find also that that the temperature, density and entropy profiles of the clusters depend on the magnetic resistivity constant, having flatter profiles in the inner regions when the magnetic resistivity increases.

preprint2011arXiv

Galactic ménage à trois: Simulating magnetic fields in colliding galaxies

We present high resolution simulations of a multiple merger of three disk galaxies including the evolution of magnetic fields performed with the N-body/SPH code Gadget. For the first time, we embed the galaxies in a magnetized, low-density medium, thus modeling an ambient IGM. The simulations include radiative cooling and a model for star formation and supernova feedback. The progenitor disks have initial magnetic seed fields in the range of 10e-9 to 10e-6 G and the IGM has initial fields of 10e-12 to 10e-9 G. The simulations are compared to a run excluding magnetic fields. We show that the propagation of interaction-driven shocks depends significantly on the initial magnetic field strength. The shocks propagate faster in simulations with stronger initial field, suggesting that the shocks are supported by magnetic pressure. The Mach numbers of the shocks range from approximately M=1.5 for the non-magnetized case up to M=6 for the highest initial magnetization, resulting in higher temperatures of the shock heated IGM gas. The magnetic field in the system saturates rapidly after the mergers at ~ 10e-6 G within the galaxies and ~ 10e-8 G in the IGM independent of the initial value. These field strengths agree with observed values and correspond to the equipartition value of the magnetic pressure with the turbulent pressure in the system. We also present synthetic radio and polarization maps for different phases of the evolution showing that shocks driven by the interaction produce a high amount of polarized emission. These idealized simulations indicate that magnetic fields play an important role for the hydrodynamics of the IGM during galactic interactions. We also show that even weak seed fields are efficiently strengthened during multiple galactic mergers. This interaction driven amplification might have been a key process for the magnetization of the Universe.

preprint2011arXiv

Matter power spectra in dynamical-Dark Energy cosmologies

(abridged) We used a suite of numerical cosmological simulations in order to investigate the effect of gas cooling and star formation on the large scale matter distribution. The simulations follow the formation of cosmic structures in five different Dark Energy models: the fiducial $Λ$CDM cosmology and four models where the Dark Energy density is allowed to have a non-trivial redshift evolution. For each cosmology we have a control run with dark matter only, in order to allow a direct assessment of the impact of baryonic processes. We found that the power spectra of gas and stars, as well as the total matter power spectrum, are in qualitative agreement with the results of previous works in the framework of the fiducial model, although several quantitative differences exist. We used the halo model in order to investigate the backreaction of gas and stars on the dark matter distribution, finding that it is very well reproduced by increasing the average dark matter halo concentration by 17%, irrespective of the mass. Moving to model universes dominated by dynamical Dark Energy, it turns out that they introduce a specific signature on the power spectra of the various matter components, that is qualitatively independent of the exact cosmology considered. This generic shape is well captured by the halo model, however the finer details of the dark matter power spectrum can be precisely captured only at the cost of a few slight modifications to the ingredients entering the model. The backreaction of baryons onto the dark matter distribution works pretty much in the same way as in the reference $Λ$CDM model. Nonetheless, the increment in average concentration is less pronounced than in the fiducial model (only $\sim 10%$), in agreement with a series of other clues pointing toward the fact that star formation is less efficient when Dark Energy displays a dynamical evolution.

preprint2011arXiv

Planck Early Results VIII: The all-sky Early Sunyaev-Zeldovich cluster sample

We present the first all-sky sample of galaxy clusters detected blindly by the Planck satellite through the Sunyaev-Zeldovich (SZ) effect from its six highest frequencies. This early SZ (ESZ) sample is comprised of 189 candidates, which have a high signal-to-noise ratio ranging from 6 to 29. Its high reliability (purity above 95%) is further ensured by an extensive validation process based on Planck internal quality assessments and by external cross-identification and follow-up observations. Planck provides the first measured SZ signal for about 80% of the 169 previously-known ESZ clusters. Planck furthermore releases 30 new cluster candidates, amongst which 20 meet the ESZ signal-to-noise selection criterion. At the submission date, twelve of the 20 ESZ candidates were confirmed as new clusters, with eleven confirmed using XMM-Newton snapshot observations, most of them with disturbed morphologies and low luminosities. The ESZ clusters are mostly at moderate redshifts (86% with z below 0.3) and span more than a decade in mass, up to the rarest and most massive clusters with masses above 10^15 Msol.

preprint2011arXiv

Planck Early Results XI: Calibration of the local galaxy cluster Sunyaev-Zeldovich scaling relations

We present precise Sunyaev-Zeldovich (SZ) effect measurements in the direction of 62 nearby galaxy clusters (z <0.5) detected at high signal-to-noise in the first Planck all-sky dataset. The sample spans approximately a decade in total mass, 10^14 < M_500 < 10^15, where M_500 is the mass corresponding to a total density contrast of 500. Combining these high quality Planck measurements with deep XMM-Newton X-ray data, we investigate the relations between D_A^2 Y_500, the integrated Compton parameter due to the SZ effect, and the X-ray-derived gas mass M_g,500, temperature T_X, luminosity L_X, SZ signal analogue Y_X,500 = M_g,500 * T_X, and total mass M_500. After correction for the effect of selection bias on the scaling relations, we find results that are in excellent agreement with both X-ray predictions and recently-published ground-based data derived from smaller samples. The present data yield an exceptionally robust, high-quality local reference, and illustrate Planck's unique capabilities for all-sky statistical studies of galaxy clusters.

preprint2011arXiv

Planck early results: Cluster Sunyaev-Zeldovich optical scaling relations

We present the Sunyaev-Zeldovich (SZ) signal-to-richness scaling relation (Y500-N200) for the MaxBCG cluster catalogue. Employing a multi-frequency matched filter on the Planck sky maps, we measure the SZ signal for each cluster by adapting the filter according to weak-lensing calibrated mass-richness relations (N200-M500). We bin our individual measurements and detect the SZ signal down to the lowest richness systems (N200=10) with high significance, achieving a detection of the SZ signal in systems with mass as low as M500~5e13 Msolar. The observed Y500-N200 relation is well modeled by a power law over the full richness range. It has a lower normalisation at given N200 than predicted based on X-ray models and published mass-richness relations. An X-ray subsample, however, does conform to the predicted scaling, and model predictions do reproduce the relation between our measured bin-average SZ signal and measured bin-average X-ray luminosities. At fixed richness, we find an intrinsic dispersion in the Y500-N200 relation of 60% rising to of order 100% at low richness. Thanks to its all-sky coverage, Planck provides observations for more than 13,000 MaxBCG clusters and an unprecedented SZ/optical data set, extending the list of known cluster scaling laws to include SZ-optical properties. The data set offers essential clues for models of galaxy formation. Moreover, the lower normalisation of the SZ-mass relation implied by the observed SZ-richness scaling has important consequences for cluster physics and cosmological studies with SZ clusters.

preprint2011arXiv

Planck early results: first assessment of the High Frequency Instrument in-flight performance

The Planck High Frequency Instrument (HFI) is designed to measure the temperature and polarization anisotropies of the Cosmic Microwave Background and galactic foregrounds in six wide bands centered at 100, 143, 217, 353, 545 and 857 GHz at an angular resolution of 10' (100 GHz), 7' (143 GHz), and 5' (217 GHz and higher). HFI has been operating flawlessly since launch on 14 May 2009. The bolometers cooled to 100 mK as planned. The settings of the readout electronics, such as the bolometer bias current, that optimize HFI's noise performance on orbit are nearly the same as the ones chosen during ground testing. Observations of Mars, Jupiter, and Saturn verified both the optical system and the time response of the detection chains. The optical beams are close to predictions from physical optics modeling. The time response of the detection chains is close to pre-launch measurements. The detectors suffer from an unexpected high flux of cosmic rays related to low solar activity. Due to the redundancy of Planck's observations strategy, the removal of a few percent of data contaminated by glitches does not affect significantly the sensitivity. The cosmic rays heat up significantly the bolometer plate and the modulation on periods of days to months of the heat load creates a common drift of all bolometer signals which do not affect the scientific capabilities. Only the high energy cosmic rays showers induce inhomogeneous heating which is a probable source of low frequency noise.

preprint2011arXiv

Planck Early Results: The Planck mission

The European Space Agency's Planck satellite was launched on 14 May 2009, and has been surveying the sky stably and continuously since 13 August 2009. Its performance is well in line with expectations, and it will continue to gather scientific data until the end of its cryogenic lifetime. We give an overview of the history of Planck in its first year of operations, and describe some of the key performance aspects of the satellite. This paper is part of a package submitted in conjunction with Planck's Early Release Compact Source Catalogue, the first data product based on Planck to be released publicly. The package describes the scientific performance of the Planck payload, and presents results on a variety of astrophysical topics related to the sources included in the Catalogue, as well as selected topics on diffuse emission.

preprint2011arXiv

Planck early results. IX. XMM-Newton follow-up for validation of Planck cluster candidates

We present the XMM-Newton follow-up for confirmation of Planck cluster candidates. Twenty-five candidates have been observed to date using snapshot (~10 ksec) exposures, ten as part of a pilot programme to sample a low range of signal-to-noise ratios (4<S/N<6), and a further 15 in a programme to observe a sample of S/N>5 candidates. The sensitivity and spatial resolution of XMM-Newton allows unambiguous discrimination between clusters and false candidates. The 4 false candidates have S/N <= 4.1. A total of 21 candidates are confirmed as extended X-ray sources. Seventeen are single clusters, the majority of which are found to have highly irregular and disturbed morphologies (about ~70%). The remaining four sources are multiple systems, including the unexpected discovery of a supercluster at z=0.45. For 20 sources we are able to derive a redshift estimate from the X-ray Fe K line (albeit of variable quality). The new clusters span the redshift range 0.09 <= z <= 0.54, with a median redshift of z~0.37. A first determination is made of their X-ray properties including the characteristic size, which is used to improve the estimate of the SZ Compton parameter, Y_SZ. The follow-up validation programme has helped to optimise the Planck candidate selection process. It has also provided a preview of the X-ray properties of these newly-discovered clusters, allowing comparison with their SZ properties, and to the X-ray and SZ properties of known clusters observed in the Planck survey. Our results suggest that Planck may have started to reveal a non-negligible population of massive dynamically perturbed objects that is under-represented in X-ray surveys. However, despite their particular properties, these new clusters appear to follow the Y_SZ-Y_X relation established for X-ray selected objects, where Y_X is the product of the gas mass and temperature.

preprint2011arXiv

Planck early results. VI. The High Frequency Instrument data processing

We describe the processing of the 336 billion raw data samples from the High Frequency Instrument (HFI) which we performed to produce six temperature maps from the first 295 days of Planck-HFI survey data. These maps provide an accurate rendition of the sky emission at 100, 143, 217, 353, 545 and 857 GHz with an angular resolution ranging from 9.9 to 4.4^2. The white noise level is around 1.5 μK degree or less in the 3 main CMB channels (100--217GHz). The photometric accuracy is better than 2% at frequencies between 100 and 353 GHz and around 7% at the two highest frequencies. The maps created by the HFI Data Processing Centre reach our goals in terms of sensitivity, resolution, and photometric accuracy. They are already sufficiently accurate and well-characterised to allow scientific analyses which are presented in an accompanying series of early papers. At this stage, HFI data appears to be of high quality and we expect that with further refinements of the data processing we should be able to achieve, or exceed, the science goals of the Planck project.

preprint2011arXiv

Planck Early Results. VII. The Early Release Compact Source Catalog

(Abridged) A brief description of the methodology of construction, contents and usage of the Planck Early Release Compact Source Catalogue (ERCSC), including the Early Cold Cores (ECC) and the Early Sunyaev-Zeldovich (ESZ) cluster catalogue is provided. The catalogue is based on data that consist of mapping the entire sky once and 60% of the sky a second time by Planck, thereby comprising the first high sensitivity radio/submillimetre observations of the entire sky. A Monte-Carlo algorithm based on the injection and extraction of artificial sources into the Planck maps was implemented to select reliable sources among all extracted candidates such that the cumulative reliability of the catalogue is >=90%. The 10sigma photometric flux density limit of the catalogue at |b|>30 deg is 0.49, 1.0, 0.67, 0.5, 0.33, 0.28, 0.25, 0.47 and 0.82 Jy at each of the nine frequencies between 30 and 857 GHz. Sources which are up to a factor of ~2 fainter than this limit, and which are present in "clean" regions of the Galaxy where the sky background due to emission from the interstellar medium is low, are included in the ERCSC if they meet the high reliability criterion. The Planck ERCSC sources have known associations to stars with dust shells, stellar cores, radio galaxies, blazars, infrared luminous galaxies and Galactic interstellar medium features. A significant fraction of unclassified sources are also present in the catalogs. In addition, two early release catalogs that contain 915 cold molecular cloud core candidates and 189 SZ cluster candidates that have been generated using multi-frequency algorithms are presented. The entire source list, with more than 15000 unique sources, is ripe for follow-up characterisation with Herschel, ATCA, VLA, SOFIA, ALMA and other ground-based observing facilities.

preprint2011arXiv

Planck Early Results. X. Statistical analysis of Sunyaev-Zeldovich scaling relations for X-ray galaxy clusters

All-sky data from the Planck survey and the Meta-Catalogue of X-ray detected Clusters of galaxies (MCXC) are combined to investigate the relationship between the thermal Sunyaev-Zeldovich (SZ) signal and X-ray luminosity. The sample comprises ~ 1600 X-ray clusters with redshifts up to ~ 1 and spans a wide range in X-ray luminosity. The SZ signal is extracted for each object individually, and the statistical significance of the measurement is maximised by averaging the SZ signal in bins of X-ray luminosity, total mass, or redshift. The SZ signal is detected at very high significance over more than two decades in X-ray luminosity (10^43 erg/s < L_500 E(z)^-7/3 < 2 X 10^45 erg/s). The relation between intrinsic SZ signal and X-ray luminosity is investigated and the measured SZ signal is compared to values predicted from X-ray data. Planck measurements and X-ray based predictions are found to be in excellent agreement over the whole explored luminosity range. No significant deviation from standard evolution of the scaling relations is detected. For the first time the intrinsic scatter in the scaling relation between SZ signal and X-ray luminosity is measured and found to be consistent with the one in the luminosity -- mass relation from X-ray studies. There is no evidence of any deficit in SZ signal strength in Planck data relative to expectations from the X-ray properties of clusters, underlining the robustness and consistency of our overall view of intra-cluster medium properties.

preprint2011arXiv

Pointing to the minimum scatter: the generalized scaling relations for galaxy clusters

We introduce a generalized scaling law, M_tot = 10^K A^a B^b, to look for the minimum scatter in reconstructing the total mass of hydrodynamically simulated X-ray galaxy clusters, given gas mass M_gas, luminosity L and temperature T. We find a locus in the plane of the logarithmic slopes $a$ and $b$ of the scaling relations where the scatter in mass is minimized. This locus corresponds to b_M = -3/2 a_M +3/2 and b_L = -2 a_L +3/2 for A=M_gas and L, respectively, and B=T. Along these axes, all the known scaling relations can be identified (at different levels of scatter), plus a new one defined as M_tot ~ (LT)^(1/2). Simple formula to evaluate the expected evolution with redshift in the self-similar scenario are provided. In this scenario, no evolution of the scaling relations is predicted for the cases (b_M=0, a_M=1) and (b_L=7/2, a_L=-1), respectively. Once the single quantities are normalized to the average values of the sample under considerations, the normalizations K corresponding to the region with minimum scatter are very close to zero. The combination of these relations allows to reduce the number of free parameters of the fitting function that relates X-ray observables to the total mass and includes the self-similar redshift evolution.

preprint2011arXiv

Resonant scattering in galaxy clusters for anisotropic gas motions on various spatial scales

The determination of characteristic amplitudes and anisotropy of hot gas motions in galaxy clusters from observations of the brightest resonance lines is discussed. Gas motions affect (i) the spectral line shape through the Doppler effect and (ii) the radial surface brightness profiles in lines during resonant scattering. Radiative transfer calculations have been performed by the Monte Carlo method in the FeXXV resonance line at 6.7 eV for the Perseus cluster (Abell426). It was shown that (a) radial motions reduce the scattering efficiency much more dramatically than purely tangential motions; (b) large-scale gas motions weakly affect the scattering efficiency. The uncertainty in measuring the characteristics of gas motions using resonant scattering has been estimated for existing and future observations of clusters.

preprint2011arXiv

The effect of primordial non-Gaussianity on the skeleton of cosmic shear maps

(abridged) We explore the imprints of deviations from Gaussian primordial density fluctuations on the skeleton of the large-scale matter distribution as mapped through cosmological weak lensing. We computed the skeleton length of simulated effective convergence maps covering $\sim 35$ sq. deg each, extracted from a suite of cosmological $n-$body runs with different levels of local primordial non-Gaussianity. The latter is expected to alter the structure formation process with respect to the fiducial Gaussian scenario, and thus to leave a signature on the cosmic web. We found that alterations of the initial conditions consistently modify both the cumulative and the differential skeleton length, although the effect is generically smaller than the cosmic variance and depends on the smoothing of the map prior to the skeleton computation. Nevertheless, the qualitative shape of these deviations is rather similar to their primordial counterparts, implying that skeleton statistics retain good memory of the initial conditions. We performed a statistical analysis in order to find out at what Confidence Level primordial non-Gaussianity could be constrained by the skeleton test on cosmic shear maps of the size we adopted. At 68.3% Confidence Level we found an error on the measured level of primordial non-Gaussianity of $Δf_\mathrm{NL}\sim 300$, while at 90% Confidence Level it is of $Δf_\mathrm{NL}\sim 500$. While these values by themselves are not competitive with the current constraints, weak lensing maps larger than those used here would have a smaller field-to-field variance, and thus would likely lead to tighter constraints. A rough estimate indicates $Δf_\mathrm{NL} \sim$ a few tens at 68.3% Confidence Level for an all-sky weak lensing survey.

preprint2011arXiv

Why are central radio relics so rare?

In this paper we address the question why cluster radio relics that are connected to shock acceleration, so-called radio gischt, have preferentially been found in the outskirts of galaxy clusters. By identifying merger shock waves in cosmological grid simulations, we explore several prescriptions for relating the energy dissipated in shocks to the energy emitted in the radio band. None of the investigated models produce detectable radio relics within 100-200 kpc from the cluster centre. All models cause > 50 per cent of the detectable relic emission at projected distances > 800 kpc. Central radio relics caused by shocks that propagate along the line-of-sight are rare events for simple geometrical reasons, and they have a low surface brightness making them elusive for current instruments. Our simulations show that the radial distribution of observed relics can be explained by the radial trend of dissipated kinetic energy in shocks, that increases with distance from the cluster centre up until half of the virial radius.

preprint2011arXiv

X-ray mass proxies from hydrodynamic simulations of galaxy clusters (paper I)

We present a detailed study of scaling relations between total cluster mass and three mass proxies based on X-ray observables: temperature of the intra-cluster medium, gas mass and the product of the two, Y_X. Our analysis is based on two sets of high-resolution hydrodynamical simulations performed with the TreePM-SPH GADGET code. The first set includes about 140 clusters with masses above 5x10^13 M_sun/h (30 having mass above 10^15 M_sun/h), that have been simulated with (i) non-radiative physics and including (ii) cooling, star formation, chemical enrichment and the effect of supernova feedback triggering galactic ejecta. This large statistics is used to quantify the robustness of the scaling relations, to determine their redshift evolution and to calibrate their intrinsic scatter and its distribution. We use a smaller set of clusters including 18 halos with masses above 5x10^13 M_sun/h to test the robustness of mass proxies against changing the physical processes included in simulations (thermal conduction, artificial viscosity, cooling and star formation, galactic winds and AGN feedback). We find the M-Y_X scaling relation to be the least sensitive one to variations of the ICM physics, with its slope and redshift evolution close to the self-similar model predictions. The distribution of the scatter around the best-fitting relations is close to a log-normal one. M_gas has the smallest scatter in mass, with values of sigma_lnM = 0.04-0.06, depending on the physics included in the simulation, and with a mild dependence on redshift. The M-T relation is the one with the largest scatter, with sigma_lnM > 0.1 at z=0, increasing to > 0.15 at z=1. The intrinsic scatter in the M-Y_X relation is slightly larger than in the M-M_gas relation. These results confirm that both Y_X and M_gas mass proxies are well suited for cosmological applications of future large X-ray surveys. [abridged]

preprint2010arXiv

Dynamical difference between the cD galaxy and the stellar diffuse component in simulated galaxy clusters

Member galaxies within galaxy clusters nowadays can be routinely identified in cosmological, hydrodynamical simulations using methods based on identifying self bound, locally over dense substructures. However, distinguishing the central galaxy from the stellar diffuse component within clusters is notoriously difficult, and in the center it is not even clear if two distinct stellar populations exist. Here, after subtracting all member galaxies, we use the velocity distribution of the remaining stars and detect two dynamically, well-distinct stellar components within simulated galaxy clusters. These differences in the dynamics can be used to apply an un-binding procedure which leads to a spatial separation of the two components into a cD and a diffuse stellar component (DSC). Applying our new algorithm to a cosmological, hydrodynamical simulation we find that -- in line with previous studies -- these two components have clearly distinguished spatial and velocity distributions as well as different star formation histories. We show that the DSC fraction -- which can broadly be associated with the observed intra cluster light -- does not depend on the virial mass of the galaxy cluster and is much more sensitive to the formation history of the cluster. We conclude that the separation of the cD and the DSC in simulations, based on our dynamical criteria, is more physically motivated than current methods which depend on implicit assumptions on a length scale associated with the cD galaxy and therefore represent a step forward in understanding the different stellar components within galaxy clusters. Our results also show the importance of analyzing the dynamics of the DSC to characterize its properties and understand its origin.

preprint2010arXiv

Gas cooling in semi-analytic models and SPH simulations: are results consistent?

We present a detailed comparison between the galaxy populations within a massive cluster, as predicted by hydrodynamical SPH simulations and by a semi-analytic model (SAM) of galaxy formation. Both models include gas cooling and a simple prescription of star formation, which consists in transforming instantaneously any cold gas available into stars, while neglecting any source of energy feedback. We find that, in general, galaxy populations from SAMs and SPH have similar statistical properties, in agreement with previous studies. However, when comparing galaxies on an object-by-object basis, we find a number of interesting differences: a) the star formation histories of the brightest cluster galaxies (BCGs) from SAM and SPH models differ significantly, with the SPH BCG exhibiting a lower level of star formation activity at low redshift, and a more intense and shorter initial burst of star formation with respect to its SAM counterpart; b) while all stars associated with the BCG were formed in its progenitors in the semi-analytic model used here, this holds true only for half of the final BCG stellar mass in the SPH simulation, the remaining half being contributed by tidal stripping of stars from the diffuse stellar component associated with galaxies accreted on the cluster halo; c) SPH satellites can loose up to 90 per cent of their stellar mass at the time of accretion, due to tidal stripping, a process not included in the semi-analytic model used in this study; d) in the SPH simulation, significant cooling occurs on the most massive satellite galaxies and this lasts for up to 1 Gyr after accretion. This physical process is not included in the semi-analytic model used in our study, as well as in most of the models discussed in the recent literature.

preprint2010arXiv

Lower limit on the strength and filling factor of extragalactic magnetic fields

High energy photons from blazars can initiate electromagnetic pair cascades interacting with the extragalactic photon background. The charged component of such cascades is deflected and delayed by extragalactic magnetic fields (EGMF), reducing thereby the observed point-like flux and leading potentially to multi degree images in the GeV energy range. We calculate the fluence of 1ES 0229+200 as seen by Fermi-LAT for different EGMF profiles using a Monte Carlo simulation for the cascade development. The non-observation of 1ES 0229+200 by Fermi-LAT suggests that the EGMF fills at least 60% of space with fields stronger than {\cal O}(10^{-16}-10^{-15})G for life times of TeV activity of {\cal O}(10^2-10^4)yr. Thus the (non-) observation of GeV extensions around TeV blazars probes the EGMF in voids and puts strong constraints on the origin of EGMFs: Either EGMFs were generated in a space filling manner (e.g. primordially) or EGMFs produced locally (e.g. by galaxies) have to be efficiently transported to fill a significant volume fraction, as e.g. by galactic outflows.

preprint2010arXiv

Measuring cosmic magnetic fields by rotation measure-galaxy cross-correlations in cosmological simulations

Using cosmological MHD simulations of the magnetic field in galaxy clusters and filaments we evaluate the possibility to infer the magnetic field strength in filaments by measuring cross-correlation functions between Faraday Rotation Measures (RM) and the galaxy density field. We also test the reliability of recent estimates considering the problem of data quality and Galactic foreground (GF) removal in current datasets. Besides the two self-consistent simulations of cosmological magnetic fields based on primordial seed fields and galactic outflows analyzed here, we also explore a larger range of models scaling up the resulting magnetic fields of one of the simulations. We find that, if an unnormalized estimator for the cross-correlation functions and a GF removal procedure is used, the detectability of the cosmological signal is only possible for future instruments (e.g. SKA and ASKAP). However, mapping of the observed RM signal to the underlying magnetization of the Universe (both in space and time) is an extremely challenging task which is limited by the ambiguities of our model parameters, as well as to the weak response of the RM signal in low density environments. Therefore, we conclude that current data cannot constrain the amplitude and distribution of magnetic fields within the large scale structure and a detailed theoretical understanding of the build up and distribution of magnetic fields within the Universe will be needed for the interpretation of future observations.

preprint2010arXiv

Photometric and clustering properties of hydrodynamical galaxies in a cosmological volume: results at z=0

In this work, we present results for the photometric and clustering properties of galaxies that arise in a LambdaCDM hydrodynamical simulation of the local universe. The present-day distribution of matter was constructed to match the observed large scale pattern of the IRAS 1.2-Jy galaxy survey. Our simulation follows the formation and evolution of galaxies in a cosmological sphere with a volume of ~130^3 (Mpc/h)^3 including supernova feedback, galactic winds, photoheating due to an uniform meta-galactic background and chemical enrichment of the gas and stellar populations. However, we do not consider AGNs. In the simulation, a total of ~20000 galaxies are formed above the resolution limit, and around 60 haloes are more massive than ~10^14 M_sun. Luminosities of the galaxies are calculated based on a stellar population synthesis model including the attenuation by dust, which is calculated from the cold gas left within the simulated galaxies. Environmental effects like colour bi-modality and differential clustering power of the hydrodynamical galaxies are qualitatively similar to observed trends. Nevertheless, the overcooling present in the simulations lead to too blue and overluminous brightest cluster galaxies (BCGs). To overcome this, we mimic the late-time suppression of star formation in massive halos by ignoring recently formed stars with the aid of a simple post-processing recipe. In this way we find luminosity functions, both for field and group/cluster galaxies, in better agreement with observations. Specifically, the BCGs then follow the observed luminosity-halo mass relation. However, in such a case, the colour bi-modality is basically lost, pointing towards a more complex interplay of late suppression of star formation than what is given by the simple scheme adopted.

preprint2010arXiv

Polarization of X-ray lines from galaxy clusters and elliptical galaxies - a way to measure tangential component of gas velocity

We study the impact of gas motions on the polarization of bright X-ray emission lines from the hot intercluster medium (ICM). The polarization naturally arises from resonant scattering of emission lines owing to a quadrupole component in the radiation field produced by a centrally peaked gas density distribution. If differential gas motions are present then a photon emitted in one region of the cluster will be scattered in another region only if their relative velocities are small enough and the Doppler shift of the photon energy does not exceed the line width. This affects both the degree and the direction of polarization. The changes in the polarization signal are in particular sensitive to the gas motions perpendicular to the line of sight. We calculate the expected degree of polarization for several patterns of gas motions, including a slow inflow expected in a simple cooling flow model and a fast outflow in an expanding spherical shock wave. In both cases, the effect of non-zero gas velocities is found to be minor. We also calculate the polarization signal for a set of clusters, taken from large-scale structure simulations and evaluate the impact of the gas bulk motions on the polarization signal. We argue that the expected degree of polarization is within reach of the next generation of space X-ray polarimeters.

preprint2010arXiv

Radio Halos From Simulations And Hadronic Models II: The Scaling Relations of Radio Halos

We use results from a constrained, cosmological MHD simulation of the Local Universe to predict radio halos and their evolution for a volume limited set of galaxy clusters and compare to current observations. The simulated magnetic field inside the clusters is a result of turbulent amplification within them, with the magnetic seed originating from star-burst driven, galactic outflows. We evaluate three models, where we choose different normalizations for the Cosmic Ray proton population within clusters. Similar to our previous analysis of the Coma cluster (Donnert et al. 2010), the radial profile and the morphological properties of observed radio halos can not be reproduced, even with a radially increasing energy fraction within the cosmic ray proton population. Scaling relations between X-ray luminosity and radio power can be reproduced by all models, however all models fail in the prediction of clusters with no radio emission. Also the evolutionary tracks of our largest clusters in all models fail to reproduce the observed bi-modality in radio luminosity. This provides additional evidence that the framework of hadronic, secondary models is disfavored to reproduce the large scale diffuse radio emission of galaxy clusters. We also provide predictions for the unavoidable emission of $γ$-rays from the hadronic models for the full cluster set. None of such secondary models is yet excluded by the observed limits in $γ$-ray emission, emphasizing that large scale diffuse radio emission is a powerful tool to constrain the amount of cosmic ray protons in galaxy clusters.

preprint2010arXiv

Rotation Measures of Radio Sources in Hot Galaxy Clusters

The goal of this work is to investigate the Faraday rotation measure (RM) of radio galaxies in hot galaxy clusters in order to establish a possible connection between the magnetic field strength and the gas temperature of the intracluster medium. We performed Very Large Array observations at 3.6 cm and 6 cm of two radio galaxies located in A401 and Ophiuchus, a radio galaxy in A2142, and a radio galaxy located in the background of A2065. All these galaxy clusters are characterized by high temperatures. We obtained detailed RM images at an angular resolution of 3'' for most of the observed radio galaxies. The RM images are patchy and reveal fine substructures of a few kpc in size. Under the assumption that the radio galaxies themselves have no effect on the measured RMs, these structures indicate that the intracluster magnetic fields fluctuate down to such small scales. These new data are compared with RM information present in the literature for cooler galaxy clusters. For a fixed projected distance from the cluster center, clusters with higher temperature show a higher dispersion of the RM distributions (sigmaRM), mostly because of the higher gas density in these clusters. Although the previously known relation between the clusters X-ray surface brightness (Sx) at the radio galaxy location and sigmaRM is confirmed, a possible connection between the sigmaRM-Sx relation and the cluster temperature, if present, is very weak. Therefore, in view of the current data, it is impossible to establish a strict link between the magnetic field strength and the gas temperature of the intracluster medium.

preprint2010arXiv

Simulating magnetic fields in the Antennae galaxies

We present self-consistent high-resolution simulations of NGC4038/4039 (the "Antennae galaxies") including star formation, supernova feedback and magnetic fields performed with the N-body/SPH code Gadget, in which magnetohydrodynamics are followed with the SPH method. We vary the initial magnetic field in the progenitor disks from 1 nG to 100 muG. At the time of the best match with the central region of the Antennae system the magnetic field has been amplified by compression and shear flows to an equilibrium field of approximately 10 muG, independent of the initial seed field. These simulations are a proof of the principle that galaxy mergers are efficient drivers for the cosmic evolution of magnetic fields. We present a detailed analysis of the magnetic field structure in the central overlap region. Simulated radio and polarization maps are in good morphological and quantitative agreement with the observations. In particular, the two cores with the highest synchrotron intensity and ridges of regular magnetic fields between the cores and at the root of the southern tidal arm develop naturally in our simulations. This indicates that the simulations are capable of realistically following the evolution of the magnetic fields in a highly non-linear environment. We also discuss the relevance of the amplification effect for present day magnetic fields in the context of hierarchical structure formation.

preprint2010arXiv

Statistical properties of SZ and X-ray cluster detections

Aims. We calibrate the number density, completeness, reliability and the lower mass limit of galaxy-cluster detections through their thermal SZ signal, and compare them to X-ray cluster detections. Methods. We simulate maps of the thermal SZ effect and the X-ray emission from light cones constructed in a large, hydrodynamical, cosmological simulation volume, including realistic noise contributions. The maps are convolved with linear, optimised, single- and multi-band filters to identify local peaks and their signal-to-noise ratios. The resulting peak catalogues are then compared to the halo population in the simulation volume to identify true and spurious detections. Results. Multi-band filtering improves the statistics of SZ cluster detections considerably compared to single-band filtering. Observations with the characteristics of ACT detect clusters with masses M>6-9e13 M_o/h, quite independent of redshift, reach 50% completeness at ~1e14 M_o/h and 100% completeness at ~2e14 M_o/h. Samples are contaminated by only a few per cent spurious detections. This is broadly comparable to X-ray cluster detections with XMM-Newton with 100 ks exposure time in the soft band, except that the mass limit for X-ray detections increases much more steeply with redshift than for SZ detections. A comparison of true and filtered signals in the SZ and X-ray maps confirms that the filters introduce at most a negligible bias.

preprint2010arXiv

The Coma cluster magnetic field from Faraday rotation measures

The aim of the present work is to constrain the Coma cluster magnetic field strength, its radial profile and power spectrum by comparing Faraday Rotation Measure (RM) images with numerical simulations of the magnetic field. We have analyzed polarization data for seven radio sources in the Coma cluster field observed with the Very Large Array at 3.6, 6 and 20 cm, and derived Faraday Rotation Measures with kiloparsec scale resolution. Random three dimensional magnetic field models have been simulated for various values of the central intensity B_0 and radial power-law slope eta, where eta indicates how the field scales with respect to the gas density profile. We derive the central magnetic field strength, and radial profile values that best reproduce the RM observations. We find that the magnetic field power spectrum is well represented by a Kolmogorov power spectrum with minimum scale ~ 2 kpc and maximum scale ~ 34 kpc. The central magnetic field strength and radial slope are constrained to be in the range (B_0=3.9 microG; eta=0.4) and (B_0=5.4 microG; eta=0.7) within 1sigma. The best agreement between observations and simulations is achieved for B_0=4.7 microG; eta=0.5. Values of B_0>7 microG and <3 microG as well as eta < 0.2 and eta > 1.0 are incompatible with RM data at 99 % confidence level.

preprint2010arXiv

The scatter in the radial profiles of X-ray luminous galaxy clusters as diagnostic of the thermodynamical state of the ICM

We study the azimuthal scatter in the radial profiles of X-ray luminous galaxy clusters, with two sets of high-resolution cosmological re-simulations obtained with the codes ENZO and GADGET2. The average gas profiles are computed for different angular sectors of the cluster projected volume, and compared with the mean cluster profiles at each radius from the center. We report that in general the level of azimuthal scatter is found to be about 10 per cent for gas density, temperature and entropy inside R200, and about 25 per cent for X-ray luminosity for the same volume. These values generally doubles going to 2 R200 from the cluster center, and are generally found to be higher (by about 20-40 percent) in the case of perturbed systems. A comparison with results from recent SUZAKU observations is discussed, showing the possibility to simply interpret the large azimuthal scatter of observables in light of our simulated results.

preprint2010arXiv

The transition from population III to population II-I star formation

We present results from the first cosmological simulations which study the onset of primordial, metal-free (population III), cosmic star formation and the transition to the present-day, metal-rich star formation (population II-I), including molecular (H$_2$, HD, etc.) evolution, tracing the injection of metals by supernovæ into the surrounding intergalactic medium and following the change in the initial stellar mass function (IMF) according to the metallicity of the corresponding stellar population. Our investigation addresses the role of a wide variety of parameters (critical metallicity for the transition, IMF slope and range, SN/pair-instability SN metal yields, star formation threshold, resolution, etc.) on the metal-enrichment history and the associated transition in the star formation mode. All simulations present common trends. Metal enrichment is very patchy, with rare, unpolluted regions surviving at all redshifts, inducing the simultaneous presence of metal-free and metal-rich star formation regimes. As a result of the rapid pollution within high-density regions due to the first SN/pair-instability SN, local metallicity is quickly boosted above the critical metallicity for the transition. The population III regime lasts for a very short period during the first stages of star formation ($\sim 10^7\,\rm yr$), and its average contribution to the total star formation rate density drops rapidly below $\sim 10^{-3}-10^{-2}$.

preprint2009arXiv

An MHD Gadget for cosmological simulations

Various radio observations have showed that the hot atmospheres of galaxy clusters are magnetized. However, our understanding of the origin of these magnetic fields, their implications on structure formation and their interplay with the dynamics of the cluster atmosphere, especially in the centers of galaxy clusters is still very limited. In preparation to the upcoming new generation of radio telescopes (like EVLA, LWA, LOFAR and SKA), a huge effort is being made to learn more about cosmological magnetic fields from the observational perspective. Here we present the implementation of magneto hydrodynamics in the cosmological SPH code GADGET. We discuss the details of the implementation and various schemes to suppress numerical instabilities as well as regularization schemes, in the context of cosmological simulations. The performance of the SPH MHD code is demonstrated in various one and two dimensional test problems, which we performed if a fully, three dimensional setup to test the code under realistic circumstances. Comparing with solutions obtained with ATHENA, we find excellent agreement with our SPH MHD implementation. Finally we apply our SPH MHD implementation to forming galaxy clusters within a large, cosmological box. Performing a resolution study we demonstrate the robustness of the predicted shape of the magnetic field profiles in galaxy clusters, which is in good agreement with previous studies.

preprint2009arXiv

Imprints of primordial non-Gaussianities in X-ray and SZ signals from galaxy clusters

Several inflationary models predict the possibility that the primordial perturbations of the density field may contain a degree of non-Gaussianity which would influence the subsequent evolution of cosmic structures at large scales. In order to study their impact, we use a set of three cosmological DM-only simulations starting from initial conditions with different levels of non-Gaussianity: f_NL=0,+/-100. More specifically, we focus on the distribution of galaxy clusters at different redshifts and, using suitable scaling relations, we determine their X-ray and SZ signals. Our analysis allows us to estimate the differences in the logN-logS and logN-logY due to the different initial conditions and to predict the cluster counts at different redshifts expected for future surveys (eROSITA and SPT). We also use a second set of simulations assuming a different cosmological scenario to estimate how the dependence on f_NL is degenerate with respect to other parameters. Our results indicate that the effects introduced by a realistic amount of primordial non-Gaussianity are small when compared to the ones connected with current uncertainties in cosmological parameters, particularly with sigma_8. However, if future surveys will be associated with optical follow-up campaigns to determine the cluster redshift, an analysis of the samples at z>1 can provide significant constraints on f_NL. In particular we predict that the SPT cluster survey will be able to detect ~1000 clusters at z>1 for the Gaussian case, with a difference of 15-20 per cent associated to f_NL=+/-100.

preprint2009arXiv

Large-scale non-Gaussian mass function and halo bias: tests on N-body simulations

The description of the abundance and clustering of halos for non-Gaussian initial conditions has recently received renewed interest, motivated by the forthcoming large galaxy and cluster surveys, which can potentially yield constraints of order unity on the non-Gaussianity parameter f_{NL}. We present tests on N-body simulations of analytical formulae describing the halo abundance and clustering for non-Gaussian initial conditions. We calibrate the analytic non-Gaussian mass function of Matarrese et al.(2000) and LoVerde et al.(2008) and the analytic description of clustering of halos for non-Gaussian initial conditions on N-body simulations. We find excellent agreement between the simulations and the analytic predictions if we make the corrections delta_c --> delta_c X sqrt{q} and delta_c --> δ_c X q where q ~ 0.75, in the density threshold for gravitational collapse and in the non-Gaussian fractional correction to the halo bias, respectively. We discuss the implications of this correction on present and forecasted primordial non-Gaussianity constraints. We confirm that the non-Gaussian halo bias offers a robust and highly competitive test of primordial non-Gaussianity.

preprint2009arXiv

Magnetic field structure due to the global velocity field in spiral galaxies

We present a set of global, self-consistent N-body/SPH simulations of the dynamic evolution of galactic discs with gas and including magnetic fields. We have implemented a description to follow the evolution of magnetic fields with the ideal induction equation in the SPH part of the Vine code. Results from a direct implementation of the field equations are compared to a representation by Euler potentials, which pose a div(B)-free description, an constraint not fulfilled for the direct implementation. All simulations are compared to an implementation of magnetic fields in the Gadget code which includes also cleaning methods for div(B). Starting with a homogeneous seed field we find that by differential rotation and spiral structure formation of the disc the field is amplified by one order of magnitude within five rotation periods of the disc. The amplification is stronger for higher numerical resolution. Moreover, we find a tight connection of the magnetic field structure to the density pattern of the galaxy in our simulations, with the magnetic field lines being aligned with the developing spiral pattern of the gas. Our simulations clearly show the importance of non-axisymmetry for the evolution of the magnetic field.

preprint2009arXiv

On the evolution of giant radio halos and their connection with cluster mergers

Giant radio halos are diffuse, Mpc-scale, synchrotron sources located in the central regions of galaxy clusters and provide the most relevant example of cluster non-thermal activity. Radio and X-ray surveys allow to investigate the statistics of halos and may contribute to constrain their origin and evolution. We investigate the distribution of clusters in the plane X-ray (thermal, L_X) vs synchrotron (P_{1.4})luminosity, where clusters hosting giant radio halos trace the P_{1.4}--L_X correlation and clusters without radio halos populate a region that is well separated from that spanned by the above correlation. The connection between radio halos and cluster mergers suggests that the cluster Mpc-scale synchrotron emission is amplified during these mergers and then suppressed when clusters become more dynamically relaxed. In this context, by analysing the distribution in the P_{1.4}--L_X plane of clusters from X-ray selected samples with adequate radio follow up, we constrain the typical time-scale of evolution of diffuse radio emission in clusters and discuss the implications for the origin of radio halos. We conclude that cluster synchrotron emission is suppressed (and amplified) in a time-scale significantly smaller than 1 Gyr. We show that this constraint appears difficult to reconcile with the hypothesis that the halo's radio power is suppressed due to dissipation of magnetic field in galaxy clusters. On the other hand, in agreement with models where turbulent acceleration plays a role, present constraints suggest that relativistic electrons are accelerated in Mpc-scale regions, in connection with cluster mergers and for a time-interval of about 1 Gyr, and then they cool in a relatively small time-scale, when the hosting cluster becomes more dynamically relaxed.

preprint2009arXiv

Radio Halos From Simulations And Hadronic Models I: The Coma cluster

We use the results from a constrained, cosmological MHD simulation of the Local Universe to predict the radio halo and the gamma-ray flux from the Coma cluster and compare it to current observations. The simulated magnetic field within the Coma cluster is the result of turbulent amplification of the magnetic field during build-up of the cluster. The magnetic seed field originates from star-burst driven, galactic outflows. The synchrotron emission is calculated assuming a hadronic model. We follow four approaches with different distributions for the cosmic-ray proton (CRp) population within galaxy clusters. The radial profile the radio halo can only be reproduced with a radially increasing energy fraction within the cosmic ray proton population, reaching $>$100% of the thermal energy content at $\approx$ 1Mpc, e.g. the edge of the radio emitting region. Additionally the spectral steepening of the observed radio halo in Coma cannot be reproduced, even when accounting for the negative flux from the thermal SZ effect at high frequencies. Therefore the hadronic models are disfavored from present analysis. The emission of $γ$-rays expected from our simulated coma is still below the current observational limits (by a factor of $\sim$6) but would be detectable in the near future.

preprint2009arXiv

Simulating the effect of AGN feedback on the metal enrichment of galaxy clusters

We present a study of the effect of AGN feedback on metal enrichment and thermal properties of the intracluster medium (ICM) in hydrodynamical simulations. The cosmological simulations are performed for a set of clusters using a version of the TreePM-SPH Gadget code that follows chemo-dynamical evolution by accounting for metal enrichment by different stellar populations. Besides runs not including any efficient form of energy feedback, we carry out simulations including: (i) kinetic feedback in the form of galactic winds triggered by supernova explosions; (ii) AGN feedback from gas accretion onto super-massive black holes (BHs); (iii) AGN feedback in which a 'radio mode' is included. We find that AGN feedback is able to quench star formation in the brightest cluster galaxies at z<4 and provides correct temperature profiles in the central regions of galaxy groups. However, its effect is not sufficient to create cool cores in massive clusters. AGN feedback creates a widespread enrichment in the outskirts of clusters, thanks to its efficiency in displacing enriched gas from galactic halos to the inter-galactic medium at relatively high redshift. Iron abundance profiles are in better agreement with observations, with a more pristine enrichment of the ICM around and beyond the cluster virial regions. From the pattern of the relative abundances of Silicon and Iron, we conclude that a significant fraction of ICM enrichment in simulations is contributed by a diffuse population of intra-cluster stars. Our simulations also predict that profiles of Z_Si/Z_Fe abundance ratio do not increase at least out to 0.5 R_vir. Our results clearly show that different sources of energy feedback leave distinct imprints in the enrichment pattern of the ICM, that are more evident when looking at cluster external regions.

preprint2009arXiv

Substructure of the galaxy clusters in the REXCESS sample: observed statistics and comparison to numerical simulations

We study the substructure statistics of a representative sample of galaxy clusters by means of two currently popular substructure characterisation methods, power ratios and centroid shifts. We use the 31 clusters from the REXCESS sample, compiled from the southern ROSAT All-Sky cluster survey REFLEX with a morphologically unbiased selection in X-ray luminosity and redshift, all of which have been reobserved with XMM-Newton. We investigate the uncertainties of the substructure parameters and examine the dependence of the results on projection effects, finding that the uncertainties of the parameters can be quite substantial. Thus while the quantification of the dynamical state of individual clusters with these parameters should be treated with extreme caution, these substructure measures provide powerful statistical tools to characterise trends of properties in large cluster samples. The centre shift parameter, w, is found to be more sensitive in general. For the REXCESS sample neither the occurence of substructure nor the presence of cool cores depends on cluster mass. There is a significant anti-correlation between the existence of substantial substructure and cool cores. The simulated clusters show on average larger substructure parameters than the observed clusters, a trend that is traced to the fact that cool regions are more pronounced in the simulated clusters, leading to stronger substructure measures in merging clusters and clusters with offset cores. Moreover, the frequency of cool regions is higher in the simulations than in the observations, implying that the description of the physical processes shaping cluster formation in the simulations requires further improvement.

preprint2009arXiv

Substructures in hydrodynamical cluster simulations

The abundance and structure of dark matter subhalos has been analyzed extensively in recent studies of dark matter-only simulations, but comparatively little is known about the impact of baryonic physics on halo substructures. We here extend the SUBFIND algorithm for substructure identification such that it can be reliably applied to dissipative hydrodynamical simulations that include star formation. This allows, in particular, the identification of galaxies as substructures in simulations of clusters of galaxies, and a determination of their content of gravitationally bound stars, dark matter, and hot and cold gas. Using a large set of cosmological cluster simulations, we present a detailed analysis of halo substructures in hydrodynamical simulations of galaxy clusters, focusing in particular on the influence both of radiative and non-radiative gas physics, and of non-standard physics such as thermal conduction and feedback by galactic outflows. We also examine the impact of numerical nuisance parameters such as artificial viscosity parameterizations. We find that diffuse hot gas is efficiently stripped from subhalos when they enter the highly pressurized cluster atmosphere. This has the effect of decreasing the subhalo mass function relative to a corresponding dark matter-only simulation. These effects are mitigated in radiative runs, where baryons condense in the central subhalo regions and form compact stellar cores. However, in all cases, only a very small fraction, of the order of one percent, of subhalos within the cluster virial radii preserve a gravitationally bound hot gaseous atmosphere. (abridged)

preprint2009arXiv

Weighing simulated galaxy clusters using lensing and X-ray

We aim at investigating potential biases in lensing and X-ray methods to measure the cluster mass profiles. We do so by performing realistic simulations of lensing and X-ray observations that are subsequently analyzed using observational techniques. The resulting mass estimates are compared among them and with the input models. Three clusters obtained from state-of-the-art hydrodynamical simulations, each of which has been projected along three independent lines-of-sight, are used for this analysis. We find that strong lensing models can be trusted over a limited region around the cluster core. Extrapolating the strong lensing mass models to outside the Einstein ring can lead to significant biases in the mass estimates, if the BCG is not modeled properly for example. Weak lensing mass measurements can be largely affected by substructures, depending on the method implemented to convert the shear into a mass estimate. Using non-parametric methods which combine weak and strong lensing data, the projected masses within R200 can be constrained with a precision of ~10%. De-projection of lensing masses increases the scatter around the true masses by more than a factor of two due to cluster triaxiality. X-ray mass measurements have much smaller scatter (about a factor of two smaller than the lensing masses) but they are generally biased low by 5-20%. This bias is ascribable to bulk motions in the gas of our simulated clusters. Using the lensing and the X-ray masses as proxies for the true and the hydrostatic equilibrium masses of the simulated clusters and averaging over the cluster sample we are able to measure the lack of hydrostatic equilibrium in the systems we have investigated.

preprint2008arXiv

Primordial non-Gaussianities in the Intergalactic Medium

We present results from the first high-resolution hydrodynamical simulations of non-Gaussian cosmological models. We focus on the statistical properties of the transmitted Lyman-alpha flux in the high redshift intergalactic medium. Imprints of non-Gaussianity are present and are larger at high redshifts. Differences larger than 20 % at z>3 in the flux probability distribution function for high transmissivity regions (voids) are expected for values of the non linearity parameter f_NL=\pm 100 when compared to a standard LCDM cosmology with f_NL=0. We investigate also the one-dimensional flux bispectrum: at the largest scales (corresponding to tens of Mpc) we expect deviations in the flux bispectrum up to 20% at z~4 (for f_NL=\pm 100), significantly larger than deviations of ~ 3% in the flux power spectrum. We briefly discuss possible systematic errors that can contaminate the signal. Although challenging, a detection of non-Gaussianities in the interesting regime of scales and redshifts probed by the Lyman-alpha forest, could be possible with future data sets.

preprint2008arXiv

The mass density field in simulated non-Gaussian scenarios

In this work we study the properties of the mass density field in the non-Gaussian world models simulated by Grossi et al. 2007. In particular we focus on the one-point density probability distribution function of the mass density field in non-Gausian models with quadratic non-linearities quantified by the usual parameter f_NL. We find that the imprints of primordial non-Gaussianity are well preserved in the negative tail of the probability function during the evolution of the density perturbation. The effect is already noticeable at redshifts as large as 4 and can be detected out to the present epoch. At z=0 we find that the fraction of the volume occupied by regions with underdensity delta < -0.9, typical of voids, is about 1.3 per cent in the Gaussian case and increases to ~2.2 per cent if f_NL=-1000 while decreases to ~0.5 per cent if f_NL=+1000. This result suggests that void-based statistics may provide a powerful method to detect non-Gaussianity even at low redshifts which is complementary to the measurements of the higher-order moments of the probability distribution function like the skewness or the kurtosis for which deviations from the Gaussian case are detected at the 25-50 per cent level.

preprint2008arXiv

UHECR observations and lensing in the magnetic field of the Virgo cluster

We discuss how lensing by magnetic fields in galaxy clusters affects ultrahigh energy cosmic ray (UHECR) observations. As specific example, we use Virgo together with the cluster magnetic fields obtained earlier in a constrained simulation of structure formation including MHD processes. We find that, if M87 is the single source of UHECRs from Virgo, the emitted flux is strongly anisotropic in the most interesting energy range, (50-100)EeV, and differs from the average value by a factor five or more for a significant fraction of observers. Since magnetic lensing is energy dependent, the external energy spectrum as seen by different observers varies strongly too. These anisotropies are averaged out in the case that all active galactic nuclei in Virgo emit UHECRs. In both cases, the anisotropies of the emitted UHECR flux may introduce an important bias in the interpretation of UHECR data like, e.g., the determination of the source density n_s and the source energy spectrum of UHECRs.

preprint2007arXiv

The actual Rees--Sciama effect from the Local Universe

Observations of the Cosmic Microwave Background (CMB) have revealed an unexpected quadrupole-octopole alignment along a preferred axis pointing toward the Virgo cluster. We here investigate whether this feature can be explained in the framework of the concordance model by secondary anisotropies produced by the non-linear evolution of the gravitational potential, the so-called Rees-Sciama (RS) effect. We focus on the effect caused by the local superclusters, which we calculate using a constrained high-resolution hydrodynamical simulation, based on the IRAS 1.2-Jy all-sky galaxy redshift survey, which reproduces the main structures of our Universe out to a distance of 110 Mpc from our Galaxy. The resulting RS effect peaks at low multipoles and has a minimum/maximum amplitude of -6.6μK 1.9μK. Even though its quadrupole is well aligned with the one measured for the CMB, its amplitude is not sufficient to explain the observed magnitude of the quadrupole/octopole alignment. In addition, we analyze the WMAP-3 data with a linear matched filter in an attempt to determine an upper limit for the RS signal amplitude on large scales. We found that it is possible to infer a weak upper limit of 30μK for its maximum amplitude.

preprint2005arXiv

A fast method for computing strong-lensing cross sections: Application to merging clusters

Strong gravitational lensing by irregular mass distributions, such as galaxy clusters, is generally not well quantified by cross sections of analytic mass models. Computationally expensive ray-tracing methods have so far been necessary for accurate cross-section calculations. We describe a fast, semi-analytic method here which is based on surface integrals over high-magnification regions in the lens plane and demonstrate that it yields reliable cross sections even for complex, asymmetric mass distributions. The method is faster than ray-tracing simulations by factors of $\sim30$ and thus suitable for large cosmological simulations, saving large amounts of computing time. We apply this method to a sample of galaxy cluster-sized dark matter haloes with simulated merger trees and show that cluster mergers approximately double the strong-lensing optical depth for lens redshifts $z_\mathrm{l}\gtrsim0.5$ and sources near $z_\mathrm{s} = 2$. We believe that this result hints at one possibility for understanding the recently detected high arcs abundance in clusters at moderate and high redshifts, and is thus worth further studies.

preprint2004arXiv

Locating the Warm-Hot Intergalactic Medium in the Simulated Local Universe

We present an analysis of mock spectral observation of warm-hot intergalactic medium (WHIM) using a constrained simulation of the local universe. The simulated map of oxygen emission lines from local WHIM reproduces well the observed structures traced by galaxies in the real local universe. We further attempt to perform mock observations of outer parts of simulated Coma cluster and A3627 adopting the expected performance of DIOS (Diffuse Intergalactic Oxygen Surveyor), which is proposed as a dedicated soft X-ray mission to search for cosmic missing baryons. We find that WHIMs surrounding nearby clusters are detectable with a typical exposure time of a day, and thus constitute realistic and promising targets for DIOS. We also find that an X-ray emitting clump in front of Coma cluster, recently reported in the XMM-Newton observation, has a counterpart in the simulated local universe, and its observed spectrum can be well reproduced in the simulated local universe if the gas temperature is set to the observationally estimated value.