Showing 1–50 of 208 results for author: Pillepich, A

The PARADIGM project I: How early merger histories shape the present-day sizes of Milky-Way-mass galaxies

Authors: Gandhali D. Joshi, Andrew Pontzen, Oscar Agertz, Martin P. Rey, Justin Read, Annalisa Pillepich

Abstract: The way in which mergers affect galaxy formation depends on both feedback processes, and on the geometry and strength of the mergers themselves. We introduce the PARADIGM project, where we study the response of a simulated Milky-Way-mass galaxy forming in a cosmological setting to differing merger histories, using genetically modified initial conditions. Each initial condition is simulated with th… ▽ More The way in which mergers affect galaxy formation depends on both feedback processes, and on the geometry and strength of the mergers themselves. We introduce the PARADIGM project, where we study the response of a simulated Milky-Way-mass galaxy forming in a cosmological setting to differing merger histories, using genetically modified initial conditions. Each initial condition is simulated with the VINTERGATAN and IllustrisTNG codes. While VINTERGATAN has been developed with an emphasis on resolving the cold interstellar medium, IllustrisTNG uses a subgrid two-phase model and consequently scales to large volume simulations, making them ideal to examine complementary views on how merger histories and feedback interact. Our genetic modifications alter the mass ratio of an important $z \approx 2$ merger while maintaining the halo's $z=0$ mass. Whether simulated with VINTERGATAN or IllustrisTNG, smaller mass ratios for this early merger result in larger galaxies at $z=0$, due to a greater build up of a kinematically cold disc. We conclude that such broad trends are robustly reproducible; however, the normalization of the resulting stellar sizes is substantially different in the two codes (ranging between $0.5-1.7\ \rm{kpc}$ for VINTERGATAN but $1.3-7.0\ \rm{kpc}$ for IllustrisTNG). The VINTERGATAN galaxies systematically form stars earlier, leading to a larger bulge component. Despite the difference in size normalization, both simulation suites lie on the observed size-mass relation for their respective morphological types. In light of these results, we discuss the interplay between internal processes and large scale gravitational interactions and gas accretion, and how the two galaxy models converge on similar emergent trends but along different evolutionary pathways. △ Less

Submitted 28 June, 2024; originally announced July 2024.

Comments: Submitted to MNRAS. 22 pages. Comments are welcome

Constraints on the in-situ and ex-situ stellar masses in nearby galaxies with Artificial Intelligence

Authors: Eirini Angeloudi, Jesús Falcón-Barroso, Marc Huertas-Company, Alina Boecker, Regina Sarmiento, Lukas Eisert, Annalisa Pillepich

Abstract: The hierarchical model of galaxy evolution suggests that the impact of mergers is substantial on the intricate processes that drive stellar assembly within a galaxy. However, accurately measuring the contribution of accretion to a galaxy's total stellar mass and its balance with in-situ star formation poses a persistent challenge, as it is neither directly observable nor easily inferred from obser… ▽ More The hierarchical model of galaxy evolution suggests that the impact of mergers is substantial on the intricate processes that drive stellar assembly within a galaxy. However, accurately measuring the contribution of accretion to a galaxy's total stellar mass and its balance with in-situ star formation poses a persistent challenge, as it is neither directly observable nor easily inferred from observational properties. Here, we present theory-motivated predictions for the fraction of stellar mass originating from mergers in a statistically significant sample of nearby galaxies, using data from MaNGA. Employing a robust machine learning model trained on mock MaNGA analogs (MaNGIA) in turn obtained from a cosmological simulation (TNG50), we unveil that in-situ stellar mass dominates almost across the entire stellar mass spectrum (1e9Msun < M* < 1e12Msun). Only in more massive galaxies (M* > 1e11Msun) does accreted mass become a substantial contributor, reaching up to 35-40% of the total stellar mass. Notably, the ex-situ stellar mass in the nearby universe exhibits significant dependence on galaxy characteristics, with higher accreted fractions favored by elliptical, quenched galaxies and slow rotators, as well as galaxies at the center of more massive dark matter halos. △ Less

Submitted 28 June, 2024; originally announced July 2024.

Comments: 16 pages, 7 Figures. Accepted for publication

The Extremely Metal Rich Knot of Stars at the Heart of the Galaxy

Authors: Hans-Walter Rix, Vedant Chandra, Gail Zasowski, Annalisa Pillepich, Sergey Khoperskov, Sofia Feltzing, Rosemary F. Wyse, Neige Frankel, Danny Horta, Juna Kollmeier, Keivan G. Stassun, Melissa Ness, Jonathan C. Bird, David L. Nidever, Jose G. Fernandez, João A. Amarante, Chervin F. Laporte, Jianhui Lian

Abstract: We show with Gaia XP spectroscopy that extremely metal-rich stars in the Milky Way (EMR; $[M/H]_{XP} > 0.5$) - but only those - are largely confined to a tight "knot" at the center of the Galaxy. This EMR knot is round in projection, has a fairly abrupt edge near $\sim 1.5$kpc, and is a dynamically hot system. This central knot also contains very metal-rich (VMR; $+0.2\le [M/H]_{XP} \le +0.4$) sta… ▽ More We show with Gaia XP spectroscopy that extremely metal-rich stars in the Milky Way (EMR; $[M/H]_{XP} > 0.5$) - but only those - are largely confined to a tight "knot" at the center of the Galaxy. This EMR knot is round in projection, has a fairly abrupt edge near $\sim 1.5$kpc, and is a dynamically hot system. This central knot also contains very metal-rich (VMR; $+0.2\le [M/H]_{XP} \le +0.4$) stars. However, in contrast to EMR stars, the bulk of VMR stars form an extended, highly flattened distribution in the inner Galaxy ($R_{\mathrm{GC}}\lesssim 5$ kpc). We draw on TNG50 simulations of Milky Way analogs for context and find that compact, metal-rich knots confined to $<1.5$kpc are a universal feature. In typical simulated analogs, the top 5-10% most metal-rich stars are confined to a central knot; however, in our Milky Way data this fraction is only 0.1%. Dust-penetrating wide-area near-infrared spectroscopy, such as SDSS-V, will be needed for a rigorous estimate of the fraction of stars in the Galactic EMR knot. Why in our Milky Way only EMR giants are confined to such a central knot remains to be explained. Remarkably, the central few kiloparsecs of the Milky Way harbor both the highest concentration of metal-poor stars (the `poor old heart') and almost all EMR stars. This highlights the stellar population diversity at the bottom of galactic potential wells. △ Less

Submitted 3 June, 2024; originally announced June 2024.

Comments: 11 pages, 7 figures, submitted to ApJ

The many colors of the TNG100 simulation

Authors: Andrea Gebek, Ana Trčka, Maarten Baes, Marco Martorano, Annalisa Pillepich, Anand Utsav Kapoor, Angelos Nersesian, Arjen van der Wel

Abstract: We apply the 3D dust radiative transfer code SKIRT to the low-redshift ($z\leq0.1$) galaxy population in the TNG100 cosmological simulation, the fiducial run of the IllustrisTNG project. We compute global fluxes and spectral energy distributions (SEDs) from the far-ultraviolet to the sub-millimeter for $\approx\,$60 000 galaxies. Our post-processing methodology follows the study of Trčka et al. (2… ▽ More We apply the 3D dust radiative transfer code SKIRT to the low-redshift ($z\leq0.1$) galaxy population in the TNG100 cosmological simulation, the fiducial run of the IllustrisTNG project. We compute global fluxes and spectral energy distributions (SEDs) from the far-ultraviolet to the sub-millimeter for $\approx\,$60 000 galaxies. Our post-processing methodology follows the study of Trčka et al. (2022) of the higher-resolution TNG50 simulation. We verify that TNG100 reproduces observational luminosity functions at low redshifts to excellent precision, unlike TNG50. Additionally, we test the realism of our TNG100 plus SKIRT fluxes by comparing various flux and color relations to data from the GAMA survey. TNG100 broadly reproduces the observed distributions, but we predict ultraviolet colors that are too blue by $\approx\,$0.4 mag, possibly related to the extinction in the star-forming regions subgrid model not being selective enough. Furthermore, we find that the simulated galaxies exhibit mid-infrared fluxes elevated by up to $\approx\,$0.5 mag that we attribute to overly effective stochastic heating of the diffuse dust. All synthetic broadband fluxes and SEDs are made publicly available in three orientations and four apertures, and can readily be used to study TNG100 galaxies in a mock observational fashion. △ Less

Submitted 30 May, 2024; v1 submitted 8 May, 2024; originally announced May 2024.

Comments: Main text 14 pages, 5 figures, accepted to MNRAS. Generated data products are publicly available at www.tng-project.org/gebek24

MAGAZ3NE: Massive, Extremely Dusty Galaxies at $z\sim2$ Lead to Photometric Overestimation of Number Densities of the Most Massive Galaxies at $3<z<4$

Authors: Ben Forrest, M. C. Cooper, Adam Muzzin, Gillian Wilson, Danilo Marchesini, Ian McConachie, Percy Gomez, Marianna Annunziatella, Z. Cemile Marsan, Joey Braspenning, Wenjun Chang, Gabriella de Lucia, Fabio Fontanot, Michaela Hirschmann, Dylan Nelson, Annalisa Pillepich, Joop Schaye, Stephanie M. Urbano Stawinski, Mauro Stefanon, Lizhi Xie

Abstract: We present rest-frame optical spectra from Keck/MOSFIRE and Keck/NIRES of 16 candidate ultramassive galaxies targeted as part of the Massive Ancient Galaxies at $z>3$ Near-Infrared (MAGAZ3NE) Survey. These candidates were selected to have photometric redshifts $3\lesssim z_{\rm phot}<4$, photometric stellar masses log($M$/M$_\odot$)$>11.7$, and well-sampled photometric spectral energy distribution… ▽ More We present rest-frame optical spectra from Keck/MOSFIRE and Keck/NIRES of 16 candidate ultramassive galaxies targeted as part of the Massive Ancient Galaxies at $z>3$ Near-Infrared (MAGAZ3NE) Survey. These candidates were selected to have photometric redshifts $3\lesssim z_{\rm phot}<4$, photometric stellar masses log($M$/M$_\odot$)$>11.7$, and well-sampled photometric spectral energy distributions (SEDs) from the UltraVISTA and VIDEO surveys. In contrast to previous spectroscopic observations of blue star-forming and post-starburst ultramassive galaxies, candidates in this sample have very red SEDs implying significant dust attenuation, old stellar ages, and/or active galactic nuclei (AGN). Of these galaxies, eight are revealed to be heavily dust-obscured $2.0<z<2.7$ galaxies with strong emission lines, some showing broad features indicative of AGN, three are Type I AGN hosts at $z>3$, one is a $z\sim1.2$ dusty galaxy, and four galaxies do not have a confirmed spectroscopic redshift. In fact, none of the sample has |$z_{\rm spec}-z_{\rm phot}$|$<0.5$, suggesting difficulties for photometric redshift programs in fitting similarly red SEDs. The prevalence of these red interloper galaxies suggests that the number densities of high-mass galaxies are overestimated at $z\gtrsim3$ in large photometric surveys, helping to resolve the `impossibly early galaxy problem' and leading to much better agreement with cosmological galaxy simulations. A more complete spectroscopic survey of ultramassive galaxies is required to pin down the uncertainties on their number densities in the early universe. △ Less

Submitted 29 April, 2024; originally announced April 2024.

Comments: Submitted to ApJ. Comments welcome!

The effects of environment on galaxies' dynamical structures: From simulations to observations

Authors: Yuchen Ding, Ling Zhu, Annalisa Pillepich, Glenn van de Ven, Enrichetta Iodice, Enrico Maria Corsini, Francesca Pinna

Abstract: We studied the effects of cluster environments on galactic structures by using the TNG50 cosmological simulation and observed galaxies in the Fornax cluster. We focused on galaxies with stellar masses of $10^{8-12}M_{\odot}$ at z=0 that reside in Fornax-like clusters with total masses of $M_{200c} = 10^{13.4-14.3}M_{\odot}$. We characterized the stellar structures by decomposing each galaxy into a… ▽ More We studied the effects of cluster environments on galactic structures by using the TNG50 cosmological simulation and observed galaxies in the Fornax cluster. We focused on galaxies with stellar masses of $10^{8-12}M_{\odot}$ at z=0 that reside in Fornax-like clusters with total masses of $M_{200c} = 10^{13.4-14.3}M_{\odot}$. We characterized the stellar structures by decomposing each galaxy into a dynamically cold disk and a hot non-disk component, and studied the evolution of both the stellar and gaseous constituents. In TNG50, we find that the cold gas is quickly removed when a galaxy falls into a Fornax-mass cluster. About 87\% of the galaxies have lost $80\%$ of their star-forming gas at 4 billion years after infall, with the remaining gas concentrating in the inner regions of the galaxy. The radius of the star-forming gaseous disk decreases to half its original size at 4 billion years after infall for 66\% of the galaxies. As a result, star formation in the extended dynamically cold disk sharply decreases, even though a low level of SF persists at the center for a few additional gigayears. This leads to a tight correlation between the average stellar age in the dynamically cold disk and the infall time of galaxies. Furthermore, the luminosity fraction of the dynamically cold disk in ancient infallers is only about 1/3 of that in recent infallers, controlling for galaxy stellar mass. This quantitatively agrees with what is observed in early-type galaxies in the Fornax cluster. Gas removal stops the possible growth of the disk, with gas removed earlier in galaxies that fell in earlier, and hence the cold-disk fraction is correlated with the infall time. The stellar disk can be significantly disrupted by tidal forces after infall, through a long-term process that enhances the difference among cluster galaxies with different infall times. △ Less

Submitted 1 April, 2024; originally announced April 2024.

Comments: Accepted to A&A, 16 pages, 12 figures

Journal ref: A&A 686, A184 (2024)

An early dark matter-dominated phase in the assembly history of Milky Way-mass galaxies suggested by the TNG50 simulation and JWST observations

Authors: Anna de Graaff, Annalisa Pillepich, Hans-Walter Rix

Abstract: Whereas well-studied galaxies at cosmic noon are found to be baryon-dominated within the effective radius, recent JWST observations of $z\sim6-7$ galaxies with stellar masses of only $M_*\sim10^{8-9}\,{\rm M_\odot}$ surprisingly indicate that they are dark matter-dominated within $r_{\rm e}\approx 1\,$kpc. Here, we place these high-redshift measurements in the context of the TNG50 galaxy formation… ▽ More Whereas well-studied galaxies at cosmic noon are found to be baryon-dominated within the effective radius, recent JWST observations of $z\sim6-7$ galaxies with stellar masses of only $M_*\sim10^{8-9}\,{\rm M_\odot}$ surprisingly indicate that they are dark matter-dominated within $r_{\rm e}\approx 1\,$kpc. Here, we place these high-redshift measurements in the context of the TNG50 galaxy formation simulation, by measuring the central (within $1\,$kpc) stellar, gas, and dark matter masses of galaxies in the simulation. The central baryon fraction varies strongly with galaxy stellar mass in TNG50, and this $M_*$-dependence is remarkably constant across $0<z<6$: galaxies of low stellar mass ($M_*\sim10^{8-9}\,{\rm M_\odot}$) are dark matter-dominated, as $f_{\rm baryon}(<1\,{\rm kpc})\sim0.25$. At $z=6$, the baryonic mass in the centers of low-mass galaxies is largely comprised of gas, exceeding the stellar mass component by a factor $\sim4$. We use the simulation to track the typical evolution of such low-mass dark matter-dominated galaxies at $z=6$, and show that these systems become baryon-dominated in their centers at cosmic noon, with high stellar-to-gas mass ratios, and grow to galaxies of $M_*\sim10^{10.5}\,{\rm M_\odot}$ at $z=0$. Comparing to the dynamical and stellar mass measurements from observations at high redshifts, these findings suggest that the inferred star formation efficiency in the early Universe is broadly in line with the established assumptions for the cosmological simulations. Moreover, our results imply that the JWST observations may indeed have reached the early, low-mass regime where the central parts of galaxies transition from being dark matter-dominated to being baryon-dominated. △ Less

Submitted 16 May, 2024; v1 submitted 1 March, 2024; originally announced March 2024.

Comments: 10 pages, 4 figures; accepted for publication in ApJL

The Hot Circum-Galactic Medium in the eROSITA All Sky Survey II. Scaling Relations between X-ray Luminosity and Galaxies' Mass

Authors: Yi Zhang, Johan Comparat, Gabriele Ponti, Andrea Merloni, Kirpal Nandra, Frank Haberl, Nhut Truong, Annalisa Pillepich, Nicola Locatelli, Xiaoyuan Zhang, Jeremy Sanders, Xueying Zheng, Ang Liu, Paola Popesso, Teng Liu, Peter Predehl, Mara Salvato

Abstract: Understanding how the hot circum-galactic medium (CGM) properties relate to the galaxy's properties can constrain galaxy evolution models. We aim to measure the scaling relations between the X-ray luminosity of the hot CGM and the fundamental properties of a galaxy, i.e., its stellar mass and halo mass. We calculate the X-ray luminosity of the hot CGM based on the surface brightness profiles of ce… ▽ More Understanding how the hot circum-galactic medium (CGM) properties relate to the galaxy's properties can constrain galaxy evolution models. We aim to measure the scaling relations between the X-ray luminosity of the hot CGM and the fundamental properties of a galaxy, i.e., its stellar mass and halo mass. We calculate the X-ray luminosity of the hot CGM based on the surface brightness profiles of central galaxy samples measured in Zhang et al. (2024a, submitted) from Spectrum Roentgen Gamma (SRG)/eROSITA all-sky survey data. We relate the X-ray luminosity to the galaxies' stellar and halo mass. We compare the observed relations to the TNG, EAGLE, and SIMBA simulations. The hot CGM X-ray luminosity correlates with the galaxy's stellar mass ($M_*$). It increases from $2.1 \pm 1.3\times10^{39} \rm erg/s$ to $2.0 \pm 0.1\times10^{41} \rm erg/s$, when $\log(M_*)$ increases from 10.0 to 11.5. A double power law describes the correlation, with a break at $\log(M_*)=11.28\pm 0.03$ and a power-law index or $1.9\pm 0.2$ ($4.2\pm0.1$) below (above) the break. The hot CGM X-ray luminosity as a function of halo mass is measured within $\log(M_{\rm 500c})=11.3-13.7$, extending our knowledge of the scaling relation by more than two orders of magnitude. $L_{\rm X,CGM}$ increases with $M_{\rm 500c}$ from $2.7 \pm 0.9\times10^{39}\ \rm erg/s$ at $\log(M_{\rm 500c})=11.3$ to $9.2 \pm 0.4\times10^{41}\ \rm erg/s$ at $\log(M_{\rm 500c})=13.7$. The relation follows a power law of $\log(L_{\rm X,CGM})= (1.35\pm 0.04)\times \log(M_{\rm 500c})+(23.8\pm0.5)$. We find a general agreement between simulations and observation. We explore, at the low mass end, the average scaling relations between the CGM X-ray luminosity and the galaxy's stellar mass or halo mass, which constitutes a new benchmark for galaxy evolution models and feedback processes. △ Less

Submitted 29 April, 2024; v1 submitted 30 January, 2024; originally announced January 2024.

Comments: 8 pages, 4 figures, submitted to A&A

The Hot Circum-Galactic Medium in the eROSITA All Sky Survey I. X-ray Surface Brightness Profiles

Authors: Yi Zhang, Johan Comparat, Gabriele Ponti, Andrea Merloni, Kirpal Nandra, Frank Haberl, Nicola Locatelli, Xiaoyuan Zhang, Jeremy Sanders, Xueying Zheng, Ang Liu, Paola Popesso, Teng Liu, Nhut Truong, Annalisa Pillepich, Peter Predehl, Mara Salvato

Abstract: We apply a stacking technique to gain enough statistics to detect the hot CGM around galaxies. We use the X-ray data from the first four SRG/eROSITA all-sky surveys (eRASS:4). We discuss how the satellite galaxies could bias the stacking and carefully build the central galaxy samples based on the SDSS spectroscopic survey and halo-based group finder algorithm and isolated galaxy sample from the ni… ▽ More We apply a stacking technique to gain enough statistics to detect the hot CGM around galaxies. We use the X-ray data from the first four SRG/eROSITA all-sky surveys (eRASS:4). We discuss how the satellite galaxies could bias the stacking and carefully build the central galaxy samples based on the SDSS spectroscopic survey and halo-based group finder algorithm and isolated galaxy sample from the ninth data release of the DESI Legacy survey (LS DR9, photometric). We mask the detected X-ray point sources and carefully model the X-ray emission from the unresolved active galactic nuclei (AGN) and X-ray binaries (XRB) to obtain the X-ray emission from the hot CGM. The X-ray surface brightness profiles are measured for $\log(M_*)>10.0$ or $\log(M_{\rm 200m})>11.5$ central galaxies, and $\log(M_*)>9.5$ isolated galaxies. We detect the X-ray emission around MW-mass ($\log(M_*)=10.5-11.0$) and more massive central galaxies extending up to the virial radius ($R_{\rm vir}$). The signal-to-noise ratio of the extended emission around the MW-mass galaxy is about $8.1σ$ within $R_{\rm vir}$. We use the $β$ model to describe the X-ray surface brightness profiles of the hot CGM. We obtain a central surface brightness of $S_{\rm X,0}\approx 3.1\times10^{35}\rm erg/s/kpc^2$ and $β\approx0.42$ for MW-mass galaxy. We estimate the baryon budget of the hot CGM and obtain a value lower than the $ΛCDM$ cosmology predicted. Our results set a firm footing for the presence of a hot CGM around MW-mass and more massive central galaxies. These measurements constitute a new benchmark for galaxy evolution models and possible implementations of feedback processes therein. △ Less

Submitted 29 April, 2024; v1 submitted 30 January, 2024; originally announced January 2024.

Comments: 21 pages, 12 figures, submitted to A&A

Preparing for low surface brightness science with the Vera C. Rubin Observatory: A Comparison of Observable and Simulated Intracluster Light Fractions

Authors: Sarah Brough, Syeda Lammim Ahad, Yannick M. Bahe, Amaël Ellien, Anthony H. Gonzalez, Yolanda Jiménez-Teja, Lucas C. Kimmig, Garreth Martin, Cristina Martínez-Lombilla, Mireia Montes, Annalisa Pillepich, Rossella Ragusa, Rhea-Silvia Remus, Chris A. Collins, Johan H. Knapen, J. Chris Mihos

Abstract: Intracluster Light (ICL) provides an important record of the interactions galaxy clusters have undergone. However, we are limited in our understanding by our measurement methods. To address this we measure the fraction of cluster light that is held in the Brightest Cluster Galaxy and ICL (BCG+ICL fraction) and the ICL alone (ICL fraction) using observational methods (Surface Brightness Threshold-S… ▽ More Intracluster Light (ICL) provides an important record of the interactions galaxy clusters have undergone. However, we are limited in our understanding by our measurement methods. To address this we measure the fraction of cluster light that is held in the Brightest Cluster Galaxy and ICL (BCG+ICL fraction) and the ICL alone (ICL fraction) using observational methods (Surface Brightness Threshold-SB, Non-Parametric Measure-NP, Composite Models-CM, Multi-Galaxy Fitting-MGF) and new approaches under development (Wavelet Decomposition-WD) applied to mock images of 61 galaxy clusters (14<log10 M_200c/M_solar <14.5) from four cosmological hydrodynamical simulations. We compare the BCG+ICL and ICL fractions from observational measures with those using simulated measures (aperture and kinematic separations). The ICL fractions measured by kinematic separation are significantly larger than observed fractions. We find the measurements are related and provide equations to estimate kinematic ICL fractions from observed fractions. The different observational techniques give consistent BCG+ICL and ICL fractions but are biased to underestimating the BCG+ICL and ICL fractions when compared with aperture simulation measures. Comparing the different methods and algorithms we find that the MGF algorithm is most consistent with the simulations, and CM and SB methods show the smallest projection effects for the BCG+ICL and ICL fractions respectively. The Ahad (CM), MGF and WD algorithms are best set up to process larger samples, however, the WD algorithm in its current form is susceptible to projection effects. We recommend that new algorithms using these methods are explored to analyse the massive samples that Rubin Observatory's Legacy Survey of Space and Time will provide. △ Less

Submitted 29 November, 2023; originally announced November 2023.

Comments: Submitted for publication in MNRAS, posted to arXiv after responding to two positive rounds of referee comments. Key results in Figs 3, 5, 6 and 11

Radio relics in massive galaxy cluster mergers in the TNG-Cluster simulation

Authors: W. Lee, A. Pillepich, J. ZuHone, D. Nelson, M. J. Jee, D. Nagai, K. Finner

Abstract: Radio relics are diffuse synchrotron sources in the outskirts of merging galaxy clusters energized by the merger shocks. In this paper, we present an overview of the radio relics in massive cluster mergers identified in the new TNG-Cluster simulation. This is a suite of magnetohydrodynamical cosmological zoom-in simulations of 352 massive galaxy clusters with… ▽ More Radio relics are diffuse synchrotron sources in the outskirts of merging galaxy clusters energized by the merger shocks. In this paper, we present an overview of the radio relics in massive cluster mergers identified in the new TNG-Cluster simulation. This is a suite of magnetohydrodynamical cosmological zoom-in simulations of 352 massive galaxy clusters with $M_{\rm 500c}= 10^{14.0-15.3}\rm~M_{\odot}$ sampled from a 1 Gpc-size cosmological box. The simulations are performed using the moving-mesh code AREPO with the galaxy formation model and high numerical resolution consistent with the TNG300 run of the IllustrisTNG series. We post-process the shock properties obtained from the on-the-fly shock finder to estimate the diffuse radio emission generated by cosmological shockwaves for a total of $\sim300$ radio relics at redshift $z=0-1$. TNG-Cluster returns a variety of radio relics with diverse morphologies, encompassing textbook examples of double radio relics, single relics, and ``inverted" radio relics that are convex to the cluster center. Moreover, the simulated radio relics reproduce both the abundance and statistical relations of observed relics. We find that extremely large radio relics ($>$ 2 Mpc) are predominantly produced in massive cluster mergers with $M_{\rm 500c}\gtrsim8\times10^{14}~\rm~M_{\odot}$. This underscores the significance of simulating massive mergers to study giant radio relics similar to those found in observations. We release a library of radio relics from the TNG-Cluster simulation, which will serve as a crucial reference for upcoming next-generation surveys. △ Less

Submitted 11 March, 2024; v1 submitted 10 November, 2023; originally announced November 2023.

Comments: 24 pages, 16 figures, 2 appendices, Accepted for publication in A&A. See the TNG-Cluster website at www.tng-project.org/cluster/

An Atlas of Gas Motions in the TNG-Cluster Simulation: from Cluster Cores to the Outskirts

Authors: Mohammadreza Ayromlou, Dylan Nelson, Annalisa Pillepich, Eric Rohr, Nhut Truong, Yuan Li, Aurora Simionescu, Katrin Lehle, Wonki Lee

Abstract: Galaxy clusters are unique laboratories for studying astrophysical processes and their impact on gas kinematics. Despite their importance, the full complexity of gas motion within and around clusters remains poorly known. This paper is part of a series presenting first results from the new TNG-Cluster simulation, a suite of 352 massive clusters including the full cosmological context, mergers, acc… ▽ More Galaxy clusters are unique laboratories for studying astrophysical processes and their impact on gas kinematics. Despite their importance, the full complexity of gas motion within and around clusters remains poorly known. This paper is part of a series presenting first results from the new TNG-Cluster simulation, a suite of 352 massive clusters including the full cosmological context, mergers, accretion, baryonic processes, feedback, and magnetic fields. Studying the dynamics and coherence of gas flows, we find that gas motions in cluster cores and intermediate regions are largely balanced between inflows and outflows, exhibiting a Gaussian distribution centered at zero velocity. In the outskirts, even the net velocity distribution becomes asymmetric, featuring a double peak where the second peak reflects cosmic accretion. Across all cluster regions, the resulting net flow distribution reveals complex gas dynamics. These are strongly correlated with halo properties: at a given total cluster mass, unrelaxed, late-forming halos with less massive black holes and lower accretion rates exhibit a more dynamic behavior. Our analysis shows no clear relationship between line-of-sight and radial gas velocities, suggesting that line-of-sight velocity alone is insufficient to distinguish between inflowing and outflowing gas. Additional properties, such as temperature, can help break this degeneracy. A velocity structure function (VSF) analysis indicates more coherent gas motion in the outskirts and more disturbed kinematics towards halo centers. In all cluster regions, the VSF shows a slope close to the theoretical models of Kolmogorov (1/3), except within 50 kpc of the cluster cores, where the slope is significantly steeper. The outcome of TNG-Cluster broadly aligns with observations of the VSF of multiphase gas across different scales in galaxy clusters, ranging from 1 kpc to Megaparsec scales. △ Less

Submitted 10 November, 2023; originally announced November 2023.

Comments: Submitted to A&A. See the TNG-Cluster website at https://www.tng-project.org/cluster/

Introducing the TNG-Cluster Simulation: overview and physical properties of the gaseous intracluster medium

Authors: Dylan Nelson, Annalisa Pillepich, Mohammadreza Ayromlou, Wonki Lee, Katrin Lehle, Eric Rohr, Nhut Truong

Abstract: We introduce the new TNG-Cluster project, an addition to the IllustrisTNG suite of cosmological magnetohydrodynamical simulations of galaxy formation. Our objective is to significantly increase the statistical sampling of the most massive and rare objects in the Universe: galaxy clusters with log(M_200c / Msun) > 14.3 - 15.4 at z=0. To do so, we re-simulate 352 cluster regions drawn from a 1 Gpc v… ▽ More We introduce the new TNG-Cluster project, an addition to the IllustrisTNG suite of cosmological magnetohydrodynamical simulations of galaxy formation. Our objective is to significantly increase the statistical sampling of the most massive and rare objects in the Universe: galaxy clusters with log(M_200c / Msun) > 14.3 - 15.4 at z=0. To do so, we re-simulate 352 cluster regions drawn from a 1 Gpc volume, thirty-six times larger than TNG300, keeping entirely fixed the IllustrisTNG physical model as well as the numerical resolution. This new sample of hundreds of massive galaxy clusters enables studies of the assembly of high-mass ellipticals and their supermassive black holes (SMBHs), brightest cluster galaxies (BCGs), satellite galaxy evolution and environmental processes, jellyfish galaxies, intracluster medium (ICM) properties, cooling and active galactic nuclei (AGN) feedback, mergers and relaxedness, magnetic field amplification, chemical enrichment, and the galaxy-halo connection at the high-mass end, with observables from the optical to radio synchrotron and the Sunyaev-Zeldovich (SZ) effect, to X-ray emission, as well as their cosmological applications. We present an overview of the simulation, the cluster sample, selected comparisons to data, and a first look at the diversity and physical properties of our simulated clusters and their hot ICM. △ Less

Submitted 6 March, 2024; v1 submitted 10 November, 2023; originally announced November 2023.

Comments: A&A accepted. See companion papers (Ayromlou, Lee, Lehle, Rohr, Truong). Additional information and visuals are available on the TNG-Cluster website at https://www.tng-project.org/cluster/

The hot circumgalactic media of massive cluster satellites in the TNG-Cluster simulation: existence and detectability

Authors: Eric Rohr, Annalisa Pillepich, Dylan Nelson, Mohammadreza Ayromlou, Elad Zinger

Abstract: The most massive galaxy clusters in the Universe host hundreds of massive satellite galaxies~$\mstar\sim10^{10-12.5} msun$, but it is unclear if these satellites are able to retain their own gaseous halos. We analyze the evolution of $\approx90,000$ satellites of stellar mass $\sim10^{9-12.5} msun$ around 352 galaxy clusters of mass $\mvir\sim10^{14.3-15.4} msun$ at $z=0$ from the new TNG-Cluster… ▽ More The most massive galaxy clusters in the Universe host hundreds of massive satellite galaxies~$\mstar\sim10^{10-12.5} msun$, but it is unclear if these satellites are able to retain their own gaseous halos. We analyze the evolution of $\approx90,000$ satellites of stellar mass $\sim10^{9-12.5} msun$ around 352 galaxy clusters of mass $\mvir\sim10^{14.3-15.4} msun$ at $z=0$ from the new TNG-Cluster suite of cosmological magneto-hydrodynamical galaxy cluster simulations. The number of massive satellites per host increases with host mass, and the mass--richness relation broadly agrees with observations. A halo of mass $\mvirhost\sim10^{14.5} (10^{15}) msun$ hosts $\sim100 (300)$ satellites today. Only a minority of satellites retain some gas, hot or cold, and this fraction increases with stellar mass. Lower-mass satellites $\sim10^{9-10} msun$ are more likely to retain part of their cold interstellar medium, consistent with ram pressure preferentially removing hot extended gas first. At higher stellar masses $\sim10^{10.5-12.5} msun$, the fraction of gas-rich satellites increases to unity, and nearly all satellites retain a portion of their hot, spatially extended circumgalactic medium (CGM), despite the ejective activity of their supermassive black holes. According to TNG-Cluster, the CGM of these gaseous satellites can be seen in soft X-ray emission (0.5-2.0 keV) that is $\gtrsim10$~times brighter than the local background. This X-ray surface brightness excess around satellites extends to $\approx30-100$ kpc, and is strongest for galaxies with higher stellar masses and larger host-centric distances. Approximately $10$~percent of the soft X-ray emission in cluster outskirts $\approx0.75-1.5\rvir$ originates from satellites. The CGM of member galaxies reflects the dynamics of cluster-satellite interactions and contributes to the observationally inferred properties of the intracluster medium. △ Less

Submitted 16 April, 2024; v1 submitted 10 November, 2023; originally announced November 2023.

Comments: Submitted to A&A, comments welcome. See companion papers by Nelson et al., Lehle et al., Truong et al., Lee et al., and Ayromlou et al., and see the TNG-Cluster website at www.tng-project.org/cluster/. Reuploaded after acceptance for publication

Journal ref: A&A 686, A86 (2024)

X-ray-inferred kinematics of the core ICM in Perseus-like clusters: Insights from the TNG-Cluster simulation

Authors: Nhut Truong, Annalisa Pillepich, Dylan Nelson, Irina Zhuravleva, Wonki Lee, Mohammadreza Ayromlou, Katrin Lehle

Abstract: The intracluster medium (ICM) of galaxy clusters encodes the impact of the physical processes that shape these massive halos, including feedback from central supermassive black holes (SMBHs). In this study, we examine the gas thermodynamics, kinematics, and the effects of SMBH feedback on the core of Perseus-like galaxy clusters with a new simulation suite: TNG-Cluster. We first make a selection o… ▽ More The intracluster medium (ICM) of galaxy clusters encodes the impact of the physical processes that shape these massive halos, including feedback from central supermassive black holes (SMBHs). In this study, we examine the gas thermodynamics, kinematics, and the effects of SMBH feedback on the core of Perseus-like galaxy clusters with a new simulation suite: TNG-Cluster. We first make a selection of simulated clusters similar to Perseus based on the total mass and inner ICM properties, such as their cool-core nature. We identify 30 Perseus-like systems among the 352 TNG-Cluster halos at $z=0$. Many exhibit thermodynamical profiles and X-ray morphologies with disturbed features such as ripples, bubbles, and shock fronts that are qualitatively similar to X-ray observations of Perseus. To study observable gas motions, we generate XRISM mock X-ray observations and conduct a spectral analysis of the synthetic data. In agreement with existing Hitomi measurements, TNG-Cluster predicts subsonic gas turbulence in the central regions of Perseus-like clusters, with a typical line-of-sight velocity dispersion of 200 km/s. This implies that turbulent pressure contributes $< 10\%$ to the dominant thermal pressure. In TNG-Cluster, such low (inferred) values of ICM velocity dispersion coexist with high-velocity outflows and bulk motions of relatively small amounts of super-virial hot gas, moving up to thousands of km/s. However, detecting these outflows in observations may prove challenging due to their anisotropic nature and projection effects. Driven by SMBH feedback, such outflows are responsible for many morphological disturbances in the X-ray maps of cluster cores. They also increase both the inferred and intrinsic ICM velocity dispersion. This effect is somewhat stronger when velocity dispersion is measured from higher-energy lines. △ Less

Submitted 28 March, 2024; v1 submitted 10 November, 2023; originally announced November 2023.

Comments: 14 pages, 8 figures. Accepted to A&A. See companion papers (Ayromlou, Lee, Lehle, Nelson, Rohr) and additional information on the TNG-Cluster website www.tng-project.org/cluster/

Journal ref: A&A 686, A200 (2024)

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

Authors: Katrin Lehle, Dylan Nelson, Annalisa Pillepich, Nhut Truong, Eric Rohr

Abstract: We analyze the physical properties of the gaseous intracluster medium (ICM) at the center of massive galaxy clusters with TNG-Cluster, a new cosmological magnetohydrodynamical simulation. Our sample contains 352 simulated clusters spanning a halo mass range of $10^{14} < {\rm M}_{\rm 500c} / M_\odot < 2 \times 10^{15}$ at $z=0$. We focus on the proposed classification of clusters into cool-core (C… ▽ More We analyze the physical properties of the gaseous intracluster medium (ICM) at the center of massive galaxy clusters with TNG-Cluster, a new cosmological magnetohydrodynamical simulation. Our sample contains 352 simulated clusters spanning a halo mass range of $10^{14} < {\rm M}_{\rm 500c} / M_\odot < 2 \times 10^{15}$ at $z=0$. We focus on the proposed classification of clusters into cool-core (CC) and non-cool-core (NCC) populations, the $z=0$ distribution of cluster central ICM properties, and the redshift evolution of the CC cluster population. We analyze resolved structure and radial profiles of entropy, temperature, electron number density, and pressure. To distinguish between CC and NCC clusters, we consider several criteria: central cooling time, central entropy, central density, X-ray concentration parameter, and density profile slope. According to TNG-Cluster and with no a-priori cluster selection, the distributions of these properties are unimodal, whereby CCs and NCCs represent the two extremes. Across the entire TNG-Cluster sample at $z=0$ and based on central cooling time, the strong CC fraction is $f_{\rm SCC} = 24\%$, compared to $f_{\rm WCC} = 60\% $ and $f_{\rm NCC} = 16\%$ for weak and non-cool-cores, respectively. However, the fraction of CCs depends strongly on both halo mass and redshift, although the magnitude and even direction of the trends vary with definition. The abundant statistics of simulated high-mass clusters in TNG-Cluster enables us to match observational samples and make a comparison with data. The CC fractions from $z=0$ to $z=2$ are in broad agreement with observations, as are radial profiles of thermodynamical quantities, globally as well as split for CC versus NCC halos. TNG-Cluster can therefore be used as a laboratory to study the evolution and transformations of cluster cores due to mergers, AGN feedback, and other physical processes. △ Less

Submitted 10 November, 2023; originally announced November 2023.

Comments: Submitted to A&A, comments welcome. See the TNG-Cluster website at www.tng-project.org/cluster/ for more details

Journal ref: A&A 687, A129 (2024)

ERGO-ML: Comparing IllustrisTNG and HSC galaxy images via contrastive learning

Authors: Lukas Eisert, Connor Bottrell, Annalisa Pillepich, Rhythm Shimakawa, Vicente Rodriguez-Gomez, Dylan Nelson, Eirini Angeloudi, Marc Huertas-Company

Abstract: Modern cosmological hydrodynamical galaxy simulations provide tens of thousands of reasonably realistic synthetic galaxies across cosmic time. However, quantitatively assessing the level of realism of simulated universes in comparison to the real one is difficult. In this paper of the ERGO-ML series (Extracting Reality from Galaxy Observables with Machine Learning), we utilize contrastive learning… ▽ More Modern cosmological hydrodynamical galaxy simulations provide tens of thousands of reasonably realistic synthetic galaxies across cosmic time. However, quantitatively assessing the level of realism of simulated universes in comparison to the real one is difficult. In this paper of the ERGO-ML series (Extracting Reality from Galaxy Observables with Machine Learning), we utilize contrastive learning to directly compare a large sample of simulated and observed galaxies based on their stellar-light images. This eliminates the need to specify summary statistics and allows to exploit the whole information content of the observations. We produce survey-realistic galaxy mock datasets resembling real Hyper Suprime-Cam (HSC) observations using the cosmological simulations TNG50 and TNG100. Our focus is on galaxies with stellar masses between $10^9$ and $10^{12} M_\odot$ at $z=0.1-0.4$. This allows us to evaluate the realism of the simulated TNG galaxies in comparison to actual HSC observations. We apply the self-supervised contrastive learning method NNCLR to the images from both simulated and observed datasets (g, r, i - bands). This results in a 256-dimensional representation space, encoding all relevant observable galaxy properties. Firstly, this allows us to identify simulated galaxies that closely resemble real ones by seeking similar images in this multi-dimensional space. Even more powerful, we quantify the alignment between the representations of these two image sets, finding that the majority ($\gtrsim 70$ per cent) of the TNG galaxies align well with observed HSC images. However, a subset of simulated galaxies with larger sizes, steeper Sersic profiles, smaller Sersic ellipticities, and larger asymmetries appears unrealistic. We also demonstrate the utility of our derived image representations by inferring properties of real HSC galaxies using simulated TNG galaxies as the ground truth. △ Less

Submitted 11 April, 2024; v1 submitted 30 October, 2023; originally announced October 2023.

Comments: 26 pages, 17 figures. Comments are welcome. Highlights in Fig. 3, 7 and 13

Journal ref: MNRAS, Volume 528, Issue 4, March 2024, Pages 7411-7439

LEM All-Sky Survey: Soft X-ray Sky at Microcalorimeter Resolution

Authors: Ildar Khabibullin, Massimiliano Galeazzi, Akos Bogdan, Jenna M. Cann, Eugene Churazov, Klaus Dolag, Jeremy J. Drake, William Forman, Lars Hernquist, Dimitra Koutroumpa, Ralph Kraft, K. D. Kuntz, Maxim Markevitch, Dan McCammon, Anna Ogorzalek, Ryan Pfeifle, Annalisa Pillepich, Paul P. Plucinsky, Gabriele Ponti, Gerrit Schellenberger, Nhut Truong, Milena Valentini, Sylvain Veilleux, Stephan Vladutescu-Zopp, Q. Daniel Wang , et al. (1 additional authors not shown)

Abstract: The Line Emission Mapper (LEM) is an X-ray Probe with with spectral resolution ~2 eV FWHM from 0.2 to 2.5 keV and effective area >2,500 cm$^2$ at 1 keV, covering a 33 arcmin diameter Field of View with 15 arcsec angular resolution, capable of performing efficient scanning observations of very large sky areas and enabling the first high spectral resolution survey of the full sky. The LEM-All-Sky Su… ▽ More The Line Emission Mapper (LEM) is an X-ray Probe with with spectral resolution ~2 eV FWHM from 0.2 to 2.5 keV and effective area >2,500 cm$^2$ at 1 keV, covering a 33 arcmin diameter Field of View with 15 arcsec angular resolution, capable of performing efficient scanning observations of very large sky areas and enabling the first high spectral resolution survey of the full sky. The LEM-All-Sky Survey (LASS) is expected to follow the success of previous all sky surveys such as ROSAT and eROSITA, adding a third dimension provided by the high resolution microcalorimeter spectrometer, with each 15 arcsec pixel of the survey including a full 1-2 eV resolution energy spectrum that can be integrated over any area of the sky to provide statistical accuracy. Like its predecessors, LASS will provide both a long-lasting legacy and open the door to the unknown, enabling new discoveries and delivering the baseline for unique GO studies. No other current or planned mission has the combination of microcalorimeter energy resolution and large grasp to cover the whole sky while maintaining good angular resolution and imaging capabilities. LASS will be able to probe the physical conditions of the hot phases of the Milky Way at multiple scales, from emission in the Solar system due to Solar Wind Charge eXchange, to the interstellar and circumgalactic media, including the North Polar Spur and the Fermi/eROSITA bubbles. It will measure velocities of gas in the inner part of the Galaxy and extract the emissivity of the Local Hot Bubble. By maintaining the original angular resolution, LASS will also be able to study classes of point sources through stacking. For classes with ~$10^4$ objects, it will provide the equivalent of 1 Ms of high spectral resolution data. We describe the technical specifications of LASS and highlight the main scientific objectives that will be addressed. (Abridged) △ Less

Submitted 24 October, 2023; originally announced October 2023.

Comments: White Paper in support of a mission concept to be submitted for the 2023 NASA Astrophysics Probes opportunity. This White Paper will be updated when required. 30 pages, 25 figures

VERTICO and IllustrisTNG: The spatially resolved effects of environment on galactic gas

Authors: Adam R. H. Stevens, Toby Brown, Benedikt Diemer, Annalisa Pillepich, Lars Hernquist, Dylan Nelson, Yannick M. Bahé, Alessandro Boselli, Timothy A. Davis, Pascal J. Elahi, Sara L. Ellison, María J. Jiménez-Donaire, Ian D. Roberts, Kristine Spekkens, Vicente Villanueva, Adam B. Watts, Christine D. Wilson, Nikki Zabel

Abstract: It has been shown in previous publications that the TNG100 simulation quantitatively reproduces the observed reduction in each of the total atomic and total molecular hydrogen gas for galaxies within massive halos, i.e.~dense environments. In this Letter, we study how well TNG50 reproduces the resolved effects of a Virgo-like cluster environment on the gas surface densities of satellite galaxies w… ▽ More It has been shown in previous publications that the TNG100 simulation quantitatively reproduces the observed reduction in each of the total atomic and total molecular hydrogen gas for galaxies within massive halos, i.e.~dense environments. In this Letter, we study how well TNG50 reproduces the resolved effects of a Virgo-like cluster environment on the gas surface densities of satellite galaxies with $m_* > \! 10^9\,{\rm M}_\odot$ and ${\rm SFR} \! > 0.05\,{\rm M}_\odot\,{\rm yr}^{-1}$. We select galaxies in the simulation that are analogous to those in the HERACLES and VERTICO surveys, and mock-observe them to the common specifications of the data. Although TNG50 does not quantitatively match the observed gas surface densities in the centers of galaxies, the simulation does qualitatively reproduce the trends of gas truncation and central density suppression seen in VERTICO in both HI and H$_2$. This result promises that modern cosmological hydrodynamic simulations can be used to reliably model the post-infall histories of cluster satellite galaxies. △ Less

Submitted 11 October, 2023; originally announced October 2023.

Comments: 9 pages, 4 figures, accepted in ApJL

Exploring chemical enrichment of the intracluster medium with the Line Emission Mapper

Authors: François Mernier, Yuanyuan Su, Maxim Markevitch, Congyao Zhang, Aurora Simionescu, Elena Rasia, Sheng-Chieh Lin, Irina Zhuravleva, Arnab Sarkar, Ralph P. Kraft, Anna Ogorzalek, Mohammadreza Ayromlou, William R. Forman, Christine Jones, Joel N. Bregman, Stefano Ettori, Klaus Dolag, Veronica Biffi, Eugene Churazov, Ming Sun, John ZuHone, Ákos Bogdán, Ildar I. Khabibullin, Norbert Werner, Nhut Truong , et al. (5 additional authors not shown)

Abstract: Synthesized in the cores of stars and supernovae, most metals disperse over cosmic scales and are ultimately deposited well outside the gravitational potential of their host galaxies. Since their presence is well visible through their X-ray emission lines in the hot gas pervading galaxy clusters, measuring metal abundances in the intracluster medium (ICM) offers us a unique view of chemical enrich… ▽ More Synthesized in the cores of stars and supernovae, most metals disperse over cosmic scales and are ultimately deposited well outside the gravitational potential of their host galaxies. Since their presence is well visible through their X-ray emission lines in the hot gas pervading galaxy clusters, measuring metal abundances in the intracluster medium (ICM) offers us a unique view of chemical enrichment of the Universe as a whole. Despite extraordinary progress in the field thanks to four decades of X-ray spectroscopy using CCD (and gratings) instruments, understanding the precise stellar origins of the bulk of metals, and when the latter were mixed on Mpc scales, requires an X-ray mission capable of spatial, non-dispersive high resolution spectroscopy covering at least the soft X-ray band over a large field of view. In this White Paper, we demonstrate how the Line Emission Mapper (LEM) probe mission concept will revolutionize our current picture of the ICM enrichment. Specifically, we show that LEM will be able to (i) spatially map the distribution of ten key chemical elements out to the virial radius of a nearby relaxed cluster and (ii) measure metal abundances in serendipitously discovered high-redshift protoclusters. Altogether, these key observables will allow us to constrain the chemical history of the largest gravitationally bound structures of the Universe. They will also solve key questions such as the universality of the initial mass function (IMF) and the initial metallicity of the stellar populations producing these metals, as well as the relative contribution of asymptotic giant branch (AGB) stars, core-collapse, and Type Ia supernovae to enrich the cosmic web over Mpc scales. Concrete observing strategies are also briefly discussed. △ Less

Submitted 6 October, 2023; originally announced October 2023.

Comments: 19 pages. White paper for a mission concept to be submitted for the 2023 NASA Astrophysics Probes opportunity

IllustrisTNG in the HSC-SSP: image data release and the major role of mini mergers as drivers of asymmetry and star formation

Authors: Connor Bottrell, Hassen M. Yesuf, Gergö Popping, Kiyoaki Christopher Omori, Shenli Tang, Xuheng Ding, Annalisa Pillepich, Dylan Nelson, Lukas Eisert, Hua Gao, Andy D. Goulding, Boris S. Kalita, Wentao Luo, Jenny E. Greene, Jingjing Shi, John D. Silverman

Abstract: At fixed galaxy stellar mass, there is a clear observational connection between structural asymmetry and offset from the star forming main sequence, $Δ$SFMS. Herein, we use the TNG50 simulation to investigate the relative roles of major mergers (stellar mass ratios $μ\geq0.25$), minor ($0.1 \leq μ< 0.25$), and mini mergers ($0.01 \leq μ< 0.1$) in driving this connection amongst star forming galaxi… ▽ More At fixed galaxy stellar mass, there is a clear observational connection between structural asymmetry and offset from the star forming main sequence, $Δ$SFMS. Herein, we use the TNG50 simulation to investigate the relative roles of major mergers (stellar mass ratios $μ\geq0.25$), minor ($0.1 \leq μ< 0.25$), and mini mergers ($0.01 \leq μ< 0.1$) in driving this connection amongst star forming galaxies (SFGs). We use dust radiative transfer post-processing with SKIRT to make a large, public collection of synthetic Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) images of simulated TNG galaxies over $0.1\leq z \leq 0.7$ with $\log M_{\star} / \mathrm{M}_{\odot}\geq9$ ($\sim750$k images). Using their instantaneous SFRs, known merger histories/forecasts, and HSC-SSP asymmetries, we show (1) that TNG50 SFGs qualitatively reproduce the observed trend between $Δ$SFMS and asymmetry and (2) a strikingly similar trend emerges between $Δ$SFMS and the time-to-coalescence for mini mergers. Controlling for redshift, stellar mass, environment, and gas fraction, we show that individual mini merger events yield small enhancements in SFRs and asymmetries that are sustained on long timescales (at least $\sim3$ Gyr after coalescence, on average) -- in contrast to major/minor merger remnants which peak at much greater amplitudes but are consistent with controls only $\sim1$ Gyr after coalescence. Integrating the boosts in SFRs and asymmetries driven by $μ\geq0.01$ mergers since $z=0.7$ in TNG50 SFGs, we show that mini mergers are responsible for (i) $55$ per cent of all merger-driven star formation and (ii) $70$ per cent of merger-driven asymmetric structure. Due to their relative frequency and prolonged boost timescales, mini mergers dominate over their minor and major counterparts in driving star formation and asymmetry in SFGs. △ Less

Submitted 7 October, 2023; v1 submitted 28 August, 2023; originally announced August 2023.

Comments: 31 pages; 19 figures; MNRAS accepted; Image data available via the TNG website: www.tng-project.org/bottrell23

On the likelihoods of finding very metal-poor (and old) stars in the Milky Way's disc, bulge, and halo

Authors: Diego Sotillo-Ramos, Maria Bergemann, Jennifer K. S. Friske, Annalisa Pillepich

Abstract: Recent observational studies have uncovered a small number of very metal-poor stars with cold kinematics in the Galactic disc and bulge. However, their origins remain enigmatic. We select a total of 138 Milky Way (MW) analogs from the TNG50 cosmological simulation based on their $z=0$ properties: disky morphology, stellar mass, and local environment. In order to make more predictive statements for… ▽ More Recent observational studies have uncovered a small number of very metal-poor stars with cold kinematics in the Galactic disc and bulge. However, their origins remain enigmatic. We select a total of 138 Milky Way (MW) analogs from the TNG50 cosmological simulation based on their $z=0$ properties: disky morphology, stellar mass, and local environment. In order to make more predictive statements for the MW, we further limit the spatial volume coverage of stellar populations in galaxies to that targeted by the upcoming 4MOST high-resolution survey of the Galactic disc and bulge. We find that across all galaxies, $\sim$20 per cent of very metal-poor (${\rm [Fe/H]} < -2$) stars belong to the disk, with some analogs reaching 30 per cent. About 50$\pm$10 per cent of the VMP disc stars are, on average, older than 12.5 Gyr and $\sim$70$\pm$10 per cent come from accreted satellites. A large fraction of the VMP stars belong to the halo ($\sim$70) and have a median age of 12 Gyr. Our results with the TNG50 cosmological simulation confirm earlier findings with simulations of fewer individual galaxies, and suggest that the stellar disc of the Milky Way is very likely to host significant amounts of very- and extremely-metal-poor stars that, although mostly of ex situ origin, can also form in situ, reinforcing the idea of the existence of a primordial Galactic disc. △ Less

Submitted 26 July, 2023; originally announced July 2023.

Comments: Accepted by MNRAS. 4 figures

X-ray metal line emission from the hot circumgalactic medium: probing the effects of supermassive black hole feedback

Authors: Nhut Truong, Annalisa Pillepich, Dylan Nelson, Ákos Bogdán, Gerrit Schellenberger, Priyanka Chakraborty, William R. Forman, Ralph Kraft, Maxim Markevitch, Anna Ogorzalek, Benjamin D. Oppenheimer, Arnab Sarkar, Sylvain Veilleux, Mark Vogelsberger, Q. Daniel Wan, Norbert Werner, Irina Zhuravleva, John Zuhone

Abstract: We derive predictions from state-of-the-art cosmological galaxy simulations for the spatial distribution of the hot circumgalactic medium (CGM, ${\rm [0.1-1]R_{200c}}$) through its emission lines in the X-ray soft band ($[0.3-1.3]$ keV). In particular, we compare IllustrisTNG, EAGLE, and SIMBA and focus on galaxies with stellar mass $10^{10-11.6}\, \MSUN$ at $z=0$. The three simulation models retu… ▽ More We derive predictions from state-of-the-art cosmological galaxy simulations for the spatial distribution of the hot circumgalactic medium (CGM, ${\rm [0.1-1]R_{200c}}$) through its emission lines in the X-ray soft band ($[0.3-1.3]$ keV). In particular, we compare IllustrisTNG, EAGLE, and SIMBA and focus on galaxies with stellar mass $10^{10-11.6}\, \MSUN$ at $z=0$. The three simulation models return significantly different surface brightness radial profiles of prominent emission lines from ionized metals such as OVII(f), OVIII, and FeXVII as a function of galaxy mass. Likewise, the three simulations predict varying azimuthal distributions of line emission with respect to the galactic stellar planes, with IllustrisTNG predicting the strongest angular modulation of CGM physical properties at radial range ${\gtrsim0.3-0.5\,R_{200c}}$. This anisotropic signal is more prominent for higher-energy lines, where it can manifest as X-ray eROSITA-like bubbles. Despite different models of stellar and supermassive black hole (SMBH) feedback, the three simulations consistently predict a dichotomy between star-forming and quiescent galaxies at the Milky-Way and Andromeda mass range, where the former are X-ray brighter than the latter. This is a signature of SMBH-driven outflows, which are responsible for quenching star formation. Finally, we explore the prospect of testing these predictions with a microcalorimeter-based X-ray mission concept with a large field-of-view. Such a mission would probe the extended hot CGM via soft X-ray line emission, determine the physical properties of the CGM, including temperature, from the measurement of line ratios, and provide critical constraints on the efficiency and impact of SMBH feedback on the CGM. △ Less

Submitted 26 August, 2023; v1 submitted 3 July, 2023; originally announced July 2023.

Comments: 22 pages, 16 figures. Accepted for publication in MNRAS

Properties of the Line-of-Sight Velocity Field in the Hot and X-ray Emitting Circumgalactic Medium of Nearby Simulated Disk Galaxies

Authors: J. A. ZuHone, G. Schellenberger, A. Ogorzalek, B. D. Oppenheimer, J. Stern, A. Bogdan, N. Truong, M. Markevitch, A. Pillepich, D. Nelson, J. N. Burchett, I. Khabibullin, C. A. Kilbourne, R. P. Kraft, P. E. J. Nulsen, S. Veilleux, M. Vogelsberger, Q. D. Wang, I. Zhuravleva

Abstract: The hot, X-ray-emitting phase of the circumgalactic medium of massive galaxies is believed to be the reservoir of baryons from which gas flows onto the central galaxy and into which feedback from AGN and stars inject mass, momentum, energy, and metals. These effects shape the velocity fields of the hot gas, which can be observed via the Doppler shifting and broadening of emission lines by X-ray IF… ▽ More The hot, X-ray-emitting phase of the circumgalactic medium of massive galaxies is believed to be the reservoir of baryons from which gas flows onto the central galaxy and into which feedback from AGN and stars inject mass, momentum, energy, and metals. These effects shape the velocity fields of the hot gas, which can be observed via the Doppler shifting and broadening of emission lines by X-ray IFUs. In this work, we analyze the gas kinematics of the hot circumgalactic medium of Milky Way-mass disk galaxies from the TNG50 simulation with synthetic observations to determine how future instruments can probe this velocity structure. We find that the hot phase is often characterized by outflows from the disk driven by feedback processes, radial inflows near the galactic plane, and rotation, though in some systems the velocity field is more disorganized and turbulent. With a spectral resolution of $\sim$1 eV, fast and hot outflows ($\sim$200-500 km s$^{-1}$) can be measured, depending on the orientation of the galaxy on the sky. The rotation velocity of the hot phase ($\sim$100-200 km s$^{-1}$) can be measured using line shifts in edge-on galaxies, and is slower than that of colder gas phases but similar to stellar rotation velocities. By contrast, the slow inflows ($\sim$50-100 km s$^{-1}$) are difficult to measure in projection with these other components, but may be detected in multi-component spectral fits. We find that the velocity measured is sensitive to which emission lines are used. Measuring these flows will constrain theories of how the gas in these galaxies evolves. △ Less

Submitted 20 May, 2024; v1 submitted 3 July, 2023; originally announced July 2023.

Comments: 46 pages, 32 figures, accepted to ApJ

Mapping the imprints of stellar and AGN feedback in the circumgalactic medium with X-ray microcalorimeters

Authors: Gerrit Schellenberger, Ákos Bogdán, John A. ZuHone, Benjamin D. Oppenheimer, Nhut Truong, Ildar Khabibullin, Fred Jennings, Annalisa Pillepich, Joseph Burchett, Christopher Carr, Priyanka Chakraborty, Robert Crain, William Forman, Christine Jones, Caroline A. Kilbourne, Ralph P. Kraft, Maxim Markevitch, Daisuke Nagai, Dylan Nelson, Anna Ogorzalek, Scott Randall, Arnab Sarkar, Joop Schaye, Sylvain Veilleux, Mark Vogelsberger , et al. (2 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 29 April, 2024; v1 submitted 3 July, 2023; originally announced July 2023.

Comments: 41 pages, 19 figures, accepted for publication in ApJ

Unusual integrated metallicity profile of our Milky Way

Authors: Jianhui Lian, Maria Bergemann, Annalisa Pillepich, Gail Zasowski, Richard R. Lane

Abstract: The heavy element abundance profiles in galaxies place stringent constraint on galaxy growth and assembly history. Low-redshift galaxies generally have a negative metallicity gradient in their gas and stars. Such gradients are thought to be a natural manifestation of galaxy inside-out formation. As the Milky Way is currently the only spiral galaxy in which we can measure temporally-resolved chemic… ▽ More The heavy element abundance profiles in galaxies place stringent constraint on galaxy growth and assembly history. Low-redshift galaxies generally have a negative metallicity gradient in their gas and stars. Such gradients are thought to be a natural manifestation of galaxy inside-out formation. As the Milky Way is currently the only spiral galaxy in which we can measure temporally-resolved chemical abundances, it enables unique insights into the origin of metallicity gradients and their correlation with the growth history of galaxies. However, until now, these unique abundance profiles had not been translated into the integrated-light measurements needed to seamlessly compare with the general galaxy population. Here we report the first measurement of the light-weighted, integrated stellar metallicity profile of our Galaxy. We find that the integrated metallicity profile of the Milky Way has a '$\wedge$'-shape broken metallicity profile, with a mildly positive gradient inside a Galactocentric radius of 7 kpc and a steep negative gradient outside. This metallicity profile appears unusual when compared to Milky Way-mass star-forming galaxies observed in the MaNGA survey and simulated in the TNG50 cosmological simulation. The analysis of the TNG50 simulated galaxies suggests that the Milky Way's positive inner gradient may be due to an inside-out quenching process. The steep negative gradient in the outer disc, however, is challenging to explain in the simulations. Our results suggest the Milky Way may not be a typical spiral galaxy for its mass regarding metallicity distribution and thus offers insight into the variety of galaxy enrichment processes. △ Less

Submitted 24 June, 2023; originally announced June 2023.

Comments: 34 pages, 6 figures, published online in Nature Astronomy with open access on 22 June 2023. This is the version prior to the peer review. The published version is available here: https://rdcu.be/dfiTf

Resonant scattering of the OVII X-ray emission line in the circumgalactic medium of TNG50 galaxies

Authors: Dylan Nelson, Chris Byrohl, Anna Ogorzalek, Maxim Markevitch, Ildar Khabibullin, Eugene Churazov, Irina Zhuravleva, Akos Bogdan, Priyanka Chakraborty, Caroline Kilbourne, Ralph Kraft, Annalisa Pillepich, Arnab Sarkar, Gerrit Schellenberger, Yuanyuan Su, Nhut Truong, Stephan Vladutescu-Zopp, Nastasha Wijers

Abstract: We study the impact of resonantly scattered X-ray line emission on the observability of the hot circumgalactic medium (CGM) of galaxies. We apply a Monte Carlo radiative transfer post-processing analysis to the high-resolution TNG50 cosmological magnetohydrodynamical galaxy formation simulation. This allows us to model the resonant scattering of OVII(r) X-ray photons within the complex, multi-phas… ▽ More We study the impact of resonantly scattered X-ray line emission on the observability of the hot circumgalactic medium (CGM) of galaxies. We apply a Monte Carlo radiative transfer post-processing analysis to the high-resolution TNG50 cosmological magnetohydrodynamical galaxy formation simulation. This allows us to model the resonant scattering of OVII(r) X-ray photons within the complex, multi-phase, multi-scale CGM. The resonant transition of the OVII He-like triplet is one of the brightest, and most promising, X-ray emission lines for detecting the hot CGM and measuring its physical properties. We focus on galaxies with stellar masses 10 < log(M*/Msun) < 11 at z ~ 0. After constructing a model for OVII(r) emission from the central galaxy as well as from CGM gas, we forward model these intrinsic photons to derive observable surface brightness maps. We find that scattering significantly boosts the observable OVII(r) surface brightness of the extended and diffuse CGM. This enhancement can be large -- an order of magnitude on average at a distance of 200 projected kpc for high-mass M* = 10^10.7 Msun galaxies. The enhancement is larger for lower mass galaxies, and can even reach a factor of 100, across the extended CGM. Galaxies with higher star formation rates, AGN luminosities, and central OVII(r) luminosities all have larger scattering enhancements, at fixed stellar mass. Our results suggest that next-generation X-ray spectroscopic missions including XRISM, LEM, ATHENA, and HUBS -- which aim to detect the hot CGM in emission -- could specifically target halos with significant enhancements due to resonant scattering. △ Less

Submitted 8 June, 2023; originally announced June 2023.

Comments: Published in MNRAS. See https://www.lem-observatory.org/ and https://www.tng-project.org/ for more details; 2023MNRAS.522.3665N

ERGO-ML: Towards a robust machine learning model for inferring the fraction of accreted stars in galaxies from integral-field spectroscopic maps

Authors: Eirini Angeloudi, Jesús Falcón-Barroso, Marc Huertas-Company, Regina Sarmiento, Annalisa Pillepich, Daniel Walo-Martín, Lukas Eisert

Abstract: Quantifying the contribution of mergers to the stellar mass of galaxies is key for constraining the mechanisms of galaxy assembly across cosmic time. However, the mapping between observable galaxy properties and merger histories is not trivial: cosmological galaxy simulations are the only tools we have for calibration. We study the robustness of a simulation-based inference of the ex-situ stellar… ▽ More Quantifying the contribution of mergers to the stellar mass of galaxies is key for constraining the mechanisms of galaxy assembly across cosmic time. However, the mapping between observable galaxy properties and merger histories is not trivial: cosmological galaxy simulations are the only tools we have for calibration. We study the robustness of a simulation-based inference of the ex-situ stellar mass fraction of nearby galaxies to different observables -- integrated and spatially-resolved -- and to different galaxy formation models -- IllustrisTNG and EAGLE -- with Machine Learning. We find that at fixed simulation, the fraction of accreted stars can be inferred with very high accuracy, with an error $\sim5$ per cent (10 per cent) from 2D integral-field spectroscopic maps (integrated quantities) throughout the considered stellar mass range. A bias (> 5 per cent) and an increase in scatter by a factor of 2 are introduced when testing with a different simulation, revealing a lack of generalization to distinct galaxy-formation models. Interestingly, upon using only stellar mass and kinematics maps in the central galactic regions for training, we find that this bias is removed and the ex-situ stellar mass fraction can be recovered in both simulations with < 15 per cent scatter, independently of the training set's origin. This opens up the door to a potential robust inference of the accretion histories of galaxies from existing Integral Field Unit surveys, such as MaNGA, covering a similar field of view (FOV) and containing spatially-resolved spectra for tens of thousands of nearby galaxies. △ Less

Submitted 1 June, 2023; originally announced June 2023.

Comments: 23 pages, 15 figures. Accepted for publication in MNRAS

A Virgo Environmental Survey Tracing Ionised Gas Emission (VESTIGE) XV. The Halpha luminosity function of the Virgo cluster

Authors: A. Boselli, M. Fossati, P. Cote, J. C. Cuillandre, L. Ferrarese, S. Gwyn, P. Amram, M. Ayromlou, M. Balogh, G. Bellusci, M. Boquien, G. Gavazzi, G. Hensler, A. Longobardi, D. Nelson, A. Pillepich, J. Roediger, R. Sanchez-Jansen, M. Sun, G. Trinchieri

Abstract: We use a complete set of deep narrow-band imaging data for 384 galaxies gathered during the VESTIGE survey to derive the first Halpha luminosity function (LF) of the Virgo cluster within R200. The data allow us to cover the whole dynamic range of the Halpha LF (10^36<LHa<10^42 erg s^-1). After they are corrected for [NII] contamination and dust attenuation, the data are used to derive the SFR func… ▽ More We use a complete set of deep narrow-band imaging data for 384 galaxies gathered during the VESTIGE survey to derive the first Halpha luminosity function (LF) of the Virgo cluster within R200. The data allow us to cover the whole dynamic range of the Halpha LF (10^36<LHa<10^42 erg s^-1). After they are corrected for [NII] contamination and dust attenuation, the data are used to derive the SFR function in the range 10^-4<SFR<10 Mo yr^-1. These LF are compared to those derived at other frequencies or using different tracers of star formation in Virgo, in other nearby and high-z clusters, in the field, and to those predicted by the IllustrisTNG cosmological hydrodynamical simulations. The Halpha LF of the Virgo cluster is fairly flat (a=-1.07) in the range 10^38.5<LHa<10^40.5 erg s^-1, and it abruptly decreases at lower luminosities. When compared to those derived for other nearby clusters and for the field, the slope and the characteristic luminosity of the Schechter function change as a function of the dynamical mass of the system, of the temperature of the X-rays gas, and of the dynamical pressure exerted on the interstellar medium of galaxies moving at high velocity within the intracluster medium. All these trends can be explained in a scenario in which the activity of SF is reduced in massive clusters due to their hydrodynamical interaction with the surrounding medium, suggesting once again that ram-pressure stripping is the dominant mechanism affecting galaxy evolution in local clusters of dynamical mass M200>10^14 Mo. The comparison with the IllustrisTNG cosmological hydrodynamical simulations shows a more pronounced decrease at the faint end of the distribution. If Virgo is representative of typical nearby clusters of similar mass, this difference suggests that the stripping process in simulated galaxies in these environments is more efficient than observed. △ Less

Submitted 25 May, 2023; originally announced May 2023.

Comments: Accepted for publication on A&A

Journal ref: A&A 675, A123 (2023)

The signature of galaxy formation models in the power spectrum of the hydrogen 21cm line during reionization

Authors: Joseph S. W. Lewis, Annalisa Pillepich, Dylan Nelson, Ralf S. Klessen, Simon C. O. Glover

Abstract: Observations of the 21cm line of neutral hydrogen are poised to revolutionize our knowledge of cosmic reionization and the high-redshift population of galaxies. However, harnessing such information requires robust and comprehensive theoretical modeling. We study the non-linear effects of hydrodynamics and astrophysical feedback processes, including stellar and AGN feedback, on the 21cm signal by p… ▽ More Observations of the 21cm line of neutral hydrogen are poised to revolutionize our knowledge of cosmic reionization and the high-redshift population of galaxies. However, harnessing such information requires robust and comprehensive theoretical modeling. We study the non-linear effects of hydrodynamics and astrophysical feedback processes, including stellar and AGN feedback, on the 21cm signal by post-processing three existing cosmological hydrodynamical simulations of galaxy formation: Illustris, IllustrisTNG, and Eagle. Overall and despite their different underlying galaxy-formation models, the three simulations return similar predictions for the global 21cm rightness temperature and its power spectrum. At fixed redshift, most differences are attributable to differences in the history of reionization, in turn driven by differences in the build-up of stellar sources of radiation. However, the impact of astrophysics is imprinted in the 21cm power spectrum through several unique signatures. First, we find significant small scale ($k \geq 10\, \rm {Mpc}^{-1}$) differences between Illustris and IllustrisTNG, where higher velocity winds generated by supernova feedback soften density peaks and lead to lower 21cm power in TNG. Second, we find more 21cm power at intermediate scales ($k \approx 0.8\, \rm {Mpc}^{-1}$) in Eagle, due to differences in ionization driven by highly effective stellar feedback, leading to lower star formation, older and redder stellar populations, and thus lower ionizing luminosities. Though subtle, these features could allow future observations of the 21cm signal, in conjunction with other reionization observables, to constrain theoretical models for galactic feedback at high redshift. △ Less

Submitted 18 December, 2023; v1 submitted 16 May, 2023; originally announced May 2023.

Comments: submitted to MNRAS (12/05/23); Accepted 07/12/23 after minor revision

The First Quiescent Galaxies in TNG300

Authors: Abigail I. Hartley, Erica J. Nelson, Katherine A. Suess, Alex M. Garcia, Minjung Park, Lars Hernquist, Rachel Bezanson, Rebecca Nevin, Annalisa Pillepich, Aimee L. Schechter, Bryan A. Terrazas, Paul Torrey, Sarah Wellons, Katherine E. Whitaker, Christina C. Williams

Abstract: We identify the first quiescent galaxies in TNG300, the largest volume of the IllustrisTNG cosmological simulation suite, and explore their quenching processes and time evolution to z=0. We find that the first quiescent galaxies with stellar masses M_* > 3 x 10^{10} M_sun and specific star formation rates sSFR < 10^{-11} yr^{-1} emerge at z~4.2 in TNG300. Suppression of star formation in these gal… ▽ More We identify the first quiescent galaxies in TNG300, the largest volume of the IllustrisTNG cosmological simulation suite, and explore their quenching processes and time evolution to z=0. We find that the first quiescent galaxies with stellar masses M_* > 3 x 10^{10} M_sun and specific star formation rates sSFR < 10^{-11} yr^{-1} emerge at z~4.2 in TNG300. Suppression of star formation in these galaxies begins with a thermal mode of AGN feedback at z~6, and a kinetic feedback mode acts in each galaxy by z~4.7 to complete the quenching process, which occurs on a time-scale of ~0.35 Gyr. Surprisingly, we find that the majority of these galaxies are not the main progenitors of their z=0 descendants; instead, four of the five galaxies fall into more massive galaxies in subsequent mergers at a range of redshifts 2.5 < z < 0.2. By z=0, these descendants are the centres of galaxy clusters with average stellar masses of 8 x 10^{11} M_sun. We make predictions for the first quenched galaxies to be located by the James Webb Space Telescope (JWST). △ Less

Submitted 18 April, 2023; originally announced April 2023.

Comments: 6 pages, 4 figures

Jellyfish galaxies with the IllustrisTNG simulations -- Citizen-science results towards large distances, low-mass hosts, and high redshifts

Authors: Elad Zinger, Gandhali Joshi, Annalisa Pillepich, Eric Rohr, Dylan Nelson

Abstract: We present the ``Cosmological Jellyfish'' project - a citizen-science classification program to identify jellyfish galaxies within the IllustrisTNG cosmological simulations. Jellyfish (JF) are satellite galaxies that exhibit long trailing gas features -- `tails' -- extending from their stellar body. Their distinctive morphology arises due to ram-pressure stripping (RPS) as they move through the ba… ▽ More We present the ``Cosmological Jellyfish'' project - a citizen-science classification program to identify jellyfish galaxies within the IllustrisTNG cosmological simulations. Jellyfish (JF) are satellite galaxies that exhibit long trailing gas features -- `tails' -- extending from their stellar body. Their distinctive morphology arises due to ram-pressure stripping (RPS) as they move through the background gaseous medium. Using the TNG50 and TNG100 simulations, we construct a sample of $\sim 80,000$ satellite galaxies spanning an unprecedented range of stellar masses, $10^{8.3-12.3}\,\mathrm{M_\odot}$, and host masses of $M_\mathrm{200,c}=10^{10.4-14.6}\,\mathrm{M_\odot}$ back to $z=2$ \citep[extending the work of][]{yun_jellyfish_2019}. Based on this sample, $\sim 90,000$ galaxy images were presented to volunteers in a citizen-science project on the Zooniverse platform who were asked to determine if each galaxy image resembles a jellyfish. Based on volunteer votes, each galaxy was assigned a score determining if it is a JF or not. This paper describes the project, the inspected satellite sample, the methodology, and the classification process that resulted in a dataset of $5,307$ visually-identified jellyfish galaxies. We find that JF galaxies are common in nearly all group- and cluster-sized systems, with the JF fraction increasing with host mass and decreasing with satellite stellar mass. We highlight JF galaxies in three relatively unexplored regimes: low-mass hosts of $M_\mathrm{200,c}\sim10^{11.5-13}\,\mathrm{M_\odot}$, radial positions within hosts exceeding the virial radius $R_\mathrm{200,c}$, and at high redshift up to $z=2$. The full dataset of our jellyfish scores is publicly available and can be used to select and study JF galaxies in the IllustrisTNG simulations. △ Less

Submitted 18 April, 2023; originally announced April 2023.

Comments: submitted to MNRAS ; See additional jellyfish companion papers today on astro-ph: Rohr et al. and Goeller et al.; Jellyfish image gallery: https://www.tng-project.org/explore/gallery/zinger23/

Jellyfish galaxies with the IllustrisTNG simulations -- No enhanced population-wide star formation according to TNG50

Authors: Junia Göller, Gandhali Joshi, Eric Rohr, Elad Zinger, Annalisa Pillepich

Abstract: Due to ram-pressure stripping, jellyfish galaxies are thought to lose large amounts, if not all, of their interstellar medium. Nevertheless, some, but not all, observations suggest that jellyfish galaxies exhibit enhanced star formation compared to control samples, even in their ram pressure-stripped tails. We use the TNG50 cosmological gravity+magnetohydrodynamical simulation, with an average spa… ▽ More Due to ram-pressure stripping, jellyfish galaxies are thought to lose large amounts, if not all, of their interstellar medium. Nevertheless, some, but not all, observations suggest that jellyfish galaxies exhibit enhanced star formation compared to control samples, even in their ram pressure-stripped tails. We use the TNG50 cosmological gravity+magnetohydrodynamical simulation, with an average spatial resolution of 50-200 pc in the star-forming regions of galaxies, to quantify the star formation activity and rates (SFRs) of more than 700 jellyfish galaxies at $z=0-1$ with stellar masses $10^{8.3-10.8}\,\mathrm{M}_\odot$ in hosts with mass $10^{10.5-14.3}\,\mathrm{M}_\odot$. We extract their global SFRs, the SFRs within their main stellar body vs. within the tails, and we follow the evolution of the star formation along their individual evolutionary tracks. We compare the findings for jellyfish galaxies to those of diversely-constructed control samples, including against satellite and field galaxies with matched redshift, stellar mass, gas fraction and host halo mass. According to TNG50, star formation and ram-pressure stripping can indeed occur simultaneously within any given galaxy, and frequently do so. Moreover, star formation can also take place within the ram pressure-stripped tails, even though the latter is typically subdominant. However, TNG50 does not predict elevated population-wide SFRs in jellyfish compared to analog satellite galaxies with the same stellar mass or gas fraction. Simulated jellyfish galaxies do undergo bursts of elevated star formation along their history but, at least according to TNG50, these do not translate into a population-wide enhancement at any given epoch. △ Less

Submitted 2 November, 2023; v1 submitted 18 April, 2023; originally announced April 2023.

Comments: 20 pages, 10 figures, 1 table, accepted to MNRAS. See additional jellyfish companion papers today on astro-ph: Zinger et al. and Rohr et al

Jellyfish galaxies with the IllustrisTNG simulations -- When, where, and for how long does ram pressure stripping of cold gas occur?

Authors: Eric Rohr, Annalisa Pillepich, Dylan Nelson, Elad Zinger, Gandhali Joshi, Mohommadreza Ayromlou

Abstract: Jellyfish galaxies are prototypical examples of satellite galaxies undergoing strong ram pressure stripping (RPS). We analyze the evolution of 512 unique, first-infalling jellyfish galaxies from the TNG50 cosmological simulation. These have been visually inspected to be undergoing RPS sometime in the past 5 billion years (since $z=0.5$), have satellite stellar masses… ▽ More Jellyfish galaxies are prototypical examples of satellite galaxies undergoing strong ram pressure stripping (RPS). We analyze the evolution of 512 unique, first-infalling jellyfish galaxies from the TNG50 cosmological simulation. These have been visually inspected to be undergoing RPS sometime in the past 5 billion years (since $z=0.5$), have satellite stellar masses $\mstarsat\sim10^{8-10.5}\,M_\odot$, and live in hosts with $\mvir\sim10^{12-14.3}\,M_\odot$ at $z=0$. We quantify the cold gas ($T\leq10^{4.5}$ K) removal using the tracer particles, confirming that for these jellyfish, RPS is the dominant driver of cold gas loss after infall. Half of these jellyfish are completely gas-less by $z=0$, and these galaxies have earlier infall times and smaller satellite-to-host mass ratios than their gaseous counterparts. RPS can act on jellyfish galaxies over long time scales of $\approx1.5-8$ Gyr. Jellyfish in more massive hosts are impacted by RPS for a shorter time span and, at a fixed host mass, jellyfish with less cold gas at infall and lower stellar masses at $z=0$ have shorter RPS time spans. While RPS may act for long periods of time, the peak RPS period -- where at least 50 pecent of the total RPS occurs -- begins within $\approx1$ Gyr of infall and lasts $\lesssim2$ Gyr. During this period, the jellyfish are at host-centric distances $\sim0.2-2\rvir$, illustrating that much of RPS occurs at large distances from the host galaxy. Interestingly, jellyfish continue forming stars until they have lost $\approx98$ percent of their cold gas. For groups and clusters in TNG50 $(\mvirhost\sim10^{13-14.3}\,M_\odot)$, jellyfish galaxies deposit more cold gas ($\sim10^{11-12}\,M_\odot$) into halos than exist in them at $z=0$, demonstrating that jellyfish, and in general satellite galaxies, are a significant source of cold gas accretion. △ Less

Submitted 16 April, 2024; v1 submitted 18 April, 2023; originally announced April 2023.

Comments: 20 pages, 11 figures + 3 appendices with 4 figures. Submitted to MNRAS. Key figures are 2, 8, 9, 11. See additional jellyfish companion papers today on astro-ph: Zinger+ and Goeller+. All data used in this publication, including the Cosmological Jellyfish Project results, are publicaly available. Reuploaded after publication

Journal ref: MNRAS, 524, 3502-3525 (2023)

Disk flaring with TNG50: diversity across Milky Way and M31 analogs

Authors: Diego Sotillo-Ramos, Martina Donnari, Annalisa Pillepich, Neige Frankel, Dylan Nelson, Volker Springel, Lars Hernquist

Abstract: We use the sample of 198 Milky Way (MW) and Andromeda (M31) analogs from TNG50 to quantify the level of disk flaring predicted by a modern, high-resolution cosmological hydrodynamical simulation. Disk flaring refers to the increase of vertical stellar disk height with galactocentric distance. The TNG50 galaxies are selected to have stellar disky morphology, a stellar mass in the range of… ▽ More We use the sample of 198 Milky Way (MW) and Andromeda (M31) analogs from TNG50 to quantify the level of disk flaring predicted by a modern, high-resolution cosmological hydrodynamical simulation. Disk flaring refers to the increase of vertical stellar disk height with galactocentric distance. The TNG50 galaxies are selected to have stellar disky morphology, a stellar mass in the range of $M_* = 10^{10.5 - 11.2}~\rm{M_{\odot}}$, and a MW-like Mpc-scale environment at $z=0$. The stellar disks of such TNG50 MW/M31 analogs exhibit a wide diversity of structural properties, including a number of galaxies with disk scalelength and thin and thick disk scaleheights that are comparable to those measured or inferred for the Galaxy and Andromeda. With one set of physical ingredients, TNG50 returns a large variety of flaring flavours and amounts, also for mono-age stellar populations. With this paper, we hence propose a non-parametric characterization of flaring. The typical MW/M31 analogs exhibit disk scaleheights that are $1.5-2$ times larger in the outer than in the inner regions of the disk for both old and young stellar populations, but with a large galaxy-to-galaxy variation. Which stellar population flares more, and by how much, also varies from galaxy to galaxy. TNG50 de facto brackets existing observational constraints for the Galaxy and all previous numerical findings. A link between the amount of flaring and the $z=0$ global galaxy structural properties or merger history is complex. However, a connection between the scaleheights and the local stellar vertical kinematics and gravitational potential is clearly in place. △ Less

Submitted 28 March, 2023; originally announced March 2023.

Comments: Submitted to MNRAS. Main figures: 5, 13. See presentation and data release of TNG50 MW/M31 analogs by Pillepich et al. and see also Ramesh et al. on astro-ph today

Milky Way and Andromeda analogs from the TNG50 simulation

Authors: Annalisa Pillepich, Diego Sotillo-Ramos, Rahul Ramesh, Dylan Nelson, Christoph Engler, Vicente Rodriguez-Gomez, Martin Fournier, Martina Donnari, Volker Springel, Lars Hernquist

Abstract: We present the properties of Milky Way- and Andromeda-like (MW/M31-like) galaxies simulated within TNG50, the highest-resolution run of the IllustrisTNG suite of $Λ$CDM magneto-hydrodynamical simulations. We introduce our fiducial selection for MW/M31 analogs, which we propose for direct usage as well as for reference in future analyses. TNG50 contains 198 MW/M31 analogs, i.e. galaxies with stella… ▽ More We present the properties of Milky Way- and Andromeda-like (MW/M31-like) galaxies simulated within TNG50, the highest-resolution run of the IllustrisTNG suite of $Λ$CDM magneto-hydrodynamical simulations. We introduce our fiducial selection for MW/M31 analogs, which we propose for direct usage as well as for reference in future analyses. TNG50 contains 198 MW/M31 analogs, i.e. galaxies with stellar disky morphology, with a stellar mass in the range of $M_* = 10^{10.5 - 11.2}$ Msun, and within a MW-like Mpc-scale environment at z=0. These are resolved with baryonic (dark matter) mass resolution of $8.5\times10^4$ Msun ($4.5\times10^5$ Msun) and $\sim150$ pc of average spatial resolution in the star-forming regions: we therefore expand by many factors (2 orders of magnitude) the sample size of cosmologically-simulated analogs with similar ($\times 10$ better) numerical resolution. The majority of TNG50 MW/M31 analogs at $z=0$ exhibit a bar, 60 per cent are star-forming, the sample includes 3 Local Group (LG)-like systems, and a number of galaxies host one or more satellites as massive as e.g. the Magellanic Clouds. Even within such a relatively narrow selection, TNG50 reveals a great diversity in galaxy and halo properties, as well as in past histories. Within the TNG50 sample, it is possible to identify several simulated galaxies whose integral and structural properties are consistent, one or more at a time, with those measured for the Galaxy and Andromeda. With this paper, we document and release a series of broadly applicable data products that build upon the IllustrisTNG public release and aim to facilitate easy access and analysis by public users. These include datacubes across snapshots ($0 \le z \le 7$) for each TNG50 MW/M31-like galaxy, and a series of value-added catalogs that will be continually expanded to provide a convenient and up to date community resource. △ Less

Submitted 28 March, 2023; originally announced March 2023.

Comments: Submitted to MNRAS. Feedback from the community encouraged. This is also a data-release paper: see visuals and documentation at https://www.tng-project.org/data/milkyway+andromeda/. See other TNG50 MW/M31-based papers also on astro-ph today: Ramesh+ and Sotillo-Ramos+

The Circumgalactic Medium of Milky Way-like Galaxies in the TNG50 Simulation -- II: Cold, Dense Gas Clouds and High-Velocity Cloud Analogs

Authors: Rahul Ramesh, Dylan Nelson, Annalisa Pillepich

Abstract: We use the TNG50 simulation of the IllustrisTNG project to study cold, dense clouds of gas in the circumgalactic media (CGM) of Milky Way-like galaxies. We find that their CGM is typically filled with of order one hundred (thousand) reasonably (marginally) resolved clouds, possible analogs of high-velocity clouds (HVCs). There is a large variation in cloud abundance from galaxy to galaxy, and the… ▽ More We use the TNG50 simulation of the IllustrisTNG project to study cold, dense clouds of gas in the circumgalactic media (CGM) of Milky Way-like galaxies. We find that their CGM is typically filled with of order one hundred (thousand) reasonably (marginally) resolved clouds, possible analogs of high-velocity clouds (HVCs). There is a large variation in cloud abundance from galaxy to galaxy, and the physical properties of clouds that we explore -- mass, size, metallicity, pressure, and kinematics -- are also diverse. We quantify the distributions of cloud properties and cloud-background contrasts, providing cosmological inputs for idealized simulations. Clouds characteristically have sub-solar metallicities, diverse shapes, small overdensities ($χ= n_{\rm cold} / n_{\rm hot} \lesssim 10$), are mostly inflowing, and have sub-virial rotation. At TNG50 resolution, resolved clouds have median masses of $\sim 10^6\,\rm{M_\odot}$ and sizes of $\sim 10$ kpc. Larger clouds are well converged numerically, while the abundance of the smallest clouds increases with resolution, as expected. In TNG50 MW-like haloes, clouds are slightly (severely) under-pressurised relative to their surroundings with respect to total (thermal) pressure, implying that magnetic fields may be important. Clouds are not distributed uniformly throughout the CGM, but are clustered around other clouds, often near baryon-rich satellite galaxies. This suggests that at least some clouds originate from satellites, via direct ram-pressure stripping or otherwise. Finally, we compare with observations of intermediate and high velocity clouds from the real Milky Way halo. TNG50 shows a similar cloud velocity distribution as observations, and predicts a significant population of currently difficult-to-detect low velocity clouds. △ Less

Submitted 28 March, 2023; originally announced March 2023.

Comments: Accepted for publication (MNRAS). Part of a set of papers based on TNG50 MW/M31-like galaxies. Additional visuals and data products at www.tng-project.org/ramesh23b

The growth of brightest cluster galaxies in the TNG300 simulation:dissecting the contributions from mergers and in situ star formation

Authors: Daniel Montenegro-Taborda, Vicente Rodriguez-Gomez, Annalisa Pillepich, Vladimir Avila-Reese, Laura V. Sales, Aldo Rodríguez-Puebla, Lars Hernquist

Abstract: We investigate the formation of brightest cluster galaxies (BCGs) in the TNG300 cosmological simulation of the IllustrisTNG project. Our cluster sample consists of 700 haloes with $M_{200} \geq 5 \times 10^{13} \, \mathrm{M}_{\odot}$ at $z=0$, along with their progenitors at earlier epochs. This includes 280 systems with $M_{200} \geq 10^{14} \, \mathrm{M}_{\odot}$ at $z=0$, as well as three haloe… ▽ More We investigate the formation of brightest cluster galaxies (BCGs) in the TNG300 cosmological simulation of the IllustrisTNG project. Our cluster sample consists of 700 haloes with $M_{200} \geq 5 \times 10^{13} \, \mathrm{M}_{\odot}$ at $z=0$, along with their progenitors at earlier epochs. This includes 280 systems with $M_{200} \geq 10^{14} \, \mathrm{M}_{\odot}$ at $z=0$, as well as three haloes with $M_{200} \geq 10^{15} \, \mathrm{M}_{\odot}$. We find that the stellar masses and star formation rates of our simulated BCGs are in good agreement with observations at $z \lesssim 0.4$, and that they have experienced, on average, $\sim$2 ($\sim$3) major mergers since $z=1$ ($z=2$). Separating the BCG from the intracluster light (ICL) by means of a fixed 30 kpc aperture, we find that the fraction of stellar mass contributed by ex situ (i.e. accreted) stars at $z=0$ is approximately 70, 80, and 90 per cent for the BCG, BCG+ICL, and ICL, respectively. Tracking our simulated BCGs back in time using the merger trees, we find that they became dominated by ex situ stars at $z \sim $1-2, and that half of the stars that are part of the BCG at $z=0$ formed early ($z \sim 3$) in other galaxies, but `assembled' onto the BCG until later times ($z \approx 0.8$ for the whole sample, $z \approx 0.5$ for BCGs in $M_{200} \geq 5 \times 10^{14} \, \mathrm{M}_{\odot}$ haloes). Finally, we show that the stellar mass profiles of BCGs are often dominated by ex situ stars at all radii, with stars from major mergers being found closer to the centre, while stars that were tidally stripped from other galaxies dominate the outer regions. △ Less

Submitted 21 February, 2023; originally announced February 2023.

Comments: 19 pages, 11 figures. Accepted for publication in MNRAS

On the nature of disks at high redshift seen by JWST/CEERS with contrastive learning and cosmological simulations

Authors: J. Vega-Ferrero, M. Huertas-Company, L. Costantin, P. G. Pérez-González, R. Sarmiento, J. S. Kartaltepe, A. Pillepich, M. B. Bagley, S. L. Finkelstein, E. J. McGrath, J. H. Knapen, P. Arrabal Haro, E. F. Bell, F. Buitrago, A. Calabrò, A. Dekel, M. Dickinson, H. Domínguez Sánchez, D. Elbaz, H. C. Ferguson, M. Giavalisco, B. W. Holwerda, D. D. Kocesvski, A. M. Koekemoer, V. Pandya , et al. (4 additional authors not shown)

Abstract: Visual inspections of the first optical rest-frame images from JWST have indicated a surprisingly high fraction of disk galaxies at high redshifts. Here, we alternatively apply self-supervised machine learning to explore the morphological diversity at $z \geq 3$. Our proposed data-driven representation scheme of galaxy morphologies, calibrated on mock images from the TNG50 simulation, is shown to… ▽ More Visual inspections of the first optical rest-frame images from JWST have indicated a surprisingly high fraction of disk galaxies at high redshifts. Here, we alternatively apply self-supervised machine learning to explore the morphological diversity at $z \geq 3$. Our proposed data-driven representation scheme of galaxy morphologies, calibrated on mock images from the TNG50 simulation, is shown to be robust to noise and to correlate well with the physical properties of the simulated galaxies, including their 3D structure. We apply the method simultaneously to F200W and F356W galaxy images of a mass-complete sample ($M_*/M_\odot>10^9$) at $ 3 \leq z \leq 6$ from the first JWST/NIRCam CEERS data release. We find that the simulated and observed galaxies do not exactly populate the same manifold in the representation space from contrastive learning. We also find that half the galaxies classified as disks -- either CNN-based or visually -- populate a similar region of the representation space as TNG50 galaxies with low stellar specific angular momentum and non-oblate structure. Although our data-driven study does not allow us to firmly conclude on the true nature of these galaxies, it suggests that the disk fraction at $z \geq 3$ remains uncertain and possibly overestimated by traditional supervised classifications. Deeper imaging and spectroscopic follow-ups as well as comparisons with other simulations will help to unambiguously determine the true nature of these galaxies, and establish more robust constraints on the emergence of disks at very high redshift. △ Less

Submitted 25 October, 2023; v1 submitted 14 February, 2023; originally announced February 2023.

Comments: 28 pages, 22 figures. Accepted for publication in ApJ

The Origin of Stars in the Inner 500 Parsecs in TNG50 Galaxies

Authors: Alina Boecker, Nadine Neumayer, Annalisa Pillepich, Neige Frankel, Rahul Ramesh, Ryan Leaman, Lars Hernquist

Abstract: We investigate the origin of stars in the innermost $500\,\mathrm{pc}$ of galaxies spanning stellar masses of $5\times10^{8-12}\,\mathrm{M}_{\odot}$ at $\mathrm{z=0}$ using the cosmological magnetohydrodynamical TNG50 simulation. Three different origins of stars comprise galactic centers: 1) in-situ (born in the center), 2) migrated (born elsewhere in the galaxy and ultimately moved to the center)… ▽ More We investigate the origin of stars in the innermost $500\,\mathrm{pc}$ of galaxies spanning stellar masses of $5\times10^{8-12}\,\mathrm{M}_{\odot}$ at $\mathrm{z=0}$ using the cosmological magnetohydrodynamical TNG50 simulation. Three different origins of stars comprise galactic centers: 1) in-situ (born in the center), 2) migrated (born elsewhere in the galaxy and ultimately moved to the center), 3) ex-situ (accreted from other galaxies). In-situ and migrated stars dominate the central stellar mass budget on average with 73% and 23% respectively. The ex-situ fraction rises above 1% for galaxies $\gtrsim10^{11}\,\mathrm{M}_{\odot}$. Yet, only 9% of all galaxies exhibit no ex-situ stars in their centers and the scatter of ex-situ mass is significant ($4-6\,\mathrm{dex}$). Migrated stars predominantly originate closely from the center ($1-2\,\mathrm{kpc}$), but if they travelled together in clumps distances reach $\sim10\,\mathrm{kpc}$. Central and satellite galaxies possess similar amounts and origins of central stars. Star forming galaxies ($\gtrsim10^{10}\,\mathrm{M}_{\odot}$) have on average more ex-situ mass in their centers than quenched ones. We predict readily observable stellar population and dynamical properties: 1) migrated stars are distinctly young ($\sim2\,\mathrm{Gyr}$) and rotationally supported, especially for Milky Way mass galaxies, 2) in-situ stars are most metal-rich and older than migrated stars, 3) ex-situ stars are on random motion dominated orbits and typically the oldest, most metal-poor and $α$-enhanced population. We demonstrate that the interaction history with other galaxies leads to diverse pathways of building up galaxy centers in a $Λ$CDM universe. Our work highlights the necessity for cosmological context in formation scenarios of central galactic components and the potential to use galaxy centers as tracers of overall galaxy assembly. △ Less

Submitted 27 January, 2023; originally announced January 2023.

Comments: 24 pages, 13 Figures, published in MNRAS

iMaNGA: mock MaNGA galaxies based on IllustrisTNG and MaStar SSPs -- II. The catalogue

Authors: Lorenza Nanni, Daniel Thomas, James Trayford, Claudia Maraston, Justus Neumann, David R. Law, Lewis Hill, Annalisa Pillepich, Renbin Yan, Yanping Chen, Dan Lazarz

Abstract: Strengthening the synergy between simulations and observations is essential to test galaxy formation and evolution theories. To achieve this goal, in the first paper of this series, we presented a method to generate mock SDSS-IV/MaNGA integral-field spectroscopic galaxy observations from cosmological simulations. In this second paper, we build the iMaNGA catalogue consisting of $\sim$1,000 unique… ▽ More Strengthening the synergy between simulations and observations is essential to test galaxy formation and evolution theories. To achieve this goal, in the first paper of this series, we presented a method to generate mock SDSS-IV/MaNGA integral-field spectroscopic galaxy observations from cosmological simulations. In this second paper, we build the iMaNGA catalogue consisting of $\sim$1,000 unique galaxies from the TNG50 cosmological simulations, selected to mimic the SDSS-IV/MaNGA-Primary sample selection. Here we present and discuss the iMaNGA sample and its comparison to the MaNGA Primary catalogue. The iMaNGA sample well recovers the MaNGA-Primary sample in terms of stellar mass vs angular size relation and spatial resolution. The Sérsic index vs angular size relation, instead, is not reproduced well by the simulations, mostly because of a paucity of high-mass elliptical galaxies in TNG50. We also investigate our ability to recover the galaxy kinematics and stellar population properties with full-spectral fitting. We demonstrate that 'intrinsic' and `recovered' stellar kinematics, stellar ages and metallicities are consistent, with residuals compatible with zero within 1$-σ$. Also`intrinsic' and `recovered' star formation histories display a great resemblance. We conclude that our mock generation and spectral fitting processes do not distort the `intrinsic' galaxy properties. Therefore, in the third paper of this series, we can meaningfully test the cosmological simulations, comparing the stellar population properties and kinematics of the iMaNGA mock galaxies and the MaNGA observational results. △ Less

Submitted 28 July, 2023; v1 submitted 23 November, 2022; originally announced November 2022.

Comments: Published on MNRAS: 03 May 2023 https://doi.org/10.1093/mnras/stad1337

MaNGIA: 10,000 mock galaxies for stellar population analysis

Authors: Regina Sarmiento, Marc Huertas-Company, Johan H. Knapen, Héctor Ibarra-Medel, Annalisa Pillepich, Sebastián F. Sánchez, Alina Boecker

Abstract: Modern astronomical observations give unprecedented access to the physical properties of nearby galaxies, including spatially resolved stellar populations. However, observations can only give a present-day view of the Universe, whereas cosmological simulations give access to the past record of the processes that galaxies have experienced in their evolution. To connect the events that happened in t… ▽ More Modern astronomical observations give unprecedented access to the physical properties of nearby galaxies, including spatially resolved stellar populations. However, observations can only give a present-day view of the Universe, whereas cosmological simulations give access to the past record of the processes that galaxies have experienced in their evolution. To connect the events that happened in the past with galactic properties as seen today, simulations must be taken to a common ground before being compared to observations. We emulate data from the MaNGA survey, which is the largest integral field spectroscopic galaxy survey to date with its 10,000 nearby galaxies of all types. For this, we use the cosmological simulations TNG50 to generate MaNGIA (Mapping Nearby Galaxies with IllustrisTNG Astrophysics), a mock MaNGA sample of similar size that emulates observations of galaxies for stellar population analysis. We choose TNG galaxies to match the MaNGA sample selection to limit the impact of selection effects. We produce MaNGA-like datacubes from all simulated galaxies, and process these with the pyPipe3D analysis code. This allows us to extract spatially resolved stellar maps. This first paper presents the approach to generate the mock sample and provides an initial exploration of its properties. We show that the stellar populations and kinematics of the simulated MaNGIA galaxies are overall in good agreement with observations. Specific discrepancies, especially in the age and metallicity gradients in low- to intermediate-mass regimes and in massive galaxies' kinematics, require further investigation. We compare our results to other attempts to mock similar observations, all of smaller data sets. Our final dataset will be released with the publication, consisting of >10,000 post-processed data-cubes analysed with pyPipe3D, along with the codes developed to create it. △ Less

Submitted 21 November, 2022; originally announced November 2022.

Comments: Submitted to A&A, 21 pages, 17 figures, 1 table

Line Emission Mapper (LEM): Probing the physics of cosmic ecosystems

Authors: Ralph Kraft, Maxim Markevitch, Caroline Kilbourne, Joseph S. Adams, Hiroki Akamatsu, Mohammadreza Ayromlou, Simon R. Bandler, Marco Barbera, Douglas A. Bennett, Anil Bhardwaj, Veronica Biffi, Dennis Bodewits, Akos Bogdan, Massimiliano Bonamente, Stefano Borgani, Graziella Branduardi-Raymont, Joel N. Bregman, Joseph N. Burchett, Jenna Cann, Jenny Carter, Priyanka Chakraborty, Eugene Churazov, Robert A. Crain, Renata Cumbee, Romeel Dave , et al. (85 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 12 April, 2023; v1 submitted 17 November, 2022; originally announced November 2022.

Comments: 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

Feedback reshapes the baryon distribution within haloes, in halo outskirts, and beyond: the closure radius from dwarfs to massive clusters

Authors: Mohammadreza Ayromlou, Dylan Nelson, Annalisa Pillepich

Abstract: We explore three sets of cosmological hydrodynamical simulations, IllustrisTNG, EAGLE, and SIMBA, to investigate the physical processes impacting the distribution of baryons in and around haloes across an unprecedented mass range of $10^8<M_{\rm 200c}/{\rm M_{\odot}}<10^{15}$, from the halo centre out to scales as large as $30\,R_{\rm 200c}$. We demonstrate that baryonic feedback mechanisms signif… ▽ More We explore three sets of cosmological hydrodynamical simulations, IllustrisTNG, EAGLE, and SIMBA, to investigate the physical processes impacting the distribution of baryons in and around haloes across an unprecedented mass range of $10^8<M_{\rm 200c}/{\rm M_{\odot}}<10^{15}$, from the halo centre out to scales as large as $30\,R_{\rm 200c}$. We demonstrate that baryonic feedback mechanisms significantly redistribute gas, lowering the baryon fractions inside haloes while simultaneously accumulating this material outside the virial radius. To understand this large-scale baryonic redistribution and identify the dominant physical processes responsible, we examine several variants of TNG that selectively exclude stellar and AGN feedback, cooling, and radiation. We find that heating from the UV background in low-mass haloes, stellar feedback in intermediate-mass haloes, and AGN feedback in groups ($10^{12} \leq M_{\rm 200c}/{\rm M_{\odot}}<10^{14}$) are the dominant processes. Galaxy clusters are the least influenced by these processes on large scales. We introduce a new halo mass-dependent characteristic scale, the closure radius $R_{\rm c}$, within which all baryons associated with haloes are found. For groups and clusters, we introduce a universal relation between this scale and the halo baryon fraction: $R_{\rm c}/R_{\rm 200c,500c}-1=β(z)(1-f_{\rm b}(<R_{\rm 200c,500c})/f_{\rm b,cosmic})$, where $β(z)=α\,(1+z)^γ$, and $α$ and $γ$ are free parameters fit using the simulations. Accordingly, we predict that all baryons associated with observed X-ray haloes can be found within $R_{\rm c}\sim 1.5-2.5 R_{\rm 200c}$. Our results can be used to constrain theoretical models, particularly the physics of supernova and AGN feedback, as well as their interplay with environmental processes, through comparison with current and future X-ray and SZ observations. △ Less

Submitted 14 November, 2022; originally announced November 2022.

Comments: Submitted to MNRAS

The survival of stellar discs in Fornax-like environments, from TNG50 to real galaxies

Authors: Pablo M. Galán-de Anta, M. Sarzi, A. Pillepich, Y. Ding, L. Zhu, L. Coccato, E. M. Corsini, K. Fahrion, J. Falcón-Barroso, D. A. Gadotti, E. Iodice, M. Lyubenova, I. Martín-Navarro, R. M. McDermid, F. Pinna, G. van de Ven, P. T. de Zeeuw

Abstract: We study the evolution of kinematically-defined stellar discs in 10 Fornax-like clusters identified in the TNG50 run from the IllustrisTNG suite of cosmological simulations. We considered disc galaxies with present-day stellar mass $M_{\star}\geq 3 \times 10^{8} M_{\odot}$ and follow their evolution since first entering their host cluster. Very few stellar discs survive since falling in such dense… ▽ More We study the evolution of kinematically-defined stellar discs in 10 Fornax-like clusters identified in the TNG50 run from the IllustrisTNG suite of cosmological simulations. We considered disc galaxies with present-day stellar mass $M_{\star}\geq 3 \times 10^{8} M_{\odot}$ and follow their evolution since first entering their host cluster. Very few stellar discs survive since falling in such dense environments, ranging from 40% surviving to all being disrupted. Such survival rates are consistent with what reported earlier for the two more massive, Virgo-like clusters in TNG50. In absolute terms, however, the low number of present-day disc galaxies in Fornax-like clusters could be at odds with the presence of three edge-on disc galaxies in the central regions of the actual Fornax cluster, as delineated by the Fornax3D survey. When looking at the Fornax analogues from random directions and with the same selection function of Fornax3D, the probability of finding three edge-on disc galaxies in any one Fornax-like cluster in TNG50 is rather low, albeit not impossible. We also compared the stellar-population properties near the equatorial plane derived from integral-field spectroscopy for the three edge-ons in Fornax to similar line-of-sight integrated values for present-day disc galaxies in TNG50. For one of these, the very old and metal-rich stellar population of its disc cannot be matched by any the disc galaxies in TNG50, including objects in the field. We discuss possible interpretations of these findings, while pointing to future studies on passive cluster spirals as a way to further test state-of-the-art cosmological simulations. △ Less

Submitted 2 November, 2022; v1 submitted 1 November, 2022; originally announced November 2022.

Comments: 12 pages, 7 figures. Accepted on 21 of October, published on 25 of October

Where are the extremely metal-poor stars in the Milky Way and Andromeda? Expectations from TNG50

Authors: Li-Hsin Chen, Annalisa Pillepich, Simon C. O. Glover, Ralf S. Klessen

Abstract: We analyse the location of extremely metal-poor stars (EMPs, [Fe/H]$ < -3$) in 198 Milky Way (MW)/M31-like galaxies at $z=0$ in the TNG50 simulation. Each system is divided into four kinematically-defined morphological stellar components based on stellar circularity and galactocentric distance, namely bulge, cold disk, warm disk, and stellar halo, in addition to satellites (with stellar mass… ▽ More We analyse the location of extremely metal-poor stars (EMPs, [Fe/H]$ < -3$) in 198 Milky Way (MW)/M31-like galaxies at $z=0$ in the TNG50 simulation. Each system is divided into four kinematically-defined morphological stellar components based on stellar circularity and galactocentric distance, namely bulge, cold disk, warm disk, and stellar halo, in addition to satellites (with stellar mass $\ge 5\times10^6\,M_\odot$). According to TNG50 and across all simulated systems, the stellar halo of the main galaxy and satellites present the highest frequency of EMPs (largest $M_{\mathrm{EMP, comp}}$-to-$M_{\mathrm{tot, comp}}$ stellar mass ratio), and thus the highest chances of finding them. Such frequency is larger in lower-mass than high-mass satellites. Moreover, TNG50 predicts that the stellar halo of the main galaxy always hosts and thus contributes the majority of the EMPs of the system. Namely, it has the highest mass ratio of EMPs in it to all the EMPs in the system (largest $M_{\mathrm{EMP, comp}}$-to-$M_\mathrm{EMP} (<300\mathrm{kpc})$). However, notably, we also find that 33 MW/M31-like galaxies in TNG50 have cold disks that contribute more than 10 per cent to the total EMP mass, each with $\gtrsim 10^{6.5-7}\, M_\odot$ of EMPs in cold circular orbits. These qualitative statements do not depend on the precise definition of EMP stars, i.e. on the adopted metallicity threshold. The results of this work provide a theoretical prediction for the location of EMP stars from both a spatial and kinematic perspective and across an unprecedented number of well-resolved MW/M31-like systems. △ Less

Submitted 1 December, 2022; v1 submitted 31 October, 2022; originally announced November 2022.

Comments: Part of a set of papers based on TNG50 MW/M31-like galaxies. 14 pages, 9 figures. Accepted manuscript by MNRAS

The merger and assembly histories of Milky Way- and M31-like galaxies with TNG50: disk survival through mergers

Authors: Diego Sotillo-Ramos, Annalisa Pillepich, Martina Donnari, Dylan Nelson, Lukas Eisert, Vicente Rodriguez-Gomez, Gandhali Joshi, Mark Vogelsberger, Lars Hernquist

Abstract: We analyze the merger and assembly histories of Milky Way (MW) and Andromeda (M31)-like galaxies to quantify how, and how often, disk galaxies of this mass can survive recent major mergers (stellar mass ratio $\ge$ 1:4). For this, we use the cosmological magneto-hydrodynamical simulation TNG50 and identify 198 analog galaxies, selected based on their $z=0$ stellar mass (… ▽ More We analyze the merger and assembly histories of Milky Way (MW) and Andromeda (M31)-like galaxies to quantify how, and how often, disk galaxies of this mass can survive recent major mergers (stellar mass ratio $\ge$ 1:4). For this, we use the cosmological magneto-hydrodynamical simulation TNG50 and identify 198 analog galaxies, selected based on their $z=0$ stellar mass ($10^{10.5-11.2} {\rm M_{\odot}}$), disky stellar morphology and local environment. Firstly, major mergers are common: 85 per cent (168) of MW/M31-like galaxies in TNG50 have undergone at least one major merger across their lifetime. In fact, 31 galaxies (16 per cent) have undergone a recent major merger, i.e. in the last 5 Gyr. The gas available during the merger suffices to either induce starbursts at pericentric passages or to sustain prolonged star formation after coalescence: in roughly half of the cases, the pre-existing stellar disk is destroyed because of the merger but reforms thanks to star formation. Moreover, higher merger mass ratios are more likely to destroy the stellar disks. In comparison to those with more ancient massive mergers, MW/M31-like galaxies with recent major mergers have, on average, somewhat thicker stellar disks, more massive and somewhat shallower stellar haloes, larger stellar ex-situ mass fractions, but similarly massive kinematically-defined bulges. All this is qualitatively consistent with the different observed properties of the Galaxy and Andromeda and with the constraints on their most recent major mergers, 8-11 and ~2 Gyr ago, respectively. According to contemporary cosmological simulations, a recent quiet merger history is not a pre-requisite for obtaining a relatively-thin stellar disk at $z=0$. △ Less

Submitted 31 October, 2022; originally announced November 2022.

Comments: Published by MNRAS. Check also Chen et al., Engler et al. and Ramesh et al., on arxiv.org today, for further analyses of Milky Way- and Andromeda-like galaxies using the cosmological simulation TNG50. Main figures: 2, 7, 16

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 516, Issue 4, November 2022, Pages 5404 to 5427

The Circumgalactic Medium of Milky Way-like Galaxies in the TNG50 Simulation -- I: Halo Gas Properties and the Role of SMBH Feedback

Authors: Rahul Ramesh, Dylan Nelson, Annalisa Pillepich

Abstract: We analyze the physical properties of gas in the circumgalactic medium (CGM) of 132 Milky Way (MW)-like central galaxies at $z=0$ from the cosmological magneto-hydrodynamical simulation TNG50, part of the IllustrisTNG project. The properties and abundance of CGM gas across the sample are diverse, and the fractional budgets of different phases (cold, warm, and hot), as well as neutral HI mass and m… ▽ More We analyze the physical properties of gas in the circumgalactic medium (CGM) of 132 Milky Way (MW)-like central galaxies at $z=0$ from the cosmological magneto-hydrodynamical simulation TNG50, part of the IllustrisTNG project. The properties and abundance of CGM gas across the sample are diverse, and the fractional budgets of different phases (cold, warm, and hot), as well as neutral HI mass and metal mass, vary considerably. Over our stellar mass range of $10^{10.5} < M_\star / \rm{M}_\odot < 10^{10.9}$, radial profiles of gas physical properties from $0.15 < R\rm{ / R_{\rm 200c}} < 1.0$ reveal great CGM structural complexity, with significant variations both at fixed distance around individual galaxies, and across different galaxies. CGM gas is multi-phase: the distributions of density, temperature and entropy are all multimodal, while metallicity and thermal pressure distributions are unimodal; all are broad. We present predictions for magnetic fields in MW-like halos: a median field strength of $|B|\sim\,1μ$G in the inner halo decreases rapidly at larger distance, while magnetic pressure dominates over thermal pressure only within $\sim0.2 \times \rm{R_{200c}}$. Virial temperature gas at $\sim 10^6\,$K coexists with a sub-dominant cool, $< 10^5\,$K component in approximate pressure equilibrium. Finally, the physical properties of the CGM are tightly connected to the galactic star formation rate, in turn dependent on feedback from supermassive black holes (SMBHs). In TNG50, we find that energy from SMBH-driven kinetic winds generates high-velocity outflows ($\gtrsim 500-2000$ km/s), heats gas to super-virial temperatures ($> 10^{6.5-7}$ K), and regulates the net balance of inflows versus outflows in otherwise quasi-static gaseous halos. △ Less

Submitted 15 January, 2023; v1 submitted 31 October, 2022; originally announced November 2022.

Comments: Minor typos fixed. Additional visuals at www.tng-project.org/ramesh22!

Satellites of Milky Way- and M31-like galaxies with TNG50: quenched fractions, gas content, and star formation histories

Authors: Christoph Engler, Annalisa Pillepich, Gandhali D. Joshi, Anna Pasquali, Dylan Nelson, Eva K. Grebel

Abstract: We analyse the quenched fractions, gas content, and star formation histories of ~1200 satellite galaxies with $M_* \geq 5 \times 10^6~{\rm M}_\odot$ around 198 Milky Way- (MW) and Andromeda-like (M31) hosts in TNG50, the highest-resolution simulation of IllustrisTNG. Satellite quenched fractions are larger for smaller masses, for smaller distances to their host galaxy, and in the more massive M31-… ▽ More We analyse the quenched fractions, gas content, and star formation histories of ~1200 satellite galaxies with $M_* \geq 5 \times 10^6~{\rm M}_\odot$ around 198 Milky Way- (MW) and Andromeda-like (M31) hosts in TNG50, the highest-resolution simulation of IllustrisTNG. Satellite quenched fractions are larger for smaller masses, for smaller distances to their host galaxy, and in the more massive M31-like compared to MW-like hosts. As satellites cross their host's virial radius, their gas content drops: most satellites within 300 kpc lack detectable gas reservoirs at $z=0$, unless they are massive like the Magellanic Clouds and M32. Nevertheless, their stellar assembly exhibits a large degree of diversity. On average, the cumulative star formation histories are more extended for brighter, more massive satellites with a later infall, and for those in less massive hosts. Based on these relationships, we can even infer infall periods for observedMWand M31 dwarfs: e.g. 0-4 Gyr ago for the Magellanic Clouds and Leo I, 4-8 and 0-2 Gyr ago for M32 and IC 10, respectively. Ram pressure stripping (in combination with tidal stripping) deprives TNG50 satellites of their gas reservoirs and ultimately quenches their star formation, even though only a few per cent of the present-day satellites around the 198 TNG50 MW/M31-like hosts appear as jellyfish. The typical time since quenching for currently quenched TNG50 satellites is $6.9^{+2.5}_{-3.3}~{\rm Gyr}$ ago. The TNG50 results are consistent with the quenched fractions and stellar assembly of observed MW and M31 satellites, however, satellites of the SAGA survey with $M_* \sim 10^{8-9}~{\rm M}_\odot$ exhibit lower quenched fractions than TNG50 and other, observed analogues. △ Less

Submitted 17 May, 2023; v1 submitted 31 October, 2022; originally announced November 2022.

Comments: Accepted for publication in MNRAS

Modeling globular clusters in the TNG50 simulation: predictions from dwarfs to giants

Authors: Jessica E. Doppel, Laura V. Sales, Dylan Nelson, Annalisa Pillepich, Mario G. Abadi, Eric W. Peng, Federico Marinacci, Jill Naiman, Paul Torrey, Mark Vogelsberger, Rainer Weinberger, Lars Hernquist

Abstract: We present a post-processing catalog of globular clusters (GCs) for the $39$ most massive groups and clusters in the TNG50 simulation of the IlllustrisTNG project (virial masses $M_{200} =[5\times 10^{12} \rm - 2 \times 10^{14}$] M$_{\odot}$). We tag GC particles to all galaxies with stellar mass $M_* \geq 5\times10^6$ M$_{\odot}$, and we calibrate their masses to reproduce the observed power-law… ▽ More We present a post-processing catalog of globular clusters (GCs) for the $39$ most massive groups and clusters in the TNG50 simulation of the IlllustrisTNG project (virial masses $M_{200} =[5\times 10^{12} \rm - 2 \times 10^{14}$] M$_{\odot}$). We tag GC particles to all galaxies with stellar mass $M_* \geq 5\times10^6$ M$_{\odot}$, and we calibrate their masses to reproduce the observed power-law relation between GC mass and halo mass for galaxies with $M_{200} \geq 10^{11}$ M$_{\odot}$ (corresponding to $M_* \sim 10^9$ $M_{\odot}$). Here we explore whether an extrapolation of this $M_{\rm GC}$-$M_{200}$ relation to lower-mass dwarfs is consistent with current observations. We find a good agreement between our predicted number and specific frequency of GCs in dwarfs with $\rm M_*=[5 \times 10^6 \rm - 10^9]$ M$_{\odot}$ and observations. Moreover, we predict a steep decline in the GC occupation fraction for dwarfs with $M_*<10^9$ M$_{\odot}$ which agrees well with current observational constraints. This declining occupation fraction is due to a combination of tidal stripping in all dwarfs plus a stochastic sampling of the GC mass function for dwarfs with $M_* < 10^{7.5}$ M$_{\odot}$. Our simulations also reproduce available constraints on the abundance of intra-cluster GCs in Virgo and Centaurus A. These successes provide support to the hypothesis that the $M_{\rm GC}$-$M_{200}$ relation holds, albeit with more scatter, all the way down to the regime of classical dwarf spheroidals in these environments. Our GC catalogs are publicly available as part of the IllustrisTNG data release. △ Less

Submitted 23 September, 2022; originally announced September 2022.

Comments: 16 pages, 8 figures (main body of paper). Resubmitted to MNRAS after positive referee report. The globular cluster catalogs for the 39 most massive groups and clusters in TNG50 are available for download from the IllustrisTNG public data release at www.tng-project.org/doppel22