Search

Travel to academic conferences -- where international flights are the norm -- is responsible for a sizeable fraction of the greenhouse gas (GHG) emissions associated with academic work. In order to provide a benchmark for comparison with other fields, as well as for future reduction strategies and assessments, we estimate the CO2-equivalent emissions for conference travel in the field of astronomy for the prepandemic year 2019. The GHG emission of the international astronomical community's 362 conferences and schools in 2019 amounted to 42,500 tCO2e, assuming a radiative-forcing index factor of 1.95 for air travel. This equates to an average of 1.0 $\pm$ 0.6 tCO2e per participant per meeting. The total travel distance adds up to roughly 1.5 Astronomical Units, that is, 1.5 times the distance between the Earth and the Sun. We present scenarios for the reduction of this value, for instance with virtual conferencing or hub models, while still prioritizing the benefits conferences bring to the scientific community., Comment: Supplementary material is available at PNAS Nexus online: https://academic.oup.com/pnasnexus/article/3/5/pgae143/7659884

After more than five years of development, we present a new version of Dark Sage, a semi-analytic model (SAM) of galaxy formation that breaks the mould for models of its kind. Included among the major changes is an overhauled treatment of stellar feedback that is derived from energy conservation, operates on local scales, affects gas gradually over time rather than instantaneously, and predicts a mass-loading factor for every galaxy. Building on the model's resolved angular-momentum structure of galaxies, we now consider the heating of stellar discs, delivering predictions for disc structure both radially and vertically. We add a further dimension to stellar discs by tracking the distribution of stellar ages in each annulus. Each annulus--age bin has its own velocity dispersion and metallicity evolved in the model. This allows Dark Sage to make structural predictions for galaxies that previously only hydrodynamic simulations could. We present the model as run on the merger trees of the highest-resolution gravity-only simulation of the MillenniumTNG suite. Despite its additional complexity relative to other SAMs, Dark Sage only has three free parameters, the least of any SAM, which we calibrate exclusively against the cosmic star formation history and the $z=0$ stellar and HI mass functions using a particle-swarm optimisation method. The Dark Sage codebase, written in C and Python, is publicly available at https://github.com/arhstevens/DarkSage, Comment: Updates post referee's report. Get in touch for data access. 34 pages before references; 14 figures, 2 tables, 130 equations

Semi-analytic models (SAMs) systematically predict higher stellar-mass scatter at a given halo mass than hydrodynamical simulations and most empirical models. Our goal is to investigate the physical origin of this scatter by exploring modifications to the physics in the SAM Dark Sage. We design two black hole formation models that approximate results from the IllustrisTNG 300-1 hydrodynamical simulation. In the first model, we assign a fixed black hole mass of $10^{6}\, \mathrm{M}_{\odot}$ to every halo that reaches $10^{10.5}\, \mathrm{M}_{\odot}$. In the second model, we disregard any black hole growth as implemented in the standard Dark Sage model. Instead, we force all black hole masses to follow the median black hole mass-halo mass relation in IllustrisTNG 300-1 with a fixed scatter. We find that each model on its own does not significantly reduce the scatter in stellar mass. To do this, we replace the native Dark Sage AGN feedback model with a simple model where we turn off cooling for galaxies with black hole masses above $10^{8}\, \mathrm{M}_{\odot}$. With this additional modification, the SMBH seeding and fixed conditional distribution models find a significant reduction in the scatter in stellar mass at halo masses between $10^{11-14}\, \mathrm{M}_{\odot}$. These results suggest that AGN feedback in SAMs acts in a qualitatively different way than feedback implemented in cosmological simulations. Either or both may require substantial modification to match the empirically inferred scatter in the Stellar Mass Halo Mass Relation (SMHMR)., Comment: 21 pages, 16 figures

It has been shown in previous publications that the TNG100 simulation quantitatively reproduces the observed reduction in each of the total atomic and total molecular hydrogen gas for galaxies within massive halos, i.e.~dense environments. In this Letter, we study how well TNG50 reproduces the resolved effects of a Virgo-like cluster environment on the gas surface densities of satellite galaxies with $m_* > \! 10^9\,{\rm M}_\odot$ and ${\rm SFR} \! > 0.05\,{\rm M}_\odot\,{\rm yr}^{-1}$. We select galaxies in the simulation that are analogous to those in the HERACLES and VERTICO surveys, and mock-observe them to the common specifications of the data. Although TNG50 does not quantitatively match the observed gas surface densities in the centers of galaxies, the simulation does qualitatively reproduce the trends of gas truncation and central density suppression seen in VERTICO in both HI and H$_2$. This result promises that modern cosmological hydrodynamic simulations can be used to reliably model the post-infall histories of cluster satellite galaxies., Comment: 9 pages, 4 figures, accepted in ApJL

We study how environment regulates the star formation cycle of 33 Virgo Cluster satellite galaxies on 720 parsec scales. We present the first resolved star-forming main sequence for cluster galaxies, dividing the sample based on their global HI properties and comparing to a control sample of field galaxies. HI-poor cluster galaxies have reduced star formation rate (SFR) surface densities with respect to both HI-normal cluster and field galaxies (0.5 dex), suggesting that mechanisms regulating the global HI content are responsible for quenching local star formation. We demonstrate that the observed quenching in HI-poor galaxies is caused by environmental processes such as ram pressure stripping (RPS) simultaneously reducing molecular gas surface density and star formation efficiency (SFE), compared to regions in HI-normal systems (by 0.38 and 0.22 dex, respectively). We observe systematically elevated SFRs that are driven by increased molecular gas surface densities at fixed stellar mass surface density in the outskirts of early-stage RPS galaxies, while SFE remains unchanged with respect to the field sample. We quantify how RPS and starvation affect the star formation cycle of inner and outer galaxy discs as they are processed by the cluster. We show both are effective quenching mechanisms with the key difference being that RPS acts upon the galaxy outskirts while starvation regulates the star formation cycle throughout disc, including within the truncation radius. For both processes, the quenching is caused by a simultaneous reduction in molecular gas surface densities and SFE at fixed stellar mass surface density., Comment: 17 pages, 1 table, 5 figures, accepted for publication in ApJ

More precise measurements of galaxy clustering will be provided by the next generation of galaxy surveys such as DESI, WALLABY and SKA. To utilize this information to improve our understanding of the Universe, we need to accurately model the distribution of galaxies in their host dark matter halos. In this work we present a new galaxy number density profile of haloes, which makes predictions for the positions of galaxies in the host halo, different to the widely adopted Navarro-Frenk-White (NFW) profile, since galaxies tend to be found more in the outskirts of halos (nearer the virial radius) than an NFW profile. The parameterised galaxy number density profile model of haloes is fit and tested using the DARKSAGE semi-analytic model of galaxy formation. We find that our galaxy number density profile model of haloes can accurately reproduce the halo occupation distribution and galaxy two-point correlation function of the DARKSAGE simulation. We also derive the analytic expressions for the circular velocity and gravitational potential energy for this profile model. We use the SDSS DR10 galaxy group catalogue to validate this galaxy number density profile model of haloes. Compared to the NFW profile, we find that our model more accurately predicts the positions of galaxies in their host halo and the galaxy two-point correlation function., Comment: 13 pages. 10 figures. Appear on ApJ

We use the EAGLE (Evolution and Assembly of GaLaxies and their Environments) and IllustrisTNG (The Next Generation) cosmological simulations to investigate the properties of the baryonic specific angular momentum (j), baryonic mass (M) and atomic gas fraction ($f_{\rm{atm}}$) plane for nearby galaxies. We find EAGLE and TNG to be in excellent agreement with each other. These simulations are also consistent with the results obtained with xGASS (eXtended GALEX Arecibo SDSS Survey) for gas fractions greater than 0.01. This implies that the disagreements previously identified between xGASS and predictions from simple analytical disc stability arguments also holds true for EAGLE and TNG. For lower gas fraction (the regime currently unconstrained by observations), both simulations deviate from the plane but still maintain good agreement with each other. Despite the challenges posed by resolution limits at low gas fractions, our findings suggest a potential disconnect between angular momentum and gas fraction in the gas-poor regime, implying that not all gas-poor galaxies have low specific angular momentum., Comment: 12 pages, 10 figures. Accepted for publication in MNRAS

We analyze cold-gas distributions in Virgo cluster galaxies using resolved CO(2-1) (tracing molecular hydrogen, H2) and HI observations from the Virgo Environment Traced In CO (VERTICO) and the VLA Imaging of Virgo in Atomic Gas (VIVA) surveys. From a theoretical perspective, it is expected that environmental processes in clusters will have a stronger influence on diffuse atomic gas compared to the relatively dense molecular gas component, and that these environmental perturbations can compress the cold interstellar medium in cluster galaxies leading to elevated star formation. In this work we observationally test these predictions for star-forming satellite galaxies within the Virgo cluster. We divide our Virgo galaxy sample into HI-normal, HI-tailed, and HI-truncated classes and show, unsurprisingly, that the HI-tailed galaxies have the largest quantitative HI asymmetries. We also compare to a control sample of non-cluster galaxies and find that Virgo galaxies, on average, have HI asymmetries that are 40 +/- 10 per cent larger than the control. There is less separation between control, HI-normal, HI-tailed, and HI-truncated galaxies in terms of H2 asymmetries, and on average, Virgo galaxies have H2 asymmetries that are only marginally (20 +/- 10 per cent) larger than the control sample. We find a weak correlation between HI and H2 asymmetries over our entire sample, but a stronger correlation for those specific galaxies being strongly impacted by environmental perturbations. Finally, we divide the discs of the HI-tailed Virgo galaxies into a leading half and trailing half according to the observed tail direction. We find evidence for excess molecular gas mass on the leading halves of the disc. This excess molecular gas on the leading half is accompanied by an excess in star formation rate such that the depletion time is, on average, unchanged., Comment: 15 pages, 8 figures, 1 table, accepted for publication in A&A

We determine the atomic hydrogen (HI) to halo mass relation (HIHM) using Arecibo Legacy Fast ALFA survey HI data at the location of optically selected groups from the Galaxy and Mass Assembly (GAMA) survey. We make direct HI detections for 37 GAMA groups. Using HI group spectral stacking of 345 groups, we study the group HI content as function of halo mass across a halo mass range of $10^{11} - 10^{14.7}\text{ M}_\odot$. We also correct our results for Eddington bias. We find that the group HI mass generally rises as a function of halo mass from $1.3\%$ of the halo mass at $10^{11.6} \text{M}_\odot$ to $0.4\%$ at $10^{13.7} \text{M}_\odot$ with some indication of flattening towards the high-mass end. Despite the differences in optical survey limits, group catalogues, and halo mass estimation methods, our results are consistent with previous group HI-stacking studies. Our results are also consistent with mock observations from SHARK and IllustrisTNG., Comment: Accepted in MNRAS; 18 pages, 12 figures

The quenching of cluster satellite galaxies is inextricably linked to the suppression of their cold interstellar medium (ISM) by environmental mechanisms. While the removal of neutral atomic hydrogen (HI) at large radii is well studied, how the environment impacts the remaining gas in the centres of galaxies, which are dominated by molecular gas, is less clear. Using new observations from the Virgo Environment traced in CO survey (VERTICO) and archival HI data, we study the HI and molecular gas within the optical discs of Virgo cluster galaxies on 1.2-kpc scales with spatially resolved scaling relations between stellar (${\Sigma}_{\star}$), HI (${\Sigma}_\mathrm{HI}$), and molecular gas (${\Sigma}_\mathrm{mol}$) surface densities. Adopting HI deficiency as a measure of environmental impact, we find evidence that, in addition to removing the HI at large radii, the cluster processes also lower the average ${\Sigma}_\mathrm{HI}$ of the remaining gas even in the central 1.2 kpc. The impact on molecular gas is comparatively weaker than on the HI, and we show that the lower ${\Sigma}_\mathrm{mol}$ gas is removed first. In the most HI-deficient galaxies, however, we find evidence that environmental processes reduce the typical ${\Sigma}_\mathrm{mol}$ of the remaining gas by nearly a factor of 3. We find no evidence for environment-driven elevation of ${\Sigma}_\mathrm{HI}$ or ${\Sigma}_\mathrm{mol}$ in HI-deficient galaxies. Using the ratio of ${\Sigma}_\mathrm{mol}$-to-${\Sigma}_\mathrm{HI}$ in individual regions, we show that changes in the ISM physical conditions, estimated using the total gas surface density and midplane hydrostatic pressure, cannot explain the observed reduction in molecular gas content. Instead, we suggest that direct stripping of the molecular gas is required to explain our results., Comment: Accepted for publication in PASA. Main text is 19 pages including 12 figures and 3 tables, plus 1 appendix. A 2.5 min, high-level summary can be found at https://youtu.be/7djMmVEpDVc

While an astronomer's job is typically to look out from Earth, the seriousness of the climate crisis has meant a shift in many astronomers' focus. Astronomers are starting to consider how our resource requirements may contribute to this crisis and how we may better conduct our research in a more environmentally sustainable fashion. Astronomers for Planet Earth is an international organisation (more than 1,700 members from over 70 countries as of November 2022) that seeks to answer the call for sustainability to be at the heart of astronomers' practices. In this article, we review the organisation's history, summarising the proactive, collaborative efforts and research into astronomy sustainability conducted by its members. We update the state of affairs with respect to the carbon footprint of astronomy research, noting an improvement in renewable energy powering supercomputing facilities in Australia, reducing that component of our footprint by a factor of 2--3. We discuss how, despite accelerated changes made throughout the pandemic, we still must address the format of our meetings. Using recent annual meetings of the Australian and European astronomical societies as examples, we demonstrate that the more online-focussed a meeting is, the greater its attendance and the lower its emissions., Comment: 11 pages, 2 figures, refereed and accepted to CAPjournal for special issue on CAP2022 conference

We employ the EAGLE hydrodynamical simulation to uncover the relationship between cluster environment and $\rm H_2$ content of star-forming galaxies at redshifts spanning $0\leq z\leq 1$. To do so, we divide the star-forming sample into those that are bound to clusters and those that are not. We find that, at any given redshift, the galaxies in clusters generally have less $\rm H_2$ than their non-cluster counterparts with the same stellar mass (corresponding to an offset of $\lesssim 0.5$ dex), but this offset varies with stellar mass and is virtually absent at $M_\star\lesssim10^{9.3}~{\rm M}_\odot$. The $\rm H_2$ deficit in star-forming cluster galaxies can be traced back to a decline in their $\rm H_2$ content that commenced after first infall into a cluster, which occurred later than a typical cluster galaxy. Evolution of the full cluster population after infall is generally consistent with `slow-then-rapid' quenching, but galaxies with $M_\star\lesssim 10^{9.5}~{\rm M}_\odot$ exhibit rapid quenching. Unlike most cluster galaxies, star-forming ones were not pre-processed in groups prior to being accreted by clusters. For both of these cluster samples, the star formation efficiency remained oblivious to the infall. We track the particles associated with star-forming cluster galaxies and attribute the drop in $\rm H_2$ mass after infall to poor replenishment, depletion due to star formation, and stripping of $\rm H_2$ in cluster environments. These results provide predictions for future surveys, along with support and theoretical insights for existing molecular gas observations that suggest there is less $\rm H_2$ in cluster galaxies., Comment: Accepted for publication in MNRAS

We measure the molecular-to-atomic gas ratio, $R_{\rm mol}$, and the star formation rate (SFR) per unit molecular gas mass, SFE$_{\rm mol}$, in 38 nearby galaxies selected from the Virgo Environment Traced in CO (VERTICO) survey. We determine their scale-lengths for the molecular and stellar components and find a roughly 3:5 ratio between them compared to $\sim$1:1 in field galaxies, indicating that the CO emission is more centrally concentrated than the stars. We compute $R_{\rm mol}$ as a function of different physical quantities. While the spatially-resolved $R_{\rm mol}$ on average decreases with increasing radius, we find that the mean molecular-to-atomic gas ratio within the stellar effective radius $R_{\rm e}$, $R_{\rm mol}(r

The next generation of galaxy surveys will provide more precise measurements of galaxy clustering than have previously been possible. The 21-cm radio signals that are emitted from neutral atomic hydrogen (HI) gas will be detected by large-area radio surveys such as WALLABY and the SKA, and deliver galaxy positions and velocities that can be used to measure galaxy clustering statistics. But, to harness this information to improve our cosmological understanding, and learn about the physics of dark matter and dark energy, we need to accurately model the manner in which galaxies detected in HI trace the underlying matter distribution of the Universe. For this purpose, we develop a new HI-based Halo Occupation Distribution (HOD) model, which makes predictions for the number of galaxies present in dark matter halos conditional on their HI mass. The parameterised HOD model is fit and validated using the Dark Sage semi-analytic model, where we show that the HOD parameters can be modelled by simple linear and quadratic functions of HI mass. However, we also find that the clustering predicted by the HOD depends sensitively on the radial distributions of the HI galaxies within their host dark matter halos, which does not follow the NFW profile in the Dark Sage simulation. As such, this work enables -- for the first time -- a simple prescription for placing galaxies of different HI mass within dark matter halos in a way that is able to reproduce the HI mass-dependent galaxy clustering and HI mass function simultaneously and without requiring knowledge of the optical properties of the galaxies. Further efforts are required to demonstrate that this model can be used to produce large ensembles of mock galaxy catalogues for upcoming surveys., Comment: 18 pages, 13 figures, 1 table. Published in ApJ

We present a rest-frame optical morphological analysis of galaxies observed with the NIRCam imager on the James Webb Space Telescope (JWST) as part of the GLASS-JWST Early Release Science program. We select 388 sources at redshifts \(0.8 < z < 5.4\) and use the seven 0.9--5\micron\ NIRCam filters to generate rest-frame $gri$ composite color images, and conduct visual morphological classification. Compared to HST-based work we find a higher incidence of disks and bulges than expected at $z>1.5$, revealed by rest frame optical imaging. We detect 123 clear disks (58 at $z>1.5$) of which 76 have bulges. No evolution of bulge fraction with redshift is evident: 61\% at \(z<2\) (\(N=110\)) versus 60\% at \(z\geq2\) (\(N=13\)). A stellar mass dependence is evident, with bulges visible in 80\% of all disk galaxies with mass \(> 10^{9.5}\, {\rm M}_{\odot}\) (\(N=41\)) but only 52\% at \(M < 10^{9.5}\, {\rm M}_{\odot}\) (\(N=82\)). We supplement visual morphologies with non-parametric measurements of Gini and Asymmetry coefficients in the rest-frame $i$-band. Our sources are more asymmetric than local galaxies, with slightly higher Gini values. When compared to high-z rest-frame ultraviolet measurements with Hubble Space Telescope, JWST shows more regular morphological types such as disks, bulges and spiral arms at $z>1.5$, with smoother (i.e. lower Gini) and more symmetrical light distributions., Comment: Accepted for publication in ApJ Letters

Using the EAGLE suite of simulations, we demonstrate that both cold gas stripping {\it and} starvation of gas inflow play an important role in quenching satellite galaxies across a range of stellar and halo masses, $M_{\star}$ and $M_{200}$. By quantifying the balance between gas inflows, outflows, and star formation rates, we show that even at $z=2$, only $\approx30\%$ of satellite galaxies are able to maintain equilibrium or grow their reservoir of cool gas - compared to $\approx50\%$ of central galaxies at this redshift. We find that the number of orbits completed by a satellite is a very good predictor of its quenching, even more so than the time since infall. On average, we show that intermediate-mass satellites with $M_{\star}$ between $10^{9}{\rm M}_{\odot}-10^{10}{\rm M}_{\odot}$ will be quenched at first pericenter in massive group environments, $M_{200}>10^{13.5}{\rm M}_{\odot}$; and will be quenched at second pericenter in less massive group environments, $M_{200}<10^{13.5}{\rm M}_{\odot}$. On average, more massive satellites ($M_{\star}>10^{10}{\rm M}_{\odot}$) experience longer depletion time-scales, being quenched between first and second pericenters in massive groups; while in smaller group environments, just $\approx30\%$ will be quenched even after two orbits. Our results suggest that while starvation alone may be enough to slowly quench satellite galaxies, direct gas stripping, particularly at pericenters, is required to produce the short quenching time-scales exhibited in the simulation., Comment: Accepted for publication in MNRAS

In this VERTICO early science paper we explore in detail how environmental mechanisms, identified in HI, affect the resolved properties of molecular gas reservoirs in cluster galaxies. The molecular gas is probed using ALMA ACA (+TP) observations of 12CO(2-1) in 51 spiral galaxies in the Virgo cluster (of which 49 are detected), all of which are included in the VIVA HI survey. The sample spans a stellar mass range of 9 < log M*/Msol < 11. We study molecular gas radial profiles, isodensity radii, and surface densities as a function of galaxy HI deficiency and morphology. There is a weak correlation between global HI and H2 deficiencies, and resolved properties of molecular gas correlate with HI deficiency: galaxies that have large HI deficiencies have relatively steep and truncated molecular gas radial profiles, which is due to the removal of low-surface density molecular gas on the outskirts. Therefore, while the environmental mechanisms observed in HI also affect molecular gas reservoirs, there is only a moderate reduction of the total amount of molecular gas., Comment: Published in ApJ. 22 pages, 6 figures, 1 table, 1 appendix. Erratum accepted for publication in ApJ

Black hole feedback has been widely implemented as the key recipe to quench star formation in massive galaxies in modern semi-analytic models and hydrodynamical simulations. As the theoretical details surrounding the accretion and feedback of black holes continue to be refined, various feedback models have been implemented across simulations, with notable differences in their outcomes. Yet, most of these simulations have successfully reproduced some observations, such as stellar mass function and star formation rate density in the local Universe. We use the recent observation on the change of neutral hydrogen gas mass (including both ${\rm H_2}$ and ${\rm HI}$) with star formation rate of massive central disc galaxies as a critical constraint of black hole feedback models across several simulations. We find that the predictions of IllustrisTNG agree with the observations much better than the other models tested in this work. This favors IllustrisTNG's treatment of active galactic nuclei - where kinetic winds are driven by black holes at low accretion rates - as more plausible amongst those we test. In turn, this also indirectly supports the idea that the massive central disc galaxy population in the local Universe was likely quenched by AGN feedback., Comment: 15 pages, 3+4 figures, accepted for publication in ApJ

We present the Virgo Environment Traced in CO (VERTICO) survey, a new effort to map $^{12}$CO($2-1$), $^{13}$CO($2-1$), and C$^{18}$O($2-1$) in 51 Virgo Cluster galaxies with the Atacama Compact Array, part of the Atacama Large Millimeter/submillimeter Array (ALMA). The primary motivation of VERTICO is to understand the physical mechanisms that perturb molecular gas disks, and therefore star formation and galaxy evolution, in dense environments. This first paper contains an overview of VERTICO's design and sample selection, $^{12}$CO($2-1$) observations, and data reduction procedures. We characterize global $^{12}$CO($2-1$) fluxes and molecular gas masses for the 49 detected VERTICO galaxies, provide upper limits for the two non-detections, and produce resolved $^{12}$CO($2-1$) data products (median resolution $= 8^{\prime\prime} \approx 640~{\rm pc}$). Azimuthally averaged $^{12}$CO($2-1$) radial intensity profiles are presented along with derived molecular gas radii. We demonstrate the scientific power of VERTICO by comparing the molecular gas size--mass scaling relation for our galaxies with a control sample of field galaxies, highlighting the strong effect that radius definition has on this correlation. We discuss the drivers of the form and scatter in the size--mass relation and highlight areas for future work. VERTICO is an ideal resource for studying the fate of molecular gas in cluster galaxies and the physics of environment-driven processes that perturb the star formation cycle. Upon public release, the survey will provide a homogeneous legacy dataset for studying galaxy evolution in our closest cluster., Comment: 68 pages, 13 Figures, 2 Figure Sets, Accepted for publication in ApJS, Online FITS versions of Tables 1, 2, and 3 are available with the journal publication

We study the shapes of spatially integrated H I emission-line profiles of galaxies in the EAGLE simulation using three separate measures of the profile's asymmetry. We show that the subset of EAGLE galaxies whose gas fractions and stellar masses are consistent with those in the xGASS survey also have similar H I line asymmetries. Central galaxies with symmetric H I line profiles typically correspond to rotationally supported H I and stellar disks, but those with asymmetric line profiles may or may not correspond to dispersion-dominated systems. Galaxies with symmetric H I emission lines are, on average, more gas rich than those with asymmetric lines, and also exhibit systematic differences in their specific star formation rates, suggesting that turbulence generated by stellar or AGN feedback may be one factor contributing to H I line asymmetry. The line asymmetry also correlates strongly with the dynamical state of a galaxy's host dark matter halo: older, more relaxed haloes host more-symmetric galaxies than those hosted by unrelaxed ones. At fixed halo mass, asymmetric centrals tend to be surrounded by a larger number of massive subhaloes than their symmetric counterparts, and also experience higher rates of gas accretion and outflow. At fixed stellar mass, central galaxies have, on average, more symmetric H I emission lines than satellites; for the latter, ram pressure and tidal stripping are significant sources of asymmetry., Comment: Accepted for publication in MNRAS

We study the clustering of HI intensity maps produced from simulations with a focus on baryonic acoustic oscillations (BAO) and the effects induced by telescope beam smoothing and foreground cleaning. We start by creating a HI catalogue at $z=1.321$ based on the Semi-Analytic Galaxy Evolution (SAGE) model applied to the UNIT simulations. With this catalogue we investigate the relation between model HI and the dark matter haloes and we also study the abundance of HI, $\Omega_{\rm HI}$, predicted by this model. We then create synthetic HI intensity maps with a Nearest-Grid-Point approach. In order to simulate the telescope beam effect, a Gaussian smoothing is applied on the plane perpendicular to the line of sight. The effect of foreground removal methods is simulated by exponentially damping the largest wavelength Fourier modes on the radial direction. We study the anisotropic 2-point correlation function (2PCF) $\xi(r_\perp,r_\parallel)$ and how it is affected by the aforementioned observational effects. In order to better isolate the BAO signal, we study several 2PCF $\mu$-wedges (with a restricted range of orientations $\mu$) tailored to address the systematics effects and we compare them with different definitions of radial 2PCFs. Finally, we discuss our findings in the context of an SKA-like survey, finding a clear BAO signal in most of the estimators here proposed., Comment: published in MNRAS; 18 pages, 10 figures

New surveys such as ESA's Euclid mission are planned to map with unprecedented precision the large-scale structure of the Universe by measuring the 3D positions of tens of millions of galaxies. It is necessary to develop theoretically modelled galaxy catalogues to estimate the expected performance and to optimise the analysis strategy of these surveys. We populate two pairs of (1 Gpc/h)^3 volume dark-matter-only simulations from the UNIT project with galaxies using the SAGE semi-analytic model of galaxy formation, coupled to the photoionisation model get_emlines to estimate their H_alpha emission. These catalogues represent a unique suite that includes galaxy formation physics and - thanks to the fixed-pair technique used - an effective volume of ~(5 Gpc/h)^3, which is several times larger than the Euclid survey. We present the performance of these data and create five additional emission-line galaxy (ELG) catalogues by applying a dust attenuation model as well as adjusting the flux threshold as a function of redshift in order to reproduce Euclid-forecast dN/dz values. As a first application, we study the abundance and clustering of those model H_alpha ELGs: for scales greater than ~5 Mpc/h, we find a scale-independent bias with a value of $b\sim 1$ at redshift z ~ 0.5, that can increase nearly linearly to $b\sim 4$ at $z\sim 2$, depending on the ELG catalogue. Model galaxy properties, including their emission-line fluxes (with and without dust extinction) are publicly available., Comment: 16 pages, revised version as accepted by MNRAS

Measuring the HI-halo mass scaling relation (HIHM) is fundamental to understanding the role of HI in galaxy formation and its connection to structure formation. While direct measurements of the HI mass in haloes are possible using HI-spectral stacking, the reported shape of the relation depends on the techniques used to measure it (e.g. monotonically increasing with mass versus flat, mass-independent). Using a simulated HI and optical survey produced with the SHARK semi-analytic galaxy formation model, we investigate how well different observational techniques can recover the intrinsic, theoretically predicted, HIHM relation. We run a galaxy group finder and mimic the HI stacking procedure adopted by different surveys and find we can reproduce their observationally derived HIHM relation. However, none of the adopted techniques recover the underlying HIHM relation predicted by the simulation. We find that systematic effects in halo mass estimates of galaxy groups modify the inferred shape of the HIHM relation from the intrinsic one in the simulation, while contamination by interloping galaxies, not associated with the groups, contribute to the inferred HI mass of a halo mass bin, when using large velocity windows for stacking. The effect of contamination is maximal at Mvir~10^(12-12.5)Msol. Stacking methods based on summing the HI emission spectra to infer the mean HI mass of galaxies of different properties belonging to a group suffer minimal contamination but are strongly limited by the use of optical counterparts, which miss the contribution of dwarf galaxies. Deep spectroscopic surveys will provide significant improvements by going deeper while maintaining high spectroscopic completeness; for example, the WAVES survey will recover ~52% of the total HI mass of the groups with Mvir~10^(14)Msol compared to ~21% in GAMA., Comment: 20 pages,16 figures. Paper submitted to MNRAS

We examine how the post-processed content of molecular hydrogen (H$_2$) in galaxies from the TNG100 cosmological, hydrodynamic simulation changes with environment at $z\!=\!0$, assessing central/satellite status and host halo mass. We make close comparisons with the carbon monoxide (CO) emission survey xCOLD GASS where possible, having mock-observed TNG100 galaxies to match the survey's specifications. For a representative sample of host haloes across $10^{11}\!\lesssim\!M_{\rm 200c}/{\rm M}_{\odot}\!<\!10^{14.6}$, TNG100 predicts that satellites with $m_*\!\geq\!10^9\,{\rm M}_{\odot}$ should have a median deficit in their H$_2$ fractions of $\sim$0.6 dex relative to centrals of the same stellar mass. Once observational and group-finding uncertainties are accounted for, the signature of this deficit decreases to $\sim$0.2 dex. Remarkably, we calculate a deficit in xCOLD GASS satellites' H$_2$ content relative to centrals of 0.2--0.3 dex, in line with our prediction. We further show that TNG100 and SDSS data exhibit continuous declines in the average star formation rates of galaxies at fixed stellar mass in denser environments, in quantitative agreement with each other. By tracking satellites from their moment of infall in TNG100, we directly show that atomic hydrogen (HI) is depleted at fractionally higher rates than H$_2$ on average. Supporting this picture, we find that the H$_2$/HI mass ratios of satellites are elevated relative to centrals in xCOLD GASS. We provide additional predictions for the effect of environment on H$_2$ -- both absolute and relative to HI -- that can be tested with spectral stacking in future CO surveys., Comment: Accepted in MNRAS. Minor proofing and reference edits from v1. Main body: 19 pages, 11 figures

Observations of the cold neutral atomic hydrogen (HI) in and around disc galaxies have revealed that spatial and kinematic asymmetries are commonplace, and are reflected in the global HI spectra. We use the TNG100 box from the IllustrisTNG suite of cosmological simulations to study the conditions under which these asymmetries may arise in current theoretical galaxy formation models. We find that more than 50% of the sample has at least a 10% difference in integrated flux between the high- and low-velocity half of the spectrum, thus the typical TNG100 galaxy has an HI profile that is not fully symmetric. We find that satellite galaxies are a more asymmetric population than centrals, consistent with observational results. Using halo mass as a proxy for environment, this trend appears to be driven by the satellite population within the virial radius of haloes more massive than $10^{13} M_{\odot}$, typical of medium/large groups. We show that, while the excess of HI asymmetry in group satellites is likely driven by ram pressure, the bulk of the asymmetric HI profiles observed in TNG100 are driven by physical processes able to affect both the central and satellite populations. Our results highlight how asymmetries are not driven solely by environment, and multiple physical processes can produce the same asymmetric shape in global HI spectra., Comment: 16 pages, 13 figures, 1 appendix, accepted for publication in MRNRAS

Determining the spatial distribution and intrinsic physical properties of neutral hydrogen on cosmological scales is one of the key goals of next-generation radio surveys. We use the EAGLE galaxy formation simulations to assess the properties of damped Lyman-alpha absorbers (DLAs) that are associated with galaxies and their underlying dark matter haloes between 0 $\leq$ z $\leq$ 2. We find that the covering fraction of DLAs increases at higher redshift; a significant fraction of neutral atomic hydrogen (HI) resides in the outskirts of galaxies with stellar mass greater than or equal to 10$^{10}$ M$_\odot$; and the covering fraction of DLAs in the circumgalactic medium (CGM) is enhanced relative to that of the interstellar medium (ISM) with increasing halo mass. Moreover, we find that the mean density of the HI in galaxies increases with increasing stellar mass, while the DLAs in high- and low-halo-mass systems have higher column densities than those in galaxies with intermediate halo masses (~ 10$^{12}$ M$_\odot$ at z = 0). These high-impact CGM DLAs in high-stellar-mass systems tend to be metal-poor, likely tracing smooth accretion. Overall, our results point to the CGM playing an important role in DLA studies at high redshift (z $\geq$ 1). However, their properties are impacted both by numerical resolution and the detailed feedback prescriptions employed in cosmological simulations, particularly that of AGN., Comment: 25 pages. Accepted for publication in MNRAS

We use SHARK, a semi-analytic galaxy formation model, to investigate the physical processes involved in dictating the shape, scatter and evolution of the HI-halo mass relation at $0\leq z \leq 2$. We compare SHARK with HI clustering and spectral stacking of the HI-halo mass relation derived from observations finding excellent agreement with the former and a deficiency of HI in SHARK at $M_{\rm vir}\approx 10^{12-13} M_{\odot}$ in the latter, but otherwise great agreement below and above that mass threshold. In SHARK, we find that the HI mass increases with the halo mass up to a critical mass of $\approx 10^{11.8} M_{\odot}$; between $\sim 10^{11.8}-10^{13}M_{\odot}$, the scatter in the relation increases by 0.7 dex and the HI mass decreases with the halo mass on average; at $M_{\rm vir} \geq 10^{13} M_{\odot}$, the HI content continues to increase with halo mass. We find that the critical halo mass of $\approx 10^{12} M_{\odot}$ is largely set by feedback from Active Galactic Nuclei (AGN), and the exact shape and scatter of the HI-halo mass relation around that mass is extremely sensitive to how AGN feedback is modelled, with other physical processes playing a less significant role. We determine the main secondary parameters responsible for the scatter of the HI-halo mass relation, namely the halo spin parameter at $M_{\rm vir}\leq 10^{11.8} M_{\odot}$, and the fractional contribution from substructure to the total halo mass for $M_{\rm vir}\geq 10^{13} M_{\odot}$. The scatter at $10^{11.8}

We apply the spectral energy distribution (SED) fitting code ProSpect to multiwavelength imaging for $\sim$7,000 galaxies from the GAMA survey at $z<0.06$, in order to extract their star formation histories. We combine a parametric description of the star formation history with a closed-box evolution of metallicity where the present-day gas-phase metallicity of the galaxy is a free parameter. We show with this approach that we are able to recover the observationally determined cosmic star formation history (CSFH), an indication that stars are being formed in the correct epoch of the Universe, on average, for the manner in which we are conducting SED fitting. We also show the contribution to the CSFH of galaxies of different present-day visual morphologies, and stellar masses. Our analysis suggests that half of the mass in present-day elliptical galaxies was in place 11 Gyr ago. In other morphological types, the stellar mass formed later, up to 6 Gyr ago for present-day irregular galaxies. Similarly, the most massive galaxies in our sample were shown to have formed half their stellar mass by 11 Gyr ago, whereas the least massive galaxies reached this stage as late as 4 Gyr ago (the well-known effect of "galaxy downsizing"). Finally, our metallicity approach allows us to follow the average evolution in gas-phase metallicity for populations of galaxies, and extract the evolution of the cosmic metal mass density in stars and in gas, producing results in broad agreement with independent, higher-redshift observations of metal densities in the Universe., Comment: 24 pages, 16 figures, accepted for publication by MNRAS

We present a comparison of galaxy atomic and molecular gas properties in three recent cosmological hydrodynamic simulations, Simba, EAGLE, and Illustris-TNG, versus observations from $z\sim 0-2$. These simulations all rely on similar sub-resolution prescriptions to model cold interstellar gas which they cannot represent directly, and qualitatively reproduce the observed $z\approx 0$ HI and H$_2$ mass functions (HIMF, H2MF), CO(1-0) luminosity functions (COLF), and gas scaling relations versus stellar mass, specific star formation rate, and stellar surface density $\mu_*$, with some quantitative differences. To compare to the COLF, we apply an H$_2$-to-CO conversion factor to the simulated galaxies based on their average molecular surface density and metallicity, yielding substantial variations in $\alpha_{\rm CO}$ and significant differences between models. Using this, predicted $z=0$ COLFs agree better with data than predicted H2MFs. Out to $z\sim 2$, EAGLE's and Simba's HIMF and COLF strongly increase, while TNG's HIMF declines and COLF evolves slowly. EAGLE and Simba reproduce high $L_{\rm CO1-0}$ galaxies at $z\sim 1-2$ as observed, owing partly to a median $\alpha_{\rm CO}(z=2)\sim 1$ versus $\alpha_{\rm CO}(z=0)\sim 3$. Examining \HI, H$_2$, and CO scaling relations, their trends with $M_*$ are broadly reproduced in all models, but EAGLE yields too little HI in green valley galaxies, TNG and Simba overproduce cold gas in massive galaxies, and Simba overproduces molecular gas in small systems. Using Simba variants that exclude individual AGN feedback modules, we find that Simba's AGN jet feedback is primarily responsible by lowering cold gas contents from $z\sim 1\to0$ by suppressing cold gas in $M_*> 10^{10}{\rm M}_\odot$ galaxies, while X-ray feedback suppresses the formation of high-$\mu_*$ systems., Comment: 22 pages, MNRAS accepted

For astronomers to make a significant contribution to the reduction of climate change-inducing greenhouse gas emissions, we first must quantify our sources of emissions and review the most effective approaches for reducing them. Here we estimate that Australian astronomers' total greenhouse gas emissions from their regular work activities are $\gtrsim$25 ktCO$_2$-e/yr (equivalent kilotonnes of carbon dioxide per year). This can be broken into $\sim$15 ktCO$_2$-e/yr from supercomputer usage, $\sim$4.2 ktCO$_2$-e/yr from flights (where individuals' flight emissions correlate with seniority), $>$3.3 ktCO$_2$-e/yr from the operation of observatories, and 2.6$\pm$0.4 ktCO$_2$-e/yr from powering office buildings. Split across faculty scientists, postdoctoral researchers, and PhD students, this averages to $\gtrsim$37 tCO$_2$-e/yr per astronomer, over 40% more than what the average Australian non-dependant emits in total, equivalent to $\sim$5$\times$ the global average. To combat these environmentally unsustainable practices, we suggest astronomers should strongly preference use of supercomputers, observatories, and office spaces that are predominantly powered by renewable energy sources. Where facilities that we currently use do not meet this requirement, their funders should be lobbied to invest in renewables, such as solar or wind farms. Air travel should also be reduced wherever possible, replaced primarily by video conferencing, which should also promote inclusivity., Comment: Accepted Perspective in Nature Astronomy. Significant additions/updates to results from v1. Text refinement post-refereeing since v2. 9 pages (excl. refs), 2 figures

We study the present-day connection between galaxy morphology and angular momentum using the {\sc Dark Sage} semi-analytic model of galaxy formation. For galaxies between $ 10^{11}-10^{12} \mathrm{M}_{\odot}$ in stellar mass, the model successfully predicts the observed trend whereby galaxies with more prominent disks exhibit higher {\em stellar} disk specific angular momentum ($j_{\rm stellar, disk}$) at fixed stellar mass. However, when we include the gas in the disk, bulge-dominated galaxies have the highest {\em total} disk specific angular momentum ($j_{\rm total, disk}$). We attribute this to a large contribution from an extended disk of cold gas in typical bulge-dominated galaxies. We find the relationship between $j_{\rm dark matter}$ and morphology to be quite complex. Surprisingly, in this stellar mass range, not only do bulge-dominated galaxies tend to live in halos with higher $j_{\rm dark matter}$ than disk-dominated galaxies, but intermediate galaxies (those with roughly equal fractions of bulge and disk mass) have the lowest $j_{\rm dark matter}$ of all. Yet, when controlling for halo mass, rather than stellar mass, the relationship between $j_{\rm dark matter}$ and morphology vanishes. Based on these results, halo mass rather than angular momentum is the main driver of the predicted morphology sequence at high masses. In fact, in our stellar mass range, disk-dominated galaxies live in dark matter halos that are roughly 1/10th the mass of their bulge-dominated counterparts., Comment: 11 pages, 9 figures, submitted to ApJ

We analytically derive the observed size-mass relation of galaxies' atomic hydrogen (HI), including limits on its scatter, based on simple assumptions about the structure of HI discs. We trial three generic profiles for HI surface density as a function of radius. Firstly, we assert that HI surface densities saturate at a variable threshold, and otherwise fall off exponentially with radius or, secondly, radius squared. Our third model assumes the total gas surface density is exponential, with the HI fraction at each radius depending on local pressure. These are tested against a compilation of 110 galaxies from the THINGS, LITTLE THINGS, LVHIS, and Bluedisk surveys, whose HI surface density profiles are well resolved. All models fit the observations well and predict consistent size-mass relations. Using an analytical argument, we explain why processes that cause gas disc truncation - such as ram-pressure stripping - scarcely affect the HI size-mass relation. This is tested with the IllustrisTNG(100) cosmological, hydrodynamic simulation and the Dark Sage semi-analytic model of galaxy formation, both of which capture radially resolved disc structure. For galaxies with m_*>10^9 M_solar and m_HI>10^8 M_solar, both simulations predict HI size-mass relations that align with observations, show no difference between central and satellite galaxies, and show only a minor, second-order dependence on host halo mass for satellites. Ultimately, the universally tight HI size-mass relation is mathematically inevitable and robust. Only by completely disrupting the structure of HI discs, e.g. through overly powerful feedback, could a simulation predict the relation poorly., Comment: Accepted in MNRAS. Post-proofing. 15 pages and 11 figures in main body (references and appendices additional)

We have recently developed a post-processing framework to estimate the abundance of atomic and molecular hydrogen (HI and H2, respectively) in galaxies in large-volume cosmological simulations. Here we compare the HI and H2 content of IllustrisTNG galaxies to observations. We mostly restrict this comparison to $z \approx 0$ and consider six observational metrics: the overall abundance of HI and H2, their mass functions, gas fractions as a function of stellar mass, the correlation between H2 and star formation rate, the spatial distribution of gas, and the correlation between gas content and morphology. We find generally good agreement between simulations and observations, particularly for the gas fractions and the HI mass-size relation. The H2 mass correlates with star formation rate as expected, revealing an almost constant depletion time that evolves up to z = 2 as observed. However, we also discover a number of tensions with varying degrees of significance, including an overestimate of the total neutral gas abundance at z = 0 by about a factor of two and a possible excess of satellites with no or very little neutral gas. These conclusions are robust to the modelling of the HI/H2 transition. In terms of their neutral gas properties, the IllustrisTNG simulations represent an enormous improvement over the original Illustris run. All data used in this paper are publicly available as part of the IllustrisTNG data release., Comment: 22 pages, 11 figures. Additional figures and link to data at http://www.benediktdiemer.com/data/hi-h2-in-illustris/

We use TheThreeHundred project, a suite of 324 resimulated massive galaxy clusters embedded in a broad range of environments, to investigate (i) how the gas content of surrounding haloes correlates with phase-space position at $z=0$, and (ii) to investigate the role that ram pressure plays in this correlation. By stacking all 324 normalised phase-space planes containing 169287 haloes and subhaloes, we show that the halo gas content is tightly correlated with phase-space position. At $\sim\,1.5-2\,\text{R}_{\text{200}}$ of the cluster dark matter halo, we find an extremely steep decline in the halo gas content of infalling haloes and subhaloes irrespective of cluster mass, possibly indicating the presence of an accretion shock. We also find that subhaloes are particularly gas-poor, even in the cluster outskirts, which could indicate active regions of ongoing pre-processing. By modelling the instantaneous ram pressure experienced by each halo and subhalo at $z=0$, we show that the ram pressure intensity is also well correlated with phase-space position, which is again irrespective of cluster mass. In fact, we show that regions in the phase-space plane with high differential velocity between a halo or subhalo and its local gas environment, are almost mutually exclusive with high halo gas content regions. This suggests a causal link between the gas content of objects and the instantaneous ram pressure they experience, where the dominant factor is the differential velocity., Comment: 17 pages, 12 figures, accepted to MNRAS

We investigate the influence of environment on the cold-gas properties of galaxies at z=0 within the TNG100 cosmological, magnetohydrodynamic simulation, part of the IllustrisTNG suite. We extend previous post-processing methods for breaking gas cells into their atomic and molecular phases, and build detailed mocks to comprehensively compare to the latest surveys of atomic hydrogen (HI) in nearby galaxies, namely ALFALFA and xGASS. We use TNG100 to explore the HI content, star formation activity, and angular momentum of satellite galaxies, each as a function of environment, and find that satellites are typically a factor of ~3 poorer in HI than centrals of the same stellar mass, with the exact offset depending sensitively on parent halo mass. Due to the large physical scales on which HI measurements are made (~45--245 kpc), contributions from gas not bound to the galaxy of interest but in the same line of sight crucially lead to larger HI mass measurements in the mocks in many cases, ultimately aligning with observations. This effect is mass-dependent and naturally greater for satellites than centrals, as satellites are never isolated by definition. We also show that HI stripping in TNG100 satellites is closely accompanied by quenching, in tension with observational data that instead favour that HI is preferentially stripped before star formation is reduced., Comment: Published in MNRAS. Main body (full paper): 18 (22) pages, 10 (11) figures. New-found bug introduced in v4 mock plots fixed. BaryMP issue fixed per footnote in Dave et al. (2020). All changes are minor and do not affect text or conclusions

We introduce the THE THREE HUNDRED project, an endeavour to model 324 large galaxy clusters with full-physics hydrodynamical re-simulations. Here we present the data set and study the differences to observations for fundamental galaxy cluster properties and scaling relations. We find that the modelled galaxy clusters are generally in reasonable agreement with observations with respect to baryonic fractions and gas scaling relations at redshift z = 0. However, there are still some (model-dependent) differences, such as central galaxies being too massive, and galaxy colours (g - r) being bluer (about 0.2 dex lower at the peak position) than in observations. The agreement in gas scaling relations down to 10^{13} h^{-1} M_{\odot} between the simulations indicates that particulars of the sub-grid modelling of the baryonic physics only has a weak influence on these relations. We also include - where appropriate - a comparison to three semi-analytical galaxy formation models as applied to the same underlying dark-matter-only simulation. All simulations and derived data products are publicly available., Comment: 20 pages, 8 figures, 7 tables. MNRAS published version

We present a comparison of the observed evolving galaxy stellar mass functions with the predictions of eight semi-analytic models and one halo occupation distribution model. While most models are able to fit the data at low redshift, some of them struggle to simultaneously fit observations at high redshift. We separate the galaxies into 'passive' and 'star-forming' classes and find that several of the models produce too many low-mass star-forming galaxies at high redshift compared to observations, in some cases by nearly a factor of 10 in the redshift range $2.5 < z < 3.0$. We also find important differences in the implied mass of the dark matter haloes the galaxies inhabit, by comparing with halo masses inferred from observations. Galaxies at high redshift in the models are in lower mass haloes than suggested by observations, and the star formation efficiency in low-mass haloes is higher than observed. We conclude that many of the models require a physical prescription that acts to dissociate the growth of low-mass galaxies from the growth of their dark matter haloes at high redshift., Comment: 16 pages, 8 figures, accepted for publication in MNRAS

We explore the connection between the atomic gas fraction, f_atm, and `global disc stability' parameter, q, of galaxies within a fully cosmological context by examining galaxies in the Dark Sage semi-analytic model. The q parameter is determined by the ratio of disc specific angular momentum to mass. Dark Sage is well suited to our study, as it includes the numerical evolution of one-dimensional disc structure, making both j_disc and q predicted quantities. We show that Dark Sage produces a clear correlation between gas fraction and j_disc at fixed disc mass, in line with recent results from observations and hydrodynamic simulations. This translates to a tight q--f_atm sequence for star-forming central galaxies, which closely tracks the analytic prediction of Obreschkow et al. The scatter in this sequence is driven by the probability distribution function of mass as a function of j (PDF of j) within discs, specifically where it peaks. We find that halo mass is primarily responsible for the peak location of the PDF of j, at least for low values of q. Two main mechanisms of equal significance are then identified for disconnecting f_atm from q. Mergers in the model can trigger quasar winds, with the potential to blow out most of the gas disc, while leaving the stellar disc relatively unharmed. Ram-pressure stripping of satellite galaxies has a similar effect, where f_atm can drop drastically with only a minimal effect to q. We highlight challenges associated with following these predictions up with observations., Comment: 13 pages, 8 figures (excluding references and appendices). Submitted to MNRAS. Revisions after referee's report

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

We apply the Jeans Anisotropic MGE (JAM) dynamical modelling method to SAGES Legacy Unifying Globulars and GalaxieS (SLUGGS) survey data of early-type galaxies in the stellar mass range $10^{10}10^{10.7}{\rm M}_{\odot}$ in the Magneticum simulations, we identify a significant anticorrelation between total-mass density profile slopes and the fraction of stellar mass formed ex situ (i.e. accreted), whereas this anticorrelation is weaker for lower stellar masses, implying that the measured total mass density slopes for low-mass galaxies are less likely to be determined by merger activity., Comment: Published in MNRAS. ArXiv paper replaced to reflect updated total-mass density slope values for model 2, as described by the published erratum to the original paper

We present a comparison of nine galaxy formation models, eight semi-analytical and one halo occupation distribution model, run on the same underlying cold dark matter simulation (cosmological box of co-moving width 125$h^{-1}$ Mpc, with a dark-matter particle mass of $1.24\times 10^9 h^{-1}$ Msun) and the same merger trees. While their free parameters have been calibrated to the same observational data sets using two approaches, they nevertheless retain some 'memory' of any previous calibration that served as the starting point (especially for the manually-tuned models). For the first calibration, models reproduce the observed z = 0 galaxy stellar mass function (SMF) within 3-{\sigma}. The second calibration extended the observational data to include the z = 2 SMF alongside the z~0 star formation rate function, cold gas mass and the black hole-bulge mass relation. Encapsulating the observed evolution of the SMF from z = 2 to z = 0 is found to be very hard within the context of the physics currently included in the models. We finally use our calibrated models to study the evolution of the stellar-to-halo mass (SHM) ratio. For all models we find that the peak value of the SHM relation decreases with redshift. However, the trends seen for the evolution of the peak position as well as the mean scatter in the SHM relation are rather weak and strongly model dependent. Both the calibration data sets and model results are publicly available., Comment: 20 pages, 14 figures (+5 supplementary figures in the Appendix), accepted for publication in MNRAS

We present the public release of the MultiDark-Galaxies: three distinct galaxy catalogues derived from one of the Planck cosmology MultiDark simulations (i.e. MDPL2, with a volume of (1 Gpc/$h$)$^{3}$ and mass resolution of $1.5 \times 10^{9} M_{\odot}/h$) by applying the semi-analytic models GALACTICUS, SAG, and SAGE to it. We compare the three models and their conformity with observational data for a selection of fundamental properties of galaxies like stellar mass function, star formation rate, cold gas fractions, and metallicities - noting that they sometimes perform differently reflecting model designs and calibrations. We have further selected galaxy subsamples of the catalogues by number densities in stellar mass, cold gas mass, and star formation rate in order to study the clustering statistics of galaxies. We show that despite different treatment of orphan galaxies, i.e. galaxies that lost their dark-matter host halo due to the finite mass resolution of the N-body simulation or tidal stripping, the clustering signal is comparable, and reproduces the observations in all three models - in particular when selecting samples based upon stellar mass. Our catalogues provide a powerful tool to study galaxy formation within a volume comparable to those probed by on-going and future photometric and redshift surveys. All model data consisting of a range of galaxy properties - including broad-band SDSS magnitudes - are publicly available., Comment: 29 pages, 16 figures, 8 tables, accepted for publication in MNRAS. All data incl. the complete galaxy catalogues for all models are publicly available from the CosmoSim database (http://www.cosmosim.org); a selected set of galaxy properties is available via the Skies & Universes website (http://www.skiesanduniverses.org)

We combine the latest spectrally stacked data of 21-cm emission from the ALFALFA survey with an updated version of the Dark Sage semi-analytic model to investigate the relative contributions of secular and environmental astrophysical processes on shaping the HI fractions and quiescence of galaxies in the local Universe. We calibrate the model to match the observed mean HI fraction of all galaxies as a function of stellar mass. Without consideration of stellar feedback, disc instabilities, and active galactic nuclei, we show how the slope and normalisation of this relation would change significantly. We find Dark Sage can reproduce the relative impact that halo mass is observed to have on satellites' HI fractions and quiescent fraction. However, the model satellites are systematically gas-poor. We discuss how this could be affected by satellite--central cross-contamination from the group-finding algorithm applied to the observed galaxies, but that it is not the full story. From our results, we suggest the anti-correlation between satellites' HI fractions and host halo mass, seen at fixed stellar mass and fixed specific star formation rate, can be attributed almost entirely to ram-pressure stripping of cold gas. Meanwhile, stripping of hot gas from around the satellites drives the correlation of quiescent fraction with halo mass at fixed stellar mass. Further detail in the modelling of galaxy discs' centres is required to solidify this result, however. We contextualise our results with those from other semi-analytic models and hydrodynamic simulations., Comment: 18 pages, 10 figures, accepted for publication in MNRAS. Correction to errorbars in Figs 3,4,9

We use EAGLE to quantify the effect galaxy mergers have on the stellar specific angular momentum of galaxies, $j_{\rm stars}$. We split mergers into: dry (gas-poor)/wet (gas-rich), major/minor, and different spin alignments and orbital parameters. Wet (dry) mergers have an average neutral gas-to-stellar mass ratio of $1.1$ ($0.02$), while major (minor) mergers are those with stellar mass ratios $\ge 0.3$ ($0.1-0.3$). We correlate the positions of galaxies in the $j_{\rm stars}$-stellar mass plane at $z=0$ with their merger history, and find that galaxies of low spins suffered dry mergers, while galaxies of normal/high spins suffered predominantly wet mergers, if any. The radial $j_{\rm stars}$ profiles of galaxies that went through dry mergers are deficient by $\approx 0.3$~dex at $r\lesssim 10\,r_{50}$ (with $r_{50}$ being the half-stellar mass radius), compared to galaxies that went through wet mergers. Studying the merger remnants reveals that dry mergers reduce $j_{\rm stars}$ by $\approx 30$\%, while wet mergers increase it by $\approx 10$\%, on average. The latter is connected to the build-up of the bulge by newly formed stars of high rotational speed. Moving from minor to major mergers accentuates these effects. When the spin vectors of the galaxies prior to the dry merger are misaligned, $j_{\rm stars}$ decreases to a greater magnitude, while in wet mergers co-rotation and high orbital angular momentum efficiently spun-up galaxies. We predict what would be the observational signatures in the $j_{\rm stars}$ profiles driven by dry mergers: (i) shallow radial profiles and (ii) profiles that rise beyond $\approx 10\,r_{50}$, both of which are significantly different from spiral galaxies., Comment: Accepted for publication in MNRAS. Compared to previous revision, Fig. 6 and A1 are new, and section 3.3 was shortened

We use the EAGLE cosmological hydrodynamic simulation suite to study the specific angular momentum of galaxies, $j$, with the aims of (i) investigating the physical causes behind the wide range of $j$ at fixed mass and (ii) examining whether simple, theoretical models can explain the seemingly complex and non-linear nature of the evolution of $j$. We find that $j$ of the stars, $j_{\rm stars}$, and baryons, $j_{\rm bar}$, are strongly correlated with stellar and baryon mass, respectively, with the scatter being highly correlated with morphological proxies such as gas fraction, stellar concentration, (u-r) intrinsic colour, stellar age and the ratio of circular velocity to velocity dispersion. We compare with available observations at $z=0$ and find excellent agreement. We find that $j_{\rm bar}$ follows the theoretical expectation of an isothermal collapsing halo under conservation of specific angular momentum to within $\approx 50$%, while the subsample of rotation-supported galaxies are equally well described by a simple model in which the disk angular momentum is just enough to maintain marginally stable disks. We extracted evolutionary tracks of the stellar spin parameter of EAGLE galaxies and found that the fate of their $j_{\rm stars}$ at $z=0$ depends sensitively on their star formation and merger histories. From these tracks, we identified two distinct physical channels behind low $j_{\rm stars}$ galaxies at $z=0$: (i) galaxy mergers, and (ii) early star formation quenching. The latter can produce galaxies with low $j_{\rm stars}$ and early-type morphologies even in the absence of mergers., Comment: 22 pages (16 without appendices). Accepted for publication in MNRAS. The only difference with previous version is a shorter abstract (to comply with the 250 words limit in MNRAS)

We use the hydrodynamic, cosmological EAGLE simulations to investigate how hot gas in haloes condenses to form and grow galaxies. We select haloes from the simulations that are actively cooling and study the temperature, distribution, and metallicity of their hot, cold, and transitioning `cooling' gas, placing these in context of semi-analytic models. Our selection criteria lead us to focus on Milky Way-like haloes. We find the hot-gas density profiles of the haloes form a progressively stronger core over time, the nature of which can be captured by a beta profile that has a simple dependence on redshift. In contrast, the hot gas that will cool over a time-step is broadly consistent with a singular isothermal sphere. We find that cooling gas carries a few times the specific angular momentum of the halo and is offset in spin direction from the rest of the hot gas. The gas loses ~60% of its specific angular momentum during the cooling process, generally remaining greater than that of the halo, and it precesses to become aligned with the cold gas already in the disc. We find tentative evidence that angular-momentum losses are slightly larger when gas cools onto dispersion-supported galaxies. We show that an exponential surface density profile for gas arriving on a disc remains a reasonable approximation, but a cusp containing ~20% of the mass is always present, and disc scale radii are larger than predicted by a vanilla Fall & Efstathiou model. These scale radii are still closely correlated with the halo spin parameter, for which we suggest an updated prescription for galaxy formation models., Comment: Accepted in MNRAS. 21 pages, 18 figures. Replacement includes fixed typos and reference updates

We present the new semi-analytic model of galaxy evolution, DARK SAGE, a heavily modified version of the publicly available SAGE code. The model is designed for detailed evolution of galactic discs. We evolve discs in a series of annuli with fixed specific angular momentum, which allows us to make predictions for the radial and angular-momentum structure of galaxies. Most physical processes, including all channels of star formation and associated feedback, are performed in these annuli. We present the surface density profiles of our model spiral galaxies, both as a function of radius and specific angular momentum, and find the discs naturally build a pseduobulge-like component. Our main results are focussed on predictions relating to the integrated mass--specific angular momentum relation of stellar discs. The model produces a distinct sequence between these properties in remarkable agreement with recent observational literature. We investigate the impact Toomre disc instabilities have on shaping this sequence and find they are crucial for regulating both the mass and spin of discs. Without instabilities, high-mass discs would be systematically deficient in specific angular momentum by a factor of ~2.5, with increased scatter. Instabilities also appear to drive the direction in which the mass--spin sequence of spiral galaxy discs evolves. With them, we find galaxies of fixed mass have higher specific angular momentum at later epochs., Comment: Published in MNRAS, 19 pages, 20 figures (both including appendices)