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. In this work, we employ a relatively simple method based on previous work to predict the submm flux densities of simulated SMGs drawn from cosmological simulations from the Illustris and IllustrisTNG projects 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. We demonstrate that IllustrisTNG hosts fewer SMGs than Illustris because in the former, high-mass ($M_{\star} \sim 10^{11} \, {\rm M}_{\odot}$) $z \sim 2-3$ galaxies have lower dust masses and star formation rates (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. Higher resolution, which would lead to stronger SFR enhancements in starbursts, could at least partially reconcile the discrepancy between the IllustrisTNG SMG number counts and those observed.

During their violent merger, two neutron stars can shed a few percent of their mass. As this ejecta expands, it collides with the surrounding interstellar gas, producing a slowly-fading {\it radio flare} that lasts for years. Radio flares uniquely probe the neutron star merger populations as many events from past decades could still be detectable. Nonetheless, no radio flare observation has been reported to date. Here we show that the radio transient FIRST J1419+3940, first observed in 1993 and still detectable, likely originated from a neutron star merger. We carry out numerical simulations of neutron star merger ejecta to demonstrate that the observed radio light curve is well reproduced by a merger model with astrophysically expected parameters. We find that a neutron star merger model is a priori more likely and provides a better fit to the data than the alternative explanation---the afterglow of an off-axis long gamma-ray burst. Near-future observations could unambiguously establish to the FIRST J1419+3940 radio transient's origin. We show that existing radio surveys likely already recorded multiple radio flares, informing us of the origin and properties of neutrons tar mergers and their role in the nucleosynthesis of the heaviest elements in the Universe.

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}_{\alpha}$ and ${\rm H}_{\beta}$ + $[\rm O \,III]$ luminosity functions up to $z=8$. The predicted ${\rm H}_{\beta}$ + $[\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}_{\alpha}$ 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$\AA (D4000) and the UV continuum slope $\beta$. 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 $\beta$ 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.

We use the IllustrisTNG cosmological hydrodynamical simulations to investigate how the specific star formation rates (sSFRs) of massive galaxies $(M_* > 10^{10} M_\odot)$ depend on the distance to their closest companions. We estimate sSFR enhancements by comparing with control samples that are matched in redshift, stellar mass, local density and isolation, and we restrict our analysis to pairs with stellar mass ratios of 0.1 to 10. At small separations (~15 kpc), the mean sSFR is enhanced by a factor of 2.0 $\pm$ 0.1 in the flagship (110.7 Mpc)$^3$ simulation (TNG100-1). Statistically significant enhancements extend out to 3D separations of 280 kpc in the (302.6 Mpc)$^3$ simulation (TNG300-1). We find similar trends in the EAGLE and Illustris simulations, although their sSFR enhancements are lower than those in TNG100-1 by about a factor of two. Enhancements in IllustrisTNG galaxies are seen throughout the redshift range explored ($0 \leq z < 1$), with the strength of the enhancements decreasing with increasing redshift for galaxies with close companions. In order to more closely compare with observational results, we separately consider 2D projected distances between galaxies in IllustrisTNG. We detect significant sSFR enhancements out to projected separations of 260 kpc in TNG300-1, with projection effects diluting the size of the enhancements by about 20 per cent below 50 kpc. We find similar sSFR enhancements in TNG100-1 and Sloan Digital Sky Survey galaxies, with enhancements extending out to projected separations of about 150 kpc for star-forming galaxies at $z < 0.2$. Finally, by summing over all separations, we estimate that the presence of closest companions boosts the average sSFR of massive galaxies in TNG100-1 by 14.5 per cent.

Galaxy mergers are a major evolutionary transformation whose effects are borne out by a plethora of observations and numerical simulations. However, most previous simulations have used idealised, isolated, binary mergers and there has not been significant progress on studying statistical samples of galaxy mergers in large cosmological simulations. We present a sample of 27,691 post-merger (PM) galaxies ($0\le z \le 1$) identified from IllustrisTNG: a cosmological, large box, magneto-hydrodynamical simulation suite. The PM sample spans a wide range of merger and galaxy properties ($M_\star$, $\mu$, $f_\mathrm{gas}$). We demonstrate that star forming (SF) PMs exhibit enhanced star formation rates (SFRs) on average by a factor of $\sim 2$, while the passive PMs show no statistical enhancement. We find that the SFR enhancements: (1) show no dependence on redshift, (2) anti-correlate with the PM's stellar mass, and (3) correlate with the gas fraction of the PM's progenitors. However, SF PMs show stronger enhancements which may indicate other processes being at play (e.g., gas phase, feedback efficiency). Although the SFR enhancement correlates mildly with the merger mass ratio, the more abundant minor mergers ($0.1 \le \mu < 0.3 $) still contribute $\sim 50\%$ of the total SFR enhancement. By tracing the PM sample forward in time, we find that galaxy mergers can drive significant SFR enhancements which decay over $\sim 0.5$ Gyr independent of the merger mass ratio, although the decay timescale is dependent on the simulation resolution. The strongest merger-driven starburst galaxies evolve to be passive/quenched on faster timescales than their controls.

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.

The correlation between galaxies' integrated stellar masses and star formation rates (the `star formation main sequence'; SFMS) is a well-established scaling relation. Recently, surveys have found a relationship between the star formation rate and stellar mass surface densities on kpc and sub-kpc scales (the `resolved SFMS'; rSFMS). In this work, we demonstrate that the rSFMS emerges naturally in FIRE-2 zoom-in simulations of Milky Way-mass galaxies. We make SFR and stellar mass maps of the simulated galaxies at a variety of spatial resolutions and star formation averaging time-scales and fit the rSFMS using multiple methods from the literature. While the absolute value of the SFMS slope depends on the fitting method, the slope is steeper for longer star formation time-scales and lower spatial resolutions regardless of the fitting method employed. We present a toy model that quantitatively captures the dependence of the simulated galaxies' rSFMS slope on spatial resolution and use it to illustrate how this dependence can be used to constrain the characteristic mass of star-forming clumps.

Recent cosmological hydrodynamical simulations are able to reproduce numerous statistical properties of galaxies that are consistent with observational data. Yet, the adopted subgrid models strongly affect the simulation outcomes, limiting the predictive power of these simulations. In this work, we perform a suite of isolated galactic disk simulations under the {\it SMUGGLE} framework and investigate how different subgrid models affect the properties of giant molecular clouds (GMCs). We employ {\sc astrodendro}, a hierarchical clump-finding algorithm, to identify GMCs in the simulations. We find that different choices of subgrid star formation efficiency, $\epsilon_{\rm ff}$, and stellar feedback channels, yield dramatically different mass and spatial distributions for the GMC populations. Without feedback, the mass function of GMCs has a shallower power-law slope and extends to higher mass ranges compared to runs with feedback. Moreover, higher $\epsilon_{\rm ff}$ results in faster molecular gas consumption and steeper mass function slopes. Feedback also suppresses power in the two-point correlation function (TPCF) of the spatial distribution of GMCs. Specifically, radiative feedback strongly reduces the TPCF on scales below 0.2~kpc, while supernova feedback reduces power on scales above 0.2~kpc. Finally, runs with higher $\epsilon_{\rm ff}$ exhibit a higher TPCF than runs with lower $\epsilon_{\rm ff}$, because the dense gas is depleted more efficiently thereby facilitating the formation of well-structured supernova bubbles. We argue that comparing simulated and observed GMC populations can help better constrain subgrid models in the next-generation of galaxy formation simulations.

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\gamma^{\prime}\rangle = 2.011 \pm 0.007$ and a scatter of $\sigma_{\gamma^{\prime}} = 0.171$ over the radial range 0.4 to 4 times the stellar half mass radius. Several correlations between $\gamma^{\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, $\gamma^{\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 $\gamma^{\prime}$ and matching observed galaxy correlations. The effects of stellar winds on $\gamma^{\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.

Over the last decades, cosmological simulations of galaxy formation have been instrumental for advancing our understanding of structure and galaxy formation in the Universe. These simulations follow the non-linear evolution of galaxies modeling a variety of physical processes over an enormous range of scales. A better understanding of the physics relevant for shaping galaxies, improved numerical methods, and increased computing power have led to simulations that can reproduce a large number of observed galaxy properties. Modern simulations model dark matter, dark energy, and ordinary matter in an expanding space-time starting from well-defined initial conditions. The modeling of ordinary matter is most challenging due to the large array of physical processes affecting this matter component. Cosmological simulations have also proven useful to study alternative cosmological models and their impact on the galaxy population. This review presents a concise overview of the methodology of cosmological simulations of galaxy formation and their different applications.

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.

We present a model for the interaction between dust and radiation fields in the radiation hydrodynamic code AREPO-RT, which solves the moment-based radiative transfer equations on an unstructured moving mesh. Dust is directly treated using live simulation particles, each of which represent a population of grains that are coupled to hydrodynamic motion through a drag force. We introduce methods to calculate radiation pressure on and photon absorption by dust grains. By including a direct treatment of dust, we are able to calculate dust opacities and update radiation fields self-consistently based on the local dust distribution. This hybrid scheme coupling dust particles to an unstructured mesh for radiation is validated using several test problems with known analytic solutions, including dust driven via spherically-symmetric flux from a constant luminosity source and photon absorption from radiation incident on a thin layer of dust. Our methods are compatible with the multifrequency scheme in AREPO-RT, which treats UV and optical photons as single-scattered and IR photons as multi-scattered. At IR wavelengths, we model heating of and thermal emission from dust. Dust and gas are not assumed to be in local thermodynamic equilibrium but transfer energy through collisional exchange. We estimate dust temperatures by balancing these dust-radiation and dust-gas energy exchange rates. This framework for coupling dust and radiation can be applied in future radiation hydrodynamic simulations of galaxy formation.

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.

We study the evolutionary trend of the total density profile of early-type galaxies (ETGs) in IllustrisTNG. To this end, we trace ETGs from $z=0$ to $z=4$ and measure the power-law slope $\gamma^{\prime}$ of the total density profile for their main progenitors. We find that their $\gamma^{\prime}$ steepen on average during $z\sim4-2$, then becoming shallower until $z=1$, after which they remain almost constant, aside from a residual trend of becoming shallower towards $z=0$. We also compare to a statistical sample of ETGs at different redshifts, selected based on their luminosity profiles and stellar masses. Due to different selection effects, the average slopes of the statistical samples follow a modified evolutionary trend. They monotonically decrease since $z=3$, and after $z\approx 1$, they remain nearly invariant with a mild increase towards $z=0$. These evolutionary trends are mass-dependent for both samples, with low-mass galaxies having in general steeper slopes than their more massive counterparts. Galaxies that transitioned to ETGs more recently have steeper mean slopes as they tend to be smaller and more compact at any given redshift. By analyzing the impact of mergers and AGN feedback on the progenitors' evolution, we conjecture a multi-phase path leading to isothermality in ETGs: dissipation associated with rapid wet mergers tends to steepen $\gamma^{\prime}$ from $z=4$ to $z=2$, whereas subsequent AGN feedback (especially in the kinetic mode) makes $\gamma^{\prime}$ shallower again from $z=2$ to $z=1$. Afterwards, passive evolution from $z=1$ to $z=0$, mainly through gas-poor mergers, mildly decreases $\gamma^{\prime}$ and maintains the overall mass distribution close to isothermal.

We present a novel framework to self-consistently model the effects of radiation fields, dust physics and molecular chemistry (H$_2$) in the interstellar medium (ISM) of galaxies. The model combines a state-of-the-art radiation hydrodynamics module with a non-equilibrium thermochemistry module that accounts for H$_2$ coupled to a realistic dust formation and destruction model, all integrated into the new stellar feedback framework SMUGGLE. We test this model on high-resolution isolated Milky-Way (MW) simulations. We show that photoheating from young stars makes stellar feedback more efficient, but this effect is quite modest in low gas surface density galaxies like the MW. The multi-phase structure of the ISM, however, is highly dependent on the strength of the interstellar radiation field. We are also able to predict the distribution of H$_2$, that allow us to match the molecular Kennicutt-Schmidt (KS) relation, without calibrating for it. We show that the dust distribution is a complex function of density, temperature and ionization state of the gas which cannot be reproduced by simple scaling relations often used in the literature. Our model is only able to match the observed dust temperature distribution if radiation from the old stellar population is considered, implying that these stars have a non-negligible contribution to dust heating in the ISM. Our state-of-the-art model is well-suited for performing next generation cosmological galaxy formation simulations, which will be able to predict a wide range of resolved ($\sim 10$ pc) properties of galaxies.

Machine learning is becoming a popular tool to quantify galaxy morphologies and identify mergers. However, this technique relies on using an appropriate set of training data to be successful. By combining hydrodynamical simulations, synthetic observations and convolutional neural networks (CNNs), we quantitatively assess how realistic simulated galaxy images must be in order to reliably classify mergers. Specifically, we compare the performance of CNNs trained with two types of galaxy images, stellar maps and dust-inclusive radiatively transferred images, each with three levels of observational realism: (1) no observational effects (idealized images), (2) realistic sky and point spread function (semi-realistic images), (3) insertion into a real sky image (fully realistic images). We find that networks trained on either idealized or semi-real images have poor performance when applied to survey-realistic images. In contrast, networks trained on fully realistic images achieve 87.1% classification performance. Importantly, the level of realism in the training images is much more important than whether the images included radiative transfer, or simply used the stellar maps (87.1% compared to 79.6% accuracy, respectively). Therefore, one can avoid the large computational and storage cost of running radiative transfer with a relatively modest compromise in classification performance. Making photometry-based networks insensitive to colour incurs a very mild penalty to performance with survey-realistic data (86.0% with r-only compared to 87.1% with gri). This result demonstrates that while colour can be exploited by colour-sensitive networks, it is not necessary to achieve high accuracy and so can be avoided if desired. We provide the public release of our statistical observational realism suite, RealSim, as a companion to this paper.

We present a new cosmological, magnetohydrodynamical simulation for galaxy formation: TNG50, the third and final installment of the IllustrisTNG project. TNG50 evolves 2x2160^3 dark-matter particles and gas cells in a volume 50 comoving Mpc across. It hence reaches a numerical resolution typical of zoom-in simulations, with a baryonic element mass of 8.5x10^4 Msun and an average cell size of 70-140 parsecs in the star-forming regions of galaxies. Simultaneously, TNG50 samples ~700 (6,500) galaxies with stellar masses above 10^10 (10^8) Msun at z=1. Here we investigate the structural and kinematical evolution of star-forming galaxies across cosmic time (0 < z < 6). We quantify their sizes, disk heights, 3D shapes, and degree of rotational vs. dispersion-supported motions as traced by rest-frame V-band light (i.e. roughly stellar mass) and by Halpha light (i.e. star-forming and dense gas). The unprecedented resolution of TNG50 enables us to model galaxies with sub-kpc half-light radii and with <300-pc disk heights. Coupled with the large-volume statistics, we characterize a diverse, redshift- and mass-dependent structural and kinematical morphological mix of galaxies all the way to early epochs. Our model predicts that for star-forming galaxies the fraction of disk-like morphologies, based on 3D stellar shapes, increases with both cosmic time and galaxy stellar mass. Gas kinematics reveal that the vast majority of 10^9-11.5 Msun star-forming galaxies are rotationally-supported disks for most cosmic epochs (Vmax/sigma>2-3, z<5), being dynamically hotter at earlier epochs (z>1.5). Despite large velocity dispersion at high redshift, cold and dense gas in galaxies predominantly arranges in disky or elongated shapes at all times and masses; these gaseous components exhibit rotationally-dominated motions far exceeding the collisionless stellar bodies.

We introduce the Stars and MUltiphase Gas in GaLaxiEs -- SMUGGLE model, an explicit and comprehensive stellar feedback model for the moving-mesh code arepo. This novel sub-resolution model resolves the multiphase gas structure of the interstellar medium and self-consistently generates gaseous outflows. The model implements crucial aspects of stellar feedback including photoionization, radiation pressure, energy and momentum injection from stellar winds and from supernovae. We explore this model in high-resolution isolated simulations of Milky Way-like disc galaxies. Stellar feedback regulates star formation to the observed level and naturally captures the establishment of a Kennicutt-Schmidt relation. This result is achieved independent of the numerical mass and spatial resolution of the simulations. Gaseous outflows are generated with average mass loading factors of the order of unity. Strong outflow activity is correlated with peaks in the star formation history of the galaxy with evidence that most of the ejected gas eventually rains down onto the disc in a galactic fountain flow that sustains late-time star formation. Finally, the interstellar gas in the galaxy shows a distinct multiphase distribution with a coexistence of cold, warm and hot phases.

We present the new TNG50 cosmological, magnetohydrodynamical simulation -- the third and final volume of the IllustrisTNG project. This simulation occupies a unique combination of large volume and high resolution, with a 50 Mpc box sampled by 2160^3 gas cells (baryon mass of 8x10^4 Msun). The median spatial resolution of star-forming ISM gas is ~100-140 parsecs. This resolution approaches or exceeds that of modern 'zoom' simulations of individual massive galaxies, while the volume contains ~20,000 resolved galaxies with M*>10^7 Msun. Herein we show first results from TNG50, focusing on galactic outflows driven by supernovae as well as supermassive black hole feedback. We find that the outflow mass loading is a non-monotonic function of galaxy stellar mass, turning over and rising rapidly above 10^10.5 Msun due to the action of the central black hole. Outflow velocity increases with stellar mass, and at fixed mass is faster at higher redshift. The TNG model can produce high velocity, multi-phase outflows which include cool, dense components. These outflows reach speeds in excess of 3000 km/s out to 20 kpc with an ejective, BH-driven origin. Critically, we show how the relative simplicity of model inputs (and scalings) at the injection scale produces complex behavior at galactic and halo scales. For example, despite isotropic wind launching, outflows exhibit natural collimation and an emergent bipolarity. Furthermore, galaxies above the star-forming main sequence drive faster outflows, although this correlation inverts at high mass with the onset of quenching, whereby low luminosity, slowly accreting, massive black holes drive the strongest outflows.

Simulating the dust content of galaxies and their surrounding gas is challenging due to the wide range of physical processes affecting the dust evolution. Here we present cosmological hydrodynamical simulations of a cluster of galaxies, $M_\text{200,crit}=6 \times 10^{14}\,{\rm M_\odot}$, including a novel dust model for the moving mesh code {\sc Arepo}. This model includes dust production, growth, supernova-shock-driven destruction, ion-collision-driven thermal sputtering, and high temperature dust cooling through far infrared re-radiation of collisionally deposited electron energies. Adopting a rather low thermal sputtering rate, we find, consistent with observations, a present-day overall dust-to-gas ratio of $\sim 2\times 10^{-5}$, a total dust mass of $\sim 2\times 10^9\,{\rm M_\odot}$, and a dust mass fraction of $\sim 3\times 10^{-6}$. The typical thermal sputtering timescales within $\sim 100\,{\rm kpc}$ are around $\sim 10\,{\rm Myr}$, and increase towards the outer parts of the cluster to $\sim 10^3\,{\rm Myr}$ at a cluster-centric distance of $1\,{\rm Mpc}$. The condensation of gas phase metals into dust grains reduces high temperature metal-line cooling, but also leads to additional dust infrared cooling. The additional infrared cooling changes the overall cooling rate in the outer parts of the cluster, beyond $\sim 1\,{\rm Mpc}$, by factors of a few. This results in noticeable changes of the entropy, temperature, and density profiles of cluster gas once dust formation is included. The emitted dust infrared emission due to dust cooling is consistent with observational constraints.

We present a machine-learning approach for estimating galaxy cluster masses from Chandra mock images. We utilize a Convolutional Neural Network (CNN), a deep machine learning tool commonly used in image recognition tasks. The CNN is trained and tested on our sample of 7,896 Chandra X-ray mock observations, which are based on 329 massive clusters from the IllustrisTNG simulation. Our CNN learns from a low resolution spatial distribution of photon counts and does not use spectral information. Despite our simplifying assumption to neglect spectral information, the resulting mass values estimated by the CNN exhibit small bias in comparison to the true masses of the simulated clusters (-0.02 dex) and reproduce the cluster masses with low intrinsic scatter, 8% in our best fold and 12% averaging over all. In contrast, a more standard core-excised luminosity method achieves 15-18% scatter. We interpret the results with an approach inspired by Google DeepDream and find that the CNN ignores the central regions of clusters, which are known to have high scatter with mass.

Motivated by the recently discovered kinematic "Hubble sequence" shown by the stellar orbit-circularity distribution of 260 CALIFA galaxies, we make use of a comparable galaxy sample at z = 0 with a stellar mass range from 5E9 to 5E11 solar masses, selected from the IllustrisTNG simulation and study their stellar orbit compositions in relation to a number of other fundamental galaxy properties.We find that the TNG100 simulation broadly reproduces the observed fractions of different orbital components and their stellar mass dependencies. In particular, the mean mass dependencies of the luminosity fractions for the kinematically warm and hot orbits are well reproduced within model uncertainties of the observed galaxies. The simulation also largely reproduces the observed peak and trough features at a stellar mass of 1-2E10 solar masses, in the mean distributions of the cold- and hot-orbit fractions, respectively, indicating fewer cooler orbits and more hotter orbits in both more- and less-massive galaxies beyond such a mass range. Several marginal disagreements are seen between the simulation and observations: the average cold-orbit (counter-rotating) fractions of the simulated galaxies below (above) a stellar mass of 6E10 solar masses, are systematically higher than the observational data by < 10% (absolute orbital fraction); the simulation also seems to produce more scatter for the cold-orbit fraction and less so for the non-cold orbits at any given galaxy mass. Possible causes that stem from the adopted heating mechanisms are discussed.

We analyze the IllustrisTNG simulations to study the mass, volume fraction and phase distribution of gaseous baryons embedded in the knots, filaments, sheets and voids of the Cosmic Web from redshift $z=8$ to redshift $z=0$. We find that filaments host more star-forming gas than knots, and that filaments also have a higher relative mass fraction of gas in this phase than knots. We also show that the cool, diffuse Intergalactic Medium (IGM; $T<10^5 \, {\rm K}$, $ n_{\rm H}<10^{-4}(1+z) \, {\rm cm^{-3}}$) and the Warm-Hot Intergalactic Medium (WHIM; $ 10^5 \, {\rm K} <T<10^7 \, {\rm K}$, $ n_{\rm H} <10^{-4}(1+z)\, {\rm cm^{-3}}$) constitute $\sim 39\%$ and $\sim 46\%$ of the baryons at redshift $z=0$, respectively. Our results indicate that the WHIM may constitute the largest reservoir of {\it missing} baryons at redshift $z=0$. Using our Cosmic Web classification, we predict the WHIM to be the dominant baryon mass contribution in filaments and knots at redshift $z=0$, but not in sheets and voids where the cool, diffuse IGM dominates. We also characterise the evolution of WHIM and IGM from redshift $z=4$ to redshift $z=0$, and find that the mass fraction of WHIM in filaments and knots evolves only by a factor $\sim 2$ from redshift $z=0$ to $z=1$, but declines faster at higher redshift. The WHIM only occupies $4-11\%$ of the volume at redshift $0\leq z \leq 1$. We predict the existence of a significant number of currently undetected OVII and NeIX absorption systems in cosmic filaments which could be detected by future X-ray telescopes like Athena.

The James Webb Space Telescop (JWST) promises to revolutionise our understanding of the early Universe, and contrasting its upcoming observations with predictions of the $\Lambda$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$\alpha$ versus star formation rate relation, and the H$\alpha$ 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}$.

We present image-based evolution of galaxy mergers from the Illustris cosmological simulation at 12 time-steps over 0.5 < z < 5. To do so, we created approximately one million synthetic deep Hubble Space Telescope and James Webb Space Telescope images and measured common morphological indicators. Using the merger tree, we assess methods to observationally select mergers with stellar mass ratios as low as 10:1 completing within +/- 250 Myr of the mock observation. We confirm that common one- or two-dimensional statistics select mergers so defined with low purity and completeness, leading to high statistical errors. As an alternative, we train redshift-dependent random forests (RFs) based on 5-10 inputs. Cross-validation shows the RFs yield superior, yet still imperfect, measurements of the late-stage merger fraction, and they select more mergers in bulge-dominated galaxies. When applied to CANDELS morphology catalogs, the RFs estimate a merger rate increasing to at least z = 3, albeit two times higher than expected by theory. This suggests possible mismatches in the feedback-determined morphologies, but affirms the basic understanding of galaxy merger evolution. The RFs achieve completeness of roughly 70% at 0.5 < z < 3, and purity increasing from 10% at z = 0.5 to 60% at z = 3. At earlier times, the training sets are insufficient, motivating larger simulations and smaller time sampling. By blending large surveys and large simulations, such machine learning techniques offer a promising opportunity to teach us the strengths and weaknesses of inferences about galaxy evolution.

We present the full public release of all data from the 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 ~750 TB, including 1200 full volume snapshots and ~80,000 high time-resolution subbox snapshots. This will increase to ~1.1 PB with the future release of TNG50. 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.

Galaxies evolve under the influence of gas flows between their interstellar medium and their surrounding gaseous halos known as the circumgalactic medium (CGM). The CGM is a major reservoir of galactic baryons and metals, and plays a key role in the long cycles of accretion, feedback, and recycling of gas that drive star formation. In order to fully understand the physical processes at work within galaxies, it is therefore essential to have a firm understanding of the composition, structure, kinematics, thermodynamics, and evolution of the CGM. In this white paper we outline connections between the CGM and galactic star formation histories, internal kinematics, chemical evolution, quenching, satellite evolution, dark matter halo occupation, and the reionization of the larger-scale intergalactic medium in light of the advances that will be made on these topics in the 2020s. We argue that, in the next decade, fundamental progress on all of these major issues depends critically on improved empirical characterization and theoretical understanding of the CGM. In particular, we discuss how future advances in spatially-resolved CGM observations at high spectral resolution, broader characterization of the CGM across galaxy mass and redshift, and expected breakthroughs in cosmological hydrodynamic simulations will help resolve these major problems in galaxy evolution.

We present a comprehensive suite of high-resolution (parsec-scale), idealised (non-cosmological) galaxy merger simulations (24 runs, stellar mass ratio ~2.5:1) to investigate the connection between interaction-induced star formation and the evolution of the interstellar medium (ISM) in various temperature-density regimes. We use the GIZMO code and the second version of the 'Feedback in Realistic Environments' model (FIRE-2), which captures the multi-phase structure of the ISM. Our simulations are designed to represent galaxy mergers in the local Universe. In this work, we focus on the 'galaxy-pair period' between first and second pericentric passage. We split the ISM into four regimes: hot, warm, cool and cold-dense, motivated by the hot, ionised, atomic and molecular gas phases observed in real galaxies. We find that, on average, interactions enhance the star formation rate of the pair (~30%, merger-suite sample average) and elevate their cold-dense gas content (~18%). This is accompanied by a decrease in warm gas (~11%), a negligible change in cool gas (~4% increase), and a substantial increase in hot gas (~400%). The amount of cold-dense gas with densities above 1000 cm^3 (the cold ultra-dense regime) is elevated significantly (~240%), but only accounts for 0.15% (on average) of the cold-dense gas budget.

We use the IllustrisTNG (TNG) cosmological simulations to provide theoretical expectations for the dark matter mass fractions (DMFs) and circular velocity profiles of galaxies. TNG predicts flat circular velocity curves for $z = 0$ Milky Way (MW)-like galaxies beyond a few kpc from the galaxy centre, in better agreement with observational constraints than its predecessor, Illustris. TNG also predicts an enhancement of the dark matter mass within the 3D stellar half-mass radius ($r_\mathrm{half}$; $M_\mathrm{200c} = 10^{10}-10^{13}\mathrm{M}_{\odot}$, $z \le2$) compared to its dark matter only and Illustris counterparts. This enhancement leads TNG present-day galaxies to be dominated by dark matter within their inner regions, with $f_\mathrm{DM}(<r_\mathrm{half})\gtrsim0.5$ at all masses and with a minimum for MW-mass galaxies. The 1$\sigma$ scatter is $\lesssim$ 10~per~cent at all apertures, which is smaller than that inferred by some observational datasets, e.g. 40 per cent from the SLUGGS survey. TNG agrees with the majority of the observationally inferred values for elliptical galaxies once a consistent IMF is adopted (Chabrier) and the DMFs are measured within the same apertures. The DMFs measured within $r_\mathrm{half}$ increase towards lower redshifts: this evolution is dominated by the increase in galaxy size with time. At $z\sim2$, the DMF in disc-like TNG galaxies decreases with increasing galaxy mass, with $f_\mathrm{DM}(<r_\mathrm{half}) \sim 0.10-0.65$ for $10^{10} \lesssim M_{\rm stars}/\mathrm{M}_{\odot} \lesssim 10^{12}$, and are two times higher than if TNG galaxies resided in Navarro-Frenk-White dark matter haloes unaffected by baryonic physics. It remains to be properly assessed whether recent observational estimates of the DMFs at $z\sim2$ rule out the contraction of the dark matter haloes predicted by the TNG model.

We study the connection between active galactic nuclei (AGN) and their host galaxies through cosmic time in the large-scale cosmological IllustrisTNG simulations. We first compare BH properties, i.e. the hard X-ray BH luminosity function, AGN galaxy occupation fraction, and distribution of Eddington ratios, to available observational constraints. The simulations produce a population of BHs in good agreement with observations, but we note an excess of faint AGN in hard X-ray (L_x ~ 10^{43-44} erg/s), and a lower number of bright AGN (L_x>10^{44} erg/s), a conclusion that varies quantitatively but not qualitatively with BH luminosity estimation method. The lower Eddington ratios of the 10^{9} Msun BHs compared to observations suggest that AGN feedback may be too efficient in this regime. We study galaxy star formation activity and structural properties, and design sample-dependent criteria to identify different galaxy types (star-forming/quiescent, extended/compact) that we apply both to the simulations and observations from the candels fields. We analyze how the simulated and observed galaxies populate the specific star formation rate - stellar mass surface density diagram. A large fraction of the z=0 M_{star}>10^{11} Msun quiescent galaxies first experienced a compaction phase (i.e. reduction of galaxy size) while still forming stars, and then a quenching event. We measure the dependence of AGN fraction on galaxies' locations in this diagram. After correcting the simulations with a redshift and AGN luminosity-dependent model for AGN obscuration, we find good qualitative and quantitative agreement with observations. The AGN fraction is the highest among compact star-forming galaxies (16-20% at z~1.5-2), and the lowest among compact quiescent galaxies (6-10% at z~1.5-2).

We have generated synthetic images of $\sim$27,000 galaxies from the IllustrisTNG and the original Illustris hydrodynamic cosmological simulations, designed to match Pan-STARRS observations of $\log_{10}(M_{\ast}/{\rm M}_{\odot}) \approx 9.8$-$11.3$ galaxies at $z \approx 0.05$. Most of our synthetic images were created with the SKIRT radiative transfer code, including the effects of dust attenuation and scattering, and performing the radiative transfer directly on the Voronoi mesh used by the simulations themselves. We have analysed both our synthetic and real Pan-STARRS images with the newly developed $\tt{statmorph}$ code, which calculates non-parametric morphological diagnostics -- including the Gini-$M_{20}$ and concentration-asymmetry-smoothness (CAS) statistics -- and performs two-dimensional S\'ersic fits. Overall, we find that the optical morphologies of IllustrisTNG galaxies are in good agreement with observations, and represent a substantial improvement compared to the original Illustris simulation. In particular, the locus of the Gini-$M_{20}$ diagram is consistent with that inferred from observations, while the median trends with stellar mass of all the morphological, size and shape parameters considered in this work lie within the $\sim$1$\sigma$ scatter of the observational trends. However, the IllustrisTNG model has some difficulty with more stringent tests, such as producing a strong morphology-colour relation. This results in a somewhat higher fraction of red discs and blue spheroids compared to observations. Similarly, the morphology-size relation is problematic: while observations show that discs tend to be larger than spheroids at a fixed stellar mass, such a trend is not present in IllustrisTNG.

The thermodynamic structure of hot gas in galaxy clusters is sensitive to astrophysical processes and typically difficult to model with galaxy formation simulations. We explore the fraction of cool-core (CC) clusters in a large sample of $370$ clusters from IllustrisTNG, examining six common CC definitions. IllustrisTNG produces continuous CC criteria distributions, the extremes of which are classified as CC and non-cool-core (NCC), and the criteria are increasingly correlated for more massive clusters. At $z=0$, the CC fractions for $2$ criteria are in reasonable agreement with the observed fractions but the other $4$ CC fractions are lower than observed. This result is partly driven by systematic differences between the simulated and observed gas fraction profiles. The simulated CC fractions with redshift show tentative agreement with the observed fractions, but linear fits demonstrate that the simulated evolution is steeper than observed. The conversion of CCs to NCCs appears to begin later and act more rapidly in the simulations. Examining the fraction of CCs and NCCs defined as relaxed we find no evidence that CCs are more relaxed, suggesting that mergers are not solely responsible for disrupting CCs. A comparison of the median thermodynamic profiles defined by different CC criteria shows that the extent to which they evolve in the cluster core is dependent on the CC criteria. We conclude that the thermodynamic structure of galaxy clusters in IllustrisTNG shares many similarities with observations, but achieving better agreement most likely requires modifications of the underlying galaxy formation model.

The Feedback In Realistic Environments (FIRE) project explores feedback in cosmological galaxy formation simulations. Previous FIRE simulations used an identical source code (FIRE-1) for consistency. Motivated by the development of more accurate numerics - including hydrodynamic solvers, gravitational softening, and supernova coupling algorithms - and exploration of new physics (e.g. magnetic fields), we introduce FIRE-2, an updated numerical implementation of FIRE physics for the GIZMO code. We run a suite of simulations and compare against FIRE-1: overall, FIRE-2 improvements do not qualitatively change galaxy-scale properties. We pursue an extensive study of numerics versus physics. Details of the star-formation algorithm, cooling physics, and chemistry have weak effects, provided that we include metal-line cooling and star formation occurs at higher-than-mean densities. We present new resolution criteria for high-resolution galaxy simulations. Most galaxy-scale properties are robust to numerics we test, provided: (1) Toomre masses are resolved; (2) feedback coupling ensures conservation, and (3) individual supernovae are time-resolved. Stellar masses and profiles are most robust to resolution, followed by metal abundances and morphologies, followed by properties of winds and circum-galactic media (CGM). Central (~kpc) mass concentrations in massive (L*) galaxies are sensitive to numerics (via trapping/recycling of winds in hot halos). Multiple feedback mechanisms play key roles: supernovae regulate stellar masses/winds; stellar mass-loss fuels late star formation; radiative feedback suppresses accretion onto dwarfs and instantaneous star formation in disks. We provide all initial conditions and numerical algorithms used.

Large-scale cosmological simulations of galaxy formation currently do not resolve the densities at which molecular hydrogen forms, implying that the atomic-to-molecular transition must be modeled either on the fly or in postprocessing. We present an improved postprocessing framework to estimate the abundance of atomic and molecular hydrogen and apply it to the IllustrisTNG simulations. We compare five different models for the atomic-to-molecular transition, including empirical, simulation-based, and theoretical prescriptions. Most of these models rely on the surface density of neutral hydrogen and the ultraviolet (UV) flux in the Lyman-Werner band as input parameters. Computing these quantities on the kiloparsec scales resolved by the simulations emerges as the main challenge. We show that the commonly used Jeans length approximation to the column density of a system can be biased and exhibits large cell-to-cell scatter. Instead, we propose to compute all surface quantities in face-on projections and perform the modeling in two dimensions. In general, the two methods agree on average, but their predictions diverge for individual galaxies and for models based on the observed midplane pressure of galaxies. We model the UV radiation from young stars by assuming a constant escape fraction and optically thin propagation throughout the galaxy. With these improvements, we find that the five models for the atomic-to-molecular transition roughly agree on average but that the details of the modeling matter for individual galaxies and the spatial distribution of molecular hydrogen. We emphasize that the estimated molecular fractions are approximate due to the significant systematic uncertainties.

We introduce the IllustrisTNG project, a new suite of cosmological magnetohydrodynamical simulations performed with the moving-mesh code AREPO employing an updated Illustris galaxy formation model. Here we focus on the general properties of magnetic fields and the diffuse radio emission in galaxy clusters. Magnetic fields are prevalent in galaxies, and their build-up is closely linked to structure formation. We find that structure formation amplifies the initial seed fields ($10^{-14}$ comoving Gauss) to the values observed in low-redshift galaxies ($1-10\,\mu{\rm G}$). The magnetic field topology is closely connected to galaxy morphology such that irregular fields are hosted by early-type galaxies, while large-scale, ordered fields are present in disc galaxies. Using two simple models for the energy distribution of relativistic electrons we predict the diffuse radio emission of $280$ clusters with a baryonic mass resolution of $1.1\times 10^{7}\,{\rm M_{\odot}}$, and generate mock observations for VLA, LOFAR, ASKAP and SKA. Our simulated clusters show extended radio emission, whose detectability correlates with their virial mass. We reproduce the observed scaling relations between total radio power and X-ray emission, $M_{500}$, and the Sunyaev-Zel'dovich $Y_{\rm 500}$ parameter. The radio emission surface brightness profiles of our most massive clusters are in reasonable agreement with VLA measurements of Coma and Perseus. Finally, we discuss the fraction of detected extended radio haloes as a function of virial mass and source count functions for different instruments. Overall our results agree encouragingly well with observations, but a refined analysis requires a more sophisticated treatment of relativistic particles in large-scale galaxy formation simulations.

The galaxy Malin 1 contains the largest stellar disk known but the formation mechanism of this structure has been elusive. In this paper, we report a Malin 1 analogue in the 100 Mpc IllustrisTNG simulation and describe its formation history. At redshift zero, this massive galaxy, having a maximum circular velocity $V_{\rm max}$ of 430 ${\rm km\ s^{-1}}$, contains a 100 kpc gas/stellar disk with morphology similar to Malin 1. The simulated galaxy reproduces well many observed features of Malin 1's vast disk, including its stellar ages, metallicities, and gas rotation curve. We trace the extended disk back in time and find that a large fraction of the cold gas at redshift zero originated from the cooling of hot halo gas, triggered by the merger of a pair of intruding galaxies. Our finding provides a novel way to form large galaxy disks as extreme as Malin 1 within the current galaxy formation framework.

Understanding how baryonic processes shape the intracluster medium (ICM) is of critical importance to the next generation of galaxy cluster surveys. However, most models of structure formation neglect potentially important physical processes, like anisotropic thermal conduction (ATC). In this letter, we explore the impact of ATC on the prevalence of cool-cores (CCs) using 12 pairs of magnetohydrodynamical galaxy cluster simulations, simulated using the IllustrisTNG model with and without ATC. Although the impact of ATC varies from cluster to cluster and with CC criterion, its inclusion produces a systematic shift to larger CC fractions at z = 0 for all CC criteria considered. Additionally, the inclusion of ATC yields a flatter CC fraction redshift evolution, easing the tension with the observed evolution. With ATC included, the energy required for the central black hole to achieve self-regulation is reduced and the gas fraction in the cluster core increases, resulting in larger CC fractions. ATC makes the ICM unstable to perturbations and the increased efficiency of AGN feedback suggests that its inclusion results in a greater level of mixing in the ICM. Therefore, ATC is potentially an important physical process in reproducing the thermal structure of the ICM.

Interstellar dust is an important component of the galactic ecosystem, playing a key role in multiple galaxy formation processes. We present a novel numerical framework for the dynamics and size evolution of dust grains implemented in the moving-mesh hydrodynamics code AREPO suited for cosmological galaxy formation simulations. We employ a particle-based method for dust subject to dynamical forces including drag and gravity. The drag force is implemented using a second-order semi-implicit integrator and validated using several dust-hydrodynamical test problems. Each dust particle has a grain size distribution, describing the local abundance of grains of different sizes. The grain size distribution is discretised with a second-order piecewise linear method and evolves in time according to various dust physical processes, including accretion, sputtering, shattering, and coagulation. We present a novel scheme for stochastically forming dust during stellar evolution and new methods for sub-cycling of dust physics time-steps. Using this model, we simulate an isolated disc galaxy to study the impact of dust physical processes that shape the interstellar grain size distribution. We demonstrate, for example, how dust shattering shifts the grain size distribution to smaller sizes resulting in a significant rise of radiation extinction from optical to near-ultraviolet wavelengths. Our framework for simulating dust and gas mixtures can readily be extended to account for other dynamical processes relevant in galaxy formation, like magnetohydrodynamics, radiation pressure, and thermo-chemical processes.

We explore the abundance, spatial distribution, and physical properties of the OVI, OVII, and OVIII ions of oxygen in circumgalactic and intergalactic media (the CGM, IGM, and WHIM). We use the TNG100 and TNG300 large volume cosmological magneto-hydrodynamical simulations. Modeling the ionization states of simulated oxygen, we find good agreement with observations of the low-redshift OVI column density distribution function (CDDF), and present its evolution for all three ions from z=0 to z=4. Producing mock quasar absorption line spectral surveys, we show that the IllustrisTNG simulations are fully consistent with constraints on the OVI content of the CGM from COS-Halos and other low redshift observations, producing columns as high as observed. We measure the total amount of mass and average column densities of each ion using hundreds of thousands of simulated galaxies spanning 10^11 < Mhalo/Msun < 10^15 corresponding to 10^9 < M*/Msun < 10^12 in stellar mass. Stacked radial profiles of OVI are computed in 3D number density and 2D projected column, decomposing into the 1-halo and 2-halo terms. Relating halo OVI to properties of the central galaxy, we find a correlation between the (g-r) color of a galaxy and the total amount of OVI in its CGM. In comparison to the COS-Halos finding, this leads to a dichotomy of columns around star-forming versus passive galaxies at fixed stellar (or halo) mass. We demonstrate that this correlation is a direct result of blackhole feedback associated with quenching and represents a causal consequence of galactic-scale baryonic feedback impacting the physical state of the circumgalactic medium.

The distribution of elements in galaxies provides a wealth of information about their production sites and their subsequent mixing into the interstellar medium. Here we investigate the distribution of elements within stars in the IllustrisTNG simulations. In particular, we analyze the abundance ratios of magnesium and europium in Milky Way-like galaxies from the TNG100 simulation (stellar masses ${\log} (M_\star / {\rm M}_\odot) \sim 9.7 - 11.2$). As abundances of magnesium and europium for individual stars in the Milky Way are observed across a variety of spatial locations and metallicities, comparison with the stellar abundances in our more than $850$ Milky Way-like galaxies provides stringent constraints on our chemical evolutionary methods. To this end we use the magnesium to iron ratio as a proxy for the effects of our SNII and SNIa metal return prescription, and a means to compare our simulated abundances to a wide variety of galactic observations. The europium to iron ratio tracks the rare ejecta from neutron star -- neutron star mergers, the assumed primary site of europium production in our models, which in turn is a sensitive probe of the effects of metal diffusion within the gas in our simulations. We find that europium abundances in Milky Way-like galaxies show no correlation with assembly history, present day galactic properties, and average galactic stellar population age. In general, we reproduce the europium to iron spread at low metallicities observed in the Milky Way, with the level of enhancement being sensitive to gas properties during redshifts $z \approx 2-4$. We show that while the overall normalization of [Eu/Fe] is susceptible to resolution and post-processing assumptions, the relatively large spread of [Eu/Fe] at low [Fe/H] when compared to that at high [Fe/H] is very robust.

We use the IllustrisTNG simulations to investigate the evolution of the mass-metallicity relation (MZR) for star-forming cluster galaxies as a function of the formation history of their cluster host. The simulations predict an enhancement in the gas-phase metallicities of star-forming cluster galaxies (10^9< M_star<10^10 M_sun) at z<1.0 in comparisons to field galaxies. This is qualitatively consistent with observations. We find that the metallicity enhancement of cluster galaxies appears prior to their infall into the central cluster potential, indicating for the first time a systematic "chemical pre-processing" signature for {\it infalling} cluster galaxies. Namely, galaxies which will fall into a cluster by z=0 show a ~0.05 dex enhancement in the MZR compared to field galaxies at z<0.5. Based on the inflow rate of gas into cluster galaxies and its metallicity, we identify that the accretion of pre-enriched gas is the key driver of the chemical evolution of such galaxies, particularly in the stellar mass range (10^9< M_star<10^10 M_sun). We see signatures of an environmental dependence of the ambient/inflowing gas metallicity which extends well outside the nominal virial radius of clusters. Our results motivate future observations looking for pre-enrichment signatures in dense environments.

The IllustrisTNG project is a new suite of cosmological magneto-hydrodynamical simulations of galaxy formation performed with the Arepo code and updated models for feedback physics. Here we introduce the first two simulations of the series, TNG100 and TNG300, and quantify the stellar mass content of about 4000 massive galaxy groups and clusters ($10^{13} \leq M_{\rm 200c}/M_{\rm sun} \leq 10^{15}$) at recent times ($z \leq 1$). The richest clusters have half of their total stellar mass bound to satellite galaxies, with the other half being associated with the central galaxy and the diffuse intra-cluster light. The exact ICL fraction depends sensitively on the definition of a central galaxy's mass and varies in our most massive clusters between 20 to 40% of the total stellar mass. Haloes of $5\times 10^{14}M_{\rm sun}$ and above have more diffuse stellar mass outside 100 kpc than within 100 kpc, with power-law slopes of the radial mass density distribution as shallow as the dark matter's ( $-3.5 < \alpha_{\rm 3D} < -3$). Total halo mass is a very good predictor of stellar mass, and vice versa: at $z=0$, the 3D stellar mass measured within 30 kpc scales as $\propto (M_{\rm 500c})^{0.49}$ with a $\sim 0.12$ dex scatter. This is possibly too steep in comparison to the available observational constraints, even though the abundance of TNG less massive galaxies ($< 10^{11}M_{\rm sun}$ in stars) is in good agreement with the measured galaxy stellar mass functions at recent epochs. The 3D sizes of massive galaxies fall too on a tight ($\sim$0.16 dex scatter) power-law relation with halo mass, with $r^{\rm stars}_{\rm 0.5} \propto (M_{\rm 500c})^{0.53}$. Even more fundamentally, halo mass alone is a good predictor for the whole stellar mass profiles beyond the inner few kpc, and we show how on average these can be precisely recovered given a single mass measurement of the galaxy or its halo.

We introduce the first two simulations of the IllustrisTNG project, a next generation of cosmological magnetohydrodynamical simulations, focusing on the optical colors of galaxies. We explore TNG100, a rerun of the original Illustris box, and TNG300, which includes 2x2500^3 resolution elements in a volume twenty times larger. Here we present first results on the galaxy color bimodality at low redshift. Accounting for the attenuation of stellar light by dust, we compare the simulated (g-r) colors of 10^9 < M*/Msun < 10^12.5 galaxies to the observed distribution from the Sloan Digital Sky Survey (SDSS). We find a striking improvement with respect to the original Illustris simulation, as well as excellent quantitative agreement in comparison to the observations, with a sharp transition in median color from blue to red at a characteristic M* ~ 10^10.5 Msun. Investigating the build-up of the color-mass plane and the formation of the red sequence, we demonstrate that the primary driver of galaxy color transition in the TNG model is supermassive blackhole feedback in its low-accretion state. Across the entire population we measure a median color transition timescale dt_green of ~1.6 Gyr, a value which drops for increasingly massive galaxies. We find signatures of the physical process of quenching: at fixed stellar mass, redder galaxies have lower SFRs, gas fractions, and gas metallicities; their stellar populations are also older and their large-scale interstellar magnetic fields weaker than in bluer galaxies. Finally, we measure the amount of stellar mass growth on the red sequence. Galaxies with M* > 10^11 Msun which redden at z<1 accumulate on average ~25% of their final z=0 mass post-reddening; at the same time, ~18% of such massive galaxies acquire half or more of their final stellar mass while on the red sequence.

Modern cosmological simulations model the universe with increasing sophistication and at higher spatial and temporal resolutions. These enhancements permit detailed comparisons between the simulation outputs and real observational data. Recent projects such as Illustris are capable of producing simulated images that are comparable to those obtained from local surveys. This paper tests how well Illustris achieves this goal across a diverse population of galaxies using visual morphologies derived from Galaxy Zoo citizen scientists. Morphological classifications provided by volunteers for simulated galaxies are compared with similar data for a compatible sample of images drawn from the SDSS Legacy Survey. This paper investigates how simple morphological characterization by human volunteers asked to distinguish smooth from featured systems differs between simulated and real galaxy images. Differences are identified, which are likely due to the limited resolution of the simulation, but which could be revealing real differences in the dynamical evolution of populations of galaxies in the real and model universes. Specifically, for stellar masses $M_{\star}\lesssim10^{11}M_{\odot}$, a larger proportion of Illustris galaxies that exhibit disk-like morphology or visible substructure, relative to their SDSS counterparts. Toward higher masses, simulated and observed galaxies converge and exhibit similar morphology distributions. The stellar mass threshold indicated by this divergent behavior confirms recent works using parametric measures of morphology from Illustris simulated images. When $M_{\star}\gtrsim10^{11}M_{\odot}$, the Illustris dataset contains fewer galaxies that classifiers regard as unambiguously featured. These results suggest that comparison between the detailed properties of observed and simulated galaxies, even when limited to reasonably massive systems, may be misleading.

Galaxies are surrounded by sizeable gas reservoirs which host a significant amount of metals: the circum-galactic medium (CGM). The CGM acts as a mediator between the galaxy and the extra-galactic medium. However, our understanding of how galaxy mergers, a major evolutionary transformation, impact the CGM remains deficient. We present a theoretical study of the effect of galaxy mergers on the CGM. We use hydrodynamical cosmological zoom-in simulations of a major merger selected from the Illustris project such that the z=0 descendant has a halo mass and stellar mass comparable to the Milky Way. To study the CGM we then re-simulated this system at a 40 times better mass resolution, and included detailed post-processing ionization modelling. Our work demonstrates the effect the merger has on the characteristic size of the CGM, its metallicity, and the predicted covering fraction of various commonly observed gas-phase species, such as H I, C IV, and O VI. We show that merger-induced outflows can increase the CGM metallicity by 0.2-0.3 dex within 0.5 Gyr post-merger. These effects last up to 6 Gyr post-merger. While the merger increases the total metal covering fractions by factors of 2-3, the covering fractions of commonly observed UV ions decrease due to the hard ionizing radiation from the active galactic nucleus, which we model explicitly. Our study of the single simulated major merger presented in this work demonstrates the significant impact that a galaxy interaction can have on the size, metallicity, and observed column densities of the CGM.

Hydrodynamical simulations of galaxy formation have now reached sufficient volume to make precision predictions for clustering on cosmologically relevant scales. Here we use our new IllustrisTNG simulations to study the non-linear correlation functions and power spectra of baryons, dark matter, galaxies and haloes over an exceptionally large range of scales. We find that baryonic effects increase the clustering of dark matter on small scales and damp the total matter power spectrum on scales up to k ~ 10 h/Mpc by 20%. The non-linear two-point correlation function of the stellar mass is close to a power-law over a wide range of scales and approximately invariant in time from very high redshift to the present. The two-point correlation function of the simulated galaxies agrees well with SDSS at its mean redshift z ~ 0.1, both as a function of stellar mass and when split according to galaxy colour, apart from a mild excess in the clustering of red galaxies in the stellar mass range 10^9-10^10 Msun/h^2. Given this agreement, the TNG simulations can make valuable theoretical predictions for the clustering bias of different galaxy samples. We find that the clustering length of the galaxy auto-correlation function depends strongly on stellar mass and redshift. Its power-law slope gamma is nearly invariant with stellar mass, but declines from gamma ~ 1.8 at redshift z=0 to gamma ~ 1.6 at redshift z ~ 1, beyond which the slope steepens again. We detect significant scale-dependencies in the bias of different observational tracers of large-scale structure, extending well into the range of the baryonic acoustic oscillations and causing nominal (yet fortunately correctable) shifts of the acoustic peaks of around ~5%.

We present final statistics from a survey for intervening MgII absorption towards 100 quasars with emission redshifts between $z=3.55$ and $z=7.08$. Using infrared spectra from Magellan/FIRE, we detect 279 cosmological MgII absorbers, and confirm that the incidence rate of $W_r>0.3 \AA$ MgII absorption per comoving path length does not evolve measurably between $z=0.25$ and $z=7$. This is consistent with our detection of seven new MgII systems at $z>6$, a redshift range that was not covered in prior searches. Restricting to relatively strong MgII systems ($W_r>1$\AA), there is significant evidence for redshift evolution. These systems roughly double in number density between $z=0$ and $z=2$-$3$, but decline by an order of magnitude from this peak by $z\sim 6$. This evolution mirrors that of the global star formation rate density, which could reflect a connection between star formation feedback and strong MgII absorbers. We compared our results to the Illustris cosmological simulation at $z=2$-$4$ by assigning absorption to catalogued dark-matter halos and by direct extraction of spectra from the simulation volume. To reproduce our results using the halo catalogs, we require circumgalactic (CGM) MgII envelopes within halos of progressively smaller mass at earlier times. This occurs naturally if we define the lower integration cutoff using SFR rather than mass. MgII profiles calculated directly from the Illustris volume yield far too few strong absorbers. We argue that this arises from unresolved phase space structure of CGM gas, particularly from turbulent velocities on sub-mesh scales. The presence of CGM MgII at $z>6$-- just $\sim 250$ Myr after the reionization redshift implied by Planck--suggests that enrichment of intra-halo gas may have begun before the presumed host galaxies' stellar populations were mature and dynamically relaxed. [abridged]

The fundamental metallicity relation (FMR) is a postulated correlation between galaxy stellar mass, star formation rate (SFR), and gas-phase metallicity. At its core, this relation posits that offsets from the mass-metallicity relation (MZR) at a fixed stellar mass are correlated with galactic SFR. In this Letter, we quantify the timescale with which galactic SFRs and metallicities evolve using hydrodynamical simulations. We find that Illustris and IllustrisTNG predict that galaxy offsets from the star formation main sequence and MZR evolve over similar timescales, are often anti-correlated in their evolution, evolve with the halo dynamical time, and produce a pronounced FMR. In fact, for a FMR to exist, the metallicity and SFR must evolve in an anti-correlated sense which requires that they evolve with similar time variability. In contrast to Illustris and IllustrisTNG, we speculate that the SFR and metallicity evolution tracks may become decoupled in galaxy formation models dominated by globally-bursty SFR histories, which could weaken the FMR residual correlation strength. This opens the possibility of discriminating between bursty and non-bursty feedback models based on the strength and persistence of the FMR -- especially at high redshift.

We analyze scaling relations and evolution histories of galaxy sizes in TNG100, part of the IllustrisTNG simulation suite. Observational qualitative trends of size with stellar mass, star-formation rate and redshift are reproduced, and a quantitative comparison of projected r-band sizes at 0~<z<~2 shows agreement to much better than 0.25dex. We follow populations of z=0 galaxies with a range of masses backwards in time along their main progenitor branches, distinguishing between main-sequence and quenched galaxies. Our main findings are as follows. (i) At M_{*,z=0}>~10^{9.5}Msun, the evolution of the median main progenitor differs, with quenched galaxies hardly growing in median size before quenching, whereas main-sequence galaxies grow their median size continuously, thus opening a gap from the progenitors of quenched galaxies. This is partly because the main-sequence high-redshift progenitors of quenched z=0 galaxies are drawn from the lower end of the size distribution of the overall population of main-sequence high-redshift galaxies. (ii) Quenched galaxies with M_{*,z=0}>~10^{9.5}Msun experience a steep size growth on the size-mass plane after their quenching time, but with the exception of galaxies with M_{*,z=0}>~10^{11}Msun, the size growth after quenching is small in absolute terms, such that most of the size (and mass) growth of quenched galaxies (and its variation among them) occurs while they are still on the main-sequence. After they become quenched, the size growth rate of quenched galaxies as a function of time, as opposed to versus mass, is similar to that of main-sequence galaxies. Hence, the size gap is retained down to z=0.

The coevolution of galaxies and their metal content serves as an important test for galaxy feedback models. We analyze the distribution and evolution of metals within the IllustrisTNG simulation suite with a focus on the gas-phase mass-metallicity relation (MZR). We find that the IllustrisTNG model broadly reproduces the slope and normalization evolution of the MZR across the redshift range $0<z<2$ and mass range $10^9 < M_*/\mathrm{M}_\odot < 10^{10.5}$. We make predictions for the high redshift ($2<z<10$) metal content of galaxies which is described by a gradual decline in the normalization of the metallicity with an average high redshift ($z>2$) evolution fit by $\mathrm{d\;log(Z)}/\mathrm{dz} \approx - 0.064$. Our simulations indicate that the metal retention efficiency of the interstellar medium (ISM) is low: a majority of gas-phase metals ($\sim$ 85 per cent at $z=0$) live outside of the ISM, either in an extended gas disk, the circumgalactic medium, or outside the halo. Nevertheless, the redshift evolution in the simulated MZR normalization is driven by the higher gas fractions of high redshift galaxies, not by changes to the metal retention efficiency. The scatter in the simulated MZR contains a clear correlation with the gas-mass or star formation rate of the system, in agreement with the observed fundamental metallicity relation. The scatter in the MZR is driven by a competition between periods of enrichment- and accretion-dominated metallicity evolution. We expect that while the normalization of the MZR declines with redshift, the strength of the correlation between metallicity and gas-mass at fixed stellar mass is not a strong function of redshift. Our results indicate that the "regulator" style models are best suited for simultaneously explaining the shape, redshift evolution, and existence of correlated scatter with gas fraction about the MZR.

The distribution of metals in the intra-cluster medium encodes important information about the enrichment history and formation of galaxy clusters. Here we explore the metal content of clusters in IllustrisTNG - a new suite of galaxy formation simulations building on the Illustris project. Our cluster sample contains 20 objects in TNG100 - a ~(100 Mpc)^3 volume simulation with 2x1820^3 resolution elements, and 370 objects in TNG300 - a ~(300 Mpc)^3 volume simulation with 2x2500^3 resolution elements. The z=0 metallicity profiles agree with observations, and the enrichment history is consistent with observational data going beyond z~1, showing nearly no metallicity evolution. The abundance profiles vary only minimally within the cluster samples, especially in the outskirts with a relative scatter of ~15%. The average metallicity profile flattens towards the center, where we find a logarithmic slope of -0.1 compared to -0.5 in the outskirts. Cool core clusters have more centrally peaked metallicity profiles (~0.8 solar) compared to non-cool core systems (~0.5 solar), similar to observational trends. Si/Fe and O/Fe radial profiles follow positive gradients. The outer abundance profiles do not evolve below z~2, whereas the inner profiles flatten towards z=0. More than ~80% of the metals in the intra-cluster medium have been accreted from the proto-cluster environment, which has been enriched to ~0.1 solar already at z~2. We conclude that the intra-cluster metal distribution is uniform among our cluster sample, nearly time-invariant in the outskirts for more than 10 Gyr, and forms through a universal enrichment history.

We introduce an updated physical model to simulate the formation and evolution of galaxies in cosmological, large-scale gravity+magnetohydrodynamical simulations with the moving mesh code AREPO. The overall framework builds upon the successes of the Illustris galaxy formation model, and includes prescriptions for star formation, stellar evolution, chemical enrichment, primordial and metal-line cooling of the gas, stellar feedback with galactic outflows, and black hole formation, growth and multi-mode feedback. In this paper we give a comprehensive description of the physical and numerical advances which form the core of the IllustrisTNG (The Next Generation) framework. We focus on the revised implementation of the galactic winds, of which we modify the directionality, velocity, thermal content, and energy scalings, and explore its effects on the galaxy population. As described in earlier works, the model also includes a new black hole driven kinetic feedback at low accretion rates, magnetohydrodynamics, and improvements to the numerical scheme. Using a suite of (25 Mpc $h^{-1}$)$^3$ cosmological boxes we assess the outcome of the new model at our fiducial resolution. The presence of a self-consistently amplified magnetic field is shown to have an important impact on the stellar content of $10^{12} M_{\rm sun}$ haloes and above. Finally, we demonstrate that the new galactic winds promise to solve key problems identified in Illustris in matching observational constraints and affecting the stellar content and sizes of the low mass end of the galaxy population.

We study the population of supermassive black holes (SMBHs) and their effects on massive central galaxies in the IllustrisTNG cosmological hydrodynamical simulations of galaxy formation. The employed model for SMBH growth and feedback assumes a two-mode scenario in which the feedback from active galactic nuclei occurs through a kinetic, comparatively efficient mode at low accretion rates relative to the Eddington limit, and in the form of a thermal, less efficient mode at high accretion rates. We show that the quenching of massive central galaxies happens coincidently with kinetic-mode feedback, consistent with the notion that active supermassive black cause the low specific star formation rates observed in massive galaxies. However, major galaxy mergers are not responsible for initiating most of the quenching events in our model. Up to black hole masses of about $10^{8.5}\,{\rm M}_\odot$, the dominant growth channel for SMBHs is in the thermal mode. Higher mass black holes stay mainly in the kinetic mode and gas accretion is self-regulated via their feedback, which causes their Eddington ratios to drop, with SMBH mergers becoming the main channel for residual mass growth. As a consequence, the quasar luminosity function is dominated by rapidly accreting, moderately massive black holes in the thermal mode. We show that the associated growth history of SMBHs produces a low-redshift quasar luminosity function and a redshift zero black hole mass-stellar bulge mass relation in good agreement with observations, whereas the simulation tends to over-predict the high-redshift quasar luminosity function.

We introduce massive black holes (BHs) in the Feedback In Realistic Environments project and perform high-resolution cosmological hydrodynamic simulations of quasar-mass halos ($M_{\rm halo}(z=2) \approx 10^{12.5}\,\rm{M}_{\odot}$) down to $z=1$. These simulations model stellar feedback by supernovae, stellar winds, and radiation, and BH growth using a gravitational torque-based prescription tied to resolved properties of galactic nuclei. We do not include BH feedback. We show that early BH growth occurs through short ($\lesssim 1\,$Myr) accretion episodes that can reach or even exceed the Eddington rate. In this regime, BH growth is limited by bursty stellar feedback continuously evacuating gas from galactic nuclei, and BHs remain under-massive relative to the local $M_{\rm BH}$-$M_{\rm bulge}$ relation. BH growth is more efficient at later times, when the nuclear stellar potential retains a significant gas reservoir, star formation becomes less bursty, and galaxies settle into a more ordered state, with BHs rapidly converging onto the scaling relation when the host reaches $M_{\rm bulge} \sim 10^{10}\,\rm{M}_{\odot}$. Our results are not sensitive to the details of the accretion model so long as BH growth is tied to the gas content within $\sim 100\,$pc of the BH. Our simulations imply that bursty stellar feedback has strong implications for BH and AGN demographics, especially in the early Universe and for low-mass galaxies.

We present full volume cosmological simulations using the moving-mesh code AREPO to study the coevolution of dust and galaxies. We extend the dust model in AREPO to include thermal sputtering of grains and investigate the evolution of the dust mass function, the cosmic distribution of dust beyond the interstellar medium, and the dependence of dust-to-stellar mass ratio on galactic properties. The simulated dust mass function is well-described by a Schechter fit and lies closest to observations at $z = 0$. The radial scaling of projected dust surface density out to distances of $10 \, \text{Mpc}$ around galaxies with magnitudes $17 < i < 21$ is similar to that seen in Sloan Digital Sky Survey data, albeit with a lower normalisation. At $z = 0$, the predicted dust density of $\Omega_\text{dust} \approx 1.3 \times 10^{-6}$ lies in the range of $\Omega_\text{dust}$ values seen in low-redshift observations. We find that dust-to-stellar mass ratio anti-correlates with stellar mass for galaxies living along the star formation main sequence. Moreover, we estimate the $850 \, \mu\text{m}$ number density functions for simulated galaxies and analyse the relation between dust-to-stellar flux and mass ratios at $z = 0$. At high redshift, our model fails to produce enough dust-rich galaxies, and this tension is not alleviated by adopting a top-heavy initial mass function. We do not capture a decline in $\Omega_\text{dust}$ from $z = 2$ to $z = 0$, which suggests that dust production mechanisms more strongly dependent on star formation may help to produce the observed number of dusty galaxies near the peak of cosmic star formation.

Gladders et al. have recently suggested that the star formation histories (SFHs) of individual galaxies are characterized by a log-normal function in time, implying a slow decline rather than rapid quenching. We test their conjecture on theoretical SFHs from the cosmological simulation Illustris and on observationally inferred SFHs. While the log-normal form necessarily ignores short-lived features such as starbursts, it fits the overall shape of the majority of SFHs very well. In particular, 85% of the cumulative SFHs are fitted to within a maximum error of 5% of the total stellar mass formed, and 99% to within 10%. The log-normal performs systematically better than the commonly used delayed-tau model, and is superseded only by functions with more than three free parameters. Poor fits are mostly found in galaxies that were rapidly quenched after becoming satellites. We explore the log-normal parameter space of normalization, peak time, and full width at half maximum, and find that the simulated and observed samples occupy similar regions, though Illustris predicts wider, later-forming SFHs on average. The ensemble of log-normal fits correctly reproduces complex metrics such as the evolution of Illustris galaxies across the star formation main sequence, but overpredicts their quenching timescales. SFHs in Illustris are a diverse population not determined by any one physical property of galaxies, but follow a tight relation, where $\mathrm{width}\propto\mathrm{(peak\ time)}^{3/2}$. We show that such a relation can be explained qualitatively (though not quantitatively) by a close connection between the growth of dark matter halos and their galaxies.

We present an analysis of the flow of metals through the circumgalactic medium (CGM) in the Feedback in Realistic Environments (FIRE) simulations of galaxy formation, ranging from isolated dwarfs to $L*$ galaxies. We find that nearly all metals produced in high-redshift galaxies are carried out in winds that reach $0.25 R_{\rm vir}$. When measured at $0.25 R_{\rm vir}$ the metallicity of outflows is slightly higher than the interstellar medium (ISM) metallicity. Many metals thus reside in the CGM. Cooling and recycling from this reservoir determine the metal budget in the ISM. The outflowing metal flux decreases by a factor of $\sim2-5$ between $0.25 R_{\rm vir}$ and $R_{\rm vir}$. Furthermore, outflow metallicity is typically lower at $R_{\rm vir}$ owing to dilution of the remaining outflow by metal-poor material swept up from the CGM. The inflow metallicity at $R_{\rm vir}$ is generally low, but outflow and inflow metallicities are similar in the inner halo. At low redshift, massive galaxies no longer generate outflows that reach the CGM, causing a divergence in CGM and ISM metallicity. Dwarf galaxies continue to generate outflows, although they preferentially retain metal ejecta. In all but the least massive galaxy considered, a majority of the metals are within the halo at $z=0$. We measure the fraction of metals in CGM, ISM and stars, and quantify the thermal state of CGM metals in each halo. The total amount of metals in the low-redshift CGM of two simulated $L*$ galaxies is consistent with estimates from the COS halos survey, while for the other two it appears to be lower.

We compare mass-selected close pairs at z > 1 with the intrinsic galaxy merger rate in the Illustris Simulations. To do so, we construct three 140 arcmin^2 lightcone catalogs and measure pair fractions, finding that they change little or decrease with increasing redshift at z > 1. Consistent with current surveys, this trend requires a decrease in the merger-pair observability time, roughly as (1 + z)^-2, in order to measure the merger rates of the same galaxies. This implies that major mergers are more common at high redshift than implied by the simplest arguments assuming a constant observability time. Several effects contribute to this trend: (1) The fraction of massive, major (4:1) pairs which merge by today increases weakly from ~0.5 at z=1 to ~0.8 at z=3. (2) The median time elapsed between an observed pair and final remnant decreases by a factor of two from z~1 to z~3. (3) An increasing specific star formation rate (sSFR) decreases the time during which common stellar-mass based pair selection criteria could identify the mergers. The average orbit of the pairs at observation time varies only weakly, suggesting that the dynamical time is not varying enough to account by itself for the pair fraction trends. Merging pairs reside in dense regions, having overdensity ~10 to ~100 times greater than the average massive galaxy. We forward model the pairs to reconstruct the merger remnant production rate, showing that it is consistent with a rapid increase in galaxy merger rates at z > 1.

The inefficiency of star formation in massive elliptical galaxies is widely believed to be caused by the interactions of an active galactic nucleus (AGN) with the surrounding gas. Achieving a sufficiently rapid reddening of moderately massive galaxies without expelling too many baryons has however proven difficult for hydrodynamical simulations of galaxy formation, prompting us to explore a new model for the accretion and feedback effects of supermassive black holes. For high accretion rates relative to the Eddington limit, we assume that a fraction of the accreted rest mass energy heats the surrounding gas thermally, similar to the `quasar mode' in previous work. For low accretion rates, we invoke a new, pure kinetic feedback model which imparts momentum into the surrounding gas in a stochastic manner. These two modes of feedback are motivated both by theoretical conjectures for the existence of different types of accretion flows as well as recent observational evidence for the importance of kinetic AGN winds in quenching galaxies. We find that a large fraction of the injected kinetic energy in this mode thermalises via shocks in the surrounding gas, thereby providing a distributed heating channel. In cosmological simulations, the resulting model produces red, non star-forming massive elliptical galaxies, and achieves realistic gas fractions, black hole growth histories and thermodynamic profiles in large haloes.

We present 90 mas (37 pc) resolution ALMA imaging of Arp 220 in the CO (1-0) line and continuum at $\lambda = 2.6$ mm. The internal gas distribution and kinematics of both galactic nuclei are well-resolved for the first time. In the West nucleus, the major gas and dust emission extends out to 0.2\arcsec radius (74 pc); the central resolution element shows a strong peak in the dust emission but a factor 3 dip in the CO line emission. In this nucleus, the dust is apparently optically thick ($\tau_{\rm 2.6mm} \sim1$) at $\lambda = 2.6$ mm with a dust brightness temperature $\sim147$ K. The column of ISM at this nucleus is $\rm N_{H2} \geq 2\times10^{26}$ cm$^{-2}$, corresponding to $\sim$900 gr cm$^{-2}$. The East nucleus is more elongated with radial extent 0.3\arcsec or $\sim111$ pc. The derived kinematics of the nuclear disks provide a good fit to the line profiles, yielding the emissivity distributions, the rotation curves and velocity dispersions. In the West nucleus, there is evidence of a central Keplerian component requiring a central mass of $8\times10^8$ \msun. The intrinsic widths of the emission lines are $\Delta \rm v (FWHM)$ = 250 (West) and 120 (East) \kms. Given the very short dissipation timescales for turbulence ($\lesssim10^5$ yrs), we suggest that the line widths may be due to semi-coherent motions within the nuclear disks. The symmetry of the nuclear disk structures is impressive -- implying the merger timescale is significantly longer than the rotation period of the disks.

Mergers and the spin of the dark matter halo are factors traditionally believed to determine the morphology of galaxies within a $\Lambda$CDM cosmology. We study this hypothesis by considering approximately 18,000 central galaxies at $z=0$ with stellar masses $M_{\ast} = 10^{9}-10^{12} \, {\rm M}_{\odot}$ selected from the Illustris cosmological hydrodynamic simulation. The fraction of accreted stars -- which measures the importance of massive, recent and dry mergers -- increases steeply with galaxy stellar mass, from less than 5 per cent in dwarfs to 80 per cent in the most massive objects, and the impact of mergers on galaxy morphology increases accordingly. For galaxies with $M_{\ast} \gtrsim 10^{11} \, {\rm M}_{\odot}$, mergers have the expected effect: if gas-poor they promote the formation of spheroidal galaxies, whereas gas-rich mergers favour the formation and survivability of massive discs. This trend, however, breaks at lower masses. For objects with $M_{\ast} \lesssim 10^{11} \, {\rm M}_{\odot}$, mergers do not seem to play any significant role in determining the morphology, with accreted stellar fractions and mean merger gas fractions that are indistinguishable between spheroidal and disc-dominated galaxies. On the other hand, halo spin correlates with morphology primarily in the least massive objects in the sample ($M_{\ast} \lesssim 10^{10} \, {\rm M}_{\odot}$), but only weakly for galaxies above that mass. Our results support a scenario where (1) mergers play a dominant role in shaping the morphology of massive galaxies, (2) halo spin is important for the morphology of dwarfs, and (3) the morphology of medium-sized galaxies -- including the Milky Way -- shows little dependence on galaxy assembly history or halo spin, at least when these two factors are considered individually.

The interpretive power of the newest generation of large-volume hydrodynamical simulations of galaxy formation rests upon their ability to reproduce the observed properties of galaxies. In this second paper in a series, we employ bulge+disc decompositions of realistic dust-free galaxy images from the Illustris simulation in a consistent comparison with galaxies from the Sloan Digital Sky Survey (SDSS). Examining the size-luminosity relations of each sample, we find that galaxies in Illustris are roughly twice as large and $0.7$ magnitudes brighter on average than galaxies in the SDSS. The trend of increasing slope and decreasing normalization of size-luminosity as a function of bulge-fraction is qualitatively similar to observations. However, the size-luminosity relations of Illustris galaxies are quantitatively distinguished by higher normalizations and smaller slopes than for real galaxies. We show that this result is linked to a significant deficit of bulge-dominated galaxies in Illustris relative to the SDSS at stellar masses $\log\mathrm{M}_{\star}/\mathrm{M}_{\odot}\lesssim11$. We investigate this deficit by comparing bulge fraction estimates derived from photometry \emph{and} internal kinematics. We show that photometric bulge fractions are systematically lower than the kinematic fractions at low masses, but with increasingly good agreement as the stellar mass increases.

We use the Illustris simulation to study the relative contributions of in situ star formation and stellar accretion to the build-up of galaxies over an unprecedentedly wide range of masses ($M_{\ast} = 10^9-10^{12} \, {\rm M_{\odot}}$), galaxy types, environments, and assembly histories. We find that the `two-phase' picture of galaxy formation predicted by some models is a good approximation only for the most massive galaxies in our simulation -- namely, the stellar mass growth of galaxies below a few times $10^{11} \, {\rm M_{\odot}}$ is dominated by in situ star formation at all redshifts. The fraction of the total stellar mass of galaxies at $z=0$ contributed by accreted stars shows a strong dependence on galaxy stellar mass, ranging from about 10 per cent for Milky Way-sized galaxies to over 80 per cent for $M_{\ast} \approx 10^{12} \, {\rm M_{\odot}}$ objects, yet with a large galaxy-to-galaxy variation. At a fixed stellar mass, elliptical galaxies and those formed at the centres of younger haloes exhibit larger fractions of ex situ stars than disc-like galaxies and those formed in older haloes. On average, $\sim$50 per cent of the ex situ stellar mass comes from major mergers (stellar mass ratio $\mu > 1/4$), $\sim$20 per cent from minor mergers ($1/10 < \mu < 1/4$), $\sim$20 per cent from very minor mergers ($\mu < 1/10$), and $\sim$10 per cent from stars that were stripped from surviving galaxies (e.g. flybys or ongoing mergers). These components are spatially segregated, with in situ stars dominating the innermost regions of galaxies, and ex situ stars being deposited at larger galactocentric distances in order of decreasing merger mass ratio.

We present an image-based method for comparing the structural properties of galaxies produced in hydrodynamical simulations to real galaxies in the Sloan Digital Sky Survey. The key feature of our work is the introduction of extensive observational realism, such as object crowding, noise and viewing angle, to the synthetic images of simulated galaxies, so that they can be fairly compared to real galaxy catalogs. We apply our methodology to the dust-free synthetic image catalog of galaxies from the Illustris simulation at $z=0$, which are then fit with bulge+disc models to obtain morphological parameters. In this first paper in a series, we detail our methods, quantify observational biases, and present publicly available bulge+disc decomposition catalogs. We find that our bulge+disc decompositions are largely robust to the observational biases that affect decompositions of real galaxies. However, we identify a significant population of galaxies (roughly 30\% of the full sample) in Illustris that are prone to internal segmentation, leading to systematically reduced flux estimates by up to a factor of 6, smaller half-light radii by up to a factor of $\sim$ 2, and generally erroneous bulge-to-total fractions of (B/T)=0.

Recent spatially resolved observations of galaxies at z=0.6-3 reveal that high-redshift galaxies show complex kinematics and a broad distribution of gas-phase metallicity gradients. To understand these results, we use a suite of high-resolution cosmological zoom-in simulations from the Feedback in Realistic Environments (FIRE) project, which include physically motivated models of the multi-phase ISM, star formation, and stellar feedback. Our simulations reproduce the observed diversity of kinematic properties and metallicity gradients, broadly consistent with observations at z=0-3. Strong negative metallicity gradients only appear in galaxies with a rotating disk, but not all rotationally supported galaxies have significant gradients. Strongly perturbed galaxies with little rotation always have flat gradients. The kinematic properties and metallicity gradient of a high-redshift galaxy can vary significantly on short time-scales, associated with starburst episodes. Feedback from a starburst can destroy the gas disk, drive strong outflows, and flatten a pre-existing negative metallicity gradient. The time variability of a single galaxy is statistically similar to the entire simulated sample, indicating that the observed metallicity gradients in high-redshift galaxies reflect the instantaneous state of the galaxy rather than the accretion and growth history on cosmological time-scales. We find weak dependence of metallicity gradient on stellar mass and specific star formation rate (sSFR). Low-mass galaxies and galaxies with high sSFR tend to have flat gradients, likely due to the fact that feedback is more efficient in these galaxies. We argue that it is important to resolve feedback on small scales in order to produce the diverse metallicity gradients observed.

The velocity function is a fundamental observable statistic of the galaxy population, similarly impor- tant as the luminosity function, but much more difficult to measure. In this work we present the first directly measured circular velocity function that is representative between 60 < v_circ < 320 km/s for galaxies of all morphological types at a given rotation velocity. For the low mass galaxy population (60 < v_circ < 170 km/s), we use the HIPASS velocity function. For the massive galaxy population (170 < v_circ < 320 km/s), we use stellar circular velocities from the Calar Alto Legacy Integral Field Area Survey (CALIFA). In earlier work we obtained the measurements of circular velocity at the 80% light radius for 226 galaxies and demonstrated that the CALIFA sample can produce volume- corrected galaxy distribution functions. The CALIFA velocity function includes homogeneous velocity measurements of both late and early-type rotation-supported galaxies and has the crucial advantage of not missing gas-poor massive ellipticals that HI surveys are blind to. We show that both velocity functions can be combined in a seamless manner, as their ranges of validity overlap. The resulting observed velocity function is compared to velocity functions derived from cosmological simulations of the z = 0 galaxy population. We find that dark matter-only simulations show a strong mismatch with the observed VF. Hydrodynamic simulations fare better, but still do not fully reproduce observations.

We use particle data from the Illustris simulation, combined with individual kinematic constraints on the mass of the Milky Way (MW) at specific distances from the Galactic center, to infer the radial distribution of the MW's dark matter halo mass. Our method allows us to convert any constraint on the mass of the MW within a fixed distance to a full circular velocity profile to the MW's virial radius. As primary examples, we take two recent (and discrepant) measurements of the total mass within 50 kpc of the Galaxy and find they imply very different mass profiles and stellar masses for the Galaxy. The dark-matter-only version of the Illustris simulation enables us to compute the effects of galaxy formation on such constraints on a halo-by-halo basis; on small scales, galaxy formation enhances the density relative to dark-matter-only runs, while the total mass density is approximately 20% lower at large Galactocentric distances. We are also able to quantify how current and future constraints on the mass of the MW within specific radii will be reflected in uncertainties on its virial mass: even a measurement of M(<50 kpc) with essentially perfect precision still results in a 20% uncertainty on the virial mass of the Galaxy, while a future measurement of M(<100 kpc) with 10% errors would result in the same level of uncertainty. We expect that our technique will become even more useful as (1) better kinematic constraints become available at larger distances and (2) cosmological simulations provide even more faithful representations of the observable Universe.

We report the discovery of 14 Lyman-alpha blobs (LABs) at z~0.3, existing at least 4-7 billion years later in the Universe than all other LABs known. Their optical diameters are 20-70 kpc, and GALEX data imply Ly-alpha luminosities of (0.4-6.3)x10^43 erg/s. Contrary to high-z LABs, they live in low-density areas. They are ionized by AGN, suggesting that cold accretion streams as a power source must deplete between z=2 and z=0.3. We also show that transient AGN naturally explain the ionization deficits observed in many LABs: Their Ly-alpha and X-ray fluxes decorrelate below 10^6 years because of the delayed escape of resonantly scattering Ly-alpha photons. High Ly-alpha luminosities do not require currently powerful AGN, independent of obscuration. Chandra X-ray data reveal intrinsically weak AGN, confirming the luminous optical nebulae as impressive ionization echoes. For the first time, we also report mid-infrared thermal echoes from the dusty tori. We conclude that the AGN have faded by 3-4 orders of magnitude within the last 10^(4-5) years, leaving fossil UV, optical and thermal radiation behind. The host galaxies belong to the group of previously discovered Green Bean galaxies (GBs). Gemini optical imaging reveals smooth spheres, mergers, spectacular outflows and ionization cones. Because of their proximity and high flux densities, GBs are perfect targets to study AGN feedback, mode switching and the Ly-alpha escape. The fully calibrated, coadded optical FITS images are publicly available.

We quantify the impact that a variety of galactic and environmental properties have on the quenching of star formation. We collate a sample of $\sim$ 400,000 central and $\sim$ 100,000 satellite galaxies from the Sloan Digital Sky Survey Data Release 7 (SDSS DR7). Specifically, we consider central velocity dispersion ($\sigma_{c}$), stellar, halo, bulge and disk mass, local density, bulge-to-total ratio, group-centric distance and galaxy-halo mass ratio. We develop and apply a new statistical technique to quantify the impact on the quenched fraction ($f_{\rm Quench}$) of varying one parameter, while keeping the remaining parameters fixed. For centrals, we find that the $f_{\rm Quench} - \sigma_{c}$ relationship is tighter and steeper than for any other variable considered. We compare to the Illustris hydrodynamical simulation and the Munich semi-analytic model (L-Galaxies), finding that our results for centrals are qualitatively consistent with their predictions for quenching via radio-mode AGN feedback, hinting at the viability of this process in explaining our observational trends. However, we also find evidence that quenching in L-Galaxies is too efficient and quenching in Illustris is not efficient enough, compared to observations. For satellites, we find strong evidence that environment affects their quenched fraction at fixed central velocity dispersion, particularly at lower masses. At higher masses, satellites behave identically to centrals in their quenching. Of the environmental parameters considered, local density affects the quenched fraction of satellites the most at fixed central velocity dispersion.

Galaxy populations at different cosmic epochs are often linked together by comoving cumulative number density in observational studies. Many theoretical works, however, have shown that the number densities of tracked galaxy populations evolve in bulk and spread out over time. We present a number density method for linking progenitor and descendant galaxy populations which takes both of these effects into account. We define probability distribution functions that capture the evolution and dispersion of galaxy populations in comoving number density space, and use these functions to assign galaxies at one redshift $z_f$ probabilities of being progenitors or descendants of a galaxy population at another redshift $z_0$. These probabilities are then used as weights for calculating distributions of physical properties such as stellar mass, star formation rate, or velocity dispersion within the progenitor/descendant population. We demonstrate that this probabilistic method provides more accurate predictions for the evolution of physical properties then either the assumption of a constant number density or the assumption of an evolving number density in a bin of fixed width by comparing the predictions against galaxy populations directly tracked through a cosmological simulation. We find that the constant number density method performs most poorly at recovering galaxy properties, the evolving number method density slightly better, and the probabilistic number density method best of all. The improvement is present for predictions of both stellar mass as well as inferred quantities such as star formation rate and velocity dispersion which were not included in the number density fits. We demonstrate that this method can also be applied robustly and easily to observational data, and provide a code package for doing so.

Galaxy comoving number density is commonly used to forge progenitor/descendant links between observed galaxy populations at different epochs. However, this method breaks down in the presence of galaxy mergers, or when galaxies experience stochastic growth rates. We present a simple analytic framework to treat the physical processes that drive the evolution and diffusion of galaxies within comoving number density space. The evolution in mass rank order of a galaxy population with time is influenced by the galaxy coagulation rate and galaxy "mass rank scatter" rate. We quantify the relative contribution of these two effects to the mass rank order evolution. We show that galaxy coagulation is dominant at lower redshifts and stellar masses, while scattered growth rates dominate the mass rank evolution at higher redshifts and stellar masses. For a galaxy population at $10^{10} M_\odot$, coagulation has been the dominant effect since $z=2.2$, but a galaxy population at $10^{11} M_\odot$ was dominated by mass rank scatter until $z=0.6$. We show that although the forward and backward median number density evolution tracks are asymmetric, the backward median number density evolution can be obtained by convolving the descendant distribution function with progenitor relative abundances. We tabulate fits for the median number density evolution and scatter which can be applied to improve the way galaxy populations are linked in multi-epoch observational datasets.

We describe a statistical approach for measuring the influence that a galaxy's closest companion has on the galaxy's properties out to arbitrarily wide separations. We begin by identifying the closest companion for every galaxy in a large spectroscopic sample of Sloan Digital Sky Survey galaxies. We then characterize the local environment of each galaxy by using the number of galaxies within 2 Mpc and by determining the isolation of the galaxy pair from other neighbouring galaxies. We introduce a sophisticated algorithm for creating a statistical control sample for each galaxy, matching on stellar mass, redshift, local density and isolation. Unlike traditional studies of close galaxy pairs, this approach is effective in a wide range of environments, regardless of how far away the closest companion is (although a very distant closest companion is unlikely to have a measurable influence on the galaxy in question). We apply this methodology to measurements of galaxy asymmetry, and find that the presence of nearby companions drives a clear enhancement in galaxy asymmetries. The asymmetry excess peaks at the smallest projected separations (< 10 kpc), where the mean asymmetry is enhanced by a factor of 2.0 $\pm$ 0.2. Enhancements in mean asymmetry decline as pair separation increases, but remain statistically significant (1-2$\sigma$) out to projected separations of at least 50 kpc.

The concept of the cosmic web, viewing the Universe as a set of discrete galaxies held together by gravity, is deeply engrained in cosmology. Yet, little is known about the most effective construction and the characteristics of the underlying network. Here we explore seven network construction algorithms that use various galaxy properties, from their location, to their size and relative velocity, to assign a network to galaxy distributions provided by both simulations and observations. We find that a model relying only on spatial proximity offers the best correlations between the physical characteristics of the connected galaxies. We show that the properties of the networks generated from simulations and observations are identical, unveiling a deep universality of the cosmic web.

It is widely believed that the large discrepancy between the observed number of satellite galaxies and the predicted number of dark subhalos can be resolved via a variety of baryonic effects which suppress star formation in low mass halos.Supporting this hypothesis, numerous high resolution simulations with star formation, and associated feedback have been shown to reproduce the satellite luminosity function around Milky Way-mass simulated galaxies at redshift zero. However, a more stringent test of these models is their ability to simultaneously match the satellite luminosity functions of a range of host halo masses and redshifts. In this work we measure the luminosity function of faint (sub-Small Magellanic Cloud luminosity) satellites around hosts with stellar masses 10.5$<\log_{10}$M$_*$/M$_\odot<11.5$ to an unprecedented redshift of 1.5. This new measurement of the satellite luminosity function provides powerful new constraining power; we compare these results with predictions from four different simulations and show that although the models perform similarly over-all, no one model reproduces the satellite luminosity function reliably at all redshifts and host stellar masses. This result highlights the continued need for improvement in understanding the fundamental physics that governs satellite galaxy evolution.

We introduce a dust model for cosmological simulations implemented in the moving-mesh code AREPO and present a suite of cosmological hydrodynamical zoom-in simulations to study dust formation within galactic haloes. Our model accounts for the stellar production of dust, accretion of gas-phase metals onto existing grains, destruction of dust through local supernova activity, and dust driven by winds from star-forming regions. We find that accurate stellar and active galactic nuclei feedback is needed to reproduce the observed dust-metallicity relation and that dust growth largely dominates dust destruction. Our simulations predict a dust content of the interstellar medium which is consistent with observed scaling relations at $z = 0$, including scalings between dust-to-gas ratio and metallicity, dust mass and gas mass, dust-to-gas ratio and stellar mass, and dust-to-stellar mass ratio and gas fraction. We find that roughly two-thirds of dust at $z = 0$ originated from Type II supernovae, with the contribution from asymptotic giant branch stars below 20 per cent for $z \gtrsim 5$. While our suite of Milky Way-sized galaxies forms dust in good agreement with a number of key observables, it predicts a high dust-to-metal ratio in the circumgalactic medium, which motivates a more realistic treatment of thermal sputtering of grains and dust cooling channels.

Observations at low redshifts thus far fail to account for all of the baryons expected in the Universe according to cosmological constraints. A large fraction of the baryons presumably resides in a thin and warm-hot medium between the galaxies, where they are difficult to observe due to their low densities and high temperatures. Cosmological simulations of structure formation can be used to verify this picture and provide quantitative predictions for the distribution of mass in different large-scale structure components. Here we study the distribution of baryons and dark matter at different epochs using data from the Illustris simulation. We identify regions of different dark matter density with the primary constituents of large-scale structure, allowing us to measure mass and volume of haloes, filaments and voids. At redshift zero, we find that 49 % of the dark matter and 23 % of the baryons are within haloes more massive than the resolution limit of $2\times 10^8$ M$_\odot$. The filaments of the cosmic web host a further 45 % of the dark matter and 46 % of the baryons. The remaining 31 % of the baryons reside in voids. The majority of these baryons have been transported there through active galactic nuclei feedback. We note that the feedback model of Illustris is too strong for heavy haloes, therefore it is likely that we are overestimating this amount. Categorizing the baryons according to their density and temperature, we find that 17.8 % of them are in a condensed state, 21.6 % are present as cold, diffuse gas, and 53.9 % are found in the state of a warm-hot intergalactic medium.

We study the interaction of feedback from active galactic nuclei (AGN) and a multi-phase interstellar medium (ISM), in simulations including explicit stellar feedback, multi-phase cooling, accretion-disk winds, and Compton heating. We examine radii ~0.1-100 pc around a black hole (BH), where the accretion rate onto the BH is determined and where AGN-powered winds and radiation couple to the ISM. We conclude: (1) The BH accretion rate is determined by exchange of angular momentum between gas and stars in gravitational instabilities. This produces accretion rates ~0.03-1 Msun/yr, sufficient to power luminous AGN. (2) The gas disk in the galactic nucleus undergoes an initial burst of star formation followed by several Myrs where stellar feedback suppresses the star formation rate (SFR). (3) AGN winds injected at small radii with momentum fluxes ~L/c couple efficiently to the ISM and have dramatic effects on ISM properties within ~100 pc. AGN winds suppress the nuclear SFR by factors ~10-30 and BH accretion rate by factors ~3-30. They increase the outflow rate from the nucleus by factors ~10, consistent with observational evidence for galaxy-scale AGN-driven outflows. (4) With AGN feedback, the predicted column density distribution to the BH is consistent with observations. Absent AGN feedback, the BH is isotropically obscured and there are not enough optically-thin sightlines to explain Type-I AGN. A 'torus-like' geometry arises self-consistently as AGN feedback evacuates gas in polar regions.

We examine the stability of feedback-regulated star formation (SF) in galactic nuclei and contrast it to SF in extended discs. In galactic nuclei the dynamical time becomes shorter than the time over which feedback from young stars evolves. We argue analytically that the balance between stellar feedback and gravity is unstable in this regime. We study this using numerical simulations with pc-scale resolution and explicit stellar feedback taken from stellar evolution models. The nuclear gas mass, young stellar mass, and SFR within the central ~100 pc (the short-timescale regime) never reach steady-state, but instead go through dramatic, oscillatory cycles. Stars form until a critical surface density of young stars is present (such that feedback overwhelms gravity), at which point they begin to expel gas from the nucleus. Since the dynamical times are shorter than the stellar evolution times, the stars do not die as the gas is expelled, but continue to push, triggering a runaway quenching of star formation in the nucleus. However the expelled gas is largely not unbound from the galaxy, but goes into a galactic fountain which re-fills the nuclear region after the massive stars from the previous burst cycle have died off (~50 Myr timescale). On large scales (>1 kpc), the galaxy-scale gas content and SFR is more stable. We examine the consequences of this episodic nuclear star formation for the Kennicutt-Schmidt (KS) relation: while a tight KS relation exists on ~1 kpc scales in good agreement with observations, the scatter increases dramatically in smaller apertures centered on galactic nuclei.

Supermassive black hole (BH) mergers produce powerful gravitational wave (GW) emission. Asymmetry in this emission imparts a recoil kick to the merged BH, which can eject the BH from its host galaxy altogether. Recoiling BHs could be observed as offset active galactic nuclei (AGN). Several candidates have been identified, but systematic searches have been hampered by large uncertainties regarding their observability. By extracting merging BHs and host galaxy properties from the Illustris cosmological simulations, we have developed a comprehensive model for recoiling AGN. Here, for the first time, we model the effects of BH spin alignment and recoil dynamics based on the gas-richness of host galaxies. We predict that if BH spins are not highly aligned, seeing-limited observations could resolve offset AGN, making them promising targets for all-sky surveys. For randomly-oriented spins, less than about 10 spatially-offset AGN may be detectable in HST-COSMOS, and > 10^3 could be found with Pan-STARRS, LSST, Euclid, and WFIRST. Nearly a thousand velocity-offset AGN are predicted within the SDSS footprint; the rarity of large broad-line offsets among SDSS quasars is likely due in part to selection effects but suggests that spin alignment plays a role in suppressing recoils. Nonetheless, in our most physically motivated model where alignment occurs only in gas-rich mergers, hundreds of offset AGN should be found in all-sky surveys. Our findings strongly motivate a dedicated search for recoiling AGN.

Massive quiescent galaxies have much smaller physical sizes at high redshift than today. The strong evolution of galaxy size may be caused by progenitor bias, major and minor mergers, adiabatic expansion, and/or renewed star formation, but it is difficult to test these theories observationally. Herein, we select a sample of 35 massive, compact galaxies ($M_* = 1-3 \times 10^{11}$ M$_\odot$, $M_*/R^{1.5} > 10^{10.5}$ M$_\odot$/kpc$^{1.5}$) at $z=2$ in the cosmological hydrodynamical simulation Illustris and trace them forward to $z=0$ to uncover their evolution and identify their descendants. By $z=0$, the original factor of 3 difference in stellar mass spreads to a factor of 20. The dark matter halo masses similarly spread from a factor of 5 to 40. The galaxies' evolutionary paths are diverse: about half acquire an ex-situ envelope and are the core of a more massive descendant, a third survive undisturbed and gain very little mass, 15% are consumed in a merger with a more massive galaxy, and a small remainder are thoroughly mixed by major mergers. The galaxies grow in size as well as mass, and only $\sim$10% remain compact by $z=0$. The majority of the size growth is driven by the acquisition of ex-situ mass. The most massive galaxies at $z=0$ are the most likely to have compact progenitors, but this trend possesses significant dispersion which precludes a direct linkage to compact galaxies at $z=2$. The compact galaxies' merger rates are influenced by their $z=2$ environments, so that isolated or satellite compact galaxies (which are protected from mergers) are the most likely to survive to the present day.

Galaxy clusters contain a large population of low mass dwarf elliptical galaxies whose exact origin is unclear: their colors, structural properties and kinematics differ substantially from those of dwarf irregulars in the field. We use the Illustris cosmological simulation to study differences in the assembly paths of dwarf galaxies (3e8 < M_*/M_sun < 1e10) according to their environment. We find that cluster dwarfs achieve their maximum total and stellar mass on average ~ 8 and ~ 4.5 Gyr ago (or redshifts z = 1.0 and z = 0.4, respectively), around the time of infall into the clusters. In contrast, field dwarfs not subjected to environmental stripping, reach their maximum mass at redshift z = 0. This different assembly history naturally produces a color bimodality, with blue isolated dwarfs and redder cluster dwarfs exhibiting negligible star-formation today. The cessation of star formation happens over median times 3.5-5 Gyr depending on stellar mass, and shows a large scatter (~ 1-8 Gyr), with the lower values associated with starburst events that occur at infall through the virial radius or pericentric passages. We argue that such starbursts together with the early assembly of cluster dwarfs can provide a natural explanation for the higher specific frequency of globular clusters (GCs) in cluster dwarfs, as found observationally. We present a simple model for the formation and stripping of GCs that supports this interpretation. The origin of dwarf ellipticals in clusters is, therefore, consistent with an environmentally-driven evolution of field dwarf irregulars. However the z = 0 field analogs of cluster dwarf progenitors have today stellar masses a factor ~ 3 larger --a difference arising from the early truncation of star formation in cluster dwarfs.

We present the full public release of all data from the Illustris simulation project. Illustris is a suite of large volume, cosmological hydrodynamical simulations run with the moving-mesh code Arepo and including a comprehensive set of physical models critical for following the formation and evolution of galaxies across cosmic time. Each simulates a volume of (106.5 Mpc)^3 and self-consistently evolves five different types of resolution elements from a starting redshift of z=127 to the present day, z=0. These components are: dark matter particles, gas cells, passive gas tracers, stars and stellar wind particles, and supermassive black holes. This data release includes the snapshots at all 136 available redshifts, halo and subhalo catalogs at each snapshot, and two distinct merger trees. Six primary realizations of the Illustris volume are released, including the flagship Illustris-1 run. These include three resolution levels with the fiducial "full" baryonic physics model, and a dark matter only analog for each. In addition, we provide four distinct, high time resolution, smaller volume "subboxes". The total data volume is ~265 TB, including ~800 full volume snapshots and ~30,000 subbox snapshots. We describe the released data products as well as tools we have developed for their analysis. All data may be directly downloaded in its native HDF5 format. Additionally, we release a comprehensive, web-based API which allows programmatic access to search and data processing tasks. In both cases we provide example scripts and a getting-started guide in several languages: currently, IDL, Python, and Matlab. This paper addresses scientific issues relevant for the interpretation of the simulations, serves as a pointer to published and on-line documentation of the project, describes planned future additional data releases, and discusses technical aspects of the release.

We study how optical galaxy morphology depends on mass and star formation rate (SFR) in the Illustris Simulation. To do so, we measure automated galaxy structures in 10808 simulated galaxies at z=0 with stellar masses 10^9.7 < M_*/M_sun < 10^12.3. We add observational realism to idealized synthetic images and measure non-parametric statistics in rest-frame optical and near-IR images from four directions. We find that Illustris creates a morphologically diverse galaxy population, occupying the observed bulge strength locus and reproducing median morphology trends versus stellar mass, SFR, and compactness. Morphology correlates realistically with rotation, following classification schemes put forth by kinematic surveys. Type fractions as a function of environment agree roughly with data. These results imply that connections among mass, star formation, and galaxy structure arise naturally from models matching global star formation and halo occupation functions when simulated with accurate methods. This raises a question of how to construct experiments on galaxy surveys to better distinguish between models. We predict that at fixed halo mass near 10^12 M_sun, disc-dominated galaxies have higher stellar mass than bulge-dominated ones, a possible consequence of the Illustris feedback model. While Illustris galaxies at M_* ~ 10^11 M_sun have a reasonable size distribution, those at M_* ~ 10^10 M_sun have half-light radii larger than observed by a factor of two. Furthermore, at M_* ~ 10^10.5-10^11 M_sun, a relevant fraction of Illustris galaxies have distinct "ring-like" features, such that the bright pixels have an unusually wide spatial extent.

Comparisons between observational surveys and galaxy formation models find that the mass of dark matter haloes can largely explain galaxies' stellar mass. However, it remains uncertain whether additional environmental variables, generally referred to as assembly bias, are necessary to explain other galaxy properties. We use the Illustris Simulation to investigate the role of assembly bias in producing galactic conformity by considering 18,000 galaxies with $M_{stellar}$ > $2 \times 10^9$ $M_{\odot}$. We find a significant signal of galactic conformity: out to distances of about 10 Mpc, the mean red fraction of galaxies around redder galaxies is higher than around bluer galaxies at fixed stellar mass. Dark matter haloes exhibit an analogous conformity signal, in which the fraction of haloes formed at earlier times (old haloes) is higher around old haloes than around younger ones at fixed halo mass. A plausible interpretation of galactic conformity can be given as a combination of the halo conformity signal with the galaxy color-halo age relation: at fixed stellar mass, particularly toward the low-mass end, Illustris' galaxy colors correlate with halo age, with the reddest galaxies (often satellites) being preferentially found in the oldest haloes. In fact, we can explain the galactic conformity effect with a simple semi-empirical model, by assigning stellar mass based on halo mass (abundance matching) and by assigning galaxy color based on halo age (age matching). We investigate other interpretations for the galactic conformity, particularly its dependence on the isolation criterion and on the central-satellite information. Regarding comparison to observations, we conclude that the adopted selection/isolation criteria, projection effects, and stacking techniques can have a significant impact on the measured amplitude of the conformity signal.

Hydrodynamical simulations of galaxy formation such as the Illustris simulations have progressed to a state where they approximately reproduce the observed stellar mass function from high to low redshift. This in principle allows self-consistent models of reionization that exploit the accurate representation of the diffuse gas distribution together with the realistic growth of galaxies provided by these simulations, within a representative cosmological volume. In this work, we apply and compare two radiative transfer algorithms implemented in a GPU-accelerated code to the $106.5\,{\rm Mpc}$ wide volume of Illustris in postprocessing in order to investigate the reionization transition predicted by this model. We find that the first generation of galaxies formed by Illustris is just about able to reionize the universe by redshift $z\sim 7$, provided quite optimistic assumptions about the escape fraction and the resolution limitations are made. Our most optimistic model finds an optical depth of $\tau\simeq 0.065$, which is in very good agreement with recent Planck 2015 determinations. Furthermore, we show that moment-based approaches for radiative transfer with the M1 closure give broadly consistent results with our angular-resolved radiative transfer scheme. In our favoured fiducial model, 20\% of the hydrogen is reionized by redshift $z=9.20$, and this rapidly climbs to 80\% by redshift $z=6.92$. It then takes until $z=6.24$ before 99\% of the hydrogen is ionized. On average, reionization proceeds `inside-out' in our models, with a size distribution of reionized bubbles that progressively features regions of ever larger size while the abundance of small bubbles stays fairly constant.

The cumulative comoving number-density of galaxies as a function of stellar mass or central velocity dispersion is commonly used to link galaxy populations across different epochs. By assuming that galaxies preserve their number-density in time, one can infer the evolution of their properties, such as masses, sizes, and morphologies. However, this assumption does not hold in the presence of galaxy mergers or when rank ordering is broken owing to variable stellar growth rates. We present an analysis of the evolving comoving number density of galaxy populations found in the Illustris cosmological hydrodynamical simulation focused on the redshift range $0\leq z \leq 3$. Our primary results are as follows: 1) The inferred average stellar mass evolution obtained via a constant comoving number density assumption is systematically biased compared to the merger tree results at the factor of $\sim$2(4) level when tracking galaxies from redshift $z=0$ out to redshift $z=2(3)$; 2) The median number density evolution for galaxy populations tracked forward in time is shallower than for galaxy populations tracked backward in time; 3) A similar evolution in the median number density of tracked galaxy populations is found regardless of whether number density is assigned via stellar mass, stellar velocity dispersion, or dark matter halo mass; 4) Explicit tracking reveals a large diversity in galaxies' assembly histories that cannot be captured by constant number-density analyses; 5) The significant scatter in galaxy linking methods is only marginally reduced by considering a number of additional physical and observable galaxy properties as realized in our simulation. We provide fits for the forward and backward median evolution in stellar mass and number density and discuss implications of our analysis for interpreting multi-epoch galaxy property observations.

We study the properties of black holes and their host galaxies across cosmic time in the Illustris simulation. Illustris is a large scale cosmological hydrodynamical simulation which resolves a (106.5 Mpc)^3 volume with more than 12 billion resolution elements and includes state-of-the-art physical models relevant for galaxy formation. We find that the black hole mass density for redshifts z = 0 - 5 and the black hole mass function at z = 0 predicted by Illustris are in very good agreement with the most recent observational constraints. We show that the bolometric and hard X-ray luminosity functions of AGN at z = 0 and 1 reproduce observational data very well over the full dynamic range probed. Unless the bolometric corrections are largely underestimated, this requires radiative efficiencies to be on average low, epsilon_r <= 0.1, noting however that in our model radiative efficiencies are degenerate with black hole feedback efficiencies. Cosmic downsizing of the AGN population is in broad agreement with the findings from X-ray surveys, but we predict a larger number density of faint AGN at high redshifts than currently inferred. We also study black hole -- host galaxy scaling relations as a function of galaxy morphology, colour and specific star formation rate. We find that black holes and galaxies co-evolve at the massive end, but for low mass, blue and star-forming galaxies there is no tight relation with either their central black hole masses or the nuclear AGN activity.

We investigate the variation of the ratio of the equivalent widths of the FeII$\lambda$2600 line to the MgII$\lambda\lambda$2796,2803 doublet as a function of redshift in a large sample of absorption lines drawn from the JHU-SDSS Absorption Line Catalog. We find that despite large scatter, the observed ratio shows a trend where the equivalent width ratio $\mathcal{R}\equiv W_{\rm FeII}/W_{\rm MgII}$ decreases monotonically with increasing redshift $z$ over the range $0.55 \le z \le 1.90$. Selecting the subset of absorbers where the signal-to-noise ratio of the MgII equivalent width $W_{\rm MgII}$ is $\ge$3 and modeling the equivalent width ratio distribution as a gaussian, we find that the mean of the gaussian distribution varies as $\mathcal{R}\propto (-0.045\pm0.005)z$. We discuss various possible reasons for the trend. A monotonic trend in the Fe/Mg abundance ratio is predicted by a simple model where the abundances of Mg and Fe in the absorbing clouds are assumed to be the result of supernova ejecta and where the cosmic evolution in the SNIa and core-collapse supernova rates is related to the cosmic star-formation rate. If the trend in $\mathcal{R}$ reflects the evolution in the abundances, then it is consistent with the predictions of the simple model.

Cold accretion is a primary growth mechanism of simulated galaxies, yet observational evidence of "cold flows" at redshifts where they should be most efficient ($z=2$-4) is scarce. In simulations, cold streams manifest as Lyman-limit absorption systems (LLSs) with low heavy-element abundances similar to those of the diffuse IGM. Here we report on an abundance survey of 17 H I-selected LLSs at $z=3.2$-4.4 which exhibit no metal absorption in SDSS spectra. Using medium-resolution spectra obtained at Magellan, we derive ionization-corrected metallicities (or limits) with a Markov-Chain Monte Carlo sampling that accounts for the large uncertainty in $N_{\rm HI}$ measurements typical of LLSs. The metal-poor LLS sample overlaps with the IGM in metallicity and is best described by a model where $71^{+13}_{-11}\%$ are drawn from the IGM chemical abundance distribution. These represent roughly half of all LLSs at these redshifts, suggesting that 28-40$\%$ of the general LLS population at $z\sim3.7$ could trace unprocessed gas. An ancillary sample of ten LLSs without any a priori metal-line selection is best fit with $48^{+14}_{-12}\%$ of metallicities drawn from the IGM. We compare these results with regions of a moving-mesh simulation; the simulation finds only half as many baryons in IGM-metallicity LLSs, and most of these lie beyond the virial radius of the nearest galaxy halo. A statistically significant fraction of all LLSs have low metallicity and therefore represent candidates for accreting gas; large-volume simulations can establish what fraction of these candidates actually lie near galaxies and the observational prospects for detecting the presumed hosts in emission.

We have constructed merger trees for galaxies in the Illustris Simulation by directly tracking the baryonic content of subhalos. These merger trees are used to calculate the galaxy-galaxy merger rate as a function of descendant stellar mass, progenitor stellar mass ratio, and redshift. We demonstrate that the most appropriate definition for the mass ratio of a galaxy-galaxy merger consists in taking both progenitor masses at the time when the secondary progenitor reaches its maximum stellar mass. Additionally, we avoid effects from `orphaned' galaxies by allowing some objects to `skip' a snapshot when finding a descendant, and by only considering mergers which show a well-defined `infall' moment. Adopting these definitions, we obtain well-converged predictions for the galaxy-galaxy merger rate with the following main features, which are qualitatively similar to the halo-halo merger rate except for the last one: a strong correlation with redshift that evolves as $\sim (1+z)^{2.4-2.8}$, a power law with respect to mass ratio, and an increasing dependence on descendant stellar mass, which steepens significantly for descendant stellar masses greater than $\sim 2 \times 10^{11} \, {\rm M_{\odot}}$. These trends are consistent with observational constraints for medium-sized galaxies ($M_{\ast} \gtrsim 10^{10} \, {\rm M_{\odot}}$), but in tension with some recent observations of the close pair fraction for massive galaxies ($M_{\ast} \gtrsim 10^{11} \, {\rm M_{\odot}}$), which report a nearly constant or decreasing evolution with redshift. Finally, we provide a fitting function for the galaxy-galaxy merger rate which is accurate over a wide range of stellar masses, progenitor mass ratios, and redshifts.

New spectral line observations, obtained with the Jansky Very Large Array (VLA), of a sample of 34 galaxies in 17 close pairs are presented in this paper. The sample of galaxy pairs is selected to contain galaxies in close, major interactions (i.e., projected separations $<$30 kpc/h, and mass ratios less extreme than 4:1), while still having a sufficiently large angular separation that the VLA can spatially resolve both galaxies in the pair. Of the 34 galaxies, 17 are detected at $> 3\sigma$. We compare the HI gas fraction of the galaxies with the triggered star formation present in that galaxy. When compared to the star formation rates (SFRs) of non-pair galaxies matched in mass, redshift, and local environment, we find that the star formation enhancement is weakly positively correlated ($\sim 2.5\sigma$) with HI gas fraction. In order to help understand the physical mechanisms driving this weak correlation, we also present results from a small suite of binary galaxy merger simulations with varying gas fractions. The simulated galaxies indicate that larger initial gas fractions are associated with lower levels of interaction-triggered star formation (relative to an identical galaxy in isolation), but also show that high gas fraction galaxies have higher absolute SFRs prior to an interaction. We show that when interaction-driven SFR enhancements are calculated relative to a galaxy with an average gas fraction for its stellar mass, the relationship between SFR and initial gas fraction dominates over the SFR enhancements driven by the interaction. Simulated galaxy interactions that are matched in stellar mass but not in gas fraction, like our VLA sample, yield the same general positive correlation between SFR enhancement and gas fraction that we observe.

We investigate how the way galaxies acquire their gas across cosmic time in cosmological hydrodynamic simulations is modified by a comprehensive physical model for baryonic feedback processes. To do so, we compare two simulations -- with and without feedback -- both evolved with the moving mesh code AREPO. The feedback runs implement the full physics model of the Illustris simulation project, including star formation driven galactic winds and energetic feedback from supermassive blackholes. We explore: (a) the accretion rate of material contributing to the net growth of galaxies and originating directly from the intergalactic medium, finding that feedback strongly suppresses the raw, as well as the net, inflow of this "smooth mode" gas at all redshifts, regardless of the temperature history of newly acquired gas. (b) At the virial radius the temperature and radial flux of inflowing gas is largely unaffected at z=2. However, the spherical covering fraction of inflowing gas at 0.25 rvir decreases substantially, from more than 80% to less than 50%, while the rates of both inflow and outflow increase, indicative of recycling across this boundary. (c) The fractional contribution of smooth accretion to the total accretion rate is lower in the simulation with feedback, by roughly a factor of two across all redshifts. Moreover, the smooth component of gas with a cold temperature history, is entirely suppressed in the feedback run at z<1. (d) The amount of time taken by gas to cross from the virial radius to the galaxy -- the "halo transit time" -- increases in the presence of feedback by a factor of ~2-3, and is notably independent of halo mass. We discuss the possible implications of this invariance for theoretical models of hot halo gas cooling.

The presence of hot gaseous coronae around present-day massive spiral galaxies is a fundamental prediction of galaxy formation models. However, our observational knowledge remains scarce, since to date only four gaseous coronae were detected around spirals with massive stellar bodies ($\gtrsim2\times10^{11} \ \rm{M_{\odot}}$). To explore the hot coronae around lower mass spiral galaxies, we utilized Chandra X-ray observations of a sample of eight normal spiral galaxies with stellar masses of $(0.7-2.0)\times10^{11} \ \rm{M_{\odot}}$. Although statistically significant diffuse X-ray emission is not detected beyond the optical radii ($\sim20$ kpc) of the galaxies, we derive $3\sigma$ limits on the characteristics of the coronae. These limits, complemented with previous detections of NGC 1961 and NGC 6753, are used to probe the Illustris Simulation. The observed $3\sigma$ upper limits on the X-ray luminosities and gas masses exceed or are at the upper end of the model predictions. For NGC 1961 and NGC 6753 the observed gas temperatures, metal abundances, and electron density profiles broadly agree with those predicted by Illustris. These results hint that the physics modules of Illustris are broadly consistent with the observed properties of hot coronae around spiral galaxies. However, a shortcoming of Illustris is that massive black holes, mostly residing in giant ellipticals, give rise to powerful radio-mode AGN feedback, which results in under luminous coronae for ellipticals.

Massive, quiescent galaxies at high redshift have been found to be considerably more compact than galaxies of similar mass in the local universe. How these compact galaxies formed has yet to be determined, though several progenitor populations have been proposed. Here we investigate the formation processes and quantify the assembly histories of such galaxies in Illustris, a suite of hydrodynamical cosmological simulations encompassing a sufficiently large volume to include rare objects, while simultaneously resolving the internal structure of galaxies. We select massive (~10^11 solar masses) and compact (stellar half-mass radius < 2 kpc) galaxies from the simulation at z=2. Within the Illustris suite, we find that these quantities are not perfectly converged, but are reasonably reliable for our purposes. The resulting population is composed primarily of quiescent galaxies, but we also find several star-forming compact galaxies. The simulated compact galaxies are similar to observed galaxies in star formation activity and appearance. We follow their evolution at high redshift in the simulation and find that there are multiple pathways to form these compact galaxies, dominated by two mechanisms: (i) intense, centrally concentrated starbursts generally triggered by gas-rich major mergers between z~2-4, reducing the galaxies' half-mass radii by a factor of a few to below 2 kpc, and (ii) assembly at very early times when the universe was much denser; the galaxies formed compact and remained so until z~2.

We employ a suite of 75 simulations of galaxies in idealised major mergers (stellar mass ratio ~2.5:1), with a wide range of orbital parameters, to investigate the spatial extent of interaction-induced star formation. Although the total star formation in galaxy encounters is generally elevated relative to isolated galaxies, we find that this elevation is a combination of intense enhancements within the central kpc and moderately suppressed activity at large galacto-centric radii. The radial dependence of the star formation enhancement is stronger in the less massive galaxy than in the primary, and is also more pronounced in mergers of more closely aligned disc spin orientations. Conversely, these trends are almost entirely independent of the encounter's impact parameter and orbital eccentricity. Our predictions of the radial dependence of triggered star formation, and specifically the suppression of star formation beyond kph-scales, will be testable with the next generation of integral-field spectroscopic surveys.

Measurements of the neutral hydrogen gas content of a sample of 93 post-merger galaxies are presented, from a combination of matches to the ALFALFA.40 data release and new Arecibo observations. By imposing completeness thresholds identical to that of the ALFALFA survey, and by compiling a mass-, redshift- and environment-matched control sample from the public ALFALFA.40 data release, we calculate gas fraction offsets (Delta f_gas) for the post-mergers, relative to the control sample. We find that the post-mergers have HI gas fractions that are consistent with undisturbed galaxies. However, due to the relative gas richness of the ALFALFA.40 sample, from which we draw our control sample, our measurements of gas fraction enhancements are likely to be conservative lower limits. Combined with comparable gas fraction measurements by Fertig et al. in a sample of galaxy pairs, who also determine gas fraction offsets consistent with zero, we conclude that there is no evidence for significant neutral gas consumption throughout the merger sequence. From a suite of 75 binary merger simulations we confirm that star formation is expected to decrease the post-merger gas fraction by only 0.06 dex, even several Gyr after the merger. Moreover, in addition to the lack of evidence for gas consumption from gas fraction offsets, the observed HI detection fraction in the complete sample of post-mergers is twice as high as the controls, which suggests that the post-merger gas fractions may actually be enhanced. We demonstrate that a gas fraction enhancement in post-mergers, relative to a stellar mass-matched control sample, would indeed be the natural result of merging randomly drawn pairs from a parent population which exhibits a declining gas fraction with increasing stellar mass.

Galactic feedback strongly affects the way galactic environments are enriched. We examine this connection by performing a suite of cosmological hydrodynamic simulations, exploring a range of parameters based on the galaxy formation model developed in Vogelsberger et al. 2013 (henceforth V13). We examine the effects of AGN feedback, wind mass loading, wind specific energy, and wind metal-loading on the properties of the circumgalactic medium (CGM) of galaxies with $M_\text{halo} > 10^{11} M_\odot$. Note that while the V13 model was tuned to match observations including the stellar mass function, no explicit tuning was done for the CGM. The wind energy per unit outflow mass has the most significant effect on the CGM enrichment. High energy winds launch metals far beyond the virial radius. AGN feedback also has a significant effect, but only at $z < 3$. We compare to high redshift HI and CIV observations. All our simulations produce the observed number of Damped Lyman-$\alpha$ Absorbers. At lower column density, several of our simulations produce enough Lyman Limit Systems (LLS) $100$ kpc from the galaxy, but in all cases the LLS abundance falls off with distance faster than observations, with too few LLS at $200$ kpc. Further, in all models the CIV abundance drops off too sharply with distance, with too little CIV $100$-$200$ kpc from the galaxy. Higher energy wind models produce more extended CIV but also produce less stars, in tension with star-formation rate density observations. This highlights the fact that circumgalactic observations are a strong constraint on galactic feedback models.

Understanding the physical processes that drive star formation is a key challenge for galaxy formation models. In this article we study the tight correlation between the star formation rate (SFR) and stellar mass of galaxies at a given redshift, how halo growth influences star formation, and star formation histories of individual galaxies. We study these topics using Illustris, a state-of-the-art cosmological hydrodynamical simulation of galaxy formation. Illustris reproduces the observed relation (the star formation main sequence; SFMS) between SFR and stellar mass at redshifts z=0 and z=4, but at intermediate redshifts of z~2, the simulated SFMS has a significantly lower normalisation than reported by observations. The scatter in the relation is consistent with the observed scatter. However, the fraction of outliers above the SFR-stellar mass relation in Illustris is less than that observed. Galaxies with halo masses of ~10^{12} solar masses dominate the SFR density of the Universe, in agreement with the results of abundance matching. Furthermore, more-massive galaxies tend to form the bulk of their stars at high redshift, which indicates that `downsizing' occurs in Illustris. We also studied the star formation histories of individual galaxies, including the use of a principal component analysis decomposition. We find that for fixed stellar mass, galaxies that form earlier have more-massive black holes at z=0, indicating that star formation and black hole growth are tightly linked processes in Illustris. While many of the properties of normal star-forming galaxies are well-reproduced in the Illustris simulation, forming a realistic population of starbursts will likely require higher resolution and probably a more sophisticated treatment of star formation and feedback from stars and black holes.

We present our methods for generating a catalog of 7,000 synthetic images and 40,000 integrated spectra of redshift z = 0 galaxies from the Illustris Simulation. The mock data products are produced by using stellar population synthesis models to assign spectral energy distributions (SED) to each star particle in the galaxies. The resulting synthetic images and integrated SEDs therefore properly reflect the spatial distribution, stellar metallicity distribution, and star formation history of the galaxies. From the synthetic data products it is possible to produce monochromatic or color-composite images, perform SED fitting, classify morphology, determine galaxy structural properties, and evaluate the impacts of galaxy viewing angle. The main contribution of this paper is to describe the production, format, and composition of the image catalog that makes up the Illustris Simulation Obsevatory. As a demonstration of this resource, we derive galactic stellar mass estimates by applying the SED fitting code FAST to the synthetic galaxy products, and compare the derived stellar masses against the true stellar masses from the simulation. We find from this idealized experiment that systematic biases exist in the photometrically derived stellar mass values that can be reduced by using a fixed metallicity in conjunction with a minimum galaxy age restriction.

Observationally, the fraction of blue satellite galaxies decreases steeply with host halo mass, and their radial distribution around central galaxies is significantly shallower in massive (M_* >10e11M_sun) than in Milky Way like systems. Theoretical models, based primarily on semi-analytical techniques, have had a long-standing problem with reproducing these trends, instead predicting too few blue satellites in general but also estimating a radial distribution that is too shallow, regardless of primary mass. In this Letter, we use the Illustris cosmological simulation to study the properties of satellite galaxies around isolated primaries. For the first time, we find good agreement between theory and observations. We identify the main source of this success relative to earlier work to be a consequence of the large gas contents of satellites at infall, a factor ~5-10 times larger than in semi-analytical models. Because of their relatively large gas reservoirs, satellites can continue to form stars long after infall, with a typical timescale for star-formation to be quenched ~2 Gyr in groups but more than ~5 Gyr for satellites around Milky Way like primaries. The gas contents we infer are consistent with z=0 observations of HI gas in galaxies, although we find large discrepancies among reported values in the literature. A testable prediction of our model is that the gas-to-stellar mass ratio of satellite progenitors should vary only weakly with cosmic time.

We introduce the Illustris Project, a series of large-scale hydrodynamical simulations of galaxy formation. The highest resolution simulation, Illustris-1, covers a volume of $(106.5\,{\rm Mpc})^3$, has a dark mass resolution of ${6.26 \times 10^{6}\,{\rm M}_\odot}$, and an initial baryonic matter mass resolution of ${1.26 \times 10^{6}\,{\rm M}_\odot}$. At $z=0$ gravitational forces are softened on scales of $710\,{\rm pc}$, and the smallest hydrodynamical gas cells have an extent of $48\,{\rm pc}$. We follow the dynamical evolution of $2\times 1820^3$ resolution elements and in addition passively evolve $1820^3$ Monte Carlo tracer particles reaching a total particle count of more than $18$ billion. The galaxy formation model includes: primordial and metal-line cooling with self-shielding corrections, stellar evolution, stellar feedback, gas recycling, chemical enrichment, supermassive black hole growth, and feedback from active galactic nuclei. At $z=0$ our simulation volume contains about $40,000$ well-resolved galaxies covering a diverse range of morphologies and colours including early-type, late-type and irregular galaxies. The simulation reproduces reasonably well the cosmic star formation rate density, the galaxy luminosity function, and baryon conversion efficiency at $z=0$. It also qualitatively captures the impact of galaxy environment on the red fractions of galaxies. The internal velocity structure of selected well-resolved disk galaxies obeys the stellar and baryonic Tully-Fisher relation together with flat circular velocity curves. In the well-resolved regime the simulation reproduces the observed mix of early-type and late-type galaxies. Our model predicts a halo mass dependent impact of baryonic effects on the halo mass function and the masses of haloes caused by feedback from supernova and active galactic nuclei.

We examine the abundance, clustering and metallicity of Damped Lyman-alpha Absorbers (DLAs) in a suite of hydrodynamic cosmological simulations using the moving mesh code AREPO. We incorporate models of supernova and AGN feedback, as well as molecular hydrogen formation. We compare our simulations to the column density distribution function at $z=3$, the total DLA abundance at $z=2-4$, the measured DLA bias at $z=2.3$ and the DLA metallicity distribution at $z=2-4$. Our preferred models produce populations of DLAs in good agreement with most of these observations. The exception is the DLA abundance at $z < 3$, which we show requires stronger feedback in $10^{11-12} \, h^{-1} M_\odot$ mass halos. While the DLA population probes a wide range of halo masses, we find the cross-section is dominated by halos of mass $10^{10} - 10^{11} \, h^{-1} M_\odot$ and virial velocities $50 - 100 \;\mathrm{km/s}$. The simulated DLA population has a linear theory bias of $1.7$, whereas the observations require $2.17 \pm 0.2$. We show that non-linear growth increases the bias in our simulations to $2.3$ at $k=1\; \mathrm{Mpc/}h$, the smallest scale observed. The scale-dependence of the bias is, however, very different in the simulations compared against the observations. We show that, of the observations we consider, the DLA abundance and column density function provide the strongest constraints on the feedback model.

We use the Illustris Simulations to gain insight into the build-up of the outer, low-surface brightness regions which surround galaxies. We characterize the stellar haloes by means of the logarithmic slope of the spherically-averaged stellar density profiles, alphaSTARS at z=0, and we relate these slopes to the properties of the underlying Dark-Matter (DM) haloes, their central galaxies, and their assembly histories. We analyze a sample of ~5,000 galaxies resolved with more than 5x10^4 particles each, and spanning a variety of morphologies and halo masses (3x10^11 < Mvir < 10^14 Msun). We find a strong trend between stellar halo slope and total halo mass, where more massive objects have shallower stellar haloes than the less massive ones (-5.5 \pm 0.5 < alphaSTARS <-3.5 \pm 0.2 in the studied mass range). At fixed halo mass, we show that disk-like, blue, young, and more massive galaxies are surrounded by significantly steeper stellar haloes than elliptical, red, older, and less massive galaxies. Overall, the stellar density profiles fall off much more steeply than the underlying DM, and no clear trend holds between stellar slope and DM halo concentration. However, DM haloes which formed more recently, or which accreted larger fractions of stellar mass from infalling satellites, exhibit shallower stellar haloes than their older analogs with similar masses, by up to Delta(alphaSTARS) ~ 0.5-0.7. Our findings, combined with the most recent measurements of the strikingly different stellar power-law indexes for M31 and the Milky Way, appear to favour a massive M31, and a Milky Way characterized by a much quieter accretion history over the past 10 Gyrs than its companion.

We present an overview of galaxy evolution across cosmic time in the Illustris Simulation. Illustris is an N-body/hydrodynamical simulation that evolves 2*1820^3 resolution elements in a (106.5Mpc)^3 box from cosmological initial conditions down to z=0 using the AREPO moving-mesh code. The simulation uses a state-of-the-art set of physical models for galaxy formation that was tuned to reproduce the z=0 stellar mass function and the history of the cosmic star-formation rate density. We find that Illustris successfully reproduces a plethora of observations of galaxy populations at various redshifts, for which no tuning was performed, and provide predictions for future observations. In particular, we discuss (a) the buildup of galactic mass, showing stellar mass functions and the relations between stellar mass and halo mass from z=7 to z=0, (b) galaxy number density profiles around massive central galaxies out to z=4, (c) the gas and total baryon content of both galaxies and their halos for different redshifts, and as a function of mass and radius, and (d) the evolution of galaxy specific star-formation rates up to z=8. In addition, we (i) present a qualitative analysis of galaxy morphologies from z=5 to z=0, for the stellar as well as the gaseous components, and their appearance in HST mock observations, (ii) follow galaxies selected at z=2 to their z=0 descendants, and quantify their growth and merger histories, and (iii) track massive z=0 galaxies to high redshift and study their joint evolution in star-formation activity and compactness. We conclude with a discussion of several disagreements with observations, and lay out possible directions for future research.

Galaxy mergers have been investigated for decades using smoothed particle hydrodynamics (SPH), but recent work highlighting inaccuracies inherent in the traditional SPH technique calls into question the reliability of previous studies. We explore this issue by comparing a suite of Gadget-3 SPH simulations of idealised (i.e., non-cosmological) isolated discs and galaxy mergers with otherwise identical calculations performed using the moving-mesh code Arepo. When black hole (BH) accretion and active galactic nucleus (AGN) feedback are not included, the star formation histories (SFHs) obtained from the two codes agree well. When BHs are included, the code- and resolution-dependent variations in the SFHs are more significant, but the agreement is still good, and the stellar mass formed over the course of a simulation is robust to variations in the numerical method. During a merger, the gas morphology and phase structure are initially similar prior to the starburst phase. However, once a hot gaseous halo has formed from shock heating and AGN feedback (when included), the agreement is less good. In particular, during the post-starburst phase, the SPH simulations feature more prominent hot gaseous haloes and spurious clumps, whereas with Arepo, gas clumps and filaments are less apparent and the hot halo gas can cool more efficiently. We discuss the origin of these differences and explain why the SPH technique yields trustworthy results for some applications (such as the idealised isolated disc and galaxy merger simulations presented here) but not others (e.g., gas flows onto galaxies in cosmological hydrodynamical simulations).

Previous simulations of the growth of cosmic structures have broadly reproduced the 'cosmic web' of galaxies that we see in the Universe, but failed to create a mixed population of elliptical and spiral galaxies due to numerical inaccuracies and incomplete physical models. Moreover, because of computational constraints, they were unable to track the small scale evolution of gas and stars to the present epoch within a representative portion of the Universe. Here we report a simulation that starts 12 million years after the Big Bang, and traces 13 billion years of cosmic evolution with 12 billion resolution elements in a volume of $(106.5\,{\rm Mpc})^3$. It yields a reasonable population of ellipticals and spirals, reproduces the distribution of galaxies in clusters and statistics of hydrogen on large scales, and at the same time the metal and hydrogen content of galaxies on small scales.

We present a new comprehensive model of the physics of galaxy formation designed for large-scale hydrodynamical simulations of structure formation using the moving mesh code AREPO. Our model includes primordial and metal line cooling with self-shielding corrections, stellar evolution and feedback processes, gas recycling, chemical enrichment, a novel subgrid model for the metal loading of outflows, black hole (BH) seeding, BH growth and merging procedures, quasar- and radio-mode feedback, and a prescription for radiative electro-magnetic (EM) feedback from active galactic nuclei (AGN). The metal mass loading of outflows can be adjusted independently of the wind mass loading. This is required to simultaneously reproduce the stellar mass content of low mass haloes and their gas oxygen abundances. Radiative EM AGN feedback is implemented assuming an average spectral energy distribution and a luminosity-dependent scaling of obscuration effects. This form of feedback suppresses star formation more efficiently than continuous thermal quasar-mode feedback alone, but is less efficient than mechanical radio-mode feedback in regulating star formation in massive haloes. We contrast simulation predictions for different variants of our galaxy formation model with key observations. Our best match model reproduces, among other things, the cosmic star formation history, the stellar mass function, the stellar mass - halo mass relation, g-, r-, i-, z-band SDSS galaxy luminosity functions, and the Tully-Fisher relation. We can achieve this success only if we invoke very strong forms of stellar and AGN feedback such that star formation is adequately reduced in both low and high mass systems. In particular, the strength of radio-mode feedback needs to be increased significantly compared to previous studies to suppress efficient cooling in massive, metal-enriched haloes.

Galactic winds are a key physical mechanism for understanding galaxy formation and evolution, yet empirical and theoretical constraints for the character of winds are limited and discrepant. Recent empirical models find that local star-forming galaxies have a deficit of oxygen that scales with galaxy stellar mass. The oxygen deficit provides unique empirical constraints on the magnitude of mass loss, composition of outflowing material and metal reaccretion onto galaxies. We formulate the oxygen deficit constraints so they may be easily implemented into theoretical models of galaxy evolution. We parameterize an effective metal loading factor which combines the uncertainties of metal outflows and metal reaccretion into a single function of galaxy virial velocity. We determine the effective metal loading factor by forward-fitting the oxygen deficit. The effective metal loading factor we derive has important implications for the implementation of mass loss in models of galaxy evolution.

We develop a dust efflux model of radiation pressure acting on dust grains which successfully reproduces the relation between stellar mass, dust opacity and star formation rate observed in local star-forming galaxies. The dust content of local star-forming galaxies is set by the competition between the physical processes of dust production and dust loss in our model. The dust loss rate is proportional to the dust opacity and star formation rate. Observations of the relation between stellar mass and star formation rate at several epochs imply that the majority of local star-forming galaxies are best characterized as having continuous star formation histories. Dust loss is a consequence of sustained interaction of dust with the radiation field generated by continuous star formation. Dust efflux driven by radiation pressure rather than dust destruction offers a more consistent physical interpretation of the dust loss mechanism. By comparing our model results with the observed relation between stellar mass, dust extinction and star formation rate in local star-forming galaxies we are able to constrain the timescale and magnitude of dust loss. The timescale of dust loss is long and therefore dust is effluxed in a "Slow Flow". Dust loss is modest in low mass galaxies but massive galaxies may lose up to 70~80% of their dust over their lifetime. Our Slow Flow model shows that mass loss driven by dust opacity and star formation may be an important physical process for understanding normal star-forming galaxy evolution.

We use the Millennium Simulation, and an abundance-matching framework, to investigate the dynamical behaviour of galaxy pairs embedded in a cosmological context. Our main galaxy-pair sample, selected to have separations under 250 kpc/h, consists of over 1.3 million pairs at redshift z = 0, with stellar masses greater than 10^9 Msun, probing mass ratios down to 1:1000. We use dark matter halo membership and energy to classify our galaxy pairs. In terms of halo membership, central-satellite pairs tend to be in isolation (in relation to external more massive galaxies), are energetically- bound to each other, and are also weakly-bound to a neighbouring massive galaxy. Satellite-satellite pairs, instead, inhabit regions in close proximity to a more massive galaxy, are energetically-unbound, and are often bound to that neighbour. We find that 60% of our paired galaxies are bound to both their companion and to a third external object. Moreover, only 9% of our pairs resemble the kind of systems described by idealised binary merger simulations in complete isolation. In sum, we demonstrate the importance of properly connecting galaxy pairs to the rest of the Universe.

We use pair and environmental classifications of ~ 211,000 star-forming galaxies from the Sloan Digital Sky Survey, along with a suite of merger simulations, to investigate the enhancement of star formation as a function of separation in galaxy pairs. Using a new technique for distinguishing between the influence of nearby neighbours and larger scale environment, we find a clear enhancement in star formation out to projected separations of ~ 150 kpc, beyond which there is no net enhancement. We find the strongest enhancements at the smallest separations (especially < 20 kpc), consistent with earlier work. Similar trends are seen in the simulations, which indicate that the strongest enhancements are produced in highly disturbed systems approaching final coalescence, whereas the more modest enhancements seen at wider separations are the result of starburst activity triggered at first pericentre passage, which persists as the galaxies move to larger separations. The absence of any net enhancement beyond 150 kpc provides reassurance that the detected enhancements are due to galaxy-galaxy interactions, rather than larger scale environmental effects or potential pair selection biases. A rough census indicates that 66 per cent of the enhanced star formation in our pair sample occurs at separations > 30 kpc. We conclude that significant interaction-induced star formation is not restricted to merger remnants or galaxies with close companions; instead, a larger population of wider separation pairs exhibit enhanced star formation due to recent close encounters.

We present a multi-epoch analysis of the galaxy populations formed within the cosmological hydrodynamical simulations presented in Vogelsberger et al. (2013). These simulations explore the performance of a recently implemented feedback model which includes primordial and metal line radiative cooling with self-shielding corrections; stellar evolution with associated mass loss and chemical enrichment; feedback by stellar winds; black hole seeding, growth and merging; and AGN quasar- and radio-mode heating with a phenomenological prescription for AGN electro-magnetic feedback. We illustrate the impact of the model parameter choices on the resulting simulated galaxy population properties at high and intermediate redshifts. We demonstrate that our scheme is capable of producing galaxy populations that broadly reproduce the observed galaxy stellar mass function extending from redshift z=0 to z=3. We also characterise the evolving galactic B-band luminosity function, stellar mass to halo mass ratio, star formation main sequence, Tully-Fisher relation, and gas-phase mass-metallicity relation and confront them against recent observational estimates. This detailed comparison allows us to validate elements of our feedback model, while also identifying areas of tension that will be addressed in future work.

We compare the structural properties of galaxies formed in cosmological simulations using the smoothed particle hydrodynamics (SPH) code GADGET with those using the moving-mesh code AREPO. Both codes employ identical gravity solvers and the same sub-resolution physics but use very different methods to track the hydrodynamic evolution of gas. This permits us to isolate the effects of the hydro solver on the formation and evolution of galactic gas disks in GADGET and AREPO haloes with comparable numerical resolution. In a matching sample of GADGET and AREPO haloes we fit simulated gas disks with exponential profiles. We find that the cold gas disks formed using the moving mesh approach have systematically larger disk scale lengths and higher specific angular momenta than their GADGET counterparts across a wide range in halo masses. For low mass galaxies differences between the properties of the simulated galaxy disks are caused by an insufficient number of resolution elements which lead to the artificial angular momentum transfer in our SPH calculation. We however find that galactic disks formed in massive halos, resolved with 10^6 particles/cells, are still systematically smaller in the GADGET run by a factor of ~2. The reasons for this are: 1) The excessive heating of haloes close to the cooling radius due to spurious dissipation of the subsonic turbulence in GADGET; and 2) The efficient delivery of low angular momentum gaseous blobs to the bottom of the potential well. While this large population of gaseous blobs in GADGET originates from the filaments which are pressure confined and fragment due to the SPH surface tension while infalling into hot halo atmospheres, it is essentially absent in the moving mesh calculation, clearly indicating numerical rather than physical origin of the blob material.

We present a sample of 1899 galaxies with a close companion taken from the SDSS DR7. The galaxy pairs are selected to have velocity differences < 300 km/s, projected separations (rp) < 80 kpc/h, mass ratios between 0.1 and 10, and robust measurements of star formation rates and gas-phase metallicities. We match the galaxies in total stellar mass, redshift, and local density to a set of 10 control galaxies per pair galaxy. For each pair galaxy we can therefore calculate the statistical change in star formation rate (SFR) and metallicity associated with the interaction process. Relative to the control sample, we find that galaxies in pairs show typical SFR enhancements that are, on average, 60% higher than the control sample at rp < 30 kpc/h. In addition, the pairs demonstrate more modest SFR enhancements of ~30% out to at least 80 kpc/h (the widest separations in our sample). Galaxies in both major and minor mergers show significant SFR enhancements at all rp, although the strongest starbursts (with SFR enhancements of a factor of ~10) appear to be found only in the major mergers. For the first time, we are also able to trace the metallicity changes in galaxy pairs as a function of projected separation. The metallicity is generally diluted in galaxy pairs by ~0.02 dex, with an average metallicity decrement of -0.03 dex at the smallest separations. The SFR and metallicity trends with projected separation are interpreted through a comparison with theoretical models. These simulations indicate that the peak in SFR enhancements at small separations is due to systems near the end of the merger process. The extended plateau in SFR enhancements out to at least 80 kpc/h is dominated by galaxies that have made a pericentric passage and are now experiencing triggered star formation on their trajectory towards apogalacticon, or on a subsequent close approach. (Abridged)

We present for the first time metallicity maps generated using data from the Wide Field Spectrograph (WiFeS) on the ANU 2.3m of 9 Luminous Infrared Galaxies (LIRGs) and discuss the abundance gradients and distribution of metals in these systems. We have carried out optical integral field spectroscopy (IFS) of several several LIRGs in various merger phases to investigate the merger process. In a major merger of two spiral galaxies with preexisting disk abundance gradients, the changing distribution of metals can be used as a tracer of gas flows in the merging system as low metallicity gas is transported from the outskirts of each galaxy to their nuclei. We employ this fact to probe merger properties by using the emission lines in our IFS data to calculate the gas-phase metallicity in each system. We create abundance maps and subsequently derive a metallicity gradient from each map. We compare our measured gradients to merger stage as well as several possible tracers of merger progress and observed nuclear abundances. We discuss our work in the context of previous abundance gradient observations and compare our results to new galaxy merger models which trace metallicity gradient. Our results agree with the observed flattening of metallicity gradients as a merger progresses. We compare our results with new theoretical predictions that include chemical enrichment. Our data show remarkable agreement with these simulations.

Nuclear inflows of metal-poor interstellar gas triggered by galaxy interactions can account for the systematically lower central oxygen abundances observed in local interacting galaxies. Here, we investigate the metallicity evolution of a large set of simulations of colliding galaxies. Our models include cooling, star formation, feedback, and a new stochastic method for tracking the mass recycled back to the interstellar medium from stellar winds and supernovae. We study the influence of merger-induced inflows, enrichment, gas consumption, and galactic winds in determining the nuclear metallicity. The central metallicity is primarily a competition between the inflow of low-metallicity gas and enrichment from star formation. An average depression in the nuclear metallicity of ~0.07 is found for gas-poor disk-disk interactions. Gas-rich disk-disk interactions, on the other hand, typically have an enhancement in the central metallicity that is positively correlated with the gas content. The simulations fare reasonably well when compared to the observed mass-metallicity and separation-metallicity relationships, but further study is warranted.