Source author record

Vicente Rodriguez-Gomez

Vicente Rodriguez-Gomez appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

32works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

32 published item(s)

preprint2022arXiv

A machine learning approach to infer the accreted stellar mass fractions of central galaxies in the TNG100 simulation

We propose a random forest (RF) machine learning approach to determine the accreted stellar mass fractions ($f_\mathrm{acc}$) of central galaxies, based on various dark matter halo and galaxy features. The RF is trained and tested using 2,710 galaxies with stellar mass $\log_{10}M_\ast/M_\odot>10.16$ from the TNG100 simulation. Galaxy size is the most important individual feature when calculated in 3-dimensions, which becomes less important after accounting for observational effects. For smaller galaxies, the rankings for features related to merger histories increase. When an entire set of halo and galaxy features are used, the prediction is almost unbiased, with root-mean-square error (RMSE) of $\sim$0.068. A combination of up to three features with different types (galaxy size, merger history and morphology) already saturates the power of prediction. If using observable features, the RMSE increases to $\sim$0.104, and a combined usage of stellar mass, galaxy size plus galaxy concentration achieves similar predictions. Lastly, when using galaxy density, velocity and velocity dispersion profiles as features, which approximately represent the maximum amount of information extracted from galaxy images and velocity maps, the prediction is not improved much. Hence the limiting precision of predicting $f_\mathrm{acc}$ is $\sim$0.1 with observables, and the multi-component decomposition of galaxy images should have similar or larger uncertainties. If the central black hole mass and the spin parameter of galaxies can be accurately measured in future observations, the RMSE is promising to be further decreased by $\sim$20%.

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

Identifying Galaxy Mergers in Simulated CEERS NIRCam Images using Random Forests

Identifying merging galaxies is an important - but difficult - step in galaxy evolution studies. We present random forest classifications of galaxy mergers from simulated JWST images based on various standard morphological parameters. We describe (a) constructing the simulated images from IllustrisTNG and the Santa Cruz SAM, and modifying them to mimic future CEERS observations as well as nearly noiseless observations, (b) measuring morphological parameters from these images, and (c) constructing and training the random forests using the merger history information for the simulated galaxies available from IllustrisTNG. The random forests correctly classify $\sim60\%$ of non-merging and merging galaxies across $0.5 < z < 4.0$. Rest-frame asymmetry parameters appear more important for lower redshift merger classifications, while rest-frame bulge and clump parameters appear more important for higher redshift classifications. Adjusting the classification probability threshold does not improve the performance of the forests. Finally, the shape and slope of the resulting merger fraction and merger rate derived from the random forest classifications match with theoretical Illustris predictions, but are underestimated by a factor of $\sim 0.5$.

preprint2022arXiv

Realistic galaxy images and improved robustness in machine learning tasks from generative modelling

We examine the capability of generative models to produce realistic galaxy images. We show that mixing generated data with the original data improves the robustness in downstream machine learning tasks. We focus on three different data sets; analytical Sérsic profiles, real galaxies from the COSMOS survey, and galaxy images produced with the SKIRT code, from the IllustrisTNG simulation. We quantify the performance of each generative model using the Wasserstein distance between the distributions of morphological properties (e.g. the Gini-coefficient, the asymmetry, and ellipticity), the surface brightness distribution on various scales (as encoded by the power-spectrum), the bulge statistic and the colour for the generated and source data sets. With an average Wasserstein distance (Fréchet Inception Distance) of $7.19 \times 10^{-2}\, (0.55)$, $5.98 \times 10^{-2}\, (1.45)$ and $5.08 \times 10^{-2}\, (7.76)$ for the Sérsic, COSMOS and SKIRT data set, respectively, our best models convincingly reproduce even the most complicated galaxy properties and create images that are visually indistinguishable from the source data. We demonstrate that by supplementing the training data set with generated data, it is possible to significantly improve the robustness against domain-shifts and out-of-distribution data. In particular, we train a convolutional neural network to denoise a data set of mock observations. By mixing generated images into the original training data, we obtain an improvement of $11$ and $45$ per cent in the model performance regarding domain-shifts in the physical pixel size and background noise level, respectively.

preprint2022arXiv

The formation of low surface brightness galaxies in the IllustrisTNG simulation

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

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.

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

MOSEL Survey: Tracking the Growth of Massive Galaxies at 2<z<4 using Kinematics and the IllustrisTNG Simulation

We use K-band spectroscopic data from the Multi-Object Spectroscopic Emission Line (MOSEL) survey to analyze the kinematic properties of galaxies at z>3. Our sample consists of 34 galaxies at 3.0<zspec<3.8 between 9.0<log(M_star)<11.0. We find that galaxies with log(M_star) > 10.2 at z > 3 have 56 +/- 21 km/s lower integrated velocity dispersion compared to galaxies at z ~ 2 of similar stellar mass. Massive galaxies at z > 3 have either a flat or declining star formation history (SFH), whereas similar stellar mass galaxies at z~2.0 exhibit a slight peak in the past 500 Myrs. Comparing with the IllustrisTNG cosmological simulation, we find that (i) the dynamical mass of massive galaxies in simulations (log(M_star) > 10.0) increases by ~0.1 dex at a fixed stellar mass between z=2.0-3.0, and (ii) dynamical mass growth is coupled with a rapid rise in the ex-situ stellar mass fraction (stars accreted from other galaxies) for massive galaxies at z < 3.5. We speculate that the rising contribution of ex-situ stellar mass to the total stellar mass growth of massive galaxies is driving the higher integrated velocity dispersion and rising SFHs of massive galaxies at z~2.0 compared to galaxies of similar stellar masses at z > 3.

preprint2020arXiv

The bivariate gas-stellar mass distributions and the mass functions of early- and late-type galaxies at $z\sim0$

We report the bivariate HI- and H$_2$-stellar mass distributions of local galaxies in addition of an inventory of galaxy mass functions, MFs, for HI, H$_2$, cold gas, and baryonic mass, separately into early- and late-type galaxies. The MFs are determined using the HI and H$_2$ conditional distributions and the galaxy stellar mass function, GSMF. For the conditional distributions we use the compilation presented in Calette et al. 2018. For determining the GSMF from $M_{\ast}\sim3\times10^{7}$ to $3\times10^{12}$ $M_{\odot}$, we combine two spectroscopic samples from the SDSS at the redshift range $0.0033<z<0.2$. We find that the low-mass end slope of the GSMF, after correcting from surface brightness incompleteness, is $α\approx-1.4$, consistent with previous determinations. The obtained HI MFs agree with radio blind surveys. Similarly, the H$_2$ MFs are consistent with CO follow-up optically-selected samples. We estimate the impact of systematics due to mass-to-light ratios and find that our MFs are robust against systematic errors. We deconvolve our MFs from random errors to obtain the intrinsic MFs. Using the MFs, we calculate cosmic density parameters of all the baryonic components. Baryons locked inside galaxies represent 5.4% of the universal baryon content, while $\sim96$% of the HI and H$_2$ mass inside galaxies reside in late-type morphologies. Our results imply cosmic depletion times of H$_2$ and total neutral H in late-type galaxies of $\sim 1.3$ and 7.2 Gyr, respectively, which shows that late type galaxies are on average inefficient in converting H$_2$ into stars and in transforming HI gas into H$_2$. Our results provide a fully self-consistent empirical description of galaxy demographics in terms of the bivariate gas--stellar mass distribution and their projections, the MFs. This description is ideal to compare and/or to constrain galaxy formation models.

preprint2020arXiv

The fate of disk galaxies in IllustrisTNG clusters

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

preprint2020arXiv

The Morphology-Density relationship in 1<z<2 clusters

The morphology-density relationship states that dense cosmic environments such as galaxy clusters have an overabundance of quiescent elliptical galaxies, but it is unclear at which redshift this relationship is first established. We study the morphology of 4 clusters with $1.2<z<1.8$ using HST imaging and the morphology computation code statmorph. By comparing median morphology of cluster galaxies to CANDELS field galaxies using Monte Carlo analysis, we find that 2 out of 4 clusters (at z=1.19 and z=1.75) have an established morphology-density relationship with more than $3σ$ significance. $\sim$50% of galaxies in these clusters are bulge-dominated compared to $\sim$30% in the field, and they are significantly more compact. This result is more significant for low-mass galaxies with $\log M/M_\odot \lessapprox 10.5$, showing that low-mass galaxies are affected the most in clusters. We also find an intriguing system of two z $\approx$ 1.45 clusters at a unusually small separation 2D separation of $3'$ and 3D separation of $\approx73$ Mpc that exhibit no morphology-density relationship but have enhanced merger signatures. We conclude that the environmental mechanism responsible for the morphology-density relationship is 1) already active as early as z=1.75, 2) forms compact, bulge-dominated galaxies and 3) affects primarily low-mass galaxies. However, there is a significant degree of intracluster variance that may depend on the larger cosmological environment in which the cluster is embedded.

preprint2020arXiv

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

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

preprint2019arXiv

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

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

preprint2019arXiv

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

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

preprint2016arXiv

Sussing Merger Trees: Stability and Convergence

Merger trees are routinely used to follow the growth and merging history of dark matter haloes and subhaloes in simulations of cosmic structure formation. Srisawat et al. (2013) compared a wide range of merger-tree-building codes. Here we test the influence of output strategies and mass resolution on tree-building. We find that, somewhat surprisingly, building the tree from more snapshots does not generally produce more complete trees; instead, it tends to short- en them. Significant improvements are seen for patching schemes which attempt to bridge over occasional dropouts in the underlying halo catalogues or schemes which combine the halo-finding and tree-building steps seamlessly. The adopted output strategy does not affec- t the average number of branches (bushiness) of the resultant merger trees. However, mass resolution has an influence on both main branch length and the bushiness. As the resolution increases, a halo with the same mass can be traced back further in time and will encounter more small progenitors during its evolutionary history. Given these results, we recommend that, for simulations intended as precursors for galaxy formation models where of order 100 or more snapshots are analysed, the tree-building routine should be integrated with the halo finder, or at the very least be able to patch over multiple adjacent snapshots.

preprint2015arXiv

An analysis of the evolving comoving number density of galaxies in hydrodynamical simulations

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

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

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

On the Assembly of Dwarf Galaxies in Clusters and their Efficient Formation of Globular Clusters

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

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

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

Periodic Radio Continuum Emission Associated with the beta Cephei Star V2187 Cyg

We present new optical time-resolved photometry and medium-resolution spectroscopy of V2187 Cyg. We confirm its classification as a beta Cephei star based on sinusoidal light variations with a period of 0.2539 days and mean amplitudes of 0.037 and 0.042 magnitudes in "i" and "V", respectively. We classified the spectrum of this star B2-3V with no evidence of variations in the profiles of its absorption lines in timescales of hours or days. The stellar spectrum is totally absent of emission lines. We detected unexpected faint radio continuum emission (between 0.4 and 0.8 mJy at 6-cm) showing a sinusoidal variation with a period of 12.8 days. The radio spectrum is thermal. We searched in the Very Large Array archive for radio continuum emission toward other 15 beta Cephei stars. None of these additional stars, some of them much closer to the Sun than V2187 Cyg, was detected, indicating that radio emission is extremely uncommon toward beta Cephei stars.

preprint2014arXiv

Sussing Merger Trees : The Impact of Halo Merger Trees on Galaxy Properties in a Semi-Analytic Model

A halo merger tree forms the essential backbone of a semi-analytic model for galaxy formation and evolution. Recent studies have pointed out that extracting merger trees from numerical simulations of structure formation is non-trivial; different tree building algorithms can give differing merger histories. These differences should be carefully understood before merger trees are used as input for models of galaxy formation. We investigate the impact of different halo merger trees on a semi-analytic model. We find that the z=0 galaxy properties in our model show differences between trees when using a common parameter set. The star formation history of the Universe and the properties of satellite galaxies can show marked differences between trees with different construction methods. Independently calibrating the semi-analytic model for each tree can reduce the discrepancies between the z=0 global galaxy properties, at the cost of increasing the differences in the evolutionary histories of galaxies. Furthermore, the underlying physics implied can vary, resulting in key quantities such as the supernova feedback efficiency differing by factors of 2. Such a change alters the regimes where star formation is primarily suppressed by supernovae. Therefore, halo merger trees extracted from a common halo catalogue using different, but reliable, algorithms can result in a difference in the semi-analytic model. Given the uncertainties in galaxy formation physics, however, these differences may not necessarily be viewed as significant.

preprint2014arXiv

Sussing Merger Trees: the influence of the halo finder

Merger tree codes are routinely used to follow the growth and merger of dark matter haloes in simulations of cosmic structure formation. Whereas in Srisawat et. al. we compared the trees built using a wide variety of such codes here we study the influence of the underlying halo catalogue upon the resulting trees. We observe that the specifics of halo finding itself greatly influences the constructed merger trees. We find that the choices made to define the halo mass are of prime importance. For instance, amongst many potential options different finders select self-bound objects or spherical regions of defined overdensity, decide whether or not to include substructures within the mass returned and vary in their initial particle selection. The impact of these decisions is seen in tree length (the period of time a particularly halo can be traced back through the simulation), branching ratio (essentially the merger rate of \subhalos) and mass evolution. We therefore conclude that the choice of the underlying halo finder is more relevant to the process of building merger trees than the tree builder itself. We also report on some built-in features of specific merger tree codes that (sometimes) help to improve the quality of the merger trees produced.

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.

preprint2013arXiv

Sussing Merger Trees: The Merger Trees Comparison Project

Merger trees follow the growth and merger of dark-matter haloes over cosmic history. As well as giving important insights into the growth of cosmic structure in their own right, they provide an essential backbone to semi-analytic models of galaxy formation. This paper is the first in a series to arise from the SUSSING MERGER TREES Workshop in which ten different tree-building algorithms were applied to the same set of halo catalogues and their results compared. Although many of these codes were similar in nature, all algorithms produced distinct results. Our main conclusions are that a useful merger-tree code should possess the following features: (i) the use of particle IDs to match haloes between snapshots; (ii) the ability to skip at least one, and preferably more, snapshots in order to recover subhaloes that are temporarily lost during merging; (iii) the ability to cope with (and ideally smooth out) large, temporary flucuations in halo mass. Finally, to enable different groups to communicate effectively, we defined a common terminology that we used when discussing merger trees and we encourage others to adopt the same language. We also specified a minimal output format to record the results.