We show that correlations between the phases of the galaxy density field in redshift space provide additional information about the growth rate of large-scale structure that is complementary to the power spectrum multipoles. In particular, we consider the multipoles of the line correlation function (LCF), which correlates phases between three collinear points, and use the Fisher forecasting method to show that the LCF multipoles can break the degeneracy between the measurement of the growth rate of structure $f$ and the amplitude of perturbations $\sigma_8$ that is present in the power spectrum multipoles at large scales. This leads to an improvement in the measurement of $f$ and $\sigma_8$ by up to 220 per cent for $k_{\rm max} = 0.15 \, h\mathrm{Mpc}^{-1}$ and up to 50 per cent for $k_{\rm max} = 0.30 \, h\mathrm{Mpc}^{-1}$ at redshift $z=0.25$, with respect to power spectrum measurements alone for the upcoming generation of galaxy surveys like DESI and Euclid. The average improvements in the constraints on $f$ and $\sigma_8$ for $k_{\rm max} = 0.15 \, h\mathrm{Mpc}^{-1}$ are $\sim 90$ per cent for the DESI BGS sample with mean redshift $\overline{z}=0.25$, $\sim 40$ per cent for the DESI ELG sample with $\overline{z}=1.25$, and $\sim 40$ per cent for the Euclid H$\alpha$ galaxies with $\overline{z}=1.3$. For $k_{\rm max} = 0.30 \, h\mathrm{Mpc}^{-1}$, the average improvements are $\sim 40$ per cent for the DESI BGS sample and $\sim 20$ per cent for both the DESI ELG and Euclid H$\alpha$ galaxies.

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

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.

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.

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.

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 present an analysis of the gas dynamics of star-forming galaxies at z~1.5 using data from the KMOS Galaxy Evolution Survey (KGES). We quantify the morphology of the galaxies using $HST$ CANDELS imaging parametrically and non-parametrically. We combine the H$\alpha$ dynamics from KMOS with the high-resolution imaging to derive the relation between stellar mass (M$_{*}$) and stellar specific angular momentum (j$_{*}$). We show that high-redshift star-forming galaxies at z~1.5 follow a power-law trend in specific stellar angular momentum with stellar mass similar to that of local late-type galaxies of the form j$_*$$\propto$M$_*^{0.53 \pm 0.10}$. The highest specific angular momentum galaxies are mostly disc-like, although generally, both peculiar morphologies and disc-like systems are found across the sequence of specific angular momentum at a fixed stellar mass. We explore the scatter within the j$_{*}$-M$_{*}$ plane and its correlation with both the integrated dynamical properties of a galaxy (e.g. velocity dispersion, Toomre Q$_{\rm g}$, H$\alpha$ star formation rate surface density $\Sigma_{\rm SFR}$) and its parameterised rest-frame UV/optical morphology (e.g. S\'ersic index, bulge to total ratio, Clumpiness, Asymmetry and Concentration). We establish that the position in the j$_{*}$-M$_{*}$ plane is strongly correlated with the star-formation surface density and the Clumpiness of the stellar light distribution. Galaxies with peculiar rest-frame UV/optical morphologies have comparable specific angular momentum to disc-dominated galaxies of the same stellar mass, but are clumpier and have higher star-formation rate surface densities. We propose that the peculiar morphologies in high--redshift systems are driven by higher star formation rate surface densities and higher gas fractions leading to a more clumpy inter-stellar medium.

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.

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

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.

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.

Redshift-space distortions are a sensitive probe of the growth of large-scale structure. In the linear regime, redshift-space distortions are fully described by the multipoles of the two-point correlation function. In the nonlinear regime, however, higher-order statistics are needed to capture the full information of the galaxy density field. In this paper, we show that the redshift-space line correlation function--which is a measure of Fourier phase correlations--is sensitive to the nonlinear growth of the density and velocity fields and to the nonlinear mapping between real and redshift space. We expand the line correlation function in multipoles, and we show that almost all of the information is encoded in the monopole, quadrupole, and hexadecapole. We argue that these multipoles are highly complementary to the multipoles of the two-point correlation function: first, because they are directly sensitive to the difference between the density and the velocity coupling kernels, which is a purely nonlinear quantity; and second, because the multipoles are proportional to different combinations of $f$ and $\sigma_8$. Measured in conjunction with the two-point correlation function and the bispectrum, the multipoles of the line correlation function could therefore allow us to disentangle efficiently these two quantities and to test modified theories of gravity.

We present adaptive optics assisted integral field spectroscopy of 34 star-forming galaxies at $z$ = 0.8-3.3 selected from the HiZELS narrow-band survey. We measure the kinematics of the ionised interstellar medium on $\sim$1 kpc scales, and show that the galaxies are turbulent, with a median ratio of rotational to dispersion support of $v$/$\sigma$=0.82$\pm$0.13. We combine the dynamics with high-resolution rest-frame optical imaging and extract emission line rotation curves. We show that high-redshift star-forming galaxies follow a similar power-law trend in specific angular momentum with stellar mass as that of local late type galaxies. We exploit the high resolution of our data and examine the radial distribution of angular momentum within each galaxy by constructing total angular momentum profiles. Although the stellar mass of a typical star-forming galaxy is expected to grow by a factor $\sim$8 in the $\sim$5 Gyrs between $z$$\sim$3.3 and $z$$\sim$0.8, we show that the internal distribution of angular momentum becomes less centrally concentrated in this period i.e the angular momentum grows outwards. To interpret our observations, we exploit the EAGLE simulation and trace the angular momentum evolution of star forming galaxies from $z$$\sim$3 to $z$$\sim$0, identifying a similar trend of decreasing angular momentum concentration. This change is attributed to a combination of gas accretion in the outer disk, and feedback that preferentially arises from the central regions of the galaxy. We discuss how the combination of the growing bulge and angular momentum stabilises the disk and gives rise to the Hubble sequence.

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

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.

We explore the scaling between the size of star-forming clumps and rotational support in massively star-forming galactic disks. The analysis relies on simulations of a clumpy galaxy at $z=2$ and the observed DYNAMO sample of rare clumpy analogs at $z\approx0.1$ to test a predictive clump size scaling proposed by \citet{Fisher2017ApJ...839L...5F} in the context of the Violent Disk Instability (VDI) theory. We here determine the clump sizes using a recently presented 2-point estimator, which is robust against resolution/noise effects, hierarchical clump substructure, clump-clump overlap and other galactic substructure. After verifying Fisher's clump scaling relation for the DYNAMO observations, we explore whether this relation remains characteristic of the VDI theory, even if realistic physical processes, such as local asymetries and stellar feedback, are included in the model. To this end, we rely on hydrodynamic zoom-simulations of a Milky Way-mass galaxy with four different feedback prescriptions. We find that, during its marginally stable epoch at $z=2$, this mock galaxy falls on the clump scaling relation, although its position on this relation depends on the feedback model. This finding implies that Toomre-like stability considerations approximately apply to large ($\sim\rm kpc$) instabilities in marginally stable turbulent disks, irrespective of the feedback model, but also emphasizes that the global clump distribution of a turbulent disk depends strongly on feedback.

We report interferometric observations tuned to the redshifted neutral hydrogen (HI) 21cm emission line in three strongly lensed galaxies at $z \sim 0.4$ with the Giant Metrewave Radio Telescope (GMRT). One galaxy spectrum (J1106+5228 at z=0.407) shows evidence of a marginal detection with an integrated signal-to-noise ratio of 3.8, which, if confirmed by follow-up observations, would represent the first strongly lensed and most distant individual galaxy detected in HI emission. Two steps are performed to transcribe the lensed integrated flux measurements into HI mass measurements for all three target galaxies. First, we calculate the HI magnification factor $\mu$ by applying general relativistic ray-tracing to a physical model of the source-lens system. The HI magnification generally differs from the optical magnification and depends largely on the intrinsic HI mass $M_{\rm HI}$ due to the HI mass-size relation. Second, we employ a Bayesian formalism to convert the integrated flux, amplified by the $M_{\rm HI}$-dependent magnification factor $\mu$, into a probability density for $M_{\rm HI}$, accounting for the asymmetric uncertainty due to the declining HI mass function (Eddington bias). In this way, we determine a value of $\log_{\rm 10} (M_{\rm HI}/M_\odot) = 10.2^{+0.3}_{-0.7}$ for J1106+5228, consistent with the estimate of $9.4\pm0.3$ from the optical properties of this galaxy. The HI mass of the other two sources are consistent with zero within a 95 per cent confidence interval however we still provide upper limits for both sources and a $1\sigma$ lower limit for J1250-0135 using the same formalism.

U-band observations with the LSST have yet to be fully optimized in cadence. The straw man survey design is a simple coverage of the medium-deep-fast survey. Here we argue that deep coverage of the four deep drilling fields (XMM-LSS, ECDFS, ELAIS-S1 and COSMOS) has a much higher scientific return, given that these are also the target of the Southern Hemisphere's Square Kilometer Array Pathfinder, the MeerKAT specifically, deep radio observations.

By analysing a sample of galaxies selected from the HI Parkes All Sky Survey (HIPASS) to contain more than 2.5 times their expected HI content based on their optical properties, we investigate what drives these HI eXtreme (HIX) galaxies to be so HI-rich. We model the HI kinematics with the Tilted Ring Fitting Code TiRiFiC and compare the observed HIX galaxies to a control sample of galaxies from HIPASS as well as simulated galaxies built with the semi-analytic model Dark Sage. We find that (1) HI discs in HIX galaxies are more likely to be warped and more likely to host HI arms and tails than in the control galaxies, (2) the average HI and average stellar column density of HIX galaxies is comparable to the control sample, (3) HIX galaxies have higher HI and baryonic specific angular momenta than control galaxies, (4) most HIX galaxies live in higher-spin haloes than most control galaxies. These results suggest that HIX galaxies are HI-rich because they can support more HI against gravitational instability due to their high specific angular momentum. The majority of the HIX galaxies inherits their high specific angular momentum from their halo. The HI content of HIX galaxies might be further increased by gas-rich minor mergers. This paper is based on data obtained with the Australia Telescope Compact Array (ATCA) through the large program C 2705.

The neutral atomic hydrogen (HI) content of spiral galaxies has been observed to vary with environment, with more HI-deficient spirals residing in high density environments. This can be attributed to environmental effects such as ram pressure stripping and tidal interactions, which remove HI from the discs of galaxies. However, some spirals in low-density environments have also been observed to have relatively low HI mass fractions. The low densities of the Intergalactic Medium and lack of nearby galaxies in such environments make ram pressure stripping and tidal interactions unlikely candidates of gas removal. What then could be making these spirals HI deficient? Obreschkow et al. introduced a parameter-free model for the neutral atomic gas fraction ($f_{atm}$), in a symmetric equilibrium disc as a function of the global atomic stability parameter ($q$), which depends on specific angular momentum. In order to examine if this model accounts for HI-deficient galaxies in low-density environments, we have used the $M_{HI} ~-$ M$_{R}$ scaling relation to select six HI-deficient spiral galaxies and observed them with the ATCA. By measuring their $f_{atm}$ and $q$ values we find that the galaxies owe their observed HI deficiencies to low specific angular momenta. Additionally, we also find that the central HI hole sizes of our sample galaxies are related to their $q$ values, following the prediction of Obreschkow et al. This result brings to light the importance of angular momentum in understanding the physics of the interstellar medium in the discs of galaxies and consequently their evolution.

In this paper we compare the molecular gas depletion times and mid-plane hydrostatic pressure in turbulent, star forming disk galaxies to internal properties of these galaxies. For this analysis we use 17 galaxies from the DYNAMO sample of nearby ($z\sim0.1$) turbulent disks. We find a strong correlation, such that galaxies with lower molecular gas depletion time (t_{dep}) have higher gas velocity dispersion ($\sigma$). Within the scatter of our data, our observations are consistent with the prediction that depletion time is inversely proportionate to velocity dispersion made in theories of feedback-regulated star formation. We also show a strong, single power-law correlation between mid-plane pressure (P) and star formation rate surface density (Sigma_{SFR}), which extends for 6 orders of magnitude in pressure. Disk galaxies with lower pressure are found to be roughly in agreement with theoretical predictions. However, in galaxies with high pressure we find P/Sigma_{SFR} values that are significantly larger than theoretical predictions. Our observations could be explained with any of the following: (1) the correlation of \Sigma_{SFR}-P is significantly sub-linear; (2) the momentum injected from star formation feedback (so-called p*/m*) is not a single, universal value; or (3) alternate sources of pressure support are important in gas-rich disk galaxies. Finally using published survey results, we find that our results are consistent with the cosmic evolution of t_{dep}(z) and velocity dispersion. Our interpretation of these results is that the cosmic evolution of t_{dep} may be regulated not just by the supply of gas, but also the internal regulation of star formation via feedback.

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.

We construct Tully-Fisher relations (TFRs), from large samples of galaxies with spatially-resolved H$\alpha$ emission maps from the K-band Multi-Object Spectrograph (KMOS) Redshift One Spectroscopic Survey (KROSS) at $z\approx1$. We compare these to data from the Sydney-Australian-Astronomical-Observatory Multi-object Integral-Field Spectrograph (SAMI) Galaxy Survey at $z\approx0$. We stringently match the data quality of the latter to the former, and apply identical analysis methods and sub-sample selection criteria to both to conduct a direct comparison of the absolute $K$-band magnitude and stellar mass TFRs at $z\approx1$ and $z\approx0$. We find that matching the quality of the SAMI data to that of KROSS results in TFRs that differ significantly in slope, zero-point and (sometimes) scatter in comparison to the corresponding original SAMI relations. These differences are in every case as large or larger than the differences between the KROSS $z\approx1$ and matched SAMI $z\approx0$ relations. Accounting for these differences, we compare the TFRs at $z\approx1$ and $z\approx0$. For disk-like, star-forming galaxies we find no significant difference in the TFR zero-points between the two epochs. This suggests the growth of stellar mass and dark matter in these types of galaxies is intimately linked over this $\approx8$ Gyr period.

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.

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.

Most statistical inference from cosmic large-scale structure relies on two-point statistics, i.e.\ on the galaxy-galaxy correlation function (2PCF) or the power spectrum. These statistics capture the full information encoded in the Fourier amplitudes of the galaxy density field but do not describe the Fourier phases of the field. Here, we quantify the information contained in the line correlation function (LCF), a three-point Fourier phase correlation function. Using cosmological simulations, we estimate the Fisher information (at redshift $z=0$) of the 2PCF, LCF and their combination, regarding the cosmological parameters of the standard $\Lambda$CDM model, as well as a Warm Dark Matter (WDM) model and the $f(R)$ and Symmetron modified gravity models. The galaxy bias is accounted for at the level of a linear bias. The relative information of the 2PCF and the LCF depends on the survey volume, sampling density (shot noise) and the bias uncertainty. For a volume of $1h^{-3}\rm Gpc^3$, sampled with points of mean density $\bar{n} = 2\times10^{-3} h^{3}\ \rm Mpc^{-3}$ and a bias uncertainty of 13\%, the LCF improves the parameter constraints by about 20\% in the $\Lambda$CDM cosmology and potentially even more in alternative models. Finally, since a linear bias only affects the Fourier amplitudes (2PCF), but not the phases (LCF), the combination of the 2PCF and the LCF can be used to break the degeneracy between the linear bias and $\sigma_8$, present in 2-point statistics.

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.

HI-Selected galaxies obey a linear relationship between their maximum detected radius Rmax and rotational velocity. This result covers measurements in the optical, ultraviolet, and HI emission in galaxies spanning a factor of 30 in size and velocity, from small dwarf irregulars to the largest spirals. Hence, galaxies behave as clocks, rotating once a Gyr at the very outskirts of their discs. Observations of a large optically-selected sample are consistent, implying this relationship is generic to disc galaxies in the low redshift Universe. A linear RV relationship is expected from simple models of galaxy formation and evolution. The total mass within Rmax has collapsed by a factor of 37 compared to the present mean density of the Universe. Adopting standard assumptions we find a mean halo spin parameter lambda in the range 0.020 to 0.035. The dispersion in lambda, 0.16 dex, is smaller than expected from simulations. This may be due to the biases in our selection of disc galaxies rather than all halos. The estimated mass densities of stars and atomic gas at Rmax are similar (~0.5 Msun/pc^2) indicating outer discs are highly evolved. The gas consumption and stellar population build time-scales are hundreds of Gyr, hence star formation is not driving the current evolution of outer discs. The estimated ratio between Rmax and disc scale length is consistent with long-standing predictions from monolithic collapse models. Hence, it remains unclear whether disc extent results from continual accretion, a rapid initial collapse, secular evolution or a combination thereof.

Stars form in molecular complexes that are visible as giant clouds ($\sim 10^{5-6} \mathrm{M}_\odot$) in nearby galaxies and as giant clumps ($\sim 10^{8-9}\mathrm{M}_\odot$) in galaxies at redshifts $z\approx1$$-$$3$. Theoretical inferences on the origin and evolution of these complexes often require robust measurements of their characteristic size, which is hard to measure at limited resolution and often ill-defined due to overlap and quasi-fractal substructure. We show that maximum and luminosity-weighted sizes of clumps seen in star formation maps (e.g.\ H$\alpha$) can be recovered statistically using the two-point correlation function (2PCF), if an approximate stellar surface density map is taken as the normalizing random field. After clarifying the link between Gaussian clumps and the 2PCF analytically, we design a method for measuring the diameters of Gaussian clumps with realistic quasi-fractal substructure. This method is tested using mock images of clumpy disk galaxies at different spatial resolutions and perturbed by Gaussian white noise. We find that the 2PCF can recover the input clump scale at $\sim20\%$ accuracy, as long as this scale is larger than the spatial resolution. We apply this method to the local spiral galaxy NGC 5194, as well as to three clumpy turbulent galaxies from the DYNAMO-HST sample. In both cases, our statistical H$\alpha$-clump size measurements agree with previous measurements and with the estimated Jeans lengths. However, the new measurements are free from subjective choices when fitting individual clumps.

We present a general modified maximum likelihood (MML) method for inferring generative distribution functions from uncertain and biased data. The MML estimator is identical to, but easier and many orders of magnitude faster to compute than the solution of the exact Bayesian hierarchical modelling of all measurement errors. As a key application, this method can accurately recover the mass function (MF) of galaxies, while simultaneously dealing with observational uncertainties (Eddington bias), complex selection functions and unknown cosmic large-scale structure. The MML method is free of binning and natively accounts for small number statistics and non-detections. Its fast implementation in the R-package "dftools" is equally applicable to other objects, such as haloes, groups and clusters, as well as observables other than mass. The formalism readily extends to multi-dimensional distribution functions, e.g. a Choloniewski function for the galaxy mass-angular momentum distribution, also handled by dftools. The code provides uncertainties and covariances for the fitted model parameters and approximate Bayesian evidences. We use numerous mock surveys to illustrate and test the MML method, as well as to emphasise the necessity of accounting for observational uncertainties in MFs of modern galaxy surveys.

We present Keck/OSIRIS adaptive optics observations with 150-400 pc spatial sampling of 7 turbulent, clumpy disc galaxies from the DYNAMO sample ($0.07<z<0.2$). DYNAMO galaxies have previously been shown to be well matched in properties to main sequence galaxies at $z\sim1.5$. Integral field spectroscopy observations using adaptive optics are subject to a number of systematics including a variable PSF and spatial sampling, which we account for in our analysis. We present gas velocity dispersion maps corrected for these effects, and confirm that DYNAMO galaxies do have high gas velocity dispersion ($\sigma=40-80$\kms), even at high spatial sampling. We find statistically significant structure in 6 out of 7 galaxies. The most common distance between the peaks in velocity dispersion and emission line peaks is $\sim0.5$~kpc, we note this is very similar to the average size of a clump measured with HST H$\alpha$ maps. This could suggest that the peaks in velocity dispersion in clumpy galaxies likely arise due to some interaction between the clump and the surrounding ISM of the galaxy, though our observations cannot distinguish between outflows, inflows or velocity shear. Observations covering a wider area of the galaxies will be needed to confirm this result.

We introduce a method for modeling disk galaxies designed to take full advantage of data from integral field spectroscopy (IFS). The method fits equilibrium models to simultaneously reproduce the surface brightness, rotation and velocity dispersion profiles of a galaxy. The models are fully self-consistent 6D distribution functions for a galaxy with a Sersic-profile stellar bulge, exponential disk and parametric dark matter halo, generated by an updated version of GalactICS. By creating realistic flux-weighted maps of the kinematic moments (flux, mean velocity and dispersion), we simultaneously fit photometric and spectroscopic data using both maximum-likelihood and Bayesian (MCMC) techniques. We apply the method to a GAMA spiral galaxy (G79635) with kinematics from the SAMI Galaxy Survey and deep $g$- and $r$-band photometry from the VST-KiDS survey, comparing parameter constraints with those from traditional 2D bulge-disk decomposition. Our method returns broadly consistent results for shared parameters, while constraining the mass-to-light ratios of stellar components and reproducing the HI-inferred circular velocity well beyond the limits of the SAMI data. While the method is tailored for fitting integral field kinematic data, it can use other dynamical constraints like central fibre dispersions and \HI circular velocities, and is well-suited for modelling galaxies with a combination of deep imaging and HI and/or optical spectra (resolved or otherwise). Our implementation (MagRite) is computationally efficient and can generate well-resolved models and kinematic maps in under a minute on modern processors.

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 present the results of a deep survey of the nearby Sculptor group and the associated Sculptor filament taken with the Parkes 64-m radio telescope in the 21-cm emission line of neutral hydrogen. We detect 31 HI sources in the Sculptor group/filament, eight of which are new HI detections. We derive a slope of the HI mass function along the Sculptor filament of $\alpha = -1.10^{+0.20}_{-0.11}$, which is significantly flatter than the global mass function and consistent with the flat slopes previously found in other low-density group environments. Some physical process, such as star formation, photoionisation or ram-pressure stripping, must therefore be responsible for removing neutral gas predominantly from low-mass galaxies. All of our HI detections have a confirmed or tentative optical counterpart and are likely associated with luminous rather than 'dark' galaxies. Despite a column density sensitivity of about $4 \times 10^{17}~\mathrm{cm}^{-2}$, we do not find any traces of extragalactic gas or tidal streams, suggesting that the Sculptor filament is, at the current time, a relatively quiescent environment that has not seen any recent major interactions or mergers.

The SKA era is set to revolutionize our understanding of neutral hydrogen (HI) in individual galaxies out to redshifts of z~0.8; and in the z > 6 intergalactic medium through the detection and imaging of cosmic reionization. Direct HI number density constraints will, nonetheless, remain relatively weak out to cosmic noon (z~2) - the epoch of peak star formation and black hole accretion - and beyond. However, as was demonstrated from the 1990s with molecular line observations, this can be overcome by utilising the natural amplification afforded by strong gravitational lensing, which results in an effective increase in integration time by the square of the total magnification (\mu^2) for an unresolved source. Here we outline how a dedicated lensed HI survey will leverage MeerKAT's high sensitivity, frequency coverage, large instantaneous bandwidth, and high dynamic range imaging to enable a lasting legacy of high-redshift HI emission detections well into the SKA era. This survey will not only provide high-impact, rapid-turnaround MeerKAT science commissioning results, but also unveil Milky Way-like systems towards cosmic noon which is not possible with any other SKA precursors/pathfinders. An ambitious lensed HI survey will therefore make a significant impact from MeerKAT commissioning all the way through to the full SKA era, and provide a more complete picture of the HI history of the Universe.

We present the first major release of data from the SAMI Galaxy Survey. This data release focuses on the emission-line physics of galaxies. Data Release One includes data for 772 galaxies, about 20% of the full survey. Galaxies included have the redshift range 0.004 < z < 0.092, a large mass range (7.6 < log(Mstellar/M$_\odot$) < 11.6), and star-formation rates of 10^-4 to 10^1\ M$_\odot$/yr. For each galaxy, we include two spectral cubes and a set of spatially resolved 2D maps: single- and multi-component emission-line fits (with dust extinction corrections for strong lines), local dust extinction and star-formation rate. Calibration of the fibre throughputs, fluxes and differential-atmospheric-refraction has been improved over the Early Data Release. The data have average spatial resolution of 2.16 arcsec (FWHM) over the 15~arcsec diameter field of view and spectral (kinematic) resolution R=4263 (sigma=30km/s) around Halpha. The relative flux calibration is better than 5\% and absolute flux calibration better than $\pm0.22$~mag, with the latter estimate limited by galaxy photometry. The data are presented online through the Australian Astronomical Observatory's Data Central.

We present molecular gas mass estimates for a sample of 13 local galaxies whose kinematic and star forming properties closely resemble those observed in $z\approx 1.5$ main-sequence galaxies. Plateau de Bure observations of the CO[1-0] emission line and Herschel Space Observatory observations of the dust emission both suggest molecular gas mass fractions of ~20%. Moreover, dust emission modeling finds $T_{dust}<$30K, suggesting a cold dust distribution compared to their high infrared luminosity. The gas mass estimates argue that $z\sim$0.1 DYNAMO galaxies not only share similar kinematic properties with high-z disks, but they are also similarly rich in molecular material. Pairing the gas mass fractions with existing kinematics reveals a linear relationship between $f_{gas}$ and $\sigma$/$v_{c}$, consistent with predictions from stability theory of a self-gravitating disk. It thus follows that high gas velocity dispersions are a natural consequence of large gas fractions. We also find that the systems with lowest depletion times ($\sim$0.5 Gyr) have the highest ratios of $\sigma$/$v_{c}$ and more pronounced clumps, even at the same high molecular gas fraction.

The SKA and its pathfinders will enable studies of HI emission at higher redshifts than ever before. In moving beyond the local Universe, this will require the use of cosmologically appropriate formulae that have traditionally been simplified to their low-redshift approximations. In this paper, we summarise some of the most important relations for tracing HI emission in the SKA era, and present an online calculator to assist in the planning and analysis of observations (hifi.icrar.org).

In this letter we study the mean sizes of Halpha clumps in turbulent disk galaxies relative to kinematics, gas fractions, and Toomre Q. We use 100~pc resolution HST images, IFU kinematics, and gas fractions of a sample of rare, nearby turbulent disks with properties closely matched to z~1.5-2 main-sequence galaxies (the DYNAMO sample). We find linear correlations of normalized mean clump sizes with both the gas fraction and the velocity dispersion-to-rotation velocity ratio of the host galaxy. We show that these correlations are consistent with predictions derived from a model of instabilities in a self-gravitating disk (the so-called "violent disk instability model"). We also observe, using a two-fluid model for Q, a correlation between the size of clumps and self-gravity driven unstable regions. These results are most consistent with the hypothesis that massive star forming clumps in turbulent disks are the result of instabilities in self-gravitating gas-rich disks, and therefore provide a direct connection between resolved clump sizes and this in situ mechanism.

We present measurements of baryonic mass Mb and specific angular momentum (sAM) jb in 14 rotating dwarf Irregular (dIrr) galaxies from the LITTLE THINGS sample. These measurements, based on 21cm kinematic data from the Very Large Array and stellar mass maps from the Spitzer Space Telescope, extend previous AM measurements by more than two orders of magnitude in Mb. The dwarf galaxies show systematically higher jb values than expected from the jb~Mb^{2/3} scaling of spiral galaxies, representative of a scale-free galaxy formation scenario. This offset can be explained by decreasing baryon mass fractions fM=Mb Mdyn (where Mdyn is the dynamical mass) with decreasing Mb (for Mb<10^{11}Msun). We find that the sAM of neutral atomic hydrogen HI alone is about 2.5 times higher than that of the stars. The M-j relation of HI is significantly steeper than that of the stars, as a direct consequence of the systematic variation of the HI fraction with Mb.

Dust attenuation in galaxies has been extensively studied nearby, however, there are still many unknowns regarding attenuation in distant galaxies. We contribute to this effort using observations of star-forming galaxies in the redshift range z = 0.05-0.15 from the DYNAMO survey. Highly star-forming DYNAMO galaxies share many similar attributes to clumpy, star-forming galaxies at high redshift. Considering integrated Sloan Digital Sky Survey observations, trends between attenuation and other galaxy properties for DYNAMO galaxies are well matched to star-forming galaxies at high redshift. Integrated gas attenuations of DYNAMO galaxies are 0.2-2.0 mags in the V-band, and the ratio of stellar E(B-V) and gas E(B-V) is 0.78-0.08 (compared to 0.44 at low redshift). Four highly star-forming DYNAMO galaxies were observed at H-alpha using the Hubble Space Telescope and at Pa-alpha using integral field spectroscopy at Keck. The latter achieve similar resolution (~0.8-1 kpc) to our HST imaging using adaptive optics, providing resolved observations of gas attenuations of these galaxies on sub-kpc scales. We find < 1.0 mag of variation in attenuation (at H-alpha) from clump to clump, with no evidence of highly attenuated star formation. Attenuations are in the range 0.3-2.2 mags in the V band, consistent with attenuations of low redshift star-forming galaxies. The small spatial variation on attenuation suggests that a majority of the star-formation activity in these four galaxies occurs in relatively unobscured regions and, thus, star-formation is well characterised by our H-alpha observations.

We present $\sim$100 pc resolution Hubble Space Telescope H$\alpha$ images of 10 galaxies from the DYnamics of Newly-Assembled Massive Objects (DYNAMO) survey of low-$z$ turbulent disk galaxies, and use these to undertake the first detailed systematic study of the effects of resolution and clump clustering on observations of clumps in turbulent disks. In the DYNAMO-{\em HST} sample we measure clump diameters spanning the range $d_{clump} \sim 100-800$~pc, and individual clump star formation rates as high as $\sim5$~M$_{\odot}$~yr$^{-1}$. DYNAMO clumps have very high SFR surface densities, $\Sigma_{SFR}\sim 15$~M$_{\odot}$~yr$^{-1}$~kpc$^{-2}$, $\sim100\times$ higher than in H{\sc ii} regions of nearby spirals. Indeed, SFR surface density provides a simple dividing line between massive star forming clumps and local star forming regions, where massive star forming clumps have $\Sigma_{SFR}> 0.5$~M$_{\odot}$~yr$^{-1}$~kpc$^{-2}$. When degraded to match the observations of galaxies in $z\sim 1-3$ surveys, DYNAMO galaxies are similar in morphology and measured clump properties to clumpy galaxies observed in the high-$z$ Universe. Emission peaks in the simulated high-redshift maps typically correspond to multiple clumps in full resolution images. This clustering of clumps systematically increases the apparent size and SFR of clumps in 1~kpc resolution maps, and decreases the measured SFR surface density of clumps by as much as a factor of 20$\times$. From these results we can infer that clump clustering is likely to strongly effect the measured properties of clumps in high-$z$ galaxies, which commonly have kiloparsec scale resolution.

The velocity function is a fundamental observable statistic of the galaxy population, similarly impor- tant as the luminosity function, but much more difficult to measure. In this work we present the first directly measured circular velocity function that is representative between 60 < v_circ < 320 km/s for galaxies of all morphological types at a given rotation velocity. For the low mass galaxy population (60 < v_circ < 170 km/s), we use the HIPASS velocity function. For the massive galaxy population (170 < v_circ < 320 km/s), we use stellar circular velocities from the Calar Alto Legacy Integral Field Area Survey (CALIFA). In earlier work we obtained the measurements of circular velocity at the 80% light radius for 226 galaxies and demonstrated that the CALIFA sample can produce volume- corrected galaxy distribution functions. The CALIFA velocity function includes homogeneous velocity measurements of both late and early-type rotation-supported galaxies and has the crucial advantage of not missing gas-poor massive ellipticals that HI surveys are blind to. We show that both velocity functions can be combined in a seamless manner, as their ranges of validity overlap. The resulting observed velocity function is compared to velocity functions derived from cosmological simulations of the z = 0 galaxy population. We find that dark matter-only simulations show a strong mismatch with the observed VF. Hydrodynamic simulations fare better, but still do not fully reproduce observations.

We investigate the relationship between stellar and gas specific angular momentum $j$, stellar mass $M_{*}$ and optical morphology for a sample of 488 galaxies extracted from the SAMI Galaxy Survey. We find that $j$, measured within one effective radius, monotonically increases with $M_{*}$ and that, for $M_{*}>$10$^{9.5}$ M$_{\odot}$, the scatter in this relation strongly correlates with optical morphology (i.e., visual classification and S\'ersic index). These findings confirm that massive galaxies of all types lie on a plane relating mass, angular momentum and stellar light distribution, and suggest that the large-scale morphology of a galaxy is regulated by its mass and dynamical state. We show that the significant scatter in the $M_{*}-j$ relation is accounted for by the fact that, at fixed stellar mass, the contribution of ordered motions to the dynamical support of galaxies varies by at least a factor of three. Indeed, the stellar spin parameter (quantified via $\lambda_R$) correlates strongly with S\'ersic and concentration indices. This correlation is particularly strong once slow-rotators are removed from the sample, showing that late-type galaxies and early-type fast rotators form a continuous class of objects in terms of their kinematic properties.

We show that the mass fraction f_atm = 1.35*MHI/M of neutral atomic gas (HI and He) in isolated local disk galaxies of baryonic mass M is well described by a straightforward stability model for flat exponential disks. In the outer disk parts, where gas at the characteristic dispersion of the warm neutral medium is stable in the sense of Toomre (1964), the disk consists of neutral atomic gas; conversely the inner part where this medium would be Toomre-unstable, is dominated by stars and molecules. Within this model, f_atm only depends on a global stability parameter q=j*sigma/(GM), where j is the baryonic specific angular momentum of the disk and sigma the velocity dispersion of the atomic gas. The analytically derived first-order solution f_atm = min{1,2.5q^1.12} provides a good fit to all plausible rotation curves. This model, with no free parameters, agrees remarkably well (+-0.2 dex) with measurements of f_atm in isolated local disk galaxies, even with galaxies that are extremely HI-rich or HI-poor for their mass. The finding that f_atm increases monotonically with q for pure stability reasons offers a powerful intuitive explanation for the mean variation of f_atm with M: in a cold dark matter universe galaxies are expected to follow j~M^(2/3), which implies the average scaling q~M^(-1/3) and hence f_atm~M^(-0.37), in agreement with observations.

There is strong evidence that cosmological N-body simulations dominated by Warm Dark Matter (WDM) contain spurious or unphysical haloes, most readily apparent as regularly spaced low-mass haloes strung along filaments. We show that spurious haloes are a feature of traditional N-body simulations of cosmological structure formation models, including WDM and Cold Dark Matter (CDM) models, in which gravitational collapse proceeds in an initially anisotropic fashion, and arises naturally as a consequence of discreteness-driven relaxation. We demonstrate this using controlled N-body simulations of plane-symmetric collapse and show that spurious haloes are seeded at shell crossing by localised velocity perturbations induced by the discrete nature of the density field, and that their characteristic separation should be approximately the mean inter-particle separation of the N-body simulation, which is fixed by the mass resolution within the volume. Using cosmological N-body simulations in which particles are split into two collisionless components of fixed mass ratio, we find that the spatial distribution of the two components show signatures of discreteness-driven relaxation in their spatial distribution on both large and small scales. Adopting a spline kernel gravitational softening that is of order the comoving mean inter-particle separation helps to suppress the effect of discreteness-driven relaxation, but cannot eliminate it completely. These results provide further motivation for recent developments of new algorithms, which include, for example, revisions of the traditional N-body approach by means of spatially adaptive anistropric gravitational softenings or explicit solutions for the evolution of dark matter in phase space.

Measurement of the evolution of both active galactic nuclei (AGN) and star-formation in galaxies underpins our understanding of galaxy evolution over cosmic time. Radio continuum observations can provide key information on these two processes, in particular via the mechanical feedback produced by radio jets in AGN, and via an unbiased dust-independent measurement of star-formation rates. In this paper we determine radio luminosity functions at 325 MHz for a sample of AGN and star-forming galaxies by matching a 138 deg sq. radio survey conducted with the Giant Metrewave Radio Telescope (GMRT), with optical imaging and redshifts from the Galaxy And Mass Assembly (GAMA) survey. We find that the radio luminosity function at 325 MHz for star-forming galaxies closely follows that measured at 1.4 GHz. By fitting the AGN radio luminosity function out to $z = 0.5$ as a double power law, and parametrizing the evolution as ${\Phi} \propto (1 + z)^{k}$ , we find evolution parameters of $k = 0.92 \pm 0.95$ assuming pure density evolution and $k = 2.13 \pm 1.96$ assuming pure luminosity evolution. We find that the Low Excitation Radio Galaxies are the dominant population in space density at lower luminosities. Comparing our 325 MHz observations with radio continuum imaging at 1.4 GHz, we determine separate radio luminosity functions for steep and flat-spectrum AGN, and show that the beamed population of flat-spectrum sources in our sample can be shifted in number density and luminosity to coincide with the unbeamed population of steep-spectrum sources, as is expected in the orientation based unification of AGN.

We measure the stellar specific angular momentum jstar=Jstar/Mstar in four nearby (z~0.1) disk galaxies that have stellar masses Mstar near the break M* of the galaxy mass function, but look like typical star-forming disks at z~2 in terms of their low stability (Q~1), clumpiness, high ionized gas dispersion (40-50 km/s), high molecular gas fraction (20-30%) and rapid star formation (~20 Msun/yr). Combining high-resolution (Keck-OSIRIS) and large-radius (Gemini-GMOS) spectroscopic maps, only available at low z, we discover that these targets have about three times less stellar angular momentum than typical local spiral galaxies of equal stellar mass and bulge fraction. Theoretical considerations show that this deficiency in angular momentum is the main cause of their low stability, while the high gas fraction plays a complementary role. Interestingly, the low jstar values of our targets are similar to those expected in the M*-population at higher z from the approximate theoretical scaling jstar~(1+z)^(-1/2) at fixed Mstar. This suggests that a change in angular momentum, driven by cosmic expansion, is the main cause for the remarkable difference between clumpy M*-disks at high z (which likely evolve into early-type galaxies) and mass-matched local spirals.

We present a new technique for the statistical evaluation of the Tully-Fisher relation (TFR) using spectral line stacking. This technique has the potential to extend TFR observations to lower masses and higher redshifts than possible through a galaxy-by-galaxy analysis. It further avoids the need for individual galaxy inclination measurements. To quantify the properties of stacked HI emission lines, we consider a simplistic model of galactic disks with analytically expressible line profiles. Using this model, we compare the widths of stacked profiles with those of individual galaxies. We then follow the same procedure using more realistic mock galaxies drawn from the S3-SAX model (a derivative of the Millennium simulation). Remarkably, when stacking the apparent HI lines of galaxies with similar absolute magnitude and random inclinations, the width of the stack is very similar to the width of the deprojected (= corrected for inclination) and dedispersed (= after removal of velocity dispersion) input lines. Therefore, the ratio between the widths of the stack and the deprojected/dedispersed input lines is approximately constant - about 0.93 - with very little dependence on the gas dispersion, galaxy mass, galaxy morphology, and shape of the rotation curve. Finally, we apply our technique to construct a stacked TFR using HIPASS data which already has a well defined TFR based on individual detections. We obtain a B-band TFR with a slope of $-8.5\pm0.4$ and a K-band relation with a slope of $-11.7\pm0.6$ for the HIPASS data set which is consistent with the existing results.

Astronomical data is often uncertain with errors that are heteroscedastic (different for each data point) and covariant between different dimensions. Assuming that a set of D-dimensional data points can be described by a (D - 1)-dimensional plane with intrinsic scatter, we derive the general likelihood function to be maximised to recover the best fitting model. Alongside the mathematical description, we also release the hyper-fit package for the R statistical language (github.com/asgr/hyper.fit) and a user-friendly web interface for online fitting (hyperfit.icrar.org). The hyper-fit package offers access to a large number of fitting routines, includes visualisation tools, and is fully documented in an extensive user manual. Most of the hyper-fit functionality is accessible via the web interface. In this paper we include applications to toy examples and to real astronomical data from the literature: the mass-size, Tully-Fisher, Fundamental Plane, and mass-spin-morphology relations. In most cases the hyper-fit solutions are in good agreement with published values, but uncover more information regarding the fitted model.

Strong gravitational lensing provides some of the deepest views of the Universe, enabling studies of high-redshift galaxies only possible with next-generation facilities without the lensing phenomenon. To date, 21 cm radio emission from neutral hydrogen has only been detected directly out to z~0.2, limited by the sensitivity and instantaneous bandwidth of current radio telescopes. We discuss how current and future radio interferometers such as the Square Kilometre Array (SKA) will detect lensed HI emission in individual galaxies at high redshift. Our calculations rely on a semi-analytic galaxy simulation with realistic HI disks (by size, density profile and rotation), in a cosmological context, combined with general relativistic ray tracing. Wide-field, blind HI surveys with the SKA are predicted to be efficient at discovering lensed HI systems, increasingly so at z > 2. This will be enabled by the combination of the magnification boosts, the steepness of the HI luminosity function at the high-mass end, and the fact that the HI spectral line is relatively isolated in frequency. These surveys will simultaneously provide a new technique for foreground lens selection and yield the highest redshift HI emission detections. More near term (and existing) cm-wave facilities will push the high redshift HI envelope through targeted surveys of known lenses.

We explore trends in galaxy properties with Mpc-scale structures using catalogues of environment and large scale structure from the Galaxy And Mass Assembly (GAMA) survey. Existing GAMA catalogues of large scale structure, group and pair membership allow us to construct galaxy stellar mass functions for different environmental types. To avoid simply extracting the known underlying correlations between galaxy properties and stellar mass, we create a mass matched sample of galaxies with stellar masses between $9.5 \leq \log{M_*/h^{-2} M_{\odot}} \leq 11$ for each environmental population. Using these samples, we show that mass normalised galaxies in different large scale environments have similar energy outputs, $u-r$ colours, luminosities, and morphologies. Extending our analysis to group and pair environments, we show galaxies that are not in groups or pairs exhibit similar characteristics to each other regardless of broader environment. For our mass controlled sample, we fail to see a strong dependence of S\'{e}rsic index or galaxy luminosity on halo mass, but do find that it correlates very strongly with colour. Repeating our analysis for galaxies that have not been mass controlled introduces and amplifies trends in the properties of galaxies in pairs, groups, and large scale structure, indicating that stellar mass is the most important predictor of the galaxy properties we examine, as opposed to environmental classifications.

We derive an analytical expression for a novel large-scale structure observable: the line correlation function. The line correlation function, which is constructed from the three-point correlation function of the phase of the density field, is a robust statistical measure allowing the extraction of information in the non-linear and non-Gaussian regime. We show that, in perturbation theory, the line correlation is sensitive to the coupling kernel F_2, which governs the non-linear gravitational evolution of the density field. We compare our analytical expression with results from numerical simulations and find a 1-sigma agreement for separations r<30 Mpc/h. Fitting formulae for the power spectrum and the non-linear coupling kernel at small scales allow us to extend our prediction into the strongly non-linear regime where we find a 1-sigma agreement with the simulations for r<2 Mpc/h. We discuss the advantages of the line correlation relative to standard statistical measures like the bispectrum. Unlike the latter, the line correlation is independent of the bias, in the regime where the bias is local and linear. Furthermore, the variance of the line correlation is independent of the Gaussian variance on the modulus of the density field. This suggests that the line correlation can probe more precisely the non-linear regime of gravity, with less contamination from the power spectrum variance.

The SAMI Galaxy Survey will observe 3400 galaxies with the Sydney-AAO Multi-object Integral-field spectrograph (SAMI) on the Anglo-Australian Telescope (AAT) in a 3-year survey which began in 2013. We present the throughput of the SAMI system, the science basis and specifications for the target selection, the survey observation plan and the combined properties of the selected galaxies. The survey includes four volume limited galaxy samples based on cuts in a proxy for stellar mass, along with low-stellar mass dwarf galaxies all selected from the Galaxy And Mass Assembly (GAMA) survey. The GAMA regions were selected because of the vast array of ancillary data available, including ultraviolet through to radio bands. These fields are on the celestial equator at 9, 12, and 14.5 hours, and cover a total of 144 square degrees (in GAMA-I). Higher density environments are also included with the addition of eight clusters. The clusters have spectroscopy from 2dFGRS and SDSS and photometry in regions covered by the Sloan Digital Sky Survey (SDSS) and/or VLT Survey Telescope/ATLAS. The aim is to cover a broad range in stellar mass and environment, and therefore the primary survey targets cover redshifts 0.004 < z < 0.095, magnitudes r$_{pet}$ < 19.4, stellar masses $10^{7} - 10^{12}$ M$_{sol}$, and environments from isolated field galaxies through groups to clusters of $10^{15}$ M$_{sol}$.

In recent years cosmology has undergone a revolution, with precise measurements of the microwave background radiation, large galaxy redshift surveys, and the discovery of the recent accelerated expansion of the Universe using observations of distant supernovae. In this light, the SKA enables us to do an ultimate test in cosmology by measuring the expansion rate of the Universe in real time. This can be done by a rather simple experiment of observing the neutral hydrogen (HI) signal of galaxies at two different epochs. The signal will encounter a change in frequency imprinted as the Universe expands over time and thus monitoring the drift in frequencies will provide a real time measure of the cosmic acceleration. Over a period of 12 years one would expected a frequency shift of the order of 0.1 Hz assuming a standard Lambda-CDM cosmology. Based on the sensitivity estimates of the SKA and the number counts of the expected HI galaxies, it is shown that the number counts are sufficiently high to compensate for the observational uncertainties of the measurements and hence allow a statistical detection of the frequency shift. [abstract abridged]

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.

The science achievable with SKA HI surveys will be greatly increased through the combination of HI data with that at other wavelengths. These multiwavelength datasets will enable studies to move beyond an understanding of HI gas in isolation to instead understand HI as an integral part of the highly complex baryonic processes that drive galaxy evolution. As they evolve, galaxies experience a host of environmental and feedback influences, many of which can radically impact their gas content. Important processes include: accretion (hot and cold mode, mergers), depletion (star formation, galactic winds, AGN), phase changes (ionised/atomic/molecular), and environmental effects (ram pressure stripping, tidal effects, strangulation). Governing all of these to various extents is the underlying dark matter distribution. In turn, the result of these processes can significantly alter the baryonic states in which material is finally observed (stellar populations, dust, chemistry) and its morphology (galaxy type, bulge/disk ratio, bars, warps, radial profile). To fully understand the evolution of HI and the role it plays in galactic evolution requires the ability to quantify each of these separate processes, and hence to coordinate SKA HI surveys with extensive multi-band photometric and spectroscopic campaigns. In addition, multiwavelength data is essential for statistical methods of HI analysis such as HI stacking and intensity mapping cross-correlations. In this chapter, we examine some of the principal science motivations for acquiring multiwavelength data to match that from the extragalactic SKA HI surveys, and review the currently planned capacity to achieve this (eg. LSST, Euclid, W-FIRST, SPICA, ALMA, and 4MOST).

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.

The interaction of galaxies with their environment, the Intergalactic Medium (IGM), is an important aspect of galaxy formation. One of the most fundamental, but unanswered questions in the evolution of galaxies is how gas circulates in and around galaxies and how it enters the galaxies to support star formation. We have several lines of evidence that the observed evolution of star formation requires gas accretion from the IGM at all times and on all cosmic scales. This gas remains largely unaccounted for and the outstanding questions are where this gas resides and what the physical mechanisms of accretion are. The gas is expected to be embedded in an extended cosmic web made of sheets and filaments. Such large-scale filaments of gas are expected by cosmological numerical simulations, which have made significant progress in recent years. Such simulations do not only model the large scale structure of the cosmic web, but also investigate the neutral gas component. To truly make significant progress in understanding the distribution of neutral hydrogen in the IGM, column densities of NHI=10^18 cm-2 and below have to be probed over large areas on the sky at sub-arcminute resolution. These are the densities of the faintest structures known today around nearby galaxies, though mostly found with single dish telescopes which do not have the resolution to resolve these structures and investigate any kinematics. Existing interferometers lack the collecting power or short baselines to achieve brightness sensitivities typically below NHI=10^19 cm-2. Reaching lower column densities with current facilities is feasible, however requires prohibitively long observing times. The SKA will for the first time break these barriers, enabling interferometric observations an order of magnitude deeper than current interferometers and with an order of magnitude better linear resolution than single-dish telescopes.

Angular momentum is one of the most fundamental physical quantities governing galactic evolution. Differences in the colours, morphologies, star formation rates and gas fractions amongst galaxies of equal stellar/baryon mass M are potentially widely explained by variations in their specific stellar/baryon angular momentum j. The enormous potential of angular momentum science is only just being realised, thanks to the emergence of the first simulations of galaxies with converged spins, paralleled by a dramatic increase in kinematic observations. Such observations are still challenged by the fact that most of the stellar/baryon angular momentum resides at large radii. In fact, the radius that maximally contributes to the angular momentum of an exponential disk (3Re-4Re) is twice as large as the radius that maximally contributes to the disk mass; thus converged measurements of angular momentum require either extremely deep IFS data or, alternatively, kinematic measurements of neutral atomic hydrogen (HI), which naturally resides at the large disk radii that dominate the angular momentum. The SKA has a unique opportunity to become the world-leading facility for angular momentum studies due to its ability to measure the resolved and/or global HI kinematics in very large and well-characterised galaxy samples. These measurements will allow, for example, (1) a very robust determination of the two-dimensional distribution of galaxies in the (M,j)-plane, (2) the largest, systematic measurement of the relationship between M, j, and tertiary galaxy properties, and (3) the most accurate measurement of the large-scale distribution and environmental dependence of angular momentum vectors, both in terms of norm and orientation. All these measurements will represent exquisite tools to build a next generation of galaxy evolution models.

The stellar and neutral hydrogen (HI) mass functions at z~0 are fundamental benchmarks for current models of galaxy evolution. A natural extension of these benchmarks is the two-dimensional distribution of galaxies in the plane spanned by stellar and HI mass, which provides a more stringent test of simulations, as it requires the HI to be located in galaxies of the correct stellar mass. Combining HI data from the ALFALFA survey, with optical data from SDSS, we find a distinct envelope in the HI-to-stellar mass distribution, corresponding to an upper limit in the HI fraction that varies monotonically over five orders of magnitude in stellar mass. This upper envelope in HI fraction does not favour the existence of a significant population of dark galaxies with large amounts of gas but no corresponding stellar population. The envelope shows a break at a stellar mass of ~10^9 Msun, which is not reproduced by modern models of galaxy populations tracing both stellar and gas masses. The discrepancy between observations and models suggests a mass dependence in gas storage and consumption missing in current galaxy evolution prescriptions. The break coincides with the transition from galaxies with predominantly irregular morphology at low masses to regular disks at high masses, as well as the transition from cold to hot accretion of gas in simulations.

We present a methodology for the regularisation and combination of sparse sampled and irregularly gridded observations from fibre-optic multi-object integral-field spectroscopy. The approach minimises interpolation and retains image resolution on combining sub-pixel dithered data. We discuss the methodology in the context of the Sydney-AAO Multi-object Integral-field spectrograph (SAMI) Galaxy Survey underway at the Anglo-Australian Telescope. The SAMI instrument uses 13 fibre bundles to perform high-multiplex integral-field spectroscopy across a one degree diameter field of view. The SAMI Galaxy Survey is targeting 3000 galaxies drawn from the full range of galaxy environments. We demonstrate the subcritical sampling of the seeing and incomplete fill factor for the integral-field bundles results in only a 10% degradation in the final image resolution recovered. We also implement a new methodology for tracking covariance between elements of the resulting datacubes which retains 90% of the covariance information while incurring only a modest increase in the survey data volume.

We use a highly complete subset of the GAMA-II redshift sample to fully describe the stellar mass dependence of close-pairs and mergers between 10^8 Msun and 10^12 Msun. Using the analytic form of this fit we investigate the total stellar mass accreting onto more massive galaxies across all mass ratios. Depending on how conservatively we select our robust merging systems, the fraction of mass merging onto more massive companions is 2.0%-5.6%. Using the GAMA-II data we see no significant evidence for a change in the close-pair fraction between redshift $z = 0.05-0.2$. However, we find a systematically higher fraction of galaxies in similar mass close-pairs compared to published results over a similar redshift baseline. Using a compendium of data and the function $\gamma_M =A(1+z)m$ to predict the major close-pair fraction, we find fitting parameters of $A = 0.021 \pm 0.001$ and $m = 1.53 \pm 0.08$, which represents a higher low-redshift normalisation and shallower power-law slope than recent literature values. We find that the relative importance of in-situ star-formation versus galaxy merging is inversely correlated, with star-formation dominating the addition of stellar material below Mstar and merger accretion events dominating beyond Mstar. We find mergers have a measurable impact on the whole extent of the GSMF, manifest as a deepening of the dip in the GSMF over the next Gyr and an increase in Mstar by as much as 0.01-0.05 dex.

We distribute an easy-to-use mock catalog of galaxies with detailed neutral atomic hydrogen (HI) and auxiliary molecular and optical properties. The catalog covers a field of 10-by-10 degrees and a redshift range of z=0-1.2. It contains galaxies with 21cm peak flux densities down to 1uJy and is, within this flux limit, complete for HI masses above 10^8 solar masses. Five random realisations of the catalog in ASCII format (~4GB/file) and subtables with HI flux limits of 10u Jy (~500MB/file) and 100uJy$ (~30MB/file) can be downloaded at http://ict.icrar.org/store/staff/do/s3sax.

In this letter we report the discovery of CO fluxes, suggesting very high gas fractions in three disk galaxies seen in the nearby Universe (z ~ 0.1). These galaxies were investigated as part of the DYnamics of Newly Assembled Massive Objects (DYNAMO) survey. High-resolution Hubble Space Telescope imaging of these objects reveals the presence of large star forming clumps in the bodies of the galaxies, while spatially resolved spectroscopy of redshifted Halpha reveals the presence of high dispersion rotating disks. The internal dynamical state of these galaxies resembles that of disk systems seen at much higher redshifts (1 < z < 3). Using CO(1-0) observations made with the Plateau de Bure Interferometer, we find gas fractions of 20-30% and depletion times of tdep ~ 0.5 Gyr (assuming a Milky Way-like CO conversion factor). These properties are unlike those expected for low- redshift galaxies of comparable specific star formation rate, but they are normal for their high-z counterparts. DYNAMO galaxies break the degeneracy between gas fraction and redshift, and we show that the depletion time per specific star formation rate for galaxies is closely tied to gas fraction, independent of redshift. We also show that the gas dynamics of two of our local targets corresponds to those expected from unstable disks, again resembling the dynamics of high-z disks. These results provide evidence that DYNAMO galaxies are local analogues to the clumpy, turbulent disks, which are often found at high redshift.

We study the spatially resolved stellar kinematics of two star-forming galaxies at z = 0.1 from the larger DYnamics of Newly Assembled Massive Objects (DYNAMO) sample. These galaxies, which have been characterized by high levels of star formation and large ionized gas velocity dispersions, are considered possible analogs to high-redshift clumpy disks. They were observed using the GMOS instrument in integral field spectroscopy (IFS) mode at the Gemini Observatory with high spectral resolution (R=5400, equivalent to 24 km/s at the observed wavelengths) and 6 hour exposure times in order to measure the resolved stellar kinematics via absorption lines. We also obtain higher-quality emission line kinematics than previous observations. The spatial resolution (1.2 kpc) is sufficient to show that the ionized gas in these galaxies (as traced by H-beta emission) is morphologically irregular, forming multiple giant clumps while stellar continuum light is smooth and well described by an exponential profile. Clumpy gas morphologies observed in IFS data are confirmed by complementary narrow band H-alpha imaging from the Hubble Space Telescope. Morphological differences between the stars and ionized gas are not reflected dynamically as stellar kinematics are found the be closely coupled to the kinematics of the ionized gas: both components are smoothly rotating with large velocity dispersions (~40 km/s) suggesting that the high gas dispersions are not primarily driven by star-formation feedback. In addition, the stellar population ages of these galaxies are estimated to be quite young (60-500 Myr). The large velocity dispersions measured for these young stars suggest that we are seeing the formation of thick disks and/or stellar bulges in support of recent models which produce these from clumpy galaxies at high redshift.

This work present high-precision measurements of the specific baryon angular momentum jb, contained in stars, atomic gas, and molecular gas, out to ~10 scale radii, in 16 nearby spiral galaxies of the THINGS sample. The accuracy of these measurements improves on existing studies by an order of magnitude, leading to the discovery of a strong correlation between the baryon mass Mb, jb, and the bulge mass fraction B/T, fitted by B/T=-(0.34+-0.03)*log(jb/Mb/[1e-7 kpc km/s/Msun])-(0.04+-0.01) on the full sample range of B/T=0.00-0.32 and Mb/Msun=1e9-1e11. The corresponding relation for the stellar quantities Ms and js is identical within the uncertainties. These M-j-B/T relations likely originate from the proportionality between j/M and the surface density of the disk that dictates its stability against (pseudo-)bulge formation. Using a CDM model, we can approximately explain classical scaling relations, such as the fundamental plane of spiral galaxies, the Tully-Fisher relation, and the mass-size relation, in terms of the M-j-B/T relation. These results advocate the use of mass and angular momentum as the most fundamental quantities of spiral galaxies.

Based on data from the Galaxy and Mass Assembly (GAMA) survey, we report on the discovery of structures that we refer to as `tendrils' of galaxies: coherent, thin chains of galaxies that are rooted in filaments and terminate in neighbouring filaments or voids. On average, tendrils contain 6 galaxies and span 10 $h^{-1}$ Mpc. We use the so-called line correlation function to prove that tendrils represent real structures rather than accidental alignments. We show that voids found in the SDSS-DR7 survey that overlap with GAMA regions contain a large number of galaxies, primarily belonging to tendrils. This implies that void sizes are strongly dependent on the number density and sensitivity limits of a survey. We caution that galaxies in low density regions, that may be defined as `void galaxies,' will have local galaxy number densities that depend on such observational limits and are likely higher than can be directly measured.

Accurate parametrization of galaxies detected in the 21-cm HI emission is of fundamental importance to the measurement of commonly used indicators of galaxy evolution, including the Tully-Fisher relation and the HI mass function. Here, we propose a new analytic function, named the 'busy function', that can be used to accurately describe the characteristic double-horn HI profile of many galaxies. The busy function is a continuous, differentiable function that consists of only two basic functions, the error function, erf(x), and a polynomial, |x|^n, of degree n >= 2. We present the basic properties of the busy function and illustrate its great flexibility in fitting a wide range of HI profiles from the Gaussian profiles of dwarf galaxies to the broad, asymmetric double-horn profiles of spiral galaxies. Applications of the busy function include the accurate and efficient parametrization of observed HI spectra of galaxies and the construction of spectral templates for simulations and matched filtering algorithms. We demonstrate the busy function's power by automatically fitting it to the HI spectra of 1000 galaxies from the HIPASS Bright Galaxy Catalog, using our own C/C++ implementation, and comparing the resulting parameters with the catalogued ones. We also present two methods for determining the uncertainties of observational parameters derived from the fit.

We predict the emission line luminosity functions (LFs) of the first 10 rotational transitions of CO in galaxies at redshift z=0 to z=10. This prediction relies on a recently presented simulation of the molecular cold gas content in ~3e7 evolving galaxies based on the Millennium Simulation. We combine this simulation with a model for the conversion between molecular mass and CO-line intensities, which incorporates the following mechanisms: (i) molecular gas is heated by the CMB, starbursts (SBs), and active galactic nuclei (AGNs); (ii) molecular clouds in dense or inclined galaxies can overlap; (iii) compact gas can attain a smooth distribution in the densest part of disks; (iv) CO-luminosities scale with metallicity changes between galaxies; (v) CO-luminosities are always detected against the CMB. We analyze the relative importance of these effects and predict the cosmic evolution of the CO-LFs. The most notable conclusion is that the detection of regular galaxies (i.e. no AGN, no massive SB) at high z>7 in CO-emission will be dramatically hindered by the weak contrast against the CMB, in contradiction to earlier claims that CMB-heating will ease the detection of high-redshift CO. The full simulation of extragalactic CO-lines and the predicted CO-LFs at any redshift can be accessed online, prior registration required} and they should be useful for the modeling of CO-line surveys with future telescopes, such as ALMA, the LMT, or the SKA.

Power-law relations between tracers of baryonic mass and rotational velocities of disk galaxies, so-called Tully-Fisher relations (TFRs), offer a wealth of applications in galaxy evolution and cosmology. However, measurements of rotational velocities require galaxy inclinations, which are difficult to measure, thus limiting the range of TFR studies. This work introduces a maximum likelihood estimation (MLE) method for recovering the TFR in galaxy samples with limited or no information on inclinations. The robustness and accuracy of this method is demonstrated using virtual and real galaxy samples. Intriguingly, the MLE reliably recovers the TFR of all test samples, even without using any inclination measurements - that is, assuming a sin(i)-distribution for galaxy inclinations. Explicitly, this 'inclination-free MLE' recovers the three TFR parameters (zero-point, slope, scatter) with statistical errors only about 1.5-times larger than the best estimates based on perfectly known galaxy inclinations with zero uncertainty. Thus, given realistic uncertainties, the inclination-free MLE is highly competitive. If inclination measurements have mean errors larger than 10 degrees, it is better not to use any inclinations, than to consider the inclination measurements to be exact. The inclination-free MLE opens interesting perspectives for future HI surveys by the SKA and its pathfinders.

We forecast the abilities of the Atacama Large Millimeter/submillimeter Array (ALMA) and the Square Kilometer Array (SKA) to detect CO and HI emission lines in galaxies at redshift z=3. A particular focus is set on Milky Way (MW) progenitors at z=3 for their detection within 24 h constitutes a key science goal of ALMA. The analysis relies on a semi-analytic model, which permits the construction of a MW progenitor sample by backtracking the cosmic history of all simulated present-day galaxies similar to the real MW. Results: (i) ALMA can best observe a MW at z=3 by looking at CO(3-2) emission. The probability of detecting a random model MW at 3-sigma in 24 h using 75 km/s channels is roughly 50%, and these odds can be increased by co-adding the CO(3-2) and CO(4-3) lines. These lines fall into ALMA band 3, which therefore represents the optimal choice towards MW detections at z=3. (ii) Higher CO transitions contained in the ALMA bands geq6 will be invisible, unless the considered MW progenitor coincidentally hosts a major starburst or an active black hole. (iii) The high-frequency array of SKA, fitted with 28.8 GHz receivers, would be a powerful instrument for observing CO(1-0) at z=3, able to detect nearly all simulated MWs in 24 h. (iv) HI detections in MWs at z=3 using the low-frequency array of SKA will be impossible in any reasonable observing time. (v) SKA will nonetheless be a supreme ha survey instrument through its enormous instantaneous field-of-view (FoV). A one year pointed HI survey with an assumed FoV of 410 sqdeg would reveal at least 10^5 galaxies at z=2.95-3.05. (vi) If the positions and redshifts of those galaxies are known from an optical/infrared spectroscopic survey, stacking allows the detection of HI at z=3 in less than 24 h.

The rich statistics of galaxy rotations as captured by the velocity function (VF) provides invaluable constraints on galactic baryon physics and the nature of dark matter (DM). However, the comparison of observed galaxy rotations against cosmological models is prone to subtle caveats that can easily lead to misinterpretations. Our analysis reveals full statistical consistency between ~5000 galaxy rotations, observed in line-of-sight projection, and predictions based on the standard cosmological model (LCDM) at the mass-resolution of the Millennium simulation (HI line-based circular velocities above ~50 km/s). Explicitly, the HI linewidths in the HI Parkes All Sky Survey (HIPASS) are found consistent with those in S3-SAX, a post-processed semi-analytic model for the Millennium simulation. Previously found anomalies in the VF can be plausibly attributed to (1) the mass-limit of the Millennium simulation, (2) confused sources in HIPASS, (3) inaccurate inclination measurements for optically faint sources, and (4) the non-detectability of gas-poor early-type galaxies. These issues can be bypassed by comparing observations and models using linewidth source counts rather than VFs. We investigate if and how well such source counts can constrain the temperature of DM.

We describe a 325-MHz survey, undertaken with the Giant Metrewave Radio Telescope (GMRT), which covers a large part of the three equatorial fields at 9, 12 and 14.5 h of right ascension from the Herschel-Astrophysical Terahertz Large Area Survey (H-ATLAS) in the area also covered by the Galaxy And Mass Assembly survey (GAMA). The full dataset, after some observed pointings were removed during the data reduction process, comprises 212 GMRT pointings covering ~90 deg^2 of sky. We have imaged and catalogued the data using a pipeline that automates the process of flagging, calibration, self-calibration and source detection for each of the survey pointings. The resulting images have resolutions of between 14 and 24 arcsec and minimum rms noise (away from bright sources) of ~1 mJy/beam, and the catalogue contains 5263 sources brighter than 5 sigma. We investigate the spectral indices of GMRT sources which are also detected at 1.4 GHz and find them to agree broadly with previously published results; there is no evidence for any flattening of the radio spectral index below S_1.4=10 mJy. This work adds to the large amount of available optical and infrared data in the H-ATLAS equatorial fields and will facilitate further study of the low-frequency radio properties of star formation and AGN activity in galaxies out to z~1.

We discover that the mass of dark matter particles mDM is imprinted in phase-correlations of the cosmic density field more significantly than in the 2-point correlation. In particular, phase-correlations trace mDM out to scales about five times larger than the 2-point correlation. This result relies on a new estimator l(r) of pure phase-information in Fourier space, which can be interpreted as a parameter-free and scale-invariant tracer of filament-like structure. Based on simulated density fields we show how mDM can, in principle, be measured using l(r), given a suitably reconstructed density field.

Spherically collapsing cavitation bubbles produce a shock wave followed by a rebound bubble. Here we present a systematic investigation of the energy partition between the rebound and the shock. Highly spherical cavitation bubbles are produced in microgravity, which suppress the buoyant pressure gradient that otherwise deteriorates the sphericity of the bubbles. We measure the radius of the rebound bubble and estimate the shock energy as a function of the initial bubble radius (2-5.6 mm) and the liquid pressure (10-80 kPa). Those measurements uncover a systematic pressure dependence of the energy partition between rebound and shock. We demonstrate that these observations agree with a physical model relying on a first-order approximation of the liquid compressibility and an adiabatic treatment of the non-condensable gas inside the bubble. Using this model we find that the energy partition between rebound and shock is dictated by a single non-dimensional parameter $\xi = \Delta p\gamma^6/[{p_{g0}}^{1/\gamma} (\rho c^2)^{1-1/\gamma}]$, where $\Delta p=p_\infty-p_v$ is the driving pressure, $p_{\infty}$ is the static pressure in the liquid, $p_v$ is the vapor pressure, $p_{g0}$ is the pressure of the non-condensable gas at the maximal bubble radius, $\gamma$ is the adiabatic index of the non-condensable gas, $\rho$ is the liquid density, and $c$ is the speed of sound in the liquid.

Cavitation bubbles collapsing and rebounding in a pressure gradient grad(p) form a "micro-jet" enveloped by a "vapor jet". This letter presents unprecedented observations of the vapor jets formed in a uniform gravity-induced grad(p), modulated aboard parabolic flights. The data uncovers that the normalized jet volume is independent of the liquid density and viscosity and proportional to zeta=grad(p)*R0/p, where R0 is the maximal bubble radius and p is the driving pressure. A derivation inspired by "Kelvin-Blake" considerations confirms this law and reveals its negligible dependence of surface tension. We further conjecture that the jet only pierces the bubble boundary if zeta>0.0004.

The unique confinement of shock waves inside isolated liquid volumes amplifies the density of shock-liquid interactions. We investigate this universal principle through an interdisciplinary study of shock-induced cavitation inside liquid volumes, isolated in 2 and 3 dimensions. By combining high-speed visualizations of ideal water drops realized in microgravity with smoothed particle simulations we evidence strong shock-induced cavitation at the focus of the confined shocks. We extend this analysis to ground-observations of jets and drops using an analytic model, and argue that cavitation caused by trapped shocks offers a distinct mechanism of erosion in high-speed impacts (>100 m/s).

[Abridged] Observations of molecular gas at all redshifts are critical for measuring the cosmic evolution in molecular gas density and understanding the star-formation history of the Universe. The 12CO molecule (J=1-0 transition = 115.27 GHz) is the best proxy for extragalactic H2, which is the gas reservoir from which star formation occurs, and has been detected out to z~6. Typically, redshifted high-J lines are observed at mm-wavelengths, the most commonly targeted systems exhibiting high SFRs (e.g. submm galaxies), and far-IR-bright QSOs. While the most luminous objects are the most readily observed, detections of more typical galaxies with modest SFRs are essential for completing the picture. ALMA will be revolutionary in terms of increasing the detection rate and pushing the sensitivity limit down to include such galaxies, however the limited FoV when observing at such high frequencies makes it difficult to use ALMA for studies of the large-scale structure traced out by molecular gas in galaxies. This article introduces a strategy for a systematic search for molecular gas during the EoR (z~7 and above), capitalizing on the fact that the J=1-0 transition of 12CO enters the upper bands of cm-wave instruments at high-z. The FoV advantage gained by observing at such frequencies, coupled with modern broadband correlators allows significant cosmological volumes to be probed on reasonable timescales. In this article we present an overview of our future observing programme which has been awarded 6,500 hours as one of the Large Survey Projects for MeerKAT, the forthcoming South African SKA pathfinder instrument. Its large FoV and correlator bandwidth, and high-sensitivity provide unprecedented survey speed for such work. An existing astrophysical simulation is coupled with instrumental considerations to demonstrate the feasibility of such observations and predict detection rates.

We present a sky simulation of the atomic HI emission line and the first ten CO rotational emission lines of molecular gas in galaxies beyond the Milky Way. The simulated sky field has a comoving diameter of 500/h Mpc, hence the actual field-of-view depends on the (user-defined) maximal redshift zmax; e.g. for zmax=10, the field of view yields ~4x4 sqdeg. For all galaxies, we estimate the line fluxes, line profiles, and angular sizes of the HI and CO emission lines. The galaxy sample is complete for galaxies with cold hydrogen masses above 10^8 Msun. This sky simulation builds on a semi-analytic model of the cosmic evolution of galaxies in a Lambda-cold dark matter (LCDM) cosmology. The evolving CDM-distribution was adopted from the Millennium Simulation, an N-body CDM-simulation in a cubic box with a side length of 500/h Mpc. This side length limits the coherence scale of our sky simulation: it is long enough to allow the extraction of the baryon acoustic oscillations (BAOs) in the galaxy power spectrum, yet the position and amplitude of the first acoustic peak will be imperfectly defined. This sky simulation is a tangible aid to the design and operation of future telescopes, such the SKA, the LMT, and ALMA. The results presented in this paper have been restricted to a graphical representation of the simulated sky and fundamental dN/dz-analyzes for peak flux density limited and total flux limited surveys of HI and CO. A key prediction is that HI will be harder to detect at redshifts z>2 than predicted by a no-evolution model. The future verification or falsification of this prediction will allow us to qualify the semi-analytic models.

The current status of the HI simulation efforts is presented, in which a self consistent simulation path is described and basic equations to calculate array sensitivities are given. There is a summary of the SKA Design Study (SKADS) sky simulation and a method for implementing it into the array simulator is presented. A short overview of HI sensitivity requirements is discussed and expected results for a simulated HI survey are presented.

The S3-Tools are a set of Python-based routines and interfaces whose purpose is to provide user-friendly access to the SKA Simulated Skies (S3) set of simulations, an effort led by the University of Oxford in the framework of the European Union's SKADS program (http://www.skads-eu.org). The databases built from the S3 simulations are hosted by the Oxford e-Research Center (OeRC), and can be accessed through a web portal at http://s-cubed.physics.ox.ac.uk. This paper focuses on the practical steps involved to make radio images from the S3-SEX and S3-SAX simulations using the S3-Map tool and should be taken as a broad overview. For a more complete description, the interested reader should look up the user's guide. The output images can then be used as input to instrument simulators, e.g. to assess technical designs and observational strategies for the SKA and SKA pathfinders.

Galaxies in the early Universe were more compact and contained more molecular gas than today. In this paper, we revisit the relation between these empirical findings, and we quantitatively predict the cosmic evolution of the surface densities of atomic (HI) and molecular (H2) hydrogen in regular galaxies. Our method uses a pressure-based model for the H2/HI-ratio of the Interstellar Medium, applied to ~3*10^7 virtual galaxies in the Millennium Simulation. We predict that, on average, the HI-surface density of these galaxies saturates at Sigma_HI<10 Msun/pc^2 at all redshifts (z), while H2-surface densities evolve dramatically as Sigma_H2(1+z)^2.4. This scaling is dominated by a (1+z)^2 surface brightness scaling originating from the (1+z)^-1 size scaling of galaxies at high z. Current measurements of Sigma_H2 at high z, derived from CO-observations, tend to have even higher values, which can be quantitatively explained by a selection bias towards merging systems. However, despite the consistency between our high-z predictions and the sparse empirical data, we emphasize that the empirical data potentially suffer from serious selection biases and that the semi-analytic models remain in many regards uncertain. As a case study, we investigate the cosmic evolution of simulated galaxies, which resemble the Milky Way at z=0. We explicitly predict their HI- and H2-distribution at z=1.5, corresponding to the CO-detected galaxy BzK-21000, and at z=3, corresponding to the primary science goal of the Atacama Large Millimeter/submillimeter Array (ALMA).

We use a pressure-based model for splitting cold hydrogen into its atomic (HI) and molecular (H2)components to tackle the co-evolution of HI, H2, and star formation rates (SFR) in ~3e7 simulated galaxies in the Millennium simulation. The main prediction is that galaxies contained similar amounts of HI at redshift z=1-5 than today, but substantially more H2, in quantitative agreement with the strong molecular line emission already detected in a few high redshift galaxies and approximately consistent with inferences from studies of the damped Lyman-alpha absorbers seen in the spectra of quasars. The cosmic H2/HI-ratio is predicted to evolve monotonically as Omega(H2)/Omega(HI) (1+z)^1.6. This decline of the H2/HI-ratio as a function of cosmic time is driven by the growth of galactic disks and the progressive reduction of the mean cold gas pressure. Finally, a comparison between the evolutions of HI, H2, and SFRs reveals two distinct cosmic epochs of star formation: an early epoch (z>3), driven by the evolution of Omega(HI+H2), and a late epoch (z<3), driven by the evolution of Omega(H2)/Omega(HI).

We revisit the mass ratio Rmol between molecular hydrogen (H2) and atomic hydrogen (HI) in different galaxies from a phenomenological and theoretical viewpoint. First, the local H2-mass function (MF) is estimated from the local CO-luminosity function (LF) of the FCRAO Extragalactic CO-Survey, adopting a variable CO-to-H2 conversion fitted to nearby observations. This implies an average H2-density Omega_H2=(6.9+-2.7) 10^5/h and Omega_H2/Omega_HI=0.26+-0.11 in the local Universe. Second, we investigate the correlations between Rmol and global galaxy properties in a sample of 245 local galaxies. Based on these correlations we introduce four phenomenological models for Rmol, which we apply to estimate H2-masses for each HI-galaxy in the HIPASS catalogue. The resulting H2-MFs (one for each model for Rmol) are compared to the reference H2-MF derived from the CO-LF, thus allowing us to determine the Bayesian evidence of each model and to identify a clear best model, in which, for spiral galaxies, Rmol negatively correlates with both galaxy Hubble type and total gas mass. Third, we derive a theoretical model for Rmol for regular galaxies based on an expression for their axially symmetric pressure profile dictating the degree of molecularization. This model is quantitatively similar to the best phenomenological one at redshift z=0, and hence represents a consistent generalization while providing a physical explanation for the dependence of Rmol on global galaxy properties. Applying the best phenomenological model for Rmol to the HIPASS sample, we derive the first integral cold gas-MF (HI+H2+helium) of the local Universe.

We present a simulation of the cosmic evolution of the atomic and molecular phases of the cold hydrogen gas in about 3e7 galaxies, obtained by post-processing the virtual galaxy catalog produced by (De Lucia et al. 2007) on the Millennium Simulation of cosmic structure (Springel et al. 2005). Our method uses a set of physical prescriptions to assign neutral atomic hydrogen (HI) and molecular hydrogen (H2) to galaxies, based on their total cold gas masses and a few additional galaxy properties. These prescriptions are specially designed for large cosmological simulations, where, given current computational limitations, individual galaxies can only be represented by simplistic model-objects with a few global properties. Our recipes allow us to (i) split total cold gas masses between HI, H2, and Helium, (ii) assign realistic sizes to both the HI- and H2-disks, and (iii) evaluate the corresponding velocity profiles and shapes of the characteristic radio emission lines. The results presented in this paper include the local HI- and H2-mass functions, the CO-luminosity function, the cold gas mass--diameter relation, and the Tully-Fisher relation (TFR), which all match recent observational data from the local Universe. We also present high-redshift predictions of cold gas diameters and the TFR, both of which appear to evolve markedly with redshift.

We have developed a semi-empirical simulation of the extragalactic radio continuum sky suitable for aiding the design of next generation radio interferometers such as the Square Kilometre Array (SKA). The emphasis is on modelling the large-scale cosmological distribution of radio sources rather than the internal details of individual galaxies. Here we provide a description of the simulation to accompany the online release of a catalogue of 320 million simulated radio sources. The simulation covers 20x20 deg^2 - a plausible upper limit to the instantaneous field of view attainable with future (e.g. SKA) aperture array technologies - out to redshift z=20, and down to flux density limits of 10 nJy at 151, 610 MHz, 1.4, 4.86 and 18 GHz. Five distinct source types are included: radio-quiet AGN, radio-loud AGN of the FRI and FRII structural classes, and star-forming galaxies, the latter split into populations of quiescent and starbursting galaxies. In our semi-empirical approach, the simulated sources are drawn from observed (or extrapolated) luminosity functions and grafted onto an underlying dark matter density field with biases which reflect their measured large-scale clustering. A numerical Press-Schechter-style filtering of the density field is used to identify and populate clusters of galaxies. Radio source structures are built from point source and elliptical sub-components, and for FRI and FRII sources an orientation-based unification and beaming model is used to partition flux between the core and extended lobes and hotspots. The simulation output can be post-processed to achieve more complete agreement with observational data in the years ahead, with the aim of using these 'idealised skies' in telescope simulators to optimise the design of the SKA itself (abridged).