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

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

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

Using 324 numerically modelled galaxy clusters we investigate the radial and galaxy-halo alignment of dark matter subhaloes and satellite galaxies orbiting within and around them. We find that radial alignment depends on distance to the centre of the galaxy cluster but appears independent of the dynamical state of the central host cluster. Furthermore, we cannot find a relation between radial alignment of the halo or galaxy shape with its own mass. We report that backsplash galaxies, i.e. objects that have already passed through the cluster radius but are now located in the outskirts, show a stronger radial alignment than infalling objects. We further find that there exists a population of well radially aligned objects passing very close to the central cluster's centre which were found to be on highly radial orbit.

Predictions for the ability of 21-cm interferometric experiments to discriminate Epoch of Reionization (EoR) signal models are typically limited by the simplicity of data models, whereby foreground signals and characteristics of the instrument are often simplified or neglected.To move towards more realistic scenarios, we explore the effects of applying more realistic foreground and instrument models to the 21cm signal, and the ability to estimate astrophysical parameters with these additional complexities. We use a highly-optimized version of \textsc{21cmFAST}, integrated into \textsc{21cmMC}, to generate lightcones of the brightness temperature fluctuation for Bayesian parameter estimation. We include a statistical point-source foreground model and an instrument model based on the Murchison Widefield Array (MWA) scaled in observation time to have an effective sensitivity similar to the future Square Kilometre Array (SKA). We also extend the current likelihood prescription to account for the presence of beam convolution and foregrounds, the 2-Dimensional Power Spectrum (PS), and the correlation of PS modes. We use frequency bands between 150 and 180 MHz to constrain the ionizing efficiency ($\zeta$), the minimum virial temperature of halos ($T_{\mathrm{vir}}$), the soft X-ray emissivity per unit Star Formation Rate (SFR) ($L_X/SFR$ ), and the X-ray energy threshold ($E_0$). We find that the inclusion of realistic foregrounds and instrumental components biases the parameter constraints due to unaccounted for cross-power between the EoR signal, foregrounds and thermal noise. This causes estimates of $\zeta$ to be biased by up to $5\sigma$ but the estimates of $T_{vir}$, L$_X$/SFR and E$_0$ remain unaffected and are all within $1\sigma$.

We present SimSpin, a new, public, software framework for generating integral field spectroscopy (IFS) data cubes from N-body/hydrodynamical simulations of galaxies, which can be compared directly with observational datasets. SimSpin provides a consistent method for studying a galaxy's stellar component. It can be used to explore how observationally inferred measurements of kinematics, such as the spin parameter $\lambda_R$, are impacted by the effects of, for example, inclination, seeing conditions, distance, etc. SimSpin is written in R and has been designed to be highly modular, flexible, and extensible. It is already being used by the astrophysics community to generate IFS-like cubes and FITS files for direct comparison of simulations to observations. In this paper, we explain the conceptual framework of SimSpin; how it is implemented in R; and we demonstrate SimSpin's current capabilities, providing as an example a brief investigation of how numerical resolution affects how reliably we can recover the intrinsic stellar kinematics of a simulated galaxy.

We present an analysis of asymmetries in global HI spectra from the extended GALEX Arecibo SDSS Survey (xGASS), a stellar mass-selected and gas fraction-limited survey which is representative of the HI properties of galaxies in the local Universe. We demonstrate that the asymmetry in a HI spectrum is strongly linked to its signal-to-noise meaning that, contrary to what was done in previous works, asymmetry distributions for different samples cannot be compared at face value. We develop a method to account for noise-induced asymmetry and find that the typical galaxy detected by xGASS exhibits higher asymmetry than what can be attributed to noise alone, with 37% of the sample showing asymmetry greater than 10% at an 80% confidence level. We find that asymmetric galaxies contain, on average, 29% less HI mass compared to their symmetric counterparts matched in both stellar mass and signal-to-noise. We also present clear evidence that satellite galaxies, as a population, exhibit more asymmetric HI spectra than centrals and that group central galaxies show a slightly higher rate of HI asymmetries compared to isolated centrals. All these results support a scenario in which environmental processes, in particular those responsible for gas removal, are the dominant driver of asymmetry in xGASS.

Hierarchical models of structure formation predict that dark matter halo assembly histories are characterised by episodic mergers and interactions with other haloes. An accurate description of this process will provide insights into the dynamical evolution of haloes and the galaxies that reside in them. Using large cosmological N-body simulations, we characterise halo orbits to study the interactions between substructure haloes and their hosts, and how different evolutionary histories map to different classes of orbits. We use two new software tools - WhereWolf, which uses halo group catalogues and merger trees to ensure that haloes are tracked accurately in dense environments, and OrbWeaver, which quantifies each halo's orbital parameters. We demonstrate how WhereWolf improves the accuracy of halo merger trees, and we use OrbWeaver to quantify orbits of haloes. We assess how well analytical prescriptions for the merger timescale from the literature compare to measured merger timescales from our simulations and find that existing prescriptions perform well, provided the ratio of substructure-to-host mass is not too small. In the limit of small substructure-to-host mass ratio, we find that the prescriptions can overestimate the merger timescales substantially, such that haloes are predicted to survive well beyond the end of the simulation. This work highlights the need for a revised analytical prescription for the merger timescale that more accurately accounts for processes such as catastrophic tidal disruption.

Galaxy cluster outskirts mark the transition region from the mildly non-linear cosmic web to the highly non-linear, virialised, cluster interior. It is in this transition region that the intra-cluster medium (ICM) begins to influence the properties of accreting galaxies and groups, as ram pressure impacts a galaxy's cold gas content and subsequent star formation rate. Conversely, the thermodynamical properties of the ICM in this transition region should also feel the influence of accreting substructure (i.e. galaxies and groups), whose passage can drive shocks. In this paper, we use a suite of cosmological hydrodynamical zoom simulations of a single galaxy cluster, drawn from the nIFTy comparison project, to study how the dynamics of substructure accreted from the cosmic web influences the thermodynamical properties of the ICM in the cluster's outskirts. We demonstrate how features evident in radial profiles of the ICM (e.g. gas density and temperature) can be linked to strong shocks, transient and short-lived in nature, driven by the passage of substructure. The range of astrophysical codes and galaxy formation models in our comparison are broadly consistent in their predictions (e.g. agreeing when and where shocks occur, but differing in how strong shocks will be); this is as we would expect of a process driven by large-scale gravitational dynamics and strong, inefficiently radiating, shocks. This suggests that mapping such shock structures in the ICM in a cluster's outskirts (via e.g. radio synchrotron emission) could provide a complementary measure of its recent merger and accretion history.

We investigate the impact of time-resolved `gradual' stellar feedback processes in high redshift dwarf spheroidal galaxies. Here `gradual' feedback refers to individual stellar feedback events which deposit energy over a period of time. We conduct high resolution hydrodynamical simulations of dwarf spheroidal galaxies with halo masses of 10$^7$ M$_{\odot}$ - 10$^8$ M$_{\odot}$, based on z = 6 progenitors of the Milky Way's dwarf spheroidal galaxies. We also include a novel feedback prescription for individual massive stars, which includes stellar winds and a HMXB (High Mass X-ray Binary) phase, on top of supernovae. We find the mass of gas unbound across a 1 Gyr starburst is uniformly lowered if gradual feedback mechanisms are included across the range of metallicities, halo concentration parameters and galaxy masses studied here. Furthermore, we find including gradual feedback in the smallest galaxies delays the unbinding of the majority of the gas and facilitates the production of `chimneys' in the dense shell surrounding the feedback generated hot, pressurised `superbubble'. These `chimneys' vent hot gas from the galaxy interior, lowering the temperature of the central 10 kpc of the gaseous halo. Additionally, we find radiative cooling has little effect on the energetics of simulations that include a short, violent starburst compared with those that have a longer, less concentrated starburst. Finally, we investigate the relative impact of HMXB feedback and stellar winds on our results, finding the ubiquity of stellar winds throughout each starburst makes them a defining factor in the final state of the interstellar medium.

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

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.

Several lines of evidence suggest that as dark matter haloes grow their scale radius increases, and that the density in their central region drops. Major mergers seem an obvious mechanism to explain both these phenomena, and the resulting patterns in the concentration--mass--redshift relation. To test this possibility, we have simulated equal-mass mergers between haloes with a variety of cosmological density profiles, placed on various different orbits. The remnants typically have higher densities than the initial conditions, but differ only slightly from self-similar scaling predictions. They are reasonably well fit by Einasto profiles, but have parameters distinct from those of the initial conditions. The net internal energy available to the merger remnant, relative to the internal energy of the initial conditions, $\kappa$, has the greatest influence on the properties of the final mass distribution. As expected, energetic encounters produce more extended remnants while mergers of strongly bound systems produce compact remnants. Surprisingly, however, the scale radius of the density profile shows the opposite trend, {\it increasing} in the remnants of low-energy encounters relative to energetic ones. Also even in the most energetic encounters, the density within the scale radius decreases only slightly (by 10--20\%), while for very low-energy systems it increases significantly after the merger. We conclude that while major mergers can produce remnants that are more diffuse at large radii, they are relatively ineffective at changing the central densities of haloes, and seem unlikely to explain the mean trends in the concentration--mass--redshift relation.

The structural properties of individual dark matter haloes, including shape, spin, concentration, and substructure, are linked to the halo's growth history, but the exact connection between the two is unclear. One open question, in particular, is the effect of major mergers on halo structure. We have performed a large set of simulations of binary equal-mass mergers between isolated haloes with various density profiles, to map out the relationship between the initial conditions and merger parameters and the structure of the final remnant. In this paper we describe our initial set-up and analysis methods, and report on the results for the size, shape, and spin of the merger remnant. The outcomes of mergers are most easily understood in terms of a scaled dimensionless energy parameter $\kappa$ and an angular momentum (or spin) parameter $\lambda$. We find that the axis ratio $c/a$ scales roughly linearly with energy $\kappa$ while the axis ratio $c/b$ scales linearly with spin $\lambda$. Qualitatively, mergers on radial orbits produce prolate remnants, while mergers on tangential orbits produce oblate remnants. The spin of the remnant can be predicted from angular momentum conservation, while the overall size changes as $\sim \kappa^{-5}$, as expected from self-similar scaling at constant mean density. We discuss potential cosmological applications for these simple patterns.

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

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

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.

WAVES is designed to study the growth of structure, mass and energy on scales of ~1 kpc to ~10 Mpc over a 7 Gyr timeline. On the largest length scales (1-10 Mpc) WAVES will measure the structures defined by groups, filaments and voids, and their emergence over recent times. Comparisons with bespoke numerical simulations will be used to confirm, refine or refute the Cold Dark Matter paradigm. At intermediate length scales (10 kpc-1 Mpc) WAVES will probe the size and mass distribution of galaxy groups, as well as the galaxy merger rates, in order to directly measure the assembly of dark matter halos and stellar mass. On the smallest length scales (1-10 kpc) WAVES will provide accurate distance and environmental measurements to complement high-resolution space-based imaging to study the mass and size evolution of galaxy bulges, discs and bars. In total, WAVES will provide a panchromatic legacy dataset of ~1.6 million galaxies, firmly linking the very low ($z < 0.1$) and intermediate ($z \sim 0.8$) redshift Universe.

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

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

Recent numerical studies of the dark matter density profiles of massive galaxy clusters ($M_{\rm halo} > 10^{15}$M$_{\odot}$) show that their median radial mass density profile remains unchanged up to $z > 1$, displaying a highly self-similar evolution. We verify this by using the data set of the THE THREE HUNDRED project, i.e. 324 cluster-sized haloes as found in full physics hydrodynamical simulations. We track the progenitors of the mass-complete sample of clusters at $z=0$, and find that their median shape is already in place by $z=2.5$. However, selecting a dynamically relaxed subsample ($\sim16$ per cent of the clusters), we observe a shift of the scale radius $r_s$ towards larger values at earlier times. Classifying the whole sample by formation time, this evolution is understood as a result of a two-phase halo mass accretion process. Early-forming clusters -- identified as relaxed today -- have already entered their slow accretion phase, hence their mass growth occurs mostly at the outskirts. Late-forming clusters -- which are still unrelaxed today -- are in their fast accretion phase, thus the central region of the clusters is still growing. We conclude that the density profile of galaxy clusters shows a profound self-similarity out to redshifts $z\sim2.5$. This result holds for both gas and total density profiles when including baryonic physics, as reported here for two rather distinct sub-grid models.

The relationship between galaxy properties and environment is a widely discussed topic within astrophysics. Here we use galaxy samples from hydrodynamical re-simulations to examine this relationship. We use the over-density ($\delta_1$) within a $1 h^{-1}{\rm Mpc}$ sphere around a galaxy to evaluate its environment. Then the relations between galaxy properties, such as specific star formation rate(sSFR), fraction of star forming galaxies, $g-r$ colour and $\delta_1$ are examined within three galactic samples formed from galaxies within large clusters, those in the vicinity of large clusters and those in the field. We find tight environmental correlations for these galaxy properties. In brief, galaxies in denser environments tend to be redder and are more likely to be quenched. This is consistent with observations. We find that although the sSFR decreases with $\delta_1$, this is mainly because that galaxies with higher stellar mass reside in environment with higher overdensity. At fixed over-density a galaxy's color is also independent of whether it lives within a cluster or within the field, but the relative fractions of the two samples varies dramatically with over-density and this drives an apparent evolution.

A primary goal of integral field spectroscopic (IFS) surveys is to provide a statistical census of galaxies classified by their internal kinematics. As a result, the observational spin parameter, $\lambda_R$, has become one of the most popular methods of quantifying the relative importance of velocity dispersion and rotation in supporting a galaxy's inner structure. The goal of this paper is to examine the relationship between the observationally deduced $\lambda_R$ and one of the most commonly used theoretical spin parameters in the literature, the Bullock et al. (2001) $\lambda'$. Using a set of $N$-body realisations of galaxies from which we construct mock IFS observations, we measure $\lambda_R$ as an observer would, incorporating the effects of beam smearing and seeing conditions. Assuming parameters typical of current IFS surveys, we confirm that there are strong positive correlations between $\lambda_R$ and measurement radius, and strong negative correlations between $\lambda_R$ and size of the PSF, for late-type galaxies; these biases can be reduced using a recently proposed empirical correction. Once observational biases are corrected for, we find that $\lambda_R$ provides a good approximation to $\sim \sqrt{3}/2 \; \lambda'(\rm R_{\rm eff})$, where $\lambda'$ is evaluated for the galactic stellar component within 1 R$_{\rm eff}$.

We present an innovative approach to the methodology of dynamical modelling, allowing practical reconstruction of the underlying dark matter mass without assuming both the density and anisotropy functions. With this, the mass-anisotropy degeneracy is reduced to simple model inference, incorporating the uncertainties inherent with observational data, statistically circumventing the mass-anisotropy degeneracy in spherical collisionless systems. We also tackle the inadequacy that the Jeans method of moments has on small datasets, with the aid of Generative Adversarial Networks: we leverage the power of artificial intelligence to reconstruct non-parametrically the projected line-of-sight velocity distribution. We show with realistic numerical simulations of dwarf spheroidal galaxies that we can distinguish between competing dark matter distributions and recover the anisotropy and mass profile of the system.

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

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.

Merger trees harvested from cosmological $N$-body simulations encode the assembly histories of dark matter halos over cosmic time, and are a fundamental component of semi-analytical models (SAMs) of galaxy formation. The ability to compare the tools used to construct merger trees, namely halo finders and tree building algorithms, in an unbiased and systematic manner is critical to assess the quality of merger trees. In this paper, we present the dendogram, a novel method to visualise merger trees, which provides a comprehensive characterisation of a halo's assembly history - tracking subhalo orbits, halo merger events, and the general evolution of halo properties. We show the usefulness of the dendogram as a diagnostic tool of merger trees by comparing halo assembly histories from a single $N$-Body simulation analysed with three different halo-finders -\textsc{VELOCIraptor}, \textsc{AHF} and \textsc{Rockstar} - and their associated tree-builders. Based on our analysis of the resulting dendograms, we highlight how they have been used to motivate improvements to \textsc{VELOCIraptor}. The dendogram software is publicly available online, at: https://github.com/rhyspoulton/MergerTree-Dendograms .

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

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

We investigate the effects of gradual heating on the evolution of turbulent molecular clouds of mass $2\times 10^6$ M$_\odot$ and virial parameters ranging between $0.7-1.2$. This gradual heating represents the energy output from processes such as winds from massive stars or feedback from High Mass X-ray binaries (HMXBs), contrasting the impulsive energy injection from supernovae (SNe). For stars with a mass high enough that their lifetime is shorter than the life of the cloud, we include a SN feedback prescription. Including both effects, we investigate the interplay between slow and fast forms of feedback and their effectiveness at triggering/suppressing star formation. We find that SN feedback can carve low density chimneys in the gas, offering a path of least resistance for the energy to escape. Once this occurs the more stable, but less energetic, gradual feedback is able to keep the chimneys open. By funneling the hot destructive gas away from the centre of the cloud, chimneys can have a positive effect on both the efficiency and duration of star formation. Moreover, the critical factor is the number of high mass stars and SNe (and any subsequent HMXBs) active within the free-fall time of each cloud. This can vary from cloud to cloud due to the stochasticity of SN delay times and in HMXB formation. However, the defining factor in our simulations is the efficiency of the cooling, which can alter the Jeans mass required for sink particle formation, along with the number of massive stars in the cloud.

We compare the evolution of voids formed under the standard cosmological model and two alternative cosmological models. The two models are a quintessence model ($\phi$CDM) and a Coupled Dark Matter-Dark Energy (CDE) model, both of which have evolving and interacting dark sectors. From $N$-body adiabatic hydrodynamical simulations of these models, we measure the statistics and quantify the properties of voids over the redshift range $z=1.5-12$: these include their population size, volumes, shapes and average densities. We find that the latter property has potential as a probe of cosmology, particularly dark energy, as significant differences in average void densities exist between the alternative models and the standard model. We postulate that this signature arises from an increased evacuation rate of particles out of voids, or an earlier start to void evacuation, in the alternative models as a direct consequence of the dynamical scalar field, which also leads to greater void merger rates. Additionally, differences between the two alternative models are likely due to the drag force arising from dark sector coupling, acting on dark matter particles in our coupled model.

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.

Analysis of Kepler mission data suggests that the Milky Way includes billions of Earth-like planets in the habitable zone of their host star. Current technology enables the detection of technosignatures emitted from a large fraction of the Galaxy. We describe a search for technosignatures that is sensitive to Arecibo-class transmitters located within ~420 ly of Earth and transmitters that are 1000 times more effective than Arecibo within ~13 000 ly of Earth. Our observations focused on 14 planetary systems in the Kepler field and used the L-band receiver (1.15-1.73 GHz) of the 100 m diameter Green Bank Telescope. Each source was observed for a total integration time of 5 minutes. We obtained power spectra at a frequency resolution of 3 Hz and examined narrowband signals with Doppler drift rates between +/-9 Hz/s. We flagged any detection with a signal-to-noise ratio in excess of 10 as a candidate signal and identified approximately 850 000 candidates. Most (99%) of these candidate signals were automatically classified as human-generated radio-frequency interference (RFI). A large fraction (>99%) of the remaining candidate signals were also flagged as anthropogenic RFI because they have frequencies that overlap those used by global navigation satellite systems, satellite downlinks, or other interferers detected in heavily polluted regions of the spectrum. All 19 remaining candidate signals were scrutinized and none were attributable to an extraterrestrial source.

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

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

The halo mass function, encoding the comoving number density of dark matter halos of a given mass, plays a key role in understanding the formation and evolution of galaxies. As such, it is a key goal of current and future deep optical surveys to constrain the mass function down to mass scales which typically host $L_\star$ galaxies. Motivated by the proven accuracy of Press-Schechter-type mass functions, we introduce a related but purely empirical form consistent with standard formulae to better than 4\% in the medium-mass regime, $10^{10}-10^{13} h^{-1}M_\odot$. In particular, our form consists of 4 parameters, each of which has a simple interpretation, and can be directly related to parameters of the galaxy distribution, such as $L_\star$. Using this form within a hierarchical Bayesian likelihood model, we show how individual mass-measurement errors can be successfully included in a typical analysis whilst accounting for Eddington bias. We apply our form to a question of survey design in the context of a semi-realistic data model, illustrating how it can be used to obtain optimal balance between survey depth and angular coverage for constraints on mass function parameters. Open-source {\tt Python} and {\tt R} codes to apply our new form are provided at \url{http://mrpy.readthedocs.org} and \url{https://cran.r-project.org/web/packages/tggd/index.html} respectively.

Using a series of cosmological simulations that includes one dark-matter-only (DM-only) run, one gas cooling-star formation-supernovae feedback (CSF) run and one that additionally includes feedback from active galactic nuclei (AGNs), we classify the large-scale structures with both a velocity-shear-tensor code (Vweb) and a tidal-tensor code (Pweb). We find that the baryonic processes have almost no impact on large-scale structures -- at least not when classified using aforementioned techniques. More importantly, our results confirm that the gas component alone can be used to infer the filamentary structure of the Universe practically un-biased, which could be applied to cosmology constrains. In addition, the gas filaments are classified with its velocity (Vweb) and density (Pweb) fields, which can theoretically connect to the radio observations, such as HI surveys. This will help us to bias-freely link the radio observations with DM distributions at large scale.

The anisotropic distribution of satellites around the central galaxy of their host halo is well-documented. However the relative impact of baryons and dark matter in shaping this distribution is still debated. Using the simulation Horizon-AGN, the angular distribution of satellite galaxies with respect to their central counterpart and halo is quantified. Below one Rvir, satellites cluster more strongly in the plane of the central, rather than merely tracing the shape of their host halo. This is due to the increased isotropy of inner haloes acquired through their inside-out assembly in vorticity-rich flows along the cosmic web. While the effect of centrals decreases with distance, halos' triaxiality increases, impacting more and more the satellite's distribution. Effects become comparable just outside one virial radius. Above this scale, the filamentary infall also impacts the satellites distribution, dominating above two virial radii. The central's morphology plays a governing role: the alignment w.r.t. the central plane is four times stronger in haloes hosting stellar discs than in spheroids. But the impact of the galactic plane decreases for lower satellite-to-central mass ratios, suggesting this might not hold for dwarf satellites of the Local group. The orientation of the Milky-Way's satellites traces their cosmic filament, their level of coplanarity is consistent with systems of similar mass and cosmic location in Horizon-AGN. However, the strong impact of galactic planes in massive groups and clusters bounds the likelihood of finding a relaxed region where satellites can be used to infer halo shape. The minor-to-major axis ratios for haloes with log(M0/Msun)>13.5 is underestimated by 10%. This error soars quickly to 30-40% for individual halo measurements.

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

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.

The planar arrangement of nearly half the satellite galaxies of M31 has been a source of mystery and speculation since it was discovered. With a growing number of other host galaxies showing these satellite galaxy planes, their stability and longevity have become central to the debate on whether the presence of satellite planes are a natural consequence of prevailing cosmological models, or represent a challenge. Given the dependence of their stability on host halo shape, we look into how a galaxy plane's dark matter environment influences its longevity. An increased number of dark matter subhalos results in increased interactions that hasten the deterioration of an already-formed plane of satellite galaxies in spherical dark halos. The role of total dark matter mass fraction held in subhalos in dispersing a plane of galaxies present non trivial effects on plane longevity as well. But any misalignments of plane inclines to major axes of flattened dark matter halos lead to their lifetimes being reduced to < 3 Gyrs. Distributing > 40% of total dark mass in subhalos in the overall dark matter distribution results in a plane of satellite galaxies that is prone to change through the 5 Gyr integration time period.

This work presents a study of galactic outflows driven by stellar feedback. We extract main sequence disc galaxies with stellar mass $10^9\le$ M$_{\star}/$M$_{\odot} \le 5.7\times10^{10}$ at redshift $z=0$ from the highest resolution cosmological simulation of the Evolution and Assembly of GaLaxies and their Environments (EAGLE) set. Synthetic gas rotation velocity and velocity dispersion ($\sigma$) maps are created and compared to observations of disc galaxies obtained with the Sydney-AAO Multi-object Integral field spectrograph (SAMI), where $\sigma$-values greater than $150$ km s$^{-1}$ are most naturally explained by bipolar outflows powered by starburst activity. We find that the extension of the simulated edge-on (pixelated) velocity dispersion probability distribution depends on stellar mass and star formation rate surface density ($\Sigma_{\rm SFR}$), with low-M$_{\star}/$low-$\Sigma_{\rm SFR}$ galaxies showing a narrow peak at low $\sigma$ ($\sim30$ km s$^{-1}$) and more active, high-M$_{\star}/$high-$\Sigma_{\rm SFR}$ galaxies reaching $\sigma>150$ km s$^{-1}$. Although supernova-driven galactic winds in the EAGLE simulations may not entrain enough gas with T $<10^5$ K compared to observed galaxies, we find that gas temperature is a good proxy for the presence of outflows. There is a direct correlation between the thermal state of the gas and its state of motion as described by the $\sigma$-distribution. The following equivalence relations hold in EAGLE: $i)$ low-$\sigma$ peak $\,\Leftrightarrow\,$ disc of the galaxy $\,\Leftrightarrow\,$ gas with T $<10^5$ K; $ii)$ high-$\sigma$ tail $\,\Leftrightarrow\,$ galactic winds $\,\Leftrightarrow\,$ gas with T $\ge 10^5$ K.

Dwarf galaxies, among the most dark matter dominated structures of our universe, are excellent test-beds for dark matter theories. Unfortunately, mass modelling of these systems suffers from the well documented mass-velocity anisotropy degeneracy. For the case of spherically symmetric systems, we describe a method for non-parametric modelling of the radial and tangential velocity moments. The method is a numerical velocity anisotropy "inversion", with parametric mass models, where the radial velocity dispersion profile, $\sigma_{\mathrm{rr}}^2$ is modeled as a B-spline, and the optimization is a three step process that consists of: (i) an Evolutionary modelling to determine the mass model form and the best B-spline basis to represent $\sigma_{\mathrm{rr}}^2$; (ii) an optimization of the smoothing parameters; (iii) a Markov chain Monte Carlo analysis to determine the physical parameters. The mass-anisotropy degeneracy is reduced into mass model inference, irrespective of kinematics. We test our method using synthetic data. Our algorithm constructs the best kinematic profile and discriminates between competing dark matter models. We apply our method to the Fornax dwarf spheroidal galaxy. Using a King brightness profile and testing various dark matter mass models, our model inference favours a simple mass-follows-light system. We find that the anisotropy profile of Fornax is tangential ($\beta(r) < 0$) and we estimate a total mass of $M_{\text{tot}} = 1.613 ^{+0.050}_{-0.075} \times 10^8 \, \text{M}_{\odot}$, and a mass-to-light ratio of $\Upsilon_V = 8.93 ^{+0.32}_{-0.47} \, (\text{M}_{\odot}/\text{L}_{\odot})$. The algorithm we present is a robust and computationally inexpensive method for non-parametric modelling of spherical clusters independent of the mass-anisotropy degeneracy.

We present a clustering comparison of 12 galaxy formation models (including Semi-Analytic Models (SAMs) and Halo Occupation Distribution (HOD) models) all run on halo catalogues and merger trees extracted from a single {\Lambda}CDM N-body simulation. We compare the results of the measurements of the mean halo occupation numbers, the radial distribution of galaxies in haloes and the 2-Point Correlation Functions (2PCF). We also study the implications of the different treatments of orphan (galaxies not assigned to any dark matter subhalo) and non-orphan galaxies in these measurements. Our main result is that the galaxy formation models generally agree in their clustering predictions but they disagree significantly between HOD and SAMs for the orphan satellites. Although there is a very good agreement between the models on the 2PCF of central galaxies, the scatter between the models when orphan satellites are included can be larger than a factor of 2 for scales smaller than 1 Mpc/h. We also show that galaxy formation models that do not include orphan satellite galaxies have a significantly lower 2PCF on small scales, consistent with previous studies. Finally, we show that the 2PCF of orphan satellites is remarkably different between SAMs and HOD models. Orphan satellites in SAMs present a higher clustering than in HOD models because they tend to occupy more massive haloes. We conclude that orphan satellites have an important role on galaxy clustering and they are the main cause of the differences in the clustering between HOD models and SAMs.

We present a comparison of void properties between the standard model of cosmology, $\Lambda$ Cold Dark Matter ($\Lambda$CDM), and two alternative cosmological models with evolving and interacting dark sectors: a quintessence model ($\phi$CDM) and a Coupled Dark Matter-Dark Energy (CDE) model. Using $N$-body simulations of these models, we derive several measures of void statistics and properties, including distributions of void volume, ellipticity, prolateness, and average density. We find that the volume distribution derived from the CDE simulation deviates from the volume distribution derived from the $\Lambda$CDM simulation in the present-day universe, suggesting that the presence of a coupled dark sector could be observable through this statistic. We also find that the distributions of void ellipticity and prolateness are practically indistinguishable among the three models over the redshift range $z=0.0-1.0$, indicating that simple void shape statistics are insensitive to small changes in dark sector physics. Interestingly, we find that the distributions of average void density measured in each of the three simulations are distinct from each other. In particular, voids on average tend to be emptiest under a quintessence model, and densest under the $\Lambda$CDM model. Our results suggest that it is the scalar field present in both alternative models that causes emptier voids to form, while the coupling of the dark sector mitigates this effect by slowing down the evacuation of matter from voids.

There are many proposed mechanisms driving the morphological transformation of disk galaxies to elliptical galaxies. In this paper, we determine if the observed transformation in low mass groups can be explained by the merger histories of galaxies. We measured the group mass-morphology relation for groups from the Galaxy and Mass Assembly group catalogue with masses from 10$^{11}$ - 10$^{15}$ M$_{\odot}$. Contrary to previous studies, the fraction of elliptical galaxies in our more complete group sample increases significantly with group mass across the full range of group mass. The elliptical fraction increases at a rate of 0.163$\pm$0.012 per dex of group mass for groups more massive than 10$^{12.5}$ M$_{\odot}$. If we allow for uncertainties in the observed group masses, our results are consistent with a continuous increase in elliptical fraction from group masses as low as 10$^{11}$M$_{\odot}$. We tested if this observed relation is consistent with merger activity using a GADGET-2 dark matter simulation of the galaxy groups. We specified that a simulated galaxy would be transformed to an elliptical morphology either if it experienced a major merger or if its cumulative mass gained from minor mergers exceeded 30 per cent of its final mass. We then calculated a group mass-morphology relation for the simulations. The position and slope of the simulated relation were consistent with the observational relation, with a gradient of 0.184$\pm$0.010 per dex of group mass. These results demonstrate a strong correlation between the frequency of merger events and disk-to-elliptical galaxy transformation in galaxy group environments.

Even as our measurements of cosmological parameters improve, the physical nature of the dark sector of the universe largely remains a mystery. Many effects of dark sector models are most prominent at very large scales and will rely on future galaxy surveys to elucidate. In this paper we compare the topological properties of the large scale dark matter distribution in a number of cosmological models using hydrodynamical simulations and the cosmological genus statistic. Genus curves are computed from z = 11 to z = 0 for {\Lambda}CDM, Quintessence and Warm Dark Matter models, over a scale range of 1 to 20 Mpc/h. The curves are analysed in terms of their Hermite spectra to describe the power contained in non-Gaussian deformations to the cosmological density field. We find that the {\Lambda}CDM and {\Lambda}WDM models produce nearly identical genus curves indicating no topological differences in structure formation. The Quintessence model, which differs solely in its expansion history, produces significant differences in the strength and redshift evolution of non-Gaussian modes associated with higher cluster abundances and lower void abundances. These effects are robust to cosmic variance and are characteristically different from those produced by tweaking the parameters of a {\Lambda}CDM model. Given the simplicity and similarity of the models, detecting these discrepancies represents a promising avenue for understanding the effect of non-standard cosmologies on large-scale structure.

We present high resolution simulations of an isolated dwarf spheroidal (dSph) galaxy between redshifts $z\sim10$ and $z\sim 4$, the epoch when several Milky Way dSph satellites experienced extended star formation, in order to understand in detail the physical processes which affect a low-mass halo's ability to retain gas. It is well-established that supernova feedback is very effective at expelling gas from a $3\times 10^7$M$_\odot$ halo, the mass of a typical redshift 10 progenitor of a redshift 0 halo with mass $\sim10^9$M$_\odot$. We investigate the conditions under which such a halo is able to retain sufficient high-density gas to support extended star formation. In particular, we explore the effects of: an increased relative concentration of the gas compared to the dark matter; a higher concentration dark matter halo; significantly lower supernova rates; enhanced metal cooling due to enrichment from earlier supernovae. We show that disk-like gas distributions retain more gas than spherical ones, primarily due to the shorter gas cooling times in the disk. However, a significant reduction in the number of supernovae compared to that expected for a standard initial mass function is still needed to allow the retention of high density gas. We conclude that the progenitors of the observed dSphs would only have retained the gas required to sustain star formation if their mass, concentration and gas morphology were already unusual for those of a dSph-mass halo progenitor by a redshift of 10.

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

The recently discovered vast thin plane of dwarf satellites orbiting the Andromeda Galaxy (M31) adds to the mystery of the small scale distribution of the Local Group's galaxy population. Such well defined planar structures are apparently rare occurrences in cold dark matter cosmological simulations, and we lack a coherent explanation of their formation and existence. In this paper, we explore the long-term survivability of thin planes of dwarfs in galactic halos, focusing, in particular, on systems mimicking the observed Andromeda distribution. The key results show that, in general, planes of dwarf galaxies are fragile, sensitive to the shape of the dark matter halo and other perturbing effects. In fact, long lived planes of satellites only exist in polar orbits in spherical dark matter halos, presenting a challenge to the observed Andromeda plane which is significantly tilted with respect to the optical disk. Our conclusion is that, in standard cosmological models, planes of satellites are generally short lived, and hence we must be located at a relatively special time in the evolution of the Andromeda Plane, lucky enough to see its coherent pattern.

We examine the properties of the galaxies and dark matter haloes residing in the cluster infall region surrounding the simulated $\Lambda$CDM galaxy cluster studied by Elahi et al. (2016) at z=0. The $1.1\times10^{15}h^{-1}\text{M}_{\odot}$ galaxy cluster has been simulated with eight different hydrodynamical codes containing a variety of hydrodynamic solvers and subgrid schemes. All models completed a dark-matter only, non-radiative and full-physics run from the same initial conditions. The simulations contain dark matter and gas with mass resolution $m_{\text{DM}}=9.01\times 10^8h^{-1}\text{M}_{\odot}$ and $m_{\text{gas}}=1.9\times 10^8h^{-1}\text{M}_{\odot}$ respectively. We find that the synthetic cluster is surrounded by clear filamentary structures that contain ~60% of haloes in the infall region with mass ~$10^{12.5} - 10^{14} h^{-1}\text{M}_{\odot}$, including 2-3 group-sized haloes ($> 10^{13}h^{-1}\text{M}_{\odot}$). However, we find that only ~10% of objects in the infall region are subhaloes residing in haloes, which may suggest that there is not much ongoing preprocessing occurring in the infall region at z=0. By examining the baryonic content contained within the haloes, we also show that the code-to-code scatter in stellar fraction across all halo masses is typically ~2 orders of magnitude between the two most extreme cases, and this is predominantly due to the differences in subgrid schemes and calibration procedures that each model uses. Models that do not include AGN feedback typically produce too high stellar fractions compared to observations by at least ~1 order of magnitude.

Using a suite of cosmology simulations of a sample of $> 120$ galaxy clusters with $\log(M_{DM, vir}) \le 14.5$. We compare clusters that form in purely dark matter run and their counterparts in hydro runs and investigate 4 independent parameters, that are normally used to classify dynamical state. We find that the virial ratio $\eta$ in hydro-dynamical runs is $\sim 10$ per cent lower than in the DM run, and there is no clear separation between the relaxed and unrelaxed clusters for any parameter. Further, using the velocity dispersion deviation parameter $\zeta$, which is defined as the ratio between cluster velocity dispersion $\sigma$ and the theoretical prediction $\sigma_t = \sqrt{G M_{total}/R}$, we find that there is a linear correlation between the virial ratio $\eta$ and this $\zeta$ parameter. We propose to use this $\zeta$ parameter, which can be easily derived from observed galaxy clusters, as a substitute of the $\eta$ parameter to quantify the cluster dynamical state.

Deep observations of galaxy outskirts reveal faint extended stellar components (ESCs) of streams, shells, and halos, which are ghostly remnants of the tidal disruption of satellite galaxies. We use cosmological galaxy formation simulations in Cold Dark Matter (CDM) and Warm Dark Matter (WDM) models to explore how the dark matter model influences the spatial, kinematic, and orbital properties of ESCs. These reveal that the spherically averaged stellar mass density at large galacto-centric radius can be depressed by up to a factor of 10 in WDM models relative to the CDM model, reflecting the anticipated suppressed abundance of satellite galaxies in WDM models. However, these differences are much smaller in WDM models that are compatible with observational limits, and are comparable in size to the system-to-system variation we find within the CDM model. This suggests that it will be challenging to place limits on dark matter using only the unresolved ESC.

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.

Warm Dark Matter (WDM) models offer an attractive alternative to the current Cold Dark Matter (CDM) cosmological model. We present a novel method to differentiate between WDM and CDM cosmologies, namely using weak lensing; this provides a unique probe as it is sensitive to all the "matter in the beam", not just dark matter haloes and the galaxies that reside in them, but also the diffuse material between haloes. We compare the weak lensing maps of CDM clusters to those in a WDM model corresponding to a thermally produced $0.5$~keV dark matter particle. Our analysis clearly shows that the weak lensing magnification, convergence and shear distributions can be used to distinguish between CDM and WDM models. WDM models {\em increase} the probability of weak magnifications, with the differences being significant to $\gtrsim5\sigma$, while leaving no significant imprint on the shear distribution. WDM clusters analysed in this work are more homogeneous than CDM ones, and the fractional decrease in the amount of material in haloes is proportional to the average increase in the magnification. This difference arises from matter that would be bound in compact haloes in CDM being smoothly distributed over much larger volumes at lower densities in WDM. Moreover, the signature does not solely lie in the probability distribution function but in the full spatial distribution of the convergence field.

Deep observations of galaxies reveal faint extended stellar components (hereafter ESCs) of streams, shells, and halos. These are a natural prediction of hierarchical galaxy formation, as accreted satellite galaxies are tidally disrupted by their host. We investigate whether or not global properties of the ESC could be used to test of dark matter, reasoning that they should be sensitive to the abundance of low-mass satellites, and therefore the underlying dark matter model. Using cosmological simulations of galaxy formation in the favoured Cold Dark Matter (CDM) and Warm Dark Matter (WDM) models ($m_{\rm WDM}$=0.5,1,2 keV/$c^2$), which suppress the abundance of low-mass satellites, we find that the kinematics and orbital structure of the ESC is consistent across models. However, we find striking differences in its spatial structure, as anticipated -- a factor of $\sim$10 drop in spherically averaged mass density between $\sim$10% and $\sim$75% of the virial radius in the more extreme WDM runs ($m_{\rm WDM}$=0.5, 1 keV/$c^2$) relative to the CDM run. These differences are consistent with the mass assembly histories of the different components, and are present across redshifts. However, even the least discrepant of the WDM models is incompatible with current observational limits on $m_{\rm WDM}$. Importantly, the differences we observe when varying the underlying dark matter are comparable to the galaxy-to-galaxy variation we expect within a fixed dark matter model. This suggests that it will be challenging to place limits on dark matter using only the unresolved spatial structure of the the ESC.

We present the Dark-ages Reionization and Galaxy-formation Observables from Numerical Simulations (DRAGONS) program and Tiamat, the collisionless N-body simulation program upon which DRAGONS is built. The primary trait distinguishing Tiamat from other large simulation programs is its density of outputs at high redshift (100 from z=35 to z=5; roughly one every 10 Myr) enabling the construction of very accurate merger trees at an epoch when galaxy formation is rapid and mergers extremely frequent. We find that the friends-of-friends halo mass function agrees well with the prediction of Watson et al. at high masses, but deviates at low masses, perhaps due to our use of a different halo finder or perhaps indicating a break from "universal" behaviour. We then analyse the dynamical evolution of galaxies during the Epoch of Reionization finding that only a small fraction (~20%) of galactic halos are relaxed. We illustrate this using standard relaxation metrics to establish two dynamical recovery time-scales: i) halos need ~1.5 dynamical times following formation, and ii) ~2 dynamical times following a major (3:1) or minor (10:1) merger to be relaxed. This is remarkably consistent across a wide mass range. Lastly, we use a phase-space halo finder to illustrate that major mergers drive long-lived massive phase-space structures which take many dynamical times to dissipate. This can yield significant differences in the inferred mass build-up of galactic halos and we suggest that care must be taken to ensure a physically meaningful match between the galaxy-formation physics of semi-analytic models and the halo finders supplying their input.

Building on the initial results of the nIFTy simulated galaxy cluster comparison, we compare and contrast the impact of baryonic physics with a single massive galaxy cluster, run with 11 state-of-the-art codes, spanning adaptive mesh, moving mesh, classic and modern SPH approaches. For each code represented we have a dark matter only (DM) and non-radiative (NR) version of the cluster, as well as a full physics (FP) version for a subset of the codes. We compare both radial mass and kinematic profiles, as well as global measures of the cluster (e.g. concentration, spin, shape), in the NR and FP runs with that in the DM runs. Our analysis reveals good consistency (<= 20%) between global properties of the cluster predicted by different codes when integrated quantities are measured within the virial radius R200. However, we see larger differences for quantities within R2500, especially in the FP runs. The radial profiles reveal a diversity, especially in the cluster centre, between the NR runs, which can be understood straightforwardly from the division of codes into classic SPH and non-classic SPH (including the modern SPH, adaptive and moving mesh codes); and between the FP runs, which can also be understood broadly from the division of codes into those that include AGN feedback and those that do not. The variation with respect to the median is much larger in the FP runs with different baryonic physics prescriptions than in the NR runs with different hydrodynamics solvers.

We examine subhaloes and galaxies residing in a simulated LCDM galaxy cluster ($M^{\rm crit}_{200}=1.1\times10^{15}M_\odot/h$) produced by hydrodynamical codes ranging from classic Smooth Particle Hydrodynamics (SPH), newer SPH codes, adaptive and moving mesh codes. These codes use subgrid models to capture galaxy formation physics. We compare how well these codes reproduce the same subhaloes/galaxies in gravity only, non-radiative hydrodynamics and full feedback physics runs by looking at the overall subhalo/galaxy distribution and on an individual objects basis. We find the subhalo population is reproduced to within $\lesssim10\%$ for both dark matter only and non-radiative runs, with individual objects showing code-to-code scatter of $\lesssim0.1$ dex, although the gas in non-radiative simulations shows significant scatter. Including feedback physics significantly increases the diversity. Subhalo mass and $V_{max}$ distributions vary by $\approx20\%$. The galaxy populations also show striking code-to-code variations. Although the Tully-Fisher relation is similar in almost all codes, the number of galaxies with $10^{9}M_\odot/h\lesssim M_*\lesssim 10^{12}M_\odot/h$ can differ by a factor of 4. Individual galaxies show code-to-code scatter of $\sim0.5$ dex in stellar mass. Moreover, strong systematic differences exist, with some codes producing galaxies $70\%$ smaller than others. The diversity partially arises from the inclusion/absence of AGN feedback. Our results combined with our companion papers demonstrate that subgrid physics is not just subject to fine-tuning, but the complexity of building galaxies in all environments remains a challenge. We argue even basic galaxy properties, such as the stellar mass to halo mass, should be treated with errors bars of $\sim0.2-0.4$ dex.

Galaxy clusters are an established and powerful test-bed for theories of both galaxy evolution and cosmology. Accurate interpretation of cluster observations often requires robust identification of the location of the centre. Using a statistical sample of clusters drawn from a suite of cosmological simulations in which we have explored a range of galaxy formation models, we investigate how the location of this centre is affected by the choice of observable - stars, hot gas, or the full mass distribution as can be probed by the gravitational potential. We explore several measures of cluster centre: the minimum of the gravitational potential, which would expect to define the centre if the cluster is in dynamical equilibrium; the peak of the density; the centre of BCG; and the peak and centroid of X-ray luminosity. We find that the centre of BCG correlates more strongly with the minimum of the gravitational potential than the X-ray defined centres, while AGN feedback acts to significantly enhance the offset between the peak X-ray luminosity and minimum gravitational potential. These results highlight the importance of centre identification when interpreting clusters observations, in particular when comparing theoretical predictions and observational data.

We have simulated the formation of a massive galaxy cluster (M$_{200}^{\rm crit}$ = 1.1$\times$10$^{15}h^{-1}M_{\odot}$) in a $\Lambda$CDM universe using 10 different codes (RAMSES, 2 incarnations of AREPO and 7 of GADGET), modeling hydrodynamics with full radiative subgrid physics. These codes include Smoothed-Particle Hydrodynamics (SPH), spanning traditional and advanced SPH schemes, adaptive mesh and moving mesh codes. Our goal is to study the consistency between simulated clusters modeled with different radiative physical implementations - such as cooling, star formation and AGN feedback. We compare images of the cluster at $z=0$, global properties such as mass, and radial profiles of various dynamical and thermodynamical quantities. We find that, with respect to non-radiative simulations, dark matter is more centrally concentrated, the extent not simply depending on the presence/absence of AGN feedback. The scatter in global quantities is substantially higher than for non-radiative runs. Intriguingly, adding radiative physics seems to have washed away the marked code-based differences present in the entropy profile seen for non-radiative simulations in Sembolini et al. (2015): radiative physics + classic SPH can produce entropy cores. Furthermore, the inclusion/absence of AGN feedback is not the dividing line -as in the case of describing the stellar content- for whether a code produces an unrealistic temperature inversion and a falling central entropy profile. However, AGN feedback does strongly affect the overall stellar distribution, limiting the effect of overcooling and reducing sensibly the stellar fraction.

We use the Gigaparsec WiggleZ (GiggleZ) simulations to characterise galaxy bias and its scale dependence for a range of redshifts and halo masses in a standard $\Lambda$LCDM cosmology. Assuming bias converges to a scale independent form at large scales, we develop a phenomenological model which fully expresses the mass and redshift dependence of bias and its scale dependence in real or redshift space. We then use this to illustrate how scale-dependent bias can systematically skew measurements of the growth-rate of cosmic structure obtained from redshift-space distortion measurements. When data is fit only to scales $k_{\rm max}{\le}0.1$ $[h^{-1} \rm{Mpc}]^{-1}$, we find that these effects are significant only for large biases ($b{\gtrsim}3$) at large redshifts ($z{\gtrsim}1$). However, when smaller scales are incorporated ($k_{\rm max}{\le}0.2$ $[h^{-1} \rm{Mpc}]^{-1}$) to increase measurement precision, the combination of reduced statistical uncertainty and increased scale dependent bias can result in highly significant systematics for most large halos across all redshifts. We identify several new interesting aspects of bias, including a significant large-scale bias boost for small halos at low-redshifts due to substructure effects ($\sim$20\% for Milky Way-like systems) and a nearly redshift-independent halo mass (corresponding to a redshift-space bias of ${\sim}1.5$) for which halo bias has little-or-no scale dependence on scales greater than $3$ $[h^{-1} {\rm Mpc}]$. This suggests an optimal strategy of targeting bias ${\sim}{1.5}$ systems for clustering studies which are dominated by systematic uncertainties in how observed halo (or galaxy) distributions map to their underlying mass distribution, such as cosmological measurements of neutrino masses. Code for generating our fitting formula is publicly available at http://gbpoole.github.io/Poole_2014a_code/ (Abridged)

Nuclear star clusters (NCs) are found to exist in the centres of many galaxies and appear to follow scaling relations similar to those of super-massive black holes. Previous analytical work has suggested that such relations are a consequence of feedback regulated growth. We explore this idea using high resolution hydrodynamical simulations, focusing on the validity of the simplifying assumptions made in analytical models. In particular, we investigate feedback emanating from multiple stellar sources rather than from a single source, as is usually assumed, and show that collisions betweens shells of gas swept up by feedback leads to momentum cancellation and the formation of high density clumps and filaments. This high density material is resistant both to expulsion from the galaxy potential and to disruption by feedback; if it falls back onto the NC, we expect the gas to be available for further star formation or for feeding a central black hole. We also note our results may have implications for the evolution of globular clusters and stellar clusters in high redshift dark matter halos.

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

The "Wide Area VISTA Extra-galactic Survey" (WAVES) is a 4MOST Consortium Design Reference Survey which will use the VISTA/4MOST facility to spectroscopically survey ~2million galaxies to $r_{\rm AB} < 22$ mag. WAVES consists of two interlocking galaxy surveys ("WAVES-Deep" and "WAVES-Wide"), providing the next two steps beyond the highly successful 1M galaxy Sloan Digital Sky Survey and the 250k Galaxy And Mass Assembly survey. WAVES will enable an unprecedented study of the distribution and evolution of mass, energy, and structures extending from 1-kpc dwarf galaxies in the local void to the morphologies of 200-Mpc filaments at $z\sim1$. A key aim of both surveys will be to compare comprehensive empirical observations of the spatial properties of galaxies, groups, and filaments, against state-of-the-art numerical simulations to distinguish between various Dark Matter models.

We study cluster mass dark matter haloes, their progenitors and surroundings in an coupled Dark Matter-Dark Energy model and compare it to quintessence and $\Lambda$CDM models with adiabatic zoom simulations. When comparing cosmologies with different expansions histories, growth functions & power spectra, care must be taken to identify unambiguous signatures of alternative cosmologies. Shared cosmological parameters, such as $\sigma_8$, need not be the same for optimal fits to observational data. We choose to set our parameters to $\Lambda$CDM $z=0$ values. We find that in coupled models, where DM decays into DE, haloes appear remarkably similar to $\Lambda$CDM haloes despite DM experiencing an additional frictional force. Density profiles are not systematically different and the subhalo populations have similar mass, spin, and spatial distributions, although (sub)haloes are less concentrated on average in coupled cosmologies. However, given the scatter in related observables ($V_{\rm max},R_{V_{\rm max}}$), this difference is unlikely to distinguish between coupled and uncoupled DM. Observations of satellites of MW and M31 indicate a significant subpopulation reside in a plane. Coupled models do produce planar arrangements of satellites of higher statistical significance than $\Lambda$CDM models, however, in all models these planes are dynamically unstable. In general, the nonlinear dynamics within and near large haloes masks the effects of a coupled dark sector. The sole environmental signature we find is that small haloes residing in the outskirts are more deficient in baryons than their $\Lambda$CDM counterparts. The lack of a pronounced signal for a coupled dark sector strongly suggests that such a phenomena would be effectively hidden from view.

We present a comparison of major methodologies of fast generating mock halo or galaxy catalogues. The comparison is done for two-point and the three-point clustering statistics. The reference catalogues are drawn from the BigMultiDark N-body simulation. Both friend-of-friends (including distinct halos only) and spherical overdensity (including distinct halos and subhalos) catalogs have been used with the typical number density of a large-volume galaxy surveys. We demonstrate that a proper biasing model is essential for reproducing the power spectrum at quasilinear and even smaller scales. With respect to various clustering statistics a methodology based on perturbation theory and a realistic biasing model leads to very good agreement with N-body simulations. However, for the quadrupole of the correlation function or the power spectrum, only the method based on semi-N-body simulation could reach high accuracy (1% level) at small scales, i.e., r<25 Mpc/h or k>0.15 h/Mpc. Full N-body solutions will remain indispensable to produce reference catalogues. Nevertheless, we have demonstrated that the far more efficient approximate solvers can reach a few percent accuracy in terms of clustering statistics at the scales interesting for the large-scale structure analysis after calibration with a few reference N-body calculations. This makes them useful for massive production aimed at covariance studies, to scan large parameter spaces, and to estimate uncertainties in data analysis techniques, such as baryon acoustic oscillation reconstruction, redshift distortion measurements, etc.

We present a comparison of 14 galaxy formation models: 12 different semi-analytical models and 2 halo-occupation distribution models for galaxy formation based upon the same cosmological simulation and merger tree information derived from it. The participating codes have proven to be very successful in their own right but they have all been calibrated independently using various observational data sets, stellar models, and merger trees. In this paper we apply them without recalibration and this leads to a wide variety of predictions for the stellar mass function, specific star formation rates, stellar-to- halo mass ratios, and the abundance of orphan galaxies. The scatter is much larger than seen in previous comparison studies primarily because the codes have been used outside of their native environment within which they are well tested and calibrated. The purpose of the `nIFTy comparison of galaxy formation models' is to bring together as many different galaxy formation modellers as possible and to investigate a common approach to model calibration. This paper provides a unified description for all participating models and presents the initial, uncalibrated comparison as a baseline for our future studies where we will develop a common calibration framework and address the extent to which that reduces the scatter in the model predictions seen here.

We have simulated the formation of a galaxy cluster in a $\Lambda$CDM universe using twelve different codes modeling only gravity and non-radiative hydrodynamics (\art, \arepo, \hydra\ and 9 incarnations of GADGET). This range of codes includes particle based, moving and fixed mesh codes as well as both Eulerian and Lagrangian fluid schemes. The various GADGET implementations span traditional and advanced smoothed-particle hydrodynamics (SPH) schemes. The goal of this comparison is to assess the reliability of cosmological hydrodynamical simulations of clusters in the simplest astrophysically relevant case, that in which the gas is assumed to be non-radiative. We compare images of the cluster at $z=0$, global properties such as mass, and radial profiles of various dynamical and thermodynamical quantities. The underlying gravitational framework can be aligned very accurately for all the codes allowing a detailed investigation of the differences that develop due to the various gas physics implementations employed. As expected, the mesh-based codes ART and AREPO form extended entropy cores in the gas with rising central gas temperatures. Those codes employing traditional SPH schemes show falling entropy profiles all the way into the very centre with correspondingly rising density profiles and central temperature inversions. We show that methods with modern SPH schemes that allow entropy mixing span the range between these two extremes and the latest SPH variants produce gas entropy profiles that are essentially indistinguishable from those obtained with grid based methods.

We describe a NOVel form of Adaptive softening (NovA) for collisionless $N$-body simulations, implemented in the Ramses adaptive mesh refinement code. We introduce a refinement criterion that the particle distribution within each cell be sufficiently isotropic, as measured by its moment of inertia tensor. In this way, collapse is only refined if it occurs along all three axes, ensuring that the softening $\epsilon$ is always of order twice the largest inter-particle spacing in a cell. This more conservative force softening criterion is designed to minimise spurious two-body effects, while maintaining high force resolution in collapsed regions of the flow. We test NovA using an antisymmetric perturbed plane wave collapse (`Valinia' test) before applying it to warm dark matter (WDM) simulations. For the Valinia test, we show that -- unlike the standard $N$-body method -- NovA produces no numerical fragmentation while still being able to correctly capture fine caustics and shells around the collapsing regions. For the WDM simulations, we find that NovA converges significantly more rapidly than standard $N$-body, producing little or no spurious halos on small scales. We show, however, that determining whether or not halos exist below the free streaming mass $M_{\rm fs}$ is complicated by the fact that our halo finder (AHF) likely incorrectly labels some caustics and criss-crossing filaments as halos, while one or two particularly massive filaments appear to fragment in any version of NovA where refinement is allowed. Such massive filaments may be physically unstable to collapse, as is the case for infinite, static, self-gravitating cylinders. We will use NovA in forthcoming papers to study the issue of halo formation below $M_{\rm fs}$; filament stability; and to obtain new constraints on the temperature of dark matter.

We present a simple method of generating approximate synthetic halo catalogues: HALOGEN. This method uses a combination of $2^{nd}$-order Lagrangian Perturbation Theory (2LPT) in order to generate the large-scale matter distribution, analytical mass functions to generate halo masses, and a single-parameter stochastic model for halo bias to position haloes. HALOGEN represents a simplification of similar recently published methods. Our method is constrained to recover the 2-point function at intermediate ($10Mpc/h<r<50Mpc/h$) scales, which we show is successful to within 2 per cent. Larger scales ($\sim100Mpc/h$) are reproduced to within 15 per cent. We compare several other statistics (e.g. power spectrum, point distribution function, redshift space distortions) with results from N-Body simulations to determine the validity of our method for different purposes. One of the benefits of HALOGEN is its flexibility, and we demonstrate this by showing how it can be adapted to varying cosmologies and simulation specifications. A driving motivation for the development of such approximate schemes is the need to compute covariance matrices and study the systematic errors for large galaxy surveys, which requires thousands of simulated realisations. We discuss the applicability of our method in this context, and conclude that it is well suited to mass production of appropriate halo catalogues. The code is publicly available at https://github.com/savila/halogen

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

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

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.

We present the Early Data Release of the Sydney-AAO Multi-object Integral field spectrograph (SAMI) Galaxy Survey. The SAMI Galaxy Survey is an ongoing integral field spectroscopic survey of ~3400 low-redshift (z<0.12) galaxies, covering galaxies in the field and in groups within the Galaxy And Mass Assembly (GAMA) survey regions, and a sample of galaxies in clusters. In the Early Data Release, we publicly release the fully calibrated datacubes for a representative selection of 107 galaxies drawn from the GAMA regions, along with information about these galaxies from the GAMA catalogues. All datacubes for the Early Data Release galaxies can be downloaded individually or as a set from the SAMI Galaxy Survey website. In this paper we also assess the quality of the pipeline used to reduce the SAMI data, giving metrics that quantify its performance at all stages in processing the raw data into calibrated datacubes. The pipeline gives excellent results throughout, with typical sky subtraction residuals in the continuum of 0.9-1.2 per cent, a relative flux calibration uncertainty of 4.1 per cent (systematic) plus 4.3 per cent (statistical), and atmospheric dispersion removed with an accuracy of 0."09, less than a fifth of a spaxel.

We investigate the influence of High Mass X-ray Binaries on their high redshift environments. Using a one-dimensional radiative transfer code, we predict the ionization and temperature profiles surrounding a coeval stellar population, composed of main sequence stars and HMXBs, at various times after its formation. We consider both uniform density surroundings, and a cluster embedded in a 10^8 solar mass NFW halo. HMXBs in a constant density environment produce negligible enhanced ionization because of their high-energy SEDs and short lifetimes. In this case, HMXBs only marginally contribute to the local heating rate. For NFW profiles, radiation from main sequence stars cannot prevent the initially ionized volume from recombining since it is unable to penetrate the high density galactic core. However, HMXB photons stall recombinations behind the front, keeping it partially ionized for longer. The increased electron density in these partially ionized regions promotes further cooling, resulting in lower IGM temperatures. In the context of this starburst model, we have shown that HMXBs do not make a major contribution to reionization or IGM heating. However, X-ray escape fractions are high in both density profile cases. Continuous star formation may result in the build up of X-rays over time, reducing the ionization timescale and potentially leading to low level ionization of the distant IGM.

Comparison of observed satellite galaxies of the Milky Way (hereafter MW) with dark matter subhaloes in cosmological $N$-body simulations of MW-mass haloes suggest that such subhaloes, if they exist, are occupied by satellites in a stochastic fashion. We examine how inefficient massive star formation and associated supernova feedback in high-redshift progenitors of present-day low-mass subhaloes might contribute to this stochasticity. Using a Monte Carlo approach to follow the assembly histories of present-day low-mass haloes with $10^7 \lesssim M \leq 10^{10}$ ${\rm M}_{\odot}$, we identify when cooling and star formation is likely to proceed, and observe that haloes with present-day masses $\lesssim 10^9 {\rm M}_{\odot}$ never grow sufficiently massive to support atomic hydrogen line cooling. Noting that the star formation timescale decreases sharply with stellar mass as $t_{\rm PMS} \propto m_{\ast}^{-2.5}$, we argue that, should the conditions for high mass star formation arise in low-mass haloes, the ensuing supernovae are likely to disrupt ongoing lower-mass star formation and unbind gas within the halo. This potentially star-forming gas is unlikely to be replenished in lower mass haloes because of, e.g. cosmological reionization, and so we expect galaxy formation to be stymied in a manner that depends on host halo assembly history and the efficiency and timing of star formation in proto-galaxies, which we illustrate using a Monte Carlo model. Based on these simple physical arguments, we assert that stochasticity of star formation and feedback is an essential but overlooked ingredient in modelling galaxy formation on the smallest scales.

We study clusters in Warm Dark Matter (WDM) models of a thermally produced dark matter particle $0.5$ keV in mass. We show that, despite clusters in WDM cosmologies having similar density profiles as their Cold Dark Matter (CDM) counterparts, the internal properties, such as the amount of substructure, shows marked differences. This result is surprising as clusters are at mass scales that are {\em a thousand times greater} than that at which structure formation is suppressed. WDM clusters gain significantly more mass via smooth accretion and contain fewer substructures than their CDM brethren. The higher smooth mass accretion results in subhaloes which are physically more extended and less dense. These fine-scale differences can be probed by strong gravitational lensing. We find, unexpectedly, that WDM clusters have {\em higher} lensing efficiencies than those in CDM cosmologies, contrary to the naive expectation that WDM clusters should be less efficient due to the fewer substructures they contain. Despite being less dense, the larger WDM subhaloes are more likely to have larger lensing cross-sections than CDM ones. Additionally, WDM subhaloes typically reside at larger distances, which radially stretches the critical lines associated with strong gravitational lensing, resulting in excess in the number of clusters with large radial cross-sections at the $\sim2\sigma$ level. Though lensing profile for an individual cluster vary significantly with the line-of-sight, the radial arc distribution based on a sample of $\gtrsim100$ clusters may prove to be the crucial test for the presence of WDM.

The nature of the dark sector of the Universe remains one of the outstanding problems in modern cosmology, with the search for new observational probes guiding the development of the next generation of observational facilities. Clues come from tension between the predictions from {\Lambda}CDM and observations of gravitationally lensed galaxies. Previous studies showed that galaxy clusters in the {\Lambda}CDM are not strong enough to reproduce the observed number of lensed arcs. This work aims to constrain the warm dark matter cosmologies by means of the lensing efficiency of galaxy clusters drawn from these alternative models. The lensing characteristics of two samples of simulated clusters in the warm dark matter ({\Lambda}WDM) and cold dark matter ({\Lambda}CDM) cosmologies have been studied. The results show that even though the CDM clusters are more centrally concentrated and contain more substructures, the WDM clusters have slightly higher lensing efficiency than their CDM counterparts. The key difference is that WDM clusters have more extended and more massive subhaloes than CDM analogues. These massive substructures significantly stretch the critical lines and caustics and hence they boost the lensing efficiency of the host halo. Despite the increase in the lensing efficiency due to the contribution of massive substructures in the WDM clusters, this is not enough to resolve the arc statistics problem.

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

In order to understand the processes that quench star formation within rich clusters, we construct a library of subhalo orbits drawn from $\Lambda$CDM cosmological N-body simulations of four rich clusters. The orbits are combined with models of star formation followed by quenching in the cluster environment. These are compared with observed bulge and disc colours and stellar absorption linestrength indices of satellite galaxies. Models in which the bulge stellar populations depend only on the galaxy subhalo mass while the disc quenching depends on the cluster environment are acceptable fits to the data. An exponential disc quenching timescale of 3 - 3.5 Gyr is preferred. Models with short ($\lesssim 1$ Gyr) quenching timescales yield cluster-centric gradients in disc colours and Balmer line indices that are too steep compared to observations. We also examine models in which there is quenching in lower mass groups prior to cluster infall ("pre-processing"), finding that such models are a better fit to the data than models without pre-processing and require similar quenching times. The data slightly prefer models where quenching occurs only for galaxies falling within about 0.5 $r_{200}$. Finally, we have examined models with short quenching timescales of 1 Gyr, but a long delay time of 3 Gyr prior to quenching. All models with short quenching timescales, even such "delayed-then-rapid" quenching models, produce excessively red galaxies near the cluster core and are strongly disfavoured by the data. These results imply that the environments of rich clusters must impact star formation rates of infalling galaxies on relatively long timescales -- several times longer than a typical halo spends within the virial radius of a cluster. This scenario favours gentler quenching mechanisms such as slow "strangulation" over more rapid ram-pressure stripping.

The characteristic prediction of the Cold Dark Matter (CDM) model of cosmological structure formation is that the Universe should contain a wealth of small-scale structure -- low-mass dark matter haloes and subhaloes. However, galaxy formation is inefficient in their shallow potential wells and so we expect these low-mass haloes and subhaloes to be dark. Can we tell the difference between a Universe in which low-mass haloes are present but dark and one in which they never formed, thereby providing a robust test of the CDM model? We address this question using cosmological simulations to examine how properties of low-mass haloes that are potentially accessible to observation, such as their spatial clustering, rate of accretions and mergers onto massive galaxies and the angular momentum content of massive galaxies, differ between a LCDM model and dark matter models in which low-mass halo formation is suppressed. Adopting an effective cut-off mass scale M_cut below which small-scale power is suppressed in the initial conditions, we study dark matter models in which M_cut varies between 5e9 Msol/h and 1e11 Msol/h, equivalent to the host haloes of dwarf and low mass galaxies. Our results show that both the clustering strength of low-mass haloes around galaxy-mass primaries and the rate at which they merge with these primaries is sensitive to the assumed value of M_cut; in contrast, suppressing low-mass halo formation has little influence on the angular momentum content of galaxy-mass haloes -- it is the quiescence or violence of a halo's assembly history that has a more marked effect. However, we expect that measuring the effect on spatial clustering or the merger rate is likely to be observationally difficult for realistic values of M_cut, and so isolating the effect of this small-scale structure would appear to be remarkably difficult to detect, at least in the present day Universe.

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.

Abridged: We simulate a massive galaxy cluster in a LCDM Universe using three different approaches to solving the equations of non-radiative hydrodynamics: `classic' Smoothed Particle Hydrodynamics (SPH); a novel SPH with a higher order dissipation switch (SPHS); and adaptive mesh refinement (AMR). We find that SPHS and AMR are in excellent agreement, with both forming a well-defined entropy core that rapidly converges with increasing mass and force resolution. By contrast, SPH exhibits rather different behaviour. At low redshift, entropy decreases systematically with decreasing cluster-centric radius, converging on ever lower central values with increasing resolution. At higher redshift, SPH is in better agreement with SPHS and AMR but shows much poorer numerical convergence. We trace these discrepancies to artificial surface tension in SPH at phase boundaries. At early times, the passage of massive substructures close to the cluster centre stirs and shocks gas to build an entropy core. At later times, artificial surface tension causes low entropy gas to sink artificially to the centre of the cluster. We use SPHS to study the contribution of numerical versus physical dissipation on the entropy core, and argue that numerical dissipation is required to ensure single-valued fluid quantities in converging flows. However, provided this dissipation occurs only at the resolution limit, and provided that it does not propagate errors to larger scales, its effect is benign. There is no requirement to build `sub-grid' models of unresolved turbulence for galaxy cluster simulations. We conclude that entropy cores in non-radiative simulations of galaxy clusters are physical, resulting from entropy generation in shocked gas during cluster assembly, putting to rest the long-standing puzzle of cluster entropy cores in AMR simulations versus their apparent absence in classic SPH simulations.

The dark matter halo mass function (HMF) is a characteristic property of cosmological structure formation models, quantifying the number density of dark matter haloes per unit mass in the Universe. A key goal of current and planned large galaxy surveys is to measure the HMF and to use it to test theories of dark matter and dark energy. We present a new web application for calculating the HMF -- the frontend HMFcalc and the engine hmf. HMFcalc has been designed to be flexible, efficient and easy to use, providing observational and theoretical astronomers alike with the means to explore standard functional forms of the HMF or to tailor their own. We outline the theoretical background needed to compute the HMF, we show how it has been implemented in hmf, and finally we provide worked examples that illustrate HMFcalc's versatility as an analysis tool.

The parameters governing the standard \Lambda Cold Dark Matter cosmological model have been constrained with unprecedented accuracy by precise measurements of the cosmic microwave background by the Wilkinson Microwave Anisotropy Probe (WMAP) and Planck satellites. Each new data release has refined further our knowledge of quantities -- such as the matter density parameter \Omega_M -- that are imprinted on the dark matter halo mass function (HMF), a powerful probe of dark matter and dark energy models. In this letter we trace how changes in the cosmological parameters over the last decade have influenced uncertainty in our knowledge of the HMF. We show that this uncertainty has reduced significantly since the 3rd WMAP data release, but the rate of this reduction is slowing. This is limited by uncertainty in the normalisation \sigma_8, whose influence is most pronounced at the high mass end of the mass function. Interestingly, we find that the accuracy with which we can constrain the HMF in terms of the cosmological parameters has now reached the point at which it is comparable to the scatter in HMF fitting functions. This suggests that the power of the HMF as a precision probe of dark matter and dark energy hinges on more accurate determination of the theoretical HMF. Finally, we assess prospects of using the HMF to differentiate between Cold and Warm Dark Matter models based on ongoing improvements in measurements of \Omega_M, and we comment briefly on optimal survey strategies for constraining dark matter and dark energy models using the HMF.

There is strong evidence that supermassive black holes reside in all galaxies that contain a stellar spheroid and their mass is tightly correlated with properties such as stellar bulge mass and velocity dispersion. There are also strong theoretical arguments that feedback from supermassive black holes plays an important role in shaping the high mass end of the galaxy mass function, hence to accurately model galaxies we also need to model the black holes. We present a comparison of two black hole growth models implemented within a large-scale, cosmological SPH simulation including star formation and feedback. One model is a modified Bondi-Hoyle prescription that grows black holes based on the smooth density of local gas, while the other is the recently proposed Accretion Disc Particle (ADP) method. This model swallows baryonic particles that pass within an accretion radius of the black hole and adds them to a subgrid accretion disc. Black holes are then grown by material from this disc. We find that both models can reproduce local scaling relations, although the ADP model is offset from the observed relations at high black hole masses. The total black hole mass density agrees between models to within a factor of three, but both struggle to reproduce the black hole mass function. The simulated mass functions are too steep and underestimate the number of intermediate and high mass black holes. In addition, the ADP model swallows excessive amounts of material at the resolution of large-scale, cosmological simulations producing unrealistically large accretion discs. Future work needs to be performed to improve the black hole mass function within simulations. This should be done through the mass growth and feedback as they are strongly coupled and should not be treated as separate entities.

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.

Massive stars at redshifts z > 6 are predicted to have played a pivotal role in cosmological reionization as luminous sources of ultra-violet (UV) photons. However, the remnants of these massive stars could be equally important as X-ray luminous (L_X 1e38 erg/s) high-mass X-ray binaries (HMXBs). Because the absorption cross section of neutral hydrogen decreases sharply with photon energy (proportional to the inverse cube), X-rays can escape more freely than UV photons from the star-forming regions in which they are produced, allowing HMXBs to make a potentially significant contribution to the ionizing X-ray background during reionization. In this paper, we explore the ionizing power of HMXBs at redshifts z > 6 using a Monte Carlo model for a coeval stellar population of main sequence stars and HMXBs. Using the archetypal Galactic HMXB Cygnus X-1 as our template, we propose a composite HMXB spectral energy distribution consisting of black-body and power-law components, whose contributions depend on the accretion state of the system. We determine the time-dependent ionizing power of a combined population of UV-luminous stars and X-ray luminous HMXBs, and deduce fitting formulae for the boost in the population's ionizing power arising from HMXBs; these fits allow for simple implementation of HMXB feedback in numerical simulations. Based on this analysis, we estimate the contribution of high redshift HMXBs to the present-day soft X-ray background, and we show that it is a factor of ~100-1000 smaller than the observed limit. Finally, we discuss the implications of our results for the role of HMXBs in reionization and in high redshift galaxy formation.

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

We investigate the means by which cold gas can accrete onto Milky Way mass galaxies from a hot corona of gas, using a new smoothed particle hydrodynamics code, 'SPHS'. We find that the 'cold clumps' seen in many classic SPH simulations in the literature are not present in our SPHS simulations. Instead, cold gas condenses from the halo along filaments that form at the intersection of supernovae-driven bubbles from previous phases of star formation. This positive feedback feeds cold gas to the galactic disc directly, fuelling further star formation. The resulting galaxies in the SPH and SPHS simulations differ greatly in their morphology, gas phase diagrams, and stellar content. We show that the classic SPH cold clumps owe to a numerical thermal instability caused by an inability for cold gas to mix in the hot halo. The improved treatment of mixing in SPHS suppresses this instability leading to a dramatically different physical outcome. In our highest resolution SPHS simulation, we find that the cold filaments break up into bound clumps that form stars. The filaments are overdense by a factor of 10-100 compared to the surrounding gas, suggesting that the fragmentation results from a physical non-linear instability driven by the overdensity. This 'fragmenting filament' mode of disc growth has important implications for galaxy formation, in particular the role of star formation in bringing cold gas into disc galaxies.

Analysing all Galaxy and Mass Assembly (GAMA) galaxies within a factor two (+/- 0.3 dex) of the stellar mass of the Milky Way (MW), there is a 11.9% chance that one of these galaxies will have a close companion (within a projected separation of 70 kpc and radial separation of 400 km/s) that is at least as massive as the Large Magellanic Cloud (LMC). Two close companions at least as massive as the Small Magellanic Cloud (SMC) are rare at the 3.4% level. Two full analogues to the MW- LMC-SMC system were found in GAMA (all galaxies late-type and star forming), suggesting such a combination of close together, late-type, star-forming galaxies is rare: only 0.4% of MW mass galaxies (in the range where we could observe both the LMC and SMC) have such a system. In summary, the MW-LMC-SMC system is a 2.7? event (when recast into Gaussian statistics). Using cross-correlation comparisons we find that there is a preference for SMC- LMC binary pair analogues to be located within 2 Mpc of a range of different lumi- nosity groups. There is a particular preference is for such binaries to be located near LG luminosity systems. When these groups are subdivided into small magnitude gap and large magnitude gap subsets, the binaries prefer to be spatially associated with the small magnitude gap systems. These systems will be dynamically less evolved, but still offer the same amount of gravitational dark matter. This suggests that binaries such as the SMC-LMC might be transient systems, usually destroyed during vigorous merger events. Details of a particularly striking analogue to the MW-SMC-LMC and M31 complex are included.

We argue that current observations of $M - \sigma$ relations for galaxies can be used to constrain theories of super-massive black holes (SMBH) feeding. In particular, assuming that SMBH mass is limited only by the feedback on the gas that feeds it, we show that SMBHs fed via a planar galaxy scale gas flow, such as a disc or a bar, should be much more massive than their counterparts fed by quasi-spherical inflows. This follows from the relative inefficiency of AGN feedback on a flattened inflow. We find that even under the most optimistic conditions for SMBH feedback on flattened inflows, the mass at which the SMBH expels the gas disc and terminates its own growth is a factor of several higher than the one established for quasi-spherical inflows. Any beaming of feedback away from the disc and any disc self-shadowing strengthens this result further. Contrary to this theoretical expectation, recent observations have shown that SMBH in pseudobulge galaxies (which are associated with barred galaxies) are typically under- rather than over-massive when compared with their classical bulge counterparts at a fixed value of $\sigma$. We conclude from this that SMBHs are not fed by large (100 pc to many kpc) scale gas discs or bars, most likely because such planar flows are turned into stars too efficiently to allow any SMBH growth. Based on this and other related observational evidence, we argue that most SMBHs grow by chaotic accretion of gas clouds with a small and nearly randomly distributed direction of angular momentum.

Accretion onto supermassive black holes (SMBHs) in galaxy formation simulations is frequently modelled by the Bondi-Hoyle formalism. Here we examine the validity of this approach analytically and numerically. We argue that the character of the flow where one evaluates the gas properties is unlikely to satisfy the simple Bondi-Hoyle model. Only in the specific case of hot virialised gas with zero angular momentum and negligible radiative cooling is the Bondi-Hoyle solution relevant. In the opposite extreme, where the gas is in a state of free-fall at the evaluation radius due to efficient cooling and the dominant gravity of the surrounding halo, the Bondi-Hoyle formalism can be erroneous by orders of magnitude in either direction. This may impose artificial trends with halo mass in cosmological simulations by being wrong by different factors for different halo masses. We propose an expression for the sub-grid accretion rate which interpolates between the free-fall regime and the Bondi-Hoyle regime, therefore taking account of the contribution of the halo to the gas dynamics.

Using a statistical sample of dark matter haloes drawn from a suite of cosmological N-body simulations of the Cold Dark Matter (CDM) model, we quantify the impact of a simulated halo's mass accretion and merging history on two commonly used measures of its dynamical state, the virial ratio eta and the centre of mass offset Delta r. Quantifying this relationship is important because the degree to which a halo is dynamically equilibrated will influence the reliability with which we can measure characteristic equilibrium properties of the structure and kinematics of a population of haloes. We begin by verifying that a halo's formation redshift zform correlates with its virial mass Mvir and we show that the fraction of its recently accreted mass and the likelihood of it having experienced a recent major merger increases with increasing Mvir and decreasing zform. We then show that both eta and Delta r increase with increasing Mvir and decreasing zform, which implies that massive recently formed haloes are more likely to be dynamically unrelaxed than their less massive and older counterparts. Our analysis shows that both eta and Delta r are good indicators of a halo's dynamical state, showing strong positive correlations with recent mass accretion and merging activity, but we argue that Delta r provides a more robust and better defined measure of dynamical state for use in cosmological N-body simulations at z~0. We find that Delta r < 0.04 is sufficient to pick out dynamically relaxed haloes at z=0. Finally, we assess our results in the context of previous studies, and consider their observational implications.

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

Due to their large dynamical mass-to-light ratios, dwarf spheroidal galaxies (dSphs) are promising targets for the indirect detection of dark matter (DM) in gamma-rays. We examine their detectability by present and future gamma-ray observatories. The key innovative features of our analysis are: (i) We take into account the angular size of the dSphs; while nearby objects have higher gamma ray flux, their larger angular extent can make them less attractive targets for background-dominated instruments. (ii) We derive DM profiles and the astrophysical J-factor (which parameterises the expected gamma-ray flux, independently of the choice of DM particle model) for the classical dSphs directly from photometric and kinematic data. We assume very little about the DM profile, modelling this as a smooth split-power law distribution, with and without sub-clumps. (iii) We use a Markov Chain Monte Carlo (MCMC) technique to marginalise over unknown parameters and determine the sensitivity of our derived J-factors to both model and measurement uncertainties. (iv) We use simulated DM profiles to demonstrate that our J-factor determinations recover the correct solution within our quoted uncertainties. Our key findings are: (i) Sub-clumps in the dSphs do not usefully boost the signal; (ii) The sensitivity of atmospheric Cherenkov telescopes to dSphs within 20 kpc with cored halos can be up to ~50 times worse than when estimated assuming them to be point-like. Even for the satellite-borne Fermi-LAT the sensitivity is significantly degraded on the relevant angular scales for long exposures, hence it is vital to consider the angular extent of the dSphs when selecting targets; (iii) No DM profile has been ruled out by current data, but using a prior on the inner dark matter cusp slope 0<=gamma<=1 provides J-factor estimates accurate to a factor of a few if an appropriate angular scale [abridged]

It has long been recognised that the main obstacle to accretion of gas onto supermassive black holes (SMBHs) is large specific angular momentum. It is feared that the gas settles in a large scale disc, and that accretion would then proceed too inefficiently to explain the masses of the observed SMBHs. Here we point out that, while the mean angular momentum in the bulge is very likely to be large, the deviations from the mean can also be significant. Indeed, cosmological simulations show that velocity and angular momentum fields of gas flows onto galaxies are very complex. Furthermore, inside bulges the gas velocity distribution can be further randomised by the velocity kicks due to feedback from star formation. We perform hydrodynamical simulations of gaseous rotating shells infalling onto an SMBH, attempting to quantify the importance of velocity dispersion in the gas at relatively large distances from the black hole. We implement this dispersion by means of a supersonic turbulent velocity spectrum. We find that, while in the purely rotating case the circularisation process leads to efficient mixing of gas with different angular momentum, resulting in a low accretion rate, the inclusion of turbulence increases this accretion rate by up to several orders of magnitude. We show that this can be understood based on the notion of "ballistic" accretion, whereby dense filaments, created by convergent turbulent flows, travel through the ambient gas largely unaffected by hydrodynamical drag. We derive a simple analytical formula that captures the numerical results to within a factor of a few. Rescaling our results to astrophysical bulges, we argue that this "ballistic" mode of accretion could provide the SMBHs with a sufficient supply of fuel without the need to channel the gas via large-scale discs or bars, and therefore that star formation in bulges can be a strong catalyst for SMBH accretion.

We discuss massive outflows in galaxy bulges, particularly ones driven by accretion episodes where the central supermassive black hole reaches the Eddington limit. We show that the quasar radiation field Compton--cools the wind shock until this reaches distances $\sim 1$ kpc from the black hole, but becomes too dilute to do this at larger radii. Radiative processes cannot cool the shocked gas within the flow time at any radius. Outflows are therefore momentum-driven at small radii (as required to explain the $M - \sigma$ relation). At large radii they are energy-driven, contrary to recent claims. We solve analytically the motion of an energy--driven shell after the central source has turned off. This shows that the thermal energy in the shocked wind can drive further expansion for a time $\sim 10$ times longer than the active time of the central source. Outflows observed at large radii with no active central source probably result from an earlier short (few Myr) active phase of this source.

We show that star formation in galaxy bulges is self-regulating through momentum feedback, limiting the stellar bulge mass to M_b ~ sigma^4. Together with a black hole mass M_BH ~ sigma^4 set by AGN momentum feedback, this produces a linear M_BH - M_b relation. At low redshift this gives M_BH/M_b ~ 0.001, close to the observed ratio. We show that AGN feedback can remove any remaining gas from the bulge and terminate star formation once the central black hole reaches the M_BH - sigma value, contrary to earlier claims. We find a mild upward deviation from the sigma^4 law at higher redshift and at higher sigma.

Black holes grow by accreting matter from their surroundings. However, angular momentum provides an efficient natural barrier to accretion and so only the lowest angular momentum material will be available to feed the black holes. The standard sub-grid model for black hole accretion in galaxy formation simulations - based on the Bondi-Hoyle method - does not account for the angular momentum of accreting material, and so it is unclear how representative the black hole accretion rate estimated in this way is likely to be. In this paper we introduce a new sub-grid model for black hole accretion that naturally accounts for the angular momentum of accreting material. Both the black hole and its accretion disc are modelled as a composite accretion disc particle. Gas particles are captured by the accretion disc particle if and only if their orbits bring them within its accretion radius R_acc, at which point their mass is added to the accretion disc and feeds the black hole on a viscous timescale t_visc. The resulting black hole accretion rate (dM/dt)_BH powers the accretion luminosity L_acc ~ (dM/dt)_BH, which drives black hole feedback. Using a series of controlled numerical experiments, we demonstrate that our new accretion disc particle method is more physically self-consistent than the Bondi-Hoyle method. We also discuss the physical implications of the accretion disc particle method for systems with a high degree of rotational support, and we argue that the M_BH-sigma relation in these systems should be offset from the relation for classical bulges and ellipticals, as appears to be observed.

With the ever increasing resolution of N-body simulations, accurate subhalo detection is becoming essential in the study of the formation of structure, the production of merger trees and the seeding of semi-analytic models. To investigate the state of halo finders, we compare two different approaches to detecting subhaloes; the first based on overdensities in a halo and the second being adaptive mesh refinement. A set of stable mock NFW dark matter haloes were produced and a subhalo was placed at different radii within a larger halo. SUBFIND (a Friends-of-Friends based finder) and AHF (an adaptive mesh based finder) were employed to recover the subhalo. As expected, we found that the mass of the subhalo recovered by SUBFIND has a strong dependence on the radial position and that neither halo finder can accurately recover the subhalo when it is very near the centre of the halo. This radial dependence is shown to be related to the subhalo being truncated by the background density of the halo and originates due to the subhalo being defined as an overdensity. If the subhalo size is instead determined using the peak of the circular velocity profile, a much more stable value is recovered. The downside to this is that the maximum circular velocity is a poor measure of stripping and is affected by resolution. For future halo finders to recover all the particles in a subhalo, a search of phase space will need to be introduced.

The distribution of cold gas in dark matter haloes is driven by key processes in galaxy formation: gas cooling, galaxy mergers, star formation and reheating of gas by supernovae. We compare the predictions of four different galaxy formation models for the spatial distribution of cold gas. We find that satellite galaxies make little contribution to the abundance or clustering strength of cold gas selected samples, and are far less important than they are in optically selected samples. The halo occupation distribution function of present-day central galaxies with cold gas mass > 10^9 h^-1 Msun is peaked around a halo mass of ~ 10^11 h^-1 Msun, a scale that is set by the AGN suppression of gas cooling. The model predictions for the projected correlation function are in good agreement with measurements from the HI Parkes All-Sky Survey. We compare the effective volume of possible surveys with the Square Kilometre Array with those expected for a redshift survey in the near-infrared. Future redshift surveys using neutral hydrogen emission will be competitive with the most ambitious spectroscopic surveys planned in the near-infrared.

(Abridged) Accurately predicting how the cosmic abundance of neutral hydrogen evolves with redshift is a challenging problem facing modellers of galaxy formation. We investigate the predictions of four currently favoured semi-analytical galaxy formation models applied to the Millennium simulation for the mass function of cold neutral gas (atomic and molecular) in galaxies as a function of redshift, and we use these predictions to construct number counts for the next generation of all-sky neutral atomic hydrogen (HI) surveys. Despite the different implementations of the physical ingredients of galaxy formation, we find that the model predictions are broadly consistent with one another; the key differences reflect how the models treat AGN feedback and how the timescale for star formation evolves with redshift. The models produce mass functions of cold gas in galaxies that are generally in good agreement with HI surveys at $z$=0. Interestingly we find that these mass functions do not evolve significantly with redshift. Adopting a simple conversion factor for cold gas mass to HI mass that we apply to all galaxies at all redshifts, we derive mass functions of HI in galaxies from the predicted mass functions of cold gas, which we use to predict the number counts of sources likely to be detected by HI surveys on next generation radio telescopes such as the Square Kilometre Array and its pathfinders. We show also how adopting a conversion factor for cold gas mass to HI mass that varies from galaxy to galaxy impacts on number counts. In addition, we examine how the typical angular sizes of galaxies vary with redshift. Taken together, these results make clear that forthcoming HI surveys will provide important and powerful tests of theoretical galaxy formation models.

The observed super-massive black hole (SMBH) mass -- galaxy velocity dispersion ($M_{\rm bh} - \sigma$) correlation may be established when winds/outflows from the SMBH drive gas out of the potential wells of classical bulges. Here we present numerical simulations of this process in a static isothermal potential. Simple spherically symmetric models of SMBH feedback at the Eddington luminosity can successfully explain the $M_{\rm bh} - \sigma$ and nuclear cluster mass $M_{\rm NC}-\sigma$ correlations, as well as why larger bulges host SMBHs while smaller ones host nuclear star clusters. However these models do not specify how SMBHs feed on infalling gas whilst simultaneously producing feedback that drives gas out of the galaxy. More complex models with rotation and/or anisotropic feedback allow SMBHs to feed via a disc or regions not exposed to SMBH winds, but in these more realistic cases it is not clear why a robust $M_{\rm bh} - \sigma$ relation should be established. In fact, some of the model predictions contradict observations. For example, an isotropic SMBH wind impacting on a disc (rather than a shell) of aspect ratio $H/R \ll 1$ requires the SMBH mass to be larger by a factor $\sim R/H$, which is opposite to what is observed. We conclude that understanding how a SMBH feeds is as important a piece of the puzzle as understanding how its feedback affects its host galaxy. Finally, we note that in aspherical cases the SMBH outflows induce differential motions in the bulge. This may pump turbulence that is known to hinder star formation in star forming regions. SMBH feedback thus may not only drive gas out of the bulge but also reduce the fraction of gas turned into stars.

We investigate the influence of the external tidal field of a dark matter halo on the dynamical evolution of star clusters using direct N-body simulations, where we assume that the halo is described by a Navarro, Frenk & White mass profile which has an inner density cusp. We assess how varying the mass and concentration of the halo affects the rate at which the star cluster loses mass and we find that increasing halo mass and concentration drives enhanced mass loss rates and in principle shorter cluster disruption timescales. In addition, we examine disruption timescales in a three-component model of a galaxy (bulge, disk and dark matter halo) and find good agreement with results based on an empirical model of the Galactic potential if we assume a halo mass of ~1e12 solar masses. In general, dark matter halos are expected to contribute significantly to the masses of galaxies and should not be ignored when modelling the evolution of star clusters. We extend our results to discuss how this can have a potentially profound effect on the disruption timescales of globular clusters, suggesting that we may underestimate the rate at which primordial globular clusters are disrupted.

Globular clusters are dense stellar systems that have typical ages of ~13 billion years, implying that they formed at redshifts of z>~6. Massive stars in newly formed or primordial globular clusters could have played an important role during the epoch of cosmological reionisation (z>~6) as sources of energetic, neutral hydrogen ionising UV photons. We investigate whether or not these stars could have been as important in death as sources of energetic X-ray photons as they were during their main sequence lives. Most massive stars are expected to form in binaries, and an appreciable fraction of these (as much as ~30%) will evolve into X-ray luminous (L_X~10^38 erg/s) high-mass X-ray binaries (HMXBs). These sources would have made a contribution to the X-ray background at z>~6. Using Monte Carlo models of a globular cluster, we estimate the total X-ray luminosity of a population of HMXBs. We compare and contrast this with the total UV luminosity of the massive stars during their main sequence lives. For reasonable estimates, we find that the bolometric luminosity of the cluster peaks at ~10^42 erg/s during the first few million years, but declines to ~10^41 erg/s after ~5 million years as the most massive stars evolve off the main sequence. From this time onwards, the total bolometric luminosity is dominated by HMXBs and falls gradually to ~10^40 erg/s after ~50 million years. Assuming a power-law spectral energy distribution for the HMXBs, we calculate the effective number of neutral hydrogen ionisations per HMXB and show that HMXBs can be as important as sources of ionising radiation as massive stars. Finally we discuss the implications of our results for modelling galaxy formation at high redshift and the prospects of using globular clusters as probes of reionisation.

Using a set of high-resolution dark matter only cosmological simulations we found a correlation between the dark matter halo mass M and its spin parameter lambda for objects forming at redshifts z > 10: the spin parameter decreases with increasing mass. However, halos forming at later times do not exhibit such a strong correlation, in agreement with the findings of previous studies. While we presented such a correlation in a previous study using the Bullock et al. (2001) spin parameter defintion we now defer to the classical definition showing that the results are independent of the definition.

We performed cosmological simulations based upon both a cold dark matter (CDM) and a warm dark matter (WDM) model. The focus of our investigations lies with selected spatial and kinematic properties of substructure halos (subhalos) orbiting within host halos, that form in both dark-matter cosmologies. We aim at using the dynamics of the subhalos as a probe of the respective cosmology.

[ABRIDGED VERSION] The future of cm and m-wave astronomy lies with the Square Kilometre Array (SKA), a telescope under development by a consortium of 17 countries. The SKA will be 50 times more sensitive than any existing radio facility. A majority of the key science for the SKA will be addressed through large-area imaging of the Universe at frequencies from 300 MHz to a few GHz. The Australian SKA Pathfinder (ASKAP) is aimed squarely in this frequency range, and achieves instantaneous wide-area imaging through the development and deployment of phase-array feed systems on parabolic reflectors. This large field-of-view makes ASKAP an unprecedented synoptic telescope poised to achieve substantial advances in SKA key science. The central core of ASKAP will be located at the Murchison Radio Observatory in inland Western Australia, one of the most radio-quiet locations on the Earth and one of the sites selected by the international community as a potential location for the SKA. Following an introductory description of ASKAP, this document contains 7 chapters describing specific science programmes for ASKAP. The combination of location, technological innovation and scientific program will ensure that ASKAP will be a world-leading radio astronomy facility, closely aligned with the scientific and technical direction of the SKA. A brief summary chapter emphasizes the point, and considers discovery space.

We identify satellites of isolated galaxies in SDSS and examine their angular distribution. Using mock catalogues generated from cosmological N-body simulations, we demonstrate that the selection criteria used to select isolated galaxies and their satellites must be very strict in order to correctly identify systems in which the primary galaxy dominates its environment. The criteria used in many previous studies instead select predominantly group members. We refine a set of selection criteria for which the group contamination is estimated to be less than 7% and present a catalogue of the resulting sample. The angular distribution of satellites about their host is biased towards the major axes for spheroidal galaxies and probably also for red disc galaxies, but is isotropic for blue disc galaxies, i.e. it is the colour of the host that determines the distribution of its satellites rather than its morphology. The similar anisotropy measured in this study as in studies that were dominated by groups implies that group-specific processes are not responsible for the angular distribution. Satellites that are most likely to have been recently accreted show a tendancy to lie along the same axis as the surrounding large scale structure. The orientations of isolated early and intermediate-type galaxies also align with the surrounding large scale structures. We discuss the origin of the anisotropic satellite distribution and consider the implications of our results, critically assessing the respective roles played by the orientation of the visible galaxy within its dark matter halo; anisotropic accretion of satellites from the larger scale environment; and the biased nature of satellites as tracers of the underlying dark matter subhalo population. (Abridged)

We explore the dependence of the radial alignment of subhalos on the mass of the host halo they orbit in. As the effect is seen on a broad range of scales including massive clusters as well as galactic systems it only appears natural to explore this phenomenon by means of cosmological simulations covering the same range in masses. We have 25 well resolved host dark matter halos at our disposal ranging from 10^15 Msun/h down to 10^12 Msun/h each consisting of order of a couple of million particles within the virial radius. We observe that subhalos tend to be more spherical than isolated objects. Both the distributions of sphericity and triaxiality of subhalos are Gaussian distributed with peak values of s approx. 0.80 and T approx. 0.56, irrespective of host mass. Interestingly we note that the radial alignment is independent of host halo mass and the distribution of \cos\theta (i.e. the angle between the major axis E_a of each subhalo and the radius vector of the subhalo in the reference frame of the host) is well fitted by a simple power law P(\cos\theta) proportional to \cos^4\theta with the same fitting parameters for all host halos.

We present a comparison of the properties of substructure halos (subhalos) orbiting within host halos that form in Cold Dark Matter (CDM) and Warm Dark Matter (WDM) cosmologies. Our study focuses on selected properties of these subhalos, namely their anisotropic spatial distribution within the hosts; the existence of a "backsplash'' population; the age-distance relation; the degree to which they suffer mass loss; and the distribution of relative (infall) velocities with respect to the hosts. We find that the number density of subhalos in our WDM model is suppressed relative to that in the CDM model, as we would expect. Interestingly, our analysis reveals that backsplash subhalos exist in both the WDM and CDM models. Indeed, there are no statistically significant differences between the spatial distributions of subhalos in the CDM and WDM models. There is evidence that subhalos in the WDM model suffer enhanced mass loss relative to their counterparts in the CDM model, reflecting their lower central densities. We note also a tendency for the (infall) velocities of subhalos in the WDM model to be higher than in the CDM model. Nevertheless, we conclude that observational tests based on either the spatial distribution or the kinematics of the subhalo population are unlikely to help us to differentiate between the CDM model and our adopted WDM model.

We report on a correlation between virial mass M and spin parameter lambda for dark matter halos forming at redshifts z > 10. We find that the spin parameter decreases with increasing halo mass. Interestingly, our analysis indicates that halos forming at later times do not exhibit such a strong correlation, in agreement with the findings of previous studies. We briefly discuss the implications of this correlation for galaxy formation at high redshifts and the galaxy population we observe today.

(abridged) We present a detailed analysis of the properties of tidally stripped material from disrupting substructure haloes or subhaloes in a sample of high resolution cosmological N-body host haloes ranging from galaxy- to cluster-mass scales. We focus on devising methods to recover the infall mass and infall eccentricity of subhaloes from the properties of their tidally stripped material (i.e. tidal streams). Our analysis reveals that there is a relation between the scatter of stream particles about the best-fit debris plane and the infall mass of the progenitor subhalo. This allows us to reconstruct the infall mass from the spread of its tidal debris in space. We also find that the spread in radial velocities of the debris material (as measured by an observer located at the centre of the host) correlates with the infall eccentricity of the subhalo, which allows us to reconstruct its orbital parameters. We devise an automated method to identify leading and trailing arms that can, in principle at least, be applied to observations of stellar streams from satellite galaxies. This method is based on the energy distribution of material in the tidal stream. Using this method, we show that the mass associated with leading and trailing arms differ. While our analysis indicates that tidal streams can be used to recover certain properties of their progenitor subhaloes (and consequently satellites), we do not find strong correlations between host halo properties and stream properties. This likely reflects the complicated relationship between the stream and the host, which in a cosmological context is characterised by a complex mass accretion history, an asymmetric mass distribution and the abundance of substructure. [...]

We explore the dependence of the central logarithmic slope of dark matter halo density profiles $\alpha$ on the spectral index $n$ of the linear matter power spectrum $P(k)$ using cosmological $N$-body simulations of scale-free models (i.e. $P(k) \propto k^n$). For each of our simulations we identify samples of well resolved haloes in dynamical equilibrium and we analyse their mass profiles. By parameterising the mass profile using a ``generalised'' Navarro, Frenk & White profile in which the central logarithmic slope $\alpha$ is allowed to vary while preserving the $r^{-3}$ asymptotic form at large radii, we obtain preferred central slopes for haloes in each of our models. There is a strong correlation between $\alpha$ and $n$, such that $\alpha$ becomes shallower as $n$ becomes steeper. However, if we normalise our mass profiles by $r_{-2}$, the radius at which the logarithmic slope of the density profile is -2, we find that these differences are no longer present. We conclude that there is no evidence for convergence to a unique central asymptotic slope, at least on the scales that we can resolve.

We present a precise estimate of the bulk virial scaling relation of halos formed via hierarchical clustering in an ensemble of simulated cold dark matter cosmologies. The result is insensitive to cosmological parameters, the presence of a trace, dissipationless gas component, and numerical resolution down to a limit of ~1000 particles. The dark matter velocity dispersion scales with total mass as log(sigma_{DM}(M,z)) = log(1082.9 +- 4.0 \kms) + (0.3361 +- 0.0026) log(h(z)M_{200}/10^{15} Msun), with h(z) the dimensionless Hubble parameter. At fixed mass, the velocity dispersion likelihood is nearly log-normal, with scatter sigma_{ln sigma} = 0.0426 +- 0.015, except for a tail to higher dispersions containing 10% of the population that are merger transients. We combine this relation with the halo mass function in LCDM models, and show that a low normalization condition, S_8 = sigma_8(Omega_m/0.3)^{0.35} \= 0.69, favored by recent WMAP and SDSS analysis requires that galaxy and gas specific energies in rich clusters be 50% larger than that of the underlying dark matter. Such large energetic biases are in conflict with the current generation of direct simulations of cluster formation. A higher normalization, S_8 = 0.80, alleviates this tension and implies that the hot gas fraction within r_{500} is (0.71 +- 0.09) h_{70}^{-3/2} Omega_b/\Omega_m, a value consistent with recent Sunyaev-Zel'dovich observations.

Galactic disks in triaxial dark matter halos become deformed by the elliptical potential in the plane of the disk in such a way as to counteract the halo ellipticity. We develop a technique to calculate the equilibrium configuration of such a disk in the combined disk-halo potential, which is based on the method of Jog (2000) but accounts for the radial variation in both the halo potential and the disk ellipticity. This crucial ingredient results in qualitatively different behavior of the disk: the disk circularizes the potential at small radii, even for a reasonably low disk mass. This effect has important implications for proposals to reconcile cuspy halo density profiles with low surface brightness galaxy rotation curves using halo triaxiality. The disk ellipticities in our models are consistent with observational estimates based on two-dimensional velocity fields and isophotal axis ratios.

We investigate the importance of interactions between dark matter substructures for the mass loss they suffer whilst orbiting within a sample of high resolution galaxy cluster mass Cold Dark Matter haloes formed in cosmological N-body simulations. We have defined a quantitative measure that gauges the degree to which interactions are responsible for mass loss from substructures. This measure indicates that interactions are more prominent in younger systems when compared to older more relaxed systems. We show that this is due to the increased number of encounters a satellite experiences and a higher mass fraction in satellites. This is in spite of the uniformity in the distributions of relative distances and velocities of encounters between substructures within the different host systems in our sample. Using a simple model to relate the net force felt by a single satellite to the mass loss it suffers, we show that interactions with other satellites account for ~30% of the total mass loss experienced over its lifetime. The relation between the age of the host and the importance of interactions increases the scatter about this mean value from ~25% for the oldest to ~45% for the youngest system we have studied. We conclude that satellite interactions play a vital role in the evolution of substructure in dark matter halos and that a significant fraction of the tidally stripped material can be attributed to these interactions.

We investigate the impact finite simulation box size has on the structural and kinematic properties of Cold Dark Matter haloes forming in cosmological simulations. Our approach involves generating a single realisation of the initial power spectrum of density perturbations and studying how truncation of this power spectrum on scales larger than L_cut affects the structure of dark matter haloes at z=0. In particular, we have examined the cases of L_cut = f_cut L_box with f_cut=1 (i.e. no truncation), 1/2, 1/3 and 1/4. In common with previous studies, we find that the suppression of long wavelength perturbations reduces the strength of clustering, as measured by a suppression of the 2-point correlation function xi(r), and reduces the numbers of the most massive haloes, as reflected in the depletion of the high mass end of the mass function n(M). Interestingly, we find that truncation has little impact on the internal properties of haloes. The masses of high mass haloes decrease in a systematic manner as L_cut is reduced, but the distribution of concentrations is unaffected. On the other hand, the median spin parameter is ~50% lower in runs with f_cut<1. We argue that this is an imprint of the linear growth phase of the halo's angular momentum by tidal torquing, and that the absence of any measurable trend in concentration and the weak trend observed in halo shape reflect the importance of virialisation and complex mass accretion histories for these quantities. These results are of interest for studies that require high mass resolution and statistical samples of simulated haloes, such as simulations of the population of first stars. Our analysis shows that large-scale tidal fields have relatively little effect on the internal properties of Cold Dark Matter haloes and hence may be ignored in such studies.

We present the proceedings from a 2-day workshop held at Swinburne University on the 24th-25th of May 2005. The workshop participants highlighted current Australian research on both theoretical and observational aspects of galaxy groups. These proceedings include short 1-page summaries of a number of the talks presented at the workshop. The talks presented ranged from reconciling N-body simulations with observations, to the HI content of galaxies in groups and the existence of ``dark galaxies''. The formation and existence of ultra-compact dwarfs in groups, and a new supergroup in Eridanus were also discussed.

The 11 known satellite galaxies within 250 kpc of the Milky Way lie close to a great circle on the sky. We use high resolution N-body simulations of galactic dark matter halos to test if this remarkable property can be understood within the context of the cold dark matter cosmology. We construct halo merger trees from the simulations and use a semianalytic model to follow the formation of satellite galaxies. We find that in all 6 of our simulations, the 11 brightest satellites are indeed distributed along thin, disk-like structures analogous to that traced by the Milky Way's satellites. This is in sharp contrast to the overall distributions of dark matter in the halo and of subhalos within it which, although triaxial, are not highly aspherical. We find that the spatial distribution of satellites is significantly different from that of the most massive subhalos but is similar to that of the subset of subhalos that had the most massive progenitors at earlier times. The elongated disk-like structure delineated by the satellites has its long axis aligned with the major axis of the dark matter halo. We interpret our results as reflecting the preferential infall of satellites along the spines of a few filaments of the cosmic web.

We investigate whether a correlation exists between the concentration and spin parameters of dark matter haloes forming in LCDM N-body simulations. In particular, we focus on haloes with virial masses characteristic of the hosts of Low Surface Brightness (LSB) and High Surface Brightness (HSB) galaxies. Our analysis demonstrates that concentration and spin are anti-correlated. If we assume that a galaxy disk forms in one of these haloes from baryons that approximately conserved angular momentum during their dissipation into the rotationally supported disk, then it is possible to estimate the disk's central surface density. For an appropriate choice of mass-to-light ratio, we can place constraints on the central surface brightness of the disk and hence identify the analogues of the host haloes of LSB and HSB galaxies. We find that our LSB galaxy analogues occupy haloes that have lower concentrations than might be expected based on consideration of the mass-concentration relation for the LCDM cosmology. The distribution of concentrations peaks at c~6, in good agreement with observational data, although there are important differences between the shapes of the simulated and observationally inferred distributions. This suggests that LSB galaxies inhabit a biased subsample of the halo population, and that this bias may be an important ingredient for resolving the current debate regarding the structure of LSB galaxy dark matter haloes.

We use N-body hydrodynamical simulations to study the structure of disks in triaxial potentials resembling CDM halos. Our analysis focuses on the accuracy of the dark mass distribution inferred from rotation curves derived from simulated long-slit spectra. We consider a massless disk embedded in a halo with axis ratios of 0.5:0.6:1.0 and with its rotation axis aligned with the minor axis of the halo. Closed orbits for the gaseous particles deviate from coplanar circular symmetry, resulting in a variety of long-slit rotation curve shapes, depending on the orientation of the disk relative to the line of sight. Rotation curves may thus differ significantly from the spherically-averaged circular velocity profile of the dark matter halo. "Solid-body" rotation curves--typically interpreted as a signature of a constant density core in the dark matter distribution--are obtained about 25% of the time for random orientations although the dark matter follows the cuspy density profile proposed by Navarro, Frenk & White (NFW). We conclude that the discrepancies reported between the shape of the rotation curve of low surface brightness galaxies and the structure of CDM halos may be resolved once the complex effects of halo triaxiality on the dynamics of the gas component is properly taken into account.

We investigate the mass profile of LambdaCDM halos using a suite of numerical simulations spanning five decades in halo mass, from dwarf galaxies to rich galaxy clusters. Our analysis confirms the proposal of Navarro, Frenk & White (NFW) that the shape of LambdaCDM halo mass profiles differs strongly from a power law and depends little on mass. The logarithmic slope of the spherically-averaged density profile, as measured by beta=-dln(rho)/dln(r), decreases monotonically towards the center and becomes shallower than isothermal (beta<2) inside a characteristic radius, r_{-2}. Although the fitting formula proposed by NFW provides a reasonably good approximation to the density and circular velocity profiles of individual halos, systematic deviations from the best NFW fits are also noticeable. Inside r_{-2}, the profile of simulated halos gets shallower with radius more gradually than predicted and, as a result, NFW fits tend to underestimate the dark matter density in these regions. This discrepancy has been interpreted as indicating a steeply divergent cusp, but our results suggest a different interpretation. We use the density and enclosed mass at our innermost resolved radii to place strong constraints on beta_{0}: density cusps as steep as r^{-1.5} are inconsistent with most of our simulations, although beta_{0}=1 is still consistent with our data. Our density profiles show no sign of converging to a well-defined asymptotic inner power law. We propose a simple formula that reproduces the radial dependence of the slope better than the NFW profile, and so may minimize errors when extrapolating our results inward to radii not yet reliably probed by numerical simulations.

We use a set of high-resolution cosmological N-body simulations to investigate the inner mass profile of galaxy-sized cold dark matter (CDM) halos. These simulations extend the thorough numerical convergence study presented in Paper I of this series (Power et al. 2003), and demonstrate that the mass profile of CDM halos can be robustly estimated beyond a minimum converged radius of order r_conv ~ 1 kpc/h in our highest resolution runs. The density profiles of simulated halos become progressively shallow from the virial radius inwards, and show no sign of approaching a well-defined power-law behaviour near the centre. At r_conv, the logarithmic slope of the density profile is steeper than the asymptotic \rho \propto r^-1 expected from the formula proposed by Navarro, Frenk, and White (1996), but significantly shallower than the steeply divergent \rho \propto r^-1.5 cusp proposed by Moore et al. (1999). We perform a direct comparison of the spherically-averaged dark matter circular velocity (V_c) profiles with rotation curves of low surface brightness (LSB) galaxies from the samples of de Blok et al. (2001), de Blok and Bosma (2002), and Swaters et al. (2003). Most (about two-thirds) LSB galaxies in this dataset are roughly consistent with CDM halo V_c profiles. However, about one third of LSBs in these samples feature a sharp transition between the rising and flat part of the rotation curve that is not seen in the V_c profiles of CDM halos. This discrepancy has been interpreted as excluding the presence of cusps, but we argue that it might simply reflect the difference between circular velocity and gas rotation speed likely to arise in gaseous disks embedded within realistic, triaxial CDM halos.

We present a comprehensive set of convergence tests which explore the role of various numerical parameters on the equilibrium structure of a simulated dark matter halo. We report results obtained with two independent, state-of-the-art, multi-stepping, parallel N-body codes: PKDGRAV and GADGET. We find that convergent mass profiles can be obtained for suitable choices of the gravitational softening, timestep, force accuracy, initial redshift, and particle number. For softenings chosen so that particle discreteness effects are negligible, convergence in the circular velocity is obtained at radii where the following conditions are satisfied: (i) the timestep is much shorter than the local orbital timescale; (ii) accelerations do not exceed a characteristic acceleration imprinted by the gravitational softening; and (iii) enough particles are enclosed so that the collisional relaxation timescale is longer than the age of the universe. The most stringent requirement for convergence is typically that imposed on the particle number by the collisional relaxation criterion, which implies that in order to estimate accurate circular velocities at radii where the density contrast may reach $\sim 10^6$, the region must enclose of order 3000 particles (or more than a few times $10^6$ within the virial radius). Applying these criteria to a galaxy-sized $\Lambda$CDM halo, we find that the spherically-averaged density profile becomes progressively shallower from the virial radius inwards, reaching a logarithmic slope shallower than -1.2 at the innermost resolved point, $r \sim 0.005 r_{200}$, with little evidence for convergence to a power-law behaviour in the inner regions.