Source author record

Shy Genel

Shy Genel appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.CO astro-ph.GA astro-ph.IM astro-ph.HE Computer Vision Machine Learning physics.flu-dyn physics.hist-ph

Catalog footprint

What is connected

59works

8topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

preprint2026arXiv

Cosmological back-reaction of baryons on dark matter in the CAMELS simulations

Baryonic processes such as radiative cooling and feedback from massive stars and active galactic nuclei (AGN) directly redistribute baryons in the Universe but also indirectly redistribute dark matter due to changes in the gravitational potential. In this work, we investigate this "back-reaction" of baryons on dark matter using thousands of cosmological hydrodynamic simulations from the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project, including parameter variations in the SIMBA, IllustrisTNG, ASTRID, and Swift-EAGLE galaxy formation models. Matching haloes to corresponding N-body (dark matter-only) simulations, we find that virial masses decrease owing to the ejection of baryons by feedback. Relative to N-body simulations, halo profiles show an increased dark matter density in the center (due to radiative cooling) and a decrease in density farther out (due to feedback), with both effects being strongest in SIMBA (> 450% increase at r < 0.01 Rvir). The clustering of dark matter strongly responds to changes in baryonic physics, with dark matter power spectra in some simulations from each model showing as much as 20% suppression or increase in power at k ~ 10 h/Mpc relative to N-body simulations. We find that the dark matter back-reaction depends intrinsically on cosmology (Omega_m and sigma_8) at fixed baryonic physics, and varies strongly with the details of the feedback implementation. These results emphasize the need for marginalizing over uncertainties in baryonic physics to extract cosmological information from weak lensing surveys as well as their potential to constrain feedback models in galaxy evolution.

preprint2022arXiv

Boundaries of chaos and determinism in the cosmos

According to the standard model of cosmology, the arrangement of matter in the cosmos on scales much larger than galaxies is entirely specified by the initial conditions laid down during inflation. But zooming in by dozens of orders of magnitude to microscopic (and human?) scales, quantum randomness reigns, independent of the initial conditions. Where is the boundary of determinism, and how does that interplay with chaos? Here, we make a first attempt at answering this question in an astronomical context, including currently understood processes. The boundary is a function, at least, of length scale, position, and matter type (dark matter being more simply predictable). In intergalactic voids, the primordial pattern of density fluctuations is largely preserved. But we argue that within galaxies, the conditions are at minimum chaotic, and may even be influenced by non-primordial information, or randomness independent of the initial conditions. Randomness could be supplied by events such as supernovae and jets from active galactic nuclei (AGN) and other accretion disks, which, with the help of chaotic dynamics, could broadcast any possible microscopic randomness to larger scales, eventually throughout a galaxy. This may be generated or amplified by a recently investigated process called spontaneous stochasticity, or effective randomness in turbulent systems arising from arbitrarily small perturbations.

preprint2022arXiv

Breaking baryon-cosmology degeneracy with the electron density power spectrum

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

preprint2022arXiv

Cosmology with one galaxy?

Galaxies can be characterized by many internal properties such as stellar mass, gas metallicity, and star-formation rate. We quantify the amount of cosmological and astrophysical information that the internal properties of individual galaxies and their host dark matter halos contain. We train neural networks using hundreds of thousands of galaxies from 2,000 state-of-the-art hydrodynamic simulations with different cosmologies and astrophysical models of the CAMELS project to perform likelihood-free inference on the value of the cosmological and astrophysical parameters. We find that knowing the internal properties of a single galaxy allow our models to infer the value of $Ω_{\rm m}$, at fixed $Ω_{\rm b}$, with a $\sim10\%$ precision, while no constraint can be placed on $σ_8$. Our results hold for any type of galaxy, central or satellite, massive or dwarf, at all considered redshifts, $z\leq3$, and they incorporate uncertainties in astrophysics as modeled in CAMELS. However, our models are not robust to changes in subgrid physics due to the large intrinsic differences the two considered models imprint on galaxy properties. We find that the stellar mass, stellar metallicity, and maximum circular velocity are among the most important galaxy properties to determine the value of $Ω_{\rm m}$. We believe that our results can be explained taking into account that changes in the value of $Ω_{\rm m}$, or potentially $Ω_{\rm b}/Ω_{\rm m}$, affect the dark matter content of galaxies. That effect leaves a distinct signature in galaxy properties to the one induced by galactic processes. Our results suggest that the low-dimensional manifold hosting galaxy properties provides a tight direct link between cosmology and astrophysics.

preprint2022arXiv

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

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

preprint2022arXiv

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

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

preprint2022arXiv

The MUSE Hubble Ultra Deep Field Survey XVI. The angular momentum of low-mass star-forming galaxies. A cautionary tale and insights from TNG50

We investigate the specific angular momentum (sAM) $ j(<r)$ profiles of intermediate redshift ($0.4<z<1.4$) star-forming galaxies (SFGs) in the relatively unexplored regime of low masses (down to $M_\star\sim 10^8$M$_{\odot}$), and small sizes (down to $R_{\rm e}\sim 1.5$ kpc) and characterize the sAM scaling relation and its redshift evolution. We have developed a 3D methodology to constrain sAM profiles of the star-forming gas using a forward modeling approach with \galpak{} that incorporates the effects of beam smearing, yielding the intrinsic morpho-kinematic properties even with limited spatial resolution data. Using mock observations from the TNG50 simulation, we find that our 3D methodology robustly recovers the star formation rate (SFR)-weighted $j(<r)$ profiles down to low effective signal-to-noise ratio (SNR) of $\gtrapprox3$. We applied our methodology blindly to a sample of 494 \OII{}-selected SFGs in the MUSE Ultra Deep Field (UDF) 9~arcmin$^2$ mosaic data, covering the unexplored $8<\log M_*/$M$_{\odot}<9$ mass range. We find that the (SFR-weighted) sAM relation follows $j\propto M_\star^α$ with an index $α$ varying from $α=0.3$ to $α=0.5$, from $\log M_\star/$M$_{\odot}=8$ to $\log M_*/$M$_{\odot}=10.5$. The UDF sample supports a redshift evolution consistent with the $(1+z)^{-0.5}$ expectation from a Universe in expansion. The scatter of the sAM sequence is a strong function of the dynamical state with $\log j|_{M_*}\propto 0.65 \times \log(V_{\rm max}/σ)$ where $σ$ is the velocity dispersion at $2 R_{\rm e}$. In TNG50, SFGs also form a $j-M_{\star}-(V/σ)$ plane but it correlates more with galaxy size than with morphological parameters. Our results suggest that SFGs might experience a dynamical transformation before their morphological transformation to becoming passive via either merging or secular evolution.

preprint2021arXiv

Finding universal relations in subhalo properties with artificial intelligence

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

preprint2021arXiv

The CAMELS Multifield Dataset: Learning the Universe's Fundamental Parameters with Artificial Intelligence

We present the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) Multifield Dataset, CMD, a collection of hundreds of thousands of 2D maps and 3D grids containing many different properties of cosmic gas, dark matter, and stars from 2,000 distinct simulated universes at several cosmic times. The 2D maps and 3D grids represent cosmic regions that span $\sim$100 million light years and have been generated from thousands of state-of-the-art hydrodynamic and gravity-only N-body simulations from the CAMELS project. Designed to train machine learning models, CMD is the largest dataset of its kind containing more than 70 Terabytes of data. In this paper we describe CMD in detail and outline a few of its applications. We focus our attention on one such task, parameter inference, formulating the problems we face as a challenge to the community. We release all data and provide further technical details at https://camels-multifield-dataset.readthedocs.io.

preprint2021arXiv

The IllustrisTNG Simulations: Public Data Release

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

preprint2021arXiv

The Role of Machine Learning in the Next Decade of Cosmology

In recent years, machine learning (ML) methods have remarkably improved how cosmologists can interpret data. The next decade will bring new opportunities for data-driven cosmological discovery, but will also present new challenges for adopting ML methodologies and understanding the results. ML could transform our field, but this transformation will require the astronomy community to both foster and promote interdisciplinary research endeavors.

preprint2020arXiv

Detecting the Cosmic Web: Lyα Emission from Simulated Filaments at z=3

The standard cosmological model ($Λ$CDM) predicts the existence of the cosmic web: a distribution of matter into sheets and filaments connecting massive halos. However, observational evidence has been elusive due to the low surface brightness of the filaments. Recent deep MUSE/VLT data and upcoming observations offer a promising avenue for Ly$α$ detection, motivating the development of modern theoretical predictions. We use hydrodynamical cosmological simulations run with the AREPO code to investigate the potential detectability of large-scale filaments, excluding contributions from the halos embedded in them. We focus on filaments connecting massive ($M_{200c}\sim(1-3)\times10^{12} M_\odot$) halos at z=3, and compare different simulation resolutions, feedback levels, and mock-image pixel sizes. We find increasing simulation resolution does not substantially improve detectability notwithstanding the intrinsic enhancement of internal filament structure. By contrast, for a MUSE integration of 31 hours, including feedback increases the detectable area by a factor of $\simeq$5.5 on average compared with simulations without feedback, implying that even the non-bound components of the filaments have substantial sensitivity to feedback. Degrading the image resolution from the native MUSE scale of (0.2")$^2$ per pixel to (5.3")$^2$ apertures has the strongest effect, increasing the detectable area by a median factor of $\simeq$200 and is most effective when the size of the pixel roughly matches the width of the filament. Finally, we find the majority of Ly$α$ emission is due to electron impact collisional excitations, as opposed to radiative recombination.

preprint2020arXiv

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

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

preprint2020arXiv

MOSEL and IllustrisTNG: Massive Extended Galaxies at z=2 Quench Later Than Normal-size Galaxies

Using the TNG100 (100 Mpc)^3 simulation of the IllustrisTNG project, we demonstrate a strong connection between the onset of star formation quenching and the stellar size of galaxies. We do so by tracking the evolutionary history of extended and normal-size galaxies selected at z=2 with log(M_star) = 10.2 - 11 and stellar-half-mass-radii above and within 1-sigma of the stellar size--stellar mass relation, respectively. We match the stellar mass and star formation rate distributions of the two populations. By z=1, only 36% of the extended massive galaxies have quenched, in contrast to a quenched fraction of 69% for the normal-size massive galaxies. We find that normal-size massive galaxies build up their central stellar mass without a significant increase in their stellar size between z=2-4, whereas the stellar size of the extended massive galaxies almost doubles in the same time. In IllustrisTNG, lower black hole masses and weaker kinetic-mode feedback appears to be responsible for the delayed quenching of star formation in the extended massive galaxies. We show that relatively gas-poor mergers may be responsible for the lower central stellar density and weaker supermassive black hole feedback in the extended massive galaxies.

preprint2020arXiv

The Angular Momentum of the Circumgalactic Medium in the TNG100 Simulation

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

preprint2020arXiv

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

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

preprint2019arXiv

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

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

preprint2019arXiv

Decoupling the rotation of stars and gas -- II: the link between black hole activity and MaNGA kinematics in TNG

We study the relationship between supermassive black hole (BH) feedback, BH luminosity and the kinematics of stars and gas for galaxies in IllustrisTNG. We use a sample of galaxies with mock MaNGA observations to identify kinematic misalignment at $z=0$ (difference in rotation of stars and gas), for which we follow the evolutionary history of BH activity and gas properties over the last 8 Gyrs. Misaligned low mass galaxies typically have boosted BH luminosity, BH growth and have had more energy injected into the gas over the last 8 Gyr in comparison to aligned galaxies. These properties likely lead to outflows and gas loss, in agreement with active low mass galaxies in observations. We show that splitting on BH luminosity at $z=0$ produces statistically consistent distributions of kinematic misalignment at $z=0$, however, splitting on the maximum BH luminosity over the last 8 Gyrs does not. While instantaneous correlation at $z=0$ is difficult due to misalignment persisting on longer timescales, the relationship between BH activity and misalignment is clear. High mass quenched galaxies with misalignment typically have similar BH luminosities, show no overall gas loss, and have typically lower gas phase metallicity over the last 8 Gyrs in comparison to those aligned; suggesting external origin.

preprint2019arXiv

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

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

preprint2016arXiv

How Environment Affects Galaxy Metallicity through Stripping and Formation History: Lessons from the Illustris Simulation

Recent studies have found higher galaxy metallicities in richer environments. It is not yet clear, however, whether metallicity-environment dependencies are merely an indirect consequence of environmentally dependent formation histories, or of environment related processes directly affecting metallicity. Here, we present a first detailed study of metallicity-environment correlations in a cosmological hydrodynamical simulation, in particular the Illustris simulation. Illustris galaxies display similar relations to those observed. Utilizing our knowledge of simulated formation histories, and leveraging the large simulation volume, we construct galaxy samples of satellites and centrals that are matched in formation histories. This allows us to find that ~1/3 of the metallicity-environment correlation is due to different formation histories in different environments. This is a combined effect of satellites (in particular, in denser environments) having on average lower z=0 star formation rates (SFRs), and of their older stellar ages, even at a given z=0 SFR. Most of the difference, ~2/3, however, is caused by the higher concentration of star-forming disks of satellite galaxies, as this biases their SFR-weighted metallicities toward their inner, more metal-rich parts. With a newly defined quantity, the 'radially averaged' metallicity, which captures the metallicity profile but is independent of the SFR profile, the metallicities of satellites and centrals become environmentally independent once they are matched in formation history. We find that circumgalactic metallicity (defined as rapidly inflowing gas around the virial radius), while sensitive to environment, has no measurable effect on the metallicity of the star-forming gas inside the galaxies.

preprint2016arXiv

Large-scale mass distribution in the Illustris simulation

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.

preprint2016arXiv

Shock finding on a moving-mesh: II. Hydrodynamic shocks in the Illustris universe

Hydrodynamical shocks are a manifestation of the non-linearity of the Euler equations and play a fundamental role in cosmological gas dynamics. In this work, we identify and analyse shocks in the Illustris simulation, and contrast the results with those of non-radiative runs. We show that simulations with more comprehensive physical models of galaxy formation pose new challenges for shock finding algorithms due to radiative cooling and star-forming processes, prompting us to develop a number of methodology improvements. We find in Illustris a total shock surface area which is about 1.4 times larger at the present epoch compared to non-radiative runs, and an energy dissipation rate at shocks which is higher by a factor of around 7. Remarkably, shocks with Mach numbers above and below $\mathcal{M}\approx10$ contribute about equally to the total dissipation across cosmic time. This is in sharp contrast to non-radiative simulations, and we demonstrate that a large part of the difference arises due to strong black hole radio-mode feedback in Illustris. We also provide an overview of the large diversity of shock morphologies, which includes complex networks of halo-internal shocks, shocks on to cosmic sheets, feedback shocks due to black holes and galactic winds, as well as ubiquitous accretion shocks. In high redshift systems more massive than $10^{12}\,\mathrm{M}_\odot$ we discover the existence of a double accretion shock pattern in haloes. They are created when gas streams along filaments without being shocked at the outer accretion shock, but then forms a second, roughly spherical accretion shock further inside.

preprint2015arXiv

Galactic Angular Momentum in the Illustris Simulation: Feedback and the Hubble Sequence

We study the stellar angular momentum of thousands of galaxies in the Illustris cosmological simulation, which captures gravitational and gas dynamics within galaxies, as well as feedback from stars and black holes. We find that the angular momentum of the simulated galaxies matches observations well, and in particular two distinct relations are found for late-type versus early-type galaxies. The relation for late-type galaxies corresponds to the value expected from full conservation of the specific angular momentum generated by cosmological tidal torques. The relation for early-type galaxies corresponds to retention of only ~30% of that, but we find that those early-type galaxies with low angular momentum at z=0 nevertheless reside at high redshift on the late-type relation. Some of them abruptly lose angular momentum during major mergers. To gain further insight, we explore the scaling relations in simulations where the galaxy formation physics is modified with respect to the fiducial model. We find that galactic winds with high mass-loading factors are essential for obtaining the high angular momentum relation typical for late-type galaxies, while AGN feedback largely operates in the opposite direction. Hence, feedback controls the stellar angular momentum of galaxies, and appears to be instrumental for establishing the Hubble sequence.

preprint2015arXiv

Galaxy Morphology and Star Formation in the Illustris Simulation at z=0

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.

preprint2015arXiv

Hot Gaseous Coronae around Spiral Galaxies: Probing the Illustris Simulation

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σ$ 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σ$ 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.

preprint2015arXiv

Hydrogen Reionization in the Illustris universe

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 $τ\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.

preprint2015arXiv

Modeling Galactic Conformity with the Color-Halo Age Relation in the Illustris Simulation

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.

preprint2015arXiv

Reducing noise in moving-grid codes with strongly-centroidal Lloyd mesh regularization

A method for improving the accuracy of hydrodynamical codes that use a moving Voronoi mesh is described. Our scheme is based on a new regularization scheme that constrains the mesh to be centroidal to high precision while still allowing the cells to move approximately with the local fluid velocity, thereby retaining the quasi-Lagrangian nature of the approach. Our regularization technique significantly reduces mesh noise that is attributed to changes in mesh topology and deviations from mesh regularity. We demonstrate the advantages of our method on various test problems, and note in particular improvements obtained in handling shear instabilities, mixing, and in angular momentum conservation. Calculations of adiabatic jets in which shear excites Kelvin Helmholtz instability show reduction of mesh noise and entropy generation. In contrast, simulations of the collapse and formation of an isolated disc galaxy are nearly unaffected, showing that numerical errors due to the choice of regularization do not impact the outcome in this case.

preprint2015arXiv

Reproducing the Kinematics of Damped Lyman-alpha Systems

We examine the kinematic structure of Damped Lyman-alpha Systems (DLAs) in a series of cosmological hydrodynamic simulations using the AREPO code. We are able to match the distribution of velocity widths of associated low ionisation metal absorbers substantially better than earlier work. Our simulations produce a population of DLAs dominated by halos with virial velocities around 70 km/s, consistent with a picture of relatively small, faint objects. In addition, we reproduce the observed correlation between velocity width and metallicity and the equivalent width distribution of SiII. Some discrepancies of moderate statistical significance remain; too many of our spectra show absorption concentrated at the edge of the profile and there are slight differences in the exact shape of the velocity width distribution. We show that the improvement over previous work is mostly due to our strong feedback from star formation and our detailed modelling of the metal ionisation state.

preprint2015arXiv

The diverse evolutionary paths of simulated high-z massive, compact galaxies to z=0

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.

preprint2015arXiv

The Formation of Massive, Compact Galaxies at z=2 in the Illustris Simulation

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.

preprint2015arXiv

The Illustris Simulation: Public Data Release

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.

preprint2015arXiv

The Illustris simulation: the evolving population of black holes across cosmic time

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.

preprint2015arXiv

The impact of feedback on cosmological gas accretion

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.

preprint2015arXiv

The Impact of Galactic Feedback on the Circumgalactic Medium

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-$α$ 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.

preprint2015arXiv

The merger rate of galaxies in the Illustris Simulation: a comparison with observations and semi-empirical models

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.

preprint2015arXiv

Zooming in on accretion - I. The structure of halo gas

We study the properties of gas in and around 10^12 solar mass halos at z=2 using a suite of high-resolution cosmological hydrodynamic 'zoom' simulations. We quantify the thermal and dynamical structure of these gaseous reservoirs in terms of their mean radial distributions and angular variability along different sightlines. With each halo simulated at three levels of increasing resolution, the highest reaching a baryon mass resolution of ~10,000 solar masses, we study the interaction of filamentary inflow and the quasi-static hot halo atmosphere. We highlight the discrepancy between the spatial resolution available in the halo gas as opposed to within the galaxy itself, and find that stream morphologies become increasingly complex at higher resolution, with large coherent flows revealing density and temperature structure at progressively smaller scales. Moreover, multiple gas components co-exist at the same radius within the halo, making radially averaged analyses misleading. This is particularly true where the hot, quasi-static, high entropy halo atmosphere interacts with cold, rapidly inflowing, low entropy accretion. We investigate the process of gas virialization and identify different regimes for the heating of gas as it accretes from the intergalactic medium. Haloes at this mass have a well-defined virial shock, associated with a sharp jump in temperature and entropy at ~1.25 r_vir. The presence, radius, and radial width of this boundary feature, however, vary not only from halo to halo, but also as a function of angular direction, covering roughly ~85% of the 4pi sphere. Our findings are relevant for the proper interpretation of observations pertaining to the circumgalactic medium, including evidence for large amounts of cold gas surrounding massive haloes at intermediate redshifts.

preprint2014arXiv

Damped Lyman-alpha absorbers as a probe of stellar feedback

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.

preprint2014arXiv

Halo Mass and Assembly History Exposed in the Faint Outskirts: the Stellar and Dark Matter Haloes of Illustris Galaxies

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.

preprint2014arXiv

Introducing the Illustris Project: Simulating the coevolution of dark and visible matter in the Universe

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.

preprint2014arXiv

Introducing the Illustris Project: the evolution of galaxy populations across cosmic time

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.

preprint2014arXiv

Properties of galaxies reproduced by a hydrodynamic simulation

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.

preprint2014arXiv

Synthetic Galaxy Images and Spectra from the Illustris Simulation

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.

preprint2014arXiv

The colors of satellite galaxies in the Illustris Simulation

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.

preprint2014arXiv

The star formation main sequence and stellar mass assembly of galaxies in the Illustris simulation

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.

preprint2013arXiv

A model for cosmological simulations of galaxy formation physics

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.

preprint2013arXiv

A physical model for cosmological simulations of galaxy formation: multi-epoch validation

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.

preprint2013arXiv

Following the flow: tracer particles in astrophysical fluid simulations

We present two numerical schemes for passive tracer particles in the hydrodynamical moving-mesh code AREPO, and compare their performance for various problems, from simple setups to cosmological simulations. The purpose of tracer particles is to allow the flow to be followed in a Lagrangian way, tracing the evolution of the fluid with time, and allowing the thermodynamical history of individual fluid parcels to be recorded. We find that the commonly-used `velocity field tracers', which are advected using the fluid velocity field, do not in general follow the mass flow correctly, and explain why this is the case. This method can result in orders-of-magnitude biases in simulations of driven turbulence and in cosmological simulations, rendering the velocity field tracers inappropriate for following these flows. We then discuss a novel implementation of `Monte Carlo tracers', which are moved along with fluid cells, and are exchanged probabilistically between them following the mass flux. This method reproduces the mass distribution of the fluid correctly. The main limitation of this approach is that it is more diffusive than the fluid itself. Nonetheless, we show that this novel approach is more reliable than what has been employed previously and demonstrate that it is appropriate for following hydrodynamical flows in mesh-based codes. The Monte Carlo tracers can also naturally be transferred between fluid cells and other types of particles, such as stellar particles, so that the mass flow in cosmological simulations can be followed in its entirety.

preprint2013arXiv

Moving mesh cosmology: tracing cosmological gas accretion

We investigate the nature of gas accretion onto haloes and galaxies at z=2 using cosmological hydrodynamic simulations run with the moving mesh code AREPO. Implementing a Monte Carlo tracer particle scheme to determine the origin and thermodynamic history of accreting gas, we make quantitative comparisons to an otherwise identical simulation run with the smoothed particle hydrodynamics (SPH) code GADGET-3. Contrasting these two numerical approaches, we find significant physical differences in the thermodynamic history of accreted gas in haloes above 10^10.5 solar masses. In agreement with previous work, GADGET simulations show a cold fraction near unity for galaxies forming in massive haloes, implying that only a small percentage of accreted gas heats to an appreciable fraction of the virial temperature during accretion. The same galaxies in AREPO show a much lower cold fraction, <20% in haloes above 10^11 solar masses. This results from a hot gas accretion rate which, at this same halo mass, is an order of magnitude larger than with GADGET, while the cold accretion rate is also lower. These discrepancies increase for more massive systems, and we explain both as due to numerical inaccuracies in the standard formulation of SPH. We also observe that the relatively sharp transition from cold to hot mode dominated accretion, at a halo mass of ~10^11, is a consequence of comparing past gas temperatures to a constant threshold value independent of virial temperature. Examining the spatial distribution of accreting gas, we find that gas filaments in GADGET tend to remain collimated and flow coherently to small radii, or artificially fragment and form a large number of purely numerical "blobs". Similar gas streams in AREPO show increased heating and disruption at 0.25-0.5 virial radii and contribute to the hot gas accretion rate in a manner distinct from classical cooling flows.

preprint2013arXiv

The Contribution of Halos with Different Mass Ratios to the Overall Growth of Cluster-Sized Halos

We provide a new observational test for a key prediction of the ΛCDM cosmological model: the contributions of mergers with different halo-to-main-cluster mass ratios to cluster-sized halo growth. We perform this test by dynamically analyzing seven galaxy clusters, spanning the redshift range $0.13 < z_c < 0.45$ and caustic mass range $0.4-1.5$ $10^{15} h_{0.73}^{-1}$ M$_{\odot}$, with an average of 293 spectroscopically-confirmed bound galaxies to each cluster. The large radial coverage (a few virial radii), which covers the whole infall region, with a high number of spectroscopically identified galaxies enables this new study. For each cluster, we identify bound galaxies. Out of these galaxies, we identify infalling and accreted halos and estimate their masses and their dynamical states. Using the estimated masses, we derive the contribution of different mass ratios to cluster-sized halo growth. For mass ratios between ~0.2 and ~0.7, we find a ~1 $σ$ agreement with ΛCDM expectations based on the Millennium simulations I and II. At low mass ratios, $\lesssim 0.2$, our derived contribution is underestimated since the detection efficiency decreases at low masses, $\sim 2 \times 10^{14}$ $h_{0.73}^{-1}$ M$_{\odot}$. At large mass ratios, $\gtrsim 0.7$, we do not detect halos probably because our sample, which was chosen to be quite X-ray relaxed, is biased against large mass ratios. Therefore, at large mass ratios, the derived contribution is also underestimated.

preprint2013arXiv

The effect of metal enrichment and galactic winds on galaxy formation in cosmological zoom simulations

We investigate the differential effects of metal cooling and galactic stellar winds on the cosmological formation of individual galaxies with three sets of cosmological, hydrodynamical zoom simulations of 45 halos in the mass range 10^11<M_halo<10^13M_sun. Models including both galactic winds and metal cooling (i) suppress early star formation at z>1 and predict reasonable star formation histories, (ii) produce galaxies with high cold gas fractions (30-60 per cent) at high redshift, (iii) significantly reduce the galaxy formation efficiencies for halos (M_halo<10^12M_sun) at all redshifts in agreement with observational and abundance matching constraints, (iv) result in high-redshift galaxies with reduced circular velocities matching the observed Tully-Fisher relation at z~2, and (v) significantly increase the sizes of low-mass galaxies (M_stellar<3x10^10M_sun) at high redshift resulting in a weak size evolution - a trend in agreement with observations. However, the low redshift (z<0.5) star formation rates of massive galaxies are higher than observed (up to ten times). No tested model predicts the observed size evolution for low-mass and high-mass galaxies simultaneously. Due to the delayed onset of star formation in the wind models, the metal enrichment of gas and stars is delayed and agrees well with observational constraints. Metal cooling and stellar winds are both found to increase the ratio of in situ formed to accreted stars - the relative importance of dissipative vs. dissipationless assembly. For halo masses below ~10^12M_sun, this is mainly caused by less stellar accretion and compares well to predictions from semi-analytical models but still differs from abundance matching models. For higher masses, the fraction of in situ stars is over-predicted due to the unrealistically high star formation rates at low redshifts.

preprint2013arXiv

The SINS/zC-SINF Survey of z~2 Galaxy Kinematics: The Nature of Dispersion Dominated Galaxies

We analyze the spectra, spatial distributions and kinematics of Ha, [NII] and [SII] emission in a sample of 42, z~2.2 UV/optically selected star forming galaxies (SFGs) from the SINS & zC-SINF surveys, 35 of which were observed in the adaptive optics mode of SINFONI. This is supplemented by kinematic data from 48 z~1-2.5 galaxies from the literature. We find that the kinematic classification of the high-z SFGs as `dispersion dominated' or `rotation dominated' correlates most strongly with their intrinsic sizes. Smaller galaxies are more likely `dispersion-dominated' for two main reasons: 1) The rotation velocity scales linearly with galaxy size but intrinsic velocity dispersion does not depend on size, and as such, their ratio is systematically lower for smaller galaxies, and 2) Beam smearing strongly decreases large-scale velocity gradients and increases observed dispersion much more for galaxies with sizes at or below the resolution. Dispersion dominated SFGs may thus have intrinsic properties similar to `rotation dominated' SFGs, but are primarily more compact, lower mass, less metal enriched and may have higher gas fractions, plausibly because they represent an earlier evolutionary state.

preprint2012arXiv

On the effect of cosmological inflow on turbulence and instability in galactic discs

We analyse the evolution of turbulence and gravitational instability of a galactic disc in a quasi-steady state governed by cosmological inflow. We focus on the possibility that the coupling between the in-streaming gas and the disc is maximal, e.g., via dense clumps, and ask whether the streams could be the driver of turbulence in an unstable disc with a Toomre parameter Q~1. Our fiducial model assumes an efficiency of ~0.5 per dynamical time for the decay of turbulence energy, and ~0.02 for each of the processes that deplete the disc gas, i.e., star formation, outflow, and inflow within the disc into a central bulge. In this case, the in-streaming drives a ratio of turbulent to rotation velocity sigma/V~0.2-0.3, which at z~2 induces an instability with Q~1, both as observed. However, in conflict with observations, this model predicts that sigma/V remains constant with time, independent of the cosmological accretion rate, because mass and turbulence have the same external source. Such strongly coupled cosmological inflow tends to stabilize the disc at low z, with Q ~ a few, which may be consistent with observations. The instability could instead be maintained for longer, with a properly declining sigma/V, if it is self-regulated to oscillations about Q~1 by a duty cycle for disc depletion. However, the 'off' phases of this duty cycle become long at low z, which may be hard to reconcile with observations. Alternatively, the coupling between the in-streaming gas and the disc may weaken in time, reflecting an evolving nature of the accretion. If, instead, that coupling is weak at all times, the likely energy source for self-regulated stirring up of the turbulence is the inflow within the disc down the potential gradient (studied in a companion paper).

preprint2012arXiv

Shocked Superwinds from the z~2 Clumpy Star-forming Galaxy, ZC406690

We have obtained high-resolution data of the z 2 ring-like, clumpy star-forming galaxy (SFG) ZC406690 using the VLT/SINFONI with AO (in K-band) and in seeing-limited mode (in H- and J-band). Our data includes all of the main strong optical emission lines: [OII], [OIII], Ha, Hb, [NII] and [SII]. We find broad, blueshifted Ha and [OIII] emission line wings in the spectra of the galaxy's massive, star-forming clumps (sigma \sim 85 km s^-1) and even broader wings (up to 70% of the total Ha flux, with sigma \sim 290 km s^-1) in regions spatially offset from the clumps by \sim 2 kpc. The broad emission likely originates from large-scale outflows with mass outflow rates from individual clumps that are 1-8x the SFR of the clumps. Based on emission line ratio diagnostics ([NII]/Ha and [SII]/Ha) and photoionization and shock models, we find that the emission from the clumps is due to a combination of photoionization from the star-forming regions and shocks generated in the outflowing component, with 5-30% of the emission deriving from shocks. In terms of the ionization parameter (6x10^7-10^8 cm/s, based on both the SFR and the O32 ratio), density (local electron densities of 300-1800 cm^-3 in and around the clumps, and ionized gas column densities of 1200-8000 Msol/pc^2), and SFR (10-40 Msol/yr), these clumps more closely resemble nuclear starburst regions of local ULIRGs and dwarf irregulars than HII regions in local galaxies. However, the star-forming clumps are not located in the nucleus as in local starburst galaxies but instead are situated in a ring several kpc from the center of their high-redshift host galaxy, and have an overall disk-like morphology. The two brightest clumps are quite different in terms of their internal properties, energetics and relative ages, and thus we are given a glimpse at two different stages in the formation and evolution of rapidly star-forming giant clumps at high-z.

preprint2012arXiv

The SINS/zC-SINF survey of z~2 galaxy kinematics: Outflow properties

Based on SINFONI Ha, [NII] and [SII] AO data of 30 z \sim 2 star-forming galaxies (SFGs) from the SINS and zcSINF surveys, we find a strong correlation of the Ha broad flux fraction with the star formation surface density of the galaxy, with an apparent threshold for strong outflows occurring at 1 Msun yr^-1 kpc^-2. Above this threshold, we find that SFGs with logm_\ast>10 have similar or perhaps greater wind mass loading factors (eta = Mdotout/SFR) and faster outflow velocities than lower mass SFGs. This trend suggests that the majority of outflowing gas at z \sim 2 may derive from high-mass SFGs, and that the z \sim 2 mass-metallicity relation is driven more by dilution of enriched gas in the galaxy gas reservoir than by the efficiency of outflows. The mass loading factor is also correlated with the SFR and inclination, such that more star-forming and face-on galaxies launch more powerful outflows. For galaxies that have evidence for strong outflows, we find that the broad emission is spatially extended to at least the half-light radius (\sim a few kpc). We propose that the observed threshold for strong outflows and the observed mass loading of these winds can be explained by a simple model wherein break-out of winds is governed by pressure balance in the disk. Using the ratio of the [SII] doublet in a broad and narrow component, we find that outflowing gas has a density of \sim10-100 cm^-3, significantly less than that of the star forming gas (600 cm^-3).

preprint2011arXiv

Evolution of violent gravitational disc instability in galaxies: Late stabilization by transition from gas to stellar dominance

We address the cosmological evolution of violent gravitational instability in high-redshift, massive, star-forming galactic discs. To this aim, we integrate in time the equations of mass and energy conservation under self-regulated instability of a two-component disc of gas and stars. The disc is assumed to be continuously fed by cold gas at the average cosmological rate. The gas forms stars and is partly driven away by stellar feedback. The gas and stars flow inward through the disc to a central bulge due to torques that drive angular momentum outwards. The gravitational energy released by the mass inflow down the gravitational potential gradient drives the disc turbulence that maintains the disc unstable with a Toomre instability parameter Q~1, compensating for the dissipative losses of the gas turbulence and raising the stellar velocity dispersion. We follow the velocity dispersion of stars and gas as they heat and cool respectively and search for disc `stabilization', to be marked by a low gas velocity dispersion comparable to the speed of sound ~10 km/s. We vary the model parameters that characterize the accreted gas fraction, turbulence dissipation rate, star-formation rate, and stellar feedback. We find that as long as the gas input roughly follows the average cosmological rate, the disc instability is a robust phenomenon at high redshift till z~1, driven by the high surface density and high gas fraction due to the intense cosmological accretion. For a broad range of model parameter values, the discs tend to `stabilize' at z~0-0.5 as they become dominated by hot stars. When the model parameters are pushed to extreme values, the discs may stabilize as early as z~2, with the gas loss by strong outflows serving as the dominant stabilizing factor.

preprint2011arXiv

Short-lived star-forming giant clumps in cosmological simulations of z~2 disks

Many observed massive star-forming z\approx2 galaxies are large disks that exhibit irregular morphologies, with \sim1kpc, \sim10^(8-10)Msun clumps. We present the largest sample to date of high-resolution cosmological SPH simulations that zoom-in on the formation of individual M*\sim10^(10.5)Msun galaxies in \sim10^(12)Msun halos at z\approx2. Our code includes strong stellar feedback parameterized as momentum-driven galactic winds. This model reproduces many characteristic features of this observed class of galaxies, such as their clumpy morphologies, smooth and monotonic velocity gradients, high gas fractions (f_g\sim50%) and high specific star-formation rates (\gtrsim1Gyr^(-1)). In accord with recent models, giant clumps (Mclump\sim(5x10^8-10^9)Msun) form in-situ via gravitational instabilities. However, the galactic winds are critical for their subsequent evolution. The giant clumps we obtain are short-lived and are disrupted by wind-driven mass loss. They do not virialise or migrate to the galaxy centers as suggested in recent work neglecting strong winds. By phenomenologically implementing the winds that are observed from high-redshift galaxies and in particular from individual clumps, our simulations reproduce well new observational constraints on clump kinematics and clump ages. In particular, the observation that older clumps appear closer to their galaxy centers is reproduced in our simulations, as a result of inside-out formation of the disks rather than inward clump migration.

preprint2010arXiv

Mergers in Lambda-CDM: Uncertainties in Theoretical Predictions and Interpretations of the Merger Rate

Different methodologies lead to order-of-magnitude variations in predicted galaxy merger rates. We examine and quantify the dominant uncertainties. Different halo merger rates and subhalo 'destruction' rates agree to within a factor ~2 given proper care in definitions. If however (sub)halo masses are not appropriately defined or are under-resolved, the major merger rate can be dramatically suppressed. The dominant differences in galaxy merger rates owe to baryonic physics. Hydrodynamic simulations without feedback and older models that do not agree with the observed galaxy mass function propagate factor ~5 bias in the resulting merger rates. However, if the model matches the galaxy mass function, properties of central galaxies are sufficiently converged to give small differences in merger rates. But variations in baryonic physics of satellites also have dramatic effects. The known problem of satellite 'over-quenching' in most semi-analytic models (SAMs), whereby SAM satellites are too efficiently stripped of gas, could lead to order-of-magnitude under-estimates of merger rates for low-mass, gas-rich galaxies. Fixing the satellite properties to observations tends to predict higher merger rates, but with factor ~2 empirical uncertainties. Choice of mass ratio definition matters: at low masses, most true major mergers (in baryonic/dynamical galaxy mass) will appear to be minor mergers in their stellar or luminosity mass ratio. Observations and models using these criteria may underestimate major merger rates by factors ~5. Orbital parameters and gas fractions also introduce factor ~3 differences in amount of bulge formed by mergers, even for fixed mass ratio encounters.

preprint2010arXiv

The growth of dark matter halos: evidence for significant smooth accretion

We study the growth of dark matter halos in the concordance LCDM cosmology using several N-body simulations of large cosmological volumes. We build merger trees from the Millennium and Millennium-II simulations, covering a range 10^9-10^15 Msun in halo mass and 1-10^5 in merger mass ratio. Our algorithm takes special care of halo fragmentation and ensures that the mass contribution of each merger to halo growth is only counted once. This way the integrated merger rate converges and we can consistently determine the contribution of mergers of different mass ratios to halo growth. We find that all resolved mergers, up to mass ratios of 10^5, contribute only ~60% of the total halo mass growth, while major mergers are subdominant, e.g. mergers with mass ratios smaller than 3:1 (10:1) contribute only ~20% (~30%). This is verified with an analysis of two additional simulation boxes, where we follow all particles individually throughout cosmic time. Our results are also robust against using several halo definitions. Under the assumption that the power-law behaviour of the merger rate at large mass ratios can be extrapolated to arbitrarily large mass ratios, it is found that, independently of halo mass, ~40% of the mass in halos comes from genuinely smooth accretion of dark matter that was never bound in smaller halos. We discuss possible implications of our findings for galaxy formation. One implication, assuming as is standard that the pristine intergalactic medium is heated and photoionized by UV photons, is that all halos accrete >40% of their baryons in smooth "cold" T>~10^4K gas, rather than as warm, enriched or clumpy gas or as stars.

Shy Genel

What is connected

Connect this record

See the researcher in context

Building this map preview

59 published item(s)

Cosmological back-reaction of baryons on dark matter in the CAMELS simulations

Boundaries of chaos and determinism in the cosmos

Breaking baryon-cosmology degeneracy with the electron density power spectrum

Cosmology with one galaxy?

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

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

The MUSE Hubble Ultra Deep Field Survey XVI. The angular momentum of low-mass star-forming galaxies. A cautionary tale and insights from TNG50

Finding universal relations in subhalo properties with artificial intelligence

The CAMELS Multifield Dataset: Learning the Universe's Fundamental Parameters with Artificial Intelligence

The IllustrisTNG Simulations: Public Data Release

The Role of Machine Learning in the Next Decade of Cosmology

Detecting the Cosmic Web: Lyα Emission from Simulated Filaments at z=3

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

MOSEL and IllustrisTNG: Massive Extended Galaxies at z=2 Quench Later Than Normal-size Galaxies

The Angular Momentum of the Circumgalactic Medium in the TNG100 Simulation

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

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

Decoupling the rotation of stars and gas -- II: the link between black hole activity and MaNGA kinematics in TNG

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

How Environment Affects Galaxy Metallicity through Stripping and Formation History: Lessons from the Illustris Simulation

Large-scale mass distribution in the Illustris simulation

Shock finding on a moving-mesh: II. Hydrodynamic shocks in the Illustris universe

Galactic Angular Momentum in the Illustris Simulation: Feedback and the Hubble Sequence

Galaxy Morphology and Star Formation in the Illustris Simulation at z=0

Hot Gaseous Coronae around Spiral Galaxies: Probing the Illustris Simulation

Hydrogen Reionization in the Illustris universe

Modeling Galactic Conformity with the Color-Halo Age Relation in the Illustris Simulation

Reducing noise in moving-grid codes with strongly-centroidal Lloyd mesh regularization

Reproducing the Kinematics of Damped Lyman-alpha Systems

The diverse evolutionary paths of simulated high-z massive, compact galaxies to z=0

The Formation of Massive, Compact Galaxies at z=2 in the Illustris Simulation

The Illustris Simulation: Public Data Release

The Illustris simulation: the evolving population of black holes across cosmic time

The impact of feedback on cosmological gas accretion

The Impact of Galactic Feedback on the Circumgalactic Medium

The merger rate of galaxies in the Illustris Simulation: a comparison with observations and semi-empirical models

Zooming in on accretion - I. The structure of halo gas

Damped Lyman-alpha absorbers as a probe of stellar feedback

Halo Mass and Assembly History Exposed in the Faint Outskirts: the Stellar and Dark Matter Haloes of Illustris Galaxies

Introducing the Illustris Project: Simulating the coevolution of dark and visible matter in the Universe

Introducing the Illustris Project: the evolution of galaxy populations across cosmic time

Properties of galaxies reproduced by a hydrodynamic simulation

Synthetic Galaxy Images and Spectra from the Illustris Simulation

The colors of satellite galaxies in the Illustris Simulation

The star formation main sequence and stellar mass assembly of galaxies in the Illustris simulation

A model for cosmological simulations of galaxy formation physics

A physical model for cosmological simulations of galaxy formation: multi-epoch validation

Following the flow: tracer particles in astrophysical fluid simulations

Moving mesh cosmology: tracing cosmological gas accretion

The Contribution of Halos with Different Mass Ratios to the Overall Growth of Cluster-Sized Halos

The effect of metal enrichment and galactic winds on galaxy formation in cosmological zoom simulations

The SINS/zC-SINF Survey of z~2 Galaxy Kinematics: The Nature of Dispersion Dominated Galaxies

On the effect of cosmological inflow on turbulence and instability in galactic discs

Shocked Superwinds from the z~2 Clumpy Star-forming Galaxy, ZC406690

The SINS/zC-SINF survey of z~2 galaxy kinematics: Outflow properties

Evolution of violent gravitational disc instability in galaxies: Late stabilization by transition from gas to stellar dominance

Short-lived star-forming giant clumps in cosmological simulations of z~2 disks

Mergers in Lambda-CDM: Uncertainties in Theoretical Predictions and Interpretations of the Merger Rate

The growth of dark matter halos: evidence for significant smooth accretion