We have used the Australian Square Kilometre Array Pathfinder (ASKAP) radio telescope to search for intervening 21 cm neutral hydrogen (HI) absorption along the line of sight to 53 bright radio continuum sources. Our observations are sensitive to HI column densities typical of Damped Lyman Alpha absorbers (DLAs) in cool gas with an HI spin temperature below about 300-500 K. The six-dish Boolardy Engineering Test Array (BETA) and twelve-antenna Early Science array (ASKAP-12) covered a frequency range corresponding to redshift $0.4<z<1.0$ and $0.37<z<0.77$ respectively for the HI line. Fifty of the 53 radio sources observed have reliable optical redshifts, giving a total redshift path $\Delta z$ = 21.37. This was a spectroscopically-untargeted survey, with no prior assumptions about the location of the lines in redshift space. Four intervening HI lines were detected, two of them new. In each case, the estimated HI column density lies above the DLA limit for HI spin temperatures above 50-80 K, and we estimate a DLA number density at redshift $z\sim0.6$ of $n(z)=0.19\substack{+0.15 \\ -0.09}$. This value lies somewhat above the general trend of $n(z)$ with redshift seen in optical DLA studies. Although the current sample is small, it represents an important proof of concept for the much larger 21cm First Large Absorption Survey in HI (FLASH) project to be carried out with the full 36-antenna ASKAP telescope, probing a total redshift path $\Delta z\sim\,50,000$.

The Widefield ASKAP L-band Legacy All-sky Blind surveY (WALLABY) is a next-generation survey of neutral hydrogen (HI) in the Local Universe. It uses the widefield, high-resolution capability of the Australian Square Kilometer Array Pathfinder (ASKAP), a radio interferometer consisting of 36 x 12-m dishes equipped with Phased-Array Feeds (PAFs), located in an extremely radio-quiet zone in Western Australia. WALLABY aims to survey three-quarters of the sky (-90 degr < Dec < +30 degr) to a redshift of z < 0.26, and generate spectral line image cubes at ~30 arcsec resolution and ~1.6 mJy/beam per 4 km/s channel sensitivity. ASKAP's instantaneous field of view at 1.4 GHz, delivered by the PAF's 36 beams, is about 30 sq deg. At an integrated signal-to-noise ratio of five, WALLABY is expected to detect over half a million galaxies with a mean redshift of z ~ 0.05 (~200 Mpc). The scientific goals of WALLABY include: (a) a census of gas-rich galaxies in the vicinity of the Local Group; (b) a study of the HI properties of galaxies, groups and clusters, in particular the influence of the environment on galaxy evolution; and (c) the refinement of cosmological parameters using the spatial and redshift distribution of low-bias gas-rich galaxies. For context we provide an overview of previous large-scale HI surveys. Combined with existing and new multi-wavelength sky surveys, WALLABY will enable an exciting new generation of panchromatic studies of the Local Universe. - First results from the WALLABY pilot survey are revealed, with initial data products publicly available in the CSIRO ASKAP Science Data Archive (CASDA).

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.

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.

Observers experience a series of limitations when measuring galaxy kinematics, such as variable seeing conditions and aperture size. These effects can be reduced using empirical corrections, but these equations are usually applicable within a restrictive set of boundary conditions (e.g. S\'ersic indices within a given range) which can lead to biases when trying to compare measurements made across a full kinematic survey. In this work, we present new corrections for two widely used kinematic parameters, $\lambda_R$ and $V/\sigma$, that are applicable across a broad range of galaxy shapes, measurement radii and ellipticities. We take a series of mock observations of N-body galaxy models and use these to quantify the relationship between the observed kinematic parameters, structural properties and different seeing conditions. Derived corrections are then tested using the full catalogue of galaxies, including hydro-dynamic models from the EAGLE simulation. Our correction is most effective for regularly-rotating systems, yet the kinematic parameters of all galaxies -- fast, slow and irregularly rotating systems -- are recovered successfully. We find that $\lambda_R$ is more easily corrected than $V/\sigma$, with relative deviations of 0.02 and 0.06 dex respectively. The relationship between $\lambda_R$ and $V/\sigma$, as described by the parameter $\kappa$, also has a minor dependence on seeing conditions. These corrections will be particularly useful for stellar kinematic measurements in current and future integral field spectroscopic (IFS) surveys of galaxies.

We use the EAGLE suite of hydrodynamical simulations to analyse accretion rates (and the breakdown of their constituent channels) onto haloes over cosmic time, comparing the behaviour of baryons and dark matter (DM). We also investigate the influence of sub-grid baryon physics on halo-scale inflow, specifically the consequences of modelling radiative cooling, as well as feedback from stars and active galactic nuclei (AGN). We find that variations in halo baryon fractions at fixed mass (particularly their circum-galactic medium gas content) are very well correlated with variations in the baryon fraction of accreting matter, which we show to be heavily suppressed by stellar feedback in low-mass haloes, $M_{\rm halo}\lesssim10^{11.5}M_{\odot}$. Breaking down accretion rates into first infall, recycled, transfer and merger components, we show that baryons are much more likely to be smoothly accreted than to have originated from mergers when compared to DM, finding (averaged across halo mass) a merger contribution of $\approx6\%$ for baryons, and $\approx15\%$ for DM at $z\approx0$. We also show that the breakdown of inflow into different channels is strongly dependent on sub-grid physics, particularly the contribution of recycled accretion (accreting matter that has been previously ejected from progenitor haloes). Our findings highlight the dual role that baryonic feedback plays in regulating the evolution of galaxies and haloes: by (i) directly removing gas from haloes, and (ii) suppressing gas inflow to haloes.

We study the effect of the gas accretion rate ($\dot M_{\rm accr}$) on the radial gas metallicity profile (RMP) of galaxies using the EAGLE cosmological hydrodynamic simulations, focusing on central galaxies of stellar mass $M_\star \gtrsim 10^9 \, {\rm M_\odot}$ at $z \le 1$. We find clear relations between $\dot M_{\rm accr}$ and the slope of the RMP (measured within an effective radius), where higher $\dot M_{\rm accr}$ are associated with more negative slopes. The slope of the RMPs depends more strongly on $\dot M_{\rm accr}$ than on stellar mass, star formation rate or gas fraction, suggesting $\dot M_{\rm accr}$ to be a more fundamental driver of the RMP slope of galaxies. We find that eliminating the dependence on stellar mass is essential for pinning down the properties that shape the slope of the RMP. Although $\dot M_{\rm accr}$ is the main property modulating the slope of the RMP, we find that it causes other correlations that are more easily testable observationally: at fixed stellar mass, galaxies with more negative RMP slopes tend to have higher gas fractions and SFRs, while galaxies with lower gas fractions and SFRs tend to have flatter metallicity profiles within an effective radius.

We apply the spectral energy distribution-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. This analysis suggests that half of the mass in present-day elliptical galaxies was in place 11 Gyr ago, whereas 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 10.5 Gyr ago, whereas the least massive galaxies formed half their stellar mass 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 observations of metal densities in the Universe.

In the local (redshift z~0) Universe, collisional ring galaxies make up only ~0.01% of galaxies and are formed by head-on galactic collisions that trigger radially propagating density waves. These striking systems provide key snapshots for dissecting galactic disks and are studied extensively in the local Universe. However, not much is known about distant (z>0.1) collisional rings. Here we present a detailed study of a ring galaxy at a look-back time of 10.8 Gyr (z=2.19). Compared with our Milky Way, this galaxy has a similar stellar mass, but has a stellar half-light radius that is 1.5-2.2 times larger and is forming stars 50 times faster. The large, diffuse stellar light outside the star-forming ring, combined with a radial velocity on the ring and an intruder galaxy nearby, provides evidence for this galaxy hosting a collisional ring. If the ring is secularly evolved, the implied large bar in a giant disk would be inconsistent with the current understanding of the earliest formation of barred spirals. Contrary to previous predictions, this work suggests that massive collisional rings were as rare 11 Gyr ago as they are today. Our discovery offers a unique pathway for studying density waves in young galaxies, as well as constraining the cosmic evolution of spiral disks and galaxy groups.

Measurements of the cosmic far-infrared background (CIB) indicate that emission from many extragalactic phenomena, including star formation and black hole accretion, in the Universe can be obscured by dust. Resolving the CIB to study the population of galaxies in which this activity takes place is a major goal of submillimetre astronomy. Here, we present interferometric 650$\mu$m submillimetre number counts. Using the Band 8 data from the ALMACAL survey, we have analysed 81 ALMA calibrator fields together covering a total area of 5.5~arcmin$^2$. The typical central rms in these fields is $\sim 100 \mu$Jy~beam$^{-1}$ with the deepest maps reaching $\sigma = 47 \mu$Jy~beam$^{-1}$ at sub-arcsec resolution. Multi-wavelength coverage from ALMACAL allows us to exclude contamination from jets associated with the calibrators. However, residual contamination by jets and lensing remain a possibility. Using a signal-to-noise threshold of $4.5\sigma$, we find 21 dusty, star-forming galaxies with 650$\mu$m flux densities of $\geq 0.7 $mJy. At the detection limit we resolve $\simeq 100$ per cent of the CIB at 650$\mu$m, a significant improvement compared to low resolution studies at similar wavelength. We have therefore identified all the sources contributing to the EBL at 650 microns and predict that the contribution from objects with flux 0.7<mJy will be small.

We present early science results from the First Large Absorption Survey in HI (FLASH), a spectroscopically blind survey for 21-cm absorption lines in cold hydrogen HI gas at cosmological distances using the Australian Square Kilometre Array Pathfinder (ASKAP). We have searched for HI absorption towards 1253 radio sources in the GAMA 23 field, covering redshifts between $z = 0.34$ and $0.79$ over a sky area of approximately 50 deg$^{2}$. In a purely blind search we did not obtain any detections of 21-cm absorbers above our reliability threshold. Assuming a fiducial value for the HI spin temperature of $T_{\rm spin}$ = 100 K and source covering fraction $c_{\rm f} = 1$, the total comoving absorption path length sensitive to all Damped Lyman $\alpha$ Absorbers (DLAs; $N_{\rm HI} \geq 2 \times 10^{20}$ cm$^{-2}$) is $\Delta{X} = 6.6 \pm 0.3$ ($\Delta{z} = 3.7 \pm 0.2$) and super-DLAs ($N_{\rm HI} \geq 2 \times 10^{21}$ cm$^{-2}$) is $\Delta{X} = 111 \pm 6$ ($\Delta{z} = 63 \pm 3$). We estimate upper limits on the HI column density frequency distribution function that are consistent with measurements from prior surveys for redshifted optical DLAs, and nearby 21-cm emission and absorption. By cross matching our sample of radio sources with optical spectroscopic identifications of galaxies in the GAMA 23 field, we were able to detect 21-cm absorption at $z = 0.3562$ towards NVSS J224500$-$343030, with a column density of $N_{\rm HI} = (1.2 \pm 0.1) \times 10^{20} (T_{\rm spin}/100~\mathrm{K})$ cm$^{-2}$. The absorber is associated with GAMA J22450.05$-$343031.7, a massive early-type galaxy at an impact parameter of 17 kpc with respect to the radio source and which may contain a massive ($M_{\rm HI} \gtrsim 3 \times 10^{9}$ M$_{\odot}$) gas disc. Such gas-rich early types are rare, but have been detected in the nearby Universe.

We introduce ProSpect, a generative galaxy spectral energy distribution (SED) package that encapsulates the best practices for SED methodologies in a number of astrophysical domains. ProSpect comes with two popular families of stellar population libraries (BC03 and EMILES), and a large variety of methods to construct star formation and metallicity histories. It models dust through the use of a Charlot & Fall attenuation model, with re-emission using Dale far-infrared templates. It also has the ability to model AGN through the inclusion of a simple AGN and hot torus model. Finally, it makes use of MAPPINGS-III photoionisation tables to produce line emission features. We test the generative and inversion utility of ProSpect through application to the Shark galaxy formation semi-analytic code, and informed by these results produce fits to the final ultraviolet to far-infrared photometric catalogues produces by the Galaxy and Mass Assembly Survey (GAMA). As part of the testing of ProSpect, we also produce a range of simple photometric stellar mass approximations covering a range of filters for both observed frame and rest frame photometry.

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's fundamental (stellar mass, angular momentum) properties. In future, PDF($s$) may also be useful as a kinematic decomposition tool.

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.

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\sigma/(GM)$, where $M$ and $j$ are the baryonic mass and specific angular momentum of the disk and $\sigma$ is the HI velocity dispersion Obreschkow et al. 2016. The log-ratio $\Delta f_q$ between this predictor and the observed atomic fraction can be seen as a physically motivated `HI deficiency'. While field disk galaxies have $\Delta f_q \approx0$, objects subject to environmental removal of HI are expected to have $\Delta 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 $\Delta f_q$ slightly increases when getting closer to the cluster centres. The $\Delta f_q$ values are similar to traditional HI deficiency estimators, but $\Delta 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 $\Delta 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 $\Delta f_q$ and quenching of star formation.

We present an analysis of morphological, kinematic and spectral asymmetries in observations of atomic neutral hydrogen (HI) gas from the Local Volume HI Survey (LVHIS), the VLA Imaging of Virgo in Atomic Gas (VIVA) survey and the Hydrogen Accretion in Local Galaxies Survey (HALOGAS). With the aim of investigating the impact of the local environment density and stellar mass on the measured HI asymmetries in future large HI surveys, we provide recommendations for the most meaningful measures of asymmetry for use in future analysis. After controlling for stellar mass, we find signs of statistically significant trends of increasing asymmetries with local density. The most significant trend we measure is for the normalised flipped spectrum residual ($A_{\mathrm{spec}}$), with mean LVHIS and VIVA values of $0.204\pm0.011$ and $0.615\pm0.068$ at average weighted $10^{\mathrm{th}}$ nearest-neighbour galaxy number densities of $\log(\rho_{10}/\mathrm{Mpc}^{-3})=-1.64$ and 0.88, respectively. Looking ahead to the WALLABY survey on the Australian Square Kilometre Array Pathfinder (ASKAP), we estimate that the number of detections will be sufficient to provide coverage over 5 orders of magnitude in both local density and stellar mass increasing the dynamic range and accuracy with which we can probe the effect of these properties on the asymmetry in the distribution of atomic gas in galaxies.

The gas content of galaxies is a key factor for their growth, starting from star formation and black hole accretion to galaxy mergers. Thus, characterising its properties via observations of tracers like the CO emission line is of big importance in order to understand the bigger picture of galaxy evolution. We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of dust continuum, CO(5-4) and CO(8-7) line emission in the quasar--star-forming companion system SMM J04135+10277 (z=2.84). Earlier low-$J$ CO studies of this system found a huge molecular gas reservoir associated to the companion galaxy, while the quasar appeared gas-poor. Our CO observations revealed that the host galaxy of the quasar is also gas-rich, with an estimated molecular gas mass of $\sim(0.7-2.3)\times10^{10}$ M$_{\odot}$. The CO line profiles of the companion galaxy are broad ($\sim1000$ km s$^{-1}$), and show signs of rotation of a compact, massive system. In contrast to previous far-infrared observations, we resolve the continuum emission and detect both sources, with the companion galaxy dominating the dust continuum and the quasar having a $\sim25\%$ contribution to the total dust emission. By fitting the infrared spectral energy distribution of the sources with \textsc{MR-MOOSE} and empirical templates, the infrared luminosities of the quasar and the companion are in the range of $L_{\rm IR, QSO}\sim(2.1-9.6)\times10^{12}$ L$_{\odot}$ and $L_{\rm IR, Comp.}\sim(2.4-24)\times10^{12}$ L$_{\odot}$, while the estimated star formation rates are $\sim210-960$ M$_{\odot}$ yr$^{-1}$ and $\sim240-2400$ M$_{\odot}$ yr$^{-1}$, respectively. Our results demonstrate that non-detection of low-$J$ CO transition lines in similar sources does not necessarily imply the absence of massive molecular gas reservoir but that the excitation conditions favour the excitation of high-$J$ transitions.

The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope Treasury program taking data from 2018-2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in WFC3/IR F105W, F125W, and F160W and ACS/WFC F606W and F814W around each of the six HFF clusters and flanking fields. This additional area has not been observed by HST but is already covered by deep multi-wavelength datasets, including Spitzer and Chandra. As with the original HFF program, BUFFALO is designed to take advantage of gravitational lensing from massive clusters to simultaneously find high-redshift galaxies which would otherwise lie below HST detection limits and model foreground clusters to study properties of dark matter and galaxy assembly. The expanded area will provide a first opportunity to study both cosmic variance at high redshift and galaxy assembly in the outskirts of the large HFF clusters. Five additional orbits are reserved for transient followup. BUFFALO data including mosaics, value-added catalogs and cluster mass distribution models will be released via MAST on a regular basis, as the observations and analysis are completed for the six individual clusters.

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 "tree entropy", to parametrise the geometry of a halo's entire mass assembly hierarchy, building on a generalisation of Shannon'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 $\delta_{\rm c}/\sigma(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.

We use the EAGLE simulations to study the effects of the intra-cluster medium (ICM) on the spatially resolved star-formation activity in galaxies. We study three cases of galaxy asymmetry dividing each galaxy in two halves using the plane (i) perpendicular to the $\rm \texttt{velocity}$ direction, differentiating the galaxy part approaching to the cluster center, hereafter dubbed as the "leading half", and the opposite one "trailing half", (ii) perpendicular to the $\rm \texttt{radial}$ position of the satellite to the centre of the cluster, (iii) that maximizes the star-formation rate ($\rm SFR$) difference between the two halves. For (i), we find an enhancement of the $\rm SFR$, star formation efficiency ($\rm SFE$), and interstellar medium pressure in the leading half with respect to the trailing one and normal star-forming galaxies in the EAGLE simulation, and a clear overabundance of gas particles in their trailing. These results suggest that ram pressure (RP) is boosting the star formation by gas compression in the leading half, and transporting the gas to the trailing half. This effect is more pronounced in satellites of intermediate stellar masses $\rm 10^{9.5-10.5} M_{\odot}$, with gas masses above $\rm 10^{9} M_{\odot}$, and located within one virial radius or in the most massive clusters. In (iii) we find an alignment between the velocity and the vector perpendicular to the plane that maximizes the $\rm SFR$ difference between the two halves. It suggests that finding this plane in real galaxies can provide insights into the velocity direction.

We use spectral stacking to measure the contribution of galaxies of different masses and in different hierarchies to the cosmic atomic hydrogen (HI) mass density in the local Universe. Our sample includes 1793 galaxies at $z < 0.11$ observed with the Westerbork Synthesis Radio Telescope, for which Sloan Digital Sky Survey spectroscopy and hierarchy information are also available. We find a cosmic HI mass density of $\Omega_{\rm HI} = (3.99 \pm 0.54)\times 10^{-4} h_{70}^{-1}$ at $\langle z\rangle = 0.065$. For the central and satellite galaxies, we obtain $\Omega_{\rm HI}$ of $(3.51 \pm 0.49)\times 10^{-4} h_{70}^{-1}$ and $(0.90 \pm 0.16)\times 10^{-4} h_{70}^{-1}$, respectively. We show that galaxies above and below stellar masses of $\sim$10$^{9.3}$ M$_{\odot}$ contribute in roughly equal measure to the global value of $\Omega_{\rm HI}$. While consistent with estimates based on targeted HI surveys, our results are in tension with previous theoretical work. We show that these differences are, at least partly, due to the empirical recipe used to set the partition between atomic and molecular hydrogen in semi-analytical models. Moreover, comparing our measurements with the cosmological semi-analytic models of galaxy formation {\sc Shark} and GALFORM reveals gradual stripping of gas via ram pressure works better to fully reproduce the properties of satellite galaxies in our sample, than strangulation. Our findings highlight the power of this approach in constraining theoretical models, and confirm the non-negligible contribution of massive galaxies to the HI mass budget of the local Universe.

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$\alpha$, H$\beta$, 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\sigma$ 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.

There is a severe tension between the observed star formation rate (SFR) - stellar mass (${\rm M}_{\star}$) relations reported by different authors at $z = 1-4$. In addition, the observations have not been successfully reproduced by state-of-the-art cosmological simulations which tend to predict a factor of 2-4 smaller SFRs at a fixed ${\rm M}_{\star}$. We examine the evolution of the SFR$-{\rm M}_{\star}$ relation of $z = 1-4 $ galaxies using the SKIRT simulated spectral energy distributions of galaxies sampled from the EAGLE simulations. We derive SFRs and stellar masses by mimicking different observational techniques. We find that the tension between observed and simulated SFR$-{\rm M}_{\star}$ relations is largely alleviated if similar methods are used to infer the galaxy properties. We find that relations relying on infrared wavelengths (e.g. 24 ${\rm \mu m}$, MIPS - 24, 70 and 160 ${\rm \mu m}$ or SPIRE - 250, 350, 500 ${\rm \mu m}$) have SFRs that exceed the intrinsic relation by 0.5 dex. Relations that rely on the spectral energy distribution fitting technique underpredict the SFRs at a fixed stellar mass by -0.5 dex at $z \sim 4$ but overpredict the measurements by 0.3 dex at $z \sim 1$. Relations relying on dust-corrected rest-frame UV luminosities, are flatter since they overpredict/underpredict SFRs for low/high star forming objects and yield deviations from the intrinsic relation from 0.10 dex to -0.13 dex at $z \sim 4$. We suggest that the severe tension between different observational studies can be broadly explained by the fact that different groups employ different techniques to infer their SFRs.

We present predictions for the evolution of radio emission from Active Galactic Nuclei (AGNs). We use a model that follows the evolution of Supermassive Black Hole (SMBH) masses and spins, within the latest version of the GALFORM semi-analytic model of galaxy formation. We use a Blandford-Znajek type model to calculate the power of the relativistic jets produced by black hole accretion discs, and a scaling model to calculate radio luminosities. First, we present the predicted evolution of the jet power distribution, finding that this is dominated by objects fuelled by hot halo accretion and an ADAF accretion state for jet powers above $10^{32}\mathrm{W}$ at $z=0$, with the contribution from objects fuelled by starbursts and in a thin disc accretion state being more important for lower jet powers at $z=0$ and at all jet powers at high redshifts ($z\geq3$). We then present the evolution of the jet power density from the model. The model is consistent with current observational estimates of jet powers from radio luminosities, once we allow for the significant uncertainties in these observational estimates. Next, we calibrate the model for radio emission to a range of observational estimates of the $z=0$ radio luminosity function. We compare the evolution of the model radio luminosity function to observational estimates for $0<z<6$, finding that the predicted evolution is similar to that observed. Finally, we explore recalibrating the model to reproduce luminosity functions of core radio emission, finding that the model is in approximate agreement with the observations.

Telescopes to be launched over the next decade-and-a-half, such as JWST, EUCLID, ATHENA and Lynx, promise to revolutionise the study of the high redshift Universe and greatly advance our understanding of the early stages of galaxy formation. We use a model that follows the evolution of the masses and spins of supermassive black holes (SMBHs) within a semi-analytic model of galaxy formation to make predictions for the Active Galactic Nucleus (AGN) luminosity function at $z\geq7$ in the broadband filters of JWST and EUCLID at near-infrared wavelengths, and ATHENA and Lynx at X-ray energies. The predictions of our model are relatively insensitive to the choice of seed black hole mass, except at the lowest luminosities ($L_{\mathrm{bol}}<10^{43}\mathrm{ergs^{-1}}$) and the highest redshifts ($z>10$). We predict that surveys with these different telescopes will select somewhat different samples of SMBHs, with EUCLID unveiling the most massive, highest accretion rate SMBHs, Lynx the least massive, lowest accretion rate SMBHs, and JWST and ATHENA covering objects inbetween. At $z=7$, we predict that typical detectable SMBHs will have masses, $M_{\mathrm{BH}}\sim10^{5-8}M_{\odot}$, and Eddington normalised mass accretion rates, $\dot{M}/\dot{M}_{\mathrm{Edd}}\sim0.6-2$. The SMBHs will be hosted by galaxies of stellar mass $M_{\star}\sim10^{8-10}M_{\odot}$, and dark matter haloes of mass $M_{\mathrm{halo}}\sim10^{11-12}M_{\odot}$. We predict that the detectable SMBHs at $z=10$ will have slightly smaller black holes, accreting at slightly higher Eddington normalised mass accretion rates, in slightly lower mass host galaxies compared to those at $z=7$, and reside in haloes of mass $M_{\mathrm{halo}}\sim10^{10-11}M_{\odot}$.

We use photometric redshifts derived from new $u$-band through 4.5$\mu$m Spitzer IRAC photometry in the 4.8\,deg$^2$ of the XMM-LSS field to construct surface density maps in the redshift range 0.1-1.5. Our density maps show evidence for large-scale structure in the form of filaments spanning several tens of Mpc. Using these maps, we identify 339 overdensities that our simulated lightcone analysis suggests are likely associated with dark matter haloes with masses, $M_{\rm halo}$, log($M_{\rm halo}/M_{\odot})>$13.7. From this list of overdensities we recover 43 of 70 known X-ray detected and spectroscopically confirmed clusters. The missing X-ray clusters are largely at lower redshifts and lower masses than our target log($M_{\rm halo}/M_{\odot})>$13.7. The bulk of the overdensities are compact, but a quarter show extended morphologies which include likely projection effects, clusters embedded in apparent filaments as well as at least one potential cluster merger (at $z\sim1.28$). The strongest overdensity in our highest redshift slice (at $z\sim1.5$) shows a compact red galaxy core potentially implying a massive evolved cluster.

We present a detailed analysis of a large sample of spectroscopically confirmed ultra-massive quiescent galaxies (${\rm{log}}(M_{\ast}/M_{\odot})\sim11.5$) at $z\gtrsim2$. This sample comprises 15 galaxies selected in the COSMOS and UDS fields by their bright K-band magnitudes and followed up with VLT/X-shooter spectroscopy and HST/WFC3 $H_{F160W}$ imaging. These observations allow us to unambiguously confirm their redshifts ascertain their quiescent nature and stellar ages, and to reliably assess their internal kinematics and effective radii. We find that these galaxies are compact, consistent with the high mass end of the mass-size relation for quiescent galaxies at $z=2$. Moreover, the distribution of the measured stellar velocity dispersions of the sample is consistent with the most massive local early-type galaxies from the MASSIVE Survey showing that evolution in these galaxies, is dominated by changes in size. The HST images reveal, as surprisingly high, that $40\ \%$ of the sample have tidal features suggestive of mergers and companions in close proximity, including three galaxies experiencing ongoing major mergers. The absence of velocity dispersion evolution from $z=2$ to $0$, coupled with a doubling of the stellar mass, with a factor of four size increase and the observed disturbed stellar morphologies support dry minor mergers as the primary drivers of the evolution of the massive quiescent galaxies over the last 10 billion years.

We study the behaviour of the spin-ellipticity radial tracks for 507 galaxies from the Sydney AAO Multi-object Integral Field (SAMI) Galaxy Survey with stellar kinematics out to $\geq1.5R_\text{e}$. We advocate for a morpho-dynamical classification of galaxies, relying on spatially-resolved photometric and kinematic data. We find the use of spin-ellipticity radial tracks is valuable in identifying substructures within a galaxy, including embedded and counter-rotating discs, that are easily missed in unilateral studies of the photometry alone. Conversely, bars are rarely apparent in the stellar kinematics but are readily identified on images. Consequently, we distinguish the spin-ellipticity radial tracks of seven morpho-dynamical types: elliptical, lenticular, early spiral, late spiral, barred spiral, embedded disc, and 2-sigma galaxies. The importance of probing beyond the inner radii of galaxies is highlighted by the characteristics of galactic features in the spin-ellipticity radial tracks present at larger radii. The density of information presented through spin-ellipticity radial tracks emphasises a clear advantage to representing galaxies as a track, rather than a single point, in spin-ellipticity parameter space.

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.

We present the first detection of mass dependent galactic spin alignments with local cosmic filaments with over 2 sigma confidence using IFS kinematics. The 3D network of cosmic filaments is reconstructed on Mpc scales across GAMA fields using the cosmic web extractor DisPerSe. We assign field galaxies from the SAMI survey to their nearest filament segment in 3D and estimate the degree of alignment between SAMI galaxies kinematic spin axis and their nearest filament in projection. Low-mass galaxies align their spin with their nearest filament while higher mass counterparts are more likely to display an orthogonal orientation. The stellar transition mass from the first trend to the second is bracketed between log stellar masses 10.4 and 10.9, with hints of an increase with filament scale. Consistent signals are found in the HorizonAGN cosmological hydrodynamic simulation. This supports a scenario of early angular momentum build-up in vorticity rich quadrants around filaments at low stellar mass followed by progressive flip of spins orthogonal to the cosmic filaments through mergers at high stellar mass. Conversely, we show that dark-matter only simulations post-processed with a semi-analytic model treatment of galaxy formation struggles to reproduce this alignment signal. This suggests that gas physics is key in enhancing the galaxy-filament alignment.

We combine the Shark semi-analytic model of galaxy formation with the ProSpect software tool for spectral energy distribution (SED) generation to study the multi-wavelength emission of galaxies from the far-ultraviolet (FUV) to the far-infrared (FIR) at $0\le z\le 10$. We produce a physical model for the attenuation of galaxies across cosmic time by combining a local Universe empirical relation to compute the dust mass of galaxies from their gas metallicity and mass, attenuation curves derived from radiative transfer calculations of galaxies in the EAGLE hydrodynamic simulation suite, and the properties of Shark galaxies. We are able to produce a wide range of galaxies, from the $z=8$ star-forming galaxies with almost no extinction, $z=2$ submillimeter galaxies, down to the normal star-forming and red sequence galaxies at $z=0$. Quantitatively, we find that Shark reproduces the observed (i) the $z=0$ FUV-to-FIR, (ii) $0\le z\le 3$ rest-frame $K$-band, and (iii) $0\le z\le 10$ rest-frame FUV luminosity functions, (iv) $z\le 8$ UV slopes, (v) the FUV-to-FIR number counts (including the widely disputed 850$\mu$m), (vi) redshift distribution of bright $850\mu$m galaxies and (vii) the integrated cosmic SED from $z=0$ to $z=1$ to an unprecedented level. This is achieved without the need to invoke changes in the stellar initial mass function, dust-to-metal mass ratio, or metal enrichment timescales. Our model predicts star formation in galaxy disks to dominate in the FUV-to-optical, while bulges dominate at the NIR at all redshifts. The FIR sees a strong evolution in which disks dominate at $z\le 1$ and starbursts (triggered by both galaxy mergers and disk instabilities, in an even mix) dominate at higher redshifts, even out to $z=10$.

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, $\Sigma_{\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.

We study the Intra-Halo Stellar Component (IHSC) of Milky Way-mass systems up to galaxy clusters in the Horizon-AGN cosmological hydrodynamical simulation. We identify the IHSC using an improved phase-space galaxy finder algorithm which provides an adaptive, physically motivated and shape-independent definition of this stellar component, that can be applied to halos of arbitrary masses. We explore the IHSC mass fraction-total halo's stellar mass, $f_{M*,IHSC}-M*$, relation and the physical drivers of its scatter. We find that on average the $f_{M*,IHSC}$ increases with $M_{*,tot}$, with the scatter decreasing strongly with mass from 2 dex at $M_{*,tot}\sim10^{11}M_\odot$ to 0.3 dex at group masses. At high masses, $M_{*,tot}>10^{11.5}M_\odot$, $f_{M*,IHSC}$ increases with the number of substructures, and with the mass ratio between the central galaxy and largest satellite, at fixed $M_{*,tot}$. From mid-size groups and systems below $M_{*,tot}<10^{12}M_\odot$, we find that the central galaxy's stellar rotation-to-dispersion velocity ratio, V/{\sigma}, displays the strongest (anti)-correlation with $f_{M*,IHSC}$ at fixed $M_{*,tot}$ of all the galaxy and halo properties explored, transitioning from $f_{M*,IHSC}$<0.1% for high V/{\sigma}, to $f_{M*,IHSC}\sim5$% for low V/{\sigma} galaxies. By studying the $f_{M*,IHSC}$ temporal evolution, we find that, in the former, mergers not always take place, but if they did, they happened early (z>1), while the high $f_{M*,IHSC}$ population displays a much more active merger history. In the case of massive groups and galaxy clusters, $M_{*,tot}>10^{12}M_\odot$, a fraction $f_{M*,IHSC}\sim$10-20% is reached at $z\sim1$ and then they evolve across lines of constant $f_{M*,IHSC}$ modulo some small perturbations. Because of the limited simulation's volume, the latter is only tentative and requires a larger sample of simulated galaxy clusters to confirm.

Current observational results show that both late-and-early-type galaxies follow tight mass-size planes, on which physical properties such as age, velocity dispersion and metallicities correlate with the scatter on the plane. We study the mass-size plane of galaxies in cosmological hydrodynamical simulations, as a function of velocity dispersion, age, chemical abundances, ellipticity and spin parameters with the aim at assessing to what extent the current cosmological paradigm can reproduce these observations and provide a physical interpretation of them. We select a sample of well-resolved galaxies from the (100 Mpc)^3 simulation of the EAGLE Project. This sample is composed by 508 spheroid-dominated galaxies and 1213 disc-dominated galaxies. The distributions of velocity dispersion, age, metallicity indicators and gradients and spin parameters across the mass-size plane are analysed. Furthermore, we study the relation between shape and kinematic parameters. The results are compared with observations. The mass-weighted ages of the EAGLE galaxies are found to vary along lines of constant velocity dispersion on the mass-size plane, except for galaxies with velocity dispersion larger than aprox 150 km s^(-1) . Negative age gradients tend to be found in extended disc galaxies in agreement with observations. However, the age distributions of early-type galaxies show a larger fraction with inverted radial profiles. The distribution of metallicity gradients does not show any clear dependence on this plane. Galaxies with similar spin parameters ({\lambda}) display larger sizes as their dynamical masses increase. Stellar-weighted ages are found to be good proxies for {\lambda} in galaxies with low ellipticity ({\epsilon}). Abridged

Using a spectral stacking technique, we measure the neutral hydrogen (HI) properties of a sample of galaxies at $z < 0.11$ across 35 pointings of the Westerbork Synthesis Radio Telescope (WSRT). The radio data contains 1,895 galaxies with redshifts and positions known from the Sloan Digital Sky Survey (SDSS). We carefully quantified the effects of sample bias, aperture used to extract spectra, sidelobes and weighting technique and use our data to provide a new estimate for the cosmic HI mass density. We find a cosmic HI mass density of $\Omega_{\rm HI} = (4.02 \pm 0.26)\times 10^{-4} h_{70}^{-1}$ at $\langle z\rangle = 0.066$, consistent with measurements from blind HI surveys and other HI stacking experiments at low redshifts. The combination of the small interferometer beam size and the large survey volume makes our result highly robust against systematic effects due to confusion at small scales and cosmic variance at large scales. Splitting into three sub-samples with $\langle z\rangle$ = 0.038, 0.067 and 0.093 shows no significant evolution of the HI gas content at low redshift.

Accurately predicting the shape of the HI velocity function of galaxies is regarded widely as a fundamental test of any viable dark matter model. Straightforward analyses of cosmological $N$-body simulations imply that the $\Lambda$CDM model predicts an overabundance of low circular velocity galaxies when compared to observed HI velocity functions. More nuanced analyses that account for the relationship between galaxies and their host haloes suggest that how we model the influence of baryonic processes has a significant impact on HI velocity function predictions. We explore this in detail by modelling HI emission lines of galaxies in the SHARK semi-analytic galaxy formation model, built on the SURFS suite of $\Lambda$CDM $N$-body simulations. We create a simulated ALFALFA survey, in which we apply the survey selection function and account for effects such as beam confusion, and compare simulated and observed HI velocity width distributions, finding differences of $\lesssim 50$%, orders of magnitude smaller than the discrepancies reported in the past. This is a direct consequence of our careful treatment of survey selection effects and, importantly, how we model the relationship between galaxy and halo circular velocity - the HI mass-maximum circular velocity relation of galaxies is characterised by a large scatter. These biases are complex enough that building a velocity function from the observed HI line widths cannot be done reliably.

Cosmological hydrodynamical simulations are rich tools to understand the build-up of stellar mass and angular momentum in galaxies, but require some level of calibration to observations. We compare predictions at $z\sim0$ from the Eagle, Hydrangea, Horizon-AGN, and Magneticum simulations with integral field spectroscopic (IFS) data from the SAMI Galaxy Survey, ATLAS3D, CALIFA and MASSIVE surveys. The main goal of this work is to simultaneously compare structural, dynamical, and stellar population measurements in order to identify key areas of success and tension. We have taken great care to ensure that our simulated measurement methods match the observational methods as closely as possible. We find that the Eagle and Hydrangea simulations reproduce many galaxy relations but with some offsets at high stellar masses. There are moderate mismatches in $R_e$ (+), $\epsilon$ (-), $\sigma_e$ (-), and mean stellar age (+), where a plus sign indicates that quantities are too high on average, and minus sign too low. The Horizon-AGN simulations qualitatively reproduce several galaxy relations, but there are a number of properties where we find a quantitative offset to observations. Massive galaxies are better matched to observations than galaxies at low and intermediate masses. Overall, we find mismatches in $R_e$ (+), $\epsilon$ (-), $\sigma_e$ (-) and $(V/\sigma)_e$ (-). Magneticum matches observations well: this is the only simulation where we find ellipticities typical for disk galaxies, but there are moderate differences in $\sigma_e$ (-), $(V/\sigma)_e$ (-) and mean stellar age (+). Our comparison between simulations and observational data has highlighted several areas for improvement, such as the need for improved modelling resulting in a better vertical disk structure, yet our results demonstrate the vast improvement of cosmological simulations in recent years.

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.

We present recalibrations of the GALFORM semi-analytical model of galaxy formation in a new N-body simulation with the Planck cosmology. The Planck Millennium simulation uses more than 128 billion particles to resolve the matter distribution in a cube of $800$ Mpc on a side, which contains more than 77 million dark matter haloes with mass greater than $2.12 \times 10^{9} h^{-1} {\rm M_{\odot}}$ at the present day. Only minor changes to a very small number of model parameters are required in the recalibration. We present predictions for the atomic hydrogen content (HI) of dark matter halos, which is a key input into the calculation of the HI intensity mapping signal expected from the large-scale structure of the Universe. We find that the HI mass $-$ halo mass relation displays a clear break at the halo mass above which AGN heating suppresses gas cooling, $\approx 3 \times 10^{11} h^{-1} M_{\rm \odot}$. Below this halo mass, the HI content of haloes is dominated by the central galaxy; above this mass it is the combined HI content of satellites that prevails. We find that the HI mass - halo mass relation changes little with redshift up to $z=3$. The bias of HI sources shows a scale dependence that gets more pronounced with increasing redshift.

Most bulge-dominated galaxies host black holes with masses that tightly correlate with the masses of their bulges. This may indicate that the black holes may regulate galaxy growth or vice versa, or that they may grow in lock-step. The quest to understand how, when, and where those black-holes formed motivates much of extragalactic astronomy. Here we focus on a population of galaxies with active black holes in their nuclei (active galactic nuclei or AGN), that are fully or partially hidden by dust and gas: the emission from the broad line region is either completely or partially obscured with a visual extinction of 1 or above. This limit, though not yet precise, appears to be the point at which the populations of AGN may evolve differently. We highlight the importance of finding and studying those dusty AGN at redshifts between 1 and 3, the epoch when the universe may have gone through its most dramatic changes. We emphasize the need for future large multiplexed spectroscopic instruments that can perform dedicated surveys in the optical and NIR to pin down the demographics of such objects and study their reddening properties, star-formation histories, and excitation conditions. These key studies will shed light on the role of black holes in galaxy evolution during the epoch of peak growth activity.

This paper reports on the atomic hydrogen gas (HI) observations of the spiral galaxy NGC 1566 using the newly commissioned Australian Square Kilometre Array Pathfinder (ASKAP) radio telescope. We measure an integrated HI flux density of $180.2$ Jy km s$^{-1}$ emanating from this galaxy, which translates to an HI mass of $1.94\times10^{10}$M$_\circ$ at an assumed distance of $21.3$ Mpc. Our observations show that NGC 1566 has an asymmetric and mildly warped HI disc. The HI-to-stellar mass fraction of NGC 1566 is $0.29$, which is high in comparison with galaxies that have the same stellar mass ($10^{10.8}$M$_\circ$). We also derive the rotation curve of this galaxy to a radius of $50$ kpc and fit different mass models to it. The NFW, Burkert and pseudo-isothermal dark matter halo profiles fit the observed rotation curve reasonably well and recover dark matter fractions of $0.62$, $0.58$ and $0.66$, respectively. Down to the column density sensitivity of our observations ($N_{HI} = 3.7\times10^{19}$ cm$^{-2}$), we detect no HI clouds connected to, or in the nearby vicinity of, the HI disc of NGC 1566 nor nearby interacting systems. We conclude that, based on a simple analytic model, ram pressure interactions with the IGM can affect the HI disc of NGC 1566 and is possibly the reason for the asymmetries seen in the HI morphology of NGC 1566.

We use the \eagle\ simulations to study the connection between the quenching timescale, $\tau_{\rm Q}$, and the physical mechanisms that transform star-forming galaxies into passive galaxies. By quantifying $\tau_{\rm Q}$ in two complementary ways - as the time over which (i) galaxies traverse the green valley on the colour-mass diagram, or (ii) leave the main sequence of star formation and subsequently arrive on the passive cloud in specific star formation rate (SSFR)-mass space - we find that the $\tau_{\rm Q}$ distribution of high-mass centrals, low-mass centrals and satellites are divergent. In the low stellar mass regime where $M_{\star}<10^{9.6}M_{\odot}$, centrals exhibit systematically longer quenching timescales than satellites ($\approx 4$~Gyr compared to $\approx 2$~Gyr). Satellites with low stellar mass relative to their halo mass cause this disparity, with ram pressure stripping quenching these galaxies rapidly. Low mass centrals are quenched as a result of stellar feedback, associated with long $\tau_{\rm Q}\gtrsim 3$~Gyr. At intermediate stellar masses where $10^{9.7}\,\rm M_{\odot}<M_{\star}<10^{10.3}\,\rm M_{\odot}$, $\tau_{\rm Q}$ are the longest for both centrals and satellites, particularly for galaxies with higher gas fractions. At $M_{\star}\gtrsim 10^{10.3}\,\rm M_{\odot}$, galaxy merger counts and black hole activity increase steeply for all galaxies. Quenching timescales for centrals and satellites decrease with stellar mass in this regime to $\tau_{\rm Q}\lesssim2$~Gyr. In anticipation of new intermediate redshift observational galaxy surveys, we analyse the passive and star-forming fractions of galaxies across redshift, and find that the $\tau_{\rm Q}$ peak at intermediate stellar masses is responsible for a peak (inflection point) in the fraction of green valley central (satellite) galaxies at $z\approx 0.5-0.7$.

We present predictions for the UV-to-mm extragalactic background light (EBL) from a recent version of the GALFORM semi-analytical model of galaxy formation which invokes a top-heavy stellar initial mass function (IMF) for galaxies undergoing dynamically-triggered bursts of star formation. We combine GALFORM with the GRASIL radiative transfer code for computing fully self-consistent UV-to-mm spectral energy distributions for each simulated galaxy, accounting for the absorption and re-emission of stellar radiation by interstellar dust. The predicted EBL is in near-perfect agreement with recent observations over the whole UV-to-mm spectrum, as is the evolution of the cosmic spectral energy distribution over the redshift range for which observations are available ($z\lesssim1$). We show that approximately 90~per~cent of the EBL is produced at $z<2$ although this shifts to higher redshifts for sub-mm wavelengths. We assess whether the top-heavy IMF in starbursts is necessary in order to reproduce the EBL at the same time as other key observables, and find that variant models with a universal solar-neighborhood IMF display poorer agreement with EBL observations over the whole UV-to-mm spectrum and fail to match the counts of galaxies in the sub-mm.

The growth channel of massive galaxies involving mergers can be studied via close pairs as putative merger progenitors, where the stellar populations of the satellite galaxies will be eventually incorporated into the massive primaries. We extend our recent analysis of the GAMA-based sample of close pairs presented in Ferreras et al. to the general spectroscopic dataset of SDSS galaxies (DR14), for which the high S/N of the data enables a detailed analysis of the differences between satellite galaxies with respect to the mass of the primary galaxy. A sample of approximately two thousand satellites of massive galaxies is carefully selected within a relatively narrow redshift range (0.07<z<0.14). Two main parameters are considered as major drivers of the star formation history of these galaxies, namely: the stellar velocity dispersion of the satellite ($\sigma$), as a proxy of "local" drivers, and the ratio between the stellar mass of the satellite and the primary, $\mu=M_{\rm SAT}/M_{\rm PRI}$, meant to serve as an indicator of environment. Consistently with the independent, GAMA-based work, we find that satellites around the most massive primaries appear older, at fixed velocity dispersion, with respect to satellites of lower mass primaries. This trend is more marked in lower mass satellites ($\sigma$~100 km/s), with SSP-equivalent age differences up to ~0.5 Gyr, and can be interpreted as a one-halo assembly bias, so that satellites corresponding to smaller values of the mass ratio $\mu$ represent older structures, akin to fossil groups.

We study the internal radial gradients of the stellar populations in a sample comprising 522 early-type galaxies (ETGs) from the SAMI (Sydney- AAO Multi-object Integral field spectrograph) Galaxy Survey. We stack the spectra of individual spaxels in radial bins, and derive basic stellar population properties: total metallicity ([Z/H]), [Mg/Fe], [C/Fe] and age. The radial gradient ($\nabla$) and central value of the fits (evaluated at R$_e$/4) are compared against a set of six possible drivers of the trends. We find that velocity dispersion ($\sigma$) - or, equivalently gravitational potential - is the dominant driver of the chemical composition gradients. Surface mass density is also correlated with the trends, especially with stellar age. The decrease of $\nabla$[Mg/Fe] with increasing $\sigma$ is contrasted by a rather shallow dependence of $\nabla$[Z/H] with $\sigma$ (although this radial gradient is overall rather steep). This result, along with a shallow age slope at the massive end, imposes stringent constraints on the progenitors of the populations that contribute to the formation of the outer envelopes of ETGs. The SAMI sample is split between a 'field' sample and a cluster sample. Only weak environment-related differences are found, most notably a stronger dependence of central total metallicity ([Z/H]$_{e4}$) with $\sigma$, along with a marginal trend of $\nabla$[Z/H] to steepen in cluster galaxies, a result that is not followed by [Mg/Fe]. The results presented here serve as constraints on numerical models of the formation and evolution of ETGs.

We present TreeFrog, a massively parallel halo merger tree builder that is capable comparing different halo catalogues and producing halo merger trees. The code is written in c++11, use the MPI and OpenMP API's for parallelisation, and includes python tools to read/manipulate the data products produced. The code correlates binding energy sorted particle ID lists between halo catalogues, determining optimal descendant/progenitor matches using multiple snapshots, a merit function that maximises the number of shared particles using pseudo-radial moments, and a scheme for correcting halo merger tree pathologies. Focusing on VELOCIraptor catalogues for this work, we demonstrate how searching multiple snapshots spanning a dynamical time significantly reduces the number of stranded halos, those lacking a descendant or a progenitor, critically correcting poorly resolved halos. We present a new merit function that improves the distinction between primary and secondary progenitors, reducing tree pathologies. We find FOF accretion rates and merger rates show similar mass ratio dependence. The model merger rates from Poole et al, (2017) agree with the measured net growth of halos through mergers.

The scatter (${\rm\sigma_{\text{sSFR}}}$) of the specific star formation rates (sSFRs) of galaxies is a measure of the diversity in their star formation histories (SFHs) at a given mass. In this paper we employ the EAGLE simulations to study the dependence of the ${\rm \sigma_{\text{sSFR}}}$ of galaxies on stellar mass (${\rm M_{\star}}$) through the ${\rm \sigma_{\text{sSFR}}}$-${\rm M_{\star}}$ relation in $ {\rm z \sim 0-4}$. We find that the relation evolves with time, with the dispersion depending on both stellar mass and redshift. The models point to an evolving U-shape form for the ${\rm \sigma_{\text{sSFR}}}$-${\rm M_{\star}}$ relation with the scatter being minimal at a characteristic mass $M^{\star}$ of ${\rm 10^{9.5}}$ ${\rm M_{\odot}}$ and increasing both at lower and higher masses. This implication is that the diversity of SFHs increases towards both at the low- and high-mass ends. We find that active galactic nuclei feedback is important for increasing the ${\rm \sigma_{\text{sSFR}}}$ for high mass objects. On the other hand, we suggest that SNe feedback increases the ${\rm \sigma_{\text{sSFR}}}$ of galaxies at the low-mass end. We also find that excluding galaxies that have experienced recent mergers does not significantly affect the ${\rm \sigma_{\text{sSFR}}}$-${\rm M_{\star}}$ relation. Furthermore, we employ the combination of the EAGLE simulations with the radiative transfer code SKIRT to evaluate the effect of SFR/stellar mass diagnostics in the ${\rm \sigma_{\text{sSFR}}}$-${\rm M_{\star}}$ relation and find that the ${\rm SFR/M_{\star}}$ methodologies (e.g. SED fitting, UV+IR, UV+IRX-$\beta$) widely used in the literature to obtain intrinsic properties of galaxies have a large effect on the derived shape and normalization of the ${\rm \sigma_{\text{sSFR}}}$-${\rm M_{\star}}$ relation.

Understanding how Active Galactic Nuclei (AGN) evolve through cosmic time allows us to probe the physical processes that control their evolution. We use an updated model for the evolution of masses and spins of supermassive black holes (SMBHs), coupled to the latest version of the semi-analytical model of galaxy formation GALFORM using the Planck cosmology and a high resolution Millennium style dark matter simulation to make predictions for AGN and SMBH properties for $0 < z < 6$. We compare the model to the observed black hole mass function and the SMBH versus galaxy bulge mass relation at $z=0$, and compare the predicted bolometric, hard X-ray, soft X-ray and optical AGN luminosity functions to observations at $z < 6$, and find that the model is in good agreement with the observations. The model predicts that at $z<2$ and $L_{\mathrm{bol}} < 10^{43} \mathrm{ergs^{-1}}$, the AGN luminosity function is dominated by objects accreting in an Advection Dominated Accretion Flow (ADAF) disc state, while at higher redshifts and higher luminosities the dominant contribution is from objects accreting via a thin disc or at super-Eddington rates. The model also predicts that the AGN luminosity function at $z<3$ and $L_{\mathrm{bol}} < 10^{44} \mathrm{ergs^{-1}}$ is dominated by the contribution from AGN fuelled by quiescent hot halo accretion, while at higher luminosities and higher redshifts, the AGN luminosity function is dominated by the contribution from AGN fuelled by starbursts triggered by disc instabilities. We employ this model to predict the evolution of SMBH masses, Eddington ratios, and spins, finding that the median SMBH spin evolves very little for $0<z<6$.

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

We present VELOCIraptor, a massively parallel galaxy/(sub)halo finder that is also capable of robustly identifying tidally disrupted objects and separate stellar halos from galaxies. The code is written in c++11, use the MPI and OpenMP API's for parallelisation, and includes python tools to read/manipulate the data products produced. We demonstrate the power of the VELOCIraptor (sub)halo finder, showing how it can identify subhalos deep within the host that have negligible density contrasts to their parent halo. We find a subhalo mass-radial distance dependence: large subhalos with mass ratios of $\gtrsim10^{-2}$ are more common in the central regions that smaller subhalos, a result of dynamical friction and low tidal mass loss rates. This dependence is completely absent in (sub)halo finders in common use, which generally search for substructure in configuration space, yet is present in codes that track particles belonging to halos as they fall into other halos, such as HBT+. VELOCIraptor largely reproduces the dependence seen without tracking, finding a similar radial dependence to HBT+ in well resolved halos from our limited resolution fiducial simulation.

Two decades of effort have been poured into both single-dish and interferometric millimeter-wave surveys of the sky to infer the volume density of dusty star-forming galaxies (DSFGs, with SFR>100M$_\odot$ yr$^{-1}$) over cosmic time. Though obscured galaxies dominate cosmic star-formation near its peak at $z\sim2$, the contribution of such heavily obscured galaxies to cosmic star-formation is unknown beyond $z\sim2.5$ in contrast to the well-studied population of Lyman-break galaxies (LBGs) studied through deep, space- and ground-based pencil beam surveys in the near-infrared. Unlocking the volume density of DSFGs beyond $z>3$, particularly within the first 1 Gyr after the Big Bang is critical to resolving key open questions about early Universe galaxy formation: (1) What is the integrated star-formation rate density of the Universe in the first few Gyr and how is it distributed among low-mass galaxies (e.g. Lyman-break galaxies) and high-mass galaxies (e.g. DSFGs and quasar host galaxies)? (2) How and where do the first massive galaxies assemble? (3) What can the most extreme DSFGs teach us about the mechanisms of dust production (e.g. supernovae, AGB stars, grain growth in the ISM) <1 Gyr after the Big Bang? We summarize the types of observations needed in the next decade to address these questions.

We present stellar specific angular momentum $j_*$ measurements of two $z\sim 1.5$ galaxies in the KGES sample and 12 DYNAMO $z\sim 0.1$ analogues of high-redshift galaxies. We combine natural seeing integral field spectroscopic data to trace line emission out to high multiples of effective radius $r_e$, with adaptive optics assisted Keck/OSIRIS observations to trace the rapid rise in rotation curve in the inner regions. Our spaxel-wise integration method gives results that are on average within measurement uncertainty of the traditional rotation curve model method. At $z\sim 0$, combining GMOS and OSIRIS datasets improves the measurement uncertainty in $j_*$ from 13\% (GMOS only) or 16\% (OSIRIS only) to 10\%. At $z\sim 1.5$, systematics allow for at best 20\% uncertainty on $j_*$. DYNAMO analogues of high-$z$ galaxies have low $j_*$ for their stellar mass $M_*$, and low bulge-to-total light ratio $\beta$ for their $j_*/M_*$. The high-$z$ galaxy COSMOS 127977 has $j_*/M_*$ consistent with normal local disk galaxies, while UDS 78317 is consistent with local analogues. However, our high-resolution OSIRIS data reveal that UDS 78317 may be a merging system. We report a relationship between distance to the $\beta-j_*/M_*$ plane and the ratio of velocity dispersion to rotational velocity $\sigma/v_{max}$, where galaxies that deviate more from the plane are more dispersion-dominated due to turbulence. Much of the scatter in $M_*-j_*$ that is not explained by variations in the bulge-to-total ratio or evolution with redshift may be driven by increased turbulence due to star formation, or by treating mergers as rotating disks.

The Sloan Digital Sky Survey (SDSS) was revolutionary because of the extraordinary breadth and ambition of its optical imaging and spectroscopy. We argue that a 'sub-millimeter SDSS' - a sensitive large-area imaging+spectroscopic survey in the sub-mm window - will revolutionize our understanding of galaxy evolution in the early Universe. By detecting the thermal dust continuum emission and atomic and molecular line emission of galaxies out to z~10 it will be possible to measure the redshifts, star formation rates, dust and gas content of hundreds of thousands of high-z galaxies down to ~L*. Many of these galaxies will have counterparts visible in the deep optical imaging of the Large Synoptic Survey Telescope. This 3D map of galaxy evolution will span the peak epoch of galaxy formation all the way back to cosmic dawn, measuring the co-evolution of the star formation rate density and molecular gas content of galaxies, tracking the production of metals and charting the growth of large-scale structure.

We present an exploration of the expected detection of the earliest Active Galactic Nuclei (AGN) in the Universe from state-of-art galaxy formation and evolution semi-analytic models and hydro-dynamical simulations. We estimate the number and radiative characteristics of Super Massive Black Holes (SMBHs) at $z\geq 6$, a redshift range that will be intensively explored by the next generation of telescopes, in particular in the radio through the Square Kilometre Array (SKA) and at high energies with ESA's Athena X-ray Observatory. We find that Athena will be able to observe over 5000 AGN/$\rm deg^2$ at the Epoch of Re-ionization (EoR), $6\leq z \leq 10$. Similarly, for the same redshift range the models/simulations suggest that SKA will detect at least 400 AGN/$\rm deg^2$. Additionally, we stress the importance of the volume of the simulation box as well as the initial physical conditions of the models/simulations on their effect on the luminosity functions (LFs) and the creation of the most massive SMBHs that we currently observe at the EoR. Furthermore, following the evolution of the accretion mode of the SMBHs in each model/simulation, we show that, while the quasar dominates over the radio mode at the EoR, detection at radio wavelengths still reaches significant numbers even at the highest redshifts. Finally, we present the effect that the radiative efficiency has on the LFs by comparing results produced with a constant value for the radiative efficiency and more complex calculations based on the spin of each SMBH.

We investigate the connection between the morphology and internal kinematics of the stellar component of central galaxies with mass $M_\star > {10}^{9.5} {\rm M}_\odot$ in the EAGLE simulations. We compare several kinematic diagnostics commonly used to describe simulated galaxies, and find good consistency between them. We model the structure of galaxies as ellipsoids and quantify their morphology via the ratios of their principal axes. We show that the differentiation of blue star-forming and red quiescent galaxies using morphological diagnostics can be achieved with similar efficacy to the use of kinematical diagnostics, but only if one is able to measure both the flattening and the triaxiality of the galaxy. Flattened oblate galaxies exhibit greater rotational support than their spheroidal counterparts, but there is significant scatter in the relationship between morphological and kinematical diagnostics, such that kinematically-similar galaxies can exhibit a broad range of morphologies. The scatter in the relationship between the flattening and the ratio of the rotation and dispersion velocities ($v/\sigma$) correlates strongly with the anisotropy of the stellar velocity dispersion: at fixed $v/\sigma$, flatter galaxies exhibit greater dispersion in the plane defined by the intermediate and major axes than along the minor axis, indicating that the morphology of simulated galaxies is influenced significantly by the structure of their velocity dispersion. The simulations reveal that this anisotropy correlates with the intrinsic morphology of the galaxy's inner dark matter halo, i.e. the halo's morphology that emerges in the absence of dissipative baryonic physics. This implies the existence of a causal relationship between the morphologies of galaxies and that of their host dark matter haloes.

At fixed stellar mass, satellite galaxies show higher passive fractions than centrals, suggesting that environment is directly quenching their star formation. Here, we investigate whether satellite quenching is accompanied by changes in stellar spin (quantified by the ratio of the rotational to dispersion velocity V/$\sigma$) for a sample of massive ($M_{*}>$10$^{10}$ M$_{\odot}$) satellite galaxies extracted from the SAMI Galaxy Survey. These systems are carefully matched to a control sample of main sequence, high $V/\sigma$ central galaxies. As expected, at fixed stellar mass and ellipticity, satellites have lower star formation rate (SFR) and spin than the control centrals. However, most of the difference is in SFR, whereas the spin decreases significantly only for satellites that have already reached the red sequence. We perform a similar analysis for galaxies in the EAGLE hydro-dynamical simulation and recover differences in both SFR and spin similar to those observed in SAMI. However, when EAGLE satellites are matched to their `true' central progenitors, the change in spin is further reduced and galaxies mainly show a decrease in SFR during their satellite phase. The difference in spin observed between satellites and centrals at $z\sim$0 is primarily due to the fact that satellites do not grow their angular momentum as fast as centrals after accreting into bigger halos, not to a reduction of $V/\sigma$ due to environmental effects. Our findings highlight the effect of progenitor bias in our understanding of galaxy transformation and they suggest that satellites undergo little structural change before and during their quenching phase.

Recently a number of studies have found a similarity between the passive fraction of central and satellite galaxies when controlled for both stellar and halo mass. These results suggest that the quenching processes that affect galaxies are largely agnostic to central/satellite status, which contradicts the traditional picture of increased satellite quenching via environmental processes such as stripping, strangulation and starvation. Here we explore this further using the Galaxy And Mass Assembly (GAMA) survey which extends to ~2dex lower in stellar mass than SDSS, is more complete for closely-separated galaxies (>95% compared to >70%), and identifies lower-halo-mass groups outside of the very local Universe (M$_{\mathrm{halo}}\sim10^{12}$M$_{\odot}$ at $0.1<z<0.2$). As far as possible we aim to replicate the selections, completeness corrections and central/satellite division of one of the previous studies but find clear differences between passive fractions of centrals and satellites. We also find that our passive fractions increase with both halo-to-satellite mass ratio and central-to-second rank mass ratio. This suggests that quenching is more efficient in satellites that are low-mass for their halo ($i.e$ at high halo-to-satellite mass ratio in comparison to low halo-to-satellite mass ratio) and are more likely to be passive in older groups - forming a consistent picture of environmental quenching of satellites. We then discuss potential explanations for the previously observed similarity, such as dependence on the group-finding method.

Recently a number of studies have proposed that the dispersion along the star formation rate - stellar mass relation ($\sigma_{\mathrm{sSFR}}$-M$_{*}$) is indicative of variations in star-formation history (SFH) driven by feedback processes. They found a 'U'-shaped dispersion and attribute the increased scatter at low and high stellar masses to stellar and active galactic nuclei feed-back respectively. However, measuring $\sigma_{\mathrm{sSFR}}$ and the shape of the $\sigma_{\mathrm{sSFR}}$-M$_{*}$ relation is problematic and can vary dramatically depending on the sample selected, chosen separation of passive/star-forming systems, and method of deriving star-formation rates ($i.e.$ H$\alpha$ emission vs spectral energy distribution fitting). As such, any astrophysical conclusions drawn from measurements of $\sigma_{\mathrm{sSFR}}$ must consider these dependencies. Here we use the Galaxy And Mass Assembly survey to explore how $\sigma_{\mathrm{sSFR}}$ varies with SFR indicator for a variety of selections for disc-like `main sequence' star-forming galaxies including colour, star-formation rate, visual morphology, bulge-to-total mass ratio, S\'{e}rsic index and mixture modelling. We find that irrespective of sample selection and/or SFR indicator, the dispersion along the sSFR-M$_{*}$ relation does follow a 'U'-shaped distribution. This suggests that the shape is physical and not an artefact of sample selection or method. We then compare the $\sigma_{\mathrm{sSFR}}$-M$_{*}$ relation to state-of-the-art hydrodynamical and semi-analytic models and find good agreement with our observed results. Finally, we find that for group satellites this 'U'-shaped distribution is not observed due to additional high scatter populations at intermediate stellar masses.

We study the origin of bulge stars and their angular momentum (AM) evolution in 10 spiral galaxies with baryonic masses above $10^{10}$M$_\odot$ in the NIHAO galaxy formation simulations. The simulated galaxies are in good agreement with observations of the relation between specific AM and mass of the baryonic component and the stellar bulge-to-total ratio ($B/T$). We divide the star particles at $z=0$ into disc and bulge components using a hybrid photometric/kinematic decomposition method that identifies all central mass above an exponential disc profile as the `bulge'. By tracking the bulge star particles back in time, we find that on average 95\% of the bulge stars formed {\it in situ}, 3\% formed {\it ex situ} in satellites of the same halo, and only 2\% formed {\it ex situ} in external galaxies. The evolution of the AM distribution of the bulge stars paints an interesting picture: the higher the final $B/T$ ratio, the more the specific AM remains preserved during the bulge formation. In all cases, bulge stars migrate significantly towards the central region, reducing their average galactocentric radius by roughly a factor 2, independently of the final $B/T$ value. However, in the higher $B/T$ ($\gtrsim0.2$) objects, the velocity of the bulge stars increases and the AM of the bulge is almost conserved, whereas at lower $B/T$ values, the velocity of the bulge stars decreases and the AM of bulge reduces. The correlation between the evolution of the AM and $B/T$ suggests that bulge and disc formation are closely linked and cannot be treated as independent processes.

We analyze 88 independent high-resolution cosmological zoom-in simulations of disk galaxies in the NIHAO simulations suite to explore the connection between the atomic gas fraction and angular momentum of baryons throughout cosmic time. The study is motivated by the analytic model of \citet{obreschkow16}, which predicts a relation between the atomic gas fraction $f_{\rm atm}$ and the global atomic stability parameter $q \equiv j\sigma / (GM)$, where $M$ and $j$ are the mass and specific angular momentum of the galaxy (stars+cold gas) and $\sigma$ is the velocity dispersion of the atomic gas. We show that the simulated galaxies follow this relation from their formation ($z\simeq4$) to present within $\sim 0.5$ dex. To explain this behavior, we explore the evolution of the local Toomre stability and find that $90\%$--$100\%$ of the atomic gas in all simulated galaxies is stable at any time. In other words, throughout the entire epoch of peak star formation until today, the timescale for accretion is longer than the timescale to reach equilibrium, thus resulting in a quasi-static equilibrium of atomic gas at any time. Hence, the evolution of $f_{\rm atm}$ depends on the complex hierarchical growth history primarily via the evolution of $q$. An exception are galaxies subject to strong environmental effects.

Until a decade ago, galaxy formation simulations were unable to simultaneously reproduce the observed angular momentum (AM) of galaxy disks and bulges. Improvements in the interstellar medium and stellar feedback modelling, together with advances in computational capabilities, have allowed the current generation of cosmological galaxy formation simulations to reproduce the diversity of AM and morphology that is observed in local galaxies. In this review I discuss where we currently stand in this area from the perspective of hydrodynamical simulations, specifically how galaxies gain their AM, and the effect galaxy mergers and gas accretion have on this process. I discuss results which suggest that a revision of the classical theory of disk formation is needed, and finish by discussing what the current challenges are.

Identifying galaxies in hydrodynamical simulations is a difficult task, particularly in regions of high density such as galaxy groups and clusters. We present a new scale-free shape-independent algorithm to robustly and accurately identify galaxies in simulation, implemented within the phase-space halo-finder code VELOCIraptor. This is achieved by using the full phase-space dispersion tensor for particle assignment and an iterative adjustment of search parameters, which help us overcome common structure finding problems. We apply our improved method to the Horizon-AGN simulation and compare galaxy stellar masses ($M_*$), star formation rates (SFR) and sizes with the elaborate configuration-space halo finder, HaloMaker. Galaxies living in halos with $> 1$ galaxy are the most affected by the shortcomings of real-space finders, with their mass, SFR, and sizes being $> 2$ times larger (smaller) in the case of host (satellite) galaxies. Thus, our ability to measure minor/major merger rates and disentangle environmental effects in simulations can be generally hindered if the identification of galaxies is not treated carefully. Though large systematic differences are obtained on a one-to-one basis, the overall Galaxy Stellar Mass Function, the Star Formation Rate Function and mass-size relations are not greatly affected. This is due to isolated galaxies being the most abundant population, dominating broad statistics.

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.

Cold neutral gas is a key ingredient for growing the stellar and central black hole mass in galaxies throughout cosmic history. We have used the Atacama Large Millimetre Array (ALMA) to detect a rare example of redshifted $^{12}$CO(2-1) absorption in PKS B1740-517, a young ($t \sim 1.6 \times 10^{3}$ yr) and luminous ($L_{\rm 5 GHz} \sim 6.6 \times 10^{43}$ erg s$^{-1}$ ) radio galaxy at $z = 0.44$ that is undergoing a tidal interaction with at least one lower-mass companion. The coincident HI 21-cm and molecular absorption have very similar line profiles and reveal a reservoir of cold gas ($M_{\rm gas} \sim 10^{7} - 10^{8}$ M$_{\odot}$), likely distributed in a disc or ring within a few kiloparsecs of the nucleus. A separate HI component is kinematically distinct and has a very narrow line width ($\Delta{v}_{\rm FWHM} \lesssim 5$ km s$^{-1}$), consistent with a single diffuse cloud of cold ($T_{\rm k} \sim 100$ K) atomic gas. The $^{12}$CO(2-1) absorption is not associated with this component, which suggests that the cloud is either much smaller than 100 pc along our sight-line and/or located in low-metallicity gas that was possibly tidally stripped from the companion. We argue that the gas reservoir in PKS B1740-517 may have accreted onto the host galaxy $\sim$50 Myr before the young radio AGN was triggered, but has only recently reached the nucleus. This is consistent with the paradigm that powerful luminous radio galaxies are triggered by minor mergers and interactions with low-mass satellites and represent a brief, possibly recurrent, active phase in the life cycle of massive early type galaxies.

We present a new, open source, free semi-analytic model (SAM) of galaxy formation, Shark, designed to be highly flexible and modular, allowing easy exploration of different physical processes and ways of modelling them. We introduce the philosophy behind Shark and provide an overview of the physical processes included in the model. Shark is written in C++11 and has been parallelized with OpenMP. In the released version (v1.1), we implement several different models for gas cooling, active galactic nuclei, stellar and photo-ionisation feedback, and star formation (SF). We demonstrate the basic performance of Shark using the Planck15 cosmology SURFS simulations, by comparing against a large set of observations, including: the stellar mass function (SMF) and stellar-halo mass relation at z=0-4; the cosmic evolution of the star formation rate density (SFRD), stellar mass, atomic and molecular hydrogen; local gas scaling relations; and structural galaxy properties, finding excellent agreement. Significant improvements over previous SAMs are seen in the mass-size relation for disks/bulges, the gas-stellar mass and stellar mass-metallicity relations. To illustrate the power of Shark in exploring the systematic effects of the galaxy formation modelling, we quantify how the scatter of the SF main sequence and the gas scaling relations changes with the adopted SF law, and the effect of the starbursts H$_2$ depletion timescale on the SFRD and $\Omega_{\rm H_2}$. We compare Shark with other SAMs and the hydrodynamical simulation EAGLE, and find that SAMs have a much higher halo baryon fractions due to large amounts of intra-halo gas, which in the case of EAGLE is in the intergalactic medium.

We study the formation and evolution of ring galaxies in the Evolution and Assembly of GaLaxies and their Environments (EAGLE) simulations. We use the largest reference model Ref-L100N1504, a cubic cosmological volume of 100 comoving megaparsecs on a side, to identify and characterise these systems through cosmic time. The number density of ring galaxies in EAGLE is in broad agreement with the observations. The vast majority of ring galaxies identified in EAGLE (83 per cent) have an interaction origin, i.e., form when one or more companion galaxies drop-through a disk galaxy. The remainder (17 per cent) have very long-lived ring morphologies (> 2 Gyr) and host strong bars. Ring galaxies are HI rich galaxies, yet display inefficient star formation activity and tend to reside in the green valley particularly at z > 0.5. This inefficiency is mainly due to the low pressure and metallicity of their interstellar medium (ISM) compared with the ISM of similar star-forming galaxies. We find that the interaction(s) is responsible for decreasing the ISM pressure by causing the ISM gas to flow from the inner regions to the outer disk, where the ring feature forms. At a fixed radius, the star formation efficiency of ring galaxies is indistinguishable from their star-forming counterparts, and thus the main reason for their integrated lower efficiency is the different gas surface density profiles. Since galaxy morphologies are not used to tune the parameters in hydrodynamical simulations, the experiment performed here demonstrates the success of the current numerical models in EAGLE.

We use the semi-analytic model (SAM) of galaxy formation and evolution SAG coupled with the MULTIDARK simulation MDPL2 to study the evolution of the stellar mass-gas metallicity relation of galaxies (MZR). We test several implementations of the dependence of the mass loading due to supernovae (SNe). We find that no evolution in the normalization of the MZR is obtained unless we introduce an explicit scaling of the reheated and ejected mass with redshift as $(1+z)^\beta$. The latter is in agreement with results from the FIRE simulations, and it should encompass small scale properties of the interstellar medium varying over time, which are not captured in SAMs, as well as other energy sources in addition to SNe. Increasing $\beta$ leads to stronger evolution of the MZR normalization; $\beta = 1.9$ reproduces the observed MZR in the range $0 < z < 3.5$. A stronger redshift dependence of outflows reduces the levels of star formation at earlier epochs with the consequent decrease of metal production. This leads to a slower increase of the gas metallicity compared to the stellar mass build-up. The cold gas can be contaminated either by receiving a direct injection of the material recycled by stellar winds and SNe or by gas cooling. The relative role of each process for a given stellar mass depends on the criterion adopted to regulate the fate of the recycled material. However, modifying the metal loading of the outflows has mild impact on the zero-point evolution and does not affect our conclusions.

The total specific angular momentum j of a galaxy disk is matched with that of its dark matter halo, but the distributions are different, in that there is a lack of both low- and high-j baryons with respect to the CDM predictions. I illustrate how the probability density function PDF(j/j_mean) can inform us of a galaxy's morphology and evolutionary history with a spanning set of examples from present-day galaxies and a galaxy at z~1.5. The shape of PDF(j/j_mean) is correlated with photometric morphology, with disk-dominated galaxies having more symmetric PDF(j/j_mean) and bulge-dominated galaxies having a strongly-skewed PDF(j/j_mean). Galaxies with bigger bulges have more strongly-tailed PDF(j/j_mean), but disks of all sizes have a similar PDF(j/j_mean). In future, PDF(j/j_mean) will be useful as a kinematic decomposition tool.

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.

We investigate the origin of the type III (up-bending) discs based on a sample of 1808 galaxies from SDSS and a sub-sample of 286 galaxies with HI data from ALFALFA. We examine how the type III fraction f3, the fraction of disc galaxies which host up-bending disc breaks, depends on other galactic properties. We confirm that f3 strongly depends on the stellar concentration of galaxies. We find that HI-rich galaxies with low spins tend to have significantly more type III disc breaks than control galaxies, which are matched in concentration and stellar mass. This trend is independent of the existence of strong bars or environment of the galaxies. This result is broadly consistent with predictions from theoretical simulations, and indicates in-situ star formation fueled by gas accretion to be an important mechanism that builds the outer discs of type III galaxies.

The Deep Extragalactic VIsible Legacy Survey (DEVILS) is a large spectroscopic campaign at the Anglo-Australian Telescope (AAT) aimed at bridging the near and distant Universe by producing the highest completeness survey of galaxies and groups at intermediate redshifts ($0.3<z<1.0$). Our sample consists of $\sim$60,000 galaxies to Y$<$21.2mag, over $\sim$6deg$^{2}$ in three well-studied deep extragalactic fields (Cosmic Origins Survey field, COSMOS, Extended Chandra Deep Field South, ECDFS and the X-ray Multi-Mirror Mission Large-Scale Structure region, XMM-LSS - all Large Synoptic Survey Telescope deep-drill fields). This paper presents the broad experimental design of DEVILS. Our target sample has been selected from deep Visible and Infrared Survey Telescope for Astronomy (VISTA) Y-band imaging (VISTA Deep Extragalactic Observations, VIDEO and UltraVISTA), with photometry measured by ProFound. Photometric star/galaxy separation is done on the basis of NIR colours, and has been validated by visual inspection. To maximise our observing efficiency for faint targets we employ a redshift feedback strategy, which continually updates our target lists, feeding back the results from the previous night's observations. We also present an overview of the initial spectroscopic observations undertaken in late 2017 and early 2018.

We use the EAGLE simulations to study the oxygen abundance gradients of gas discs in galaxies within the stellar mass range [10^9.5, 10^10.8]Mo at z=0. The estimated median oxygen gradient is -0.011 (0.002) dex kpc^-1, which is shallower than observed. No clear trend between simulated disc oxygen gradient and galaxy stellar mass is found when all galaxies are considered. However, the oxygen gradient shows a clear correlation with gas disc size so that shallower abundance slopes are found for increasing gas disc sizes. Positive oxygen gradients are detected for ~40 per cent of the analysed gas discs, with a slight higher frequency in low mass galaxies. Galaxies that have quiet merger histories show a positive correlation between oxygen gradient and stellar mass, so that more massive galaxies tend to have shallower metallicity gradients. At high stellar mass, there is a larger fraction of rotational-dominated galaxies in low density regions. At low stellar mass, non-merger galaxies show a large variety of oxygen gradients and morphologies. The normalization of the disc oxygen gradients in non-merger galaxies by the effective radius removes the trend with stellar mass. Conversely, galaxies that experienced mergers show a weak relation between oxygen gradient and stellar mass. Additionally, the analysed EAGLE discs show no clear dependence of the oxygen gradients on local environment, in agreement with current observational findings.

Spatially resolved HI studies of dwarf galaxies have provided a wealth of precision data. However these high-quality, resolved observations are only possible for handful of dwarf galaxies in the Local Volume. Future HI surveys are unlikely to improve the current situation. We therefore explore a method for estimating the surface density of the atomic gas from global HI parameters, which are conversely widely available. We perform empirical tests using galaxies with resolved HI maps, and find that our approximation produces values for the surface density of atomic hydrogen within typically 0.5dex of the true value. We apply this method to a sample of 147 galaxies drawn from modern near-infrared stellar photometric surveys. With this sample we confirm a strict correlation between the atomic gas surface density and the star formation rate surface density, that is vertically offset from the Kennicutt-Schmidt relation by a factor of 10-30, and significantly steeper than the classical N=1.4 of Kennicutt (1998). We further infer the molecular fraction in the sample of low surface brightness, predominantly dwarf galaxies by assuming that the star formation relationship with molecular gas observed for spiral galaxies also holds in these galaxies, finding a molecular-to-atomic gas mass fraction within the range of 5-15%. Comparison of the data to available models shows that a model in which the thermal pressure balances the vertical gravitational field captures better the shape of the $\Sigma_{\rm SFR}-\Sigma_{\rm gas}$ relationship. However, such models fail to reproduce the data completely, suggesting that thermal pressure plays an important role in the disks of dwarf galaxies.

A decade of study has established that the molecular gas properties of star-forming galaxies follow coherent scaling relations out to z~3, suggesting remarkable regularity of the interplay between molecular gas, star formation, and stellar growth. Passive galaxies, however, are expected to be gas-poor and therefore faint, and thus little is known about molecular gas in passive galaxies beyond the local universe. Here we present deep Atacama Large Millimeter/submillimeter Array (ALMA) observations of CO(2-1) emission in 8 massive (Mstar ~ 10^11 Msol) galaxies at z~0.7 selected to lie a factor of 3-10 below the star-forming sequence at this redshift, drawn from the Large Early Galaxy Astrophysics Census (LEGA-C) survey. We significantly detect half the sample, finding molecular gas fractions <~0.1. We show that the molecular and stellar rotational axes are broadly consistent, arguing that the molecular gas was not accreted after the galaxies became quiescent. We find that scaling relations extrapolated from the star-forming population over-predict both the gas fraction and gas depletion time for passive objects, suggesting the existence of either a break or large increase in scatter in these relations at low specific star formation rate. Finally, we show that the gas fractions of the passive galaxies we have observed at intermediate redshifts are naturally consistent with evolution into local massive early-type galaxies by continued low-level star formation, with no need for further gas accretion or dynamical stabilization of the gas reservoirs in the intervening 6 billion years.

We present optical photometry of Hubble Space Telescope (HST) ACS/WFC data of the resolved stellar populations in the outer disc of the dwarf irregular galaxy DDO 154. The photometry reveals that young main sequence stars are almost absent from the outermost HI disc. Instead, most are clustered near the main stellar component of the galaxy. We constrain the stellar initial mass function (IMF) by comparing the luminosity function of the main sequence stars to simulated stellar populations assuming a constant star formation rate over the dynamical timescale. The best-fitting IMF is deficient in high mass stars compared to a canonical Kroupa IMF, with a best-fit slope $\alpha = -2.45$ and upper mass limit $M_U = 16\ M_{\odot}$. This top-light IMF is consistent with predictions of the Integrated Galaxy-wide IMF theory. Combining the HST images with HI data from The HI Nearby Galaxy Survey Treasury (THINGS) we determine the star formation law (SFL) in the outer disc. The fit has a power law exponent $N = 2.92 \pm0.22$ and zero point $A=4.47 \pm 0.65 \times 10^{-7} \ M_{\odot} \ \text{yr}^{-1} \ \text{kpc}^{-2}$. This is depressed compared to the Kennicutt-Schmidt Star Formation Law, but consistent with weak star formation observed in diffuse HI environments. Extrapolating the SFL over the outer disc implies that there could be significant star formation occurring that is not detectable in H$\alpha$. Last, we determine the Toomre stability parameter $Q$ of the outer disc of DDO 154 using the THINGS HI rotation curve and velocity dispersion map. 72% of the HI in our field has $Q\leq 4$ and this incorporates 96% of the observed MS stars. Hence 28% of the HI in the field is largely dormant.

We present the relation between stellar specific angular momentum $j_*$, stellar mass $M_*$, and bulge-to-total light ratio $\beta$ for THINGS, CALIFA and Romanowsky \& Fall datasets, exploring the existence of a fundamental plane between these parameters as first suggested by Obreschkow \& Glazebrook. Our best-fit $M_*-j_*$ relation yields a slope of $\alpha = 1.03 \pm 0.11$ with a trivariate fit including $\beta$. When ignoring the effect of $\beta$, the exponent $\alpha = 0.56 \pm 0.06$ is consistent with $\alpha = 2/3$ predicted for dark matter halos. There is a linear $\beta - j_*/M_*$ relation for $\beta \lesssim 0.4$, exhibiting a general trend of increasing $\beta$ with decreasing $j_*/M_*$. Galaxies with $\beta \gtrsim 0.4$ have higher $j_*$ than predicted by the relation. Pseudobulge galaxies have preferentially lower $\beta$ for a given $j_*/M_*$ than galaxies that contain classical bulges. Pseudobulge galaxies follow a well-defined track in $\beta - j_*/M_*$ space, consistent with Obreschkow \& Glazebrook, while galaxies with classical bulges do not. These results are consistent with the hypothesis that while growth in either bulge type is linked to a decrease in $j_*/M_*$, the mechanisms that build pseudobulges seem to be less efficient at increasing bulge mass per decrease in specific angular momentum than those that build classical bulges.

Taipan is a multi-object spectroscopic galaxy survey starting in 2017 that will cover 2pi steradians over the southern sky, and obtain optical spectra for about two million galaxies out to z<0.4. Taipan will use the newly-refurbished 1.2m UK Schmidt Telescope at Siding Spring Observatory with the new TAIPAN instrument, which includes an innovative 'Starbugs' positioning system capable of rapidly and simultaneously deploying up to 150 spectroscopic fibres (and up to 300 with a proposed upgrade) over the 6-deg diameter focal plane, and a purpose-built spectrograph operating from 370 to 870nm with resolving power R>2000. The main scientific goals of Taipan are: (i) to measure the distance scale of the Universe (primarily governed by the local expansion rate, H_0) to 1% precision, and the structure growth rate of structure to 5%; (ii) to make the most extensive map yet constructed of the mass distribution and motions in the local Universe, using peculiar velocities based on improved Fundamental Plane distances, which will enable sensitive tests of gravitational physics; and (iii) to deliver a legacy sample of low-redshift galaxies as a unique laboratory for studying galaxy evolution as a function of mass and environment. The final survey, which will be completed within 5 years, will consist of a complete magnitude-limited sample (i<17) of about 1.2x10^6 galaxies, supplemented by an extension to higher redshifts and fainter magnitudes (i<18.1) of a luminous red galaxy sample of about 0.8x10^6 galaxies. Observations and data processing will be carried out remotely and in a fully-automated way, using a purpose-built automated 'virtual observer' software and an automated data reduction pipeline. The Taipan survey is deliberately designed to maximise its legacy value, by complementing and enhancing current and planned surveys of the southern sky at wavelengths from the optical to the radio.

Satellite galaxies are commonly used as tracers to measure the line-of-sight velocity dispersion ($\sigma_{\rm LOS}$) of the dark matter halo associated with their central galaxy, and thereby to estimate the halo's mass. Recent observational dispersion estimates of the Local Group, including the Milky Way and M31, suggest $\sigma\sim$50 km/s, which is surprisingly low when compared to the theoretical expectation of $\sigma\sim$100s km/s for systems of their mass. Does this pose a problem for $\Lambda$CDM? We explore this tension using the {\small{SURFS}} suite of $N$-body simulations, containing over 10000 (sub)haloes with well tracked orbits. We test how well a central galaxy's host halo velocity dispersion can be recovered by sampling $\sigma_{\rm LOS}$ of subhaloes and surrounding haloes. Our results demonstrate that $\sigma_{\rm LOS}$ is biased mass proxy. We define an optimal window in $v_{\rm LOS}$ and projected distance ($D_p$) -- $0.5\lesssim D_p/R_{\rm vir}\lesssim1.0$ and $v_{\rm LOS} \lesssim0.5V_{\rm esc}$, where $R_{\rm vir}$ is the virial radius and $V_{\rm esc}$ is the escape velocity -- such that the scatter in LOS to halo dispersion is minimised - $\sigma_{\rm LOS}=(0.5\pm0.1)\sigma_{v,{\rm H}}$. We argue that this window should be used to measure line-of-sight dispersions as a proxy for mass, as it minimises scatter in the $\sigma_{\rm LOS}-M_{\rm vir}$ relation. This bias also naturally explains the results from \cite{mcconnachie2012a}, who used similar cuts when estimating $\sigma_{\rm LOS,LG}$, producing a bias of $\sigma_{\rm LG}=(0.44\pm0.14)\sigma_{v,{\rm H}}$. We conclude that the Local Group's velocity dispersion does not pose a problem for $\Lambda$CDM and has a mass of $\log M_{\rm LG, vir}/M_\odot=12.0^{+0.8}_{-2.0}$.

Recent observations from integral field spectroscopy (IFS) indicate that the fraction of galaxies that are slow rotators, $F_{\rm SR}$, depends primarily on stellar mass, with no significant dependence on environment. We investigate these trends and the formation paths of slow rotators (SRs) using the EAGLE and Hydrangea hydro-dynamical simulations. EAGLE consists of several cosmological boxes of volumes up to $(100\,\rm Mpc)^3$, while Hydrangea consists of $24$ cosmological simulations of galaxy clusters and their environment. Together they provide a statistically significant sample in the stellar mass range $10^{9.5}\,\rm M_{\odot}-10^{12.3}\,\rm M_{\odot}$, of $16,358$ galaxies. We construct IFS-like cubes and measure stellar spin parameters, $\lambda_{\rm R}$, and ellipticities, allowing us to classify galaxies into slow/fast rotators as in observations. The simulations display a primary dependence of $F_{\rm SR}$ on stellar mass, with a weak dependence on environment. At fixed stellar mass, satellite galaxies are more likely to be SRs than centrals. $F_{\rm SR}$ shows a dependence on halo mass at fixed stellar mass for central galaxies, while no such trend is seen for satellites. We find that $\approx 70$% of SRs at $z=0$ have experienced at least one merger with mass ratio $\ge 0.1$, with dry mergers being at least twice more common than wet mergers. Individual dry mergers tend to decrease $\lambda_{\rm R}$, while wet mergers mostly increase it. However, $30$% of SRs at $z=0$ have not experienced mergers, and those inhabit halos with median spins twice smaller than the halos hosting the rest of the SRs. Thus, although the formation paths of SRs can be varied, dry mergers and/or halos with small spins dominate.

We present the Synthetic UniveRses For Surveys ({\sc surfs}) simulations, a set of N-body/Hydro simulations of the concordance $\Lambda$ Cold Dark Matter (\LCDM) cosmology. These simulations use Planck cosmology, contain up to 10 billion particles and sample scales & halo masses down to $1~$kpc & $10^8{\rm M}_\odot$. We identify and track haloes from $z=24$ to today using a state-of-the-art 6D halo finder and merger tree builder. We demonstrate that certain properties of halo merger trees are numerically converged for haloes composed of $\gtrsim100$ particles. Haloes smoothly grow in mass, $V_{\rm max}$, with the mass history characterised by $\log M(a)\propto\exp\left[-(a/\beta)^\alpha\right]$ where $a$ is the scale factor, $\alpha(M)\approx0.8$ \& $\beta(M)\approx0.024$, with these parameters decreasing with decreasing halo mass. Subhaloes follow power-law cumulative mass and velocity functions, i.e. $n(>f)\propto f^{-\alpha}$ with $\alpha_{M}=0.83\pm0.01$ and $\alpha_{V_{\rm max}}=2.13\pm0.03$ for mass \& velocity respectively, independent of redshift, as seen in previous studies. The halo-to-halo scatter in amplitude is $0.9$~dex. The number of subhaloes in a halo weakly correlates with a halo's concentration $c$ \& spin $\lambda$:haloes of high $c$ \& low $\lambda$ have $60\%$ more subhaloes than similar mass haloes of low $c$ \& high $\lambda$. High cadence tracking shows subhaloes are dynamic residents, with $25\%$ leaving their host halo momentarily, becoming a backsplash subhalo, and another $20\%$ changing hosts entirely, in agreement with previous studies. In general, subhaloes have elliptical orbits, $e\approx0.6$, with periods of $2.3^{+2.1}_{-1.7}$~Gyrs. Subhaloes lose most of their mass at pericentric passage with mass loss rates of $\sim40\%$~Gyr$^{-1}$. These catalogues will be made publicly available.

We present a phenomological model of the Cosmic Spectral Energy Distribution (CSED) and the integrated galactic light (IGL) over all cosmic time. This model, based on an earlier model by Driver et al. (2013), attributes the cosmic star formation history to two processes -- firstly, chaotic clump accretion and major mergers, resulting in the early-time formation of bulges and secondly, cold gas accretion, resulting in late-time disc formation. Under the assumption of a Universal Chabrier initial mass function, we combine the Bruzual & Charlot (2003) stellar libraries, the Charlot & Fall (2000) dust attenuation prescription and template spectra for emission by dust and active galactic nuclei to predict the CSED -- pre- and post-dust attenuation -- and the IGL throughout cosmic time. The phenomological model, as constructed, adopts a number of basic axioms and empirical results and has minimal free parameters. We compare the model output, as well as predictions from the semi-analytic model GALFORM to recent estimates of the CSED out to $z=1$. By construction, our empirical model reproduces the full energy output of the Universe from the ultraviolet to the far-infrared extremely well. We use the model to derive predictions of the stellar and dust mass densities, again finding good agreement. We find that GALFORM predicts the CSED for $z < 0.3$ in good agreement with the observations. This agreement becomes increasingly poor towards $z = 1$, when the model CSED is $\sim$50 per cent fainter. The latter is consistent with the model underpredicting the cosmic star formation history. As a consequence, GALFORM predicts a $\sim$30 per cent fainter IGL.

We report the study of far-IR sizes of submillimeter galaxies (SMGs) in relation to their dust-obscured star formation rate (SFR) and active galactic nuclei (AGN) presence, determined using mid-IR photometry. We determined the millimeter-wave ($\lambda_{\rm obs}=1100 \mu$m) sizes of 69 ALMA-identified SMGs, selected with $\geq10$$\sigma$ confidence on ALMA images ($F_{\rm 1100 \mu m}=1.7$--7.4 mJy). We found that all the SMGs are located above an avoidance region in the millimeter size-flux plane, as expected by the Eddington limit for star formation. In order to understand what drives the different millimeter-wave sizes in SMGs, we investigated the relation between millimeter-wave size and AGN fraction for 25 of our SMGs at $z=1$--3. We found that the SMGs for which the mid-IR emission is dominated by star formation or AGN have extended millimeter-sizes, with respective median $R_{\rm c,e} = 1.6^{+0.34}_{-0.21}$ and 1.5$^{+0.93}_{-0.24}$ kpc. Instead, the SMGs for which the mid-IR emission corresponds to star-forming/AGN composites have more compact millimeter-wave sizes, with median $R_{\rm c,e}=1.0^{+0.20}_{-0.20}$ kpc. The relation between millimeter-wave size and AGN fraction suggests that this size may be related to the evolutionary stage of the SMG. The very compact sizes for composite star-forming/AGN systems could be explained by supermassive black holes growing rapidly during the SMG coalescing, star-formation phase.

We use the energy-balance code MAGPHYS to determine stellar and dust masses, and dust corrected star-formation rates for over 200,000 GAMA galaxies, 170,000 G10-COSMOS galaxies and 200,000 3D-HST galaxies. Our values agree well with previously reported measurements and constitute a representative and homogeneous dataset spanning a broad range in stellar mass (10^8---10^12 Msol), dust mass (10^6---10^9 Msol), and star-formation rates (0.01---100 Msol per yr), and over a broad redshift range (0.0 < z < 5.0). We combine these data to measure the cosmic star-formation history (CSFH), the stellar-mass density (SMD), and the dust-mass density (DMD) over a 12 Gyr timeline. The data mostly agree with previous estimates, where they exist, and provide a quasi-homogeneous dataset using consistent mass and star-formation estimators with consistent underlying assumptions over the full time range. As a consequence our formal errors are significantly reduced when compared to the historic literature. Integrating our cosmic star-formation history we precisely reproduce the stellar-mass density with an ISM replenishment factor of 0.50 +/- 0.07, consistent with our choice of Chabrier IMF plus some modest amount of stripped stellar mass. Exploring the cosmic dust density evolution, we find a gradual increase in dust density with lookback time. We build a simple phenomenological model from the CSFH to account for the dust mass evolution, and infer two key conclusions: (1) For every unit of stellar mass which is formed 0.0065---0.004 units of dust mass is also formed; (2) Over the history of the Universe approximately 90 to 95 per cent of all dust formed has been destroyed and/or ejected.

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.

We use three different cosmological dark matter simulations to study how the orientation of the angular momentum vector (AM) in dark matter haloes evolve with time. We find that haloes in this kind of simulations are constantly affected by a spurious change of mass, which translates into an artificial change in the orientation of the AM. After removing the haloes affected by artificial mass change, we found that the change in the orientation of the AM vector is correlated with time. The change in its angle and direction (i.e. the angle subtended by the AM vector in two consecutive timesteps) that affect the AM vector has a dependence on the change of mass that affects a halo, the time elapsed in which the change of mass occurs and the halo mass. We create a Monte-Carlo simulation that reproduces the change of angle and direction of the AM vector. We reproduce the angular separation of the AM vector since a look back time of 8.5 Gyrs to today ( $\rm \alpha$) with an accuracy of approximately 0.05 in $\rm cos(\alpha)$. We are releasing this Monte-Carlo simulation together with this publication. We also create a Monte Carlo simulation that reproduces the change of the AM modulus. We find that haloes in denser environments display the most dramatic evolution in their AM direction, as well as haloes with a lower specific AM modulus. These relations could be used to improve the way we follow the AM vector in low-resolution simulations.

It is now possible for hydrodynamical simulations to reproduce a representative galaxy population. Accordingly, it is timely to assess critically some of the assumptions of traditional semi-analytic galaxy formation models. We use the Eagle simulations to assess assumptions built into the Galform semi-analytic model, focussing on those relating to baryon cycling, angular momentum and feedback. We show that the assumption in Galform that newly formed stars have the same specific angular momentum as the total disc leads to a significant overestimate of the total stellar specific angular momentum of disc galaxies. In Eagle, stars form preferentially out of low specific angular momentum gas in the interstellar medium (ISM) due to the assumed gas density threshold for stars to form, leading to more realistic galaxy sizes. We find that stellar mass assembly is similar between Galform and Eagle but that the evolution of gas properties is different, with various indications that the rate of baryon cycling in Eagle is slower than is assumed in Galform. Finally, by matching individual galaxies between Eagle and Galform, we find that an artificial dependence of AGN feedback and gas infall rates on halo mass doubling events in Galform drives most of the scatter in stellar mass between individual objects. Put together our results suggest that the Galform semi-analytic model can be significantly improved in light of recent advances.

We examine the z = 0 group-integrated stellar and cold baryonic (stars + cold atomic gas) mass functions (group SMF and CBMF) and the baryonic collapse efficiency (group cold baryonic to dark matter halo mass ratio) using the RESOLVE and ECO survey galaxy group catalogs and a galform semi-analytic model (SAM) mock catalog. The group SMF and CBMF fall off more steeply at high masses and rise with a shallower low-mass slope than the theoretical halo mass function (HMF). The transition occurs at group-integrated cold baryonic mass M_coldbary ~ 10^11 Msun. The SAM, however, has significantly fewer groups at the transition mass ~ 10^11 Msun and a steeper low-mass slope than the data, suggesting that feedback is too weak in low-mass halos and conversely too strong near the transition mass. Using literature prescriptions to include hot halo gas and potential unobservable galaxy gas produces a group BMF with slope similar to the HMF even below the transition mass. Its normalization is lower by a factor of ~2, in agreement with estimates of warm-hot gas making up the remaining difference. We compute baryonic collapse efficiency with the halo mass calculated two ways, via halo abundance matching (HAM) and via dynamics (extended all the way to three-galaxy groups using stacking). Using HAM, we find that baryonic collapse efficiencies reach a flat maximum for groups across the halo mass range of M_halo ~ 10^11.4-12 Msun, which we label "nascent groups." Using dynamics, however, we find greater scatter in baryonic collapse efficiencies, likely indicating variation in group hot-to-cold baryon ratios. Similarly, we see higher scatter in baryonic collapse efficiencies in the SAM when using its true groups and their group halo masses as opposed to friends-of-friends groups and HAM masses.

We investigate measuring the peculiar motions of galaxies up to $z=0.5$ using Type Ia supernovae (SNe Ia) from LSST, and predict the subsequent constraints on the growth rate of structure. We consider two cases. Our first is based on measurements of the volumetric SNe Ia rate and assumes we can obtain spectroscopic redshifts and light curves for varying fractions of objects that are detected pre-peak luminosity by LSST (some of which may be obtained by LSST itself and others which would require additional follow-up). We find that these measurements could produce growth rate constraints at $z<0.5$ that significantly outperform those using Redshift Space Distortions (RSD) with DESI or 4MOST, even though there are $\sim4\times$ fewer objects. For our second case, we use semi-analytic simulations and a prescription for the SNe Ia rate as a function of stellar mass and star formation rate to predict the number of LSST SNe IA whose host redshifts may already have been obtained with the Taipan+WALLABY surveys, or with a future multi-object spectroscopic survey. We find $\sim 18,000$ and $\sim 160,000$ SN Ia with host redshifts for these cases respectively. Whilst this is only a fraction of the total LSST-detected SNe Ia, they could be used to significantly augment and improve the growth rate constraints compared to only RSD. Ultimately, we find that combining LSST SNe Ia with large numbers of galaxy redshifts will provide the most powerful probe of large scale gravity in the $z<0.5$ regime over the coming decades.

We built a multi-wavelength dataset for galaxies from the Local Volume HI Survey (LVHIS), which comprises 82 galaxies. We also select a sub-sample of ten large galaxies for investigating properties in the galactic outskirts. The LVHIS sample covers nearly four orders of magnitude in stellar mass and two orders of magnitude in HI mass fraction (fHI). The radial distribution of HI gas with respect to the stellar disc is correlated with fHI but with a large scatter. We confirm the previously found correlations between the total HI mass and star formation rate (SFR), and between HI surface densities and SFR surface densities beyond R25. However, the former correlation becomes much weaker when the average surface densities rather than total mass or rate are considered, and the latter correlation also becomes much weaker when the effect of stellar mass is removed or controlled. Hence the link between SFR and HI is intrinsically weak in these regions, consistent with what was found on kpc scales in the galactic inner regions. We find a strong correlation between the SFR surface density and the stellar mass surface density, which is consistent with the star formation models where the gas is in quasi-equilibrium with the mid-plane pressure. We find no evidence for HI warps to be linked with decreasing star forming efficiencies.

We investigate the evolution of the galaxy Star Formation Rate Function (SFRF) and Cosmic Star Formation Rate Density (CSFRD) of $z\sim 0-8 $ galaxies in the Evolution and Assembly of GaLaxies and their Environments (EAGLE) simulations. In addition, we present a compilation of UV, IR and H$\alpha$ SFRFs and compare these with the predictions from the EAGLE suite of cosmological hydrodynamic simulations. We find that the constraints implied by different indicators are inconsistent with each other for the highest star-forming objects at z < 2, a problem that is possibly related to selection biases and the uncertainties of dust attenuation effects. EAGLE's feedback parameters were calibrated to reproduce realistic galaxy sizes and stellar masses at z = 0.1. In this work we test if and why those choices yield realistic Star Formation Rates (SFRs) for $z \sim 0-8$ as well. We demonstrate that SNe feedback plays a major role at setting the abundance of galaxies at all star-forming regimes, especially at high redshifts. On the contrary, Active Galactic Nuclei (AGN) feedback becomes more prominent at lower redshifts and is a major mechanism that affects only the highest star-forming systems. Furthermore, we find that galaxies with SFR $\sim 1-10 \, {\rm M_{\odot} \, yr^{-1}}$ dominate the CSFRD at redshifts z < 5, while rare high star-forming galaxies (SFR $\sim 10-100 \,{\rm M_{\odot} \, yr^{-1}}$) contribute significantly only briefly around the peak era ($z \sim 2$) and then are quenched by AGN feedback. In the absence of this prescription objects with SFR $\sim 10-100 \,{\rm M_{\odot} \, yr^{-1}}$ would dominate the CSFRD, while the cosmic budget of star formation would be extremely high. Finally, we demonstrate that the majority of the cosmic star formation occurs in relatively rare high mass halos ($ {\rm M_{Halo}} \sim 10^{11-13} \, {\rm M_{\odot}}$) even at the earliest epochs.

Cosmic voids are promising tools for cosmological tests due to their sensitivity to dark energy, modified gravity and alternative cosmological scenarios. Most previous studies in the literature of void properties use cosmological N-body simulations of dark matter (DM) particles that ignore the potential effect of baryonic physics. Using a spherical underdensity finder, we analyse voids using the mass field and subhalo tracers in the EAGLE simulations, which follow the evolution of galaxies in a $\rm{\Lambda}$ cold dark matter Universe with state-of-the-art subgrid models for baryonic processes in a $(100 \rm{cMpc})^3$ volume. We study the effect of baryons on void statistics by comparing results with DM-only simulations that use the same initial conditions as EAGLE. When identifying voids in the mass field, we find that a DM-only simulation produces 24 per cent more voids than a hydrodynamical one due to the action of galaxy feedback polluting void regions with hot gas, specially for small voids with $r_{\rm{void}} \le 10\ \rm{Mpc}$. We find that the way in which galaxy tracers are selected has a strong impact on the inferred void properties. Voids identified using galaxies selected by their stellar mass are larger and have cuspier density profiles than those identified by galaxies selected by their total mass. Overall, baryons have minimal effects on void statistics, as void properties are well captured by DM-only simulations, but it is important to account for how galaxies populate DM haloes to estimate the observational effect of different cosmological models on the statistics of voids.

We explore the multiwavelength properties of AGN host galaxies for different classes of radio-selected AGN out to z$\lesssim$6 via a multiwavelength analysis of about 7700 radio sources in the COSMOS field. The sources were selected with the Very Large Array (VLA) at 3 GHz (10 cm) within the VLA-COSMOS 3 GHz Large Project, and cross-matched with multiwavelength ancillary data. This is the largest sample of high-redshift (z$\lesssim$6) radio sources with exquisite photometric coverage and redshift measurements available. We constructed a sample of moderate-to-high radiative luminosity AGN (HLAGN) via spectral energy distribution (SED) decomposition combined with standard X-ray and mid-infrared diagnostics. Within the remainder of the sample we further identified low-to-moderate radiative luminosity AGN (MLAGN) via excess in radio emission relative to the star formation rates in their host galaxies. We show that AGN power in HLAGN occurs predominantly in radiative form, while MLAGN display a substantial mechanical AGN luminosity component. We found significant differences in the host properties of the two AGN classes, as a function of redshift. At z$<$1.5, MLAGN appear to reside in significantly more massive and less star-forming galaxies compared to HLAGN. At z$>$1.5, we observed a reversal in the behaviour of the stellar mass distributions with the HLAGN populating the higher stellar mass tail. We interpret this finding as a possible hint of the downsizing of galaxies hosting HLAGN, with the most massive galaxies triggering AGN activity earlier than less massive galaxies, and then fading to MLAGN at lower redshifts. Our conclusion is that HLAGN and MLAGN samples trace two distinct galaxy and AGN populations in a wide range of redshifts, possibly resembling the radio AGN types often referred to as radiative- and jet-mode (or high- and low-excitation), respectively.

We present a MUSE and KMOS dynamical study 405 star-forming galaxies at redshift z=0.28-1.65 (median redshift z=0.84). Our sample are representative of star-forming, main-sequence galaxies, with star-formation rates of SFR=0.1-30Mo/yr and stellar masses M=10^8-10^11Mo. For 49+/-4% of our sample, the dynamics suggest rotational support, 24+/-3% are unresolved systems and 5+/-2% appear to be early-stage major mergers with components on 8-30kpc scales. The remaining 22+/-5% appear to be dynamically complex, irregular (or face-on systems). For galaxies whose dynamics suggest rotational support, we derive inclination corrected rotational velocities and show these systems lie on a similar scaling between stellar mass and specific angular momentum as local spirals with j*=J/M*\propto M^(2/3) but with a redshift evolution that scales as j*\propto M^{2/3}(1+z)^(-1). We identify a correlation between specific angular momentum and disk stability such that galaxies with the highest specific angular momentum, log(j*/M^(2/3))>2.5, are the most stable, with Toomre Q=1.10+/-0.18, compared to Q=0.53+/-0.22 for galaxies with log(j*/M^(2/3))<2.5. At a fixed mass, the HST morphologies of galaxies with the highest specific angular momentum resemble spiral galaxies, whilst those with low specific angular momentum are morphologically complex and dominated by several bright star-forming regions. This suggests that angular momentum plays a major role in defining the stability of gas disks: at z~1, massive galaxies that have disks with low specific angular momentum, appear to be globally unstable, clumpy and turbulent systems. In contrast, galaxies with high specific angular have evolved in to stable disks with spiral structures.

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.

We compare observed far infra-red/sub-millimetre (FIR/sub-mm) galaxy spectral energy distributions (SEDs) of massive galaxies ($M_{\star}\gtrsim10^{10}$ $h^{-1}$M$_{\odot}$) derived through a stacking analysis with predictions from a new model of galaxy formation. The FIR SEDs of the model galaxies are calculated using a self-consistent model for the absorption and re-emission of radiation by interstellar dust based on radiative transfer calculations and global energy balance arguments. Galaxies are selected based on their position on the specific star formation rate (sSFR) - stellar mass ($M_{\star}$) plane. We identify a main sequence of star-forming galaxies in the model, i.e. a well defined relationship between sSFR and $M_\star$, up to redshift $z\sim6$. The scatter of this relationship evolves such that it is generally larger at higher stellar masses and higher redshifts. There is remarkable agreement between the predicted and observed average SEDs across a broad range of redshifts ($0.5\lesssim z\lesssim4$) for galaxies on the main sequence. However, the agreement is less good for starburst galaxies at $z\gtrsim2$, selected here to have elevated sSFRs$>10\times$ the main sequence value. We find that the predicted average SEDs are robust to changing the parameters of our dust model within physically plausible values. We also show that the dust temperature evolution of main sequence galaxies in the model is driven by star formation on the main sequence being more burst-dominated at higher redshifts.

We present the detailed characterization of two extremely red submillimeter galaxies (SMGs), ASXDF1100.053.1 and 231.1, with the Atacama Large Millimeter/submillimeter Array (ALMA) and the Jansky Very Large Array (VLA). These SMGs were selected originally using AzTEC at 1100 micron, and are observed by Herschel to be faint at 100--500 micron. Their (sub)millimeter colors are as red as -- or redder -- than known z>~5 SMGs; indeed, ASXDF1100.053.1 is redder than HFLS 3, which lies at z=6.3. They are also faint and red in the near-/mid-infrared: ~1 microJy at IRAC 4.5 micron and <0.2 microJy in the Ks filter. These SMGs are also faint in the radio waveband, where F_6GHz=4.5 microJy for ASXDF1100.053.1 and F_1.4GHz=28 microJy for ASXDF1100.231.1, suggestive of z=6.5^{+1.4}_{-1.1} and z=4.1^{+0.6}_{-0.7} for ASXDF1100.053.1 and 231.1, respectively. ASXDF1100.231.1 has a flux excess in the 3.6-micron filter, probably due to H$\alpha$ emission at z=4--5. Derived properties of ASXDF1100.053.1 for z=5.5--7.5 and 231.1 for z=3.5--5.5 are as follows: their infrared luminosities are [6.5-7.4]x10^{12} and [4.2-4.5]x10^{12} L_sun; their stellar masses are [0.9-2]x10^{11} and [0.4-3]x10^{10} M_sun; their circularized half-light radii in the ALMA maps are ~1 and <~0.2 kpc (~2--3 kpc for 90% of the total flux). Lastly, their surface infrared luminosity densities, Sigma_IR, are ~1x10^{12} and >~1.5x10^{13} L_sun kpc^{-2}, similar to values seen for local (U)LIRGs. These data suggest that ASXDF1100.053.1 and 231.1 are compact SMGs at z>~4 and can plausibly evolve into z>~3 compact quiescent galaxies.

In this paper we investigate environment driven gas depletion in satellite galaxies, taking full advantage of the atomic hydrogen (HI) spectral stacking technique to quantify the gas content for the entire gas-poor to -rich regime. We do so using a multi-wavelength sample of 10,600 satellite galaxies, selected according to stellar mass (log M$_{\star}$/M$_{\odot}$ $\geq$ 9) and redshift (0.02 $\leq$ z $\leq$ 0.05) from the Sloan Digital Sky Survey, with HI data from the Arecibo Legacy Fast ALFA (ALFALFA) survey. Using key HI-to-stellar mass scaling relations, we present evidence that the gas content of satellite galaxies is, to a significant extent, dependent on the environment in which a galaxy resides. For the first time, we demonstrate that systematic environmental suppression of gas content at both fixed stellar mass and fixed specific star formation rate (sSFR) in satellite galaxies begins in halo masses typical of the group regime (log M$_{h}$/M$_{\odot}$ < 13.5), well before galaxies reach the cluster environment. We also show that environment driven gas depletion is more closely associated to halo mass than local density. Our results are then compared with state-of-the-art semi-analytic models and hydrodynamical simulations and discussed within this framework, showing that more work is needed if models are to reproduce the observations. We conclude that the observed decrease of gas content in the group and cluster environments cannot be reproduced by starvation of the gas supply alone and invoke fast acting processes such as ram-pressure stripping of cold gas to explain this.

The formation and evolution of galaxies with low neutral atomic hydrogen (HI) masses, M$_{\rm HI}$$<$10$^{8}h^{-2}$M$_{\odot}$, are affected by host dark matter halo mass and photoionisation feedback from the UV background after the end of reionization. We study how the physical processes governing the formation of galaxies with low HI mass are imprinted on the distribution of neutral hydrogen in the Universe using the hierarchical galaxy formation model, GALFORM. We calculate the effect on the correlation function of changing the HI mass detection threshold at redshifts $0 \le z \le 0.5$. We parameterize the clustering as $\xi(r)=(r/r_{0})^{-\gamma}$ and we find that including galaxies with M$_{\rm HI}$$<$10$^{8}h^{-2}$M$_{\odot}$ increases the clustering amplitude $r_{0}$ and slope $\gamma$ compared to samples of higher HI masses. This is due to these galaxies with low HI masses typically being hosted by haloes with masses greater than 10$^{12}{h}^{-1}$M$_{\odot}$, and is in contrast to optically selected surveys for which the inclusion of faint, blue galaxies lowers the clustering amplitude. We show the HI mass function for different host dark matter halo masses and galaxy types (central or satellite) to interpret the values of $r_{0}$ and $\gamma$ of the clustering of HI-selected galaxies. We also predict the contribution of low HI mass galaxies to the 21cm intensity mapping signal. We calculate that a dark matter halo mass resolution better than $\sim$10$^{10}{h}^{-1}$M$_{\odot}$ at redshifts higher than 0.5 is required in order to predict converged 21cm brightness temperature fluctuations.

We examine the properties of atomic hydrogen (HI) associated with galaxies in the EAGLE simulations of galaxy formation. EAGLE's feedback parameters were calibrated to reproduce the stellar mass function and galaxy sizes at $z=0.1$, and we assess whether this calibration also yields realistic HI properties. We estimate the self-shielding density with a fitting function calibrated using radiation transport simulations, and correct for molecular hydrogen with empirical or theoretical relations. The `standard-resolution' simulations systematically underestimate HI column densities, leading to an HI deficiency in low-mass ($M_\star < 10^{10}M_\odot$) galaxies and poor reproduction of the observed HI mass function. These shortcomings are largely absent from EAGLE simulations featuring a factor of 8 (2) better mass (spatial) resolution, within which the HI mass of galaxies evolves more mildly from $z=1$ to $0$ than in the standard-resolution simulations. The largest-volume simulation reproduces the observed clustering of HI systems, and its dependence on HI-richness. At fixed $M_\star$, galaxies acquire more HI in simulations with stronger feedback, as they become associated with more massive haloes and higher infall rates. They acquire less HI in simulations with a greater star formation efficiency, since the star formation and feedback necessary to balance the infall rate is produced by smaller gas reservoirs. The simulations indicate that the HI of present-day galaxies was acquired primarily by the smooth accretion of ionized, intergalactic gas at $z\simeq1$, which later self-shields, and that only a small fraction is contributed by the reincorporation of gas previously heated strongly by feedback. HI reservoirs are highly dynamic: over $40$ percent of HI associated with $z=0.1$ galaxies is converted to stars or ejected by $z=0$.

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.

(Abridged) In the last decade several massive molecular gas reservoirs were found <100 kpc distance from active galactic nuclei (AGNs), residing in gas-rich companion galaxies. The study of AGN-gas-rich companion systems opens the opportunity to determine whether the stellar mass of massive local galaxies was formed in their host after a merger event or outside of their host galaxy in a close starbursting companion and later incorporated via mergers. We study the quasar-companion galaxy system of SMM J04135+10277 (z=2.84) and investigate the expected frequency of quasar-starburst galaxy pairs at high redshift using a cosmological galaxy formation model. We use archive data and new APEX ArTeMiS data to construct and model the spectral energy distribution of SMM J04135. We also carry out a comprehensive analysis of the cosmological galaxy formation model GALFORM with the aim of characterising how typical the system of SMM J04135 is and whether quasar-star-forming galaxy pairs may constitute an important stage in galaxy evolution. The companion galaxy of SMM J04135 is a heavily dust-obscured starburst galaxy with a median star formation rate (SFR) of $700\,\rm{M_{\odot}\,yr^{-1}}$, median dust mass of $5.1\times 10^9\,\rm{M_{\odot}}$ and median dust luminosity of $\textrm 9.3 \times 10^{12}\,\rm{L_{\odot}}$. Our simulations, performed at z=2.8, suggest that SMM J04135 is not unique. In fact, at a distance of <100 kpc, 22% of our simulated quasar sample have at least one companion galaxy of a stellar mass $>10^8\, \rm{M_{\odot}}$, and 0.3% have at least one highly star-forming companion ($\rm{SFR}>100\,\rm{M_{\odot}\,yr^{-1}}$). Our results suggest that quasar-gas-rich companion galaxy systems are common phenomena in the early Universe and the high incidence of companions makes the study of such systems crucial to understand the growth and hierarchical build-up of galaxies and black holes.

We investigate the contentious issue of the presence, or lack thereof, of satellites mass segregation in galaxy groups using the Galaxy And Mass Assembly (GAMA) survey, the GALFORM semi-analytic and the EAGLE cosmological hydrodynamical simulation catalogues of galaxy groups. We select groups with halo mass $12 \leqslant \log(M_{\text{halo}}/h^{-1}M_\odot) <14.5$ and redshift $z \leqslant 0.32$ and probe the radial distribution of stellar mass out to twice the group virial radius. All the samples are carefully constructed to be complete in stellar mass at each redshift range and efforts are made to regularise the analysis for all the data. Our study shows negligible mass segregation in galaxy group environments with absolute gradients of $\lesssim0.08$ dex and also shows a lack of any redshift evolution. Moreover, we find that our results at least for the GAMA data are robust to different halo mass and group centre estimates. Furthermore, the EAGLE data allows us to probe much fainter luminosities ($r$-band magnitude of 22) as well as investigate the three-dimensional spatial distribution with intrinsic halo properties, beyond what the current observational data can offer. In both cases we find that the fainter EAGLE data show a very mild spatial mass segregation at $z \leqslant 0.22$, which is again not apparent at higher redshift. Interestingly, our results are in contrast to some earlier findings using the Sloan Digital Sky Survey. We investigate the source of the disagreement and suggest that subtle differences between the group finding algorithms could be the root cause.

We present a new version of the GALFORM semi-analytical model of galaxy formation. This brings together several previous developments of GALFORM into a single unified model, including a different initial mass function (IMF) in quiescent star formation and in starbursts, feedback from active galactic nuclei supressing gas cooling in massive halos, and a new empirical star formation law in galaxy disks based on their molecular gas content. In addition, we have updated the cosmology, introduced a more accurate treatment of dynamical friction acting on satellite galaxies, and updated the stellar population model. The new model is able to simultaneously explain both the observed evolution of the K-band luminosity function and stellar mass function, and the number counts and redshift distribution of sub-mm galaxies selected at 850 mu. This was not previously achieved by a single physical model within the LambdaCDM framework, but requires having an IMF in starbursts that is somewhat top-heavy. The new model is tested against a wide variety of observational data covering wavelengths from the far-UV to sub-mm, and redshifts from z=0 to z=6, and is found to be generally successful. These observations include the optical and near-IR luminosity functions, HI mass function, fraction of early type galaxies, Tully-Fisher, metallicity-luminosity and size-luminosity relations at z=0, as well as far-IR number counts, and far-UV luminosity functions at z ~ 3-6. [abridged]

We perform automated bulge + disc decomposition on a sample of $\sim$7500 galaxies from the Galaxy And Mass Assembly (GAMA) survey in the redshift range of 0.002$<$z$<$0.06 using SIGMA, a wrapper around GALFIT3. To achieve robust profile measurements we use a novel approach of repeatedly fitting the galaxies, varying the input parameters to sample a large fraction of the input parameter space. Using this method we reduce the catastrophic failure rate significantly and verify the confidence in the fit independently of $\chi^2$. Additionally, using the median of the final fitting values and the 16$^{th}$ and 84$^{th}$ percentile produces more realistic error estimates than those provided by GALFIT, which are known to be underestimated. We use the results of our decompositions to analyse the stellar mass - half-light radius relations of bulges, discs and spheroids. We further investigate the association of components with a parent disc or elliptical relation to provide definite $z=0$ disc and spheroid $\mathcal{M_\star}-R_{\rm e}$ relations. We conclude by comparing our local disc and spheroid $\mathcal{M_\star}-R_{\rm e}$ to simulated data from EAGLE and high redshift data from CANDELS-UDS. We show the potential of using the mass-size relation to study galaxy evolution in both cases but caution that for a fair comparison all data sets need to be processed and analysed in the same manner.

We examine the evolution of intrinsic u-r colours of galaxies in the EAGLE cosmological hydrodynamical simulations, which has been shown to reproduce the observed redshift z=0.1 colour-magnitude distribution well. The median u-r of star-forming ('blue cloud') galaxies reddens by 1 mag from z=2 to 0 at fixed stellar mass, as their specific star formation rates decrease with time. A red sequence starts to build-up around z=1, due to the quenching of low-mass satellite galaxies at the faint end, and due to the quenching of more massive central galaxies by their active galactic nuclei (AGN) at the bright end. This leaves a dearth of intermediate-mass red sequence galaxies at z=1, which is mostly filled in by z=0. We quantify the time-scales of colour transition due to satellite and AGN quenching, finding that most galaxies spend less than 2 Gyr in the 'green valley'. On examining the trajectories of galaxies in a colour-stellar mass diagram, we identify three characteristic tracks that galaxies follow (quiescently star-forming, quenching and rejuvenating galaxies) and quantify the fraction of galaxies that follow each track.

We compare global predictions from the EAGLE hydrodynamical simulation, and two semi-analytic (SA) models of galaxy formation, L-GALAXIES and GALFORM. All three models include the key physical processes for the formation and evolution of galaxies and their parameters are calibrated against a small number of observables at $z\approx 0$. The two SA models have been applied to merger trees constructed from the EAGLE dark matter only simulation. We find that at $z\leq 2$, both the galaxy stellar mass functions for stellar masses $M_*<10^{10.5} M_{\odot}$ and the median specific star formation rates (sSFRs) in the three models agree to better than $0.4$~dex. The evolution of the sSFR predicted by the three models closely follows the mass assembly history of dark matter haloes. In both EAGLE and L-GALAXIES there are more central passive galaxies with $M_*<10^{9.5} M_{\odot}$ than in L-GALAXIES. This difference is related to galaxies that have entered and then left a larger halo and which are treated as satellites in GALFORM. In the range $0<z<1$, the slope of the evolution of the star formation rate density in EAGLE is a factor of $\approx 1.5$ steeper than for the two SA models. The median sizes for galaxies with $M_*>10^{9.5} M_{\odot}$ differ in some instances by an order of magnitude, while the stellar mass-size relation in EAGLE is a factor of $\approx 2$ tighter than for the two SA models. Our results suggest the need for a revision of how SA models treat the effect of baryonic self-gravity on the underlying dark matter. The treatment of gas flows in the models needs to be revised based on detailed comparison with observations to understand in particular the evolution of the stellar mass-metallicity relation.

We investigate correlations between different physical properties of star-forming galaxies in the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) cosmological hydrodynamical simulation suite over the redshift range $0\le z\le 4.5$. A principal component analysis reveals that neutral gas fraction ($f_{\rm gas, neutral}$), stellar mass ($M_{\rm stellar}$) and star formation rate (SFR) account for most of the variance seen in the population, with galaxies tracing a two-dimensional, nearly flat, surface in the three-dimensional space of $f_{\rm gas, neutral}-M_{\rm stellar}-\rm SFR$ with little scatter. The location of this plane varies little with redshift, whereas galaxies themselves move along the plane as their $f_{\rm gas, neutral}$ and SFR drop with redshift. The positions of galaxies along the plane are highly correlated with gas metallicity. The metallicity can therefore be robustly predicted from $f_{\rm gas, neutral}$, or from the $M_{\rm stellar}$ and SFR. We argue that the appearance of this "fundamental plane of star formation" is a consequence of self-regulation, with the plane's curvature set by the dependence of the SFR on gas density and metallicity. We analyse a large compilation of observations spanning the redshift range $0\lesssim \rm z\lesssim 2.5$, and find that such a plane is also present in the data. The properties of the observed fundamental plane of star formation are in good agreement with EAGLE's predictions.

High-redshift galaxy clusters allow us to examine galaxy formation in extreme environments. Here we compile data for 15 $z>1$ galaxy clusters to test the predictions from a state-of-the-art semi-analytical model of galaxy formation. The model gives a good match to the slope and zero-point of the cluster red sequence. The model is able to match the cluster galaxy luminosity function at faint and bright magnitudes, but underestimates the number of galaxies around the break in the cluster luminosity function. We find that simply assuming a weaker dust attenuation improves the model predictions for the cluster galaxy luminosity function, but worsens the predictions for the red sequence at bright magnitudes. Examination of the properties of the bright cluster galaxies suggests that the default dust attenuation is large due to these galaxies having large reservoirs of cold gas as well as small radii. We find that matching the luminosity function and colours of high-redshift cluster galaxies, whilst remaining consistent with local observations, poses a challenge for galaxy formation models.

We compare the mass and internal distribution of atomic hydrogen (HI) in 2200 present-day central galaxies with M_star > 10^10 M_Sun from the 100 Mpc EAGLE Reference simulation to observational data. Atomic hydrogen fractions are corrected for self-shielding using a fitting formula from radiative transfer simulations and for the presence of molecular hydrogen using an empirical or a theoretical prescription from the literature. The resulting neutral hydrogen fractions, M_(HI+H2)/M_star, agree with observations to better than 0.1 dex for galaxies with M_star between 10^10 and 10^11 M_Sun. Our fiducial, empirical H2 model based on gas pressure results in galactic HI mass fractions, M_HI/M_star, that agree with observations from the GASS survey to better than 0.3 dex, but the alternative theoretical H2 formula leads to a negative offset in M_HI/M_star of up to 0.5 dex. Visual inspection reveals that most HI disks in simulated HI-rich galaxies are vertically disturbed, plausibly due to recent accretion events. Many galaxies (up to 80 per cent) contain spuriously large HI holes, which are likely formed as a consequence of the feedback implementation in EAGLE. The HI mass-size relation of all simulated galaxies is close to (but 16 per cent steeper than) observed, and when only galaxies without large holes in the HI disc are considered, the agreement becomes excellent (better than 0.1 dex). The presence of large HI holes also makes the radial HI surface density profiles somewhat too low in the centre, at \Sigma_HI > 1 M_Sun pc^-2 (by a factor of <~ 2 compared to data from the Bluedisk survey). In the outer region (\Sigma_HI < 1 M_Sun pc^-2), the simulated profiles agree quantitatively with observations. Scaled by HI size, the simulated profiles of HI-rich (M_HI > 10^9.8 M_Sun) and control galaxies (10^9.1 M_Sun > M_HI > 10^9.8 M_Sun) follow each other closely, as observed. (Abridged)

In this work we present IRAM-30m telescope observations of a sample of bulge-dominated galaxies with large dust lanes, which have had a recent minor merger. We find these galaxies are very gas rich, with H2 masses between 4x10^8 and 2x10^10 Msun. We use these molecular gas masses, combined with atomic gas masses from an accompanying paper, to calculate gas-to-dust and gas-to-stellar mass ratios. The gas-to-dust ratios of our sample objects vary widely (between ~50 and 750), suggesting many objects have low gas-phase metallicities, and thus that the gas has been accreted through a recent merger with a lower mass companion. We calculate the implied minor companion masses and gas fractions, finding a median predicted stellar mass ratio of ~40:1. The minor companion likely had masses between ~10^7 - 10^10 Msun. The implied merger mass ratios are consistent with the expectation for low redshift gas-rich mergers from simulations. We then go on to present evidence that (no matter which star-formation rate indicator is used) our sample objects have very low star-formation efficiencies (star-formation rate per unit gas mass), lower even than the early-type galaxies from ATLAS3D which already show a suppression. This suggests that minor mergers can actually suppress star-formation activity. We discuss mechanisms that could cause such a suppression, include dynamical effects induced by the minor merger.

We explore the galaxy formation physics governing the low mass end of the HI mass function in the local Universe. Specifically, we predict the effects on the HI mass function of varying i) the strength of photoionisation feedback and the redshift of the end of the epoch of reionization, ii) the cosmology, iii) the supernovae feedback prescription, and iv) the efficiency of star formation. We find that the shape of the low-mass end of the HI mass function is most affected by the critical halo mass below which galaxy formation is suppressed by photoionisation heating of the intergalactic medium. We model the redshift dependence of this critical dark matter halo mass by requiring a match to the low-mass end of the HI mass function. The best fitting critical dark matter halo mass decreases as redshift increases in this model, corresponding to a circular velocity of $\sim 50 \, {\rm km \,s}^{-1}$ at $z=0$, $\sim 30 \, {\rm km\, s}^{-1}$ at $z \sim 1$ and $\sim 12 \, {\rm km \, s}^{-1}$ at $z=6$. We find that an evolving critical halo mass is required to explain both the shape and abundance of galaxies in the HI mass function below $M_{\rm HI} \sim 10^{8} h^{-2} {\rm M_{\odot}}$. The model makes specific predictions for the clustering strength of HI-selected galaxies with HI masses > $10^{6} h^{-2} {\rm M_{\odot}}$ and $> 10^{7} h^{-2} {\rm M_{\odot}}$ and for the relation between the HI and stellar mass contents of galaxies which will be testable with upcoming surveys with the Square Kilometre Array and its pathfinders. We conclude that measurements of the HI mass function at $z \ge 0$ will lead to an improvement in our understanding of the net effect of photoionisation feedback on galaxy formation and evolution.

We report the source size distribution, as measured by ALMA millimetric continuum imaging, of a sample of 13 AzTEC-selected submillimeter galaxies (SMGs) at z_photo ~ 3-6. Their infrared luminosities and star-formation rates (SFR) are L_IR ~ 2-6 x 10^12 L_sun and ~ 200-600 M_sun yr-1, respectively. The size of z ~ 3-6 SMGs ranges from 0".10 to 0".38 with a median of 0".20+0".03-0".05 (FWHM), corresponding to a median circularized effective radius (Rc,e) of 0.67+0.13-0.14 kpc, comparable to the typical size of the stellar component measured in compact quiescent galaxies at z ~ 2 (cQGs) --- R ~ 1 kpc. The median surface SFR density of our z ~ 3-6 SMGs is 100+42-26 M_sun yr-1 kpc-2, comparable to that seen in local merger-driven (U)LIRGsrather than in extended disk galaxies at low and high redshifts. The discovery of compact starbursts in z >~ 3 SMGs strongly supports a massive galaxy formation scenario wherein z ~ 3-6 SMGs evolve into the compact stellar components of z ~ 2 cQGs. These cQGs are then thought to evolve into the most massive ellipticals in the local Universe, mostly via dry mergers. Our results thus suggest that z >~ 3 SMGs are the likely progenitors of massive local ellipticals, via cQGs, meaning that we can now trace the evolutionary path of the most massive galaxies over a period encompassing ~ 90% of the age of the Universe.

We investigate the abundance of galactic molecular hydrogen (H$_2$) in the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) cosmological hydrodynamic simulations. We assign H$_2$ masses to gas particles in the simulations in post-processing using two different prescriptions that depend on the local dust-to-gas ratio and the interstellar radiation field. Both result in H$_2$ galaxy mass functions that agree well with observations in the local and high-redshift Universe. The simulations reproduce the observed scaling relations between the mass of H$_2$ and the stellar mass, star formation rate and stellar surface density. Towards high edshifts, galaxies in the simulations display larger H$_2$ mass fractions, and correspondingly lower H$_2$ depletion timescales, also in good agreement with observations. The comoving mass density of H$_2$ in units of the critical density, $\Omega_{\rm H_2}$, peaks at $z\approx 1.2-1.5$, later than the predicted peak of the cosmic star formation rate activity, at $z\approx 2$. This difference stems from the decrease in gas metallicity and increase in interstellar radiation field with redshift, both of which hamper H$_2$ formation. We find that the cosmic H$_2$ budget is dominated by galaxies with $M_{\rm H_2}>10^9\,\rm M_{\odot}$, star formation rates $>10\,\rm M_{\odot}\,\rm yr^{-1}$ and stellar masses $M_{\rm stellar}>10^{10}\,\rm M_{\odot}$, which are readily observable in the optical and near-IR. The match between the H$_2$ properties of galaxies that emerge in the simulations and observations is remarkable, particularly since H$_2$ observations were not used to adjust parameters in EAGLE.

We present predictions for the two-point correlation function of galaxy clustering as a function of stellar mass, computed using two new versions of the GALFORM semi-analytic galaxy formation model. These models make use of a high resolution, large volume N-body simulation, set in the WMAP7 cosmology. One model uses a universal stellar initial mass function (IMF), while the other assumes different IMFs for quiescent star formation and bursts. Particular consideration is given to how the assumptions required to estimate the stellar masses of observed galaxies (such as the choice of IMF, stellar population synthesis model and dust extinction) influence the perceived dependence of galaxy clustering on stellar mass. Broad-band spectral energy distribution fitting is carried out to estimate stellar masses for the model galaxies in the same manner as in observational studies. We show clear differences between the clustering signals computed using the true and estimated model stellar masses. As such, we highlight the importance of applying our methodology to compare theoretical models to observations. We introduce an alternative scheme for the calculation of the merger timescales for satellite galaxies in GALFORM, which takes into account the dark matter subhalo information from the simulation. This reduces the amplitude of small-scale clustering. The new merger scheme offers improved or similar agreement with observational clustering measurements, over the redshift range 0 < z < 0.7. We find reasonable agreement with clustering measurements from GAMA, but find larger discrepancies for some stellar mass ranges and separation scales with respect to measurements from SDSS and VIPERS, depending on the GALFORM model used.

A major goal of the Atacama Large Millimeter/submillimeter Array (ALMA) is to make accurate images with resolutions of tens of milliarcseconds, which at submillimeter (submm) wavelengths requires baselines up to ~15 km. To develop and test this capability, a Long Baseline Campaign (LBC) was carried out from September to late November 2014, culminating in end-to-end observations, calibrations, and imaging of selected Science Verification (SV) targets. This paper presents an overview of the campaign and its main results, including an investigation of the short-term coherence properties and systematic phase errors over the long baselines at the ALMA site, a summary of the SV targets and observations, and recommendations for science observing strategies at long baselines. Deep ALMA images of the quasar 3C138 at 97 and 241 GHz are also compared to VLA 43 GHz results, demonstrating an agreement at a level of a few percent. As a result of the extensive program of LBC testing, the highly successful SV imaging at long baselines achieved angular resolutions as fine as 19 mas at ~350 GHz. Observing with ALMA on baselines of up to 15 km is now possible, and opens up new parameter space for submm astronomy.

Theoretical galaxy formation models are an established and powerful tool for interpreting the astrophysical significance of observational data, particularly galaxy surveys. Such models have been utilised with great success by optical surveys such as 2dFGRS and SDSS, but their application to radio surveys of cold gas in galaxies has been limited. In this chapter we describe recent developments in the modelling of the cold gas properties in the models, and how these developments are essential if they are to be applied to cold gas surveys of the kind that will be carried out with the SKA. By linking explicitly a galaxy's star formation rate to the abundance of molecular hydrogen in the galaxy rather than cold gas abundance, as was assumed previously, the latest models reproduce naturally many of the global atomic and molecular hydrogen properties of observed galaxies. We review some of the key results of the latest models and highlight areas where further developments are necessary. We discuss also how model predictions can be most accurately compared with observational data, what challenges we expect when creating synthetic galaxy surveys in the SKA era, and how the SKA can be used to test models of dark matter.

One of the key science drivers for the development of the SKA is to observe the neutral hydrogen, HI, in galaxies as a means to probe galaxy evolution across a range of environments over cosmic time. Over the past decade, much progress has been made in theoretical simulations and observations of HI in galaxies. However, recent HI surveys on both single dish radio telescopes and interferometers, while providing detailed information on global HI properties, the dark matter distribution in galaxies, as well as insight into the relationship between star formation and the interstellar medium, have been limited to the local universe. Ongoing and upcoming HI surveys on SKA pathfinder instruments will extend these measurements beyond the local universe to intermediate redshifts with long observing programmes. We present here an overview of the HI science which will be possible with the increased capabilities of the SKA and which will build upon the expected increase in knowledge of HI in and around galaxies obtained with the SKA pathfinder surveys. With the SKA1 the greatest improvement over our current measurements is the capability to image galaxies at reasonable linear resolution and good column density sensitivity to much higher redshifts (0.2 < z < 1.7). So one will not only be able to increase the number of detections to study the evolution of the HI mass function, but also have the sensitivity and resolution to study inflows and outflows to and from galaxies and the kinematics of the gas within and around galaxies as a function of environment and cosmic time out to previously unexplored depths. The increased sensitivity of SKA2 will allow us to image Milky Way-size galaxies out to redshifts of z=1 and will provide the data required for a comprehensive picture of the HI content of galaxies back to z~2 when the cosmic star formation rate density was at its peak.

We study the origin of the wide distribution of angles between the angular momenta of the stellar and gas components, $\alpha_{\rm G,S}$, in early-type galaxies (ETGs). We use the GALFORM model of galaxy formation, set in the $\Lambda$ cold dark matter framework, and coupled it with a Monte-Carlo simulation to follow the angular momenta flips driven by matter accretion onto haloes and galaxies. We consider a gas disk to be misaligned with respect to the stellar body if $\alpha_{\rm G,S}>30$~degrees. By assuming that the only sources of misaligments in galaxies are galaxy mergers, we place a lower limit of $2-5$ per cent on the fraction of ETGs with misaligned gas/stellar components. These low fractions are inconsistent with the observed value of $\approx 42\pm 6$ per cent in ATLAS$^{\rm 3D}$. In the more general case, in which smooth gas accretion in addition to galaxy mergers can drive misalignments, our calculation predicts that $\approx 46$ per cent of ETGs have $\alpha_{\rm G,S}>30$~degrees. In this calculation, we find correlations between $\alpha_{\rm G,S}$ and stellar mass, cold gas fraction and star formation rate, such that ETGs with high masses, low cold gas fractions and low star formation rates are more likely to display aligned cold gas and stellar components. We confirm these trends observationally for the first time using ATLAS$^{\rm 3D}$ data. We argue that the high fraction of misaligned gas discs observed in ETGs is mostly due to smooth gas accretion (e.g. cooling from the hot halo of galaxies) which takes place after most of the stellar mass of the galaxy is in place and comes misaligned with respect to the stellar component. Galaxies that have accreted most of their cold gas content prior to the time where most of the stellar mass was in place show aligned components.

Using Hubble Space Telescope (HST) ACS/WFC data we present the photometry and spatial distribution of resolved stellar populations in the outskirts of NGC 2915, a blue compact dwarf with an extended HI disc. These observations reveal an elliptical distribution of red giant branch stars, and a clumpy distribution of main-sequence stars that correlate with the HI gas distribution. We constrain the upper-end initial mass function (IMF) and determine the star formation law (SFL) in this field, using the observed main-sequence stars and an assumed constant star formation rate. Previously published H{\alpha} observations of the field, which show one faint HII region, are used to provide further constraints on the IMF. We find that the main-sequence luminosity function analysis alone results in a best-fitting IMF with a power-law slope {\alpha}=-2.85 and upper-mass limit M$_\rm{u}$ = 60 M$_\odot$. However, if we assume that all H{\alpha} emission is confined to HII regions then the upper-mass limit is restricted to M$_\rm{u}$ $\le$20 M$_\odot$. For the luminosity function fit to be correct we have to discount the H{\alpha} observations implying significant diffuse ionized gas or escaping ionizing photons. Combining the HST photometry with HI imaging we find the SFL has a power law index $N=1.53 \pm 0.21$. Applying these results to the entire outer HI disc indicates that it contributes 11--28% of the total recent star formation in NGC 2915, depending on whether the IMF is constant within the disc or varies from the centre to the outer region.

We aim to measure the average dust and molecular gas content of massive star-forming galaxies ($\rm > 3 \times 10^{10}\,M_\odot$) up to z=4 in the COSMOS field to determine if the intense star formation observed at high redshift is induced by major mergers or caused by large gas reservoirs. Firstly, we measured the evolution of the average spectral energy distributions as a function of redshift using a stacking analysis of Spitzer, Herschel, LABOCA, and AzTEC data for two samples of galaxies: normal star-forming objects and strong starbursts, as defined by their distance to the main sequence. We found that the mean intensity of the radiation field $< U >$ heating the dust (strongly correlated with dust temperature) increases with increasing redshift up to z$\sim$4 in main-sequence galaxies. We can reproduce this evolution with simple models that account for the decrease of the gas metallicity with redshift. No evolution of $< U >$ with redshift is found in strong starbursts. We then deduced the evolution of the molecular gas fraction (defined here as $\rm M_{\rm mol}/(M_{\rm mol}+M_\star)$) with redshift and found a similar, steeply increasing trend for both samples. At z$\sim$4, this fraction reaches $\sim$60%. The average position of the main-sequence galaxies is on the locus of the local, normal star-forming disks in the integrated Schmidt-Kennicutt diagram (star formation rate versus mass of molecular gas), suggesting that the bulk of the star formation up to z=4 is dominated by secular processes.

We study the atomic (HI) and molecular hydrogen (H2) contents of early-type galaxies (ETGs) and their gas sources using the GALFORM model of galaxy formation. This model uses a self-consistent calculation of the star formation rate (SFR), which depends on the H2 content of galaxies. We first present a new analysis of HIPASS and ATLAS3D surveys, with special emphasis on ETGs. The model predicts HI and H2 contents of ETGs in agreement with the observations from these surveys only if partial ram pressure stripping of the hot gas is included, showing that observations of neutral gas in `quenched' galaxies place stringent constraints on the treatment of the hot gas in satellites. We find that ~90% of ETGs at z=0 have neutral gas contents supplied by radiative cooling from their hot halos, 8% were supplied by gas accretion from minor mergers that took place in the last 1Gyr, while 2% were supplied by mass loss from old stars. The model predicts neutral gas fractions strongly decreasing with increasing bulge fraction. This is due to the impeded disk regeneration in ETGs, resulting from both active galactic nuclei feedback and environmental quenching by partial ram pressure stripping of the hot gas.

We exploit ALMA 870um (345GHz) observations of submillimetre sources in the Extended Chandra Deep Field South to investigate the far-infrared properties of high-redshift submillimetre galaxies (SMGs). Using the precisely located 870um ALMA positions of 99 SMGs, together with 24um and radio imaging of this field, we deblend the Herschel/SPIRE imaging of this region to extract their far-infrared fluxes and colours. The median photometric redshifts for ALMA LESS (ALESS) SMGs which are detected in at least two SPIRE bands increases with wavelength of the peak in their SEDs, with z=2.3+/-0.2, 2.5+/-0.3 and 3.5+/-0.5 for the 250, 350 and 500-um peakers respectively. We find that 34 ALESS SMGs do not have a >3-sigma counterpart at 250, 350 or 500-um. These galaxies have a median photometric redshift of z=3.3+/-0.5, which is higher than the full ALESS SMG sample; z=2.5+/-0.2. Using the photometric redshifts together with the 250-870um photometry, we estimate the far-infrared luminosities and characteristic dust temperature of each SMG. The median infrared luminosity of the S_870um>2mJy SMGs is L_IR=(3.0+/-0.3)x10^{12}Lo(SFR=300+/-30Mo/yr). At a fixed luminosity, the characteristic dust temperature of these high-redshift SMGs is 2-3K lower than comparably luminous galaxies at z=0, reflecting the more extended star formation occurring in these systems. By extrapolating the 870um number counts to S_ 870um=1mJy, we show that the contribution of S_870um>1mJy SMGs to the cosmic star formation budget is 20% of the total over the redshift range z~1-4. We derive a median dust mass for these SMGs of M_d=(3.6+/-0.3)x10^8Mo and by adopting an appropriate gas-to-dust ratio, we estimate an average molecular mass of M_H2=(4.2+/-0.4)x10^{10}Mo. Finally, we use our estimates of the H2 masses to show that SMGs with S_870um>1mJy contain ~10% of the z~2 volume-averaged H2 mass density at this epoch.

We study the contribution of galaxies with different properties to the global densities of star formation rate (SFR), atomic (HI) and molecular hydrogen (H2) as a function of redshift. We use the GALFORM model of galaxy formation, which is set in the LCDM framework. This model includes a self-consistent calculation of the SFR, which depends on the H2 content of galaxies. The predicted SFR density and how much of this is contributed by galaxies with different stellar masses and infrared luminosities are in agreement with observations. The model predicts a modest evolution of the HI density at z<3, which is also in agreement with the observations. The HI density is predicted to be always dominated by galaxies with SFR<1Msun/yr. This contrasts with the H2 density, which is predicted to be dominated by galaxies with SFR>10Msun/yr. Current high-redshift galaxy surveys are limited to detect carbon monoxide in galaxies with SFR>30Msun/yr, which in our model make up, at most, 20% of the H2 in the universe. In terms of stellar mass, the predicted H2 density is dominated by massive galaxies, Mstellar>10^10Msun, while the HI density is dominated by low mass galaxies, Mstellar<10^9Msun. In the context of upcoming neutral gas surveys, we suggest that the faint nature of the galaxies dominating the HI content of the Universe will hamper the identification of optical counterparts, while for H2, we expect follow up observations of molecular emission lines of already existing galaxy catalogues to be able to uncover the H2 density of the Universe.

We present a new release of the GALFORM semi-analytical model of galaxy formation and evolution, which exploits a Millennium Simulation-class N-body run performed with the WMAP7 cosmology. We use this new model to study the impact of the choice of stellar population synthesis (SPS) model on the predicted evolution of the galaxy luminosity function. The semi-analytical model is run using seven different SPS models. In each case we obtain the rest-frame luminosity function in the far-ultra-violet, optical and near-infrared (NIR) wavelength ranges. We find that both the predicted rest-frame ultra-violet and optical luminosity function are insensitive to the choice of SPS model. However, we find that the predicted evolution of the rest-frame NIR luminosity function depends strongly on the treatment of the thermally pulsating asymptotic giant branch (TP-AGB) stellar phase in the SPS models, with differences larger than a factor of 2 for model galaxies brighter than $M_{\rm AB}(K)-5$log$h<-22$ ($\sim$L$_*$ for $0\leq z\leq 1.5$). We have also explored the predicted number counts of galaxies, finding remarkable agreement between the results with different choices of SPS model, except when selecting galaxies with very red optical-NIR colours. The predicted number counts of these extremely red galaxies appear to be more affected by the treatment of star formation in disks than by the treatment of TP-AGB stars in the SPS models.

The HRS is a complete volume-limited sample of nearby objects including Virgo cluster and isolated objects. Using a recent compilation of HI and CO data we study the effects of the cluster on the molecular gas content of spiral galaxies. We first identify M* as the scaling variable that traces the total H2 mass of galaxies better. We show that, on average, HI-deficient galaxies are significantly offset from the M(H2) vs. M* relation for HI-normal galaxies. We use the M(H2) vs. M* scaling relation to define the H2-deficiency parameter. This parameter shows a weak and scattered relation with the HI-def, here taken as a proxy for galaxy interactions with the cluster environment. We also show that, as for the HI, the extent of the H2 disc decreases with increasing HI-deficiency. These results show that cluster galaxies have, on average, a lower H2 content than similar objects in the field. The slope of the H2-def vs. HI-def relation is less than 1, while the D(HI)/D(i) vs. HI-def relation is steeper than the D(CO)/D(i) vs. HI-def relation, thereby indicating that the H2 gas is removed less efficiently than the HI. This result can be understood if the HI is distributed on a flat disc more extended than the stellar disc, thus less anchored to the gravitational potential well of the galaxy than the H2. There is a clear trend between the NUV-i colour and H2-def, which suggests that H2 removal quenches the activity of star formation. This causes galaxies migrate from the blue cloud to the green valley and, eventually, to the red sequence. The total gas-consumption timescale of gas deficient cluster galaxies is comparable to that of isolated systems, and is significantly larger than the typical timescale for total gas removal in a ram pressure stripping process, thus suggesting that ram pressure, rather than starvation, is the dominant process driving the evolution of these cluster galaxies.

We study the properties of the cold gas component of the interstellar medium of the Herschel Reference Survey, a complete volume-limited (15<D<25 Mpc), K-band-selected sample of galaxies spanning a wide range in morphological type (from E to Im) and stellar mass (10^9<M*<10^11 Mo). The multifrequency data in our hands are used to trace the molecular gas mass distribution and the main scaling relations of the sample, which put strong constraints on galaxy formation simulations. We extend the main scaling relations concerning the total and the molecular gas component determined for massive galaxies (M* > 10^10 Mo) from the COLD GASS survey down to stellar masses M* ~ 10^9 Mo. As scaling variables we use M*, the stellar surface density mu*, the specific star formation rate SSFR, and the metallicity of the target galaxies. By comparing molecular gas masses determined using a constant or a luminosity dependent conversion factor, we estimate the robustness of these scaling relations on the very uncertain assumptions used to transform CO line intensities into molecular gas masses. The molecular gas distribution of a K-band-selected sample is different from that of a far-infrared-selected sample since it includes a significantly smaller number of objects with M(H2) < 6 10^9 Mo. In spiral galaxies the molecular gas phase is only 25-30% of the atomic gas. The analysis also indicates that the slope of the main scaling relations depends on the adopted conversion factor. Among the sampled relations, all those concerning M(gas)/M* are statistically significant and show little variation with X_CO. We observe a significant correlation between M(H2)/M* and SSFR, M(H2)/M(HI) and mu*, M(H2)/M(HI), and 12+log(O/H) regardless of the adopted X_CO. The total and molecular gas consumption timescales are anticorrelated with the SSFR.

We present a dynamical model of supernova feedback which follows the evolution of pressurised bubbles driven by supernovae in a multi-phase interstellar medium (ISM). The bubbles are followed until the point of break-out into the halo, starting from an initial adiabatic phase to a radiative phase. We show that a key property which sets the fate of bubbles in the ISM is the gas surface density, through the work done by the expansion of bubbles and its role in setting the gas scaleheight. The multi-phase description of the ISM is essential, and neglecting it leads to order of magnitude differences in the predicted outflow rates. We compare our predicted mass loading and outflow velocities to observations of local and high-redshift galaxies and find good agreement over a wide range of stellar masses and velocities. With the aim of analysing the dependence of the mass loading of the outflow, beta (i.e. the ratio between the outflow and star formation rates), on galaxy properties, we embed our model in the galaxy formation simulation, GALFORM, set in the LCDM framework. We find that a dependence of beta solely on the circular velocity, as is widely assumed in the literature, is actually a poor description of the outflow rate, as large variations with redshift and galaxy properties are obtained. Moreover, we find that below a circular velocity of 80km/s the mass loading saturates. A more fundamental relation is that between beta and the gas scaleheight of the disk, hg, and the gas fraction, fgas, as beta hg^(1.1) fgas^(0.4), or the gas surface density, \Sigma_g, and the gas fraction, as beta \Sigma_g^(-0.6) fgas^(0.8). We find that using the new mass loading model leads to a shallower faint-end slope in the predicted optical and near-IR galaxy luminosity functions.

We introduce a method for constructing end-to-end mock galaxy catalogues using a semi-analytical model of galaxy formation, applied to the halo merger trees extracted from a cosmological N-body simulation. The mocks that we construct are lightcone catalogues, in which a galaxy is placed according to the epoch at which it first enters the past lightcone of the observer, and incorporate the evolution of galaxy properties with cosmic time. We determine the position between the snapshot outputs at which a galaxy enters the observer's lightcone by interpolation. As an application, we consider the effectiveness of the BzK colour selection technique, which was designed to isolate galaxies in the redshift interval 1.4<z<2.5. The mock catalogue is in reasonable agreement with the observed number counts of all BzK galaxies, as well as with the observed counts of the subsample of BzKs that are star-forming galaxies. We predict that over 75 per cent of the model galaxies with K_{AB}<=23, and 1.4<z<2.5, are selected by the BzK technique. Interloper galaxies, outside the intended redshift range, are predicted to dominate bright samples of BzK galaxies (i.e. with K_{AB}<=21). Fainter K-band cuts are necessary to reduce the predicted interloper fraction. We also show that shallow B-band photometry can lead to confusion in classifying BzK galaxies as being star-forming or passively evolving. Overall, we conclude that the BzK colour selection technique is capable of providing a sample of galaxies that is representative of the 1.4<z<2.5 galaxy population.

We investigate the impact of feedback - from supernovae (SNe), active galactic nuclei (AGN) and a photo-ionizing background at high redshifts - on the neutral atomic hydrogen (HI) mass function, the $b_{\rm J}$ band luminosity function, and the spatial clustering of these galaxies at $z$=0. We use a version of the semi-analytical galaxy formation model GALFORM that calculates self-consistently the amount of HI in a galaxy as a function of cosmic time and links its star formation rate to its mass of molecular hydrogen (H$_2$). We find that a systematic increase or decrease in the strength of SNe feedback leads to a systematic decrease or increase in the amplitudes of the luminosity and HI mass functions, but has little influence on their overall shapes. Varying the strength of AGN feedback influences only the numbers of the brightest or most HI massive galaxies, while the impact of varying the strength of photo-ionization feedback is restricted to changing the numbers of the faintest or least HI massive galaxies.Our results suggest that the HI mass function is a more sensitive probe of the consequences of cosmological reionization for galaxy formation than the luminosity function. We find that increasing the strength of any of the modes of feedback acts to weaken the clustering strength of galaxies, regardless of their HI-richness. In contrast, weaker AGN feedback has little effect on the clustering strength whereas weaker SNe feedback increases the clustering strength of HI-poor galaxies more strongly than HI-rich galaxies. These results indicate that forthcoming HI surveys on next generation radio telescopes such as the Square Kilometre Array and its pathfinders will be exploited most fruitfully as part of multiwavelength survey campaigns.

The latest observations of molecular gas and the atomic hydrogen content of local and high-redshift galaxies, coupled with how these correlate with star formation activity, have revolutionized our ideas about how to model star formation in a galactic context. A successful theory of galaxy formation has to explain some key facts: (i) high-redshift galaxies have higher molecular gas fractions and star formation rates than local galaxies, (ii) scaling relations show that the atomic-to-stellar mass ratio decreases with stellar mass in the local Universe, and (iii) the global abundance of atomic hydrogen evolves very weakly with time. We review how modern cosmological simulations of galaxy formation attempt to put these pieces together and highlight how approaches simultaneously solving dark matter and gas physics, and approaches first solving the dark matter N-body problem and then dealing with gas physics using semi-analytic models, differ and complement each other. We review the observable predictions, what we think we have learned so far and what still needs to be done in the simulations to allow robust testing by the new observations expected from telescopes such as ALMA, PdBI, LMT, JVLA, ASKAP, MeerKAT, SKA.

We combine the galaxy formation model GALFORM with the Photon Dominated Region code UCL_PDR to study the emission from the rotational transitions of 12CO (CO) in galaxies from z=0 to z=6 in the Lambda CDM framework. GALFORM is used to predict the molecular (H2) and atomic hydrogen (HI) gas contents of galaxies using the pressure-based empirical star formation relation of Blitz & Rosolowsky. From the predicted H2 mass and the conditions in the interstellar medium, we estimate the CO emission in the rotational transitions 1-0 to 10-9 by applying the UCL_PDR model to each galaxy. We find that deviations from the Milky-Way CO-H2 conversion factor come mainly from variations in metallicity, and in the average gas and star formation rate surface densities. In the local universe, the model predicts a CO(1-0) luminosity function (LF), CO-to-total infrared (IR) luminosity ratios for multiple CO lines and a CO spectral line energy distribution (SLED) which are in good agreement with observations of luminous and ultra-luminous IR galaxies. At high redshifts, the predicted CO SLED of the brightest IR galaxies reproduces the shape and normalization of the observed CO SLED. The model predicts little evolution in the CO-to-IR luminosity ratio for different CO transitions, in good agreement with observations up to z~5. We use this new hybrid model to explore the potential of using colour selected samples of high-redshift star-forming galaxies to characterise the evolution of the cold gas mass in galaxies through observations with the Atacama Large Millimeter Array.

Using a sample of BzK-selected galaxies at z~2 identified from the CFHT/WIRCAM near-infrared survey of GOODS-North, we discuss the relation between star formation rate (SFR), specific star formation rate (SSFR), and stellar mass (M_{*}), and the clustering of galaxies as a function of these parameters. For star-forming galaxies (sBzKs), the UV-based SFR, corrected for extinction, scales with the stellar mass as SFR ~ M_{*}^{alpha} with alpha = 0.74+/-0.20 down to M_{*} ~ 10^{9} M_{solar}, indicating a weak dependence on the stellar mass of the SSFR. We also measure the angular correlation function and hence infer the correlation length for sBzK galaxies as a function of M_{*}, SFR, and SSFR, as well as K-band apparent magnitude. We show that passive galaxies (pBzKs) are more strongly clustered than sBzK galaxies at a given stellar mass, mirroring the color-density relation seen at lower redshifts. We also find that the correlation length of sBzK galaxies ranges from 4 to 20 h^{-1}Mpc, being a strong function of M_{K}, M_{*}, and SFR. On the other hand, the clustering dependence on SSFR changes abruptly at 2x10^{-9} yr^{-1}, which is the typical value for "main sequence" star-forming galaxies at z~2. We show that the correlation length reaches a minimum at this characteristic value, and is larger for galaxies with both smaller and larger SSFRs; a dichotomy that is only marginally implied from the predictions of the semi-analytical models. Our results suggest that there are two types of environmental effects at work at z~2. Stronger clustering for relatively quiescent galaxies implies that the environment has started to play a role in quenching star formation. At the same time, stronger clustering for galaxies with elevated SSFRs ("starbursts") might be attributed to an increased efficiency for galaxy interactions and mergers in dense environments.

We study the evolution of the cold gas content of galaxies by splitting the interstellar medium into its atomic and molecular hydrogen components, using the galaxy formation model GALFORM in the LCDM framework. We calculate the molecular-to-atomic hydrogen mass ratio, H2/HI, in each galaxy using two different approaches; the pressure-based empirical relation of Blitz & Rosolowsky and the theoretical model of Krumholz, McKeee & Tumlinson, and apply them to consistently calculate the star formation rates of galaxies. We find that the model based on the Blitz & Rosolowsky law predicts an HI mass function, CO(1-0) luminosity function, correlations between the H2/HI ratio and stellar and cold gas mass, and infrared-CO luminosity relation in good agreement with local and high redshift observations. The HI mass function evolves weakly with redshift, with the number density of high mass galaxies decreasing with increasing redshift. In the case of the H2 mass function, the number density of massive galaxies increases strongly from z=0 to z=2, followed by weak evolution up to z=4. We also find that the H2/HI ratio of galaxies is strongly dependent on stellar and cold gas mass, and also on redshift. The slopes of the correlations between H2/HI and stellar and cold gas mass hardly evolve, but the normalisation increases by up to two orders of magnitude from z=0-8. The strong evolution in the H2 mass function and the H2/HI ratio is primarily due to the evolution in the sizes of galaxies and secondarily, in the gas fractions. The predicted cosmic density evolution of HI agrees with the observed evolution inferred from DLAs, and is dominated by low/intermediate mass halos. We find that previous theoretical studies have largely overestimated the redshift evolution of the global H2/HI ratio due to limited resolution. We predict a maximum of rho_H2/rho_HI~1.2 at z~3.5.

We investigate the consequences of applying different star formation laws in the galaxy formation model GALFORM. Three broad star formation laws are implemented: the empirical relations of Kennicutt and Schmidt and Blitz & Rosolowsky and the theoretical model of Krumholz, McKee & Tumlinson. These laws have no free parameters once calibrated against observations of the star formation rate (SFR) and gas surface density in nearby galaxies. We start from published models, and investigate which observables are sensitive to a change in the star formation law, without altering any other model parameters. We show that changing the star formation law (i) does not significantly affect either the star formation history of the universe or the galaxy luminosity functions in the optical and near-IR, due to an effective balance between the quiescent and burst star formation modes; (ii) greatly affects the cold gas contents of galaxies; (iii) changes the location of galaxies in the SFR versus stellar mass plane, so that a second sequence of "passive" galaxies arises, in addition to the known "active" sequence. We show that this plane can be used to discriminate between the star formation laws.

We present IRAM Plateau de Bure interferometric detections of CO(1-0) emission from a 24um-selected sample of star-forming galaxies at z=0.4. The galaxies have PAH 7.7um-derived star formation rates of SFR~30-60 M_Sun/yr and stellar masses M*~10^{11} M_Sun. The CO(1-0) luminosities of the galaxies imply that the disks still contain a large reservoir of molecular gas, contributing ~20% of the baryonic mass, but have star-formation 'efficiencies' similar to local quiescent disks and gas-dominated disks at z~1.5-2. We reveal evidence that the average molecular gas fraction has undergone strong evolution since z~2, with f_gas ~ (1+z)^{2 +/- 0.5}. The evolution of f_gas encodes fundamental information about the relative depletion/replenishment of molecular fuel in galaxies, and is expected to be a strong function of halo mass. We show that the latest predictions for the evolution of the molecular gas fraction in semi-analytic models of galaxy formation within a LCDM Universe are supported by these new observations.

We determine the intrinsic shapes and orientations of 27,450 type I and II active galactic nucleus (AGN) galaxies in the spectroscopic sample of the SDSS DR7, by studying the distribution of projected axis ratios of AGN hosts. Our aim is to study possible alignments between the AGN and host galaxy systems (e.g. the accretion disc and the galaxy angular momentum) and the effect of dust obscuration geometry on the AGN type. We define control samples of non-AGN galaxies that mimic the morphology, colour, luminosity and concentration distributions of the AGN population, taking into account the effects of dust extinction and reddening. Assuming that AGN galaxies have the same underlying three-dimensional shape distribution as their corresponding control samples, we find that the spiral and elliptical type I AGN populations are strongly skewed toward face-on galaxies, while ellipticals and spirals type II AGN are skewed toward edge-on orientations. These findings rule out random orientations for AGN hosts at high confidence for type I spirals (delta chi^2~230) and type II ellipticals (delta chi^2~15), while the signal for type I ellipticals and type II spirals is weaker (delta chi^2~3 and delta chi^2~6, respectively). We obtain a much stronger tendency for the type II spirals to be edge-on when just high [OIII] equivalent width (EW) AGN are considered, suggesting that >20% of low [OIII] EW edge-on type II AGN may be missing from the optical sample. Galactic dust absorption of the broad-line region alone cannot explain the observed inclination angle and projected axis ratio distributions of type I and II Seyfert types, implying that obscuration by a small-scale circumnuclear torus is necessary. These results favour a scenario in which the angular momentum of the material which feeds the black hole retains a memory of its original gas source at least to some small, non-negligible degree.

We use a coupled model of the formation and evolution of galaxies and black holes (BH) to study the evolution of active galactic nuclei (AGN) in a cold dark matter universe. The model predicts the BH mass, spin and mass accretion history. BH mass grows via accretion triggered by discs becoming dynamically unstable or galaxy mergers (called the starburst mode) and accretion from quasi-hydrostatic hot gas haloes (called the hot-halo mode). By taking into account AGN obscuration, we obtain a very good fit to the observed luminosity functions (LF) of AGN (optical, soft and hard X-ray, and bolometric) for a wide range of redshifts (0<z<6). The model predicts a hierarchical build up of BH mass, with the typical mass of actively growing BHs increasing with decreasing redshift. Remarkably, despite this, we find downsizing in the AGN population, in terms of the differential growth with redshift of the space density of faint and bright AGN. This arises naturally from the interplay between the starburst and hot-halo accretion modes. The faint end of the LF is dominated by massive BHs experiencing quiescent accretion via a thick disc, primarily during the hot-halo mode. The bright end of the LF, on the other hand, is dominated by AGN which host BHs accreting close to or in excess of the Eddington limit during the starburst mode. The model predicts that the comoving space density of AGN peaks at z~3, similar to the star formation history. However, when taking into account obscuration, the space density of faint AGN peaks at lower redshift (z<2) than that of bright AGN (z~2-3). This implies that the cosmic evolution of AGN is shaped in part by obscuration.

We develop a new numerical approach to describe the action of ram pressure stripping (RPS) within a semi-analytic model of galaxy formation and evolution which works in combination with non-radiative hydrodynamical simulations of galaxy clusters. The new feature in our method is the use of the gas particles to obtain the kinematical and thermodynamical properties of the intragroup and intracluster medium (ICM). This allows a self-consistent estimation of the RPS experienced by satellite galaxies. We find that the ram pressure (RP) in the central regions of clusters increases approximately one order of magnitude between z = 1 and 0, consistent with the increase in the density of the ICM. The mean RP experienced by galaxies within the virial radius increases with decreasing redshift. In clusters with virial masses M_vir ~10^15 h^-1 M_Sun, over 50 per cent of satellite galaxies have experienced RP ~10^(-11) h^2 dyn cm^-2 or higher for z <= 0.5. In smaller clusters (M_vir ~10^14 h^-1 M_Sun) the mean RP are approximately one order of magnitude lower at all redshifts. RPS has a strong effect on the cold gas content of galaxies for all cluster masses. At z = 0, over 70 per cent of satellite galaxies within the virial radius are completely depleted of cold gas. For the more massive clusters the fraction of depleted galaxies is already established at z ~ 1, whereas for the smaller clusters this fraction increases appreciably between z = 1 and 0. This indicates that the rate at which the cold gas is stripped depends on the virial mass of the host cluster. Compared to our new approach, the use of an analytic profile to describe the ICM results in an overestimation of the RP larger than 50 per cent for z >= 0.5.

We study the properties of building blocks (BBs, i.e. accreted satellites) and surviving satellites of present-day galaxies using the SAG semi-analytic model of galaxy formation in the context of a concordance Lambda Cold Dark Matter (LCDM) cosmology. We consider large numbers of DM halo merger trees spanning a wide range of masses (~1x10^10 - 2.14x10^15 Msun). We find higher metallicities for BBs with respect to surviving satellites, an effect produced by the same processes behind the build-up of the mass-metallicity relation. We prove that these metallicity differences arise from the higher peak height in the density fluctuation field occupied by BBs and central galaxies which have collapsed into a single object earlier than surviving satellites. BBs start to form stars earlier, during the peak of the merger activity in LCDM, and build-up half of their final stellar mass (measured at the moment of disruption) up to four times faster than surviving satellites. Surviving satellites keep increasing their stellar masses rather quiescently down to z~1. The difference between the metallicities of satellites, BBs and central galaxies depends on the host DM halo mass, in a way that can be used as a further test for the concordance cosmology.

We use a combination of a cosmological N-body simulation of the concordance Lambda cold dark matter (LCDM) paradigm and a semi-analytic model of galaxy formation to investigate the spin development of central supermassive black holes (BHs) and its relation to the BH host galaxy properties. In order to compute BH spins, we use the alpha-model of Shakura & Sunyaev and consider the King et al. warped disc alignment criterion. The orientation of the accretion disc is inferred from the angular momentum of the source of accreted material, which bears a close relationship to the large-scale structure in the simulation. We find that the final BH spin depends almost exclusively on the accretion history and only weakly on the warped disc alignment. The main mechanisms of BH spin-up are found to be gas cooling processes and disc instabilities, a result that is only partially compatible with Monte-Carlo models where the main spin-up mechanisms are major mergers and disc instabilities; the latter results are reproduced when implementing randomly oriented accretion discs in our model. Regarding the BH population, we find that more massive BHs, which are hosted by massive ellipticals, have higher spin values than less-massive BHs, hosted by spiral galaxies. We analyse whether gas accretion rates and BH spins can be used as tracers of the radio loudness of active galactic nuclei (AGN). We find that the current observational indications of an increasing trend of radio-loud AGN fractions with stellar and BH mass can be easily obtained when placing lower limits on the BH spin, with a minimum influence from limits on the accretion rates; a model with random accretion disc orientations is unable to reproduce this trend. (ABRIDGED)

We study the effects of Active Galactic Nuclei (AGN) feedback on the formation and evolution of galaxies in a semi-analytic model of galaxy formation. This model is an improved version of the one described by Cora (2006), which now considers the growth of black holes (BHs) as driven by (i) gas accretion during merger-driven starbursts and mergers with other BHs, (ii) accretion during starbursts triggered by disc instabilities, and (iii) accretion of gas cooled from quasi-hydrostatic hot gas haloes. It is assumed that feedback from AGN operates in the later case. The model has been calibrated in order to reproduce observational correlations between BH mass and mass, velocity dispersion, and absolute magnitudes of the galaxy bulge. AGN feedback has a strong impact on reducing or even suppressing gas cooling, an effect that becomes important at lower redshifts. This phenomenon helps to reproduce the observed galaxy luminosity function (LF) in the optical and near IR bands at z=0, and the cosmic star formation rate and stellar mass functions over a wide redshift range (0<z<5). It also allows to have a population of massive galaxies already in place at z>1, which are mostly early-type and have older and redder stellar populations than lower mass galaxies, reproducing the observed bimodality in the galaxy colour distribution, and the morphological fractions. The evolution of the optical QSO LF is also reproduced, provided that the presence of a significant fraction of obscured QSOs is assumed. We explore the effects of AGN feedback during starbursts and new recent prescriptions for dynamical friction time-scales. (ABRIDGED)