Observations Of Galaxy Clusters Research Articles

Cosmological perturbations on the phantom brane

We obtain a closed system of equations for scalar perturbations in a multi-component braneworld. Our braneworld possesses a phantom-like equation of state at late times, weff < −1, but no big-rip future singularity. In addition to matter and radiation, the braneworld possesses a new effective degree of freedom—the `Weyl fluid' or `dark radiation'. Setting initial conditions on super-Hubble spatial scales at the epoch of radiation domination, we evolve perturbations of radiation, pressureless matter and the Weyl fluid until the present epoch. We observe a gradual decrease in the amplitude of the Weyl-fluid perturbations after Hubble-radius crossing, which results in a negligible effect of the Weyl fluid on the evolution of matter perturbations on spatial scales relevant for structure formation. Consequently, the quasi-static approximation of Koyama and Maartens provides a good fit to the exact results during the matter-dominated epoch. We find that the late-time growth of density perturbations on the brane proceeds at a faster rate than in ΛCDM. Additionally, the gravitational potentials Φ and Ψ evolve differently on the brane than in ΛCDM, for which Φ = Ψ. On the brane, by contrast, the ratio Φ/Ψ exceeds unity during the late matter-dominated epoch (z ≲ 50). These features emerge as smoking gun tests of phantom brane cosmology and allow predictions of this scenario to be tested against observations of galaxy clustering and large-scale structure.

Observing dynamical friction in galaxy clusters

We present a novel method to detect the effects of dynamical friction in observed galaxy clusters. Following accretion into clusters, massive satellite galaxies will backsplash to systematically smaller radii than less massive satellites, an effect that may be detected by stacking the number density profiles of galaxies around clusters. We show that this effect may be understood using a simple toy model which reproduces the trends with halo properties observed in simulations. We search for this effect using SDSS redMaPPer clusters with richness 010 < λ < 2, and find that bright (Mi < −21.5) satellites have smaller splashback radii than fainter (0Mi > −2) satellites at 99% confidence.

Covariance in the thermal SZ–weak lensing mass scaling relation of galaxy clusters

The thermal Sunyaev-Zel'dovich (tSZ) effect signal is widely recognized as a robust mass proxy of galaxy clusters with small intrinsic scatter. However, recent observational calibration of the tSZ scaling relation using weak lensing (WL) mass exhibits considerably larger scatter than the intrinsic scatter predicted from numerical simulations. This raises a question as to whether we can realize the full statistical power of ongoing and upcoming tSZ-WL observations of galaxy clusters. In this work, we investigate the origin of observed scatter in the tSZ-WL scaling relation, using mock maps of galaxy clusters extracted from cosmological hydrodynamic simulations. We show that the inferred intrinsic scatter from mock tSZ-WL analyses is considerably larger than the intrinsic scatter measured in simulations, and comparable to the scatter in the observed tSZ-WL relation. We show that this enhanced scatter originates from the combination of the projection of correlated structures along the line of sight and the uncertainty in the cluster radius associated with WL mass estimates, causing the amplitude of the scatter to depend on the covariance between tSZ and WL signals. We present a statistical model to recover the unbiased cluster scaling relation and cosmological parameter by taking into account the covariance in the tSZ-WL mass relation from multi-wavelength cluster surveys.

Observational capabilities of the new medium- and low-resolution spectrograph at the 1.6-m telescope of the Sayan Observatory

The capabilities of the new medium- and low-resolution spectrograph installed recently on the 1.6-m AZT-33IK telescope at the Sayan Observatory of the Institute of Solar–Terrestrial Physics to solve the problems of ground-based optical support for the future all-skyX-ray survey of the SRGobservatory are discussed. Results of the test observations of galaxy clusters, active galactic nuclei (AGNs) and quasars, and cataclysmic variables performed immediately after the installation of the spectrograph on the telescope are presented. The results of these observations show that the AZT-33IK telescope equipped with the new medium- and low-resolution spectrograph can provide a substantial fraction of the necessary optical observations in the program of ground-based optical support for the all-sky survey of the SRG observatory.

The galaxy correlation function as a constraint on galaxy formation physics

We introduce methods which allow observed galaxy clustering to be used together with observed luminosity or stellar mass functions to constrain the physics of galaxy formation. We show how the projected two-point correlation function of galaxies in a large semi-analytic simulation can be estimated to better than ~10% using only a very small subsample of the subhalo merger trees. This allows measured correlations to be used as constraints in a Monte Carlo Markov Chain exploration of the astrophysical and cosmological parameter space. An important part of our scheme is an analytic profile which captures the simulated satellite distribution extremely well out to several halo virial radii. This is essential to reproduce the correlation properties of the full simulation at intermediate separations. As a first application, we use low-redshift clustering and abundance measurements to constrain a recent version of the Munich semi-analytic model. The preferred values of most parameters are consistent with those found previously, with significantly improved constraints and somewhat shifted "best" values for parameters that primarily affect spatial distributions. Our methods allow multi-epoch data on galaxy clustering and abundance to be used as joint constraints on galaxy formation. This may lead to significant constraints on cosmological parameters even after marginalising over galaxy formation physics.

Gaussian covariance matrices for anisotropic galaxy clustering measurements

Measurements of the redshift-space galaxy clustering have been a prolific source of cosmological information in recent years. Accurate covariance estimates are an essential step for the validation of galaxy clustering models of the redshift-space two-point statistics. Usually, only a limited set of accurate N-body simulations is available. Thus, assessing the data covariance is not possible or only leads to a noisy estimate. Further, relying on simulated realisations of the survey data means that tests of the cosmology dependence of the covariance are expensive. With these points in mind, this work presents a simple theoretical model for the linear covariance of anisotropic galaxy clustering observations with synthetic catalogues. Considering the Legendre moments (`multipoles') of the two-point statistics and projections into wide bins of the line-of-sight parameter (`clustering wedges'), we describe the modelling of the covariance for these anisotropic clustering measurements for galaxy samples with a trivial geometry in the case of a Gaussian approximation of the clustering likelihood. As main result of this paper, we give the explicit formulae for Fourier and configuration space covariance matrices. To validate our model, we create synthetic HOD galaxy catalogues by populating the haloes of an ensemble of large-volume N-body simulations. Using linear and non-linear input power spectra, we find very good agreement between the model predictions and the measurements on the synthetic catalogues in the quasi-linear regime.

Comparison of weak lensing by NFW and Einasto halos and systematic errors

Recent N-body simulations have shown that Einasto radial profiles provide the most accurate description of dark matter halos. Predictions based on the traditional NFW functional form may fail to describe the structural properties of cosmic objects at the percent level required by precision cosmology. We computed the systematic errors expected for weak lensing analyses of clusters of galaxies if one wrongly models the lens density profile. Even though the NFW fits of observed tangential shear profiles can be excellent, viral masses and concentrations of very massive halos (≳ 1015M⊙/h) can be over- and underestimated by 0∼ 1 per cent, respectively. Misfitting effects also steepen the observed mass-concentration relation, as observed in multi-wavelength observations of galaxy groups and clusters. Based on shear analyses, Einasto and NFW halos can be set apart either with deep observations of exceptionally massive structures (≳ 2×1015M⊙/h) or by stacking the shear profiles of thousands of group-sized lenses (≳ 1014M⊙/h).

SPECTRAL IMAGING OF GALAXY CLUSTERS WITHPLANCK

The Sunyaev-Zeldovich (SZ) effect is a promising tool for detecting the presence of hot gas out to the galaxy cluster peripheries. We developed a spectral imaging algorithm dedicated to the SZ observations of nearby galaxy clusters with Planck, with the aim of revealing gas density anisotropies related to the filamentary accretion of materials, or pressure discontinuities induced by the propagation of shock fronts. To optimize an unavoidable trade-off between angular resolution and precision of the SZ flux measurements, the algorithm performs a multiscale analysis of the SZ maps as well as of other extended components, such as the cosmic microwave background (CMB) anisotropies and the Galactic thermal dust. The demixing of the SZ signal is tackled through kernel weighted likelihood maximizations. The CMB anisotropies are further analyzed through a wavelet analysis, while the Galactic foregrounds and SZ maps are analyzed via a curvelet analysis that best preserves their anisotropic details. The algorithm performance has been tested against mock observations of galaxy clusters obtained by simulating the Planck High Frequency Instrument and by pointing a few characteristic positions in the sky. These tests suggest that Planck should easily allow us to detect filaments in the cluster peripheries and detect large-scale shocks in colliding galaxy clusters that feature favorable geometry.

Galaxy clustering, photometric redshifts and diagnosis of systematics in the DES Science Verification data

We study the clustering of galaxies detected at $i<22.5$ in the Science Verification observations of the Dark Energy Survey (DES). Two-point correlation functions are measured using $2.3\times 10^6$ galaxies over a contiguous 116 deg$^2$ region in five bins of photometric redshift width $\Delta z = 0.2$ in the range $0.2 < z < 1.2.$ The impact of photometric redshift errors are assessed by comparing results using a template-based photo-$z$ algorithm (BPZ) to a machine-learning algorithm (TPZ). A companion paper (Leistedt et al 2015) presents maps of several observational variables (e.g. seeing, sky brightness) which could modulate the galaxy density. Here we characterize and mitigate systematic errors on the measured clustering which arise from these observational variables, in addition to others such as Galactic dust and stellar contamination. After correcting for systematic effects we measure galaxy bias over a broad range of linear scales relative to mass clustering predicted from the Planck $\Lambda$CDM model, finding agreement with CFHTLS measurements with $\chi^2$ of 4.0 (8.7) with 5 degrees of freedom for the TPZ (BPZ) redshifts. We test a "linear bias" model, in which the galaxy clustering is a fixed multiple of the predicted non-linear dark-matter clustering. The precision of the data allow us to determine that the linear bias model describes the observed galaxy clustering to $2.5\%$ accuracy down to scales at least $4$ to $10$ times smaller than those on which linear theory is expected to be sufficient.

Is there scale-dependent bias in single-field inflation?

Scale-dependent halo bias due to local primordial non-Gaussianity provides a strong test of single-field inflation. While it is universally understood that single-field inflation predicts negligible scale-dependent bias compared to current observational uncertainties, there is still disagreement on the exact level of scale-dependent bias at a level that could strongly impact inferences made from future surveys. In this paper, we clarify this confusion and derive in various ways that there is exactly zero scale-dependent bias in single-field inflation. Much of the current confusion follows from the fact that single-field inflation does predict a mode coupling of matter perturbations at the level of fNLlocal; ≈ −5/3, which naively would lead to scale-dependent bias. However, we show explicitly that this mode coupling cancels out when perturbations are evaluated at a fixed physical scale rather than fixed coordinate scale. Furthermore, we show how the absence of scale-dependent bias can be derived easily in any gauge. This result can then be incorporated into a complete description of the observed galaxy clustering, including the previously studied general relativistic terms, which are important at the same level as scale-dependent bias of order fNLlocal ∼ 1. This description will allow us to draw unbiased conclusions about inflation from future galaxy clustering data.

The large-scale properties of simulated cosmological magnetic fields

We perform uniformly sampled large-scale cosmological simulations including magnetic fields with the moving mesh code AREPO. We run two sets of MHD simulations: one including adiabatic gas physics only; the other featuring the fiducial feedback model of the Illustris simulation. In the adiabatic case, the magnetic field amplification follows the $B \propto \rho^{2/3}$ scaling derived from `flux-freezing' arguments, with the seed field strength providing an overall normalization factor. At high baryon overdensities the amplification is enhanced by shear flows and turbulence. Feedback physics and the inclusion of radiative cooling change this picture dramatically. In haloes, gas collapses to much larger densities and the magnetic field is amplified strongly and to the same maximum intensity irrespective of the initial seed field of which any memory is lost. At lower densities a dependence on the seed field strength and orientation, which in principle can be used to constrain models of cosmic magnetogenesis, is still present. Inside the most massive haloes magnetic fields reach values of $\sim 10-100\,\,{\rm \mu G}$, in agreement with galaxy cluster observations. The topology of the field is tangled and gives rise to rotation measure signals in reasonable agreement with the observations. However, the rotation measure signal declines too rapidly towards larger radii as compared to observational data.

NON-UNIVERSALITY OF DARK-MATTER HALOS: CUSPS, CORES, AND THE CENTRAL POTENTIAL

Dark-matter halos grown in cosmological simulations appear to have central NFW-like density cusps with mean values of $d\log\rho/d\log r \approx -1$, and some dispersion, which is generally parametrized by the varying index $\alpha$ in the Einasto density profile fitting function. Non-universality in profile shapes is also seen in observed galaxy clusters and possibly dwarf galaxies. Here we show that non-universality, at any given mass scale, is an intrinsic property of DARKexp, a theoretically derived model for collisionless self-gravitating systems. We demonstrate that DARKexp - which has only one shape parameter, $\phi_0$ - fits the dispersion in profile shapes of massive simulated halos as well as observed clusters very well. DARKexp also allows for cored dark-matter profiles, such as those found for dwarf spheroidal galaxies. We provide approximate analytical relations between DARKexp $\phi_0$, Einasto $\alpha$, or the central logarithmic slope in the Dehnen-Tremaine analytical $\gamma$-models. The range in halo parameters reflects a substantial variation in the binding energies per unit mass of dark-matter halos.

Heavy concerns about the light axino explanation of the 3.5 keV X-ray line

An unidentified 3.5 keV line from X-ray observations of galaxy clusters has been reported recently. Although still under scrutiny, decaying dark matter could be responsible for this signal. We investigate whether an axino with a mass of 7 keV could explain the line, keeping the discussion as model independent as possible. We point out several obstacles, which were overlooked in the literature, and which make the axino an unlikely candidate. The only viable scenario predicts a light metastable neutralino, with a mass between 0.1 and 10 GeV and a lifetime between 10−3 and 104 s.

MASS ACCRETION AND ITS EFFECTS ON THE SELF-SIMILARITY OF GAS PROFILES IN THE OUTSKIRTS OF GALAXY CLUSTERS

Galaxy clusters exhibit remarkable self-similar behavior which allows us to establish simple scaling relationships between observable quantities and cluster masses, making galaxy clusters useful cosmological probes. Recent X-ray observations suggest that self-similarity may be broken in the outskirts of galaxy clusters. In this work, we analyze a mass-limited sample of massive galaxy clusters from the Omega500 cosmological hydrodynamic simulation to investigate the self-similarity of the diffuse X-ray emitting intracluster medium (ICM) in the outskirts of galaxy clusters. We find that the self-similarity of the outer ICM profiles is better preserved if they are normalized with respect to the mean density of the universe, while the inner profiles are more self-similar when normalized using the critical density. However, the outer ICM profiles as well as the location of accretion shock around clusters are sensitive to their mass accretion rate, which causes the apparent breaking of self-similarity in cluster outskirts. We also find that the collisional gas does not follow the distribution of collisionless dark matter perfectly in the infall regions of galaxy clusters, leading to 10% departures in the gas-to-dark matter density ratio from the cosmic mean value. Our results have a number implications for interpreting observations of galaxy clusters in X-ray and through the Sunyaev-Zel'dovich effect and their application to cluster cosmology.

THE SHAPE OF X-RAY CAVITIES IN GALAXY CLUSTERS: PROBING JET PROPERTIES AND VISCOSITY

X-ray observations of galaxy clusters have detected numerous X-ray cavities, evolved from the interaction of AGN jets with the intracluster medium (ICM) and providing compelling evidence for the importance of jet-mode AGN feedback. Here we argue for the physical importance of the cavity shape, which we characterize with two geometric parameters: radial elongation $\tau$ and top wideness $b$. We study the cavity shape with 16 hydrodynamic jet simulations in two representative clusters, and find that the shapes of young cavities are mainly determined by various jet properties. Our simulations successfully reproduce two observed types of young cavities elongated along either the jet ($\tau>1$; type-II) or perpendicular ($\tau\leq1$; type-I) direction. Bottom-wide type-I cavities are produced by very light internally-subsonic jets, while top-wide type-II cavities are produced by heavier, internally-supersonic jets, which may also produce center-wide cavities with $\tau\sim 1$ if the jets are only slightly supersonic. Bottom-wide type-II cavities can be produced by very light jets with very long durations and cylindrical cavities are produced by very light internally-supersonic jets. While not appreciably affecting the shapes of young cavities, viscosity significantly affects the long-term cavity evolution, suppressing both interface instabilities and the formation of torus-like morphology. We encourage observers to study the shapes of young and old X-ray cavities separately, the former probing the properties of AGN jets and the latter potentially probing the ICM viscosity level.

SALT spectroscopic observations of galaxy clusters detected by ACT and a type II quasar hosted by a brightest cluster galaxy

We present Southern African Large Telescope (SALT) follow-up observations of seven massive clusters detected by the Atacama Cosmology Telescope (ACT) on the celestial equator using the Sunyaev-Zel'dovich (SZ) effect. We conducted multi-object spectroscopic observations with the Robert Stobie Spectrograph in order to measure galaxy redshifts in each cluster field, determine the cluster line-of-sight velocity dispersions, and infer the cluster dynamical masses. We find that the clusters, which span the redshift range 0.3 < z < 0.55, range in mass from (5 -- 20) x 10$^{14}$ solar masses (M200c). Their masses, given their SZ signals, are similar to those of southern hemisphere ACT clusters previously observed using Gemini and the VLT. We note that the brightest cluster galaxy in one of the systems studied, ACT-CL J0320.4+0032 at z = 0.38, hosts a Type II quasar. Only a handful of such systems are currently known, and therefore ACT-CL J0320.4+0032 may be a rare example of a very massive halo in which quasar-mode feedback is actively taking place.

Multiplicative errors in the galaxy power spectrum: self-calibration of unknown photometric systematics for precision cosmology

We develop a general method to ‘self-calibrate’ observations of galaxy clustering with respect to systematics associated with photometric calibration errors. We first point out the danger posed by the multiplicative effect of calibration errors, where large-angle error propagates to small scales and may be significant even if the large-scale information is cleaned or not used in the cosmological analysis. We then propose a method to measure the arbitrary large-scale calibration errors and use these measurements to correct the small-scale (high-multipole) power which is most useful for constraining the majority of cosmological parameters. We demonstrate the effectiveness of our approach on synthetic examples and briefly discuss how it may be applied to real data.

An algorithm to build mock galaxy catalogues using MICE simulations

We present a method to build mock galaxy catalogues starting from a halo catalogue that uses halo occupation distribution (HOD) recipes as well as the subhalo abundance matching (SHAM) technique. Combining both prescriptions we are able to push the absolute magnitude of the resulting catalogue to fainter luminosities than using just the SHAM technique and can interpret our results in terms of the HOD modelling. We optimize the method by populating with galaxies friends-of-friends dark matter haloes extracted from the Marenostrum Institut de Ci\`{e}ncies de l'Espai (MICE) dark matter simulations and comparing them to observational constraints. Our resulting mock galaxy catalogues manage to reproduce the observed local galaxy luminosity function and the colour-magnitude distribution as observed by the Sloan Digital Sky Survey. They also reproduce the observed galaxy clustering properties as a function of luminosity and colour. In order to achieve that, the algorithm also includes scatter in the halo mass - galaxy luminosity relation derived from direct SHAM and a modified NFW mass density profile to place satellite galaxies in their host dark matter haloes. Improving on general usage of the HOD that fits the clustering for given magnitude limited samples, our catalogues are constructed to fit observations at all luminosities considered and therefore for any luminosity subsample. Overall, our algorithm is an economic procedure of obtaining galaxy mock catalogues down to faint magnitudes that are necessary to understand and interpret galaxy surveys.

Discovery of O vii line emitting gas in elliptical galaxies

In the cores of ellipticals, clusters, and groups of galaxies, the gas has a cooling time shorter than 1 Gyr. It is possible to probe cooling flows through the detection of Fe XVII and O VII emission lines, but so far O VII has not been detected in any individual object. The Reflection Grating Spectrometers (RGS) aboard XMM-Newton are currently the only instruments able to detect O VII in extended objects such as elliptical galaxies and galaxy clusters. We searched for evidence of O VII through all the archival RGS observations of galaxy clusters, groups of galaxies, and elliptical galaxies focusing on those with core temperatures below 1 keV. We have discovered O VII resonance (21.6A) and forbidden (22.1A) lines for the first time in the spectra of individual objects. O VII was detected at a level higher than three sigma in six elliptical galaxies: M 84, M 86, M 89, NGC 1316, NGC 4636, and NGC 5846. M 84, M 86 and M 89 are members of the Virgo Cluster, the others are central dominant galaxies of groups, and most them show evidence of O VI in UV spectra. We detect no significant trend between the Fe XVII and O VII resonance-to-forbidden line ratios, possibly because of the limited statistics. The observed line ratios <Fe(r/f), O(r/f)>= (0.52+/-0.02, 0.9+/-0.2) indicate that the spectra of all these ellipticals are affected by resonance scattering, suggesting low turbulence. Deeper exposures will help to understand whether the Fe XVII and O VII lines are both produced by the same cooling gas or by multiphase gas. Our O VII luminosities correspond to 0.2-2 Msun/yr, which agree with the predictions for ellipticals. Such weak cooling rates would not be detected in clusters because their spectra are dominated by the emission of hotter gas, and owing to their greater distance, the expected O VII line flux would be undetectable.

Bayes-x: a Bayesian inference tool for the analysis of X-ray observations of galaxy clusters

We present the first public release of our Bayesian inference tool, Bayes-X, for the analysis of X-ray observations of galaxy clusters. We illustrate the use of Bayes-X by analysing a set of four simulated clusters at z=0.2-0.9 as they would be observed by a Chandra-like X-ray observatory. In both the simulations and the analysis pipeline we assume that the dark matter density follows a spherically-symmetric Navarro, Frenk and White (NFW) profile and that the gas pressure is described by a generalised NFW (GNFW) profile. We then perform four sets of analyses. By numerically exploring the joint probability distribution of the cluster parameters given simulated Chandra-like data, we show that the model and analysis technique can robustly return the simulated cluster input quantities, constrain the cluster physical parameters and reveal the degeneracies among the model parameters and cluster physical parameters. We then analyse Chandra data on the nearby cluster, A262, and derive the cluster physical profiles. To illustrate the performance of the Bayesian model selection, we also carried out analyses assuming an Einasto profile for the matter density and calculated the Bayes factor. The results of the model selection analyses for the simulated data favour the NFW model as expected. However, we find that the Einasto profile is preferred in the analysis of A262. The Bayes-X software, which is implemented in Fortran 90, is available at http://www.mrao.cam.ac.uk/facilities/software/bayesx/.