Search

The physical properties of the intracluster medium (ICM) reflect signatures of the underlying gravitational potential, mergers and strong interactions with other halos and satellite galaxies, as well as galactic feedback from supernovae and supermassive black holes (SMBHs). Traditionally, clusters have been characterized in terms of summary statistics, such as halo mass, X-ray luminosity, cool-core state, luminosity of AGN, and number of merging components. In this paper of the Extracting Reality from Galaxy Observables with Machine Learning series (ERGO-ML), we instead consider the full information content available in maps of X-ray emission from the ICM. We employ Nearest Neighbour Contrastive Learning (NNCLR) to identify and populate a low-dimensional representation space of such images. Using idealized X-ray maps of the 352 clusters of the TNG-Cluster cosmological magnetohydrodynamical simulation suite, we take three orthogonal projections of each cluster at eight snapshots within the redshift range $0\leq z<1$, resulting in a dataset of $\sim$8,000 images. Our findings reveal that this representation space forms a continuous distribution from relaxed to merging objects, and from centrally-peaked to flat emission profiles. The representation also exhibits clear trends with redshift, with halo, gas, stellar, and SMBH mass, with time since a last major merger, and with indicators of dynamical state. We show that an 8-dimensional representation can be used to predict a variety of cluster properties, find analogs, and identify correlations between physical properties, thereby suggesting causal relationships. Our analysis demonstrates that contrastive learning is a powerful tool for characterizing galaxy clusters from their images alone, allowing us to derive constraints on their physical properties and formation histories using cosmological hydrodynamical galaxy simulations., Comment: Please see more results from TNG-Cluster on astro-ph this week from Rohr+ and Prunier+

ChatGPT and other state-of-the-art large language models (LLMs) are rapidly transforming multiple fields, offering powerful tools for a wide range of applications. These models, commonly trained on vast datasets, exhibit human-like text generation capabilities, making them useful for research tasks such as ideation, literature review, coding, drafting, and outreach. We conducted a study involving 13 astronomers at different career stages and research fields to explore LLM applications across diverse tasks over several months and to evaluate their performance in research-related activities. This work was accompanied by an anonymous survey assessing participants' experiences and attitudes towards LLMs. We provide a detailed analysis of the tasks attempted and the survey answers, along with specific output examples. Our findings highlight both the potential and limitations of LLMs in supporting research while also addressing general and research-specific ethical considerations. We conclude with a series of recommendations, emphasizing the need for researchers to complement LLMs with critical thinking and domain expertise, ensuring these tools serve as aids rather than substitutes for rigorous scientific inquiry., Comment: Paper submitted to RASTI. We share our experience, ethical and legal concerns (5.3), and recommendations for individuals and journals (6.). We welcome feedback

Abell 2744, also known as Pandora's Cluster, is a complex merging galaxy cluster. While a major merger is clear along the north-south axis, the dynamical state of the northwest subcluster has been highly uncertain. We present ultra-deep ($\approx$2.1 Ms) X-ray observations of Abell 2744 obtained with the Chandra X-ray Observatory and reinterpret the multi-wavelength picture with a suite of idealised simulations of galaxy cluster mergers. The new data reveal in unprecedented detail the disruption of cool cores in the three X-ray luminous subclusters and confirm the presence of a shock to the NW. A position-velocity clustering of the cluster member galaxies shows a clearly separated S2 component, with a $\Delta z$ implying a separation of 53 Mpc or a line-of-sight velocity of $4500\ \rm{km \ s^{-1}}$, or likely some combination of the two. While binary simulations allow NW to have undergone a gravitational slingshot after the first pericenter passage, triple merger simulations rule out this scenario, because the two mergers would have had to occur $\sim$0.5 Gyr apart, and the joint impact of the shocks from the two mergers would completely disrupt the SE and NW cool cores; they only reform after 1-2 Gyr, by which point the core separations greatly exceed observations. The scenario that best describes Abell 2744 is a head-on N-S merger $0.5-0.6$ Gyrs ago followed by a first infall of the NW subcluster. Furthermore, we note that a model with three cluster-size halos, with masses consistent with gravitational lensing constraints, nevertheless produces a lensing convergence and surface brightness lower than observed in most of the field of view, whereas the temperatures are consistent with observations. This suggests the presence of a large-scale overdensity, which contributes to the diffuse emission and total surface density without heating the densest gas., Comment: Submitted to the Open Journal of Astrophysics

Although supermassive black holes (SMBHs) reside in the heart of virtually every massive galaxy, it remains debated whether dwarf galaxies also commonly host SMBHs. Because low-mass galaxies may retain a memory of the assembly history of their black holes, probing the black hole occupation fraction of local dwarf galaxies might offer insights into the growth and seeding mechanisms of the first black holes. In this work, we exploit the Western half of the eROSITA all-sky survey (covering $20,000~\rm{deg^2}$) and compile a catalog of accreting SMBHs in local ($D<200$~Mpc) dwarf galaxies. After cleaning our sample from cosmic X-ray background sources, X-ray binaries, and ultraluminous X-ray sources, we identify 74 AGN-dwarf galaxy pairs. Using this large and uniform sample, we derive a luminosity function of dwarf galaxy AGN, fitting it with a power law function and obtaining ${\rm d}N/{\rm d}L_{\rm X} = (15.9\pm2.2)\times L_{\rm X}^{-1.63\pm0.05}$. Measuring the offset between the centroid of dwarf galaxies and the X-ray sources, we find that about $50\%$ of the AGN are likely off-nuclear, in agreement with theoretical predictions. We also compare the black hole-to-stellar mass relation of the AGN in our sample with the local and high-redshift relations, finding that our sources better adhere to the former. This suggests that local AGN across different mass scales underwent a similar growth history. Finally, we compare our sources with semi-analytical models: while our sample is too shallow to distinguish between different seeding models, it favors a growth mechanism linked to the star-formation rate of the host galaxy., Comment: 17 pages, 12 figures, 1 table. Submitted to ApJ. Comments are welcome

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

While supermassive black holes (BHs) are widely observed in the nearby and distant universe, their origin remains debated with two viable formation scenarios with light and heavy seeds. In the light seeding model, the first BHs form from the collapse of massive stars with masses of $10-100 \ \rm{M_{\odot}}$, while the heavy seeding model posits the formation of $10^{4-5} \ \rm{M_{\odot}}$ seeds from direct collapse. The detection of BHs at redshifts $z\gtrsim10$, edging closer to their formation epoch, provides critical observational discrimination between these scenarios. Here, we focus on the JWST-detected galaxy, GHZ 9, at $z\approx10$ that is lensed by the foreground cluster, Abell 2744. Based on 2.1 Ms deep Chandra observations, we detect a candidate X-ray AGN, which is spatially coincident with the high-redshift galaxy, GHZ 9. The BH candidate is inferred to have a bolometric luminosity of $(1.0^{+0.5}_{-0.4})\times10^{46} \ \rm{erg \ s^{-1}}$, which corresponds to a BH mass of $(8.0^{+3.7}_{-3.2})\times10^7 \ \rm{M_{\odot}}$ assuming Eddington-limited accretion. This extreme mass at such an early cosmic epoch suggests the heavy seed origin for this BH candidate. Based on the Chandra and JWST discoveries of extremely high-redshift quasars, we have constructed the first simple AGN luminosity function extending to $z\approx10$. Comparison of this luminosity function with theoretical models indicates an over-abundant $z\approx10$ BH population, consistent with a higher-than-expected seed formation efficiency., Comment: 9 pages, 4 figures, accepted for publication in The Astrophysical Journal Letters

JKCS041 ($z=1.8$) is one of the most distant galaxy cluster systems known, seen when the Universe was less than 4 billion years old. Recent Sunyaev-Zeldovich (SZ) observations show a temperature decrement that is less than expected based on mass estimates of the system from X-ray, weak gravitational lensing and galaxy richness measurements. In this paper we seek to explain the observables - in particular the low SZ decrement and single SZ peak, the projected offset between the X-ray and SZ peaks of $\approx$220 kpc, the gas mass measurements and the lensing mass estimate. We use the GAMER-2 hydrodynamic code to carry out idealized numerical simulations of cluster mergers and compare resulting synthetic maps with the observational data. The observations are not well reproduced by an isolated cluster, while instead they are when considering cluster mergers viewed a few tenths of a Gyr after first core passage. A range of merger scenarios is consistent with the observations, but parts of parameter space can be ruled out, and generically some kind of merger process is necessary to reproduce the offset between the SZ and X-ray peaks. In particular, a total mass of $\approx$2$\times 10^{14} M_\odot$, mass ratio of $\approx$2:3, gas fraction of $0.05-0.1$ and Navarro, Frenk and White (NFW) mass density profile concentration $c$$\approx$5 for both components are scenarios that are consistent with the observational data., Comment: 11 Pages and 11 figures

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

Galaxy cluster mass functions are a function of cosmology, but mass is not a direct observable, and systematic errors abound in all its observable proxies. Mass-free inference can bypass this challenge, but it requires large suites of simulations spanning a range of cosmologies and models for directly observable quantities. In this work, we devise a U-net - an image-to-image machine learning algorithm - to ``paint'' the IllustrisTNG model of baryons onto dark-matter-only simulations of galaxy clusters. Using 761 galaxy clusters with $M_{200c} \gtrsim 10^{14}M_\odot$ from the TNG-300 simulation at $z<1$, we train the algorithm to read in maps of projected dark matter mass and output maps of projected gas density, temperature, and X-ray flux. The models train in under an hour on two GPUs, and then predict baryonic images for $\sim2700$ dark matter maps drawn from the TNG-300 dark-matter-only (DMO) simulation in under two minutes. Despite being trained on individual images, the model reproduces the true scaling relation and scatter for the $M_{DM}-L_X$, as well as the distribution functions of the cluster X-ray luminosity and gas mass. For just one decade in cluster mass, the model reproduces three orders of magnitude in $L_X$. The model is biased slightly high when using dark matter maps from the DMO simulation. The model performs well on inputs from TNG-300-2, whose mass resolution is 8 times coarser; further degrading the resolution biases the predicted luminosity function high. We conclude that U-net-based baryon painting is a promising technique to build large simulated cluster catalogs which can be used to improve cluster cosmology by combining existing full-physics and large $N$-body simulations., Comment: Accepted to MNRAS

Observations of quasars reveal that many supermassive black holes (BHs) were in place less than 700 million years after the Big Bang. However, the origin of the first BHs remains a mystery. Seeds of the first BHs are postulated to be either light (i.e., $10-100~\rm{M_{\odot}})$, remnants of the first stars or heavy (i.e., $10^4-10^5~\rm{M_{\odot}})$, originating from the direct collapse of gas clouds. Harnessing recent data from the Chandra X-ray Observatory, we report the detection of an X-ray-luminous massive BH in a gravitationally-lensed galaxy identified by JWST at $z\approx10.3$ behind the cluster lens Abell 2744. This heavily-obscured quasar with a bolometric luminosity of $L_{\rm bol}\sim5\times10^{45}~\rm{erg\ s^{-1}}$ harbors a $M_{\rm BH}\sim10^7-10^8~\rm{M_{\odot}}$ BH assuming accretion at the Eddington limit. This mass is comparable to the inferred stellar mass of its host galaxy, in contrast to what is found in the local Universe wherein the BH mass is $\sim0.1\%$ of the host galaxy's stellar mass. The combination of such a high BH mass and large BH-to-galaxy stellar mass ratio just $\sim$500 Myrs after the Big Bang was theoretically predicted and is consistent with a picture wherein BHs originated from heavy seeds., Comment: 27 pages, 6 figures, accepted

We conduct a systematic tidal disruption event (TDE) demographics analysis using the largest sample of optically selected TDEs. A flux-limited, spectroscopically complete sample of 33 TDEs is constructed using the Zwicky Transient Facility over three years (from October 2018 to September 2021). We infer the black hole (BH) mass ($M_{\rm BH}$) with host galaxy scaling relations, showing that the sample $M_{\rm BH}$ ranges from $10^{5.1}\,M_\odot$ to $10^{8.2}\,M_\odot$. We developed a survey efficiency corrected maximum volume method to infer the rates. The rest-frame $g$-band luminosity function (LF) can be well described by a broken power-law of $\phi (L_g)\propto [(L_g / L_{\rm bk})^{0.3} + (L_g / L_{\rm bk})^{2.6}]^{-1}$, with $L_{\rm bk}=10^{43.1}\,{\rm erg\,s^{-1}}$. In the BH mass regime of $10^{5.3}\lesssim (M_{\rm BH}/M_\odot) \lesssim 10^{7.3}$, the TDE mass function follows $\phi(M_{\rm BH})\propto M_{\rm BH}^{-0.25}$, which favors a flat local BH mass function ($dn_{\rm BH}/d{\rm log}M_{\rm BH}\approx{\rm constant}$). We confirm the significant rate suppression at the high-mass end ($M_{\rm BH}\gtrsim 10^{7.5}\,M_\odot$), which is consistent with theoretical predictions considering direct capture of hydrogen-burning stars by the event horizon. At a host galaxy mass of $M_{\rm gal}\sim 10^{10}\,M_\odot$, the average optical TDE rate is $\approx 3.2\times 10^{-5}\,{\rm galaxy^{-1}\,yr^{-1}}$. We constrain the optical TDE rate to be [3.7, 7.4, and 1.6$]\times 10^{-5}\,{\rm galaxy^{-1}\,yr^{-1}}$ in galaxies with red, green, and blue colors., Comment: Replaced following peer-review process. 38 pages, 23 figures. Accepted for publication in ApJL

Dwarf galaxies are promising test beds for constraining models of supermassive and intermediate-mass black holes (MBH) via their black hole occupation fraction (BHOF). Disentangling seeding from the confounding effects of mass assembly over a Hubble time is a challenging problem, that we tackle in this study with a suite of semi-analytical models (SAMs). We show how measured BHOF depends on the lowest black hole mass or AGN luminosity achieved by a survey. To tell seeding models apart, we need to detect or model all AGN brighter than $10^{37}\ \rm{erg \ s^{-1}}$ in galaxies of $M_* \sim 10^{8-10} \ \rm{M_{\odot}}$. Shallower surveys, like eRASS, cannot distinguish between seed models even with the compensation of a much larger survey volume. We show that the AMUSE survey, with its inference of the MBH population underlying the observed AGN, strongly favors heavy seed models, growing with either a power-law Eddington Ratio Distribution Function (ERDF) or one in which black hole accretion is tagged to the star-formation rate (AGN-MS). These two growth channels can then be distinguished by the AGN luminosity function at $> 10^{40}\ \rm{erg \ s^{-1}}$, with the AGN-MS model requiring more accretion than observed at z $\sim$ 0. Thus, current X-ray observations favour heavy seeds whose Eddington ratios follow a power-law distribution. The different models also predict different radio scaling relations, which we quantify using the fundamental plane of black hole activity. We close with recommendations for the design of upcoming multi-wavelength campaigns that can optimally detect MBHs in dwarf galaxies., Comment: Accepted to ApJ

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

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

The hot ($>10^6$ K) phase of the circumgalactic medium (CGM) contains a large fraction of baryons in galaxies. It also retains signatures of the processes that shaped the galaxies, such as feedback from active galactic nuclei (AGNs) and supernovae, and offers a uniquely powerful way to constrain theoretical models of feedback. It is, however, notoriously difficult to detect. By stacking 2643 optically selected galaxies in the eROSITA Final Equatorial Depth Survey (eFEDS), we present spatially resolved properties of the extended CGM in both star-forming and quiescent galaxies spanning an order of magnitude in stellar mass. We mask out resolved point sources and galaxy groups/clusters and model the contribution from X-ray binaries and the hot ISM, producing accurate radial profiles. We compare the profiles to mock X-ray observations of galaxy stacks in the IllustrisTNG100 (TNG) and EAGLE cosmological simulations. We detect extended emission from both the high-mass ($10.7<\log(M_*/M_\odot)<11.2$) and low-mass ($10.2<\log(M_*/M_\odot)<10.7$) galaxy stacks. Galaxies have somewhat more luminous CGM between $10-100$~kpc if they are more massive or star-forming. However, the luminosity increases slower with stellar mass than predicted in simulations. Simulated quenched galaxies are far dimmer than observed, suggesting that they rely too heavily on CGM ejection for quenching. Star-forming galaxies are observed to have flatter and more extended profiles than in simulations, suggesting under-efficient stellar feedback models. Our results highlight the need to modify future prescriptions of galaxy feedback models., Comment: Accepted to ApJ Letters

Kelvin-Helmholtz Instabilities (KHI) along contact discontinuities in galaxy clusters have been used to constrain the strength of magnetic fields in galaxy clusters, following the assumption that, as magnetic field lines drape around the interface between the cold and hot phases, their magnetic tension resists the growth of perturbations. This has been observed in simulations of rigid objects moving through magnetised media and sloshing galaxy clusters, and then applied in interpreting observations of merger cold fronts. Using a suite of MHD simulations of binary cluster mergers, we show that even magnetic field strengths stronger than yet observed ($\beta = P_{\rm th}/P_B = 50$) show visible KHI features. This is because our initial magnetic field is tangled, producing Alfven waves and associated velocity fluctuations in the ICM; stronger initial fields therefore seed larger fluctuations, so that even a reduced growth rate due to magnetic tension produces a significant KHI. The net result is that a stronger initial magnetic field produces more dramatic fluctuations in surface brightness and temperature, not the other way around. We show that this is hard to distinguish from the evolution of turbulent perturbations of the same initial magnitude. Therefore, in order to use observations of KHI in the ICM to infer magnetic field strengths by comparing to idealized simulations, the perturbations which seed the KHI must be well-understood and (if possible) carefully controlled.

Galaxy cluster mergers are a powerful laboratory for testing cosmological and astrophysical models. However, interpreting individual merging clusters depends crucially on their merger configuration, defined by the masses, velocities, impact parameters, and orientation of the merger axis with respect to the plane of the sky. In this work, we investigate the impact of merger parameters on the X-ray emitting intracluster medium and gravitational lensing maps using a suite of idealised simulations of binary cluster mergers performed using the GAMER-2 code. As a test case, we focus on modeling the Bullet Cluster-like merging system Abell 2146, in which deep \textit{Chandra} X-ray and lensing observations revealed prominent merger shocks as well as the mass distribution and substructures associated with this merging cluster. We identify the most interesting parameter combinations, and evaluate the effects of various parameters on the properties of merger shocks observed by deep \textit{Chandra} and lensing observations. We show that due gravitational compression of the cluster halos during the merger, previous mass estimates from weak lensing are too high. The plane of the merger is tilted further from the plane of the sky than estimated previously, up to $30^\circ$ from the plane of the sky. We discuss the applicability of our results to multi-wavelength observations of merging galaxy clusters and their use as probes of cosmology and plasma physics.

The intracluster medium (ICM) is a multi-phase environment, dynamically regulated by Active Galactic Nuclei (AGN), the motions of galaxies through it, and mergers with other clusters. AGN as a central heating source are key to preventing runaway cooling flows, but their role in heating cores in a cosmological context is still poorly understood. The activity of the AGN is strongly linked to star formation, especially in the Brightest Cluster Galaxy (BCG), likely because both rely on cold phase gas. A self-consistent model for AGN and star formation in galaxy clusters thus requires cosmological context, higher resolution, and a careful modeling of cooling and heating balance. In this paper, we use the high-resolution hydrodynamical cosmological simulation of the RomulusC galaxy cluster to study in detail the role of AGN and a major, head-on merger in shaping the cluster core. The unprecedented resolution of the RomulusC simulation captures the multiphase structure of the ICM. The realistic large-scale outflows launched by very small-scale thermal injections, the improved modeling of turbulent diffusion and mixing, and the particle nature of the simulation allow us to carefully separate different heating channels. We show that AGN activity, while efficient at regulating star formation, is incapable of destroying a CC. Instead, that process is facilitated by a head-on, 1:8 mass ratio merger. The merger generates bulk and turbulent motions, which in turn mix high entropy gas generated by AGN and merger driven shocks, turbulent dissipation and sloshing of the ICM by infalling substructures. While central cooling times remain shorter than the Hubble time, restoring a CC is made more difficult by the reduced precipitation rates at larger radii, emphasizing that the AGN-ICM connection is truly a multi-scale problem., Comment: 16 pages, 9 figures

Clusters of galaxies, the largest collapsed structures in the Universe, are located at the intersection of extended filaments of baryons and dark matter. Cosmological accretion onto clusters through large scale filaments adds material at cluster outskirts. Kinetic energy in the form of bulk motions and turbulence due to this accretion provides a form of pressure support against gravity, supplemental to thermal pressure. Significant amount of non-thermal pressure support could bias cluster masses derived assuming hydrostatic equilibrium, the primary proxy for cluster cosmology studies. Sensitive measurements of Doppler broadening and shift of astrophysical lines, and the relative fluctuations in thermodynamical quantities (e.g., density, pressure, and entropy) are primary diagnostic tools. Forthcoming planned and proposed X-ray (with large etendue, throughput, and high spectral resolution) and SZ observatories will provide crucial information on the assembly and virialisation processes of clusters, involving turbulent eddies cascading at various spatial scales and larger gas bulk motions in their external regions to the depth or their potential wells., Comment: 5 pages, 3 figures, Science white paper submitted to the Astro2020 Decadal Survey

We present the first results from RomulusC, the highest resolution cosmological hydrodynamic simulation of a galaxy cluster run to date. RomulusC, a zoom-in simulation of a halo with $z=0$ mass $10^{14}$ M$_{\odot}$, is run with the same sub-grid physics and resolution as Romulus25 (Tremmel et al. 2017). With unprecedented mass and spatial resolution, RomulusC represents a unique opportunity to study the evolution of galaxies in dense environments down to dwarf masses. We demonstrate that RomulusC results in an intracluster medium (ICM) consistent with observations. The star formation history and stellar mass of the brightest cluster galaxy (BCG) is consistent with observations and abundance matching results, indicating that our sub-grid models, optimized only to reproduce observations of field dwarf and Milky Way mass galaxies, are able to produce reasonable galaxy masses and star formation histories in much higher mass systems. Feedback from supermassive black holes (SMBHs) regulates star formation by driving large-scale, collimated outflows that coexist with a low entropy core. We find that non-BCG cluster member galaxies are substantially quenched compared to the field down to dwarf galaxy masses and, at low masses, quenching is seen to have no dependence on mass or distance from the cluster center. This enhanced quenched population extends beyond $R_{200}$ and is in place at high redshift. Similarly, we predict that SMBH activity is significantly suppressed within clusters outside of the BCG, but show how the effect could be lost when only focusing on the brightest AGN in the most massive galaxies., Comment: Submitted to MNRAS, 29 pages, 19 figures, Accepted to MNRAS

We map the lensing-inferred substructure in the first three clusters observed by the Hubble Space Telescope Frontier Fields Initiative (HSTFF): Abell 2744 (z = 0.308), MACSJ0416, (z = 0.396) and MACSJ1149 (z = 0.543). Statistically resolving dark-matter subhaloes down to ~10^{9.5} solar masses, we compare the derived subhalo mass functions (SHMFs) to theoretical predictions from analytical models and with numerical simulations in a Lambda Cold Dark Matter (LCDM) cosmology. Mimicking our observational cluster member selection criteria in the HSTFF, we report excellent agreement in both amplitude and shape of the SHMF over four decades in subhalo mass (10^{9-13} solar masses). Projection effects do not appear to introduce significant errors in the determination of SHMFs from simulations. We do not find evidence for a substructure crisis, analogous to the missing satellite problem in the Local Group, on cluster scales, but rather excellent agreement of the count-matched HSTFF SHMF down to M_{sub halo}/M_{halo} ~ 10^{-5}. However, we do find discrepancies in the radial distribution of sub haloes inferred from HSTFF cluster lenses compared to determinations from simulated clusters. This suggests that although the selected simulated clusters match the HSTFF sample in mass, they do not adequately capture the dynamical properties and complex merging morphologies of these observed cluster lenses. Therefore, HSTFF clusters are likely observed in a transient evolutionary stage that is presently insufficiently sampled in cosmological simulations. The abundance and mass function of dark matter substructure in cluster lenses continues to offer an important test of the LCDM paradigm, and at present we find no tension between model predictions and observations., Comment: 20 pages, 29 figures, Accepted for publication, in press, MNRAS, replaced with added author affiliations and added reference