Researcher profile

Lars Hernquist

Lars Hernquist contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
66works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

66 published item(s)

preprint2025arXiv

The THESAN-ZOOM project: Star-formation efficiencies in high-redshift galaxies

Recent JWST observations hint at unexpectedly intense cosmic star-formation in the early Universe, often attributed to enhanced star-formation efficiencies (SFEs). Here, we analyze the SFE in THESAN-ZOOM, a novel zoom-in radiation-hydrodynamic simulation campaign of high-redshift ($z \gtrsim 3$) galaxies employing a state-of-the-art galaxy formation model resolving the multiphase interstellar medium (ISM). The halo-scale SFE ($ε^{\ast}_{\rm halo}$) - the fraction of baryons accreted by a halo that are converted to stars - follows a double power-law dependence on halo mass, with a mild redshift evolution above $M_{\rm halo} \gtrsim 10^{9.5}\,{\rm M}_{\odot}$. The power-law slope is roughly $1/3$ at large halo masses, consistent with expectations when gas outflows are momentum-driven. At lower masses, the slope is roughly $2/3$ and is more aligned with the energy-driven outflow scenario. $ε^{\ast}_{\rm halo}$ is a factor of $2-3$ larger than commonly assumed in empirical galaxy-formation models at $M_{\rm halo} \lesssim 10^{11}\,{\rm M}_{\odot}$. On galactic (kpc) scales, the Kennicutt-Schmidt (KS) relation of neutral gas is universal in THESAN-ZOOM, following $Σ_{\rm SFR} \propto Σ_{\rm gas}^2$, indicative of a turbulent energy balance in the ISM maintained by stellar feedback. The rise of $ε^{\ast}_{\rm halo}$ with halo mass can be traced primarily to increasing gas surface densities in massive galaxies, while the underlying KS relation and neutral, star-forming gas fraction remain unchanged. Although the increase in $ε^{\ast}_{\rm halo}$ with redshift is relatively modest, it is sufficient to explain the large observed number density of UV-bright galaxies at $z \gtrsim 12$. However, reproducing the brightest sources at $M_{\rm UV} \lesssim -21$ may require extrapolating the SFE beyond the halo mass range directly covered by THESAN-ZOOM.

preprint2024arXiv

Gas Morphology of Milky Way-like Galaxies in the TNG50 Simulation: Signals of Twisting and Stretching

We present an in-depth analysis of gas morphologies for a sample of 25 Milky Way-like galaxies from the IllustrisTNG TNG50 simulation. We constrain the morphology of cold, warm, hot gas, and gas particles as a whole using a Local Shell Iterative Method (LSIM) and explore its observational implications by computing the hard-to-soft X-ray ratio, which ranges between $10^{-3}$-$10^{-2}$ in the inner $\sim 50 \rm kpc$ of the distribution and $10^{-5}$-$10^{-4}$ at the outer portion of the hot gas distribution. We group galaxies into three main categories: simple, stretched, and twisted. These categories are based on the radial reorientation of the principal axes of the reduced inertia tensor. We find that a vast majority ($77\%$) of the galaxies in our sample exhibit twisting patterns in their radial profiles. Additionally, we present detailed comparisons between 1) the gaseous distributions belonging to individual temperature regimes, 2) the cold gas distributions and stellar distributions, and 3) the gaseous distributions and dark matter (DM) halos. We find a strong correlation between the morphological properties of the cold gas and stellar distributions. Furthermore, we find a correlation between gaseous distributions with DM halo that increases with gas temperature, implying that we may use the warm-hot gaseous morphology as a tracer to probe the DM morphology. Finally, we show gaseous distributions exhibit significantly more prolate morphologies than the stellar distributions and DM halos, which we hypothesize is due to stellar and AGN feedback.

preprint2022arXiv

Breaking baryon-cosmology degeneracy with the electron density power spectrum

Uncertain feedback processes in galaxies affect the distribution of matter, currently limiting the power of weak lensing surveys. If we can identify cosmological statistics that are robust against these uncertainties, or constrain these effects by other means, then we can enhance the power of current and upcoming observations from weak lensing surveys such as DES, Euclid, the Rubin Observatory, and the Roman Space Telescope. In this work, we investigate the potential of the electron density auto-power spectrum as a robust probe of cosmology and baryonic feedback. We use a suite of (magneto-)hydrodynamic simulations from the CAMELS project and perform an idealized analysis to forecast statistical uncertainties on a limited set of cosmological and physically-motivated astrophysical parameters. We find that the electron number density auto-correlation, measurable through either kinematic Sunyaev-Zel'dovich observations or through Fast Radio Burst dispersion measures, provides tight constraints on $Ω_{m}$ and the mean baryon fraction in intermediate-mass halos, $\bar{f}_{\mathrm{bar}}$. By obtaining an empirical measure for the associated systematic uncertainties, we find these constraints to be largely robust to differences in baryonic feedback models implemented in hydrodynamic simulations. We further discuss the main caveats associated with our analysis, and point out possible directions for future work.

preprint2022arXiv

Dynamics of Intermediate-Mass Black Holes Wandering in the Milky Way Galaxy Using the Illustris TNG50 Simulation

The detection of Intermediate-Mass Black Holes (IMBHs) in dwarf galaxies is crucial to closing the gap in the wide mass distribution of black holes ($\sim 3 \, \rm M_{\odot}$ to $\sim 5 \times 10^{10} \, \rm M_{\odot}$). IMBHs originally located at the center of dwarfs that later collide with the Milky Way (MW) could be wandering, undetected, in our Galaxy. We used TNG50, the highest-resolution run of the IllustrisTNG project, to study the kinematics and dynamics of star clusters, in the appropriate mass range, acting as IMBH proxies in a MW analog galaxy. We showed that $\sim 87\%$ of our studied IMBHs drift inward. The radial velocity of these sinking IMBHs has a median magnitude of $\sim 0.44 \, \rm ckpc \, h^{-1} \, Gyr^{-1}$ and no dependence on the black hole mass. The central $1 \, \rm ckpc \, h^{-1}$ has the highest number density of IMBHs in the galaxy. A physical toy model with linear drag forces was developed to explain the orbital circularization with time. These findings constrain the spatial distribution of IMBHs, suggesting that future searches should focus on the central regions of the Galaxy. Additionally, we found that the 3D velocity distribution of IMBHs with respect to the galactic center has a mean of $\sim 180 \, \rm km \, s^{-1}$ and larger variance with decreasing radius. Remarkably, the velocity distribution relative to the local gas shows significantly lower values, with a mean of $\sim 88 \, \rm km \, s^{-1}$. These results are instrumental for predicting the accretion and radiation properties of IMBHs, facilitating their detection with future surveys.

preprint2022arXiv

Early-type galaxy density profiles from IllustrisTNG: III. Effects on outer kinematic structure

Early-type galaxies (ETGs) possess total density profiles close to isothermal, which can lead to non-Gaussian line-of-sight velocity dispersion (LOSVD) under anisotropic stellar orbits. However, recent observations of local ETGs in the MASSIVE Survey reveal outer kinematic structures at $1.5 R_{\mathrm{eff}}$ (effective radius) that are inconsistent with fixed isothermal density profiles; the authors proposed varying density profiles as an explanation. We aim to verify this conjecture and understand the influence of stellar assembly on these kinematic features through mock ETGs in IllustrisTNG. We create mock Integral-Field-Unit observations to extract projected stellar kinematic features for 207 ETGs with stellar mass $M_{\ast}\geqslant 10^{11} \mathrm{M_{\odot}}$ in TNG100-1. The mock observations reproduce the key outer ($1.5R_{\mathrm{eff}}$) kinematic structures in the MASSIVE ETGs, including the puzzling positive correlation between velocity dispersion profile outer slope $γ_{\mathrm{outer}}$ and the kurtosis $h_{4}$'s gradient. We find that $h_{4}$ is uncorrelated with stellar orbital anisotropy beyond $R_{\mathrm{eff}}$; instead we find that the variations in $γ_{\mathrm{outer}}$ and outer $h_{4}$ (a good proxy for $h_{4}$ gradient) are both driven by variations of the density profile at the outskirts across different ETGs. These findings corroborate the proposed conjecture and rule out velocity anisotropy as the origin of non-Gaussian outer kinematic structure in ETGs. We also find that the outer kurtosis and anisotropy correlate with different stellar assembly components, with the former related to minor mergers or flyby interactions while the latter is mainly driven by major mergers, suggesting distinct stellar assembly origins that decorrelates the two quantities.

preprint2022arXiv

Fast, Slow, Early, Late: Quenching Massive Galaxies at z~0.8

We investigate the stellar populations for a sample of 161 massive, mainly quiescent galaxies at $\langle z_{\rm obs} \rangle=0.8$ with deep Keck/DEIMOS rest-frame optical spectroscopy (HALO7D survey). With the fully Bayesian framework Prospector, we simultaneously fit the spectroscopic and photometric data with an advanced physical model (including non-parametric star-formation histories, emission lines, variable dust attenuation law, and dust and AGN emission) together with an uncertainty and outlier model. We show that both spectroscopy and photometry are needed to break the dust-age-metallicity degeneracy. We find a large diversity of star-formation histories: although the most massive ($M_{\star}>2\times10^{11}~M_{\odot}$) galaxies formed the earliest (formation redshift of $z_{\rm f}\approx5-10$ with a short star-formation timescale of $τ_{\rm SF}\lesssim1~\mathrm{Gyr}$), lower-mass galaxies have a wide range of formation redshifts, leading to only a weak trend of $z_{\rm f}$ with $M_{\star}$. Interestingly, several low-mass galaxies with have formation redshifts of $z_{\rm f}\approx5-8$. Star-forming galaxies evolve about the star-forming main sequence, crossing the ridgeline several times in their past. Quiescent galaxies show a wide range and continuous distribution of quenching timescales ($τ_{\rm quench}\approx0-5~\mathrm{Gyr}$) with a median of $\langleτ_{\rm quench}\rangle=1.0_{-0.9}^{+0.8}~\mathrm{Gyr}$ and of quenching epochs of $z_{\rm quench}\approx0.8-5.0$ ($\langle z_{\rm quench}\rangle=1.3_{-0.4}^{+0.7}$). This large diversity of quenching timescales and epochs points toward a combination of internal and external quenching mechanisms. In our sample, rejuvenation and "late bloomers" are uncommon. In summary, our analysis supports the "grow & quench" framework and is consistent with a wide and continuously-populated diversity of quenching timescales.

preprint2022arXiv

Formation and fate of low metallicity stars in IllustrisTNG50

Low metallicity stars give rise to unique spectacular transients and are of immense interest for understanding stellar evolution. Their importance has only grown further with the recent detections of mergers of stellar mass black holes that likely originate mainly from low metallicity progenitor systems. Moreover, the formation of low metallicity stars is intricately linked to galaxy evolution, in particular to early enrichment and to later accretion and mixing of lower metallicity gas. Because low metallicity stars are difficult to observe directly, cosmological simulations are crucial for understanding their formation. Here we quantify the rates and locations of low metallicity star formation using the high-resolution TNG50 magnetohydrodynamical cosmological simulation, and we examine where low metallicity stars end up at $z=0$. We find that $20\%$ of stars with $Z_*<0.1\,\mathrm{Z_\odot}$ form after $z=2$, and that such stars are still forming in galaxies of all masses at $z=0$ today. Moreover, most low-metallicity stars at $z=0$ reside in massive galaxies. We analyse the radial distribution of low metallicity star formation, and discuss the curious case of seven galaxies in TNG50 that form stars from primordial gas even at $z=0$.

preprint2022arXiv

Galactic angular momentum in the IllustrisTNG simulation -- I. Connection to morphology, halo spin, and black hole mass

We use the TNG100 simulation of the IllustrisTNG project to investigate the stellar specific angular momenta ($j_{\ast}$) of $\sim$12,000 central galaxies at $z=0$ in a full cosmological context, with stellar masses ($M_{\ast}$) ranging from $10^{9}$ to $10^{12} \, {\rm M}_{\odot}$. We find that the $j_{\ast}$-$M_{\ast}$ relations for early-type and late-type galaxies in IllustrisTNG are in good overall agreement with observations, and that these galaxy types typically `retain&#39; $\sim$10-20 and $\sim$50-60 per cent of their host haloes&#39; specific angular momenta, respectively, with some dependence on the methodology used to measure galaxy morphology. We present results for kinematic as well as visual-like morphological measurements of the simulated galaxies. Next, we explore the scatter in the $j_{\ast}$-$M_{\ast}$ relation with respect to the spin of the dark matter halo and the mass of the supermassive black hole (BH) at the galactic centre. We find that galaxies residing in faster spinning haloes, as well as those hosting less massive BHs, tend to have a higher specific angular momentum. We also find that, at fixed galaxy or halo mass, halo spin and BH mass are anticorrelated with each other, probably as a consequence of more efficient gas flow toward the galactic centre in slowly rotating systems. Finally, we show that halo spin plays an important role in determining galaxy sizes - larger discs form at the centres of faster-rotating haloes - although the trend breaks down for massive galaxies with $M_{\ast} \gtrsim 10^{11} \, {\rm M}_{\odot}$, roughly the mass scale at which a galaxy&#39;s stellar mass becomes dominated by accreted stars.

preprint2022arXiv

Galaxy Formation in the Santa Cruz semi-analytic model compared with IllustrisTNG -- I. Galaxy scaling relations, dispersions, and residuals at z=0

We present the first results from applying the Santa Cruz semi-analytic model (SAM) for galaxy formation on merger trees extracted from a dark matter only version of the IllustrisTNG (TNG) simulations. We carry out a statistical comparison between the predictions of the Santa Cruz SAM and TNG for a subset of central galaxy properties at z = 0, with a focus on stellar mass, cold and hot gas mass, star formation rate (SFR), and black hole (BH) mass. We find fairly good agreement between the mean predictions of the two methods for stellar mass functions and the stellar mass vs. halo mass (SMHM) relation, and qualitatively good agreement between the SFR or cold gas mass vs. stellar mass relation and quenched fraction as a function of stellar mass. There are greater differences between the predictions for hot (circumgalactic) gas mass and BH mass as a function of halo mass. Going beyond the mean relations, we also compare the dispersion in the predicted scaling relations, and the correlation in residuals on a halo-by-halo basis between halo mass and galaxy property scaling relations. Intriguingly, we find similar correlations between residuals in SMHM in the SAM and in TNG, suggesting that these relations may be shaped by similar physical processes. Other scaling relations do not show significant correlations in the residuals, indicating that the physics implementations in the SAM and TNG are significantly different.

preprint2022arXiv

Gas-phase metallicity break radii of star-forming galaxies in IllustrisTNG

We present radial gas-phase metallicity profiles, gradients, and break radii at redshift $z = 0 - 3$ from the TNG50-1 star-forming galaxy population. These metallicity profiles are characterized by an emphasis on identifying the steep inner gradient and flat outer gradient. From this, the break radius, $r_{\rm Break}$, is defined as the region where the transition occurs. We observe the break radius having a positive trend with mass that weakens with redshift. When normalized by the stellar half-mass radius, the break radius has a weaker relation with both mass and redshift. To test if our results are dependent on the resolution or adopted physics of TNG50-1, the same analysis is performed in TNG50-2 and Illustris-1. We find general agreement between each of the simulations in their qualitative trends; however, the adopted physics between TNG and Illustris differ and therefore the breaks, normalized by galaxy size, deviate by a factor of $\sim$2. In order to understand where the break comes from, we define two relevant time-scales: an enrichment time-scale and a radial gas mixing time-scale. We find that $r_{\rm Break}$ occurs where the gas mixing time-scale is $\sim$10 times as long as the enrichment time-scale in all three simulation runs, with some weak mass and redshift dependence. This implies that galactic disks can be thought of in two-parts: a star-forming inner disk with a steep gradient and a mixing-dominated outer disk with a flat gradient, with the break radius marking the region of transition between them.

preprint2022arXiv

H-alpha emission in local galaxies: star formation, time variability and the diffuse ionized gas

The nebular recombination line H$α$ is widely used as a star-formation rate (SFR) indicator in the local and high-redshift Universe. We present a detailed H$α$ radiative transfer study of high-resolution isolated Milky-Way and Large Magellanic Cloud simulations that include radiative transfer, non-equilibrium thermochemistry, and dust evolution. We focus on the spatial morphology and temporal variability of the H$α$ emission, and its connection to the underlying gas and star formation properties. The H$α$ and H$β$ radial and vertical surface brightness profiles are in excellent agreement with observations of nearby galaxies. We find that the fraction of H$α$ emission from collisional excitation amounts to $f_{\rm col}\sim5-10\%$, only weakly dependent on radius and vertical height, and that scattering boosts the H$α$ luminosity by $\sim40\%$. The dust correction via the Balmer decrement works well (intrinsic H$α$ emission recoverable within $25\%$), though the dust attenuation law depends on the amount of attenuation itself both on spatially resolved and integrated scales. Important for the understanding of the H$α$-SFR connection is the dust and helium absorption of ionizing radiation (Lyman continuum [LyC] photons), which are about $f_{\rm abs}\approx28\%$ and $f_{\rm He}\approx9\%$, respectively. Together with an escape fraction of $f_{\rm esc}\approx6\%$, this reduces the available budget for hydrogen line emission by nearly half ($f_{\rm H}\approx57\%$). We discuss the impact of the diffuse ionized gas, showing - among other things - that the extraplanar H$α$ emission is powered by LyC photons escaping the disc. Future applications of this framework to cosmological (zoom-in) simulations will assist in the interpretation of spectroscopy of high-redshift galaxies with the upcoming James Webb Space Telescope.

preprint2022arXiv

Mass of the dynamically hot inner stellar halo predicts the ancient accreted stellar mass

Galactic dynamical structures are fossil records of the assembly histories of galaxies. By analyzing the cosmological hydrodynamical simulation TNG50, we find that a dynamical structure that we call the &#34;hot inner stellar halo,&#34; defined by stars on dynamically hot orbits with circularity $λ_z < 0.5$ at $3.5\,{\rm kpc}<r \lesssim 2\,R_e$, is a strong indicator of the mass of accreted satellite galaxies. We find a strong correlation between the mass of this hot inner stellar halo and the total ex situ stellar mass. There is a similarly strong correlation with the stellar mass of the most massive secondary galaxy ever merged. These TNG50 correlations are compatible with those predicted by other simulations, for example by TNG100 across the whole mass range under study (galaxy stellar masses, $M_*$, in the $10^{10.3-11.6}$\,\Msun\, range) and by EAGLE for $M_* \gtrsim 10^{10.6} $\,\Msun\, galaxies.\ This shows that our predictions are robust across different galaxy formation and feedback models and hold across a wide range of numerical resolution. The hot inner stellar halo is a product of massive and typically ancient mergers, with inner-halo stars exhibiting three main physical origins: accreted and stripped from massive satellites, dynamically heated by mergers from the bulge and/or disk in the main progenitor, and formed from star formation triggered during mergers. The mass of the hot inner stellar halo defined in this paper is a quantity that can be robustly obtained for real galaxies by applying a population-orbit superposition method to integral-field-unit spectroscopy data, out to a distance of $\sim2\,R_e$, which is possible with current observations. Hence, this paper shows that integral-field-unit observations and dynamical models of the inner regions of galaxies provide a way to quantitatively determine the mass of ancient accreted satellites.

preprint2022arXiv

Morphological decomposition of TNG50 galaxies: methodology and catalogue

We present MORDOR (MORphological DecOmposeR, a new algorithm for structural decomposition of simulated galaxies based on stellar kinematics. The code measures the properties of up to five structural components (a thin/cold and a thick/warm disc, a classical and a secular bulge, and a spherical stellar halo), and determines the properties of a stellar bar (if present). A comparison with other algorithms presented in the literature yields overall good agreement, with MORDOR displaying a higher flexibility in correctly decomposing systems and identifying bars in crowded environments (e.g. with ongoing fly-bys, often observable in cosmological simulations). We use MORDOR to analyse galaxies in the TNG50 simulation and find the following: ($i$) the thick disc component undergoes the strongest evolution in the binding energy-circularity plane, as expected when disc galaxies decrease their turbulent-rotational support with cosmic time; ($ii$) smaller galaxies (with stellar mass, $10^{9} \lesssim M_{*} / {\rm M_{\odot}} \leq 5 \times 10^{9}$) undergo a major growth in their disc components after $z\sim 1$, whereas ($iii$) the most massive galaxies ($5 \times 10^{10} < M_{*} / {\rm M_{\odot}} \leq 5\times10^{11}$) evolve toward more spheroidal dominated objects down to $z=0$ due to frequent gravitational interactions with satellites; ($iv$) the fraction of barred galaxies grows rapidly at high redshift and stabilizes below $z\sim 2$, except for the most massive galaxies that show a decrease in the bar occupation fraction at low redshift; ($v$) galaxies with $M_{*} \sim 10^{11}~{\rm M_{\odot}}$ exhibit the highest relative occurrence of bars at $z=0$, in agreement with observational studies. We publicly release MORDOR and the morphological catalogue of TNG50 galaxies.

preprint2022arXiv

Percent-level constraints on baryonic feedback with spectral distortion measurements

High-significance measurements of the monopole thermal Sunyaev-Zel&#39;dovich CMB spectral distortions have the potential to tightly constrain poorly understood baryonic feedback processes. The sky-averaged Compton-y distortion and its relativistic correction are measures of the total thermal energy in electrons in the observable universe and their mean temperature. We use the CAMELS suite of hydrodynamic simulations to explore possible constraints on parameters describing the subgrid implementation of feedback from active galactic nuclei and supernovae, assuming a PIXIE-like measurement. The small 25 Mpc/h CAMELS boxes present challenges due to the significant cosmic variance. We utilize machine learning to construct interpolators through the noisy simulation data. Using the halo model, we translate the simulation halo mass functions into correction factors to reduce cosmic variance where required. Our results depend on the subgrid model. In the case of IllustrisTNG, we find that the best-determined parameter combination can be measured to ~2% and corresponds to a product of AGN and SN feedback. In the case of SIMBA, the tightest constraint is ~0.2% on a ratio between AGN and SN feedback. A second orthogonal parameter combination can be measured to ~8%. Our results demonstrate the significant constraining power a measurement of the late-time spectral distortion monopoles would have for baryonic feedback models.

preprint2022arXiv

Probing the $z\gtrsim6$ quasars in a universe with IllustrisTNG physics: Impact of gas-based black hole seeding models

We explore implications of a range of black hole (BH) seeding prescriptions on the formation of the brightest $z\gtrsim6$ quasars in cosmological hydrodynamic simulations. The underlying galaxy formation model is the same as in IllustrisTNG. Using constrained initial conditions, we study the growth of BHs in rare overdense regions (forming $\gtrsim10^{12}M_{\odot}/h$ halos by $z=7$) using a $(9~\mathrm{Mpc}/h)^3$ simulated volume. BH growth is maximal within halos that are compact and have a low tidal field. For these halos, we consider an array of gas-based seeding prescriptions wherein $M_{\mathrm{seed}}=10^4-10^6~M_{\odot}/h$ seeds are inserted in halos above critical thresholds for halo mass and dense, metal-poor gas mass (defined as $\tilde{M}_{\mathrm{h}}$ and $\tilde{M}_{\mathrm{sf,mp}}$, respectively, in units of $M_{\mathrm{seed}}$). We find that a seed model with $\tilde{M}_{\mathrm{sf,mp}}=5$ and $\tilde{M}_{\mathrm{h}}=3000$ successfully produces a $z\sim6$ quasar with $\sim10^9~M_{\odot}$ mass and $\sim10^{47}~\mathrm{ergs~s^ {-1}}$ luminosity. BH mergers play a crucial role at $z\gtrsim9$, causing an early boost in BH mass at a time when accretion-driven BH growth is negligible. When more stringent seeding conditions are applied (for e.g., $\tilde{M}_{\mathrm{sf,mp}}=1000$), the relative paucity of BH seeds results in a much lower merger rate. In this case, $z\gtrsim6$ quasars can only be formed if we enhance the maximum allowed BH accretion rates (by factors $\gtrsim10$) compared to the accretion model used in IllustrisTNG. This can be achieved either by allowing for super-Eddington accretion, or by reducing the radiative efficiency. Our results show that progenitors of $z\sim6$ quasars have distinct BH merger histories for different seeding models, which will be distinguishable with LISA observations.

preprint2022arXiv

Sunyaev-Zel&#39;dovich effect and X-ray scaling relations of galaxies, groups and clusters in the IllustrisTNG simulations

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

preprint2022arXiv

The buildup of strongly barred galaxies in the TNG100 simulation

We analyse the properties of strongly barred disc galaxies using the TNG100 simulation, a cosmological hydrodynamical realisation of the IllustrisTNG suite. We identify 270 disc galaxies at $z=0$ in the stellar mass range $M_{*}=10^{10.4-11}M_{\odot}$, of which 40 per cent are barred. Of the detected bars, more than half are strong. We find that the fraction of barred galaxies increases with stellar mass, in agreement with observational results. Strongly barred galaxies exhibit, overall, lower gas-to-stellar mass ratios compared to unbarred galaxies. The majority of barred galaxies are quenched (sSFR $\sim10^{-11.7} $yr$^{-1}$), whereas unbarred galaxies continue to be active (sSFR $\sim10^{-10.3}$yr$^{-1}$) on the main sequence of star-forming galaxies. We explore the evolution of strongly barred and unbarred galaxies to investigate their formation and quenching histories. We find that strong bars form between $0.5< z< 1.5$, with more massive galaxies hosting older bars. Strong bars form in galaxies with an early-established prominent disc component, undergoing periods of enhanced star formation and black hole accretion, possibly assisted by cosmological inflows. Unbarred galaxies, on the other hand, assemble most of their mass and disc component at late times. The nuclear region of strongly barred galaxies quenches shortly after bar formation, while unbarred galaxies remain active across time. Our findings are indicative of bar quenching, possibly assisted by nuclear feedback processes. We conclude that the cosmological environment, together with small scale feedback processes, determine the chances of a galaxy to form a bar and to rapidly quench its central region.

preprint2022arXiv

The Circumgalactic Medium from the CAMELS Simulations: Forecasting Constraints on Feedback Processes from Future Sunyaev-Zeldovich Observations

The cycle of baryons through the circumgalactic medium (CGM) is important to understand in the context of galaxy formation and evolution. In this study we forecast constraints on the feedback processes heating the CGM with current and future Sunyaev-Zeldovich (SZ) observations. To constrain these processes, we use a suite of cosmological simulations, the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS), that varies four different feedback parameters of two previously existing hydrodynamical simulations, IllustrisTNG and SIMBA. We capture the dependencies of SZ radial profiles on these feedback parameters with an emulator, calculate their derivatives, and forecast future constraints on these feedback parameters from upcoming experiments. We find that for a DESI-like (Dark Energy Spectroscopic Instrument) galaxy sample observed by the Simons Observatory all four feedback parameters are able to be constrained (some within the $10\%$ level), indicating that future observations will be able to further restrict the parameter space for these sub-grid models. Given the modeled galaxy sample and forecasted errors in this work, we find that the inner SZ profiles contribute more to the constraining power than the outer profiles. Finally, we find that, despite the wide range of AGN feedback parameter variation in the CAMELS simulation suite, we cannot reproduce the tSZ signal of galaxies selected by the Baryon Oscillation Spectroscopic Survey as measured by the Atacama Cosmology Telescope.

preprint2022arXiv

The distinct stellar-to-halo mass relations of satellite and central galaxies: insights from the IllustrisTNG simulations

We study the stellar-to-halo mass relation (SHMR) for central and satellite galaxies with total dynamical masses above 10^10.5 Msun using the suite of cosmological magneto-hydrodynamical simulations IllustrisTNG. In particular, we quantify environmental effects on satellite populations from TNG50, TNG100, and TNG300 located within the virial radius of group- and cluster-like hosts with total masses of 10^12-15.2 Msun. At fixed stellar mass, the satellite SHMR exhibits a distinct shift towards lower dynamical mass compared to the SHMR of centrals. Conversely, at fixed dynamical mass, satellite galaxies appear to have larger stellar-to-total mass fractions than centrals by up to a factor of a few. The systematic deviation from the central SHMR is larger for satellites in more massive hosts, at smaller cluster-centric distances, with earlier infall times, and that inhabit higher local density environments; moreover, it is in place already at early times (z < 2). Systematic environmental effects might contribute to the perceived galaxy-to-galaxy variation in the measured SHMR when galaxies cannot be separated into satellites and centrals. The SHMR of satellites exhibits a larger scatter than centrals, over the whole range of dynamical mass (by up to 0.8 dex). The shift of the satellite SHMR results mostly from tidal stripping of their dark matter, which affects satellites in an outside-in fashion: the departure of the satellite SHMR from the centrals&#39; relation diminishes for measurements of dynamical mass in progressively smaller apertures. Finally, we provide a family of fitting functions for the SHMR predicted by IllustrisTNG.

preprint2022arXiv

The evolution of the barred galaxy population in the TNG50 simulation

We use the magnetic-hydrodynamical simulation TNG50 to study the evolution of barred massive disc galaxies. Massive spiral galaxies are already present as early as $z=4$, and bar formation takes place already at those early times. The bars grow longer and stronger as the host galaxies evolve, with the bar sizes increasing at a pace similar to that of the disc scale lengths. The bar fraction mildly evolves with redshift for galaxies with $M_{*}\geq10^{10}M\odot$, being greater than $\sim40\%$ at $0.5<z<3$ and $\sim30\%$ at $z=0$. When bars larger than a given physical size ($\geq 2\,\rm kpc$) or the angular resolution limit of twice the I-band angular PSF FWHM of the HST are considered, the bar fraction dramatically decreases with increasing redshift, reconciling the theoretical predictions with observational data. We find that barred galaxies have an older stellar population, lower gas fractions and star formation rates than unbarred galaxies. In most cases, the discs of barred galaxies assembled earlier and faster than the discs of unbarred galaxies. We also find that barred galaxies are typical in haloes with larger concentrations and smaller spin parameters than unbarred galaxies. Furthermore, the inner regions of barred galaxies are more baryon-dominated than those of unbarred galaxies but have comparable global stellar mass fractions. Our findings suggest that the bar population could be used as a potential tracer of the buildup of disc galaxies and their host haloes. With this paper, we release a catalogue of barred galaxies in TNG50 at $6$ redshifts between $z=4$ and $z=0$.

preprint2022arXiv

The formation of low surface brightness galaxies in the IllustrisTNG simulation

We explore the nature of low surface brightness galaxies (LSBGs) in the hydrodynamic cosmological simulation TNG100 of the IllustrisTNG project, selecting a sample of LSBGs ($r$-band effective surface brightness $μ_r > 22.0$ mag arcsec$^{-2}$) at $z=0$ over a wide range of stellar masses ($M_{\ast} = 10^{9}$-$10^{12}$ M$_\odot$). We find LSBGs of all stellar masses, although they are particularly prevalent at $M_{\ast} < 10^{10}$ M$_\odot$. We show that the specific star formation rates of LSBGs are not significantly different from those of high surface brightness galaxies (HSBGs) but, as a population, LSBGs are systematically less massive and more extended than HSBGs, and tend to display late-type morphologies according to a kinematic criterion. At fixed stellar mass, we find that haloes hosting LSBGs are systematically more massive and have a higher baryonic fraction than those hosting HSBGs. We find that LSBGs have higher stellar specific angular momentum and halo spin parameter values compared to HSBGs, as suggested by previous works. We track the evolution of these quantities back in time, finding that the spin parameters of the haloes hosting LSBGs and HSBGs exhibit a clear bifurcation at $z \sim 2$, which causes a similar separation in the evolutionary tracks of other properties such as galactic angular momentum and effective radius, ultimately resulting in the values observed at $z =$ 0. The higher values of specific stellar angular momentum and halo spin in LSBGs seem to be responsible for their extended nature, preventing material from collapsing into the central regions of the galaxies, also causing LSBGs to host less massive black holes at their centres.

preprint2022arXiv

The observability of galaxy merger signatures in nearby gas-rich spirals

Galaxy mergers are crucial to understanding galaxy evolution, therefore we must determine their observational signatures to select them from large IFU galaxy samples such as MUSE and SAMI. We employ 24 high-resolution idealised hydrodynamical galaxy merger simulations based on the &#34;Feedback In Realistic Environment&#34; (FIRE-2) model to determine the observability of mergers to various configurations and stages using synthetic images and velocity maps. Our mergers cover a range of orbital configurations at fixed 1:2.5 stellar mass ratio for two gas rich spirals at low redshift. Morphological and kinematic asymmetries are computed for synthetic images and velocity maps spanning each interaction. We divide the interaction sequence into three: (1) the pair phase; (2) the merging phase; and (3) the post-coalescence phase. We correctly identify mergers between first pericentre passage and 500 Myr after coalescence using kinematic asymmetry with 66% completeness, depending upon merger phase and the field-of-view of the observation. We detect fewer mergers in the pair phase (40%) and many more in the merging and post-coalescence phases (97%). We find that merger detectability decreases with field-of-view, except in retrograde mergers, where centrally concentrated asymmetric kinematic features enhances their detectability. Using a cut-off derived from a combination of photometric and kinematic asymmetry, we increase these detections to 89% overall, 79% in pairs, and close to 100% in the merging and post-coalescent phases. By using this combined asymmetry cut-off we mitigate some of the effects caused by smaller fields-of-view subtended by massively multiplexed integral field spectroscopy programmes.

preprint2022arXiv

The THESAN project: predictions for multi-tracer line intensity mapping in the epoch of reionization

Line intensity mapping (LIM) is rapidly emerging as a powerful technique to study galaxy formation and cosmology in the high-redshift Universe. We present LIM estimates of select spectral lines originating from the interstellar medium (ISM) of galaxies and 21 cm emission from neutral hydrogen gas in the Universe using the large volume, high resolution THESAN reionization simulations. A combination of sub-resolution photo-ionization modelling for HII regions and Monte Carlo radiative transfer calculations is employed to estimate the dust-attenuated spectral energy distributions (SEDs) of high-redshift galaxies ($z\gtrsim5.5$). We show that the derived photometric properties such as the ultraviolet (UV) luminosity function and the UV continuum slopes match observationally inferred values, demonstrating the accuracy of the SED modelling. We provide fits to the luminosity--star formation rate relation (L-SFR) for the brightest emission lines and find that important differences exist between the derived scaling relations and the widely used low-$z$ ones because the interstellar medium of reionization era galaxies is generally less metal-enriched than in their low redshift counterparts. We use these relations to construct line intensity maps of nebular emission lines and cross correlate with the 21 cm emission. Interestingly, the wavenumber at which the correlation switches sign ($k_\mathrm{transition}$) depends heavily on the reionization model and to a lesser extent on the targeted emission line, which is consistent with the picture that $k_\mathrm{transition}$ probes the typical sizes of ionized regions. The derived scaling relations and intensity maps represent a timely state-of-the-art framework for forecasting and interpreting results from current and upcoming LIM experiments.

preprint2022arXiv

The THESAN project: properties of the intergalactic medium and its connection to Reionization-era galaxies

The high-redshift intergalactic medium (IGM) and the primeval galaxy population are rapidly becoming the new frontier of extra-galactic astronomy. We investigate the IGM properties and their connection to galaxies at $z\geq5.5$ under different assumptions for the ionizing photon escape and the nature of dark matter, employing our novel THESAN radiation-hydrodynamical simulation suite, designed to provide a comprehensive picture of the emergence of galaxies in a full reionization context. Our simulations have realistic `late&#39; reionization histories, match available constraints on global IGM properties and reproduce the recently-observed rapid evolution of the mean free path of ionizing photons. We additionally examine high-z Lyman-$α$ transmission. The optical depth evolution is consistent with data, and its distribution suggests an even-later reionization than simulated, although with a strong sensitivity to the source model. We show that the effects of these two unknowns can be disentangled by characterising the spectral shape and separation of Lyman-$α$ transmission regions, opening up the possibility to observationally constrain both. For the first time in simulations, THESAN reproduces the modulation of the Lyman-$α$ flux as a function of galaxy distance, demonstrating the power of coupling a realistic galaxy formation model with proper radiation-hydrodynamics. We find this feature to be extremely sensitive on the timing of reionization, while being relatively insensitive to the source model. Overall, THESAN produces a realistic IGM and galaxy population, providing a robust framework for future analysis of the high-z Universe.

preprint2022arXiv

Unifying Sunyaev-Zel&#39;dovich and X-ray predictions from clusters to galaxy groups: the impact of X-ray mass estimates on the $Y-M$ scaling relation

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

preprint2022arXiv

UV to submillimetre luminosity functions of TNG50 galaxies

We apply the radiative transfer (RT) code SKIRT on a sample of ~14000 low-redshift (z<= 0.1) galaxies extracted from the TNG50 simulation to enable an apples-to-apples comparison with observations. The RT procedure is calibrated via comparison of a subsample of TNG50 galaxies with the DustPedia observational sample: we compare several luminosity and colour scaling relations and spectral energy distributions in different specific SFR bins. We consistently derive galaxy luminosity functions for the TNG50 simulation in 14 broadband filters from UV to submillimetre wavelengths and investigate the effects of the aperture, orientation, radiative transfer recipe, and numerical resolution. We find that, while our TNG50+RT fiducial model agrees well with the observed luminosity functions at the knee (+/- 0.04 dex typical agreement), the TNG50+RT luminosity functions evaluated within 5R_1/2 are generally higher than observed at both the faint and bright ends, by 0.004 (total IR)-0.27 (UKIDSS H) dex and 0.12 (SPIRE250)-0.8 (GALEX FUV) dex, respectively. A change in the aperture does affect the bright end of the luminosity function, easily by up to 1 dex depending on the choice. However, we also find that the galaxy luminosity functions of a worse-resolution run of TNG50 (TNG50-2, with 8 times worse mass resolution than TNG50, similar to TNG100) are in better quantitative agreement with observational constraints. Finally, we publicly release the photometry for the TNG50 sample in 53 broadbands from FUV to submillimetre, in three orientations and four apertures, as well as galaxy spectral energy distributions.

preprint2021arXiv

Finding universal relations in subhalo properties with artificial intelligence

We use a generic formalism designed to search for relations in high-dimensional spaces to determine if the total mass of a subhalo can be predicted from other internal properties such as velocity dispersion, radius, or star-formation rate. We train neural networks using data from the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project and show that the model can predict the total mass of a subhalo with high accuracy: more than 99% of the subhalos have a predicted mass within 0.2 dex of their true value. The networks exhibit surprising extrapolation properties, being able to accurately predict the total mass of any type of subhalo containing any kind of galaxy at any redshift from simulations with different cosmologies, astrophysics models, subgrid physics, volumes, and resolutions, indicating that the network may have found a universal relation. We then use different methods to find equations that approximate the relation found by the networks and derive new analytic expressions that predict the total mass of a subhalo from its radius, velocity dispersion, and maximum circular velocity. We show that in some regimes, the analytic expressions are more accurate than the neural networks. We interpret the relation found by the neural network and approximated by the analytic equation as being connected to the virial theorem.

preprint2021arXiv

High and low Sérsic index bulges in Milky Way- and M31-like galaxies: origin and connection to the bar with TNG50

We study bulge formation in MW/M31-like galaxies in a $Λ$-cold dark matter scenario, focusing on the origin of high- and low-Sersic index bulges. For this purpose we use TNG50, a simulation of the IllustrisTNG project that combines a resolution of $\sim 8 \times 10^4 M_{\odot}$ in stellar particles with a cosmological volume 52 cMpc in extent. We parametrize bulge surface brightness profiles by the Sérsic index and the bulge-to-total (B/T) ratio obtained from two-component photometric decompositions. In our sample of 287 MW/M31-like simulated galaxies, $17.1\%$ of photometric bulges exhibit high-Sérsic indices and $82.9\%$ show low-Sérsic indices. We study the impact that the environment, mergers and bars have in shaping the surface brightness profiles. We explore two different definitions for local environment and find no correlation between bulge properties and the environment where they reside. Simulated galaxies with higher Sérsic indices show, on average, a higher fraction of ex-situ stars in their kinematically selected bulges. For this bulge population the last significant merger (total mass ratio $m_{\rm sat}/m_{\rm host} > 0.1$) occurs, on average, at later times. However, a substantial fraction of low-Sérsic index bulges also experience a late significant merger. We find that bars play an important role in the development of the different types of photometric bulges. We show that the fraction of simulated galaxies with strong bars is smaller for the high- than for the low-Sérsic index population, reaching differences of $20\%$ at $z > 1$. Simulated galaxies with high fractions of ex-situ stars in the bulge do not develop strong bars. Conversely, simulated galaxies with long-lived strong bars have bulges with ex-situ fractions, $f_{\rm ex-situ} < 0.2$.

preprint2021arXiv

Hot and counter-rotating star-forming disk galaxies in IllustrisTNG and their real-world counterparts

A key feature of a large population of low-mass, late-type disk galaxies are star-forming disks with exponential light distributions. They are typically also associated with thin and flat morphologies, blue colors, and dynamically cold stars moving along circular orbits within co-planar thin gas disks. However, the latter features do not necessarily always imply the former, in fact, a variety of different kinematic configurations do exist. In this work, we use the cosmological hydrodynamical IllustrisTNG Simulation to study the nature and origin of dynamically hot, sometimes even counter-rotating, star-forming disk galaxies in the lower stellar mass range (between $5\times 10^9\,\mathrm{M_{\odot}}$ and $2\times 10^{10}\,\mathrm{M_{\odot}}$). We find that being dynamically hot arises in most cases as an induced transient state, for example due to galaxy interactions and merger activities, rather than as an age-dependent evolutionary phase of star-forming disk galaxies. The dynamically hot but still actively star-forming disks show a common feature of hosting kinematically misaligned gas and stellar disks, and centrally concentrated on-going star formation. The former is often accompanied by disturbed gas morphologies, while the latter is reflected in low gas and stellar spins in comparison to their dynamically cold, normal disk counterparts. Interestingly, observed galaxies from MaNGA with kinematic misalignment between gas and stars show remarkably similar general properties as the IllustrisTNG galaxies, and therefore are plausible real-world counterparts. In turn, this allows us to make predictions for the stellar orbits and gas properties of these misaligned galaxies.

preprint2021arXiv

On the formation of massive quiescent galaxies with diverse morphologies in the TNG50 simulation

Observations have shown that the star-formation activity and the morphology of galaxies are closely related, but the underlying physical connection is not well understood. Using the TNG50 simulation, we explore the quenching and the morphological evolution of the 102 massive quiescent galaxies in the mass range of $10.5<\log(M_{\rm stellar}/M_{\odot})<11.5$ selected at $z=0$. The morphology of galaxies is quantified based on their kinematics, and we measure the quenching timescale of individual galaxies directly from star formation history. We show that galaxies tend to be quenched more rapidly if they: (i) are satellites in massive halos, (ii) have lower star-forming gas fractions, or (iii) inject a larger amount of black hole kinetic feedback energy. By following the global evolutionary pathways, we conclude that quiescent discs are mainly disc galaxies that are recently and slowly quenched. Approximately half of the quiescent ellipticals at $z=0$ are rapidly quenched at higher redshifts while still disc-like. While being quiescent, they gradually become more elliptical mostly by disc heating, yet these ellipticals still retain some degree of rotation. The other half of quiescent ellipticals with the most random motion-dominated kinematics build up large spheroidal components before quenching primarily by mergers, or in some cases, misaligned gas accretion. However, the mergers that contribute to morphological transformation do not immediately quench galaxies in many cases. In summary, we find that quenching and morphological transformation are decoupled. We conclude that the TNG black hole feedback -- in combination with the stochastic merger history of galaxies -- leads to a large diversity of quenching timescales and a rich morphological landscape.

preprint2021arXiv

Submillimetre galaxies in cosmological hydrodynamical simulations -- an opportunity for constraining feedback models

Submillimetre galaxies (SMGs) have long posed a challenge for theorists, and self-consistently reproducing the properties of the SMG population in a large-volume cosmological hydrodynamical simulation has not yet been achieved. We use a scaling relation derived from previous simulations plus radiative transfer calculations to predict the submm flux densities of simulated SMGs drawn from cosmological simulations from the Illustris and IllustrisTNG projects based on the simulated galaxies&#39; star formation rates (SFRs) and dust masses and compare the predicted number counts with observations. We find that the predicted SMG number counts based on IllustrisTNG are significantly less than observed (more than 1 dex at $S_{850} \gtrsim 4$ mJy). The simulation from the original Illustris project yields more SMGs than IllustrisTNG: the predicted counts are consistent with those observed at both $S_{850} \lesssim 5$ mJy and $S_{850} \gtrsim 9$ mJy and only a factor of $\sim 2$ lower than observed at intermediate flux densities. The redshift distribution of SMGs with $S_{850} > 3$ mJy in IllustrisTNG is consistent with the observed distribution, whereas the Illustris redshift distribution peaks at significantly lower redshift (1.5 vs. 2.8). We demonstrate that IllustrisTNG hosts fewer SMGs than Illustris because in the former, high-mass ($M_{\star} \sim 10^{11} \, \text{M}_{\odot}$) $z \sim 2-3$ galaxies have lower dust masses and SFRs than in Illustris owing to differences in the sub-grid models for stellar or/and active galactic nucleus (AGN) feedback between the two simulations (we unfortunately cannot isolate the specific cause(s) post hoc). Our results demonstrate that because our method enables predicting SMG number counts in post-processing with a negligible computational expense, SMGs can provide useful constraints for tuning sub-grid models in future large-volume cosmological simulations.

preprint2021arXiv

The IllustrisTNG Simulations: Public Data Release

We present the full public release of all data from the TNG50, TNG100 and TNG300 simulations of the IllustrisTNG project. IllustrisTNG is a suite of large volume, cosmological, gravo-magnetohydrodynamical simulations run with the moving-mesh code Arepo. TNG includes a comprehensive model for galaxy formation physics, and each TNG simulation self-consistently solves for the coupled evolution of dark matter, cosmic gas, luminous stars, and supermassive blackholes from early time to the present day, z=0. Each of the flagship runs -- TNG50, TNG100, and TNG300 -- are accompanied by lower-resolution and dark-matter only counterparts, and we discuss scientific and numerical cautions and caveats relevant when using TNG. Full volume snapshots are available at 100 redshifts; halo and subhalo catalogs at each snapshot and merger trees are also released. The data volume now directly accessible online is ~1.1 PB, including 2,000 full volume snapshots and ~110,000 high time-resolution subbox snapshots. Data access and analysis examples are available in IDL, Python, and Matlab. We describe improvements and new functionality in the web-based API, including on-demand visualization and analysis of galaxies and halos, exploratory plotting of scaling relations and other relationships between galactic and halo properties, and a new JupyterLab interface. This provides an online, browser-based, near-native data analysis platform which supports user computation with fully local access to TNG data, alleviating the need to download large simulated datasets.

preprint2021arXiv

The impact of galactic feedback on the shapes of dark-matter haloes

We quantify the impact of galaxy formation on dark matter halo shapes using cosmological simulations at redshift $z=0$. The haloes are drawn from the IllustrisTNG project, a suite of magneto-hydrodynamic simulations of galaxies. We focus on haloes of mass $10^{10-14} M_\odot$ from the 50-Mpc (TNG50) and 100-Mpc (TNG100) boxes, and compare them to dark matter-only (DMO) analogues and other simulations e.g. NIHAO and Eagle. We further quantify the prediction uncertainty by varying the baryonic feedback models in a series of smaller 25 Mpc $h^{-1}$ boxes. We find that: (i) galaxy formation results in rounder haloes compared to the DMO simulations, in qualitative agreement with past hydrodynamic models. Haloes of mass $\approx 2\times 10^{12} M_\odot$ are most spherical, with an average minor-to-major axis ratio of $\left< s \right> \approx 0.75$ in the inner halo, an increase of 40 per cent compared to their DMO counterparts. No significant change in halo shape is found for low-mass $10^{10} M_\odot$ haloes; (ii) stronger feedback, e.g. increasing galactic wind speed, reduces the impact of baryons; (iii) the inner halo shape correlates with the stellar mass fraction, which can explain the dependence of halo shapes on different feedback models; (iv) the fiducial and weaker feedback models are most consistent with observational estimates of the Milky Way halo shape. Yet, at fixed halo mass, very diverse and possibly unrealistic feedback models all predict inner halo shapes that are closer to one another than to the DMO results. This implies that a larger observational sample would be required to statistically distinguish between different baryonic prescriptions due to large halo-to-halo variation in halo shapes.

preprint2020arXiv

A missing outskirts problem? Comparisons between stellar halos in the Dragonfly Nearby Galaxies Survey and the TNG100 simulation

Low surface brightness galactic stellar halos provide a challenging but promising path towards unraveling the past assembly histories of individual galaxies. Here, we present detailed comparisons between the stellar halos of Milky Way-mass disk galaxies observed as part of the Dragonfly Nearby Galaxies Survey (DNGS) and stellar mass-matched galaxies in the TNG100 run of the IllustrisTNG project. We produce stellar mass maps as well as mock $g$ and $r$-band images for randomly oriented simulated galaxies, convolving the latter with the Dragonfly PSF and taking care to match the background noise, surface brightness limits and spatial resolution of DNGS. We measure azimuthally averaged stellar mass density and surface brightness profiles, and find that the DNGS galaxies generally have less stellar mass (or light) at large radii (>20 kpc) compared to their mass-matched TNG100 counterparts, and that simulated galaxies with similar surface density profiles tend to have low accreted mass fractions for their stellar mass. We explore potential solutions to this apparent &#34;missing outskirts problem&#34; by implementing several ad-hoc adjustments within TNG100 at the stellar particle level. Although we are unable to identify any single adjustment that fully reconciles the differences between the observed and simulated galaxy outskirts, we find that artificially delaying the disruption of satellite galaxies and reducing the spatial extent of in-situ stellar populations result in improved matches between the outer profile shapes and stellar halo masses, respectively. Further insight can be achieved with higher resolution simulations that are able to better resolve satellite accretion, and with larger samples of observed galaxies.

preprint2020arXiv

A redshift-dependent IRX-$β$ dust attenuation relation for TNG50 galaxies

We study the relation between the UV-slope, $β$, and the ratio between the infrared- and UV-luminosities (IRX) of galaxies from TNG50, the latest installment of the IllustrisTNG galaxy formation simulations. We select 7280 star-forming main-sequence (SFMS) galaxies with stellar mass $10^9M_\odot$ at redshifts $0 \leq z \leq 4$ and perform radiative transfer with SKIRT to model effects of interstellar medium dust on the emitted stellar light. Assuming a Milky Way (MW) dust type and a dust-to-metal ratio of 0.3, we find that TNG50 SFMS galaxies generally agree with observationally-derived IRX-$β$ relations at $z \lesssim 1$. However, we find a redshift-dependent systematic offset with respect to empirically-derived local relations, with the TNG50 IRX-$β$ relation shifting towards lower $β$ and steppening at higher redshifts. This is partially driven by variations in the dust-uncorrected UV-slope of galaxies, due to different star-formation histories of galaxies selected at different cosmic epochs; we suggest the remainder of the effect is caused by differences in the effective dust attenuation curves (EDACs) of galaxies as a function of redshift. We find a typical galaxy-to-galaxy variation of 0.3 dex in IRX at fixed $β$, correlated with intrinsic galaxy properties: galaxies with higher star-formation rates, star-formation efficiencies, gas metallicities and stellar masses exhibit larger IRX values. We demonstrate a degeneracy between stellar age, dust geometry and dust composition: $z=4$ galaxies with a Small Magellanic Cloud dust type follow the same IRX-$β$ relation as low-redshift galaxies with MW dust. We provide a redshift-dependent fitting function for the IRX-$β$ relation for MW dust based on our models.

preprint2020arXiv

Atomic hydrogen in IllustrisTNG galaxies: the impact of environment parallelled with local 21-cm surveys

We investigate the influence of environment on the cold-gas properties of galaxies at z=0 within the TNG100 cosmological, magnetohydrodynamic simulation, part of the IllustrisTNG suite. We extend previous post-processing methods for breaking gas cells into their atomic and molecular phases, and build detailed mocks to comprehensively compare to the latest surveys of atomic hydrogen (HI) in nearby galaxies, namely ALFALFA and xGASS. We use TNG100 to explore the HI content, star formation activity, and angular momentum of satellite galaxies, each as a function of environment, and find that satellites are typically a factor of ~3 poorer in HI than centrals of the same stellar mass, with the exact offset depending sensitively on parent halo mass. Due to the large physical scales on which HI measurements are made (~45--245 kpc), contributions from gas not bound to the galaxy of interest but in the same line of sight crucially lead to larger HI mass measurements in the mocks in many cases, ultimately aligning with observations. This effect is mass-dependent and naturally greater for satellites than centrals, as satellites are never isolated by definition. We also show that HI stripping in TNG100 satellites is closely accompanied by quenching, in tension with observational data that instead favour that HI is preferentially stripped before star formation is reduced.

preprint2020arXiv

Cosmological insights into the assembly of the radial and compact stellar halo of the Milky Way

Recent studies using Gaia DR2 have identified a massive merger in the history of the Milky Way (MW) whose debris is markedly radial and counterrotating. This event, known as the Gaia-Enceladus/Gaia-Sausage (GE/GS), is also hypothesized to have built the majority of the inner stellar halo. We use the cosmological hydrodynamic simulation Illustris to place this merger in the context of galaxy assembly within $Λ$CDM. From $\sim$150 MW analogs, $\sim 80 \%$ have experienced at least one merger of similar mass and infall time as GE/GS. Within this sample, 37 have debris as radial as that of the GE/GS, which we dub the Ancient Radial Mergers (ARMs). Counterrotation is not rare among ARMs, with $43 \%$ having $>40 \%$ of their debris in counterrotating orbits. However, the compactness inferred for the GE/GS debris, given its large $β$ and its substantial contribution to the stellar halo, is more difficult to reproduce. The median radius of ARM debris is r$_{*,deb}\simeq 45$kpc, while GE/GS is thought to be mostly contained within $r\sim 30$ kpc. For most MW analogs, a few mergers are required to build the inner stellar halo, and ARM debris only accounts for $\sim 12 \%$ of inner accreted stars. Encouragingly, we find one ARM that is both compact and dominates the inner halo of its central, making it our best GE/GS analog. Interestingly, this merger deposits a significant number of stars (M$_*\simeq1.5 \times 10^9 M_\odot$) in the outer halo, suggesting that an undiscovered section of GE/GS may await detection.

preprint2020arXiv

Dust entrainment in galactic winds

Winds driven by stellar feedback are an essential part of the galactic ecosystem and are the main mechanism through which low-mass galaxies regulate their star formation. These winds are generally observed to be multi-phase with detections of entrained neutral and molecular gas. They are also thought to enrich the circum-galactic medium around galaxies with metals and dust. This ejected dust encodes information about the integrated star formation and outflow history of the galaxy. It is therefore, important to understand how much dust is entrained and driven out of the disc by galactic winds. Here we demonstrate that stellar feedback is efficient in driving dust-enriched winds and eject enough material to account for the amount of extraplanar dust observed in nearby galaxies. The amount of dust in the wind depends on the sites from where they are launched, with dustier galaxies launching more dust enriched outflows. Moreover, the outflowing cold-dense gas is significantly more dust-enriched than the volume filling hot tenuous material, naturally reproducing the complex multiphase structure of the outflowing wind observed in nearby galaxies. These results provide an important new insight into the dynamics, structure, and composition of galactic winds and their role in determining the dust content of the extragalactic gas in galaxies.

preprint2020arXiv

Early-type galaxy density profiles from IllustrisTNG: I. Galaxy correlations and the impact of baryons

We explore the isothermal total density profiles of early-type galaxies (ETGs) in the IllustrisTNG simulation. For the selected 559 ETGs at $z = 0$ with stellar mass $10^{10.7}\mathrm{M}_{\odot} \leqslant M_{\ast} \leqslant 10^{11.9}\mathrm{M}_{\odot}$, the total power-law slope has a mean of $\langleγ^{\prime}\rangle = 2.011 \pm 0.007$ and a scatter of $σ_{γ^{\prime}} = 0.171$ over the radial range 0.4 to 4 times the stellar half mass radius. Several correlations between $γ^{\prime}$ and galactic properties including stellar mass, effective radius, stellar surface density, central velocity dispersion, central dark matter fraction and in-situ-formed stellar mass ratio are compared to observations and other simulations, revealing that IllustrisTNG reproduces many correlation trends, and in particular, $γ^{\prime}$ is almost constant with redshift below $z = 2$. Through analyzing IllustrisTNG model variations we show that black hole kinetic winds are crucial to lowering $γ^{\prime}$ and matching observed galaxy correlations. The effects of stellar winds on $γ^{\prime}$ are subdominant compared to AGN feedback, and differ due to the presence of AGN feedback from previous works. The density profiles of the ETG dark matter halos are well-described by steeper-than-NFW profiles, and they are steeper in the full physics (FP) run than their counterparts in the dark matter only (DMO) run. Their inner density slopes anti-correlates (remain constant) with the halo mass in the FP (DMO) run, and anti-correlates with the halo concentration parameter $c_{200}$ in both types of runs. The dark matter halos of low-mass ETGs are contracted whereas high-mass ETGs are expanded, suggesting that variations in the total density profile occur through the different halo responses to baryons.

preprint2020arXiv

Ejective and preventative: the IllustrisTNG black hole feedback and its effects on the thermodynamics of the gas within and around galaxies

Supermassive black holes (SMBHs) which reside at the centres of galaxies can inject vast amounts of energy into the surrounding gas and are thought to be a viable mechanism to quench star-formation in massive galaxies. Here we study the $10^{9\textrm{--}12.5}\,\mathrm{M_\odot}$ stellar mass central galaxy population of the IllustrisTNG simulation, specifically the TNG100 and TNG300 volumes at \zeq{0}, and show how the three components -- SMBH, galaxy, and circumgalactic medium (CGM) -- are interconnected in their evolution. We find that gas entropy is a sensitive diagnostic of feedback injection. In particular, we demonstrate how the onset of the low-accretion BH feedback mode, realised in the IllustrisTNG model as a kinetic, BH-driven wind, leads not only to star-formation quenching at stellar masses $\gtrsim10^{10.5}\mathrm{M_\odot}$ but also to a change in thermodynamic properties of the (\emph{non}-star-forming) gas, both within the galaxy and beyond. The IllustrisTNG kinetic feedback from SMBHs increases the average gas entropy, within the galaxy and in the CGM, lengthening typical gas cooling times from $10\textrm{--}100\,\mathrm{Myr}$ to $1\textrm{--}10\,\mathrm{Gyr}$, effectively ceasing ongoing star-formation and inhibiting radiative cooling and future gas accretion. In practice, the same AGN feedback channel is simultaneously `ejective&#39; and `preventative&#39; and leaves an imprint on the temperature, density, entropy, and cooling times also in the outer reaches of the gas halo, up to distances of several hundred kiloparsecs. In the IllustrisTNG model, a long-lasting quenching state can occur for a heterogeneous CGM, whereby the hot and dilute CGM gas of quiescent galaxies contains regions of low-entropy gas with short cooling times.

preprint2020arXiv

First results from SMAUG: Characterization of Multiphase Galactic Outflows from a Suite of Local Star-Forming Galactic Disk Simulations

Large scale outflows in star-forming galaxies are observed to be ubiquitous, and are a key aspect of theoretical modeling of galactic evolution in a cosmological context, the focus of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project. Gas blown out from galactic disks, similar to gas within galaxies, consists of multiple phases with large contrasts of density, temperature, and other properties. To study multiphase outflows as emergent phenomena, we run a suite of ~pc-resolution local galactic disk simulations using the TIGRESS framework. Explicit modeling of the interstellar medium (ISM), including star formation and self-consistent radiative heating plus supernova feedback, regulates ISM properties and drives the outflow. We investigate the scaling of outflow mass, momentum, energy, and metal loading factors with galactic disk properties, including star formation rate (SFR) surface density (Σ_SFR~10^{-4}-1 M_sun/kpc^2/yr), gas surface density (~1-100 M_sun/pc^2), and total midplane pressure (or weight) (~10^3-10^6 k_B cm^{-3} K). The main components of outflowing gas are mass-delivering cool gas (T~10^4 K) and energy/metal-delivering hot gas (T~10^6 K). Cool mass outflow rates measured at outflow launch points (one or two scale heights) are 1-100 times the SFR (decreasing with Σ_SFR), although in massive galaxies most mass falls back due to insufficient outflow velocity. The hot galactic outflow carries mass comparable to 10% of the SFR, together with 10-20% of the energy and 30-60% of the metal mass injected by SN feedback. The characteristic outflow velocities of both phases scale very weakly with SFR, as v_out \propto Σ_SFR^{0.1~0.2}, consistent with observations. Importantly, our analysis demonstrates that in any physically-motivated cosmological wind model, it is crucial to include at least two distinct thermal wind components.

preprint2020arXiv

First results from SMAUG: Insights into star formation conditions from spatially-resolved ISM properties in TNG50

Physical and chemical properties of the interstellar medium (ISM) at sub-galactic ($\sim$kpc) scales play an indispensable role in controlling the ability of gas to form stars. As part of the SMAUG (Simulating Multiscale Astrophysics to Understand Galaxies) project, in this paper, we use the TNG50 cosmological simulation to explore the physical parameter space of 8 resolved ISM properties in star-forming regions to constrain the areas of this hyperspace over which most star-forming environments exist. We deconstruct our simulated galaxies spanning a wide range of mass (M$_\star = 10^{7-11}$ M$_\odot$) and redshift ($0 \leq z \leq 3$) into kpc-sized regions, and statistically analyze the gas/stellar surface densities, gas metallicity, vertical stellar velocity dispersion, epicyclic frequency and dark-matter volumetric density representative of each region in the context of their star formation activity and galactic environment (radial galactocentric location). By examining the star formation rate (SFR) weighted distributions of these properties, we show that stars primarily form in two spatially distinct environmental regimes, which are brought about by an underlying bi-component radial SFR surface density profile in galaxies. We examine how the relative prominence of these two regimes depends on host galaxy mass and cosmic time. We also compare our findings with those from integral field spectroscopy observations and achieve a good overall agreement. Further, using dimensionality reduction, we characterise the aforementioned hyperspace to reveal a high-degree of multicollinearity in relationships amongst ISM properties that drive the distribution of star formation at kpc-scales. Based on this, we show that a reduced 3D representation underpinned by a multi-variate radius relationship is sufficient to capture most of the variance in the original 8D space.

preprint2020arXiv

High redshift JWST predictions from IllustrisTNG: II. Galaxy line and continuum spectral indices and dust attenuation curves

We present predictions for high redshift ($z=2-10$) galaxy populations based on the IllustrisTNG simulation suite and a full Monte Carlo dust radiative transfer post-processing. Specifically, we discuss the ${\rm H}_α$ and ${\rm H}_β$ + $[\rm O \,III]$ luminosity functions up to $z=8$. The predicted ${\rm H}_β$ + $[\rm O \,III]$ luminosity functions are consistent with present observations at $z\lesssim 3$ with $\lesssim 0.1\,{\rm dex}$ differences in luminosities. However, the predicted ${\rm H}_α$ luminosity function is $\sim 0.3\,{\rm dex}$ dimmer than the observed one at $z\simeq 2$. Furthermore, we explore continuum spectral indices, the Balmer break at $4000$Å(D4000) and the UV continuum slope $β$. The median D4000 versus sSFR relation predicted at $z=2$ is in agreement with the local calibration despite a different distribution pattern of galaxies in this plane. In addition, we reproduce the observed $A_{\rm UV}$ versus $β$ relation and explore its dependence on galaxy stellar mass, providing an explanation for the observed complexity of this relation. We also find a deficiency in heavily attenuated, UV red galaxies in the simulations. Finally, we provide predictions for the dust attenuation curves of galaxies at $z=2-6$ and investigate their dependence on galaxy colors and stellar masses. The attenuation curves are steeper in galaxies at higher redshifts, with bluer colors, or with lower stellar masses. We attribute these predicted trends to dust geometry. Overall, our results are consistent with present observations of high redshift galaxies. Future JWST observations will further test these predictions.

preprint2020arXiv

Kinematic decomposition of IllustrisTNG disk galaxies: morphology and relation with morphological structures

We recently developed an automated method, auto-GMM to decompose simulated galaxies. It extracts kinematic structures in an accurate, efficient, and unsupervised way. We use auto-GMM to study the stellar kinematic structures of disk galaxies from the TNG100 run of IllustrisTNG. We identify four to five structures that are commonly present among the diverse galaxy population. Structures having strong to moderate rotation are defined as cold and warm disks, respectively. Spheroidal structures dominated by random motions are classified as bulges or stellar halos, depending on how tightly bound they are. Disky bulges are structures that have moderate rotation but compact morphology. Across all disky galaxies and accounting for the stellar mass within 3 half-mass radii, the kinematic spheroidal structures, obtained by summing up stars of bulges and halos, contribute ~45% of the total stellar mass, while the disky structures constitute 55%. This study also provides important insights about the relationship between kinematically and morphologically derived galactic structures. Comparing the morphology of kinematic structures with that of traditional bulge+disk decomposition, we conclude: (1) the morphologically decomposed bulges are composite structures comprised of a slowly rotating bulge, an inner halo, and a disky bulge; (2) kinematically disky bulges, akin to what are commonly called pseudo bulges in observations, are compact disk-like components that have rotation similar to warm disks; (3) halos contribute almost 30% of the surface density of the outer part of morphological disks when viewed face-on; and (4) both cold and warm disks are often truncated in central regions.

preprint2020arXiv

Limitations to the &#34;basic&#39;&#39; HOD model and beyond

We make use of the IllustrisTNG cosmological, hydrodynamical simulations to test fundamental assumptions of the mass-based Halo Occupation Distribution (HOD) approach to modelling the galaxy-halo connection. By comparing the clustering of galaxies measured in the 300 Mpc TNG box (TNG300) with that predicted by the standard (``basic&#39;&#39;) HOD model, we find that, on average, the ``basic&#39;&#39; HOD model underpredicts the real-space correlation function in the TNG300 box by $\sim$ 15\% on scales of $1 \ {\rm Mpc}/h < r < 20 \ {\rm Mpc}/h$, which is well beyond the target precision demanded of next-generation galaxy redshift surveys. We perform several tests to establish the robustness of our findings to systematic effects, including the effect of finite box size and the choice of halo finder. In our exploration of ``secondary&#39;&#39; parameters with which to augment the ``basic&#39;&#39; HOD, we find that the local environment of the halo, the velocity dispersion anisotropy, $β$, and the product of the half-mass radius and the velocity dispersion, $σ^2 R_{\rm halfmass}$, are the three most effective measures of assembly bias that help reconcile the ``basic&#39;&#39; HOD-predicted clustering with that in TNG300. In addition, we test other halo properties such as halo spin, formation epoch and halo concentration. We also find that at fixed halo mass, galaxies in one type of environment cluster differently from galaxies in another. We demonstrate that a more complete model of the galaxy-halo connection can be constructed if we combine both mass and local environment information about the halo.

preprint2020arXiv

Redshift Evolution of the Fundamental Plane Relation in the IllustrisTNG Simulation

We investigate the fundamental plane (FP) evolution of early-type galaxies in the IllustrisTNG-100 simulation (TNG100) from redshift $z=0$ to $z=2$. We find that a tight plane relation already exists as early as $z=2$. Its scatter stays as low as $\sim 0.08$ dex across this redshift range. Both slope parameters $b$ and $c$ (where $R \propto σ^b I^c$ with $R$, $σ$, and $I$ being the typical size, velocity dispersion, and surface brightness) of the plane evolve mildly since $z=2$, roughly consistent with observations. The FP residual $\rm Res$ ($\equiv\,a\,+\,b\log σ\,+\,c\log I\,-\,\log R$, where $a$ is the zero point of the FP) is found to strongly correlate with stellar age, indicating that stellar age can be used as a crucial fourth parameter of the FP. However, we find that $4c+b+2=δ$, where $δ\sim 0.8$ for FPs in TNG, rather than zero as is typically inferred from observations. This implies that a tight power-law relation between the dynamical mass-to-light ratio $M_{\rm dyn}/L$ and the dynamical mass $M_{\rm dyn}$ (where $M_{\rm dyn}\equiv 5σ^2R/G$, with $G$ being the gravitational constant) is not present in the TNG100 simulation. Recovering such a relation requires proper mixing between dark matter and baryons, as well as star formation occurring with correct efficiencies at the right mass scales. This represents a powerful constraint on the numerical models, which has to be satisfied in future hydrodynamical simulations.

preprint2020arXiv

Resolving small-scale cold circumgalactic gas in TNG50

We use the high-resolution TNG50 cosmological magnetohydrodynamical simulation to explore the properties and origin of cold circumgalactic medium (CGM) gas around massive galaxies (M* > 10^11 Msun) at intermediate redshift (z~0.5). We discover a significant abundance of small-scale, cold gas structure in the CGM of &#39;red and dead&#39; elliptical systems, as traced by neutral HI and MgII. Halos can host tens of thousands of discrete absorbing cloudlets, with sizes of order a kpc or smaller. With a Lagrangian tracer analysis, we show that cold clouds form due to strong drho/rho >> 1 gas density perturbations which stimulate thermal instability. These local overdensities trigger rapid cooling from the hot virialized background medium at ~10^7 K to radiatively inefficient ~10^4 K clouds, which act as cosmologically long-lived, &#39;stimulated cooling&#39; seeds in a regime where the global halo does not satisfy the classic tcool/tff < 10 criterion. Furthermore, these small clouds are dominated by magnetic rather than thermal pressure, with plasma beta << 1, suggesting that magnetic fields may play an important role. The number and total mass of cold clouds both increase with resolution, and the ~8x10^4 Msun cell mass of TNG50 enables the ~few hundred pc, small-scale CGM structure we observe to form. Finally, we make a preliminary comparison against observations from the COS-LRG, LRG-RDR, COS-Halos, and SDSS LRG surveys. We broadly find that our recent, high-resolution cosmological simulations produce sufficiently high covering fractions of extended, cold gas as observed to surround massive galaxies.

preprint2020arXiv

Separate Universe Simulations with IllustrisTNG: baryonic effects on power spectrum responses and higher-order statistics

We measure power spectrum response functions in the presence of baryonic physical processes using separate universe simulations with the IllustrisTNG galaxy formation model. The response functions describe how the small-scale power spectrum reacts to long-wavelength perturbations and they can be efficiently measured with the separate universe technique by absorbing the effects of the long modes into a modified cosmology. Specifically, we focus on the total first-order matter power spectrum response to an isotropic density fluctuation $R_1(k,z)$, which is fully determined by the logarithmic derivative of the nonlinear matter power spectrum ${\rm dln}P_m(k,z)/{\rm dln}k$ and the growth-only response function $G_1(k,z)$. We find that $G_1(k,z)$ is not affected by the baryonic physical processes in the simulations at redshifts $z < 3$ and on all scales probed ($k \lesssim 15h/{\rm Mpc}$, i.e. length scales $\gtrsim 0.4 {\rm Mpc}/h$). In practice, this implies that the power spectrum fully specifies the baryonic dependence of its response function. Assuming an idealized lensing survey setup, we evaluate numerically the baryonic impact on the squeezed-lensing bispectrum and the lensing super-sample power spectrum covariance, which are given in terms of responses. Our results show that these higher-order lensing statistics can display varying levels of sensitivity to baryonic effects compared to the power spectrum, with the squeezed-bispectrum being the least sensitive. We also show that ignoring baryonic effects on lensing covariances slightly overestimates the error budget (and is therefore conservative from the point of view of parameter error bars) and likely has negligible impact on parameter biases in inference analyses.

preprint2020arXiv

Simulating the &#34;hidden giant&#34; in cold and self-interacting dark matter models

We perform a series of controlled N-body simulations to study realizations of the recently discovered Antlia 2 galaxy in cold dark matter (CDM) and self-interacting dark matter (SIDM) scenarios. Our simulations contain six benchmark models, where we vary the initial halo concentration and the self-scattering cross section. We adopt well-motivated initial stellar and halo masses, and our fiducial orbit has a small pericenter. After evolving in the Milky Way&#39;s tidal field, the simulated galaxies experience significant mass loss and their stellar distributions expand accordingly. These tidal effects are more prominent if the initial halo concentration is lower and if the self-scattering cross section is larger. Our results show that Antlia 2-like galaxies could be realized in CDM if the halo concentration is low and the stellar distribution is diffuse at the infall time, while these conditions could be relaxed in SIDM. We also find all the simulated galaxies predict approximately the same stellar velocity dispersion after imposing selection criteria for stellar particles. This has important implications for testing dark matter models using tidally disturbed systems.

preprint2020arXiv

The Angular Momentum of the Circumgalactic Medium in the TNG100 Simulation

We present an analysis of the angular momentum content of the circumgalactic medium (CGM) using TNG100, one of the flagship runs of the IllustrisTNG project. We focus on Milky Way-mass halos ($\sim 10^{12} \; M_{\odot}$) at $z=0$ but also analyze other masses and redshifts up to $z=5$. We find that the CGM angular momentum properties are strongly correlated with the stellar angular momentum of the corresponding galaxy: the CGM surrounding high-angular momentum galaxies has a systematically higher angular momentum and is better aligned to the rotational axis of the galaxy itself than the CGM surrounding low-angular momentum galaxies. Both the hot and cold phases of the CGM show this dichotomy, though it is stronger for colder gas. The CGM of high-angular momentum galaxies is characterized by a large wedge of cold gas with rotational velocities at least $\sim1/2$ of the halo&#39;s virial velocity, extending out to $\sim 1/2$ of the virial radius, and by biconical polar regions dominated by radial velocities suggestive of galactic fountains; both of these features are absent from the CGM of low-angular momentum galaxies. These conclusions are general to halo masses $\lesssim 10^{12} \; M_{\odot}$ and for $z \lesssim 2$, but they do not apply for more massive halos or at the highest redshift studied. By comparing simulations run with alterations to the fiducial feedback model, we identify the better alignment of the CGM to high-angular momentum galaxies as a feedback-independent effect and the galactic winds as a dominant influence on the CGM&#39;s angular momentum.

preprint2020arXiv

The Diversity and Variability of Star Formation Histories in Models of Galaxy Evolution

Understanding the variability of galaxy star formation histories (SFHs) across a range of timescales provides insight into the underlying physical processes that regulate star formation within galaxies. We compile the SFHs of galaxies at $z=0$ from an extensive set of models, ranging from cosmological hydrodynamical simulations (Illustris, IllustrisTNG, Mufasa, Simba, EAGLE), zoom simulations (FIRE-2, g14, and Marvel/Justice League), semi-analytic models (Santa Cruz SAM) and empirical models (UniverseMachine), and quantify the variability of these SFHs on different timescales using the power spectral density (PSD) formalism. We find that the PSDs are well described by broken power-laws, and variability on long timescales ($\gtrsim1$ Gyr) accounts for most of the power in galaxy SFHs. Most hydrodynamical models show increased variability on shorter timescales ($\lesssim300$ Myr) with decreasing stellar mass. Quenching can induce $\sim0.4-1$ dex of additional power on timescales $>1$ Gyr. The dark matter accretion histories of galaxies have remarkably self-similar PSDs and are coherent with the in-situ star formation on timescales $>3$ Gyr. There is considerable diversity among the different models in their (i) power due to SFR variability at a given timescale, (ii) amount of correlation with adjacent timescales (PSD slope), (iii) evolution of median PSDs with stellar mass, and (iv) presence and locations of breaks in the PSDs. The PSD framework is a useful space to study the SFHs of galaxies since model predictions vary widely. Observational constraints in this space will help constrain the relative strengths of the physical processes responsible for this variability.

preprint2020arXiv

The effect of differential accretion on the Gravitational Wave Background and the present day MBH Binary population

Massive black hole binaries (MBHBs) form as a consequence of galaxy mergers. However, it is still unclear whether they typically merge within a Hubble time, and how accretion may affect their evolution. These questions will be addressed by pulsar timing arrays (PTAs), which aim to detect the GW background (GWB) emitted by MBHBs during the last Myrs of inspiral. Here we investigate the influence of differential accretion on MBHB merger rates, chirp masses and the resulting GWB spectrum. We evolve a MBHB sample from the Illustris hydrodynamic cosmological simulation using semi-analytic models and for the first time self-consistently evolve their masses with binary accretion models. In all models, MBHBs coalesce with median total masses up to $1.5 \times 10^8 M_{\odot}$, up to $3-4$ times larger than in models neglecting accretion. In our model with the largest plausible impact, the median mass ratio of coalescing MBHBs increases by a factor $3.6$, the coalescence rate by $52.3\%$, and the GWB amplitude by a factor $4.0$, yielding a dimensionless GWB strain $A_{yr^{-1}} = 1 \times 10^{-15}$. Our model that favours accretion onto the primary MBH reduces the median mass ratio of coalescing MBHBs by a factor of $2.9$, and yields a GWB amplitude $A_{yr^{-1}} = 3.1 \times 10^{-16}$. This is nearly indistinguishable from our model neglecting accretion, despite higher MBHB masses at coalescence. \textbf{We further predict binary separation and mass ratio distributions of stalled MBHBs in the low-redshift universe, and find that these depend sensitively on binary accretion models. This presents the potential for combined EM and GW observational constraints on merger rates and accretion models of MBHB populations.}

preprint2020arXiv

The fate of disk galaxies in IllustrisTNG clusters

We study the stellar morphological evolution of disc galaxies within clusters in the TNG50 and TNG100 runs from the IllustrisTNG simulation suite. We select satellites of masses $10^{9.7} \leq M_{*,z=0}/\text{M}_{\odot} \leq 10^{11.6}$ residing in clusters of masses $10^{14} \lesssim M_{\text{200c,z=0}}/\text{M}_{\odot} \leq 10^{14.6}$ at $z=0$ and that were discs at accretion according to a kinematic morphology indicator (the circularity fraction). These are traced from the time of accretion to $z=0$ and compared to a control sample of central galaxies mass-matched at accretion. Most cluster discs become non-discy by $z=0$, in stark contrast with the control discs, of which a significant fraction remains discy over the same timescales. Cluster discs become non-discy accompanied by gas removal and star formation quenching, loss of dark matter and little growth or a loss of stellar mass. In contrast, control discs transform while also losing gas mass and quenching, but growing significantly in dark matter and stellar mass. Most cluster satellites change morphologies on similar timescales regardless of stellar mass, in $\sim0.5-4$ Gyr after accretion. Cluster discs that experienced more numerous and closer pericentric passages show the largest change in morphology. Morphological change in all cases requires the presence of a gravitational perturbation to drive stellar orbits to non-discy configurations, along with gas removal/heating to prevent replenishment of the disc through continued star-formation. For cluster discs, the perturbation is impulsive tidal shocking at pericentres and not tidal stripping of outer disc stellar material, whereas for control discs, a combination of mergers and AGN feedback appears to be the key driving force behind morphological transformations.

preprint2020arXiv

The Formation History of Subhalos and the Evolution of Satellite Galaxies

Satellites constitute an important fraction of the overall galaxy population and are believed to form in dark matter subhalos. Here we use the cosmological hydrodynamic simulation TNG100 to investigate how the formation histories of subhalos affect the properties and evolution of their host galaxies. We use a scaled formation time ($a_{\rm nf}$) to characterize the mass assembly histories of the subhalos before they are accreted by massive host halos. We find that satellite galaxies in young subhalos (low $a_{\rm nf}$) are less massive and more gas rich, and have stronger star formation and a higher fraction of ex situ stellar mass than satellites in old subhalos (high $a_{\rm nf}$). Furthermore, these low $a_{\rm nf}$ satellites require longer timescales to be quenched as a population than the high $a_{\rm nf}$ counterparts. We find very different merger histories between satellites in fast accretion (FA, $a_{\rm nf}<1.3$) and slow accretion (SA, $a_{\rm nf}>1.3$) subhalos. For FA satellites, the galaxy merger frequency dramatically increases just after accretion, which enhances the star formation at accretion. While, for SA satellites, the mergers occur smoothly and continuously across the accretion time. Moreover, mergers with FA satellites happen mainly after accretion, while a contrary trend is found for SA satellites. Our results provide insight into the evolution and star formation quenching of the satellite population.

preprint2020arXiv

The Formation of Ultra-Diffuse Galaxies in Clusters

We use the IllustrisTNG cosmological hydrodynamical simulation to study the formation of ultra-diffuse galaxies (UDGs) in galaxy clusters. We supplement the simulations with a realistic mass-size relation for galaxies at the time of infall into the cluster, as well as an analytical model to describe the tidally-induced evolution of their stellar mass, velocity dispersion and size. The model assumes &#34;cuspy&#34; NFW halos and, contrary to recent claims, has no difficulty reproducing the observed number of UDGs in clusters. Our results further suggest that the UDG population consists of a mixture of &#34;normal&#34; low surface brightness galaxies such as those found in the field (&#34;born&#34; UDGs, or B-UDGs), as well as a distinct population that owe their large size and low surface brightness to the effects of cluster tides (&#34;tidal&#34;, or T-UDGs). The simulations indicate that T-UDGs entered the cluster earlier and should be more prevalent than B-UDGs near the cluster centres. T-UDGs should also have, at given stellar mass, lower velocity dispersion, higher metallicities, and lower dark matter content than B-UDGs. Our results suggest that systems like DF-44 are consistent with having been born as UDGs, while others such as DF2, DF4 and VLSB-D are possibly extreme T-UDG examples.

preprint2020arXiv

The relationship between fine galaxy stellar morphology and star formation activity in cosmological simulations: a deep learning view

Hydrodynamical simulations of galaxy formation and evolution attempt to fully model the physics that shapes galaxies. The agreement between the morphology of simulated and real galaxies, and the way the morphological types are distributed across galaxy scaling relations are important probes of our knowledge of galaxy formation physics. Here we propose an unsupervised deep learning approach to perform a stringent test of the fine morphological structure of galaxies coming from the Illustris and IllustrisTNG (TNG100 and TNG50) simulations against observations from a subsample of the Sloan Digital Sky Survey. Our framework is based on PixelCNN, an autoregressive model for image generation with an explicit likelihood. We adopt a strategy that combines the output of two PixelCNN networks in a metric that isolates the fine morphological details of galaxies from the sky background. We are able to \emph{quantitatively} identify the improvements of IllustrisTNG, particularly in the high-resolution TNG50 run, over the original Illustris. However, we find that the fine details of galaxy structure are still different between observed and simulated galaxies. This difference is driven by small, more spheroidal, and quenched galaxies which are globally less accurate regardless of resolution and which have experienced little improvement between the three simulations explored. We speculate that this disagreement, that is less severe for quenched disky galaxies, may stem from a still too coarse numerical resolution, which struggles to properly capture the inner, dense regions of quenched spheroidal galaxies.

preprint2020arXiv

Tidally induced warps of spiral galaxies in IllustrisTNG

Warps are common features in both stellar and gaseous disks of nearby spiral galaxies with the latter usually easier to detect. Several theories have been proposed in the literature to explain their formation and prevalence, including tidal interactions with external galaxies. Observational correlations also suggest the importance of tides for warp formation. Here, we use the TNG100 run from the magnetohydrodynamical cosmological simulation suite IllustrisTNG to investigate the connection between interactions and the formation of gas warps. We find that in the sample of well-resolved gas-rich spiral galaxies ($10^{10}\lesssim\mathrm{M_{*}/M_{\odot}}\lesssim10^{11}$ at $z=0$) from the simulation TNG100-1, about $16\%$ possess the characteristic S-shaped warp. Around one third of these objects have their vertical morphology induced by interactions with other galaxies. Half of these interactions end with the perturber absorbed by the host. Warps formed in interactions are more asymmetrical than the remaining sample, however after the interaction the asymmetry decreases with time. We find that warps induced by interactions survive on average for $<1$ Gyr. The angle between the orbital angular momentum of the perturber and the angular momentum of the host&#39;s disk that most likely leads to warp formation is around 45 degrees. While our main goal is to investigate tidally induced warps, we find that during interactions in addition to tides, new gas that is accreted from infalling satellites also can contribute to warp formation.

preprint2020arXiv

X-ray Signatures of Black Hole Feedback: Hot Galactic Atmospheres in IllustrisTNG and X-ray Observations

Hot gaseous atmospheres that permeate galaxies and extend far beyond their stellar distribution, where they are commonly referred to as the circumgalactic medium (CGM), imprint important information about feedback processes powered by the stellar populations of galaxies and their central supermassive black holes (SMBH). In this work we study the properties of this hot X-ray emitting medium using the IllustrisTNG cosmological simulations. We analyse their mock X-ray spectra, obtained from the diffuse and metal-enriched gas in TNG100 and TNG50, and compare the results with X-ray observations of nearby early-type galaxies. The simulations reproduce the observed X-ray luminosities ($L_{\rm X}$) and temperature ($T_{\rm X})$ at small ($<R_{\rm e}$) and intermediate ($<5R_{\rm e}$) radii reasonably well. We find that the X-ray properties of lower mass galaxies depend on their star formation rates. In particular, in the magnitude range where the star-forming and quenched populations overlap, $M_{\rm K}\sim-24$ $ (M_*\sim10^{10.7}M_\odot)$, we find that the X-ray luminosities of star-forming galaxies are on average about an order of magnitude higher than those of their quenched counterparts. We show that this diversity in $L_{\rm X}$ is a direct manifestation of the quenching mechanism in the simulations, where the galaxies are quenched due to gas expulsion driven by SMBH kinetic feedback. The observed dichotomy in $L_{\rm X}$ is thus an important observable prediction for the SMBH feedback-based quenching mechanisms implemented in state-of-the-art cosmological simulations. While the current X-ray observations of star forming galaxies are broadly consistent with the predictions of the simulations, the observed samples are small and more decisive tests are expected from the sensitive all-sky X-ray survey with eROSITA.

preprint2019arXiv

Antlia2&#39;s role in driving the ripples in the outer gas disk of the Galaxy

We employ the earlier published proper motions of the newly discovered Antlia 2 dwarf galaxy derived from Gaia data to calculate its orbital distribution in the cosmologically recent past. Using these observationally motivated orbits, we calculate the effect of the Antlia 2 dwarf galaxy on the outer HI disk of the Milky Way, using both test particle and Smoothed Particle Hydrodynamics simulations. We find that orbits with low pericenters, $\sim$ 10 kpc, produce disturbances that match the observed outer HI disk perturbations. We have independently recalculated the proper motion of the Antlia 2 dwarf from Gaia data and found a proper motion of $(μ_αcosδ, μ_δ) = (-0.068,0.032) \pm (0.023,-0.031)~\rm mas/yr$, which agrees with results from Torrealba et al. (2019) within the errors, but gives lower mean pericenters, e.g., $\sim$ 15 kpc for our fiducial model of the Milky Way. We also show that the Sagittarius dwarf galaxy interaction does not match the observed perturbations in the outer gas disk. Thus, Antlia 2 may be the driver of the observed large perturbations in the outer gas disk of the Galaxy. The current location of the Antlia 2 dwarf galaxy closely matches that predicted by an earlier dynamical analysis (Chakrabarti \& Blitz 2009) of the dwarf that drove ripples in the outer Galaxy, and, in particular, its orbit is nearly coplanar to the Galactic disk. If the Antlia 2 dwarf galaxy is responsible for the perturbations in the outer Galactic disk, it would have a specific range of proper motions that we predict here; this can be tested soon with Gaia DR-3 and Gaia DR-4 data.

preprint2019arXiv

Baryons in the CosmicWeb of IllustrisTNG -- II: the Connection among Galaxies, Halos, their Formation Time and their Location in the Cosmic Web

The connections among galaxies, the dark matter halos where they form and the properties of the large-scale Cosmic Web still need to be completely disentangled. We use the cosmological hydrodynamical simulation TNG100 of the IllustrisTNG suite to quantify the effects played by the large-scale density field and the Cosmic Web morphology on the relation between halo mass and galaxy stellar mass. We select objects with total dynamical mass in the range $\geq 6.3\times 10^{10} h ^{-1}\, M_{\odot}$ up to a few $10^{14} h^{-1} \, M_{\odot}$ between redshift $z=4$ and redshift $z=0$. A Cosmic Web class (knot, filament, sheet, void) is assigned to each region of the volume using a density field deformation tensor-based method. We find that galaxy stellar mass strongly correlates with total dynamical mass and formation time, and more weakly with large-scale overdensity and Cosmic Web class. The latter two quantities correlate with each other, but are not entirely degenerate. Furthermore, we find that at fixed halo mass, galaxies with stellar mass lower than the median value are more likely to be found in voids and sheets, whereas galaxies with stellar mass higher than the median are more likely to be found in filaments and knots. Finally, we find that the dependence on environment is stronger for satellites than for centrals, and discuss the physical implications of these results.

preprint2019arXiv

Correlations Between Black Holes and Host Galaxies in the Illustris and IllustrisTNG Simulations

We study black hole - host galaxy correlations, and the relation between the over-massiveness (the distance from the average $M_{BH}-σ$ relation) of super-massive black holes (SMBHs) and star formation histories of their host galaxies in the Illustris and TNG100 simulations. We find that both simulations are able to produce black hole scaling relations in general agreement with observations at $z=0$, but with noticeable discrepancies. Both simulations show an offset from the observations for the $M_{BH}-σ$ relation, and the relation between $M_{BH}$ and the Sersic index. The relation between $M_{BH}$ and stellar mass $M_*$ is tighter than the observations, especially for TNG100. For massive galaxies in both simulations, the hosts of over-massive SMBHs (those above the mean $M_{BH}-σ$ relation) tend to have larger Sersic indices and lower baryon conversion efficiency, suggesting a multidimensional link between SMBHs and properties of their hosts. In Illustris, the hosts of over-massive SMBHs have formed earlier and have lower present-day star formation rates, in qualitative agreement with the observations for massive galaxies with $σ>100 \rm km/s$. For low-mass galaxies, such a correlation still holds in Illustris but does not exist in the observed data. For TNG100, the correlation between SMBH over-massiveness and star formation history is much weaker. The hosts of over-massive SMBHs generally have consistently larger star formation rates throughout history. These galaxies have higher stellar mass as well, due to the strong $M_{BH}-M_*$ correlation. Our findings show that simulated SMBH scaling relations and correlations are sensitive to features in the modeling of SMBHs.

preprint2019arXiv

Galaxy Formation with BECDM -- II. Cosmic Filaments and First Galaxies

Bose-Einstein Condensate Dark Matter (BECDM; also known as Fuzzy Dark Matter) is motivated by fundamental physics and has recently received significant attention as a serious alternative to the established Cold Dark Matter (CDM) model. We perform cosmological simulations of BECDM gravitationally coupled to baryons and investigate structure formation at high redshifts ($z \gtrsim 5$) for a boson mass $m=2.5\cdot 10^{-22}~{\rm eV}$, exploring the dynamical effects of its wavelike nature on the cosmic web and the formation of first galaxies. Our BECDM simulations are directly compared to CDM as well as to simulations where the dynamical quantum potential is ignored and only the initial suppression of the power spectrum is considered -- a Warm Dark Matter-like (&#34;WDM&#34;) model often used as a proxy for BECDM. Our simulations confirm that &#34;WDM&#34; is a good approximation to BECDM on large cosmological scales even in the presence of the baryonic feedback. Similarities also exist on small scales, with primordial star formation happening both in isolated haloes and continuously along cosmic filaments; the latter effect is not present in CDM. Global star formation and metal enrichment in these first galaxies are delayed in BECDM/&#34;WDM&#34; compared to the CDM case: in BECDM/&#34;WDM&#34; first stars form at $z\sim 13$/$13.5$ while in CDM star formation starts at $z\sim 35$. The signature of BECDM interference, not present in &#34;WDM&#34;, is seen in the evolved dark matter power spectrum: although the small scale structure is initially suppressed, power on kpc scales is added at lower redshifts. Our simulations lay the groundwork for realistic simulations of galaxy formation in BECDM.

preprint2019arXiv

Galaxy interactions in IllustrisTNG-100, I: The power and limitations of visual identification

We present a sample of 446 galaxy pairs constructed using the cosmological simulation IllustrisTNG-100 at z = 0, with M$_{FoF, dm}$ = 10$^{11}$-10$^{13.5}$ M$_{\odot}$. We produce ideal mock SDSS g-band images of all pairs to test the reliability of visual classification schema employed to produce samples of interacting galaxies. We visually classify each image as interacting or not based on the presence of a close neighbour, the presence of stellar debris fields, disturbed discs, and/or tidal features. By inspecting the trajectories of the pairs, we determine that these indicators correctly identify interacting galaxies $\sim45\%$ of the time. We subsequently split the sample into the visually identified interacting pairs (VIP; 38 pairs) and those which are interacting but are not visually identified (nonVIP; 47 pairs). We find that VIP have undergone a close passage nearly twice as recently as the nonVIP, and typically have higher stellar masses. Further, the VIP sit in dark matter haloes that are approximately 2.5 times as massive, in environments nearly 2 times as dense, and are almost a factor of 10 more affected by the tidal forces of their surroundings than the nonVIP. These factors conspire to increase the observability of tidal features and disturbed morphologies, making the VIP more likely to be identified. Thus, merger rate calculations which rely on stellar morphologies are likely to be significantly biased toward massive galaxy pairs which have recently undergone a close passage.

preprint2019arXiv

High redshift JWST predictions from IllustrisTNG: Dust modelling and galaxy luminosity functions

The James Webb Space Telescop (JWST) promises to revolutionise our understanding of the early Universe, and contrasting its upcoming observations with predictions of the $Λ$CDM model requires detailed theoretical forecasts. Here, we exploit the large dynamic range of the IllustrisTNG simulation suite, TNG50, TNG100, and TNG300, to derive multi-band galaxy luminosity functions from $z=2$ to $z=10$. We put particular emphasis on the exploration of different dust attenuation models to determine galaxy luminosity functions for the rest-frame ultraviolet (UV), and apparent wide NIRCam bands. Our most detailed dust model is based on continuum Monte Carlo radiative transfer calculations employing observationally calibrated dust properties. This calibration results in constraints on the redshift evolution of the dust attenuation normalisation and dust-to-metal ratios yielding a stronger redshift evolution of the attenuation normalisation compared to most previous theoretical studies. Overall we find good agreement between the rest-frame UV luminosity functions and observational data for all redshifts, also beyond the regimes used for the dust-model calibrations. Furthermore, we also recover the observed high redshift ($z=4-6$) UV luminosity versus stellar mass relation, the H$α$ versus star formation rate relation, and the H$α$ luminosity function at $z=2$. The bright end ($M_{\rm UV}>-19.5$) cumulative galaxy number densities are consistent with observational data. For the F200W NIRCam band, we predict that JWST will detect $\sim 80$ ($\sim 200$) galaxies with a signal-to-noise ratio of $10$ ($\sim 5$) within the NIRCam field of view, $2.2\times2.2 \,{\rm arcmin}^{2}$, for a total exposure time of $10^5{\rm s}$ in the redshift range $z=8 \pm 0.5$. These numbers drop to $\sim 10$ ($\sim 40$) for an exposure time of $10^4{\rm s}$.

preprint2019arXiv

The relationship between black hole mass and galaxy properties: Examining the black hole feedback model in IllustrisTNG

Supermassive black hole feedback is thought to be responsible for the lack of star formation, or quiescence, in a significant fraction of galaxies. We explore how observable correlations between the specific star formation rate (sSFR), stellar mass (M$_{\rm{star}}$), and black hole mass (M$_{\rm{BH}}$) are sensitive to the physics of black hole feedback in a galaxy formation model. We use the IllustrisTNG simulation suite, specifically the TNG100 simulation and ten model variations that alter the parameters of the black hole model. Focusing on central galaxies at $z = 0$ with M$_{\rm{star}} > 10^{10}$ M$_{\odot}$, we find that the sSFR of galaxies in IllustrisTNG decreases once the energy from black hole kinetic winds at low accretion rates becomes larger than the gravitational binding energy of gas within the galaxy stellar radius. This occurs at a particular M$_{\rm{BH}}$ threshold above which galaxies are found to sharply transition from being mostly star-forming to mostly quiescent. As a result of this behavior, the fraction of quiescent galaxies as a function of M$_{\rm{star}}$ is sensitive to both the normalization of the M$_{\rm{BH}}$-M$_{\rm{star}}$ relation and the M$_{\rm{BH}}$ threshold for quiescence in IllustrisTNG. Finally, we compare these model results to observations of 91 central galaxies with dynamical M$_{\rm{BH}}$ measurements with the caveat that this sample is not representative of the whole galaxy population. While IllustrisTNG reproduces the observed trend that quiescent galaxies host more massive black holes, the observations exhibit a broader scatter in M$_{\rm{BH}}$ at a given M$_{\rm{star}}$ and show a smoother decline in sSFR with M$_{\rm{BH}}$.

preprint2018arXiv

Efficacy of early stellar feedback in low gas surface density environments

We present a suite of high resolution radiation hydrodynamic simulations of a small patch ($1 \ {\rm kpc}^2$) of the inter-stellar medium (ISM) performed with Arepo-RT, with the aim to quantify the efficacy of various feedback processes like supernovae explosions (SNe), photoheating and radiation pressure in low gas surface density galaxies ($Σ_{\rm gas} \simeq 10 \ {\rm M}_\odot \ {\rm pc}^{-2}$). We show that radiation fields decrease the star formation rate and therefore the total stellar mass formed by a factor of $\sim 2$. This increases the gas depletion timescale and brings the simulated Kennicutt-Schmidt relation closer to the observational estimates. Radiation feedback coupled with SNe is more efficient at driving outflows with the mass and energy loading increasing by a factor of $\sim 10$. This increase is mainly driven by the additional entrainment of medium density ($10^{-2} \leq n< 1 \ {\rm cm}^{-3}$), warm ($300 \ {\rm K}\leq T<8000 \ {\rm K}$) material. Therefore including radiation fields tends to launch colder, denser and higher mass and energy loaded outflows. This is because photoheating of the high density gas around a newly formed star over-pressurises the region, causing it to expand. This reduces the ambient density in which the SNe explode by a factor of $10-100$ which in turn increases their momentum output by a factor of $\sim 1.5-2.5$. Finally, we note that in these low gas surface density environments, radiation fields primarily impact the ISM via photoheating and radiation pressure has only a minimal role in regulating star formation.