Search

We introduce the new cosmological simulation project cosmosTNG, a first-of-its-kind suite of constrained galaxy formation simulations for the universe at Cosmic Noon ($z\sim 2$). cosmosTNG simulates a $0.2$ deg$^2$ patch of the COSMOS field at $z \simeq 2.0-2.2$ using an initial density field inferred from galaxy redshift surveys and the CLAMATO Lyman-alpha forest tomography survey, reconstructed by the TARDIS algorithm. We evolve eight different realizations of this volume to capture small-scale variations. All runs use the IllustrisTNG galaxy formation model with a baryonic mass resolution of $10^6$ M$_\odot$, equal to TNG100-1. In this initial study, we demonstrate qualitative agreement between the evolved large-scale structure and the spatial distribution of observed galaxy populations in COSMOS, emphasizing the zFIRE protocluster region. We then compare the statistical properties and scaling relations of the galaxy population, covering stellar, gaseous, and supermassive black hole (SMBH) components, between cosmosTNG, observations in COSMOS, and $z \sim 2$ observational data in general. We find that galaxy quenching and environmental effects in COSMOS are modulated by its specific large-scale structure, particularly the collapsing protoclusters in the region. With respect to a random region of the universe, the abundance of high-mass galaxies is larger, and the quenched fraction of galaxies is significant lower at fixed mass. This suggests an accelerated growth of stellar mass, as reflected in a higher cosmic star formation rate density, due to the unique assembly histories of galaxies in the simulated COSMOS subvolume. The cosmosTNG suite will be a valuable tool for studying galaxy formation at cosmic noon, particularly when interpreting extragalactic observations with HST, JWST, and other large multi-wavelength survey programs of the COSMOS field.

According to the well established hierarchical framework for galaxy evolution, galaxies grow through mergers with other galaxies and the LambdaCDM cosmological model predicts that the stellar halos of galaxies are rich in remnants from minor mergers. The Stellar Streams Legacy Survey has provided a first release of a catalogue with a statistically significant sample of stellar streams in the Local Universe that can be used to study minor mergers and test the cosmological models. The main objective is to compare the results of the observations of stellar tidal streams with the predictions of state-of-the-art cosmological simulations regarding the formation of stellar streams up to a redshift z < 0.02, according to the LambdaCDM model. We use the predictions of the cosmological simulations Copernicus Complexio, TNG50 of the IllustrisTNG project and Auriga to generate 225 mock-images of nearby halos at a distance of 70 Mpc, and search for stellar streams. We compare the obtained stream frequency and characteristics with those obtained from the Stellar Streams Legacy Survey. We find good agreement between the results of analysing real images from the Dark Energy Survey and mock-images from cosmological simulations. We obtained predictions for the detection rate of stellar streams to a surface brightness limit of 35 mag arcsec^-2. The cosmological simulations predict that for a surface brightness limit of 32 mag arcsec^-2 a frequency of almost 70% in the detection of streams around galaxies can be achieved., Comment: 18 pages, 17 figures

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., Comment: Submitted to MNRAS. 22 pages. Comments are welcome

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., Comment: 16 pages, 7 Figures. Accepted for publication

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., Comment: 11 pages, 7 figures, submitted to ApJ

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\v{c}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., Comment: Main text 14 pages, 5 figures, accepted to MNRAS. Generated data products are publicly available at www.tng-project.org/gebek24

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.03$, 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., Comment: Accepted for publication in ApJ

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., Comment: Accepted to A&A, 16 pages, 12 figures

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

Understanding how the properties of galaxies relate to the properties of the hot circum-galactic medium (CGM) around them can constrain galaxy evolution models. We measured the X-ray luminosity of the hot CGM based on the surface brightness profiles of central galaxy samples measured from Spectrum Roentgen Gamma (SRG)/eROSITA all-sky survey data. We related the X-ray luminosity to the galaxies' stellar and halo mass, and we compared the observed relations to the self-similar model and intrinsic (i.e., not forward-modeled) output of the IllustrisTNG, EAGLE, and SIMBA simulations. The average hot CGM X-ray luminosity ($L_{\rm X,CGM}$) correlates with the galaxy's stellar mass ($M_*$). It increases from $(1.6 \pm 2.1)\times10^{39} \rm erg\,s^{-1}$ to $(3.4 \pm 0.3)\times10^{41} \rm erg\,s^{-1}$, when $\log(M_*)$ increases from 10.0 to 11.5. A power law describes the correlation as $\log(L_{\rm X,CGM})= (2.4\pm 0.1)\times \log(M_*)+(14.6\pm1.5)$. 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 $(3.0 \pm 1.6)\times10^{39}\ \rm erg\,s^{-1}$ at $\log(M_{\rm 500c})=11.3$ to $(1.3 \pm 0.1)\times10^{42}\ \rm erg\,s^{-1}$ at $\log(M_{\rm 500c})=13.7$. The relation follows a power law of $\log(L_{\rm X,CGM})= (1.32\pm 0.05)\times \log(M_{\rm 500c})+(24.1\pm0.7)$. Our observations highlight the necessity of non-gravitational processes at the galaxy group scale while suggesting these processes are sub-dominant at the galaxy scale. We show that the outputs of current cosmological galaxy simulations generally align with the observational results uncovered here but with possibly important deviations in selected mass ranges., Comment: Abstract abridged for arXiv submission, accepted by A&A

The circumgalactic medium (CGM) provides the material needed for galaxy formation and influences galaxy evolution. The hot ($T>10^6K$) CGM is poorly detected around galaxies with stellar masses ($M_*$) lower than $3\times10^{11}M_\odot$ due to the low surface brightness. We used the X-ray data from the first four SRG/eROSITA All-Sky Surveys (eRASS:4). Based on the SDSS spectroscopic survey and halo-based group finder algorithm, we selected central galaxies with spectroscopic redshifts of $z_{\rm spec}<0.2$ and stellar masses of $10.0<\log(M_*/M_\odot)<11.5$ (85,222 galaxies) -- or halo masses of $11.5<\log(M_{\rm 200m}/M_\odot)<14.0$ (125,512 galaxies). By stacking the X-ray emission around galaxies, masking the detected X-ray point sources and carefully modeling the X-ray emission from the unresolved active galactic nuclei (AGN) and X-ray binaries (XRB), we obtain the X-ray emission from the hot CGM. We detected the X-ray emission around MW-mass and more massive central galaxies extending up to the virial radius ($R_{\rm vir}$). We used a $\beta$ model to describe the X-ray surface brightness profile and found $\beta =0.43^{+0.10}_{-0.06}\,(0.37^{+0.04}_{-0.02})$ for MW-mass (M31-mass) galaxies.We estimated the baryon budget of the hot CGM and obtained a value that is lower than the prediction of $\Lambda$CDM cosmology, indicating significant gas depletion in these halos. We extrapolated the hot CGM profile measured within $R_{\rm vir}$ to larger radii and found that within $\approx 3 R_{\rm vir}$, the baryon budget is close to the $\Lambda$CDM cosmology prediction. Our results set a firm footing for the presence of the hot CGM around such galaxies. These measurements constitute a new benchmark for galaxy evolution models and possible implementations of feedback processes therein., Comment: Abstract abridged for arXiv submission, accepted by A&A

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., Comment: 24 pages, 16 figures, 2 appendices, Accepted for publication in A&A. See the TNG-Cluster website at www.tng-project.org/cluster/

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., Comment: 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

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., Comment: Submitted to A&A. See the TNG-Cluster website at https://www.tng-project.org/cluster/

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., Comment: Submitted to A&A, comments welcome. See the TNG-Cluster website at www.tng-project.org/cluster/ for more details

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., Comment: 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 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., Comment: 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/

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., Comment: 26 pages, 17 figures. Comments are welcome. Highlights in Fig. 3, 7 and 13

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), Comment: 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

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., Comment: 9 pages, 4 figures, accepted in ApJL

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., Comment: 19 pages. White paper for a mission concept to be submitted for the 2023 NASA Astrophysics Probes opportunity

At fixed galaxy stellar mass, there is a clear observational connection between structural asymmetry and offset from the star forming main sequence, $\Delta$SFMS. Herein, we use the TNG50 simulation to investigate the relative roles of major mergers (stellar mass ratios $\mu\geq0.25$), minor ($0.1 \leq \mu < 0.25$), and mini mergers ($0.01 \leq \mu < 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 $\Delta$SFMS and asymmetry and (2) a strikingly similar trend emerges between $\Delta$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 $\mu\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., Comment: 31 pages; 19 figures; MNRAS accepted; Image data available via the TNG website: www.tng-project.org/bottrell23

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., Comment: Accepted by MNRAS. 4 figures

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., Comment: 22 pages, 16 figures. Accepted for publication in MNRAS

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., Comment: 46 pages, 32 figures, accepted to ApJ

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

The 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., Comment: 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

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., Comment: Published in MNRAS. See https://www.lem-observatory.org/ and https://www.tng-project.org/ for more details; 2023MNRAS.522.3665N

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., Comment: 23 pages, 15 figures. Accepted for publication in MNRAS

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^3610^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., Comment: Accepted for publication on A&A

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., Comment: submitted to MNRAS (12/05/23); Accepted 07/12/23 after minor revision

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)., Comment: 6 pages, 4 figures

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., Comment: 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/

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., Comment: 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 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., Comment: 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

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., Comment: 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

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 $\Lambda$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., Comment: 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+

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 ($\chi = 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., Comment: 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

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., Comment: 19 pages, 11 figures. Accepted for publication in MNRAS

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., Comment: 28 pages, 22 figures. Accepted for publication in ApJ

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 $\alpha$-enhanced population. We demonstrate that the interaction history with other galaxies leads to diverse pathways of building up galaxy centers in a $\Lambda$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., Comment: 24 pages, 13 Figures, published in MNRAS

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\'ersic 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$-{\sigma}$. 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., Comment: Published on MNRAS: 03 May 2023 https://doi.org/10.1093/mnras/stad1337

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., Comment: Submitted to A&A, 21 pages, 17 figures, 1 table

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

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

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., Comment: 12 pages, 7 figures. Accepted on 21 of October, published on 25 of October

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., Comment: Part of a set of papers based on TNG50 MW/M31-like galaxies. 14 pages, 9 figures. Accepted manuscript by MNRAS

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$., Comment: 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

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\mu$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., Comment: Minor typos fixed. Additional visuals at www.tng-project.org/ramesh22!