Showing papers by "Klaus Dolag published in 2016"

Planck 2015 results - XIII. Cosmological parameters

Abstract: This paper presents cosmological results based on full-mission Planck observations of temperature and polarization anisotropies of the cosmic microwave background (CMB) radiation. Our results are in very good agreement with the 2013 analysis of the Planck nominal-mission temperature data, but with increased precision. The temperature and polarization power spectra are consistent with the standard spatially-flat 6-parameter ΛCDM cosmology with a power-law spectrum of adiabatic scalar perturbations (denoted “base ΛCDM” in this paper). From the Planck temperature data combined with Planck lensing, for this cosmology we find a Hubble constant, H0 = (67.8 ± 0.9) km s-1Mpc-1, a matter density parameter Ωm = 0.308 ± 0.012, and a tilted scalar spectral index with ns = 0.968 ± 0.006, consistent with the 2013 analysis. Note that in this abstract we quote 68% confidence limits on measured parameters and 95% upper limits on other parameters. We present the first results of polarization measurements with the Low Frequency Instrument at large angular scales. Combined with the Planck temperature and lensing data, these measurements give a reionization optical depth of τ = 0.066 ± 0.016, corresponding to a reionization redshift of . These results are consistent with those from WMAP polarization measurements cleaned for dust emission using 353-GHz polarization maps from the High Frequency Instrument. We find no evidence for any departure from base ΛCDM in the neutrino sector of the theory; for example, combining Planck observations with other astrophysical data we find Neff = 3.15 ± 0.23 for the effective number of relativistic degrees of freedom, consistent with the value Neff = 3.046 of the Standard Model of particle physics. The sum of neutrino masses is constrained to ∑ mν < 0.23 eV. The spatial curvature of our Universe is found to be very close to zero, with | ΩK | < 0.005. Adding a tensor component as a single-parameter extension to base ΛCDM we find an upper limit on the tensor-to-scalar ratio of r0.002< 0.11, consistent with the Planck 2013 results and consistent with the B-mode polarization constraints from a joint analysis of BICEP2, Keck Array, and Planck (BKP) data. Adding the BKP B-mode data to our analysis leads to a tighter constraint of r0.002 < 0.09 and disfavours inflationarymodels with a V(φ) ∝ φ2 potential. The addition of Planck polarization data leads to strong constraints on deviations from a purely adiabatic spectrum of fluctuations. We find no evidence for any contribution from isocurvature perturbations or from cosmic defects. Combining Planck data with other astrophysical data, including Type Ia supernovae, the equation of state of dark energy is constrained to w = −1.006 ± 0.045, consistent with the expected value for a cosmological constant. The standard big bang nucleosynthesis predictions for the helium and deuterium abundances for the best-fit Planck base ΛCDM cosmology are in excellent agreement with observations. We also constraints on annihilating dark matter and on possible deviations from the standard recombination history. In neither case do we find no evidence for new physics. The Planck results for base ΛCDM are in good agreement with baryon acoustic oscillation data and with the JLA sample of Type Ia supernovae. However, as in the 2013 analysis, the amplitude of the fluctuation spectrum is found to be higher than inferred from some analyses of rich cluster counts and weak gravitational lensing. We show that these tensions cannot easily be resolved with simple modifications of the base ΛCDM cosmology. Apart from these tensions, the base ΛCDM cosmology provides an excellent description of the Planck CMB observations and many other astrophysical data sets.

Planck 2015 results

Abstract: ESA; CNES (France); CNRS/INSU-IN2P3-INP (France); ASI (Italy); CNR (Italy); INAF (Italy); NASA (USA); DoE (USA); STFC (UK); UKSA (UK); CSIC (Spain); MINECO (Spain); JA (Spain); RES (Spain); Tekes (Finland); AoF (Finland); CSC (Finland); DLR (Germany); MPG (Germany); CSA (Canada); DTU Space (Denmark); SER/SSO (Switzerland); RCN (Norway); SFI (Ireland); FCT/MCTES (Portugal); ERC (EU); PRACE (EU); Higher Education Funding Council for England; Science and Technology Facilities Council; Alfred P. Sloan Foundation; National Science Foundation; US Department of Energy Office of Science

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

Abstract: This paper presents the Planck 2015 likelihoods, statistical descriptions of the 2-point correlationfunctions of the cosmic microwave background (CMB) temperature and polarization fluctuations that account for relevant uncertainties, both instrumental and astrophysical in nature. They are based on the same hybrid approach used for the previous release, i.e., a pixel-based likelihood at low multipoles (l< 30) and a Gaussian approximation to the distribution of cross-power spectra at higher multipoles. The main improvements are the use of more and better processed data and of Planck polarization information, along with more detailed models of foregrounds and instrumental uncertainties. The increased redundancy brought by more than doubling the amount of data analysed enables further consistency checks and enhanced immunity to systematic effects. It also improves the constraining power of Planck, in particular with regard to small-scale foreground properties. Progress in the modelling of foreground emission enables the retention of a larger fraction of the sky to determine the properties of the CMB, which also contributes to the enhanced precision of the spectra. Improvements in data processing and instrumental modelling further reduce uncertainties. Extensive tests establish the robustness and accuracy of the likelihood results, from temperature alone, from polarization alone, and from their combination. For temperature, we also perform a full likelihood analysis of realistic end-to-end simulations of the instrumental response to the sky, which were fed into the actual data processing pipeline; this does not reveal biases from residual low-level instrumental systematics. Even with the increase in precision and robustness, the ΛCDM cosmological model continues to offer a very good fit to the Planck data. The slope of the primordial scalar fluctuations, n_s, is confirmed smaller than unity at more than 5σ from Planck alone. We further validate the robustness of the likelihood results against specific extensions to the baseline cosmology, which are particularly sensitive to data at high multipoles. For instance, the effective number of neutrino species remains compatible with the canonical value of 3.046. For this first detailed analysis of Planck polarization spectra, we concentrate at high multipoles on the E modes, leaving the analysis of the weaker B modes to future work. At low multipoles we use temperature maps at all Planck frequencies along with a subset of polarization data. These data take advantage of Planck’s wide frequency coverage to improve the separation of CMB and foreground emission. Within the baseline ΛCDM cosmology this requires τ = 0.078 ± 0.019 for the reionization optical depth, which is significantly lower than estimates without the use of high-frequency data for explicit monitoring of dust emission. At high multipoles we detect residual systematic errors in E polarization, typically at the μK^2 level; we therefore choose to retain temperature information alone for high multipoles as the recommended baseline, in particular for testing non-minimal models. Nevertheless, the high-multipole polarization spectra from Planck are already good enough to enable a separate high-precision determination of the parameters of the ΛCDM model, showing consistency with those established independently from temperature information alone.

XXIV. Cosmology from Sunyaev-Zeldovich cluster counts

Abstract: We present cluster counts and corresponding cosmological constraints from the Planck full mission data set. Our catalogue consists of 439 clusters detected via their Sunyaev-Zeldovich (SZ) signal d ...

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

Abstract: The polarized thermal emission from diffuse Galactic dust is the main foreground present in measurements of the polarization of the cosmic microwave background (CMB) at frequencies above 100 GHz. In this paper we exploit the uniqueness of the Planck HFI polarization data from 100 to 353 GHz to measure the polarized dust angular power spectra C_l^(EE) and C_l^(BB) over the multipole range 40

Planck 2015 results: XIX. Constraints on primordial magnetic fields

Abstract: We predict 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 temperature and polarization spectra, related to their contribution to cosmological perturbations; the effect on CMB polarization induced by Faraday rotation; magnetically-induced non-Gaussianities and related non-zero bispectra; and the magnetically-induced breaking of statistical isotropy. We present constraints on the amplitude of PMFs derived from different combinations of Planck data products, depending on the specific effect that is analysed. Overall, Planck data constrain the amplitude of PMFs to less than a few nanogauss, with different bounds depending on the considered model. In particular, individual limits coming from the analysis of the CMB angular power spectra, using the Planck likelihood, are B1Mpc < 4:4 nG (where B1Mpc is the comoving field amplitude at a scale of 1 Mpc) at 95 % confidence level, assuming zero helicity, and B1Mpc < 5:6 nG when we consider a maximally helical field. For nearly scaleinvariant PMFs we obtain B1Mpc < 2:1 nG and B1Mpc < 0:7 nG if the impact of PMFs on the ionization history of the Universe is included in the analysis. 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 B1Mpc < 2:8 nG. A search for preferred directions in the magnetically-induced passive bispectrum yields B1Mpc < 4:5 nG, whereas the the compensated-scalar bispectrum gives B1Mpc < 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 B1Mpc < 1380 nG. In our final analysis, we consider the harmonic-space correlations produced by Alfv´ en 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.

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

Abstract: We have constructed all-sky Compton parameters maps, 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 These reconstructed y-maps are delivered as part of the Planck 2015 release The y-maps are characterized in terms of noise properties and residual foreground contamination, mainly thermal dust emission at large angular scales, and cosmic infrared background 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

Halo mass function: baryon impact, fitting formulae, and implications for cluster cosmology

Abstract: We use a set of hydrodynamical and dark matter-only (DMonly) simulations to calibrate the halo mass function (HMF). We explore the impact of baryons, propose an improved parametrization for spherical overdensity masses, and identify differences between our DMonly HMF and previously published HMFs. We use the Magneticum simulations, which are well suited because of their accurate treatment of baryons, high resolution, and large cosmological volumes of up to (3818 Mpc)(3). Baryonic effects globally decrease the masses of galaxy clusters, which, at a given mass, results in a decrease of their number density. This effect vanishes at high redshift z similar to 2 and for high masses M-200m greater than or similar to 10(14)M(circle dot). We perform cosmological analyses of three idealized approximations to the cluster surveys by the South Pole Telescope (SPT), Planck, and eROSITA. We pursue two main questions. (1) What is the impact of baryons? - for the SPT-like and the Planck-like samples, the impact of baryons on cosmological results is negligible. In the eROSITA-like case, however, neglecting the baryonic impact leads to an underestimate of Omega(m) by about 0.01, which is comparable to the expected uncertainty from eROSITA. (2) How does our DMonly HMF compare with previous work? - for the Planck-like sample, results obtained using our DMonly HMF are shifted by Delta(sigma(8)) similar or equal to (sigma(8)(Omega(m)/0.27)(0.3)) similar or equal to 0.02 with respect to results obtained using the Tinker et al. fit. This suggests that using our HMF would shift results from Planck clusters towards better agreement with cosmic-microwave-background anisotropy measurements. Finally, we discuss biases that can be introduced through inadequate HMF parametrizations that introduce false cosmological sensitivity.

An improved SPH scheme for cosmological simulations

Abstract: We present an implementation of smoothed particle hydrodynamics (SPH) with improved accuracy for simulations of galaxies and the large-scale structure. In particular, we implement and test a vast majority of SPH improvement in the developer version of GADGET-3. We use the Wendland kernel functions, a particle wake-up time-step limiting mechanism and a time-dependent scheme for artificial viscosity including 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 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 GADGET-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 discs to be colder and more extended and galaxy clusters showing 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 GADGET-SPH, thus becoming the core of an efficient code for large cosmological simulations.

Planck 2015 results - XII. Full focal plane simulations

Abstract: 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 18 144 maps, together with the most massive suite of Monte Carlo realizations of instrument noise and CMB ever generated, comprising 104 mission realizations reduced to about 106 maps. The resulting maps incorporate the dominant instrumental, scanning, and data analysis effects, and the 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 to validate and verify analysis algorithms and their implementations, and to remove biases from and quantify uncertainties in the results of analyses of the real data.

SZ effects in the Magneticum Pathfinder Simulation: Comparison with the Planck, SPT, and ACT results

Abstract: We calculate the one-point probability density distribution functions (PDF) and the power spectra of the thermal and kinetic Sunyaev-Zeldovich (tSZ and kSZ) effects and the mean Compton Y parameter using the Magneticum Pathfinder simulations, state-of-the-art cosmological hydrodynamical simulations of a large cosmological volume of (896 Mpc h(-1))(3). These simulations follow in detail the thermal and chemical evolution of the intracluster medium as well as the evolution of supermassive black holes and their associated feedback processes. We construct full-sky maps of tSZ and kSZ from the light-cones out to z = 0.17, and one realization of 8.degrees 8 x 8.degrees 8 deep light-cone out to z = 5.2. The local universe at z = 1.18 x 10(-6) for Omega(m) = 0.272 and sigma(8) = 0.809. Nearly half (approximate to 5 x 10(-7)) of the signal comes from haloes below a virial mass of 10(13) M-circle dot h(-1). Scaling this to the Planck 2015 parameters, we find (Y) over bar = 1.57 x 10(-6).

On the Nature of Hydrostatic Equilibrium in Galaxy Clusters

Abstract: In this paper, we investigate the level of hydrostatic equilibrium (HE) in the intracluster 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 are. stellar and active galactic nucleus 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 and. 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 the deviation from HE. very well, 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 (greater than or similar to 30%), relative to thermal ones. As a general result, the HE-deviation and mass bias, at a. given distance from the cluster center, are not very sensitive to the temperature inhomogeneities in the gas.

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

Abstract: By looking at the kinetic Sunyaev-Zeldovich effect (kSZ) in Planck nominal mission data, we present a significant detection of b aryons participating in large-scale bulk flows around central galaxies (CG s) at redshift z≈ 0.1. We estimate the pairwise momentum of the kSZ temperature fluctuations at the positions of the Central Galaxy Catalogu e (CGC) samples extracted from Sloan Digital Sky Survey (SDSS-DR7) data. For the foreground-cleanedSEVEM,SMICA,NILC, andCOMMANDER 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 comp ute for the first time the cross-correlation function betwee n 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 the kSZ temperature-velocity field correlation find evidence for kSZ signatures out to aper tures of 8 arcmin and beyond, corresponding to a physical radius of> 1 Mpc, more than twice the mean virial radius of halos. This is consisten t with the predictions from hydrodynamical simulations that most of the baryons are outside the virialized halos. We fit a simple model, in which t he temperature-velocity cross-correlation is proportion al to the signal seen in a semi-analytic model built upon N-body simulations, and interpret the proportionality cons tant as an effective optical depth to Thomson scattering. We find τT = (1.4± 0.5)× 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 Sunyaev-Zeldovich observations.

On the nature of hydrostatic equilibrium in galaxy clusters

Abstract: 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.

Origin and properties of dual and offset active galactic nuclei in a cosmological simulation at z=2

Abstract: In the last few years, it became possible to observationally resolve galaxies with two distinct nuclei in their centre. For separations smaller than 10 kpc, dual and offset active galactic nuclei (AGN) are distinguished: in dual AGN, both nuclei are active, whereas in offset AGN only one nucleus is active. To study the origin of such AGN pairs, we employ a cosmological, hydrodynamic simulation with a large volume of (182 Mpc)(3) from the set of Magneticum Pathfinder Simulations. The simulation self-consistently produces 35 resolved black hole (BH) pairs at redshift z = 2, with a comoving distance smaller than 10 kpc. 14 of them are offset AGN and nine are dual AGN, resulting in a fraction of (1.2 +/- 0.3) per cent AGN pairs with respect to the total number of AGN. In this paper, we discuss fundamental differences between the BH and galaxy properties of dual AGN, offset AGN and inactive BH pairs and investigate their different triggering mechanisms. We find that in dual AGN the BHs have similar masses and the corresponding BH from the less massive progenitor galaxy always accretes with a higher Eddington ratio. In contrast, in offset AGN the active BH is typically more massive than its non-active counterpart. Furthermore, dual AGN in general accrete more gas from the intergalactic medium than offset AGN and non-active BH pairs. This highlights that merger events, particularly minor mergers, do not necessarily lead to strong gas inflows and thus, do not always drive strong nuclear activity.

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

Abstract: 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.

DEMNUni: ISW, Rees-Sciama, and weak-lensing in the presence of massive neutrinos

Abstract: We present, for the first time in the literature, a full reconstruction of the total (linear and non-linear) ISW/Rees-Sciama effect in the presence of massive neutrinos, together with its cross-correlations with CMB-lensing and weak-lensing signals. The present analyses make use of all-sky maps extracted via ray-tracing across the gravitational potential distribution provided by the "Dark Energy and Massive Neutrino Universe" (DEMNUni) project, a set of large-volume, high-resolution cosmological N-body simulations, where neutrinos are treated as separate collisionless particles. We correctly recover, at 1-2% accuracy, the linear predictions from CAMB. Concerning the CMB-lensing and weak-lensing signals, we also recover, with similar accuracy, the signal predicted by Boltzmann codes, once non-linear neutrino corrections to HALOFIT are accounted for. Interestingly, in the ISW/Rees-Sciama signal, and its cross correlation with lensing, we find an excess of power with respect to the massless case, due to free streaming neutrinos, roughly at the transition scale between the linear and non-linear regimes. The excess is similar to 5-10% at l similar to 100 for the ISW/Rees-Sciama auto power spectrum, depending on the total neutrino mass M-v, and becomes a factor of similar to 4 for M-v = 0.3 eV, at l similar to 600, for the ISW/Rees-Sciama cross power with CMB-lensing. This effect should be taken into account for the correct estimation of the CMB temperature bispectrum in the presence of massive neutrinos.

New constraints on modelling the random magnetic field of the MW

Abstract: We extend the description of the isotropic and anisotropic random component of the small-scale magnetic field within the existing magnetic field model of the Milky Way from Jansson Farrar, by including random realizations of the small-scale component. Using a magnetic-field power spectrum with Gaussian random fields, the NE2001 model for the thermal electrons and the Galactic cosmic-ray electron distribution from the current GALPROP model we derive full-sky maps for the total and polarized synchrotron intensity as well as the Faraday rotation-measure distribution. While previous work assumed that small-scale fluctuations average out along the line-of-sight or which only computed ensemble averages of random fields, we show that these fluctuations need to be carefully taken into account. Comparing with observational data we obtain not only good agreement with 408 MHz total and WMAP7 22 GHz polarized intensity emission maps, but also an improved agreement with Galactic foreground rotation-measure maps and power spectra, whose amplitude and shape strongly depend on the parameters of the random field. We demonstrate that a correlation length of approximate to 220 pc (50 pc being a 5 sigma lower limit) is needed to match the slope of the observed power spectrum of Galactic foreground rotation-measure maps. Using multiple realizations allows us also to infer errors on individual observables. We find that previously used amplitudes for random and anisotropic random magnetic field components need to be rescaled by factors of approximate to 0.3 and 0.6 to account for the new small-scale contributions. Our model predicts a rotation measure of -2.8 +/- 7.1 rad/m(2) and 4.4 +/- 11.0 rad/m(2) for the north and south Galactic poles respectively, in good agreement with observations. Applying our model to deflections of ultra-high-energy cosmic rays we infer a mean deflection of approximate to 3.5 +/- 1.1 degree for 60 EeV protons arriving from CenA.

Neutral hydrogen in galaxy clusters: impact of AGN feedback and implications for intensity mapping

Abstract: By means of zoom-in hydrodynamic simulations, we quantify the amount of neutral hydrogen (H I) hosted by groups and clusters of galaxies. Our simulations, which are based on an improved formulation of smoothed particle hydrodynamics, include radiative cooling, star formation, metal enrichment and supernova feedback, and can be split into two different groups, depending on whether feedback from active galactic nuclei (AGN) is turned on or off. Simulations are analysed to account for HI self-shielding and the presence of molecular hydrogen. We find that the mass in neutral hydrogen of dark matter haloes monotonically increases with the halomass and can be well described by a power law of the form M-H I (M, z) proportional to M-3/4. Our results point out that AGN feedback reduces both the total halo mass and its HI mass, although it is more efficient in removing HI. We conclude that AGN feedback reduces the neutral hydrogen mass of a given halo by similar to 50 per cent, with a weak dependence on halo mass and redshift. The spatial distribution of neutral hydrogen within haloes is also affected by AGN feedback, whose effect is to decrease the fraction of HI that resides in the halo inner regions. By extrapolating our results to haloes not resolved in our simulations, we derive astrophysical implications from the measurements of Omega(H) (I)(z): haloes with circular velocities larger than similar to 25 km s(-1) are needed to host HI in order to reproduce observations. We find that only the model with AGN feedback is capable of reproducing the value of Omega(HI)b(HI) derived from available 21 cm intensity mapping observations.

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

Abstract: 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, enable a detailed study of this big object through the SZ effect. Virgo is well resolved by Planck , showing an elongated structure that correlates well with the morphology observed from X-rays, but extends beyond the observed X-ray signal. We find good agreement between the SZ signal (or Compton parameter, y c ) observed by Planck and the expected signal inferred from X-ray observations and simple analytical models. Owing to its proximity to us, the gas beyond the virial radius in Virgo 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 elusivewarm/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. Finally, based on the combination of the same SZ and X-ray data, we constrain the total amount of gas in Virgo. Under the hypothesis that the abundance of baryons in Virgo is representative of the cosmic average, we also infer a distance for Virgo of approximately 18 Mpc, in good agreement with previous estimates.

A web portal for hydrodynamical, cosmological simulations

Abstract: This article describes a data center hosting a web portal for accessing and sharing the output of large, cosmological, hydro-dynamical simulations with a broad scientific community. It also allows users to receive related scientific data products by directly processing the raw simulation data on a remote computing cluster. The data center has a multi-layer structure: a web portal, a job control layer, a computing cluster and a HPC storage system. The outer layer enables users to choose an object from the simulations. Objects can be selected by visually inspecting 2D maps of the simulation data, by performing highly compounded and elaborated queries or graphically by plotting arbitrary combinations of properties. The user can run analysis tools on a chosen object. These services allow users to run analysis tools on the raw simulation data. The job control layer is responsible for handling and performing the analysis jobs, which are executed on a computing cluster. The innermost layer is formed by a HPC storage system which hosts the large, raw simulation data. The following services are available for the users: (I) {\sc ClusterInspect} visualizes properties of member galaxies of a selected galaxy cluster; (II) {\sc SimCut} returns the raw data of a sub-volume around a selected object from a simulation, containing all the original, hydro-dynamical quantities; (III) {\sc Smac} creates idealised 2D maps of various, physical quantities and observables of a selected object; (IV) {\sc Phox} generates virtual X-ray observations with specifications of various current and upcoming instruments.

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

Abstract: 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.

DEMNUni: ISW, Rees-Sciama, and weak-lensing in the presence of massive neutrinos

Abstract: We present, for the first time in the literature, a full reconstruction of the total (linear and non-linear) ISW/Rees-Sciama effect in the presence of massive neutrinos, together with its cross-correlations with CMB-lensing and weak-lensing signals. The present analyses make use of all-sky maps extracted via ray-tracing across the gravitational potential distribution provided by the "Dark Energy and Massive Neutrino Universe" (DEMNUni) project, a set of large-volume, high-resolution cosmological N-body simulations, where neutrinos are treated as separate collisionless particles. We correctly recover, at $1-2\%$ accuracy, the linear predictions from CAMB. Concerning the CMB-lensing and weak-lensing signals, we also recover, with similar accuracy, the signal predicted by Boltzmann codes, once non-linear neutrino corrections to Halofit are accounted for. Interestingly, in the ISW/Rees-Sciama signal, and its cross correlation with lensing, we find an excess of power with respect to the massless case, due to free streaming neutrinos, roughly at the transition scale between the linear and non-linear regimes. The excess is $\sim 5-10\%$ at $l\sim 100$ for the ISW/Rees-Sciama auto power spectrum, depending on the total neutrino mass $M_ u$, and becomes a factor of $\sim 4$ for $M_ u=0.3$ eV, at $l\sim 600$, for the ISW/Rees-Sciama cross power with CMB-lensing. This effect should be taken into account for the correct estimation of the CMB temperature bispectrum in the presence of massive neutrinos.

Geometrical on-the-fly shock detection in smoothed particle hydrodynamics

Abstract: 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 setups. 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.

Extreme Scale-out SuperMUC Phase 2 - lessons learned

Abstract: In spring 2015, the Leibniz Supercomputing Centre (Leibniz-Rechenzentrum, LRZ), installed their new Peta-Scale System SuperMUC Phase2. Selected users were invited for a 28 day extreme scale-out block operation during which they were allowed to use the full system for their applications. The following projects participated in the extreme scale-out workshop: BQCD (Quantum Physics), SeisSol (Geophysics, Seismics), GPI-2/GASPI (Toolkit for HPC), Seven-League Hydro (Astrophysics), ILBDC (Lattice Boltzmann CFD), Iphigenie (Molecular Dynamic), FLASH (Astrophysics), GADGET (Cosmological Dynamics), PSC (Plasma Physics), waLBerla (Lattice Boltzmann CFD), Musubi (Lattice Boltzmann CFD), Vertex3D (Stellar Astrophysics), CIAO (Combustion CFD), and LS1-Mardyn (Material Science). The projects were allowed to use the machine exclusively during the 28 day period, which corresponds to a total of 63.4 million core-hours, of which 43.8 million core-hours were used by the applications, resulting in a utilization of 69%. The top 3 users were using 15.2, 6.4, and 4.7 million core-hours, respectively.

An online theoretical virtual observatory for hydrodynamical, cosmological simulations

Abstract: This article describes a virtual observatory hosting a web portal for accessing and sharing the output of large, cosmological, hydro-dynamical simulations with a broad scientific community. It also allows users to receive related scientific data products by directly processing the raw simulation data on a remote computing cluster. The virtual observatory is a multi-layer structure: a web portal, a job control layer, a computing cluster and a HPC storage system. The outer layer enables users to choose an object from the simulations. Objects can be selected by visually inspecting 2D maps of the simulation data, by performing highly compounded and elaborated queries or graphically from plotting arbitrary combinations of properties. The user can apply several services to a chosen object. These services allow users to run analysis tools on the raw simulation data. The job control layer is responsible for handling and performing the analysis jobs, which are executed on a computing cluster. The inner most layer is formed by a HPC storage system which host the large, raw simulation data. The virtual observatory provides the following services for the users: (I) ClusterInspect visualizes properties of member galaxies of a selected galaxy cluster; (II) SimCut returns the raw data of a sub-volume around a selected object from a simulation, containing all the original, hydro-dynamical quantities; (III) Smac creates idealised 2D maps of various, physical quantities and observable of a selected object; (IV) Phox generates virtual X-ray observations with specifications of various current and upcoming instruments.

Extreme Scale-out SuperMUC Phase 2 - lessons learned

Abstract: In spring 2015, the Leibniz Supercomputing Centre (Leibniz-Rechenzentrum, LRZ), installed their new Peta-Scale System SuperMUC Phase2. Selected users were invited for a 28 day extreme scale-out block operation during which they were allowed to use the full system for their applications. The following projects participated in the extreme scale-out workshop: BQCD (Quantum Physics), SeisSol (Geophysics, Seismics), GPI-2/GASPI (Toolkit for HPC), Seven-League Hydro (Astrophysics), ILBDC (Lattice Boltzmann CFD), Iphigenie (Molecular Dynamic), FLASH (Astrophysics), GADGET (Cosmological Dynamics), PSC (Plasma Physics), waLBerla (Lattice Boltzmann CFD), Musubi (Lattice Boltzmann CFD), Vertex3D (Stellar Astrophysics), CIAO (Combustion CFD), and LS1-Mardyn (Material Science). The projects were allowed to use the machine exclusively during the 28 day period, which corresponds to a total of 63.4 million core-hours, of which 43.8 million core-hours were used by the applications, resulting in a utilization of 69%. The top 3 users were using 15.2, 6.4, and 4.7 million core-hours, respectively.

Hydrodynamic Simulations of Galaxy Clusters: Scaling Relations and Evolution

Abstract: We analyze hydrodynamical and cosmological simulations of galaxy clusters to study scaling relations between the cluster total masses and observable quantities such as gas luminosity, gas mass, temperature, and YX , i.e., the product of the last two properties. Our simulations are performed with the Smoothed-Particle-Hydrodynamic GADGET-3 code and include different physical processes. The twofold aim of our study is to compare our simulated scaling relations with observations at low (z~0) and intermediate (z~0.5) redshifts and to explore their evolution over the redshift range z=0-2. The result of the comparative study shows a good agreement between our numerical models and real data. We find that AGN feedback significantly affects low-mass haloes at the highest redshifts resulting in a reduction of the slope of the mass-gas mass relation (~13%) and the mass-YX relation (~10%) at z=2 in comparison to z=0. The drop of the slope of the mass-temperature relation at z=2 (~14%) is, instead, caused by early mergers. We investigate the impact of the slope variation on the study of the evolution of the normalization. We conclude that the observed scaling relations should be limited to the redshift range z=0-1 for cosmological studies because in that redshift range the slope, the scatter, and the covariance matrix of the relations do not exhibit significant evolution. The mass-YX relation continues to be the most suitable relation for this goal. Extending the analysis to the redshift range between 1 and 2 will be crucial to evaluate the impact generated by the AGN feedback.

The effect of AGN feedback on Sunyaev-Zeldovich properties of simulated galaxy clusters

Abstract: Abstract We studied the imprints that feedback from Active Galactic Nuclei (AGN) leaves on the intracluster plasma during the assembly history of galaxy clusters. To this purpose we used state-of-the-art cosmological hydrodynamical simulations based on an updated version of the Tree-PM SPH GADGET-3 code, comparing three sets of simulations with different prescriptions for the physics of baryons (including AGN and/or stellar feedback). We explore the effect of these different physics, in particular AGN feedback, on IntraCluster medium (ICM) properties observed via Sunyaev-Zel’dovich (SZ) effect using an extended set of galaxy clusters (~100 clusters with M 500 masses above 5 × 1013 M ⊙/h). Some of the main findings are that the scaling relation between the integrated SZ flux and the galaxy cluster total mass is in good accordance with several observed samples, especially for massive clusters, and does not show any clear redshift evolution, with the slope of the relation close to the theoretical one in the AGN feedback case. As for the scatter of this relation, we obtain a mild dependence on the cluster dynamical state.