Galaxy Clusters Research Articles

Detecting relativistic Doppler in galaxy clustering with tailored galaxy samples

We present a method to obtain a high-significance detection of relativistic effects on cosmological scales. Measurements of such effects would be instrumental for our understanding of the Universe, as they would provide a further confirmation of the validity of general relativity as the correct description of the gravitational interaction, in a regime very far from that of strong gravity, where it has been tested to exquisite accuracy. Despite its relevance, the detection of relativistic effects has hitherto eluded us, mainly because they are stronger on the largest cosmic scales, plagued by cosmic variance. Our work focuses on the cosmological probe of galaxy clustering, describing the excess probability of finding pairs of galaxies at a given separation due to them being part of the same underlying cosmic large-scale structure. We focus on the two-point correlation function of the distribution of galaxies in Fourier space — the power spectrum — where relativistic effects appear as an imaginary contribution to the real power spectrum. By carefully tailoring cuts in magnitude/luminosity, we are able to obtain two samples (bright and faint) of the same galaxy population, whose cross-correlation power spectrum allows for a detection of the relativistic contribution. In particular, we optimise the definition of the samples to maximise the detection significance of the relativistic Doppler term for both a low-z Bright Galaxy Sample and a high-z Hα emission line galaxy population.

Vortex-p: A Helmholtz-Hodge and Reynolds decomposition algorithm for particle-based simulations

Astrophysical turbulent flows display an intrinsically multi-scale nature, making their numerical simulation and the subsequent analyses of simulated data a complex problem. In particular, two fundamental steps in the study of turbulent velocity fields are the Helmholtz-Hodge decomposition (compressive+solenoidal; HHD) and the Reynolds decomposition (bulk+turbulent; RD). These problems are relatively simple to perform numerically for uniformly-sampled data, such as the one emerging from Eulerian, fix-grid simulations; but their computation is remarkably more complex in the case of non-uniformly sampled data, such as the one stemming from particle-based or meshless simulations. In this paper, we describe, implement and test vortex-p, a publicly available tool evolved from the vortex code, to perform both these decompositions upon the velocity fields of particle-based simulations, either from smoothed particle hydrodynamics (SPH), moving-mesh or meshless codes. The algorithm relies on the creation of an ad-hoc adaptive mesh refinement (AMR) set of grids, on which the input velocity field is represented. HHD is then addressed by means of elliptic solvers, while for the RD we adapt an iterative, multi-scale filter. We perform a series of idealised tests to assess the accuracy, convergence and scaling of the code. Finally, we present some applications of the code to various SPH and meshless finite-mass (MFM) simulations of galaxy clusters performed with OpenGadget3, with different resolutions and physics, to showcase the capabilities of the code.

LtU-ILI: An All-in-One Framework for Implicit Inference in Astrophysics and Cosmology

This paper presents the Learning the Universe Implicit Likelihood Inference (LtU-ILI) pipeline, a codebase for rapid, user-friendly, and cutting-edge machine learning (ML) inference in astrophysics and cosmology. The pipeline includes software for implementing various neural architectures, training schema, priors, and density estimators in a manner easily adaptable to any research workflow. It includes comprehensive validation metrics to assess posterior estimate coverage, enhancing the reliability of inferred results. Additionally, the pipeline is easily parallelizable, designed for efficient exploration of modeling hyperparameters. To demonstrate its capabilities, we present real applications across a range of astrophysics and cosmology problems, such as: estimating galaxy cluster masses from X-ray photometry; inferring cosmology from matter power spectra and halo point clouds; characterising progenitors in gravitational wave signals; capturing physical dust parameters from galaxy colors and luminosities; and establishing properties of semi-analytic models of galaxy formation. We also include exhaustive benchmarking and comparisons of all implemented methods as well as discussions about the challenges and pitfalls of ML inference in astronomical sciences. All code and examples are made publicly available at https://github.com/maho3/ltu-ili.

Fast supermassive black hole growth in the SPT2349--56 protocluster at z=4.3

Large-scale environment is one of the main physical drivers of galaxy evolution. The densest regions at high redshifts (i.e. $z>2$ protoclusters) are gas-rich regions characterised by high star formation activity. The same physical properties that enhance star formation in protoclusters are also thought to boost the growth of supermassive black holes (SMBHs), most likely in heavily obscured conditions. We aim to test this scenario by probing the active galactic nucleus (AGN) content of SPT2349--56:\ a massive, gas-rich, and highly star-forming protocluster core at $z=4.3$ discovered as an overdensity of dusty star-forming galaxies (DSFGs). We compare our results with data on the field environment and other protoclusters We observed SPT2349--56 with Chandra (200 ks) and searched for X-ray emission from the known galaxy members. We also performed a spectral energy distribution fitting procedure to derive the physical properties of the discovered AGNs. In the X-ray band, we detected two protocluster members:\ C1 and C6, corresponding to an AGN fraction among DSFGs in the structure of $ This value is consistent with other protoclusters at $z=2-4$, but higher than the AGN incidence among DSFGs in the field environment. Both AGNs are heavily obscured sources, hosted in star-forming galaxies with $ M_ odot $ stellar masses. We estimate that the intergalactic medium in the host galaxies contributes to a significant fraction (or even entirely) to the nuclear obscuration. In particular, C1 is a highly luminous ($L_X=2 and Compton-thick ($N_H=2 cm^ $) AGN, likely powered by a $M_ BH M_ odot $ SMBH, assuming Eddington-limited accretion. Its high accretion rate suggests that it is in the phase of efficient growth that is generally required to explain the presence of extremely massive SMBHs in the centres of local galaxy clusters. Considering SPT2349--56 and DRC, a similar protocuster at $z=4$, and under different assumptions on their volumes, we find that gas-rich protocluster cores at $z enhance the triggering of luminous (log$ L_X lunit =45-46$) AGNs by three to five orders of magnitude with respect to the predictions from the AGN X-ray luminosity function at a similar redshift in the field environment. We note that this result is not solely driven by the overdensity of the galaxy population in the structures. Our results indicate that gas-rich protoclusters at high redshift boost the growth of SMBHs, which will likely impact the subsequent evolution of the structures. Therefore, they stand as key science targets to obtain a complete understanding of the relation between the environment and galaxy evolution. Dedicated investigations of similar protoclusters are required to definitively confirm this conclusion with a higher statistical significance.

MeerKAT observations of starburst galaxies and AGNs within the core of XMMXCS J2215.9−1738 at z = 1.46

ABSTRACT We present the first detailed radio study of the galaxy cluster XMMXCS J2215.9−1738 at z = 1.46 using MeerKAT L-band (1.3 GHz) observations. We combine our radio observation with archival optical and infrared data to investigate the star formation and active galactic nucleus (AGN) population within $R_{200}$ ($R =$ 0.8 Mpc) of the cluster centre. Using three selection criteria; the radio luminosity, the far-infrared radio ratio ($q_{\rm {IR}}$), and the mid-infrared colour, we distinguish galaxies with radio emission predominantly powered by star formation from that powered by AGNs. We selected 24 cluster members within $R_{\rm {200}}$ in the MeerKAT image based on either their photometric or spectroscopic redshift. We classified 12/24 (50 per cent) as galaxies whose radio emission is dominated by star-formation activity, 6/24 (25 per cent) as intermediate star-forming galaxies, and 6/24 (25 per cent) as AGN-dominated galaxies. Using the radio continuum luminosities of the star-forming cluster galaxies, we estimated an integrated star formation rate (SFR) value of 1700 $\pm$ 330 M$_{\odot }$ yr$^{-1}$ within $R_{200}$. We derived a mass-normalized integrated SFR value of $(570 \pm 110) \times 10^{-14}$ yr$^{-1}$. This supports previous observational and theoretical studies that indicated a rapid increase in star formation activity within the core of high-redshift clusters. We also show that the high-AGN fraction within the cluster core is consistent with previous cluster observations at $z \gt $ 1.5.

The infall region as a complementary probe to cluster abundance

ABSTRACT Galaxy cluster abundance measurements provide a classic test of cosmology. They are most sensitive to the evolved amplitude of fluctuations, usually expressed as $S_8 = \sigma _8\sqrt{\Omega _{\rm m}/0.3}$. Thus, abundance constraints exhibit a strong degeneracy between $\sigma _8$ and $\Omega _{\rm m}$, as do other similar low-redshift tests such as cosmic shear. The mass distribution in the infall region around galaxy clusters, where material is being accreted from the surrounding field, also exhibits a cosmological dependence, but in this case it is nearly orthogonal to the $S_8$ direction in the $\Omega _{\rm m}$–$\sigma _8$ plane, making it highly complementary to halo abundance or cosmic shear studies. We explore how weak-lensing measurements of the infall region might be used to complement abundance studies, considering three different tests. The splashback radius is a prominent feature of the infall region; we show that detection of this feature in lensing data from the Euclid survey could independently constrain $\Omega _{\rm m}$ and $\sigma _8$ to $\pm 0.05$. Another feature, the depletion radius where the bias reaches a minimum, also shows cosmological dependence, though it is challenging to observe in practice. The strongest constraints come from direct measurements of the shear profile in the infall region at 2–$4\, r_{200{\rm c}}$. Combining the latter with abundance constraints such as those reported from SRG$/$eROSITA should reduce the area of the error contours by an estimated factor of 1.2 using a sample of clusters observed by the UNIONS survey, or a factor of 3 using clusters observed by the Euclid Wide survey over a broader range of redshift.

Inferring intrahalo light from stellar kinematics. A deep learning approach

In the context of structure formation, disentangling the central galaxy stellar population from the stellar intrahalo light can help us shed light on the formation history of the halo as a whole, as the properties of the stellar components are expected to retain traces of the formation history. Many approaches are adopted to assess the task, depending on different physical assumptions (e.g. the light profile, chemical composition, and kinematical differences) and depending on whether the full six-dimensional phase-space information is known (much like in simulations) or whether one analyses projected quantities (i.e. observations). This paper paves the way for a new approach to bridge the gap between observational and simulation methods. We propose the use of projected kinematical information from stars in simulations in combination with deep learning to create a robust method for identifying intrahalo light in observational data to enhance understanding and consistency in studying the process of galaxy formation. Using deep learning techniques, particularly a convolutional neural network called U-Net, we developed a methodology for predicting these contributions in simulated galaxy cluster images. We created a sample of mock images from hydrodynamical simulations (including masking of the interlopers) to train, validate and test the network. Reinforced training (Attention U-Net) was used to improve the first results, as the innermost central regions of the mock images consistently overestimate the stellar intrahalo contribution. Our work shows that adequate training over a representative sample of mock images can lead to good predictions of the intrahalo light distribution. The model is mildly dependent on the training size and its predictions are less accurate when applied to mock images from different simulations. However, the main features (spatial scales and gradients of the stellar fractions) are recovered for all tests. While the method presented here should be considered as a proof of concept, future work (e.g. generating more realistic mock observations) is required to enable the application of the proposed model to observational data.

Merge and strip: Dark matter-free dwarf galaxies in clusters can be formed by galaxy mergers

Context. Recent observations of galaxy mergers inside galaxy cluster environments, such as NGC 5291 in the vicinity of Abell 3574, report high star formation rates in the ejected tidal tails, which point towards currently developing tidal dwarf galaxies. This prompts the intriguing question whether these newly formed stellar structures could get stripped from the galaxy potential by the cluster and thus populate it with dwarf galaxies. Aims. We verify whether environmental stripping of tidal dwarf galaxies from galaxy mergers inside galaxy cluster environments is a possible evolutionary channel to populate a galaxy cluster with low-mass and low surface brightness galaxies. Methods. We performed three high-resolution hydrodynamical simulations of mergers between spiral galaxies in a cluster environment, implementing a stellar mass ratio of 2:1 with M200 = 9.5 × 1011 M⊙ for the more massive galaxy. Between the three different simulations, we varied the initial orbit of the infalling galaxies with respect to the cluster center. Results. We demonstrate that cluster environments are capable of stripping tidal dwarf galaxies from the host potential independently of the infall orbit of the merging galaxy pair, without instantly destroying the tidal dwarfs. Starting to evolve separately from their progenitor, these newly formed dwarf galaxies reach total masses of Mtot ≈ 107 − 9 M⊙ within the limits of our resolution. In the three tested orbit scenarios, we find three, seven, and eight tidal dwarf galaxies per merger, respectively, which survive longer than 1 Gyr after the merger event. Exposed to ram pressure, these gas dominated dwarf galaxies exhibit high star formation rates while also losing gas to the environment. Experiencing a strong headwind due to their motion through the intracluster medium, they quickly lose momentum and start spiraling towards the cluster center, reaching distances on the order of 1 Mpc from their progenitor. About 4 Gyr after the merger event, we still find three and four intact dwarf galaxies in two of the tested scenarios, respectively. The other stripped tidal dwarf galaxies either evaporate in the hostile cluster environment due to their low initial mass, or are disrupted as soon as they reach the cluster center. Conclusions. The dwarf production rate due to galaxy mergers is elevated when the interaction with a cluster environment is taken into account. Comparing their contribution to the observed galaxy mass function in clusters, our results indicate that ∼30% of dwarf galaxies in clusters could have been formed by stripping from galaxy mergers.

Simulating the LOcal Web (SLOW)

Context. This is the second paper in a series presenting the results from a 500 h−1Mpc large constrained simulation of the local Universe (SLOW). The initial conditions for this cosmological hydro-dynamical simulation are based on peculiar velocities derived from the CosmicFlows-2 catalog. The simulation follows cooling, star formation, and the evolution of super-massive black holes. This allows one to directly predict observable properties of the intracluster medium (ICM) within galaxy clusters, including X-ray luminosity, temperatures, and the Compton-y signal. Aims. Comparing the properties of observed galaxy clusters within the local Universe with the properties of their simulated counterparts enables us to assess the effectiveness of the initial condition constraints in accurately replicating the mildly nonlinear properties of the largest, collapsed objects within the simulation. Methods. Based on the combination of several, publicly available surveys we compiled a sample of galaxy clusters within the local Universe, of which we were able to cross-identify 46 of them with an associated counterpart within the SLOW simulation. We then derived the probability of the cross identification based on mass, X-ray luminosity, temperature, and Compton-y by comparing it to a random selection. Results. Our set of 46 cross-identified local Universe clusters contains the 13 most massive clusters from the Planck SZ catalog as well as 70% of clusters with M500 larger than 2 × 1014 M⊙. Compared to previous constrained simulations of the local volume, we found in SLOW a much larger amount of replicated galaxy clusters, where their simulation-based mass prediction falls within the uncertainties of the observational mass estimates. Comparing the median observed and simulated masses of our cross-identified sample allows us to independently deduce a hydrostatic mass bias of (1 − b)≈0.87. Conclusions. The SLOW constrained simulation of the local Universe faithfully reproduces numerous fundamental characteristics of a sizable number of galaxy clusters within our local neighborhood, opening a new avenue for studying the formation and evolution of a large set of individual galaxy clusters as well as testing our understanding of physical processes governing the ICM.

Sensitive 3 mm Imaging of Discrete Sources in the Fields of Thermal Sunyaev–Zel’dovich Effect–Selected Galaxy Clusters

In this paper, we present the results of a blind survey for compact sources in 243 Galaxy clusters that were identified using the thermal Sunyaev–Zel'dovich effect (tSZ). The survey was carried out at 90 GHz using MUSTANG2 on the Green Bank Telescope and achieved a 5σ detection limit of 1 mJy in the center of each cluster. We detected 24 discrete sources. The majority (18) of these correspond to known radio sources, and of these, five show signs of significant variability, either with time or in spectral index. The remaining sources have no clear counterparts at other wavelengths. Searches for galaxy clusters via the tSZ strongly rely on observations at 90 GHz, and the sources found have the potential to bias mass estimates of clusters. We compare our results to the Websky simulation that can be used to estimate the source contamination in galaxy cluster catalogs. While the simulation shows a good match to our observations at the clusters’ centers, it does not match our source distribution further out. Sources over 104″ from a cluster’s center bias the tSZ signal high, for some of the sources found, by over 50%. When averaged over the whole cluster population, the effect is smaller but still at a level of 1%–2%. We also discovered that unlike previous measurements and simulations, we see an enhancement of source counts in the outer regions of the clusters and fewer sources than expected in the centers of this tSZ-selected sample.

Dianoga SIDM: Galaxy cluster self-interacting dark matter simulations

Context. Self-interacting dark matter (SIDM) can tackle or alleviate small-scale issues within the cosmological standard model ΛCDM, and diverse flavours of SIDM can produce unique astrophysical predictions, resulting in different possible signatures which can be used to test these models with dedicated observations of galaxy clusters. Aims. This work aims to assess the impact of dark matter self-interactions on the properties of galaxy clusters. In particular, the goal is to study the angular dependence of the cross section by testing rare (large angle scattering) and frequent (small angle scattering) SIDM models with velocity-dependent cross sections. Methods. We re-simulated six galaxy cluster zoom-in initial conditions with a dark matter-only run and with full-physics set-up simulations that include a self-consistent treatment of baryon physics. We tested the dark matter-only setup and the full physics setup with the collisionless cold dark matter, rare self-interacting dark matter, and frequent self-interacting dark matter models. We then studied their matter density profiles as well as their subhalo population. Results. Our dark matter-only SIDM simulations agree with theoretical models, and when baryons are included in simulations, our SIDM models substantially increase the central density of galaxy cluster cores compared to full-physics simulations using collisionless dark matter. SIDM subhalo suppression in full-physics simulations is milder compared to the one found in the dark matter-only simulations because of the cuspier baryonic potential that prevents subhalo disruption. Moreover, SIDM with small-angle scattering significantly suppresses a larger number of subhaloes compared to large-angle scattering SIDM models. Additionally, SIDM models generate a broader range of subhalo concentration values, including a tail of more diffuse subhaloes in the outskirts of galaxy clusters and a population of more compact subhaloes in the cluster cores.

The SRG/eROSITA All-Sky Survey

Aims. Characterising galaxy cluster populations from a catalogue of sources selected in astronomical surveys requires knowledge of sample incompleteness, known as the selection function. The first All-Sky Survey (eRASS1) by eROSITA on board Spectrum Roentgen Gamma (SRG) has enabled the collection of large samples of galaxy clusters detected in the soft X-ray band over the western Galactic hemisphere. The driving goal consists in constraining cosmological parameters, which puts stringent requirements on the accuracy and flexibility of explainable selection function models. Methods. We used a large set of mock observations of the eRASS1 survey and we processed simulated data identically to the real eRASS1 events. We matched detected sources to simulated clusters and we associated detections to intrinsic cluster properties. We trained a series of models to build selection functions depending only on observable surface brightness data. We developed a second series of models relying on global cluster characteristics such as X-ray luminosity, flux, and the expected instrumental count rate as well as on morphological properties. We validated our models using our simulations and we ranked them according to selected performance metrics. We validated the models with datasets of clusters detected in X-rays and via the Sunyaev–Zeldovich effect. We present the complete Bayesian population modelling framework developed for this purpose. Results. Our results reveal the surface brightness characteristics most relevant to cluster selection in the eRASS1 sample, in particular the ambiguous role of central surface brightness at the scale of the instrument resolution. We have produced a series of user-friendly selection function models and demonstrated their validity and their limitations. Our selection function for bright sources reproduces the catalogue matches with external datasets well. We discuss potential inconsistencies in the selection models at a low signal-to-noise revealed by comparison with a deep X-ray sample acquired by eROSITA during its performance verification phase. Conclusions. Detailed modelling of the eRASS1 galaxy cluster selection function is made possible by reformulating selection into a classification problem. Our models are used in the first eRASS1 cosmological analysis and in sample studies of eRASS1 cluster and groups. These models are crucial for science with eROSITA cluster samples and our new methods pave the way for further investigation of faint cluster selection effects.

Properties of Nuclear Star Clusters in Low Surface Brightness Galaxies

Using the Systematically Measuring Ultra-Diffuse Galaxies and Sloan Digital Sky Survey catalogs and our own reprocessing of the Legacy Survey imaging, we investigate the properties of nuclear star clusters (NSCs) in galaxies having central surface brightnesses as low as 27 mag arcsec−2. We identify 273 (123 with known redshift) and 32 NSC-bearing galaxies in the two samples, respectively, where we require candidate NSCs to have a separation of less than 0.10r e from the galaxy center. We find that galaxies with low central surface brightness (μ 0,g > 24 mag arcsec−2) are more likely to contain an NSC if they (1) have a higher stellar mass, (2) have a higher stellar-to-total mass ratio, (3) have a brighter central surface brightness, (4) have a larger axis ratio, or (5) lie in a denser environment. Because of the correlations among these various quantities, it is likely that only one or two are true physical drivers. We also find scaling relations for the NSC mass with stellar mass (M NSC/ M⊙=106.02±0.03(M*,gal/108 M⊙)0.77±0.04 ) and halo mass (M NSC/ M⊙=106.11±0.05(Mh,gal/1010 M⊙)0.92±0.05 ), although it is the scaling with halo mass that is consistent with a direct proportionality. In galaxies with an NSC, M NSC ≈ 10−4 M h,gal. This proportionality echoes the finding of a direct proportionality between the mass (or number) of globular clusters (GCs) in galaxies and the galaxy’s total mass. These findings favor a related origin for GCs and NSCs.

The interface of gravity and dark energy

At sufficiently large radii dark energy modifies the behavior of (a) bound orbits around a galaxy and (b) virialized gas in a cluster of galaxies. Dark energy also provides a natural cutoff to a cluster’s dark matter halo. In (a) there exists a maximum circular orbit beyond which periodic motion is no longer possible, and orbital evolution near critical binding is analytically calculable using an adiabatic invariant integral. The finding implicates the study of wide galaxy pairs. In (b), dark energy necessitates the use of a generalized Virial Theorem to describe gas at the outskirts of a cluster. When coupled to the baryonic escape condition, aided by dark energy, the results is a radius beyond which the continued establishment of a hydrostatic halo of thermalized baryons is untenable. This leads to a theoretically motivated virial radius. We use this theory to probe the structure of a cluster’s baryonic halo and apply it to X-ray and weak-lensing data collected on cluster Abell 1835. We find that gas in its outskirts deviates significantly from hydrostatic equilibrium beginning at ∼1.3 Mpc , the ‘inner’ virial radius. We also define a model dependent dark matter halo cutoff radius to A1835. The dark matter cutoff gives an upper limit to the cluster’s total mass of ∼7×1015M⊙ . Moreover, it is possible to derive an ‘outer’ hydrostatic equilibrium cutoff radius given a dark matter cutoff radius. A region of cluster gas transport and turbulence occurs between the inner and outer cutoff radii.

Can the angular scale of cosmic homogeneity be used as a cosmological test?

In standard cosmology, the cosmic homogeneity scale is the transition scale above which the patterns arising from non-uniformities – such as groups and clusters of galaxies, voids, and filaments – become indistinguishable from a random distribution of sources. Recently, different groups have investigated the feasibility of using such a scale as a cosmological test and arrived at different conclusions. In this paper, we complement and extend these studies by exploring the evolution of the spatial (RH\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$R_{\ extrm{H}}$$\\end{document}) and angular (θH\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\ heta _{\ extrm{H}}$$\\end{document}) homogeneity scales with redshift, assuming a spatially flat, Λ\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\varLambda $$\\end{document}-Cold Dark Matter universe and linear cosmological perturbation theory. We confirm previous results concerning the non-monotonicity of RH\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$R_{\ extrm{H}}$$\\end{document} with the matter density parameter Ωm0\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\varOmega _{\ extrm{m0}}$$\\end{document} but also show that it exhibits a monotonical behavior with the Hubble constant H0\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$H_{0}$$\\end{document} within a large redshift interval. More importantly, we find that, for z≳0.6\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$z \\gtrsim 0.6$$\\end{document}, θH\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\ heta _{\ extrm{H}}$$\\end{document} presents a monotonical behavior with Ωm0\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\varOmega _{\ extrm{m0}}$$\\end{document}, as well as for H0\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$H_0$$\\end{document} the entire redshift interval analyzed. We find also that the angular homogeneity scale is sensitive to H0\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$H_{0}$$\\end{document}, especially at higher redshifts. Using the currently available θH\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\ heta _{\ extrm{H}}$$\\end{document} measurements, we illustrate the constraints on the Ωm0\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\varOmega _{\ extrm{m0}}$$\\end{document}–H0\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$H_{0}$$\\end{document} plane through a MCMC analysis and show the feasibility of using the angular homogeneity scale as a new, model-independent way to constrain cosmological parameters.

Deciphering baryonic feedback with galaxy clusters

Upcoming cosmic shear analyses will precisely measure the cosmic matter distribution at low redshifts. At these redshifts, the matter distribution is affected by galaxy formation physics, primarily baryonic feedback from star formation and active galactic nuclei. Employing measurements from the Magneticum and IllustrisTNG simulations and a dark matter + baryon (DMB) halo model, this paper demonstrates that Sunyaev-Zel'dovich (SZ) effect observations of galaxy clusters, whose masses have been calibrated using weak gravitational lensing, can constrain the baryonic impact on cosmic shear with statistical and systematic errors subdominant to the measurement errors of DES-Y3 and LSST-Y1, with systematic errors on S8 and Ω m reaching 10% and 50% of the statistical errors, respectively. For LSST-Y6 and Roman surveys, these systematic errors increase to 150% and 100% of the statistical errors, indicating the necessity for further model developments for future surveys. We further dissect the contributions from different scales and halos with different masses to cosmic shear, highlighting the dominant role of SZ clusters at scales critical for cosmic shear analyses. These findings suggest a promising avenue for future joint analyses of Cosmic Microwave Background (CMB) and lensing surveys.

The heart of galaxy clusters: Demographics and physical properties of cool-core and non-cool-core halos in the TNG-Cluster simulation

We analyzed the physical properties of the gaseous intracluster medium (ICM) at the center of massive galaxy clusters with TNG-Cluster, a new cosmological magnetohydrodynamical simulation. Our sample contains 352 simulated clusters spanning a halo mass range of 1014 < M500c/M⊙ < 2 × 1015 at z = 0. We focused on the proposed classification of clusters into cool-core (CC) and non-cool-core (NCC) populations, the z = 0 distribution of cluster central ICM properties, and the redshift evolution of the CC cluster population. We analyzed the resolved structure and radial profiles of entropy, temperature, electron number density, and pressure. To distinguish between CC and NCC clusters, we considered several criteria: central cooling time, central entropy, central density, X-ray concentration parameter, and density profile slope. According to TNG-Cluster and with no a priori cluster selection, the distributions of these properties are unimodal, whereby CCs and NCCs represent the two extremes. Across the entire TNG-Cluster sample at z = 0 and based on the central cooling time, the strong CC fraction is fSCC = 24%, compared to fWCC = 60% and fNCC = 16% for weak and NCCs, respectively. However, the fraction of CCs depends strongly on both halo mass and redshift, although the magnitude and even direction of the trends vary with definition. The abundant statistics of simulated high-mass clusters in TNG-Cluster enabled us to match observational samples and make a comparison with data. The CC fractions from z = 0 to z = 2 are in broad agreement with observations, as are the radial profiles of thermodynamical quantities, globally as well as when divided as CC versus NCC halos. TNG-Cluster can therefore be used as a laboratory to study the evolution and transformations of cluster cores due to mergers, AGN feedback, and other physical processes.

The SRG/eROSITA All-Sky Survey

Context. Number counts of galaxy clusters across redshift are a powerful cosmological probe if a precise and accurate reconstruction of the underlying mass distribution is performed – a challenge called mass calibration. With the advent of wide and deep photometric surveys, weak gravitational lensing (WL) by clusters has become the method of choice for this measurement. Aims. We measured and validated the WL signature in the shape of galaxies observed in the first three years of the Dark Energy Survey (DES Y3) caused by galaxy clusters and groups selected in the first all-sky survey performed by SRG (Spectrum Roentgen Gamma)/eROSITA (eRASS1). These data were then used to determine the scaling between the X-ray photon count rate of the clusters and their halo mass and redshift. Methods. We empirically determined the degree of cluster member contamination in our background source sample. The individual cluster shear profiles were then analyzed with a Bayesian population model that self-consistently accounts for the lens sample selection and contamination and includes marginalization over a host of instrumental and astrophysical systematics. To quantify the accuracy of the mass extraction of that model, we performed mass measurements on mock cluster catalogs with realistic synthetic shear profiles. This allowed us to establish that hydrodynamical modeling uncertainties at low lens redshifts (z < 0.6) are the dominant systematic limitation. At high lens redshift, the uncertainties of the sources’ photometric redshift calibration dominate. Results. With regard to the X-ray count rate to halo mass relation, we determined its amplitude, its mass trend, the redshift evolution of the mass trend, the deviation from self-similar redshift evolution, and the intrinsic scatter around this relation. Conclusions. The mass calibration analysis performed here sets the stage for a joint analysis with the number counts of eRASS1 clusters to constrain a host of cosmological parameters. We demonstrate that WL mass calibration of galaxy clusters can be performed successfully with source galaxies whose calibration was performed primarily for cosmic shear experiments, opening the way for the cluster cosmological exploitation of future optical and NIR surveys like Euclid and LSST.

Energetic Particles and High-Energy Processes in Cosmological Filaments and Their Astronomical Implications

Large-scale cosmic filaments connect galaxies, clusters, and voids. They are permeated by magnetic fields with a variety of topologies. Cosmic rays with energies up to 1020eV can be produced in astrophysical environments associated with star-formation and AGN activities. The fate of these cosmic rays in filaments, which cannot be directly observed on Earth, are rarely studied. We investigate the high-energy processes associated with energetic particles (cosmic rays) in filaments, adopting an ecological approach that includes galaxies, clusters/superclusters, and voids as key cosmological structures in the filament ecosystem. We derive the phenomenology for modelling interfaces between filaments and these structures, and investigate how the transfer and fate of energetic cosmic ray protons are affected by the magnetism of the interfaces. We consider different magnetic field configurations in filaments and assess the implications for cosmic ray confinement and survival against hadronic pion-producing and photo-pair interactions. Our analysis shows that the fate of the particles depends on the location of their origin within a filament ecosystem, and that filaments act as ‘highways’, channelling cosmic rays between galaxies, galaxy clusters, and superclusters. Filaments can also operate as cosmic ‘fly paper’, capturing cosmic ray protons with energies up to 1018eV from cosmic voids. Our analysis predicts the presence of a population of ∼1012–1016eV cosmic ray protons in filaments and voids accumulated continually over cosmic time. These protons do not suffer significant energy losses through photo-pair or pion production, nor can they be cooled efficiently. Instead, they form a cosmic ray fossil record of the power generation history of the Universe.

A diffused background from axion-like particles in the microwave sky

The nature of dark matter is an unsolved cosmological problem and axions are one of the weakly interacting cold dark matter candidates. Axions or ALPs (Axion-like particles) are pseudo-scalar bosons predicted by beyond-standard model theories. The weak coupling of ALPs with photons leads to the conversion of CMB photons to ALPs in the presence of a transverse magnetic field. If they have the same mass as the effective mass of a photon in a plasma, the resonant conversion would cause a polarized spectral distortion leading to temperature fluctuations with the distortion spectrum. The probability of resonant conversion depends on the properties of the cluster such as the magnetic field, electron density, and its redshift. We show that this kind of conversion can happen in numerous unresolved galaxy clusters up to high redshifts, which will lead to a diffused polarised anisotropy signal in the microwave sky. The spectrum of the signal and its shape in the angular scale will be different from the lensed CMB polarization signal. This new polarised distortion spectrum will be correlated with the distribution of clusters in the universe and hence, with the large-scale structure. The spectrum can then be probed using its spectral and spatial variation with respect to the CMB and various foregrounds. An SNR of ~ 4.36 and ~ 93.87 are possible in the CMB-S4 145 GHz band and CMB-HD 150 GHz band respectively for a photon-ALPs coupling strength ofgaγ = 10-12GeV-1using galaxy clusters beyond redshiftz= 1. The same signal would lead to additional RMS fluctuations of ~7.5 × 10-2 μK at 145 GHz. In the absence of any signal, future CMB experiments such as Simons Observatory (SO), CMB-S4, and CMB-HD can put constraints on the coupling strength better than current bounds from particle physics experiment CERN Axion Solar Telescope (CAST).