Researcher profile

Danail Obreschkow

Danail Obreschkow contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
19works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

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

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

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

Building this graph slice

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

Published work

19 published item(s)

preprint2022arXiv

An empirical measurement of the Halo Mass Function from the combination of GAMA DR4, SDSS DR12, and REFLEX II data

We construct the halo mass function (HMF) from the GAMA galaxy group catalogue over the mass range 10^12.7M_sol to 10^15.5M_sol, and find good agreement with the expectation from LambdaCDM. In comparison to previous studies, this result extends the mass range over which the HMF has now been measured over by an order of magnitude. We combine the GAMA DR4 HMF with similar data from the SDSS DR12 and REFLEX II surveys, and fit a four-parameter Murray-Robotham-Power (MRP) function, valid at z~0.1, yielding: a density normalisation of: log10 (phi Mpc^3)=-3.96[+0.55,-0.82], a high mass turn-over of: log10(M/M_sol)=14.13[+0.43,-0.40], a low mass power law slope of: alpha=-1.68[+0.21,-0.24] , and a high mass softening parameter of: beta= 0.63[+0.25,-0.11]. If we fold in the constraint on Omega_M from Planck 2018 Cosmology, we are able to reduce these uncertainties further, but this relies on the assumption that the power-law trend can be extrapolated from 10^12.7M_sol to zero mass. Throughout, we highlight the effort needed to improve on our HMF measurement: improved halo mass estimates that do not rely on calibration to simulations; reduced halo mass uncertainties needed to mitigate the strong Eddington Bias that arises from the steepness of the HMF low mass slope; and deeper wider area spectroscopic surveys. To our halo mass limit of 10^12.7 M_sol, we are directly resolving (`seeing') 41+/-5 per cent of the total mass density, i.e. Omega_[M>12.7]=0.128+/-0.016, opening the door for the direct construction of 3D dark matter mass maps at Mpc resolution.

preprint2022arXiv

Galaxy And Mass Assembly (GAMA): Data Release 4 and the z < 0.1 total and z < 0.08 morphological galaxy stellar mass functions

In Galaxy And Mass Assembly Data Release 4 (GAMA DR4), we make available our full spectroscopic redshift sample. This includes 248682 galaxy spectra, and, in combination with earlier surveys, results in 330542 redshifts across five sky regions covering ~250deg^2. The redshift density, is the highest available over such a sustained area, has exceptionally high completeness (95 per cent to r_KIDS=19.65mag), and is well suited for the study of galaxy mergers, galaxy groups, and the low redshift (z<0.25) galaxy population. DR4 includes 32 value-added tables or Data Management Units (DMUs) that provide a number of measured and derived data products including GALEX, ESO KiDS, ESO VIKING, WISE and Herschel Space Observatory imaging. Within this release, we provide visual morphologies for 15330 galaxies to z<0.08, photometric redshift estimates for all 18million objects to r_KIDS~25mag, and stellar velocity dispersions for 111830 galaxies. We conclude by deriving the total galaxy stellar mass function (GSMF) and its sub-division by morphological class (elliptical, compact-bulge and disc, diffuse-bulge and disc, and disc only). This extends our previous measurement of the total GSMF down to 10^6.75 M_sol h^-2_70 and we find a total stellar mass density of rho_*=(2.97+/-0.04)x10^8 M_sol h_70 Mpc^-3 or Omega_*=(2.17+/-0.03)x10^-3 h^-1_70. We conclude that at z<0.1, the Universe has converted 4.9+/-0.1 per cent of the baryonic mass implied by Big Bang Nucleosynthesis into stars that are gravitationally bound within the galaxy population.

preprint2022arXiv

LADUMA: Discovery of a luminous OH megamaser at $z > 0.5$

In the local Universe, OH megamasers (OHMs) are detected almost exclusively in infrared-luminous galaxies, with a prevalence that increases with IR luminosity, suggesting that they trace gas-rich galaxy mergers. Given the proximity of the rest frequencies of OH and the hyperfine transition of neutral atomic hydrogen (HI), radio surveys to probe the cosmic evolution of HI in galaxies also offer exciting prospects for exploiting OHMs to probe the cosmic history of gas-rich mergers. Using observations for the Looking At the Distant Universe with the MeerKAT Array (LADUMA) deep HI survey, we report the first untargeted detection of an OHM at $z > 0.5$, LADUMA J033046.20$-$275518.1 (nicknamed &#34;Nkalakatha&#34;). The host system, WISEA J033046.26$-$275518.3, is an infrared-luminous radio galaxy whose optical redshift $z \approx 0.52$ confirms the MeerKAT emission line detection as OH at a redshift $z_{\rm OH} = 0.5225 \pm 0.0001$ rather than HI at lower redshift. The detected spectral line has 18.4$σ$ peak significance, a width of $459 \pm 59\,{\rm km\,s^{-1}}$, and an integrated luminosity of $(6.31 \pm 0.18\,{\rm [statistical]}\,\pm 0.31\,{\rm [systematic]}) \times 10^3\,L_\odot$, placing it among the most luminous OHMs known. The galaxy&#39;s far-infrared luminosity $L_{\rm FIR} = (1.576 \pm 0.013) \times 10^{12}\,L_\odot$ marks it as an ultra-luminous infrared galaxy; its ratio of OH and infrared luminosities is similar to those for lower-redshift OHMs. A comparison between optical and OH redshifts offers a slight indication of an OH outflow. This detection represents the first step towards a systematic exploitation of OHMs as a tracer of galaxy growth at high redshifts.

preprint2022arXiv

Spin transfer from dark matter to gas during halo formation

In the protogalactic density field, diffuse gas and collision-less cold dark matter (DM) are often assumed sufficiently mixed that both components experience identical tidal torques. However, haloes in cosmological simulations consistently end up with a higher specific angular momentum (sAM) in gas, even in simulations without radiative cooling and galaxy formation physics. We refine this result by analysing the spin distributions of gas and DM in $\sim$50,000 well-resolved haloes in a non-radiative cosmological simulation from the SURFS suite. The sAM of the halo gas on average ends up $\sim$40\% above that of the DM. This can be pinned down to an excess AM in the inner halo ($<$50\% virial radius), paralleled by a more coherent rotation pattern in the gas. We uncover the leading driver for this AM difference through a series of control simulations of a collapsing ellipsoidal top-hat, where gas and DM are initially well mixed. These runs reveal that the pressurised inner gas shells collapse more slowly, causing the DM ellipsoid to spin ahead of the gas ellipsoid. The arising torque generally transfers AM from the DM to the gas. The amount of AM transferred via this mode depends on the initial spin, the initial axes ratios and the collapse factor. These quantities can be combined in a single dimensionless parameter, which robustly predicts the AM transfer of the ellipsoidal collapse. This simplistic model can quantitatively explain the average AM excess of the gas found in the more complex non-radiative cosmological simulation.

preprint2022arXiv

Stellar Masses of Clumps in Gas-rich, Turbulent Disk Galaxies

In this paper we use HST/WFC3 observations of 6 galaxies from the DYNAMO survey, combined with stellar population modelling of the SED, to determine the stellar masses of DYNAMO clumps. The DYNAMO sample has been shown to have properties similar to $z\approx1.5$ turbulent, clumpy disks. DYNAMO sample clump masses offer a useful comparison for studies of $z>1$ in that the galaxies have the same properties, yet the observational biases are significantly different. Using DYNAMO we can more easily probe rest-frame near-IR wavelengths and also probe finer spatial scales. We find that the stellar mass of DYNAMO clumps is typically $10^{7}-10^8 \mathrm{M}_\odot$. We employ a technique that makes non-parametric corrections in removal of light from nearby clumps, and carries out a locally determined disk subtraction. The process of disk subtraction is the dominant effect, and can alter clump masses at the 0.3~dex level. Using these masses, we investigate the stellar mass function of clumps in DYNAMO galaxies. DYNAMO stellar mass functions follow a declining power law with slope $α\approx -1.4$, which is slightly shallower than, but similar to what is observed in $z>1$ lensed galaxies. We compare DYNAMO clump masses to results of simulations. The masses and galactocentric position of clumps in DYNAMO galaxies are more similar to long-lived clumps in simulations. Similar to recent DYNAMO results on the stellar population gradients, these results are consistent with simulations that do not employ strong &#34;early&#34; radiative feedback prescriptions.

preprint2022arXiv

The variation of the gas content of galaxy groups and pairs compared to isolated galaxies

We measure how the atomic gas (HI) fraction ($f_{HI}={\rm \frac{M_{HI}}{M_{*}}}$) of groups and pairs taken as single units vary with average stellar mass ($\langle {\rm M_*} \rangle$) and average star-formation rate ($\langle {\rm SFR} \rangle$), compared to isolated galaxies. The HI 21 cm emission observation are from (i) archival ALFALFA survey data covering three fields from the GAMA survey (provides environmental and galaxy properties), and (ii) DINGO pilot survey data of one of those fields. The mean $f_{HI}$ for different units (groups/pairs/isolated galaxies) are measured in regions of the log($\langle {\rm M_*} \rangle$) -- log($\langle {\rm SFR} \rangle$) plane, relative to the z $\sim 0$ star-forming main sequence (SFMS) of individual galaxies, by stacking $f_{HI}$ spectra of individual units. For ALFALFA, $f_{HI}$ spectra of units are measured by extracting HI spectra over the full groups/pair areas and dividing by the total stellar mass of member galaxies. For DINGO, $f_{HI}$ spectra of units are measured by co-adding HI spectra of individual member galaxies, followed by division by their total stellar mass. For all units the mean $f_{HI}$ decreases as we move to higher $\langle {\rm M_*} \rangle$ along the SFMS, and as we move from above the SFMS to below it at any $\langle {\rm M_*} \rangle$. From the DINGO-based study, mean $f_{HI}$ in groups appears to be lower compared to isolated galaxies for all $\langle {\rm M_*} \rangle$ along the SFMS. From the ALFALFA-based study we find substantially higher mean $f_{HI}$ in groups compared to isolated galaxies (values for pairs being intermediate) for ${\langle{\rm M_*}\rangle}\lesssim10^{9.5}~{\rm M_{\odot}}$, indicating the presence of substantial amounts of HI not associated with cataloged member galaxies in low mass groups.

preprint2022arXiv

xGASS: The connection between angular momentum, mass and atomic gas fraction in nearby galaxies

We use a sample of 559 disc galaxies extracted from the eXtended GALEX Arecibo SDSS Survey (xGASS) to study the connection between baryonic angular momentum, mass and atomic gas fraction in the local Universe. Baryonic angular momenta are determined by combining HI and H$_{2}$ integrated profiles with two-dimensional stellar mass surface density profiles. In line with previous work, we confirm that specific angular momentum and atomic gas fraction are tightly correlated, but we find a larger scatter than previously observed. This is most likely due to the wider range of galaxy properties covered by our sample. We compare our findings with the predictions of the analytical stability model developed by Obreschkow et al. and find that, while the model provides a very good first-order approximation for the connection between baryonic angular momentum, mass and gas fraction, it does not fully match our data. Specifically, we find that at fixed baryonic mass, the dependence of specific angular momentum on gas fraction is significantly weaker, and at fixed gas fraction, the slope of the angular momentum vs. mass relation is shallower than what was predicted by the model. The reasons behind this tension remain unclear, but we speculate that multiple factors may simultaneously play a role, all related to the fact that the model is not able to encapsulate the full diversity of galaxy properties in our sample.

preprint2021arXiv

GAMA/DEVILS: Constraining the cosmic star-formation history from improved measurements of the 0.3-2.2 micron Extragalactic Background Light

We present a revised measurement of the optical extragalactic background light (EBL), based on the contribution of resolved galaxies to the integrated galaxy light (IGL). The cosmic optical background radiation (COB), encodes the light generated by star-formation, and provides a wealth of information about the cosmic star formation history (CSFH). We combine wide and deep galaxy number counts from the Galaxy And Mass Assembly survey (GAMA) and Deep Extragalactic VIsible Legacy Survey (DEVILS), along with the Hubble Space Telescope (HST) archive and other deep survey datasets, in 9 multi-wavelength filters to measure the COB in the range from 0.35 micron to 2.2 micron. We derive the luminosity density in each band independently and show good agreement with recent and complementary estimates of the optical-EBL from very high-energy (VHE) experiments. Our error analysis suggests that the IGL and Gamma-ray measurements are now fully consistent to within ~10%, suggesting little need for any additional source of diffuse light beyond the known galaxy population. We use our revised IGL measurements to constrain the cosmic star-formation history, and place amplitude constraints on a number of recent estimates. As a consistency check, we can now demonstrate convincingly, that the CSFH, stellar mass growth, and the optical-EBL provide a fully consistent picture of galaxy evolution. We conclude that the peak of star-formation rate lies in the range 0.066-0.076 Msol/yr/Mpc^3 at a lookback time of 9.1 to 10.9 Gyrs.

preprint2021arXiv

K-Band Imaging of the Nearby, Clumpy Turbulent Disk Galaxy DYNAMO G04-1

We present a case study of stellar clumps in G04-1, a clumpy, turbulent disk galaxy located at $z$ = 0.13 from the DYNAMO sample, using adaptive optics enabled K-band imaging ($\sim2.25$ kpc/arcsec) with Keck/NIRC2. We identify 15 stellar clumps in G04-1 with a range of masses from $3.6 \times 10^{6}$ to $2.7 \times 10^{8}\ \rm M_{\odot}$, and with a median mass of $\sim2.9 \times 10^{7}\ \rm M_{\odot}$. Note that these masses decrease by about one-half when we apply a light correction for the underlying stellar disk. A majority (12 of 15) of clumps observed in the $K_{P}$-band imaging have associated components in H$α$ maps ($\sim2.75$ kpc/arcsec; $<$R$_{clump}> \sim$500 pc) and appear co-located ($\overline{Δx} \sim 0.1$ arcsec). Using Hubble Space Telescope WFC/ACS observations with the F336W and F467M filters, we also find evidence of radial trends in clump stellar properties: clumps closer to the centre of G04-1 are more massive (consistent with observations at high-$z$) and appear more red, suggesting they may be more evolved. Using our high-resolution data, we construct a star forming main sequence for G04-1 in terms of spatially-resolved quantities and find that all regions (both clump and intra-clump) within the galaxy are experiencing an enhanced mode of star formation routinely observed in galaxies at high-$z$. In comparison to recent simulations, our observation of a number of clumps with masses $10^{7}-10^{8}\ \rm M_{\odot}$ is not consistent with strong radiative feedback in this galaxy.

preprint2021arXiv

The SAMI Galaxy Survey: the third and final data release

We have entered a new era where integral-field spectroscopic surveys of galaxies are sufficiently large to adequately sample large-scale structure over a cosmologically significant volume. This was the primary design goal of the SAMI Galaxy Survey. Here, in Data Release 3 (DR3), we release data for the full sample of 3068 unique galaxies observed. This includes the SAMI cluster sample of 888 unique galaxies for the first time. For each galaxy, there are two primary spectral cubes covering the blue (370-570nm) and red (630-740nm) optical wavelength ranges at spectral resolving power of R=1808 and 4304 respectively. For each primary cube, we also provide three spatially binned spectral cubes and a set of standardized aperture spectra. For each galaxy, we include complete 2D maps from parameterized fitting to the emission-line and absorption-line spectral data. These maps provide information on the gas ionization and kinematics, stellar kinematics and populations, and more. All data are available online through Australian Astronomical Optics (AAO) Data Central.

preprint2020arXiv

Angular momentum-related probe of cold gas deficiencies

Recent studies of neutral atomic hydrogen (HI) in nearby galaxies found that all field disk galaxies are HI saturated, in that they carry roughly as much HI as permitted before this gas becomes gravitationally unstable. By taking this HI saturation for granted, the atomic gas fraction $f_{\rm atm}$ of galactic disks can be predicted as a function of the stability parameter $q=jσ/(GM)$, where $M$ and $j$ are the baryonic mass and specific angular momentum of the disk and $σ$ is the HI velocity dispersion Obreschkow et al. 2016. The log-ratio $Δf_q$ between this predictor and the observed atomic fraction can be seen as a physically motivated `HI deficiency&#39;. While field disk galaxies have $Δf_q \approx0$, objects subject to environmental removal of HI are expected to have $Δf_q>0$. Within this framework, we revisit the HI deficiencies of satellite galaxies in the Virgo cluster and in clusters of the EAGLE simulation. We find that observed and simulated cluster galaxies are HI deficient and that $Δf_q$ slightly increases when getting closer to the cluster centres. The $Δf_q$ values are similar to traditional HI deficiency estimators, but $Δf_q$ is more directly comparable between observations and simulations than morphology-based deficiency estimators. By tracking the simulated HI deficient cluster galaxies back in time, we confirm that $Δf_q\approx0$ until the galaxies first enter a halo with $M_{\rm halo}>10^{13} {\rm M_{\odot}}$, at which moment they quickly lose HI by environmental effects. Finally, we use the simulation to investigate the links between $Δf_q$ and quenching of star formation.

preprint2020arXiv

Characterising the Structure of Halo Merger Trees Using a Single Parameter: The Tree Entropy

Linking the properties of galaxies to the assembly history of their dark matter haloes is a central aim of galaxy evolution theory. This paper introduces a dimensionless parameter $s\in[0,1]$, the &#34;tree entropy&#34;, to parametrise the geometry of a halo&#39;s entire mass assembly hierarchy, building on a generalisation of Shannon&#39;s information entropy. By construction, the minimum entropy ($s=0$) corresponds to smoothly assembled haloes without any mergers. In contrast, the highest entropy ($s=1$) represents haloes grown purely by equal-mass binary mergers. Using simulated merger trees extracted from the cosmological $N$-body simulation SURFS, we compute the natural distribution of $s$, a skewed bell curve peaking near $s=0.4$. This distribution exhibits weak dependences on halo mass $M$ and redshift $z$, which can be reduced to a single dependence on the relative peak height $δ_{\rm c}/σ(M,z)$ in the matter perturbation field. By exploring the correlations between $s$ and global galaxy properties generated by the SHARK semi-analytic model, we find that $s$ contains a significant amount of information on the morphology of galaxies $-$ in fact more information than the spin, concentration and assembly time of the halo. Therefore, the tree entropy provides an information-rich link between galaxies and their dark matter haloes.

preprint2020arXiv

From rest-frame luminosity functions to observer-frame colour distributions: tackling the next challenge in cosmological simulations

Galaxy spectral energy distributions (SEDs) remain among the most challenging yet informative quantities to reproduce in simulations due to the large and complex mixture of physical processes that shape the radiation output of a galaxy. With the increasing number of surveys utilising broadband colours as part of their target selection criteria, the production of realistic SEDs in simulations is necessary for assisting in survey design and interpretation of observations. The recent success in reproducing the observed luminosity functions (LF) from far-UV to far-IR, using the state-of-the-art semi-analytic model \shark\ and the SED generator \prospect, represents a critical step towards better galaxy colour predictions. We show that with \shark\ and \prospect\ we can closely reproduce the optical colour distributions observed in the panchromatic GAMA survey. The treatment of feedback, star formation, central-satellite interactions and radiation re-processing by dust are critical for this achievement. The first three processes create a bimodal distribution, while dust attenuation defines the location and shape of the blue and red populations. While a naive comparison between observation and simulations displays the known issue of over-quenching of satellite galaxies, the introduction of empirically-motivated observational errors and classification from the same group finder used in GAMA greatly reduces this tension. The introduction of random re-assignment of $\sim 15\%$ of centrals/satellites as satellites/centrals on the simulation classification closely resembles the outcome of the group finder, providing a computationally less intensive method to compare simulations with observations.

preprint2020arXiv

Multi-wavelength consensus of large-scale linear bias

We model the large-scale linear galaxy bias $b_g(x,z)$ as a function of redshift $z$ and observed absolute magnitude threshold $x$ for broadband continuum emission from the far infrared to ultra-violet, as well as for prominent emission lines, such as the H$α$, H$β$, Lya and [OII] lines. The modelling relies on the semi-analytic galaxy formation model GALFORM, run on the state-of-the-art $N$-body simulation SURFS with the Planck 2015 cosmology. We find that both the differential bias at observed absolute magnitude $x$ and the cumulative bias for magnitudes brighter than $x$ can be fitted with a five-parameter model: $b_g(x,z)=a + b(1+z)^e(1 + \exp{[(x-c)d]})$. We also find that the bias for the continuum bands follows a very similar form regardless of wavelength due to the mixing of star-forming and quiescent galaxies in a magnitude limited survey. Differences in bias only become apparent when an additional colour separation is included, which suggest extensions to this work could look at different colours at fixed magnitude limits. We test our fitting formula against observations, finding reasonable agreement with some measurements within $1σ$ statistical uncertainties, and highlighting areas of improvement. We provide the fitting parameters for various continuum bands, emission lines and intrinsic galaxy properties, enabling a quick estimation of the linear bias in any typical survey of large-scale structure.

preprint2020arXiv

Stellar angular momentum distribution linked to galaxy morphology

We study the spatially-resolved stellar specific angular momentum $j_*$ in a high-quality sample of 24 CALIFA galaxies covering a broad range of visual morphology, accounting for stellar velocity and velocity dispersion. The shape of the spaxel-wise probability density function of normalised $s=j_*/j_{*mean}$, PDF($s$), deviates significantly from the near-universal initial distribution expected of baryons in a dark matter halo and can be explained by the expected baryonic effects in galaxy formation that remove and redistribute angular momentum. Further we find that the observed shape of the PDF($s$) correlates significantly with photometric morphology, where late-type galaxies have PDF($s$) that is similar to a normal distribution, whereas early types have a strongly-skewed PDF($s$) resulting from an excess of low-angular momentum material. Galaxies that are known to host pseudobulges (bulge Sersic index $n_b <2.2$) tend to have less skewed bulge PDF($s$), with skewness $(b_{1rb})\lesssim0.8$. The PDF($s$) encodes both kinematic and photometric information and appears to be a robust tracer of morphology. Its use is motivated by the desire to move away from traditional component-based classifications which are subject to observer bias, to classification on a galaxy&#39;s fundamental (stellar mass, angular momentum) properties. In future, PDF($s$) may also be useful as a kinematic decomposition tool.

preprint2020arXiv

The influence of angular momentum and environment on the HI gas of late-type galaxies

We use high-resolution HI data from the WHISP survey to study the HI and angular momentum properties of a sample of 114 late-type galaxies. We explore the specific baryonic angular momentum -- baryonic mass ($j_b - M_b$) relation, and find that an unbroken power law of the form $j_b \propto M_b^{0.55 \pm 0.02}$ fits the data well, with an intrinsic scatter of $\sim 0.13 \pm 0.01$ dex. We revisit the relation between the atomic gas fraction, $f_{atm}$, and the integrated atomic stability parameter $q$ (the $f_{atm} - q$ relation), originally introduced by Obreschkow et al., and probe this parameter space by populating it with galaxies from different environments, in order to study the influence of the environment on their $j_b$, $f_{atm}$ and $q$ values. We find evidence that galaxies with close neighbours show a larger intrinsic scatter about the $f_{atm} - q$ relation compared to galaxies without close-neighbours. We also find enhanced SFR among the deviating galaxies with close neighbours. In addition, we use the bulge-to-total (B/T) ratio as a morphology proxy, and find a general trend of decreasing B/T values with increasing disc stability and HI fraction in the $f_{atm} - q$ plane, indicating a fundamental link between mass, specific angular momentum, gas fraction and morphology of galaxies.

preprint2020arXiv

The physical drivers of the atomic hydrogen-halo mass relation

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}<M_{\rm vir}<10^{13} M_{\odot}$ is best described by the black-hole-to-stellar mass ratio of the central galaxy, reflecting the AGN feedback relevance. We present a numerical model to populate dark matter-only simulations with HI at $0\leq z \leq 2$ based solely on halo parameters that are measurable in such simulations.

preprint2019arXiv

Fade to grey: systematic variation of the galaxy attenuation curves with galaxy properties in EAGLE

We present a simple model for galaxy attenuation by distilling SKIRT radiative transfer calculations for ~100,000 EAGLE galaxies at redshifts z=2-0. Our model adapts the two component screen model of Charlot & Fall (2000), parametrising the optical depth and slope of the ISM screen using the average dust surface density, $Σ_{\rm dust}$. We recover relatively tight relations between these parameters for the EAGLE sample, but also provide the scatter in these parameter owing to the morphological variation and orientation of galaxies. We also find that these relations are nearly independent of redshift in the EAGLE model. By pairing our model with an empirical prescription for birth clouds below the resolution scale of the simulation, we reproduce the observed relation between attenuation slope and optical depth for the first time in a cosmological simulation. We demonstrate that this result is remarkably independent of the attenuation properties assumed for birth cloud screen, merely requiring a boosted attenuation for infant stars. We present this model with a view to interpreting observations, as well as processing semi-analytic models and other hydrodynamic simulations.

preprint2019arXiv

The SAMI Galaxy Survey: The contribution of different kinematic classes to the stellar mass function of nearby galaxies

We use the complete Sydney-AAO Multi-object Integral field spectrograph (SAMI) Galaxy Survey to determine the contribution of slow rotators, as well as different types of fast rotators, to the stellar mass function of galaxies in the local Universe. We use stellar kinematics not only to discriminate between fast and slow rotators, but also to distinguish between dynamically cold systems (i.e., consistent with intrinsic axis ratios$<0.3$) and systems including a prominent dispersion-supported bulge. We show that fast rotators account for more than $80\%$ of the stellar mass budget of nearby galaxies, confirming that their number density overwhelms that of slow rotators at almost all masses from $10^{9}$ to $10^{11.5}{\rm M_\odot}$. Most importantly, dynamically cold disks contribute to at least $25\%$ of the stellar mass budget of the local Universe, significantly higher than what is estimated from visual morphology alone. For stellar masses up to $10^{10.5}{\rm M_\odot}$, this class makes up $>=30\%$ of the galaxy population in each stellar mass bin. The fact that many galaxies that are visually classified as having two-components have stellar spin consistent with dynamically cold disks suggests that the inner component is either rotationally-dominated (e.g., bar, pseudo-bulge) or has little effect on the global stellar kinematics of galaxies.