Search

The intracluster medium (ICM) in the far outskirts (r $>$ 2-3 R$_{200}$) of galaxy clusters interfaces with the intergalactic medium (IGM) and is theorized to comprise diffuse, multiphase gas. This medium may hold vital clues to clusters' thermodynamic evolution and far-reaching impacts on infalling, future cluster galaxies. The diffuse outskirts of clusters are well-suited for quasar absorption line observations, capable of detecting gas to extremely low column densities. We analyze 18 QSO spectra observed with the Cosmic Origins Spectrograph aboard the Hubble Space Telescope whose lines of sight trace the gaseous environments of 26 galaxy clusters from within R$_{200}$ to 6 R$_{200}$ in projection. We measure the dN/dz and covering fraction of H I and O VI associated with the foreground clusters as a function of normalized impact parameter. We find the dN/dz for H I is consistent with the IGM field value for all impact parameter bins, with an intriguing slight elevation between 2 and 3 R$_{200}$. The dN/dz for O VI is also consistent with the field value (within 3$\sigma$) for all impact parameter bins, with potential elevations in dN/dz both within 1-2 R$_{200}$ and beyond 4 R$_{200}$ at $>2\sigma$. We propose physical scenarios that may give rise to these tentative excesses, such as a buildup of neutral gas at the outer accretion shock front and a signature of the warm-hot IGM. We do not find a systematic excess of potentially associated galaxies near the sightlines where O VI is detected; thus, the detected O VI does not have a clear circumgalactic origin., Comment: Submitted to ApJ, 33 pages, 11 primary figures, 2 tables, 25 supplementary figures

We present the Baryon Pasted (BP) X-ray and thermal Sunyaev-Zel'dovich (tSZ) maps derived from the half-sky Uchuu Lightcone simulation. These BP-Uchuu maps are constructed using more than $75$ million dark matter halos with masses $M_{500c} \geq 10^{13} M_\odot$ within the redshift range $0 \leq z \leq 2$. A distinctive feature of our BP-Uchuu Lightcone maps is their capability to assess the influence of both extrinsic and intrinsic scatter caused by triaxial gaseous halos and internal gas characteristics, respectively, at the map level. We show that triaxial gas drives substantial scatter in X-ray luminosities of clusters and groups, accounting for nearly half of the total scatter in core-excised measurements. Additionally, scatter in the thermal pressure and gas density profiles of halos enhances the X-ray and SZ power spectra, leading to biases in cosmological parameter estimates. These findings are statistically robust due to the extensive sky coverage and large halo sample in the BP-Uchuu maps., Comment: 16 pages, 12 figures, submitted to ApJ. The BP-Uchuu maps are publicly available

The recent tension in the value of the cosmological parameter $S_8 \equiv \sigma_8(\Omega_M/0.3)^{1/2}$, which represents the amplitude of the matter density fluctuations of the universe, has not been resolved. In this work, we present constraints on $S_8$ with the X-ray angular power spectra of clusters and groups measured with the half-sky map from the eROSITA All Sky Survey data release 1 (eRASS1). Thanks to the extensive sky coverage of eRASS1, its power spectrum achieves unprecedented precision compared to previous measurements. Using a well-calibrated, physical halo gas model that includes astrophysics of feedback and non-thermal pressure support, we obtain $S_8 = 0.80^{+0.02}_{-0.01}$ with 1$\sigma$ uncertainty that is competitive against other cosmological probes. Our derived $S_8$ value is smaller than the primary CMB measurements from Planck, but still consistent to within $1\sigma$. We also obtain constraints on the astrophysics of feedback, non-thermal pressure, equation of state in cluster cores, and gas boundaries in clusters and groups. We discuss how additional X-ray observations, and cosmological surveys in microwave and optical, will further improve the cosmological constraints with the angular power spectrum. Our work demonstrates that the angular power spectrum of clusters and groups is a promising probe of both cosmology and astrophysics., Comment: 12 pages, 4 figures, submitted to ApJ

The baryonic physics shaping galaxy formation and evolution are complex, spanning a vast range of scales and making them challenging to model. Cosmological simulations rely on subgrid models that produce significantly different predictions. Understanding how models of stellar and active galactic nuclei (AGN) feedback affect baryon behavior across different halo masses and redshifts is essential. Using the SIMBA and IllustrisTNG suites from the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project, we explore the effect of parameters governing the subgrid implementation of stellar and AGN feedback. We find that while IllustrisTNG shows higher cumulative feedback energy across all halos, SIMBA demonstrates a greater spread of baryons, quantified by the closure radius and circumgalactic medium (CGM) gas fraction. This suggests that feedback in SIMBA couples more effectively to baryons and drives them more efficiently within the host halo. There is evidence that different feedback modes are highly interrelated in these subgrid models. Parameters controlling stellar feedback efficiency significantly impact AGN feedback, as seen in the suppression of black hole mass growth and delayed activation of AGN feedback to higher mass halos with increasing stellar feedback efficiency in both simulations. Additionally, AGN feedback efficiency parameters affect the CGM gas fraction at low halo masses in SIMBA, hinting at complex, non-linear interactions between AGN and SNe feedback modes. Overall, we demonstrate that stellar and AGN feedback are intimately interwoven, especially at low redshift, due to subgrid implementation, resulting in halo property effects that might initially seem counterintuitive., Comment: 28 pages, 8 figures, submitted to ApJ

The relationship between galaxy size and environment has remained enigmatic, with over a decade of conflicting results. We present one of the first comprehensive studies of the variation of galaxy radius with environment beyond the local Universe and demonstrate that large-scale environmental density is correlated with galaxy radius independent of stellar mass and galaxy morphology. We confirm with $>5\sigma$ confidence that galaxies in denser environments are up to $\sim25\%$ larger than their equally massive counterparts with similar morphology in less dense regions of the Universe. We achieve this result by correlating projected two-dimensional densities over $\sim360$ deg$^2$ with the structural parameters of $\sim3$ million Hyper Suprime-Cam galaxies at $0.3 \leq z < 0.7$ with $\log M/M_{\odot} \geq 8.9$. Compared to most previous studies, this sample is $\sim100-10,000$ times larger and goes $\sim1$ dex deeper in mass-completeness. We demonstrate that past conflicting results have been driven by small sample sizes and a lack of robust measurement uncertainties. We verify the presence of the above correlation separately for disk-dominated, bulge-dominated, star-forming, and quiescent subpopulations. We find the strength of the correlation to be dependent on redshift, stellar mass, and morphology. The correlation is strongest at lower redshifts and systematically weakens or disappears beyond $z \geq 0.5$. At $z\geq0.5$, more massive galaxies still display a statistically significant correlation. Although some existing theoretical frameworks can be selectively invoked to explain some of the observed correlations, our work demonstrates the need for more comprehensive theoretical investigations of the correlation between galaxy size and environment., Comment: 29 pages, 13 figures. Published in The Astrophysical Journal. We welcome comments and constructive criticism

We present a stacked lensing analysis of 96 galaxy clusters selected by the thermal Sunyaev-Zel'dovich (SZ) effect in maps of the cosmic microwave background (CMB). We select foreground galaxy clusters with a $5\sigma$-level SZ threshold in CMB observations from the Atacama Cosmology Telescope, while we define background source galaxies for the lensing analysis with secure photometric redshift cuts in Year 3 data of the Subaru Hyper Suprime Cam survey. We detect the stacked lensing signal in the range of $0.1 < R\, [h^{-1}\mathrm{Mpc}] < 100$ in each of three cluster redshift bins, $0.092

A systematic comparison of the models of the circumgalactic medium (CGM) and their observables is crucial to understanding the predictive power of the models and constraining physical processes that affect the thermodynamics of CGM. This paper compares four analytic CGM models: precipitation, isentropic, cooling flow, and baryon pasting models for the hot, volume-filling CGM phase, all assuming hydrostatic or quasi-hydrostatic equilibrium. We show that for fiducial parameters of the CGM of a Milky-Way (MW) like galaxy ($\rm M_{vir} \sim 10^{12}~M_{\odot}$ at $z\sim 0$), the thermodynamic profiles -- entropy, density, temperature, and pressure -- show most significant differences between different models at small ($r\lesssim 30$ kpc) and large scales ($r\gtrsim 100$ kpc) while converging at intermediate scales. The slope of the entropy profile, which is one of the most important differentiators between models, is $\approx 0.8$ for the precipitation and cooling flow models, while it is $\approx0.6$ and 0 for the baryon pasting and isentropic models, respectively. We make predictions for various observational quantities for an MW mass halo for the different models, including the projected Sunyaev-Zeldovich (SZ) effect, soft X-ray emission (0.5--2 keV), dispersion measure, and column densities of oxygen ions (OVI, OVII, and OVIII) observable in absorption., Comment: Comments: 15 pages, 4 figures, 3 tables, accepted for publication in MNRAS

Galaxy cluster mergers are representative of a wide range of physics, making them an excellent probe of the properties of dark matter and the ionized plasma of the intracluster medium. To date, most studies have focused on mergers occurring in the plane of the sky, where morphological features can be readily identified. To allow study of mergers with arbitrary orientation, we have assembled multi-probe data for the eight-cluster ICM-SHOX sample sensitive to both morphology and line of sight velocity. The first ICM-SHOX paper (Silich+2023) provided an overview of our methodology applied to one member of the sample, MACS J0018.5+1626, in order to constrain its merger geometry. That work resulted in an exciting new discovery of a velocity space decoupling of its gas and dark matter distributions. In this work, we describe the availability and quality of multi-probe data for the full ICM-SHOX galaxy cluster sample. These datasets will form the observational basis of an upcoming full ICM-SHOX galaxy cluster sample analysis., Comment: 6 pages, 6 figures; published in Proc. of the mm Universe 2023 conference, EPJ Web of conferences, EDP Sciences

Galaxy formation models within cosmological hydrodynamical simulations contain numerous parameters with non-trivial influences over the resulting properties of simulated cosmic structures and galaxy populations. It is computationally challenging to sample these high dimensional parameter spaces with simulations, particularly for halos in the high-mass end of the mass function. In this work, we develop a novel sampling and reduced variance regression method, CARPoolGP, which leverages built-in correlations between samples in different locations of high dimensional parameter spaces to provide an efficient way to explore parameter space and generate low variance emulations of summary statistics. We use this method to extend the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) to include a set of 768 zoom-in simulations of halos in the mass range of $10^{13} - 10^{14.5} M_\odot\,h^{-1}$ that span a 28-dimensional parameter space in the IllustrisTNG model. With these simulations and the CARPoolGP emulation method, we explore parameter trends in the Compton $Y-M$, black hole mass-halo mass, and metallicity-mass relations, as well as thermodynamic profiles and quenched fractions of satellite galaxies. We use these emulations to provide a physical picture of the complex interplay between supernova and active galactic nuclei feedback. We then use emulations of the $Y-M$ relation of massive halos to perform Fisher forecasts on astrophysical parameters for future Sunyaev-Zeldovich observations and find a significant improvement in forecasted constraints. We publicly release both the simulation suite and CARPoolGP software package., Comment: The manuscript was submitted to arxiv after receiving and responding to comments from the first referee report

Most diffuse baryons, including the circumgalactic medium (CGM) surrounding galaxies and the intergalactic medium (IGM) in the cosmic web, remain unmeasured and unconstrained. Fast Radio Bursts (FRBs) offer an unparalleled method to measure the electron dispersion measures (DMs) of ionized baryons. Their distribution can resolve the missing baryon problem, and constrain the history of feedback theorized to impart significant energy to the CGM and IGM. We analyze the Cosmology and Astrophysics in Machine Learning (CAMEL) Simulations, using three suites: IllustrisTNG, SIMBA, and Astrid, each varying 6 parameters (2 cosmological & 4 astrophysical feedback), for a total of 183 distinct simulation models. We find significantly different predictions between the fiducial models of the suites, owing to their different implementations of feedback. SIMBA exhibits the strongest feedback, leading to the smoothest distribution of baryons, reducing the sightline-to-sightline variance in DMs between z=0-1. Astrid has the weakest feedback and the largest variance. We calculate FRB CGM measurements as a function of galaxy impact parameter, with SIMBA showing the weakest DMs due to aggressive AGN feedback and Astrid the strongest. Within each suite, the largest differences are due to varying AGN feedback. IllustrisTNG shows the most sensitivity to supernova feedback, but this is due to the change in the AGN feedback strengths, demonstrating that black holes, not stars, are most capable of redistributing baryons in the IGM and CGM. We compare our statistics directly to recent observations, paving the way for the use of FRBs to constrain the physics of galaxy formation and evolution., Comment: 15 pages, 7 figures, Accepted to ApJ

We explore the evolution of cold streams from the cosmic web that feed galaxies through their shock-heated circumgalactic medium (CGM) at cosmic noon, $z\simeq 1-5$. In addition to the hydrodynamical instabilities and radiative cooling that we have incorporated in earlier works, we embed the stream and the hot CGM in the gravitational potential of the host dark-matter halo, deriving equilibrium profiles for both. Self-gravity within the stream is tentatively ignored. We find that the cold streams gradually entrain a large mass of initially hot CGM gas that cools in the mixing layer and condenses onto the stream. This entrainment, combined with the acceleration down the gravitational potential well, typically triples the inward cold inflow rate into the central galaxy, compared to the original rate at the virial radius, which makes the entrained gas the dominant source of gas supply to the galaxy. The potential sources for the hot gas to be entrained are recycled enriched gas that has been previously ejected from the galaxy, and fresh virial-shock-heated gas that has accumulated in the CGM. This can naturally elevate the star formation rate in the galaxy by a factor of $\sim 3$ compared to the gas accretion rate onto the halo, thus explaining the otherwise puzzling observed excess of star formation at cosmic noon. When accounting for self-shielding of dense gas from the UV background, we find that the energy radiated from the streams, originating predominantly from the cooling of the entrained gas, is consistent with observed Lyman-$\alpha$ blobs around galaxies., Comment: 27 pages, 11 figures, accepted in MNRAS

We present a novel method for identifying cosmic web filaments using the IllustrisTNG (TNG100) cosmological simulations and investigate the impact of filaments on galaxies. We compare the use of cosmic density field estimates from the Delaunay Tessellation Field Estimator (DTFE) and the Monte Carlo Physarum Machine (MCPM), which is inspired by the slime mold organism, in the DisPerSE structure identification framework. The MCPM-based reconstruction identifies filaments with higher fidelity, finding more low-prominence/diffuse filaments and better tracing the true underlying matter distribution than the DTFE-based reconstruction. Using our new filament catalogs, we find that most galaxies are located within 1.5-2.5 Mpc of a filamentary spine, with little change in the median specific star formation rate and the median galactic gas fraction with distance to the nearest filament. Instead, we introduce the filament line density, Sigma_fil(MCPM), as the total MCPM overdensity per unit length of a local filament segment, and find that this parameter is a superior predictor of galactic gas supply and quenching. Our results indicate that most galaxies are quenched and gas-poor near high-line density filaments at z<=1. At z=0, quenching in log(M*/Msun)>10.5 galaxies is mainly driven by mass, while lower-mass galaxies are significantly affected by the filament line density. In high-line density filaments, satellites are strongly quenched, whereas centrals have reduced star formation, but not gas fraction, at z<=0.5. We discuss the prospect of applying our new filament identification method to galaxy surveys with SDSS, DESI, Subaru PFS, etc. to elucidate the effect of large-scale structure on galaxy formation., Comment: Accepted for publication in ApJ. Data available at https://doi.org/10.5281/zenodo.10929036

We present an investigation into a hitherto unexplored systematic that affects the accuracy of galaxy cluster mass estimates with weak gravitational lensing. Specifically, we study the covariance between the weak lensing signal, $\Delta\Sigma$, and the "true" cluster galaxy number count, $N_{\rm gal}$, as measured within a spherical volume that is void of projection effects. By quantifying the impact of this covariance on mass calibration, this work reveals a significant source of systematic uncertainty. Using the MDPL2 simulation with galaxies traced by the SAGE semi-analytic model, we measure the intrinsic property covariance between these observables within the 3D vicinity of the cluster, spanning a range of dynamical mass and redshift values relevant for optical cluster surveys. Our results reveal a negative covariance at small radial scales ($R \lesssim R_{\rm 200c}$) and a null covariance at large scales ($R \gtrsim R_{\rm 200c}$) across most mass and redshift bins. We also find that this covariance results in a $2-3\%$ bias in the halo mass estimates in most bins. Furthermore, by modeling $N_{\rm gal}$ and $\Delta\Sigma$ as multi-(log)-linear equations of secondary halo properties, we provide a quantitative explanation for the physical origin of the negative covariance at small scales. Specifically, we demonstrate that the $N_{\rm gal}$-$\Delta\Sigma$ covariance can be explained by the secondary properties of halos that probe their formation history. We attribute the difference between our results and the positive bias seen in other works with (mock)-cluster finders to projection effects. These findings highlight the importance of accounting for the covariance between observables in cluster mass estimation, which is crucial for obtaining accurate constraints on cosmological parameters.

High-resolution spectroscopy in soft X-rays will open a new window to map multiphase gas in galaxy clusters and probe physics of the intracluster medium (ICM), including chemical enrichment histories, circulation of matter and energy during large-scale structure evolution, stellar and black hole feedback, halo virialization, and gas mixing processes. An eV-level spectral resolution, large field-of-view, and effective area are essential to separate cluster emissions from the Galactic foreground and efficiently map the cluster outskirts. Several mission concepts that meet these criteria have been proposed recently, e.g., LEM, HUBS, and SuperDIOS. This theoretical study explores what information on ICM physics could be recovered with such missions and the associated challenges. We emphasize the need for a comprehensive comparison between simulations and observations to interpret the high-resolution spectroscopic observations correctly. Using Line Emission Mapper (LEM) characteristics as an example, we demonstrate that it enables the use of soft X-ray emission lines (e.g., O VII/VIII and Fe-L complex) from the cluster outskirts to measure the thermodynamic, chemical, and kinematic properties of the gas up to $r_{200}$ and beyond. By generating mock observations with full backgrounds, analysing their images/spectra with observational approaches, and comparing the recovered characteristics with true ones from simulations, we develop six key science drivers for future missions, including the exploration of multiphase gas in galaxy clusters (e.g., temperature fluctuations, phase-space distributions), metallicity, ICM gas bulk motions and turbulence power spectra, ICM-cosmic filament interactions, and advances for cluster cosmology., Comment: 24 pages, 26 figures, submitted to MNRAS. Comments are welcome

Galaxy cluster mergers are rich sources of information to test cluster astrophysics and cosmology. However, cluster mergers produce complex projected signals that are difficult to interpret physically from individual observational probes. Multi-probe constraints on the gas and dark matter cluster components are necessary to infer merger parameters that are otherwise degenerate. We present ICM-SHOX (Improved Constraints on Mergers with SZ, Hydrodynamical simulations, Optical, and X-ray), a systematic framework to jointly infer multiple merger parameters quantitatively via a pipeline that directly compares a novel combination of multi-probe observables to mock observables derived from hydrodynamical simulations. We report a first application of the ICM-SHOX pipeline to MACS J0018.5+1626, wherein we systematically examine simulated snapshots characterized by a wide range of initial parameters to constrain the MACS J0018.5+1626 merger geometry. We constrain the epoch of MACS J0018.5+1626 to the range $0$--$60$ Myr post-pericenter passage, and the viewing angle is inclined $\approx 27$--$40$ degrees from the merger axis. We obtain constraints for the impact parameter ($\lesssim 250$ kpc), mass ratio ($\approx 1.5$--$3.0$), and initial relative velocity when the clusters are separated by 3 Mpc ($\approx 1700$--3000 km s$^{-1}$). The primary and secondary clusters initially (at 3 Mpc) have gas distributions that are moderately and strongly disturbed, respectively. We discover a velocity space decoupling of the dark matter and gas distributions in MACS J0018.5+1626, traced by cluster-member galaxy velocities and the kinematic Sunyaev-Zel'dovich effect, respectively. Our simulations indicate this decoupling is dependent on the different collisional properties of the two distributions for particular merger epochs, geometries, and viewing angles., Comment: 40 pages, 18 figures; accepted for publication in ApJ

A narrow linear object extending $\sim 60 \,{\rm kpc}$ from the centre of a galaxy at redshift $z \sim 1$ has recently been discovered and interpreted as shocked gas filament forming stars. The host galaxy presents an irregular morphology, implying recent merger events. Supposing that each of the progenitor galaxies has a central supermassive black hole (SMBH) and the SMBHs are accumulated at the centre of the merger remnant, a fraction of them can be ejected from the galaxy with a high velocity due to interactions between SMBHs. When such a runaway SMBH (RSMBH) passes through the circumgalactic medium (CGM), converging flows are induced along the RSMBH path, and star formation could eventually be ignited. We show that the CGM temperature prior to the RSMBH perturbation should be below the peak temperature in the cooling function to trigger filament formation. While the gas is temporarily heated due to compression, the cooling efficiency increases, and gas accumulation becomes allowed along the path. When the CGM density is sufficiently high, the gas can cool down and develop a dense filament by $z = 1$. The mass and velocity of the RSMBH determine the scale of filament formation. Hydrodynamical simulations validate the analytical expectations. Therefore, we conclude that the perturbation by RSMBHs is a viable channel to form the observed linear object. Using the analytic model validated by simulations, we show that the CGM around the linear object to be warm ($T < 2 \times 10^5 \, K$) and dense ($n > 2 \times 10^{-5} (T/2 \times 10^5 \, K)^{-1} \, {\rm cm^{-3}}$)., Comment: corrected typos, accepted for publication in MNRAS

The circum-galactic medium (CGM) can feasibly be mapped by multiwavelength surveys covering broad swaths of the sky. With multiple large datasets becoming available in the near future, we develop a likelihood-free Deep Learning technique using convolutional neural networks (CNNs) to infer broad-scale physical properties of a galaxy's CGM and its halo mass for the first time. Using CAMELS (Cosmology and Astrophysics with MachinE Learning Simulations) data, including IllustrisTNG, SIMBA, and Astrid models, we train CNNs on Soft X-ray and 21-cm (HI) radio 2D maps to trace hot and cool gas, respectively, around galaxies, groups, and clusters. Our CNNs offer the unique ability to train and test on ''multifield'' datasets comprised of both HI and X-ray maps, providing complementary information about physical CGM properties and improved inferences. Applying eRASS:4 survey limits shows that X-ray is not powerful enough to infer individual halos with masses $\log(M_{\rm{halo}}/M_{\odot}) < 12.5$. The multifield improves the inference for all halo masses. Generally, the CNN trained and tested on Astrid (SIMBA) can most (least) accurately infer CGM properties. Cross-simulation analysis -- training on one galaxy formation model and testing on another -- highlights the challenges of developing CNNs trained on a single model to marginalize over astrophysical uncertainties and perform robust inferences on real data. The next crucial step in improving the resulting inferences on physical CGM properties hinges on our ability to interpret these deep-learning models.

We introduce a new model of the evolution of the concentration of dark matter halos, c(t). For individual halos, our model approximates c(t) as a power law with a time-dependent index, such that at early times, concentration has a nearly constant value of c=3-4, and as cosmic time progresses, c(t) smoothly increases. Using large samples of halo merger trees taken from the Bolshoi-P and MDPL2 cosmological simulations, we demonstrate that our 3-parameter model can approximate the evolution of the concentration of individual halos with a typical accuracy of 0.1 dex for t>2 Gyr for all Bolshoi-P and MDPL2 halos of present-day mass greater than 10^11.5 Msun. We additionally present a new model of the evolution of the concentration of halo populations, which we show faithfully reproduces both average concentration growth, as well as the diversity of smooth trajectories of c(t), including capturing correlations with halo mass and halo assembly history. Our publicly available source code, Diffprof, can be used to generate Monte Carlo realizations of the concentration histories of cosmologically representative halo populations; Diffprof is differentiable due to its implementation in the JAX autodiff library, which facilitates the incorporation of our model into existing analytical halo model frameworks., Comment: 11 pages, 4 appendices, v2 updated to match MNRAS publication

We quantify the cosmological spread of baryons relative to their initial neighboring dark matter distribution using thousands of state-of-the-art simulations from the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project. We show that dark matter particles spread relative to their initial neighboring distribution owing to chaotic gravitational dynamics on spatial scales comparable to their host dark matter halo. In contrast, gas in hydrodynamic simulations spreads much further from the initial neighboring dark matter owing to feedback from supernovae (SNe) and Active Galactic Nuclei (AGN). We show that large-scale baryon spread is very sensitive to model implementation details, with the fiducial \textsc{SIMBA} model spreading $\sim$40\% of baryons $>$1\,Mpc away compared to $\sim$10\% for the IllustrisTNG and \textsc{ASTRID} models. Increasing the efficiency of AGN-driven outflows greatly increases baryon spread while increasing the strength of SNe-driven winds can decrease spreading due to non-linear coupling of stellar and AGN feedback. We compare total matter power spectra between hydrodynamic and paired $N$-body simulations and demonstrate that the baryonic spread metric broadly captures the global impact of feedback on matter clustering over variations of cosmological and astrophysical parameters, initial conditions, and galaxy formation models. Using symbolic regression, we find a function that reproduces the suppression of power by feedback as a function of wave number ($k$) and baryonic spread up to $k \sim 10\,h$\,Mpc$^{-1}$ while highlighting the challenge of developing models robust to variations in galaxy formation physics implementation., Comment: 17 pages, 15 figures

The Astro2020 Decadal Survey has identified the mapping of the circumgalactic medium (CGM, gaseous plasma around galaxies) as a key objective. We explore the prospects for characterizing the CGM in and around nearby galaxy halos with a future, large grasp X-ray microcalorimeter. We create realistic mock observations from hydrodynamical simulations (EAGLE, IllustrisTNG, and Simba) that demonstrate a wide range of potential measurements, which will address the open questions in galaxy formation and evolution. By including all background and foreground components in our mock observations, we show why it is impossible to perform these measurements with current instruments, such as X-ray CCDs, and only microcalorimeters will allow us to distinguish the faint CGM emission from the bright Milky Way (MW) foreground emission lines. We find that individual halos of MW mass can, on average and depending on star formation rate, be traced out to large radii, around R500, and for larger galaxies even out to R200, using prominent emission lines, such as OVII, or OVIII. Furthermore, we show that emission line ratios for individual halos can reveal the radial temperature structure. Substructure measurements show that it will be possible to relate azimuthal variations to the feedback mode of the galaxy. We demonstrate the ability to construct temperature, velocity, and abundance ratio maps from spectral fitting for individual galaxy halos, which reveal rotation features, AGN outbursts, and enrichment., Comment: 41 pages, 19 figures, accepted for publication in ApJ

The circumgalactic medium (CGM) plays a crucial role in galaxy evolution as it fuels star formation, retains metals ejected from the galaxies, and hosts gas flows in and out of galaxies. For Milky Way-type and more massive galaxies, the bulk of the CGM is in hot phases best accessible at X-ray wavelengths. However, our understanding of the CGM remains largely unconstrained due to its tenuous nature. A promising way to probe the CGM is via X-ray absorption studies. Traditional absorption studies utilize bright background quasars, but this method probes the CGM in a pencil beam, and, due to the rarity of bright quasars, the galaxy population available for study is limited. Large-area, high spectral resolution X-ray microcalorimeters offer a new approach to exploring the CGM in emission and absorption. Here, we demonstrate that the cumulative X-ray emission from cosmic X-ray background sources can probe the CGM in absorption. We construct column density maps of major X-ray ions from the Magneticum simulation and build realistic mock images of nine galaxies to explore the detectability of X-ray absorption lines arising from the large-scale CGM. We conclude that the OVII absorption line is detectable around individual massive galaxies at the $3\sigma-6\sigma$ confidence level. For Milky Way-type galaxies, the OVII and OVIII absorption lines are detectable at the $\sim\,6\sigma$ and $\sim\,3\sigma$ levels even beyond the virial radius when co-adding data from multiple galaxies. This approach complements emission studies, does not require additional exposures, and will allow probing of the baryon budget and the CGM at the largest scales., Comment: 16 pages, 8 figures, accepted for publication in ApJ

We analyse the internal structure and dynamics of cosmic-web filaments that connect massive high-$z$ haloes. Our analysis is based on a high-resolution AREPO cosmological simulation zooming-in on a volume encompassing three ${\rm Mpc}$-scale filaments feeding three massive haloes of $\sim 10^{12}\,\text{M}_\odot$ at $z \sim 4$, embedded in a large-scale sheet. Each filament is surrounded by a cylindrical accretion shock of radius $r_{\rm shock} \sim 50 \,{\rm kpc}$. The post-shock gas is in virial equilibrium with the potential well set by an isothermal dark-matter filament. The filament line-mass is $\sim 9\times 10^8\,\text{M}_\odot\,{\rm kpc}^{-1}$, the gas fraction within $r_{\rm shock}$ is the universal baryon fraction, and the virial temperature is $\sim 7\times 10^5 {\rm K}$. In the outer ''thermal'' (T) zone, $r \geq 0.65 \, r_{\rm shock}$, inward gravity and ram-pressure forces are over-balanced by outwards thermal pressure forces, decelerating the inflowing gas expanding the shock outward. In the intermediate ''vortex'' (V) zone, $0.25 \leq r/ r_{\rm shock} \leq 0.65$, the velocity field is dominated by a quadrupolar vortex structure due to offset inflow along the sheet through the post-shock gas. The outwards force is dominated by centrifugal forces associated with these vortices, with additional contributions from global rotation and thermal pressure. The shear and turbulent forces associated with the vortices act inward. The inner ''stream'' (S) zone, $r < 0.25 \, r_{\rm shock}$, is a dense isothermal core, $T\sim 3 \times 10^4 \, {\rm K}$ and $n_{\rm H}\sim 0.01 \,{\rm cm^{-3}}$, defining the cold streams that feed galaxies. The core is formed by an isobaric cooling flow and is associated with a decrease in outwards forces, though it exhibits both inflows and outflows. [abridged], Comment: 32 pages, 22 figures, submitted to MNRAS

We investigate how cosmic web structures affect galaxy quenching in the IllustrisTNG (TNG100) cosmological simulations by reconstructing the cosmic web within each snapshot using the DisPerSE framework. We measure the comoving distance from each galaxy with stellar mass $\log(M_{\ast}/\mathrm{M}_{\odot}) \geq 8$ to the nearest node ($d_{\mathrm{node}}$) and the nearest filament spine ($d_{\mathrm{fil}}$) to study the dependence of both median specific star formation rate () and median gas fraction (<$f_{\mathrm{gas}}$>) on these distances. We find that the of galaxies is only dependent on cosmic web environment at $z<2$, with the dependence increasing with time. At $z\leq0.5$, $8 \leq \log(M_{\ast}/\mathrm{M}_{\odot}) < 9$ galaxies are quenched at $d_{\mathrm{node}}\lesssim1$~Mpc, and have significantly-suppressed star formation at $d_{\mathrm{fil}}\lesssim1$~Mpc, trends driven mostly by satellite galaxies. At $z\leq1$, in contrast to the monotonic drop in of $\log(M_{\ast}/\mathrm{M}_{\odot}) <10$ galaxies with decreasing $d_{\mathrm{node}}$ and $d_{\mathrm{fil}}$, $\log(M_{\ast}/\mathrm{M}_{\odot}) \geq 10$ galaxies - both centrals and satellites - experience an upturn in at $d_{\mathrm{node}}\lesssim0.2$~Mpc. Much of this cosmic web dependence of star formation activity can be explained by an evolution in $$. Our results suggest that in the past $\sim$10 Gyr, low-mass satellites are quenched by rapid gas stripping in dense environments near nodes and gradual gas starvation in intermediate-density environments near filaments, while at earlier times cosmic web structures efficiently channeled cold gas into most galaxies. State-of-the-art ongoing spectroscopic surveys such as SDSS and DESI, as well as those planned with the Subaru Prime Focus Spectrograph, JWST and Roman, are required to test our predictions against observations., Comment: Accepted for publication in ApJ. 5 Figures

We evaluate the effectiveness of deep learning (DL) models for reconstructing the masses of galaxy clusters using X-ray photometry data from next-generation surveys. We establish these constraints using a catalogue of realistic mock eROSITA X-ray observations which use hydrodynamical simulations to model realistic cluster morphology, background emission, telescope response, and AGN sources. Using bolometric X-ray photon maps as input, DL models achieve a predictive mass scatter of $\sigma_{\ln M_\mathrm{500c}} = 17.8\%$, a factor of two improvements on scalar observables such as richness $N_\mathrm{gal}$, 1D velocity dispersion $\sigma_\mathrm{v,1D}$, and photon count $N_\mathrm{phot}$ as well as a $32\%$ improvement upon idealised, volume-integrated measurements of the bolometric X-ray luminosity $L_X$. We then show that extending this model to handle multichannel X-ray photon maps, separated in low, medium, and high energy bands, further reduces the mass scatter to $16.2\%$. We also tested a multimodal DL model incorporating both dynamical and X-ray cluster probes and achieved marginal gains at a mass scatter of $15.9\%$. Finally, we conduct a quantitative interpretability study of our DL models and find that they greatly down-weight the importance of pixels in the centres of clusters and at the location of AGN sources, validating previous claims of DL modelling improvements and suggesting practical and theoretical benefits for using DL in X-ray mass inference., Comment: 14 pages, 9 figures, 3 tables, accepted in MNRAS

The interaction of a runaway supermassive black hole (SMBH) with the circumgalactic medium (CGM) can lead to the formation of a wake of shocked gas and young stars behind it. Here we report the serendipitous discovery of an extremely narrow linear feature in HST/ACS images that may be an example of such a wake. The feature extends 62 kpc from the nucleus of a compact star-forming galaxy at z=0.964. Keck LRIS spectra show that the [OIII]/H$\beta$ ratio varies from ~1 to ~10 along the feature, indicating a mixture of star formation and fast shocks. The feature terminates in a bright [OIII] knot with a luminosity of 1.9x10$^{41}$ ergs/s. The stellar continuum colors vary along the feature, and are well-fit by a simple model that has a monotonically increasing age with distance from the tip. The line ratios, colors, and the overall morphology are consistent with an ejected SMBH moving through the CGM at high speed while triggering star formation. The best-fit time since ejection is ~39 Myr and the implied velocity is v~1600 km/s. The feature is not perfectly straight in the HST images, and we show that the amplitude of the observed spatial variations is consistent with the runaway SMBH interpretation. Opposite the primary wake is a fainter and shorter feature, marginally detected in [OIII] and the rest-frame far-ultraviolet. This feature may be shocked gas behind a binary SMBH that was ejected at the same time as the SMBH that produced the primary wake., Comment: Accepted for publication in ApJ Letters. The key data are in Figure 1: a really odd thin streak in HST images, with a complex emission line spectrum. Figure 7 is an illustration of our proposed interpretation

The "cosmic web", the filamentary large-scale structure in a cold dark matter Universe, is readily apparent via galaxy tracers in spectroscopic surveys. However, the underlying dark matter structure is as of yet unobservable and mapping the diffuse gas permeating it lies beyond practical observational capabilities. A recently developed technique, inspired by the growth and movement of Physarum polycephalum "slime mold", has been used to map the cosmic web of a low redshift sub-sample of the SDSS spectroscopic galaxy catalog. This model, the Monte Carlo Physarum Machine (MCPM) was shown to promisingly reconstruct the cosmic web. Here, we improve the formalism used in calibrating the MCPM to better recreate the Bolshoi-Planck cosmological simulation's density distributions and apply them to a significantly larger cosmological volume than previous works using the Sloan Digital Sky Survey (SDSS, $z < 0.1$) and the Extended Baryon Oscillation Spectroscopic Survey (eBOSS) Luminous Red Galaxy (LRG, $z \lesssim 0.5$) spectroscopic catalogs. We present the "Cosmic Slime Value Added Catalog" which provides estimates for the cosmic overdensity for the sample of galaxies probed spectroscopically by the above SDSS surveys. In addition, we provide the fully reconstructed 3D density cubes of these volumes. These data products were released as part of Sloan Digital Sky Survey Data Release 17 and are publicly available. We present the input catalogs and the methodology for constructing these data products. We also highlight exciting potential applications to galaxy evolution, cosmology, the intergalactic and circumgalactic medium, and transient phenomenon localization., Comment: Submitted to ApJS, 15 pages, 9 figures

We use the Galaxy Morphology Posterior Estimation Network (GaMPEN) to estimate morphological parameters and associated uncertainties for $\sim 8$ million galaxies in the Hyper Suprime-Cam (HSC) Wide survey with $z \leq 0.75$ and $m \leq 23$. GaMPEN is a machine learning framework that estimates Bayesian posteriors for a galaxy's bulge-to-total light ratio ($L_B/L_T$), effective radius ($R_e$), and flux ($F$). By first training on simulations of galaxies and then applying transfer learning using real data, we trained GaMPEN with $<1\%$ of our dataset. This two-step process will be critical for applying machine learning algorithms to future large imaging surveys, such as the Rubin-Legacy Survey of Space and Time (LSST), the Nancy Grace Roman Space Telescope (NGRST), and Euclid. By comparing our results to those obtained using light-profile fitting, we demonstrate that GaMPEN's predicted posterior distributions are well-calibrated ($\lesssim 5\%$ deviation) and accurate. This represents a significant improvement over light profile fitting algorithms which underestimate uncertainties by as much as $\sim60\%$. For an overlapping sub-sample, we also compare the derived morphological parameters with values in two external catalogs and find that the results agree within the limits of uncertainties predicted by GaMPEN. This step also permits us to define an empirical relationship between the S\'ersic index and $L_B/L_T$ that can be used to convert between these two parameters. The catalog presented here represents a significant improvement in size ($\sim10 \times $), depth ($\sim4$ magnitudes), and uncertainty quantification over previous state-of-the-art bulge+disk decomposition catalogs. With this work, we also release GaMPEN's source code and trained models, which can be adapted to other datasets., Comment: 39 pages, 31 figures. Published in The Astrophysical Journal. We welcome comments and constructive criticism. Public Data Release at http://gampen.ghosharitra.com/

Known as the "Missing Baryon Problem", about one-third of baryons in the local universe remain unaccounted for. The missing baryons are thought to reside in the warm-hot intergalactic medium (WHIM) of the cosmic web filaments, which are challenging to detect. Recent Chandra X-ray observations used a novel stacking analysis and detected an OVII absorption line toward the sightline of a luminous quasar, hinting that the missing baryons may reside in the WHIM. To explore how the properties of the OVII absorption line depend on feedback physics, we compare the observational results with predictions obtained from the Cosmology and Astrophysics with MachinE Learning (CAMEL) Simulation suite. CAMELS consists of cosmological simulations with state-of-the-art supernova (SN) and active galactic nuclei (AGN) feedback models from the IllustrisTNG and SIMBA simulations, with varying strengths. We find that the simulated OVII column densities are higher in the outskirts of galaxies than in the large-scale WHIM, but they are consistently lower than those obtained in the Chandra observations, for all feedback runs. We establish that the OVII distribution is primarily sensitive to changes in the SN feedback prescription, whereas changes in the AGN feedback prescription have minimal impact. We also find significant differences in the OVII column densities between the IllustrisTNG and SIMBA runs. We conclude that the tension between the observed and simulated OVII column densities cannot be explained by the wide range of feedback models implemented in CAMELS., Comment: 11 pages, 9 figures. Matching to the published version on MNRAS

The Line Emission Mapper (LEM) is an X-ray Probe for the 2030s that will answer the outstanding questions of the Universe's structure formation. It will also provide transformative new observing capabilities for every area of astrophysics, and to heliophysics and planetary physics as well. LEM's main goal is a comprehensive look at the physics of galaxy formation, including stellar and black-hole feedback and flows of baryonic matter into and out of galaxies. These processes are best studied in X-rays, and emission-line mapping is the pressing need in this area. LEM will use a large microcalorimeter array/IFU, covering a 30x30' field with 10" angular resolution, to map the soft X-ray line emission from objects that constitute galactic ecosystems. These include supernova remnants, star-forming regions, superbubbles, galactic outflows (such as the Fermi/eROSITA bubbles in the Milky Way and their analogs in other galaxies), the Circumgalactic Medium in the Milky Way and other galaxies, and the Intergalactic Medium at the outskirts and beyond the confines of galaxies and clusters. LEM's 1-2 eV spectral resolution in the 0.2-2 keV band will make it possible to disentangle the faintest emission lines in those objects from the bright Milky Way foreground, providing groundbreaking measurements of the physics of these plasmas, from temperatures, densities, chemical composition to gas dynamics. While LEM's main focus is on galaxy formation, it will provide transformative capability for all classes of astrophysical objects, from the Earth's magnetosphere, planets and comets to the interstellar medium and X-ray binaries in nearby galaxies, AGN, and cooling gas in galaxy clusters. In addition to pointed observations, LEM will perform a shallow all-sky survey that will dramatically expand the discovery space., Comment: 18 pages. White paper for a mission concept to be submitted for the 2023 NASA Astrophysics Probes opportunity. v2: All-sky survey figure expanded, references fixed. v3: Added energy resolution measurements for prototype detector array. v4: Author list and reference fixes

Although weak lensing (WL) is a powerful method to estimate a galaxy cluster mass without any dynamical assumptions, a model bias can arise when the cluster density profile departs from the assumed model profile. In a merging system, the bias is expected to become most severe because the constituent halos undergo significant structural changes. In this study, we investigate WL mass bias in binary cluster mergers using a suite of idealized hydrodynamical simulations. Realistic WL shear catalogs are generated by matching the source galaxy properties, such as intrinsic shape dispersion, measurement noise, source densities, etc., to those from Subaru and {\it Hubble Space Telescope} observations. We find that, with the typical mass-concentration ($M$-$c$) relation and the Navarro-Frenk-White (NFW) profile, the halo mass bias depends on the time since the first pericenter passage and increases with the mass of the companion cluster. The time evolution of the mass bias is similar to that of the concentration, indicating that, to first order, the mass bias is modulated by the concentration change. For a collision between two $\sim10^{15}~M_{\odot}$ clusters, the maximum bias amounts to $\sim60\%$. This suggests that previous WL studies may have significantly overestimated the mass of the clusters in some of the most massive mergers. Finally, we apply our results to three merger cases: Abell 2034, MACS J1752.0+4440, and ZwCl 1856.8+6616, and report their mass biases at the observed epoch, as well as their pre-merger masses, utilizing their merger shock locations as tracers of the merger phases., Comment: 14 pages, 11 figures, accepted to ApJ

While awaiting direct velocity measurement of gas motions in the hot intracluster medium, we rely on indirect probes, including gas perturbations in galaxy clusters. Using a sample of $\sim 80$ clusters in different dynamic states from Omega500 cosmological simulations, we examine scaling relations between the fluctuation amplitudes of gas density, $\delta\rho/\rho$, pressure, $\delta P/P$, X-ray surface brightness, Sunyaev-Zeldovich (SZ) y-parameter, and the characteristic Mach number of gas motions, $M_{\rm 1d}$. In relaxed clusters, accounting for halo ellipticities reduces $\delta\rho/\rho$ or $\delta P/P$ by a factor of up to 2 within $r_{500c}$. We confirm a strong linear correlation between $\delta\rho/\rho$ (or $\delta P/P$) and $M_{\rm 1d}$ in relaxed clusters, with the proportionality coefficient $\eta \approx 1$. For unrelaxed clusters, the correlation is less strong and has a larger $\eta\approx 1.3\pm 0.5$ ($1.5\pm0.5$) for $\delta\rho/\rho$ ($\delta P/P$). Examination of the power-law scaling of $M_{\rm 1d}$ with $\delta\rho/\rho$ shows that it is almost linear for relaxed clusters, while for the unrelaxed ones, it is closer to $\delta\rho/\rho\propto M_{\rm 1d}^2$, supporting an increasing role of non-linear terms and compressive modes. In agreement with previous studies, we observe a strong correlation of $M_{\rm 1d}$ with radius. Correcting for these correlations leaves a residual scatter in $M_{\rm 1d}$ of $\sim 4 (7)$ per cent for relaxed (perturbed) clusters. Hydrostatic mass bias correlates with $M_{\rm 1d}$ as strongly as with $\delta\rho/\rho$ in relaxed clusters. The residual scatters after correcting for derived trends is $\sim 6-7$ per cent. These predictions can be verified with existing X-ray and SZ observations of galaxy clusters combined with forthcoming velocity measurements with X-ray microcalorimeters., Comment: 15 pages, 11 figures, submitted

We present the public data release of the Uchuu-UM galaxy catalogues by applying the UniverseMachine algorithm to assign galaxies to the dark matter halos in the Uchuu $N$-body cosmological simulation. It includes a variety of baryonic properties for all galaxies down to $\sim 5\times10^8 M_{\odot}$ with halos in a mass range of $10^{10}

Feedback from active galactic nuclei (AGN) and supernovae can affect measurements of integrated SZ flux of halos ($Y_\mathrm{SZ}$) from CMB surveys, and cause its relation with the halo mass ($Y_\mathrm{SZ}-M$) to deviate from the self-similar power-law prediction of the virial theorem. We perform a comprehensive study of such deviations using CAMELS, a suite of hydrodynamic simulations with extensive variations in feedback prescriptions. We use a combination of two machine learning tools (random forest and symbolic regression) to search for analogues of the $Y-M$ relation which are more robust to feedback processes for low masses ($M\lesssim 10^{14}\, h^{-1} \, M_\odot$); we find that simply replacing $Y\rightarrow Y(1+M_*/M_\mathrm{gas})$ in the relation makes it remarkably self-similar. This could serve as a robust multiwavelength mass proxy for low-mass clusters and galaxy groups. Our methodology can also be generally useful to improve the domain of validity of other astrophysical scaling relations. We also forecast that measurements of the $Y-M$ relation could provide percent-level constraints on certain combinations of feedback parameters and/or rule out a major part of the parameter space of supernova and AGN feedback models used in current state-of-the-art hydrodynamic simulations. Our results can be useful for using upcoming SZ surveys (e.g., SO, CMB-S4) and galaxy surveys (e.g., DESI and Rubin) to constrain the nature of baryonic feedback. Finally, we find that the an alternative relation, $Y-M_*$, provides complementary information on feedback than $Y-M$, Comment: Version appearing in MNRAS. Minor change to Fig.6 and added Fig. A5 compared to the previous version. 7+5 pages. The code and data associated with this paper are available at https://github.com/JayWadekar/ScalingRelations_ML

The eROSITA X-ray telescope, launched in 2019, is predicted to observe roughly 100,000 galaxy clusters. Follow-up observations of these clusters from Chandra, for example, will be needed to resolve outstanding questions about galaxy cluster physics. Deep Chandra cluster observations are expensive and follow-up of every eROSITA cluster is infeasible, therefore, objects chosen for follow-up must be chosen with care. To address this, we have developed an algorithm for predicting longer duration, background-free observations based on mock eROSITA observations. We make use of the hydrodynamic cosmological simulation Magneticum, have simulated eROSITA instrument conditions using SIXTE, and have applied a novel convolutional neural network to output a deep Chandra-like "super observation" of each cluster in our simulation sample. Any follow-up merit assessment tool should be designed with a specific use case in mind; our model produces observations that accurately and precisely reproduce the cluster morphology, which is a critical ingredient for determining cluster dynamical state and core type. Our model will advance our understanding of galaxy clusters by improving follow-up selection and demonstrates that image-to-image deep learning algorithms are a viable method for simulating realistic follow-up observations., Comment: 21 pages, 11 figures, 3 tables. Minor changes upon revision. Corrected caption of Figure 3. Added discussion of alternative asymmetry metrics. To be published in the Astrophysical Journal

Dark matter haloes have long been recognized as one of the fundamental building blocks of large scale structure formation models. Despite their importance -- or perhaps because of it! -- halo definitions continue to evolve towards more physically motivated criteria. Here, we propose a new definition that is physically motivated, and effectively unique and parameter-free: ''A dark matter halo is comprised of the collection of particles orbiting in their own self-generated potential.'' This definition is enabled by the fact that, even with as few as $\approx 300$ particles per halo, nearly every particle in the vicinity of a halo can be uniquely classified as either orbiting or infalling based on its dynamical history. For brevity, we refer to haloes selected in this way as physical haloes. We demonstrate that: 1) the mass function of physical haloes is Press-Schechter, provided the critical threshold for collapse is allowed to vary slowly with peak height; and 2) the peak-background split prediction of the clustering amplitude of physical halos is statistically consistent with the simulation data, with an accuracy no worse than $\approx 5\%$., Comment: 16 pages, 14 figures

The Sunyaev-Zeldovich (SZ) effect is a powerful tool in modern cosmology. With future observations promising ever improving SZ measurements, the relativistic corrections to the SZ signals from galaxy groups and clusters are increasingly relevant. As such, it is important to understand the differences between three temperature measures: (a) the average relativistic SZ (rSZ) temperature, (b) the mass-weighted temperature relevant for the thermal SZ (tSZ) effect, and (c) the X-ray spectroscopic temperature. In this work, we compare these cluster temperatures, as predicted by the {\sc Bahamas} \& {\sc Macsis}, {\sc Illustris-TNG}, {\sc Magneticum}, and {\sc The Three Hundred Project} simulations. Despite the wide range of simulation parameters, we find the SZ temperatures are consistent across the simulations. We estimate a $\simeq 10\%$ level correction from rSZ to clusters with $Y\simeq10^{-4}$~Mpc$^{-2}$. Our analysis confirms a systematic offset between the three temperature measures; with the rSZ temperature $\simeq 20\%$ larger than the other measures, and diverging further at higher redshifts. We demonstrate that these measures depart from simple self-similar evolution and explore how they vary with the defined radius of haloes. We investigate how different feedback prescriptions and resolution affect the observed temperatures, and discover the SZ temperatures are rather insensitive to these details. The agreement between simulations indicates an exciting avenue for observational and theoretical exploration, determining the extent of relativistic SZ corrections. We provide multiple simulation-based fits to the scaling relations for use in future SZ modelling., Comment: Accepted by MNRAS

Refereeing is a crucial component of publishing astronomical research, but few professional astronomers receive formal training on how to effectively referee a manuscript. In this article, we lay out considerations and best practices for referees. This document is intended as a tool for early career researchers to develop a fair, effective, and efficient approach to refereeing., Comment: Submitted to the Bulletin of the AAS

One of the main limitations in precision cluster cosmology arises from systematic errors and uncertainties in estimating cluster masses. Using the Mock-X pipeline, we produce synthetic X-ray images and derive cluster and galaxy group X-ray properties for a sample of over 30,000 simulated galaxy groups and clusters with $M_{\rm 500crit}$ between $10^{12}$ and $2\times 10^{15}$ M$_{\odot}$ in IllustrisTNG. We explore the similarities and differences between IllustrisTNG predictions of the Sunyaev-Zel'dovich and X-ray scaling relations with mass. We find a median hydrostatic mass bias $b = 0.125 \pm 0.003$ for $M_{\rm 500crit}$ $>10^{13}$ M$_{\odot}$. The bias increases to $b = 0.17 \pm 0.004$ when masses are derived from synthetic X-ray observations. We model how different underlying assumptions about the dependence of $Y_{\rm X}$ on halo mass can generate biases in the observed $Y_{\rm SZ} - M_{Y_{\rm X}}$ scaling relation. In particular, the simplifying assumption that $Y_{\rm X} - M_{\rm tot}$ is self-similar at all mass scales largely hides the break in $Y_{\rm SZ} - M_{\rm tot}$ and overestimates $Y_{\rm SZ}$ at galaxy and groups scales. We show that calibrating the $Y_{\rm X}-$mass proxy using a new model for a smoothly broken power law reproduces the true underlying $Y_{\rm SZ} - M_{\rm tot}$ scaling relation with high accuracy. Moreover, $M_{Y_{\rm X}}$ estimates calibrated with this method lead to $Y_{\rm SZ} - M_{Y_{\rm X}}$ predictions that are not biased by the presence of lower mass clusters or galaxy groups in the sample. Finally, we show that our smoothly broken power law model provides a robust way to derive the $Y_{\rm X}-$mass proxy, significantly reducing the level of mass bias for clusters, groups, and galaxies., Comment: 12 pages, 5 figures; submitted to MNRAS. This is a companion paper to Pop et al. 2022: "Sunyaev-Zel'dovich effect and X-ray scaling relations of galaxies, groups and clusters in the IllustrisTNG simulations"

Observable thermodynamical properties of the intracluster medium (ICM) reflect the complex interplay between AGN feedback and the gravitational collapse of haloes. Using the large volume TNG300 simulation of the IllustrisTNG project we provide predictions for X-ray and Sunyaev-Zel'dovich (SZ) scaling relations for a sample of over 30,000 haloes that cover a wide mass range from galaxies to massive galaxy clusters ($M_{\rm 500crit}$ $\in [10^{12}$ M$_{\odot} - 2\times 10^{15}$ M$_{\odot}$]). We produce mock X-ray observations of simulated haloes using methods that are consistent with observational techniques. Thus, we investigate the scaling relations between the soft-band X-ray luminosity, spectroscopic temperature, gas mass fraction, $Y_{\rm X}$ and $Y_{\rm SZ}$ as a function of halo mass, and we find broad agreement between IllustrisTNG and the observed relations. Our results highlight the scatter and bias introduced by estimated masses, and thus the importance of converting simulated ICM properties to the observable space when comparing simulations to current X-ray observations. The wide range of halo masses in our sample provides new insights into the shape of the X-ray and SZ scaling relations across three orders of magnitude in mass. Our findings show strong evidence for a break in $z=0$ scaling relations. We introduce a smoothly broken power law model which robustly captures the location of this break, the width of the transition region around the break, as well as the slope dependence on halo mass. Our results inform the next generation of subgrid black hole feedback models and provide predictions for ongoing and future observational surveys., Comment: 29 pages, 15 figures, 7 tables, 2 appendices; submitted to MNRAS. This is a companion paper to Pop et al. 2022: "Unifying Sunyaev-Zel'dovich and X-ray predictions from clusters to galaxy groups: the impact of X-ray mass estimates on the Y-M scaling relation"

The ultra-diffuse galaxies DF2 and DF4 in the NGC1052 group share several unusual properties: they both have large sizes, rich populations of overluminous and large globular clusters, and very low velocity dispersions indicating little or no dark matter. It has been suggested that these galaxies were formed in the aftermath of high velocity encounters of gas rich galaxies, events that resemble the collision that created the bullet cluster but on much smaller scales. The gas separates from the dark matter in the collision and subsequent star formation leads to the formation of one or more dark matter-free galaxies. Here we show that the present-day line-of-sight distances and radial velocities of DF2 and DF4 are consistent with their joint formation in the aftermath of a single bullet-dwarf collision, around eight billion years ago. Moreover, we find that DF2 and DF4 are part of an apparent linear substructure of 7-11 large, low-luminosity objects. We propose that these all originated in the same event, forming a trail of dark matter-free galaxies that is more than 2 Mpc long and angled 7 +- 2 degrees from the line of sight. We also tentatively identify the highly dark matter-dominated remnants of the two progenitor galaxies that are expected at the leading edges of the trail., Comment: Nature, May 19 edition. See https://www.nature.com/articles/s41586-022-04665-6 . A pretty 40-orbit HST/ACS image of DF2 (used in Shen et al. 2021) can be found at https://www.pietervandokkum.com/publications

The rotation curves of some star forming massive galaxies at redshift two decline over the radial range of a few times the effective radius, indicating a significant deficit of dark matter (DM) mass in the galaxy centre. The DM mass deficit is interpreted as the existence of a DM density core rather than the cuspy structure predicted by the standard cosmological model. A recent study proposed that a galaxy merger, in which the smaller satellite galaxy is significantly compacted by dissipative contraction of the galactic gas, can heat the centre of the host galaxy and help make a large DM core. By using an $N$-body simulation, we find that a large amount of DM mass is imported to the centre by the merging satellite, making this scenario an unlikely solution for the DM mass deficit. In this work, we consider giant baryonic clumps in high redshift galaxies as alternative heating source for creating the baryon dominated galaxies with a DM core. Due to dynamical friction, the orbit of clumps decays in a few Gyr and the baryons condensate at the galactic centre. As a back-reaction, the halo centre is heated up and the density cusp is flattened out. The combination of the baryon condensation and core formation makes the galaxy baryon dominated in the central 2-5 kpc, comparable to the effective radius of the observed galaxies. Thus, the dynamical heating by giant baryonic clumps is a viable mechanism for explaining the observed dearth of DM in high redshift galaxies., Comment: 11 pages, 7 figures, 1 table and 1 appendix, accepted for publication in MNRAS

We present Classification of Cluster GAlaxy MEmbers (C$^2$-GaMe), a classification algorithm based on a suite of machine learning models that differentiates galaxies into orbiting, infalling, and background (interloper) populations, using phase space information as input. We train and test C$^2$-GaMe with the galaxies from UniverseMachine mock catalog based on Multi-Dark Planck 2 N-body simulations. We show that probabilistic classification is superior to deterministic classification in estimating the physical properties of clusters, including density profiles and velocity dispersion. We propose a set of estimators to get an unbiased estimation of cluster properties. We demonstrate that C$^2$-GaMe can recover the distribution of orbiting and infalling galaxies' position and velocity distribution with $<1\%$ statistical error when using probabilistic predictions in the presence of interlopers in the projected phase space. Additionally, we demonstrate the robustness of trained models by applying them to a different simulation. Finally, adding a specific star formation rate and the ratio of the galaxy's halo mass to the cluster's halo mass as additional features improves the classification performance. We discuss potential applications of this technique to enhance cluster cosmology and galaxy quenching., Comment: Accepted in Astronomy and Computing. Some figures are downsampled from journal version to appear correctly in Chrome

Halos of similar mass and redshift exhibit a large degree of variability in their differential properties, such as dark matter, hot gas, and stellar mass density profiles. This variability is an indicator of diversity in the formation history of these dark matter halos that is reflected in the coupling of scatters about the mean relations. In this work, we show that the strength of this coupling depends on the scale at which halo profiles are measured. By analyzing the outputs of the IllustrisTNG hydrodynamical cosmological simulations we report the radial- and mass-dependent couplings between the dark matter, hot gas, and stellar mass radial density profiles utilizing the population diversity in dark matter halos. We find that for the same mass halos the scatters in the density of baryons and dark matter are strongly coupled at large scales ($r>R_{200}$); but the coupling between gas and dark matter density profiles fades near the core of halos ($r < 0.3 R_{200}$). We then show that the correlation between halo profile and integrated quantities induces a radius-dependent additive bias in the profile observables of halos when halos are selected on properties other than their mass. We discuss the impact of this effect on cluster abundance and cross-correlations cosmology with multi-wavelength cosmological surveys., Comment: 8 figures, 4 appendices, data in Table 2 (Accepted)

We present a percent-level accurate model of the line-of-sight velocity distribution of galaxies around dark matter halos as a function of projected radius and halo mass. The model is developed and tested using synthetic galaxy catalogs generated with the UniverseMachine run on the Multi-Dark Planck 2 N-body simulations. The model decomposes the galaxies around a cluster into three kinematically distinct classes: orbiting, infalling, and interloping galaxies. We demonstrate that: 1) we can statistically distinguish between these three types of galaxies using only projected line-of-sight velocity information; 2) the halo edge radius inferred from the line-of-sight velocity dispersion is an excellent proxy for the three-dimensional halo edge radius; 3) we can accurately recover the full velocity dispersion profile for each of the three populations of galaxies. Importantly, the velocity dispersion profiles of the orbiting and infalling galaxies contain five independent parameters -- three distinct radial scales and two velocity dispersion amplitudes -- each of which is correlated with mass. Thus, the velocity dispersion profile of galaxy clusters has inherent redundancies that allow us to perform nontrivial systematics check from a single data set. We discuss several potential applications of our new model for detecting the edge radius and constraining cosmology and astrophysics using upcoming spectroscopic surveys., Comment: 10 pages, 11 figures, accepted in MNRAS

The eROSITA Final Equatorial Depth Survey (eFEDS), with a sky area of 140 square degrees with depth equivalent to the equatorial patch of the final eROSITA all-sky survey, represents the largest continuous non-full-sky X-ray fields to-date, making it the premier data set for measuring the angular power spectrum. In this work, we measure the X-ray angular power spectrum of galaxy clusters and groups in the eFEDS field. We show that the measured power spectrum is consistent with past observations, including the ROSAT All Sky Survey, and the Chandra COSMOS and Bootes fields. The predictions of cluster gas halo model that is calibrated from Chandra observations is also consistent with the eFEDS power spectrum. While the eFEDS does not have large enough sky coverage to provide meaningful cosmological constraints, we predict that the X-ray power spectrum from the cycle 4 of the eROSITA All-Sky Survey (eRASS4) will provide constraints on $\Omega_M$ and $\sigma_8$ at the 10% level., Comment: 10 pages, 4 Figures, accepted for publication in MNRAS. Updated arxiv version matched to the published version. Major change: tensions at small-scale power are resolved after point source mask correction

This is a solicited whitepaper for the Snowmass 2021 community planning exercise. The paper focuses on measurements and science with the Cosmic Microwave Background (CMB). The CMB is foundational to our understanding of modern physics and continues to be a powerful tool driving our understanding of cosmology and particle physics. In this paper, we outline the broad and unique impact of CMB science for the High Energy Cosmic Frontier in the upcoming decade. We also describe the progression of ground-based CMB experiments, which shows that the community is prepared to develop the key capabilities and facilities needed to achieve these transformative CMB measurements., Comment: contribution to Snowmass 2021

This Snowmass 2021 White Paper describes the Cosmic Microwave Background Stage 4 project CMB-S4, which is designed to cross critical thresholds in our understanding of the origin and evolution of the Universe, from the highest energies at the dawn of time through the growth of structure to the present day. We provide an overview of the science case, the technical design, and project plan., Comment: Contribution to Snowmass 2021. arXiv admin note: substantial text overlap with arXiv:1908.01062, arXiv:1907.04473

Cosmological data in the next decade will be characterized by high-precision, multi-wavelength measurements of thousands of square degrees of the same patches of sky. By performing multi-survey analyses that harness the correlated nature of these datasets, we will gain access to new science, and increase the precision and robustness of science being pursued by each individual survey. However, effective application of such analyses requires a qualitatively new level of investment in cross-survey infrastructure, including simulations, associated modeling, coordination of data sharing, and survey strategy. The scientific gains from this new level of investment are multiplicative, as the benefits can be reaped by even present-day instruments, and can be applied to new instruments as they come online., Comment: Contribution to Snowmass 2021

CMB-HD is a proposed millimeter-wave survey over half the sky that would be ultra-deep (0.5 uK-arcmin) and have unprecedented resolution (15 arcseconds at 150 GHz). Such a survey would answer many outstanding questions about the fundamental physics of the Universe. Major advances would be 1.) the use of gravitational lensing of the primordial microwave background to map the distribution of matter on small scales (k~10 h Mpc^(-1)), which probes dark matter particle properties. It will also allow 2.) measurements of the thermal and kinetic Sunyaev-Zel'dovich effects on small scales to map the gas density and velocity, another probe of cosmic structure. In addition, CMB-HD would allow us to cross critical thresholds: 3.) ruling out or detecting any new, light (< 0.1 eV) particles that were in thermal equilibrium with known particles in the early Universe, 4.) testing a wide class of multi-field models that could explain an epoch of inflation in the early Universe, and 5.) ruling out or detecting inflationary magnetic fields. CMB-HD would also provide world-leading constraints on 6.) axion-like particles, 7.) cosmic birefringence, 8.) the sum of the neutrino masses, and 9.) the dark energy equation of state. The CMB-HD survey would be delivered in 7.5 years of observing 20,000 square degrees of sky, using two new 30-meter-class off-axis crossed Dragone telescopes to be located at Cerro Toco in the Atacama Desert. Each telescope would field 800,000 detectors (200,000 pixels), for a total of 1.6 million detectors., Comment: Contribution to Snowmass 2021. Note some text overlap with CMB-HD Astro2020 APC and RFI (arXiv:1906.10134, arXiv:2002.12714). Science case further broadened and updated