We use cosmological simulations of isolated Milky Way-mass galaxies, as well as Local Group analogues, to define the "edge" -- a caustic manifested in a drop in density or radial velocity -- of Galactic-sized haloes, both in dark matter and in stars. In the dark matter, we typically identify two caustics: the outermost caustic located at ~1.4r_200m corresponding to the "splashback" radius, and a second caustic located at ~0.6r_200m which likely corresponds to the edge of the virialized material which has completed at least two pericentric passages. The splashback radius is ill defined in Local Group type environments where the halos of the two galaxies overlap. However, the second caustic is less affected by the presence of a companion, and is a more useful definition for the boundary of the Milky Way halo. Curiously, the stellar distribution also has a clearly defined caustic, which, in most cases, coincides with the second caustic of the dark matter. This can be identified in both radial density and radial velocity profiles, and should be measurable in future observational programmes. Finally, we show that the second caustic can also be identified in the phase-space distribution of dwarf galaxies in the Local Group. Using the current dwarf galaxy population, we predict the edge of the Milky Way halo to be 292 +/- 61 kpc.

We determine the inner density profiles of massive galaxy clusters (M$_{200}$ > $5 \times 10^{14}$ M$_{\odot}$) in the Cluster-EAGLE (C-EAGLE) hydrodynamic simulations, and investigate whether the dark matter density profiles can be correctly estimated from a combination of mock stellar kinematical and gravitational lensing data. From fitting mock stellar kinematics and lensing data generated from the simulations, we find that the inner density slopes of both the total and the dark matter mass distributions can be inferred reasonably well. We compare the density slopes of C-EAGLE clusters with those derived by Newman et al. for 7 massive galaxy clusters in the local Universe. We find that the asymptotic best-fit inner slopes of "generalized" NFW (gNFW) profiles, ${\gamma}_{\rm gNFW}$, of the dark matter haloes of the C-EAGLE clusters are significantly steeper than those inferred by Newman et al. However, the mean mass-weighted dark matter density slopes of the simulated clusters are in good agreement with the Newman et al. estimates. We also find that the estimate of ${\gamma}_{\rm gNFW}$ is very sensitive to the constraints from weak lensing measurements in the outer parts of the cluster and a bias can lead to an underestimate of ${\gamma}_{\rm gNFW}$.

We present a cosmological hydrodynamical simulation of a 10^13 Msun galaxy group and its environment (out to 10 times the virial radius) carried out using the EAGLE model of galaxy formation. Exploiting a novel technique to increase the resolution of the dark matter calculation independently of that of the gas, the simulation resolves dark matter haloes and subhaloes of mass 5x10^6 Msun . It is therefore useful for studying the abundance and properties of the haloes and subhaloes targeted in strong lensing tests of the cold dark matter model. We estimate the halo and subhalo mass functions and discuss how they are affected both by the inclusion of baryons in the simulation and by the environment. We find that the halo and subhalo mass functions have lower amplitude in the hydrodynamical simulation than in its dark matter only counterpart. This reflects the reduced growth of haloes in the hydrodynamical simulation due to the early loss of gas by reionisation and galactic winds and, additionally, in the case of subhaloes, disruption by enhanced tidal effects within the host halo due to the presence of a massive central galaxy. The distribution of haloes is highly anisotropic reflecting the filamentary character of mass accretion onto the cluster. As a result, there is significant variation in the number of structures with viewing direction. The median number of structures near the centre of the halo, when viewed in projection, is reduced by a factor of two when baryons are included.

The existence of hot, accreted gaseous coronae around massive galaxies is a long-standing central prediction of galaxy formation models in the $\Lambda$CDM cosmology. While observations now confirm that extraplanar hot gas is present around late-type galaxies, the origin of the gas is uncertain with suggestions that galactic feedback could be the dominant source of energy powering the emission. We investigate the origin and X-ray properties of the hot gas that surrounds galaxies of halo mass, $(10^{11}-10^{14}) \mathrm{M}_\odot$, in the cosmological hydrodynamical EAGLE simulations. We find that the central X-ray emission, $\leq 0.10 R_{\mathrm{vir}}$, of halos of mass $\leq 10^{13} \mathrm{M}_\odot$ originates from gas heated by supernovae (SNe). However, beyond this region, a quasi-hydrostatic, accreted atmosphere dominates the X-ray emission in halos of mass $\geq 10^{12} \mathrm{M}_\odot$. We predict that a dependence on halo mass of the hot gas to dark matter mass fraction can significantly change the slope of the $L_{\mathrm{X}}-M_{\mathrm{vir}}$ relation (which is typically assumed to be $4/3$ for clusters) and we derive the scaling law appropriate to this case. As the gas fraction in halos increases with halo mass, we find a steeper slope for the $L_{\mathrm{X}}-M_{\mathrm{vir}}$ in lower mass halos, $\leq 10^{14} \mathrm{M}_\odot$. This varying gas fraction is driven by active galactic nuclei (AGN) feedback. We also identify the physical origin of the so-called "missing feedback" problem, the apparently low X-ray luminosities observed from high star-forming, low-mass galaxies. This is explained by the ejection of SNe-heated gas from the central regions of the halo.

We study the properties of the dark matter component of the radially anisotropic stellar population recently identified in the Gaia data, using magneto-hydrodynamical simulations of Milky Way-like halos from the Auriga project. We identify 10 simulated galaxies that approximately match the rotation curve and stellar mass of the Milky Way. Four of these have an anisotropic stellar population reminiscent of the Gaia structure. We find an anti-correlation between the dark matter mass fraction of this population in the Solar neighbourhood and its orbital anisotropy. We estimate the local dark matter density and velocity distribution for halos with and without the anisotropic stellar population, and use them to simulate the signals expected in future xenon and germanium direct detection experiments. We find that a generalized Maxwellian distribution fits the dark matter halo integrals of the Milky Way-like halos containing the radially anisotropic stellar population. For dark matter particle masses below approximately 10 GeV, direct detection exclusion limits for the simulated halos with the anisotropic stellar population show a mild shift towards smaller masses compared to the commonly adopted Standard Halo Model.

In the $\Lambda$CDM model of structure formation galactic haloes build up by accretion of mass and mergers of smaller subunits. The most recent massive merger event experienced by the Milky Way (MW) halo was the accretion of the Large Magellanic Cloud (LMC; which has a stellar mass of $\sim \; 10^9M_\odot$). Recent analyses of galactic stellar data from the \textit{Gaia} satellite have uncovered an earlier massive accretion event, the Gaia-Enceladus-Sausage (GES), which merged with the MW around 10 Gyr ago. Here, we use the EAGLE cosmological hydrodynamics simulation to study properties of simulated MW-mass haloes constrained to have accretion histories similar to that of the MW, specifically the recent accretion of an "LMC" galaxy and a "GES'' merger, with a quiescent period between the GES merger and the infall of the LMC (the "GES \& LMC" class). We find that $\sim 16$ per cent of MW-mass haloes have an LMC; $\sim 5$ per cent have a GES event and no further merger with an equally massive object since $z=1$; and only $0.65$ per cent belong to the LMC \& GES category. The progenitors of the MWs in this last category are much less massive than average at early times but eventually catch up with the mean. The LMC \& GES class of galaxies naturally end up in the "blue cloud" in the colour-magnitude diagram at $z=0$, tend to have a disc morphology and have a larger than average number of satellite galaxies.

We analyse systems analogous to the Milky Way (MW) in the EAGLE cosmological hydrodynamics simulation in order to deduce the likely structure of the MW's dark matter halo. We identify MW-mass haloes in the simulation whose satellite galaxies have similar kinematics and spatial distribution to those of the bright satellites of the MW, specifically systems in which the majority of the satellites (8 out of 11) have nearly co-planar orbits that are also perpendicular to the central stellar disc. We find that the normal to the common orbital plane of the co-planar satellites is well aligned with the minor axis of the host dark matter halo, with a median misalignment angle of only $17.3^\circ$. Based on this result, we infer that the minor axis of the Galactic dark matter halo points towards $(l,b)=(182^\circ,-2^\circ)$, with an angular uncertainty at the 68 and 95 percentile confidence levels of 22$^\circ$ and 43$^\circ$ respectively. Thus, the inferred minor axis of the MW halo lies in the plane of the stellar disc. The halo, however, is not homologous and its flattening and orientation vary with radius. The inner parts of the halo are rounder than the outer parts and well-aligned with the stellar disc (that is the minor axis of the halo is perpendicular to the disc). Further out, the halo twists and the minor axis changes direction by $90^\circ$. This twist occurs over a very narrow radial range and reflects variations in the filamentary network along which mass was accreted into the MW.

We use a compilation of disc galaxy rotation curves to assess the role of the luminous component ("baryons") in the rotation curve diversity problem. As in earlier work, we find that rotation curve shape correlates with baryonic surface density: high surface density galaxies have rapidly-rising rotation curves consistent with cuspy cold dark matter halos; slowly-rising rotation curves (characteristic of galaxies with inner mass deficits or "cores") occur only in low surface density galaxies. The correlation, however, seems too weak to be the main driver of the diversity. In addition, dwarf galaxies exhibit a clear trend, from "cuspy" systems where baryons are unimportant in the inner mass budget to "cored" galaxies where baryons actually dominate. This trend constrains the various scenarios proposed to explain the diversity, such as (i) baryonic inflows and outflows during galaxy formation; (ii) dark matter self-interactions; (iii) variations in the baryonic mass structure coupled to rotation velocities through the "mass discrepancy-acceleration relation" (MDAR); or (iv) non-circular motions in gaseous discs. Together with analytical modeling and cosmological hydrodynamical simulations, our analysis shows that each of these scenarios has promising features, but none seems to fully account for the observed diversity. The MDAR, in particular, is inconsistent with the observed trend between rotation curve shape and baryonic importance; either the trend is caused by systematic errors in the data or the MDAR does not apply. The origin of the dwarf galaxy rotation curve diversity and its relation to the structure of cold dark matter halos remains an open issue.

We present a publicly-available toolkit of flight-proven hardware and software to retrieve 5 TB of data or small physical samples from a stratospheric balloon platform. Before launch, a capsule is attached to the balloon, and rises with it. Upon remote command, the capsule is released and descends via parachute, continuously transmitting its location. Software to predict the trajectory can be used to select a safe but accessible landing site. We dropped two such capsules from the SuperBIT telescope, in September 2019. The capsules took ~37 minutes to descend from ~30 km altitude. They drifted 32 km and 19 km horizontally, but landed within 300 m and 600 m of their predicted landing sites. We found them easily, and successfully recovered the data. We welcome interest from other balloon teams for whom the technology would be useful.

We determine the Milky Way (MW) mass profile inferred from fitting physically motivated models to the Gaia DR2 Galactic rotation curve and other data. Using various hydrodynamical simulations of MW-mass haloes, we show that the presence of baryons induces a contraction of the dark matter (DM) distribution in the inner regions, r<20 kpc. We provide an analytic expression that relates the baryonic distribution to the change in the DM halo profile. For our galaxy, the contraction increases the enclosed DM halo mass by factors of roughly 1.3, 2 and 4 at radial distances of 20, 8 and 1 kpc, respectively compared to an uncontracted halo. Ignoring this contraction results in systematic biases in the inferred halo mass and concentration. We provide a best-fitting contracted NFW halo model to the MW rotation curve that matches the data very well. The best-fit has a DM halo mass, $M_{200}^{\rm DM}=0.97_{-0.19}^{+0.24}\times10^{12} M_\odot$, and concentration before baryon contraction of $9.4_{-2.6}^{+1.9}$, which lie close to the median halo mass--concentration relation predicted in $\Lambda$CDM. The inferred total mass, $M_{200}^{\rm total}=1.08_{-0.14}^{+0.20} \times 10^{12} M_\odot$, is in good agreement with recent measurements. The model gives a MW stellar mass of $5.04_{-0.52}^{+0.43}\times10^{10} M_\odot$ and infers that the DM density at the Solar position is $\rho_{\odot}^{\rm DM}=8.8_{-0.5}^{+0.5}\times10^{-3} M_\odot \rm{pc}^{-3}\equiv0.33_{-0.02}^{+0.02}~\rm{GeV}~\rm{cm}^{-3}$. The rotation curve data can also be fitted with an uncontracted NFW halo model, but with very different DM and stellar parameters. The observations prefer the physically motivated contracted NFW halo, but the measurement uncertainties are too large to rule out the uncontracted NFW halo.

We present a model for the formation of the first galaxies before and after the reionization of hydrogen in the early universe. In this model, galaxy formation can only take place in dark matter haloes whose mass exceeds a redshift-dependent critical value, which, before reionization, is equal (in the simplest case) to the mass at which atomic hydrogen cooling becomes effective and, after reionization, is equal to the mass above which gas cannot remain in hydrostatic equilibrium. We define the Halo Occupation Fraction (HOF) as the fraction of haloes that host a luminous galaxy as a function of halo mass. The HOF is established by the interplay between the evolution of the critical mass and the assembly history of haloes and depends on three factors: the minimum halo mass for galaxy formation before reionization, the redshift of reionization, and the intensity of the (evolving) external photoheating rate. Our fiducial model predicts a cutoff in the galaxy mass function at a present-day halo mass, $M_{200} \sim 3\times 10^{8} M_{\odot}$; 100\% occupation at $M_{200} > 5\times 10^9 M_{\odot}$; and a population of starless gaseous haloes of present-day mass in the range $10^{6} \lesssim M_{200} / M_{\odot}\lesssim 5\times 10^{9}$, in which the gas is in thermal equilibrium with the ultraviolet background radiation and in hydrostatic equilibrium in the gravitational potential of the halo. The transition between HOF = 0 and HOF=1 reflects the stochastic nature of halo mass growth. We explore how these characteristic masses vary with model assumptions and parameter values. The results of our model are in excellent agreement with cosmological hydrodynamic simulations of galaxy formation.

We study structure formation in a set of cosmological simulations to uncover the scales in the initial density field that gave rise to the formation of present-day structures. Our simulations share a common primordial power spectrum (here Lambda-CDM), but the introduction of hierarchical variations of the phase information allows us to systematically study the scales that determine the formation of structure at later times. We consider the variance in z=0 statistics such as the matter power spectrum and halo mass function. We also define a criterion for the existence of individual haloes across simulations, and determine what scales in the initial density field contain sufficient information for the non-linear formation of unique haloes. We study how the characteristics of individual haloes such as the mass and concentration, as well as the position and velocity, are affected by variations on different scales, and give scaling relations for haloes of different mass. Finally, we use the example of a cluster-mass halo to show how our hierarchical parametrisation of the initial density field can be used to create variants of particular objects. With properties such as mass, concentration, kinematics and substructure of haloes set on distinct and well-determined scales, and its unique ability to introduce variations localised in real space, our method is a powerful tool to study structure formation in cosmological simulations.

Cosmology now has a standard model - a remarkably simple description of the universe, its contents and its history. A symposium held last September in Cambridge, UK, gave this model a 'health check' and discussed fascinating questions that lie beyond it.

We use magneto-hydrodynamical simulations of Milky Way-mass haloes from the Auriga project to examine the properties of surviving and destroyed dwarf galaxies that are accreted by these haloes over cosmic time. We show that the combined luminosity function of surviving and destroyed dwarfs at infall is similar in the various Auriga haloes, and is dominated by the destroyed dwarfs. There is, however, a strong dependence on infall time: destroyed dwarfs have typically early infall times, $t_{infall}<6$ Gyr, whereas the majority of dwarfs accreted at $t_{infall}>10$ Gyr have survived to the present day. Because of their late infall the surviving satellites today had higher metallicites at infall than their destroyed counterparts of similar infall mass; the difference is even more pronounced for the present-day metallicites of satellites, many of which continue to form stars after infall. In agreement with previous work, we find that a small number of relatively massive destroyed dwarf galaxies dominate the mass of the stellar haloes. However, there is a significant radial dependence: while 90 per cent of the mass in the inner regions ($<\,20\,$kpc) is contributed, on average, by only 3 massive progenitors, the outer regions ($>\,100\,$kpc) typically have $\sim8$ main progenitors of relatively lower mass. Finally, we show that a few massive progenitors dominate the metallicity distribution of accreted stars, even at the metal poor end. Contrary to common assumptions in the literature, dwarf galaxies of mass $M_{*}<10^7 \, M_{\odot}$ make up less than 10 per cent of the accreted, metal poor stars ([Fe/H] $<\,-3$) in the inner $50\,$kpc.

The nature of the dark matter can affect the collapse time of dark matter haloes, and can therefore be imprinted in observables such as the stellar population ages and star formation histories of dwarf galaxies. In this paper we use high resolution hydrodynamical simulations of Local Group-analogue (LG) volumes in cold dark matter (CDM), sterile neutrino warm dark matter (WDM) and self-interacting dark matter (SIDM) models with the EAGLE galaxy formation code to study how galaxy formation times change with dark matter model. We are able to identify the same haloes in different simulations, since they share the same initial density field phases. We show that the stellar mass varies systematically with resolution by over a factor of two, in a manner that depends on the final stellar mass. The evolution of the stellar populations in SIDM is largely identical to that of CDM, but in WDM early star formation is instead suppressed. The time at which LG haloes can begin to form stars through atomic cooling is delayed by $\sim$200~Myr in WDM models compared to CDM. 70~per~cent of WDM haloes of mass $>10^{8}M_{\odot}$ collapse early enough to form stars before $z=6$, compared to 90~per~cent of CDM and SIDM galaxies. It will be necessary to measure stellar ages of old populations to a precision of better than 100~Myr, and to address degeneracies with the redshift of reionization, in order to use these observables to distinguish between dark matter models.

We study the orbital phase-space of dark matter (DM) halos in the AURIGA suite of cosmological hydrodynamics simulations of Milky Way analogues. We characterise halos by their spherical action distribution, a function of the specific angular momentum, and the radial action, of the DM particles. By comparing DM-only and hydrodynamical simulations of the same halos, we investigate the contraction of DM halos caused by the accumulation of baryons at the centre. We find a small systematic suppression of the radial action in the DM halos of the hydrodynamical simulations, suggesting that the commonly used adiabatic contraction approximation can result in an underestimate of the density by ~ 8%. We apply an iterative algorithm to contract the AURIGA DM halos given a baryon density profile and halo mass, recovering the true contracted DM profiles with an accuracy of ~15%, that reflects halo-to-halo variation. Using this algorithm, we infer the total mass profile of the Milky Way's contracted DM halo. We derive updated values for the key astrophysical inputs to DM direct detection experiments: the DM density and velocity distribution in the Solar neighbourhood.

Dark matter haloes are the basic units of all cosmic structure. They grew by gravitational amplification of weak initial density fluctuations that are still visible on large scales in the cosmic microwave background radiation. Galaxies formed within relatively massive haloes as gas cooled and condensed at their centres, but many hypotheses for the nature of dark matter imply that the halo population should extend to masses many orders of magnitude below those where galaxies can form. Here, we use a novel, multi-zoom technique to create the first consistent simulation of the formation of present-day haloes over the full mass range populated when dark matter is aWeakly Interacting Massive Particle (WIMP) of mass ~100 GeV. The simulation has a dynamic range of 30 orders of magnitude in mass, resolving the internal structure of hundreds of Earth-mass haloes just as well as that of hundreds of rich galaxy clusters. Remarkably, halo density profiles are universal over the entire mass range and are well described by simple two-parameter fitting formulae. Halo mass and concentration are tightly related in a way which depends on cosmology and on the nature of the dark matter. At fixed mass, concentration is independent of local environment for haloes less massive than those of typical galaxies. These results are important for predicting annihilation radiation signals from dark matter, since these should be dominated by contributions from the smallest structures.

A non-parametric higher-order Jeans analysis method, GravSphere, was recently used to constrain the density profiles of Local Group dwarf galaxies, classifying them into those that are more likely to have an inner dark matter cusp and those that are likely to have a core (Read et al.). In this work we test this method using 31 simulated galaxies, comparable to Fornax, selected from the APOSTLE suite of cosmological hydrodynamics simulations, which include CDM and Self-Interacting Dark Matter (SIDM) cosmologies. We find that the mass profiles of the simulated dwarfs are often, but not always, well recovered by GravSphere. The less successful cases may be identified using a chi-squared diagnostic. Although the uncertainties are large in the inner regions, the inferred mass profiles are unbiased and exhibit smaller scatter than comparable Jeans methods. We find that GravSphere recovers the density profiles of simulated dwarfs below the half-light radius and down to the resolution limit of our simulations with better than 10+-30 per cent accuracy, making it a promising Jeans-based approach for modelling dark matter distributions in dwarf galaxies.

We use N-body simulations to revisit the globular cluster (GC) ``timing problem'' in the Fornax dwarf spheroidal (dSph). In agreement with earlier work, we find that, due to dynamical friction, GCs sink to the center of dark matter halos with a cuspy inner density profile but ``stall'' at roughly 1/3 of the core radius ($r_{\rm core}$) in halos with constant-density cores. The timescales to sink or stall depend strongly on the mass of the GC and on the initial orbital radius, but are essentially the same for either cuspy (NFW) or cored halos normalized to have the same total mass within $r_{\rm core}$. Arguing against a cusp on the basis that GCs have not sunk to the center is thus no different from arguing against a core, unless all clusters are today at $\sim (1/3)\, r_{\rm core}$. This would imply a core radius exceeding $\sim 3$ kpc, much larger than seems plausible in any core-formation scenario. (The average projected distance of Fornax GCs is $\langle R_{\rm GC,Fnx}\rangle\sim 1$ kpc and its effective radius is $\sim 700$ pc.) A simpler explanation is that Fornax GCs have only been modestly affected by dynamical friction, as expected if clusters started orbiting at initial radii of order $\sim 1$-$2$ kpc, just outside Fornax's present-day half-light radius but well within the tidal radius imprinted by Galactic tides. This is not entirely unexpected. Fornax GCs are significantly older and more metal-poor than most Fornax stars, and such populations in dSphs tend to be more spatially extended than their younger and more metal-rich counterparts. Contrary to some earlier claims, our simulations further suggest that GCs do not truly ``stall'' at $\sim 0.3\, r_{\rm core}$, but rather continue decaying toward the center, albeit at reduced rates. We conclude that dismissing the presence of a cusp in Fornax based on the spatial distribution of its GC population is unwarranted.

We investigate the formation of ultra-diffuse galaxies (UDGs) using the Auriga high-resolution cosmological magneto-hydrodynamical simulations of Milky Way-sized galaxies. We identify a sample of $92$ UDGs in the simulations that match a wide range of observables such as sizes, central surface brightness, S\'{e}rsic indices, colors, spatial distribution and abundance. Auriga UDGs have dynamical masses similar to normal dwarfs. In the field, the key to their origin is a strong correlation present in low-mass dark matter haloes between galaxy size and halo spin parameter. Field UDGs form in dark matter haloes with larger spins compared to normal dwarfs in the field, in agreement with previous semi-analytical models. Satellite UDGs, on the other hand, have two different origins: $\sim 55\%$ of them formed as field UDGs before they were accreted; the remaining $\sim 45\%$ were normal field dwarfs that subsequently turned into UDGs as a result of tidal interactions.

We explore the predictive power of cosmological, hydrodynamical simulations for stellar phase space substructure and velocity correlations with the Auriga simulations and Aurigaia mock-Gaia catalogues. We show that at the solar circle the Auriga simulations commonly host phase space structures in the stellar component that have constant orbital energies and arise from accreted subhaloes. These structures can persist for a few Gyrs, even after coherent streams in position space have been erased. We also explore velocity two-point correlation functions and find this diagnostic is not deterministic for particular clustering patterns in phase space. Finally, we explore these structure diagnostics with the Aurigaia catalogues and show that the current catalogues have the ability to recover some structures in phase space but careful consideration is required to separate physical structures from numerical structures arising from catalogue generation methods.

Ultrafaint dwarf galaxies ($M_\star\lesssim10^{5-6}\,{\rm M}_\odot$) are relics of an early phase of galaxy formation. They contain some of the oldest and most metal-poor stars in the Universe which likely formed before the epoch of hydrogen reionisation. These galaxies are so faint that they can only be detected as satellites of the Milky Way. They are so small that they are just barely resolved in current cosmological hydrodynamics simulations. Here we combine very high resolution cosmological $N$-body simulations with a semi-analytic model of galaxy formation to study the demographics and spatial distribution of ultrafaint satellites in Milky Way-mass haloes. We show that the abundance of these galaxies is correlated with the assembly history of the host halo: at fixed mass, haloes assembled earlier contain, on average, more ultrafaint satellites today than haloes assembled later. We identify simulated galactic haloes that experience an ancient Gaia-Enceladus-Sausage-like and a recent LMC-like accretion event and find that the former occurs in 33% of the sample and the latter in 9%. Only 3% experience both events and these are especially rich in ultrafaint satellites, most acquired during the ancient accretion event. Our models predict that the radial distribution of satellites is more centrally concentrated in early-forming haloes. Accounting for the depletion of satellites by tidal interactions with the central disc, we find a very good match to the observed radial distribution of satellites in the Milky Way over the entire radial range. This agreement is mainly due to the ability of our model to track 'orphan' galaxies after their subhaloes fall below the resolution limit of the simulation.

We examine the formation of dark matter (DM) cores in dwarf galaxies simulated with the EAGLE model of galaxy formation. As in earlier work, we find that the star formation (SF) gas density threshold ($\rho_{\rm th}$) plays a critical role. At low thresholds (LT), gas is unable to reach densities high enough to dominate the gravitational potential before being dispersed by feedback from supernovae. LT runs show little effect on the inner DM profile, even in systems with extended and bursty SF, two ingredients often cited as critical for core formation. For higher thresholds, gas is able to dominate the gravitational potential before being ejected by feedback. This can lead to a substantial reduction in the inner DM content, but only if the gas is gravitationally important over an extended period of time, allowing the halo to contract before gas removal. Rapid assembly and removal of gas in short SF bursts is less effective at altering the inner DM content. Subsequent gas accretion may draw DM back in and reform a cusp, unless SF is bursty enough to prevent it, preserving the core. Thus, for the EAGLE SF+feedback model, there is no simple relation between core formation and SF history, contrary to recent claims. The dependence of the inner DM content of dwarfs on $\rho_{\rm th}$ hinders robust predictions and the interpretation of observations. A simulation of a $(12 \rm \ Mpc)^3$ volume with high $\rho_{\rm th}$ results in dwarfs with sizeable cores over a limited halo mass range, but with insufficient variety in mass profiles to explain the observed diversity of dwarf galaxy rotation curves.

The observed Lyman-$\alpha$ flux power spectrum (FPS) is suppressed on scales below $\sim~ 30~{\rm km~s}^{-1}$. This cutoff could be due to the high temperature, $T_0$, and pressure, $p_0$, of the absorbing gas or, alternatively, it could reflect the free streaming of dark matter particles in the early universe. We perform a set of very high resolution cosmological hydrodynamic simulations in which we vary $T_0$, $p_0$ and the amplitude of the dark matter free streaming, and compare the FPS of mock spectra to the data. We show that the location of the dark matter free-streaming cutoff scales differently with redshift than the cutoff produced by thermal effects and is more pronounced at higher redshift. We, therefore, focus on a comparison to the observed FPS at $z>5$. We demonstrate that the FPS cutoff can be fit assuming cold dark matter, but it can be equally well fit assuming that the dark matter consists of $\sim 7$ keV sterile neutrinos in which case the cutoff is due primarily to the dark matter free streaming.

We analyse the orbital kinematics of the Milky Way (MW) satellite system utilizing the latest systemic proper motions for 38 satellites based on data from Gaia Data Release 2. Combining these data with distance and line-of-sight velocity measurements from the literature, we use a likelihood method to model the velocity anisotropy, $\beta$, as a function of Galactocentric distance and compare the MW satellite system with those of simulated MW-mass haloes from the APOSTLE and Auriga simulation suites. The anisotropy profile for the MW satellite system increases from $\beta\sim -2$ at $r\sim20$ kpc to $\beta\sim 0.5$ at $r\sim200$ kpc, indicating that satellites closer to the Galactic centre have tangentially-biased motions while those farther out have radially-biased motions. The motions of satellites around APOSTLE host galaxies are nearly isotropic at all radii, while the $\beta(r)$ profiles for satellite systems in the Auriga suite, whose host galaxies are substantially more massive in baryons than those in APOSTLE, are more consistent with that of the MW satellite system. This shape of the $\beta(r)$ profile may be attributed to the central stellar disc preferentially destroying satellites on radial orbits, or intrinsic processes from the formation of the Milky Way system.

The development of methods and algorithms to solve the $N$-body problem for classical, collisionless, non-relativistic particles has made it possible to follow the growth and evolution of cosmic dark matter structures over most of the Universe's history. In the best studied case $-$ the cold dark matter or CDM model $-$ the dark matter is assumed to consist of elementary particles that had negligible thermal velocities at early times. Progress over the past three decades has led to a nearly complete description of the assembly, structure and spatial distribution of dark matter haloes, and their substructure in this model, over almost the entire mass range of astronomical objects. On scales of galaxies and above, predictions from this standard CDM model have been shown to provide a remarkably good match to a wide variety of astronomical data over a large range of epochs, from the temperature structure of the cosmic background radiation to the large-scale distribution of galaxies. The frontier in this field has shifted to the relatively unexplored subgalactic scales, the domain of the central regions of massive haloes, and that of low-mass haloes and subhaloes, where potentially fundamental questions remain. Answering them may require: (i) the effect of known but uncertain baryonic processes (involving gas and stars), and/or (ii) alternative models with new dark matter physics. Here we present a review of the field, focusing on our current understanding of dark matter structure from $N$-body simulations and on the challenges ahead.

We present the status of the Dark Energy Spectroscopic Instrument (DESI) and its plans and opportunities for the coming decade. DESI construction and its initial five years of operations are an approved experiment of the US Department of Energy and is summarized here as context for the Astro2020 panel. Beyond 2025, DESI will require new funding to continue operations. We expect that DESI will remain one of the world's best facilities for wide-field spectroscopy throughout the decade. More about the DESI instrument and survey can be found at https://www.desi.lbl.gov.

A number of Local Group dwarf galaxies are known to have two spatially segregated stellar metallicity populations, a centrally concentrated metal-rich population and a more extended metal-poor population. In this work we discuss mechanisms that lead to the formation of two spatially segregated metallicity populations. Using a set of high-resolution hydrodynamical simulations of Local Group-like environments, we select a sample of satellite and field galaxies, spanning the stellar mass range $10^6-10^9$M$_{\odot}$, that exhibit bimodality in their metallicity distributions. Among those, we identify a subsample with a strong spatial segregation in the two populations. We find three distinct mechanisms for their formation. In field dwarfs and in a small fraction of satellites, a merger causes the metal-poor stars to migrate to larger radii and encourages the available gas to sink to the centre of the dwarf. Most of the gas is subsequently blown out of the halo through star formation feedback, but the remaining gas is consumed in the formation of a metal-rich population. In the exclusive case of satellites that have retained some of their gas at infall, it is the compression of this gas by ram pressure near pericentre that triggers the formation of metal-rich stars, whilst simultaneously preventing star formation at larger radii through stripping. Additionally, in a small number of field and satellite dwarfs, interactions with gaseous filaments and other galaxies can result in the formation of a metal-rich population. Regardless of the formation mechanism, a history of mergers typically enhances the spatial segregation.

N-body simulations make unambiguous predictions for the abundance of substructures within dark matter halos. However, the inclusion of baryons in the simulations changes the picture because processes associated with the presence of a large galaxy in the halo can destroy subhalos and substantially alter the mass function and velocity distribution of subhalos. We compare the effect of galaxy formation on subhalo populations in two state-of-the-art sets of hydrodynamical CDM simulations of Milky Way mass halos, APOSTLE and AURIGA. We introduce a new method for tracking the orbits of subhalos between simulation snapshots that gives accurate results down to a few kiloparsecs from the centre of the halo. Relative to a dark matter-only simulation, the abundance of subhalos in APOSTLE is reduced by 50% near the centre and by 10% within r200. In AURIGA the corresponding numbers are 80% and 40%. The velocity distributions of subhalos are also affected by the presence of the galaxy, much more so in AURIGA than in APOSTLE . The differences on subhalo properties in the two simulations can be traced back to the mass of the central galaxies, which in AURIGA are typically twice as massive as those in APOSTLE . We show that some of the results from previous studies are inaccurate due to systematic errors in the modelling of subhalo orbits near the centre of halos.

We study the formation of planes of dwarf galaxies around Milky Way (MW)-mass haloes in the EAGLE galaxy formation simulation. We focus on satellite systems similar to the one in the MW: spatially thin or with a large fraction of members orbiting in the same plane. To characterise the latter, we introduce a robust method to identify the subsets of satellites that have the most co-planar orbits. Out of the 11 MW classical dwarf satellites, 8 have highly clustered orbital planes whose poles are contained within a $22^\circ$ opening angle centred around $(l,b)=(182^\circ,-2^\circ)$. This configuration stands out when compared to both isotropic and typical $\Lambda$CDM satellite distributions. Purely flattened satellite systems are short-lived chance associations and persist for less than $1~\rm{Gyr}$. In contrast, satellite subsets that share roughly the same orbital plane are longer lived, with half of the MW-like systems being at least $4~\rm{Gyrs}$ old. On average, satellite systems were flatter in the past, with a minimum in their minor-to-major axes ratio about $9~\rm{Gyrs}$ ago, which is the typical infall time of the classical satellites. MW-like satellite distributions have on average always been flatter than the overall population of satellites in MW-mass haloes and, in particular, they correspond to systems with a high degree of anisotropic accretion of satellites. We also show that torques induced by the aspherical mass distribution of the host halo channel some satellite orbits into the host's equatorial plane, enhancing the fraction of satellites with co-planar orbits. In fact, the orbital poles of co-planar satellites are tightly aligned with the minor axis of the host halo.

We study the Local Group (LG) dwarf galaxy population using cosmological simulations of volumes selected to resemble the surroundings of the Milky Way (MW) and Andromeda (M31) galaxies. The simulations, from the APOSTLE project, indicate that the total mass within 3 Mpc of the pair's midpoint ($M_{\rm 3Mpc}$) typically exceeds $\sim 3$ times the sum of the virial masses of the two primaries and that the dwarf galaxy formation efficiency per unit mass is uniform throughout the volume. This suggests that MW and M31 satellites should make up fewer than one third of all LG dwarfs within 3 Mpc. This is consistent with the fraction (28 per cent) of galaxies with stellar mass $M_*>10^7\,{\rm M_{\odot}})$ that are satellites. In addition, the total number of such galaxies (42) suggests, for the APOSTLE galaxy mass-halo mass relation, a total LG mass of $M_{\rm 3Mpc}\sim 10^{13} \,{\rm M_{\odot}}$. At lower galaxy masses, however, the satellite fraction is substantially higher (42 per cent for $M_*>10^5\,{\rm M_{\odot}}$, or 42 out of 99). If this is due to incompleteness in the field sample, then $\sim 50$ dwarf galaxies at least as massive as the Draco dwarf spheroidal are missing from our current LG field dwarf inventory. The incompleteness interpretation is supported by the pronounced flattening of the LG luminosity function below $M_*\sim 10^7\,{\rm M_{\odot}}$, and by the scarcity of low-surface brightness LG field galaxies: whereas $\sim 40$ per cent of all $M_*>10^5\,{\rm M_{\odot}}$ MW/M31 satellites have effective surface brightness below $\Sigma_{\rm eff}=26.5$ m$_V$/arcsec$^2$, only 1 out of all field dwarfs does. The simulations indicate that most missing dwarfs should lie near the virial boundaries of the two LG primaries, and predict a trove of nearby dwarfs that await discovery by upcoming wide-field imaging surveys.

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

The dark matter velocity distribution in the Solar neighbourhood is an important astrophysical input which enters in the predicted event rate of dark matter direct detection experiments. It has been recently suggested that the local dark matter velocity distribution can be inferred from that of old or metal-poor stars in the Milky Way. We investigate this potential relation using six high resolution magneto-hydrodynamical simulations of Milky Way-like galaxies of the Auriga project. We do not find any correlation between the velocity distributions of dark matter and old stars in the Solar neighbourhood. Likewise, there are no strong correlations between the local velocity distributions of dark matter and metal-poor stars selected by applying reasonable cuts on metallicity. In some simulated galaxies, extremely metal-poor stars have a velocity distribution that is statistically consistent with that of the dark matter, but the sample of such stars is so small that we cannot draw any strong conclusions.

Measurements of the rotation curves of dwarf galaxies are often interpreted as requiring a constant density core at the centre, at odds with the "cuspy" inner profiles predicted by $N$-body simulations of cold dark matter (CDM) haloes. It has been suggested that this conflict could be resolved by fluctuations in the inner gravitational potential caused by the periodic removal of gas following bursts of star formation. Earlier work has suggested that core formation requires a bursty and extended star formation history (SFH). Here we investigate the structure of CDM haloes of dwarf galaxies ($M_{{\rm DM}} \sim 10^9-5\times10^{10}\,{\rm M}_\odot$) formed in the APOSTLE ('A Project of Simulating the Local Environment') and AURIGA cosmological hydrodynamic simulations. Our simulations have comparable or better resolution than others that make cores ($M_{{\rm gas}} \sim 10^4\,{\rm M}_\odot$, gravitational softening $\sim 150$ pc). Yet, we do not find evidence of core formation at {\it any} mass or any correlation between the inner slope of the DM density profile and temporal variations in the SFH. APOSTLE and AURIGA dwarfs display a similar diversity in their cumulative SFHs to available data for Local Group dwarfs. Dwarfs in both simulations are DM-dominated on all resolved scales at all times, likely limiting the ability of gas outflows to alter significantly the central density profiles of their haloes. We conclude that recurrent bursts of star formation are not sufficient to cause the formation of cores, and that other conditions must also be met for baryons to be able to modify the central DM cusp.

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

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.

The nearby Perseus galaxy cluster is a key target for indirect detection searches for decaying dark matter. We use the C-EAGLE simulations of galaxy clusters to predict the flux, width and shape of a dark matter decay line, paying particular attention to the unexplained 3.55keV line detected in the spectra of some galaxies and clusters, and the upcoming XRISM X-ray observatory mission. We show that the line width in C-EAGLE clusters similar to Perseus is typically [600-800]$\mathrm{kms^{-1}}$, and therefore narrower than the amplitude of the velocity dispersion of galaxies in the cluster. Haloes that are significantly disturbed can, however, exhibit galaxy velocity dispersions higher than $1000\mathrm{kms^{-1}}$, and in this case will show a large difference between the line profiles of on- and off-center observations. We show that the line profile is likely to be slightly asymmetric, but still well approximated by a Gaussian at the 10% level, and that the halo asymmetry can lead to fluxes that vary by a factor of two. In summary, we predict that, if the previously reported 3.55keV line detections do originate from dark matter decay, the XRISM mission will detect a line with a roughly Gaussian profile at a rest frame energy of 3.55keV, with a width $>600\mathrm{kms^{-1}}$ and flux approximately in the range $[4-9]\times10^{-8}\mathrm{counts/sec/cm^{2}}$.

We investigate the alignment of galaxies and haloes relative to cosmic web filaments using the EAGLE hydrodynamical simulation. We identify filaments by applying the NEXUS+ method to the mass distribution and the Bisous formalism to the galaxy distribution. Both web finders return similar filamentary structures that are well aligned and that contain comparable galaxy populations. EAGLE haloes have an identical spin alignment with filaments as their counterparts in dark matter only simulations: a complex mass dependent trend with low mass haloes spinning preferentially parallel to and high mass haloes spinning preferentially perpendicular to filaments. In contrast, galaxy spins do not show such a spin transition and have a propensity for perpendicular alignments at all masses, with the degree of alignment being largest for massive galaxies. This result is valid for both NEXUS+ and Bisous filaments. When splitting by morphology, we find that elliptical galaxies show a stronger orthogonal spin--filament alignment than spiral galaxies of similar mass. The same is true of their haloes, with the host haloes of elliptical galaxies having a larger degree of orthogonal alignment than the host haloes of spirals. Due to the misalignment between galaxy shape and spin, galaxy minor axes are oriented differently with filaments than galaxy spins. We find that the galaxies whose minor axis is perpendicular to a filament are much better aligned with their host haloes. This suggests that many of the same physical processes determine both the galaxy--filament and the galaxy--halo alignments.

We compare the mass cooling rates and cumulative cooled-down masses predicted by several semi-analytical (SA) cooling models with cosmological hydrodynamical simulations performed using the AREPO code (ignoring processes such as feedback and chemical enrichment). The SA cooling models are the new GALFORM cooling model introduced in Hou et al. (2017), along with two earlier GALFORM cooling models and the L-GALAXIES and MORGANA cooling models. We find that the predictions of the new GALFORM cooling model are generally in best agreement with the simulations. For halos with $M_{\rm halo}\lesssim 3\times 10^{11}\,{\rm M}_{\odot}$, the SA models predict that the timescale for radiative cooling is shorter than or comparable to the gravitational infall timescale. Even though SA models assume that gas falls onto galaxies from a spherical gas halo, while the simulations show that the cold gas is accreted through filaments, both methods predict similar mass cooling rates, because in both cases the gas accretion occurs on similar timescales. For halos with $M_{\rm halo}\gtrsim 10^{12}\,{\rm M}_{\odot}$, gas in the simulations typically cools from a roughly spherical hot gas halo, as assumed in the SA models, but the halo gas gradually contracts during cooling, leading to compressional heating. SA models ignore this heating, and so overestimate mass cooling rates by factors of a few. At low redshifts halo major mergers or a sequence of successive smaller mergers are seen in the simulations to strongly heat the halo gas and suppress cooling, while mergers at high redshifts do not suppress cooling, because the gas filaments are difficult to heat up. The new SA cooling model best captures these effects.

We use the APOSTLE and Auriga cosmological simulations to study the star formation histories (SFHs) of field and satellite dwarf galaxies. Despite sizeable galaxy-to-galaxy scatter, the SFHs of APOSTLE and Auriga dwarfs exhibit robust average trends with galaxy stellar mass: faint field dwarfs ($10^5<M_{\rm star}/M_\odot<10^{6.5}$) have, on average, steadily declining SFHs, whereas brighter dwarfs ($10^{7.5}<M_{\rm star}/M_\odot<10^{9}$) show the opposite trend. Intermediate-mass dwarfs have roughly constant SFHs. Satellites exhibit similar average trends, but with substantially suppressed star formation in the most recent $\sim 5$ Gyr, likely as a result of gas loss due to tidal and ram-pressure stripping after entering the haloes of their primaries. These simple mass and environmental trends are in good agreement with the derived SFHs of Local Group (LG) dwarfs whose photometry reaches the oldest main sequence turnoff. SFHs of galaxies with less deep data show deviations from these trends, but this may be explained, at least in part, by the large galaxy-to-galaxy scatter, the limited sample size, and the large uncertainties of the inferred SFHs. Confirming the predicted mass and environmental trends will require deeper photometric data than currently available, especially for isolated dwarfs.

We examine the stellar haloes of the Auriga simulations, a suite of thirty cosmological magneto-hydrodynamical high-resolution simulations of Milky Way-mass galaxies performed with the moving-mesh code AREPO. We study halo global properties and radial profiles out to $\sim 150$ kpc for each individual galaxy. The Auriga haloes are diverse in their masses and density profiles; mean metallicity and metallicity gradients; ages; and shapes, reflecting the stochasticity inherent in their accretion and merger histories. A comparison with observations of nearby late-type galaxies shows very good agreement between most observed and simulated halo properties. However, Auriga haloes are typically too massive. We find a connection between population gradients and mass assembly history: galaxies with few significant progenitors have more massive haloes, possess large negative halo metallicity gradients and steeper density profiles. The number of accreted galaxies, either disrupted or under disruption, that contribute 90% of the accreted halo mass ranges from 1 to 14, with a median of 6.5, and their stellar masses span over three orders of magnitude. The observed halo mass--metallicity relation is well reproduced by Auriga and is set by the stellar mass and metallicity of the dominant satellite contributors. This relationship is found not only for the accreted component but also for the total (accreted + in-situ) stellar halo. Our results highlight the potential of observable halo properties to infer the assembly history of galaxies.

Pan-STARRS1 has carried out a set of distinct synoptic imaging sky surveys including the $3\pi$ Steradian Survey and the Medium Deep Survey in 5 bands ($grizy_{P1}$). The mean 5$\sigma$ point source limiting sensitivities in the stacked 3$\pi$ Steradian Survey in $grizy_{P1}$ are (23.3, 23.2, 23.1, 22.3, 21.4) respectively. The upper bound on the systematic uncertainty in the photometric calibration across the sky is 7-12 millimag depending on the bandpass. The systematic uncertainty of the astrometric calibration using the Gaia frame comes from a comparison of the results with Gaia: the standard deviation of the mean and median residuals ($ \Delta ra, \Delta dec $) are (2.3, 1.7) milliarcsec, and (3.1, 4.8) milliarcsec respectively. The Pan-STARRS system and the design of the PS1 surveys are described and an overview of the resulting image and catalog data products and their basic characteristics are described together with a summary of important results. The images, reduced data products, and derived data products from the Pan-STARRS1 surveys are available to the community from the Mikulski Archive for Space Telescopes (MAST) at STScI.

We present and apply a method to infer the mass of the Milky Way (MW) by comparing the dynamics of MW satellites to those of model satellites in the EAGLE cosmological hydrodynamics simulations. A distribution function (DF) for galactic satellites is constructed from EAGLE using specific angular momentum and specific energy, which are scaled so as to be independent of host halo mass. In this 2-dimensional space, the orbital properties of satellite galaxies vary according to the host halo mass. The halo mass can be inferred by calculating the likelihood that the observed satellite population is drawn from this DF. Our method is robustly calibrated on mock EAGLE systems. We validate it by applying it to the completely independent suite of 30 AURIGA high-resolution simulations of MW-like galaxies: the method accurately recovers their true mass and associated uncertainties. We then apply it to ten classical satellites of the MW with 6D phase-space measurements, including updated proper motions from the GAIA satellite. The mass of the MW is estimated to be $M_{200}^{\textnormal{MW}}=1.17_{-0.15}^{+0.21}\times10^{12}M_{\odot}$ (68\% confidence limits). We combine our total mass estimate with recent mass estimates in the inner regions of the Galaxy to infer an inner dark matter (DM) mass fraction $M^\textnormal{DM}(<20~\rm{kpc})/M^\textnormal{DM}_{200}=0.12$ which is typical of ${\sim}10^{12}M_{\odot}$ $\Lambda$CDM haloes in hydrodynamical galaxy formation simulations. Assuming an NFW profile, this is equivalent to a halo concentration of $c_{200}^{\textnormal{MW}}=10.9^{+2.6}_{-2.0}$.

The Milky Way (MW) offers a uniquely detailed view of galactic structure and is often regarded as a prototypical spiral galaxy. But recent observations indicate that the MW is atypical: it has an undersized supermassive black hole at its centre; it is surrounded by a very low mass, excessively metal-poor stellar halo; and it has an unusually large nearby satellite galaxy, the Large Magellanic Cloud (LMC). Here we show that the LMC is on a collision course with the MW with which it will merge in $2.4^{+1.2}_{-0.8}$ Gyrs (68% confidence level). This catastrophic and long-overdue event will restore the MW to normality. Using the EAGLE galaxy formation simulation, we show that, as a result of the merger, the central supermassive black hole will increase in mass by up to a factor of 8. The Galactic stellar halo will undergo an equally impressive transformation, becoming 5 times more massive. The additional stars will come predominantly from the disrupted LMC, but a sizeable number will be ejected onto the halo from the stellar disc. The post-merger stellar halo will have the median metallicity of the LMC, [Fe/H]=-0.5 dex, which is typical of other galaxies of similar mass to the MW. At the end of this exceptional event, the MW will become a true benchmark for spiral galaxies, at least temporarily.

Deep observations of faint surface brightness stellar tidal streams in external galaxies with LSST are addressed in this White Paper contribution. We propose using the Wide--Fast--Deep survey that contains several nearby galaxies (at distances where the stars themselves are not resolved, i.e., beyond 20 Mpc). In the context of hierarchical galaxy formation, it is necessary to understand the prevalence and properties of tidal substructure around external galaxies based on integrated (i.e., unresolved) diffuse light. This requires collecting observations on much larger samples of galaxies than the Milky Way and M31. We will compare the observed structures to the predictions of cosmological models of galactic halo formation that inform us about the number and properties of streams around Milky Way-like galaxies. The insight gained from these comparisons will allow us to infer the properties of stream progenitors (masses, dynamics, metallicities, stellar populations). The changes in the host galaxies caused by the interactions with the dissolving companion galaxies will be another focus of our studies. We conclude by discussing synergies with WFIRST and Euclid, and also provide concrete suggestions for how the effects of scattered light could be minimized in LSST images to optimize the search for low surface brightness features, such as faint unresolved stellar tidal streams.

We use cosmological hydrodynamical simulations of the APOSTLE project along with high-quality rotation curve observations to examine the fraction of baryons in {\Lambda}CDM haloes that collect into galaxies. This 'galaxy formation efficiency' correlates strongly and with little scatter with halo mass, dropping steadily towards dwarf galaxies. The baryonic mass of a galaxy may thus be used to place a lower limit on total halo mass and, consequently, on its asymptotic maximum circular velocity. A number of observed dwarfs seem to violate this constraint, having baryonic masses up to ten times higher than expected from their rotation speeds, or, alternatively, rotating at only half the speed expected for their mass. Taking the data at face value, either these systems have formed galaxies with extraordinary efficiency - highly unlikely given their shallow potential wells - or their dark matter content is much lower than expected from {\Lambda}CDM haloes. This 'missing dark matter' is reminiscent of the inner mass deficit of galaxies with slowly-rising rotation curves, but cannot be explained away by star formation-induced 'cores' in the dark mass profile, since the anomalous deficit applies to regions larger than the luminous galaxies themselves. We argue that explaining the structure of these galaxies would require either substantial modification of the standard Lambda cold dark matter paradigm or else significant revision to the uncertainties in their inferred mass profiles, which should be much larger than reported. Systematic errors in inclination may provide a simple resolution to what would otherwise be a rather intractable problem for the current paradigm.

We use mock interferometric HI measurements and a conventional tilted-ring modelling procedure to estimate circular velocity curves of dwarf galaxy discs from the APOSTLE suite of {\Lambda}CDM cosmological hydrodynamical simulations. The modelling yields a large diversity of rotation curves for an individual galaxy at fixed inclination, depending on the line-of-sight orientation. The diversity is driven by non-circular motions in the gas; in particular, by strong bisymmetric fluctuations in the azimuthal velocities that the tilted-ring model is ill-suited to account for and that are difficult to detect in model residuals. Large misestimates of the circular velocity arise when the kinematic major axis coincides with the extrema of the fluctuation pattern, in some cases mimicking the presence of kiloparsec-scale density 'cores', when none are actually present. The thickness of APOSTLE discs compounds this effect: more slowly-rotating extra-planar gas systematically reduces the average line-of-sight speeds. The recovered rotation curves thus tend to underestimate the true circular velocity of APOSTLE galaxies in the inner regions. Non-circular motions provide an appealing explanation for the large apparent cores observed in galaxies such as DDO 47 and DDO 87, where the model residuals suggest that such motions might have affected estimates of the inner circular velocities. Although residuals from tilted ring models in the simulations appear larger than in observed galaxies, our results suggest that non-circular motions should be carefully taken into account when considering the evidence for dark matter cores in individual galaxies.

Using the astrometry from the ESA's Gaia mission, previous works have shown that the Milky Way stellar halo is dominated by metal-rich stars on highly eccentric orbits. To shed light on the nature of this prominent halo component, we have analysed 28 Galaxy analogues in the Auriga suite of cosmological hydrodynamics zoom-in simulations. Some three quarters of the Auriga galaxies contain significant components with high radial velocity anisotropy, beta > 0.6. However, only in one third of the hosts do the high-beta stars contribute significantly to the accreted stellar halo overall, similar to what is observed in the Milky Way. For this particular subset we reveal the origin of the dominant stellar halo component with high metallicity, [Fe/H]~-1, and high orbital anisotropy, beta>0.8, by tracing their stars back to the epoch of accretion. It appears that, typically, these stars come from a single dwarf galaxy with a stellar mass of order of 10^9-10^10 Msol that merged around 6-10 Gyr ago, causing a sharp increase in the halo mass. Our study therefore establishes a firm link between the excess of radially anisotropic stellar debris in the Milky Way halo and an ancient head-on collision between the young Milky Way and a massive dwarf galaxy

Carbon enhanced metal poor (CEMP) stars with low abundances of neutron capture elements (CEMP-no stars) are ubiquitous among metal poor stars in the Milky Way. Recent observations have uncovered their two subgroups that differ in the carbon to magnesium ([C/Mg]) abundance ratio. Here we demonstrate that similar abundance patterns are also present in Milky Way-like galaxies in the EAGLE cosmological hydrodynamical simulation, where these patterns originate from the fact that stars may form from gas enriched predominantly by asymptotic giant branch (AGB) stars or by type-II supernovae. This occurs when stars form in the poorly mixed interstellar medium of Milky Way progenitor galaxies that are undergoing bursty star formation. The CEMP-no stars with lower [C/Mg] form at the onset of a starburst from gas enriched by low metallicity type-II supernovae that power a strong outflow, quenching further star formation. When star formation resumes following cosmological gas accretion, the CEMP-no stars with higher [C/Mg] form, with enrichment by AGB ejecta evident in their higher abundance of barium and lower abundance of magnesium. This suggests that bursty star formation in the progenitors of the Galaxy leaves a permanent imprint in the abundance patterns of CEMP stars.

Dark matter particles may decay, emitting photons. Drawing on the EAGLE family of hydrodynamic simulations of galaxy formation -- including the APOSTLE and C-EAGLE simulations -- we assess the systematic uncertainties and scatter on the decay flux from different galaxy classes, from Milky Way satellites to galaxy clusters, and compare our results to studies of the 3.55~keV line. We demonstrate that previous detections and non-detections of this line are consistent with a dark matter interpretation. For example, in our simulations the width of the the dark matter decay line for Perseus-analogue galaxy clusters lies in the range 1300-1700~\kms. Therefore, the non-detection of the 3.55~keV line in the centre of the Perseus cluster by the {\it Hitomi} collaboration is consistent with detections by other instruments. We also consider trends with stellar and halo mass and evaluate the scatter in the expected fluxes arising from the anisotropic halo mass distribution and from object-to-object variations. We provide specific predictions for observations with {\it XMM-Newton} and with the planned X-ray telescopes {\it XRISM} and {\it ATHENA}. If future detections of unexplained X-ray lines match our predictions, including line widths, we will have strong evidence that we have discovered the dark matter.

The Milky Way's (MW) satellite population is a powerful probe of warm dark matter (WDM) models as the abundance of small substructures is very sensitive to the properties of the WDM particle. However, only a partial census of the MW's complement of satellite galaxies exists because surveys of the MW's close environs are incomplete both in depth and in sky coverage. We present a new Bayesian analysis that combines the sample of satellites recently discovered by the Dark Energy Survey (DES) with those found in the Sloan Digital Sky Survey (SDSS) to estimate the total satellite galaxy luminosity function down to $M_{\rm V}=0$. We find that there should be at least $124^{+40}_{-27}$ ($68\%$ CL, statistical error) satellites as bright or brighter than $M_{\rm V}=0$ within 300 kpc of the Sun, with only a weak dependence on MW halo mass. When it comes online the Large Synoptic Survey Telescope should detect approximately half of this population. We also show that WDM models infer the same number of satellites as in $\Lambda$CDM, which will allow us to rule out those models that produce insufficient substructure to be viable.

We investigate the alignment of haloes with the filaments of the cosmic web using an unprecedently large sample of dark matter haloes taken from the P-Millennium $\Lambda$CDM cosmological N-body simulation. We use the state-of-the-art NEXUS morphological formalism which, due to its multiscale nature, simultaneously identifies structures at all scales. We find strong and highly significant alignments, with both the major axis of haloes and their peculiar velocity tending to orient along the filament. However, the spin - filament alignment displays a more complex trend changing from preferentially parallel at low masses to preferentially perpendicular at high masses. This "spin flip" occurs at an average mass of $5\times10^{11}~h^{-1}M_\odot$. This mass increases with increasing filament diameter, varying by more than an order of magnitude between the thinnest and thickest filament samples. We also find that the inner parts of haloes have a spin flip mass that is several times smaller than that of the halo as a whole. These results confirm that recent accretion is responsible for the complex behaviour of the halo spin - filament alignment. Low-mass haloes mainly accrete mass along directions perpendicular to their host filament and thus their spins tend to be oriented along the filaments. In contrast, high-mass haloes mainly accrete along their host filaments and have their spins preferentially perpendicular to them. Furthermore, haloes located in thinner filaments are more likely to accrete along their host filaments than haloes of the same mass located in thicker filaments.

We present and analyse mock stellar catalogues that match the selection criteria and observables (including uncertainties) of the Gaia satellite data release 2 (DR2). The source are six cosmological high-resolution magneto-hydrodynamic $\Lambda$CDM zoom simulations of the formation of Milky Way analogues from the AURIGA project. Mock data are provided for stars with $V < 16$ mag, and $V < 20$ mag at $|b|>20$ degrees. The mock catalogues are made using two different methods: the public SNAPDRAGONS code, and a method based on that of Lowing et al. that preserves the phase-space distribution of the model stars. These publicly available catalogues contain 5-parameter astrometry, radial velocities, multi-band photometry, stellar parameters, dust extinction values, and uncertainties in all these quantities. In addition, we provide the gravitational potential and information on the origin of each star. By way of demonstration, we apply the mock catalogues to analyses of the young stellar disc and the stellar halo. We show that: i) the young outer stellar disc exhibits a flared distribution that is detectable in the height and vertical velocity distribution of A- and B-dwarf stars up to radii of ~15 kpc; and ii) the spin of the stellar halo out to 100 kpc can be accurately measured with Gaia DR2 RR Lyrae stars. These catalogues are well suited for comparisons with observations and should help to: i) develop and test analysis methods for the Gaia DR2 data; ii) gauge the limitations and biases of the data and iii) interpret the data in the light of theoretical predictions from realistic $ab$ $initio$ simulations of galaxy formation in the $\Lambda$CDM cosmological model.

The (re)ionisation of hydrogen in the early universe has a profound effect on the formation of the first galaxies: by raising the gas temperature and pressure, it prevents gas from cooling into small haloes thus affecting the abundance of present-day small galaxies. Using the Galform semi-analytic model of galaxy formation, we show that two key aspects of the reionisation process -- when reionisation takes place and the characteristic scale below which it suppresses galaxy formation -- are imprinted in the luminosity function of dwarf galaxies. We focus on the luminosity function of satellites of galaxies like the Milky Way and the LMC, which is easier to measure than the luminosity function of the dwarf population as a whole. Our results show that the details of these two characteristic properties of reionisation determine the shape of the luminosity distribution of satellites in a unique way, and is largely independent of the other details of the galaxy formation model. Our models generically predict a bimodality in the distribution of satellites as a function of luminosity: a population of faint satellites and population of bright satellites separated by a 'valley' forged by reionisation. We show that this bimodal distribution is present at high statistical significance in the combined satellite luminosity function of the Milky Way and M31. We make predictions for the expected number of satellites around LMC-mass dwarfs where the bimodality may also be measurable in future observational programmes. Our preferred model predicts a total of $26 \pm 10$ (68 per cent confidence) satellites brighter than ${\rm M}_V=0$ in LMC-mass systems.

The total number and luminosity function of the population of dwarf galaxies of the Milky Way (MW) provide important constraints on the nature of the dark matter and on the astrophysics of galaxy formation at low masses. However, only a partial census of this population exists because of the flux limits and restricted sky coverage of existing Galactic surveys. We combine the sample of satellites recently discovered by the Dark Energy Survey (DES) with the satellites found in Sloan Digital Sky Survey (SDSS) Data Release 9 (together these surveys cover nearly half the sky) to estimate the total luminosity function of satellites down to $M_{\rm V}=0$. We apply a new Bayesian inference method in which we assume that the radial distribution of satellites independently of absolute magnitude follows that of subhaloes selected according to their peak maximum circular velocity. We find that there should be at least $124^{+40}_{-27}$ (68 per cent CL, statistical error) satellites brighter than $M_{\rm V}=0$ within $300$ kpc of the Sun. As a result of our use of new data and better simulations, and a more robust statistical method, we infer a much smaller population of satellites than reported in previous studies using earlier SDSS data only; we also address an underestimation of the uncertainties in earlier work by accounting for stochastic effects. We find that the inferred number of faint satellites depends only weakly on the assumed mass of the MW halo and we provide scaling relations to extend our results to different assumed halo masses and outer radii. We predict that half of our estimated total satellite population of the MW should be detected by the Large Synoptic Survey Telescope. The code implementing our estimation method is available online.

We investigate the population of dwarf galaxies with stellar masses similar to the Large Magellanic Cloud (LMC) and M33 in the EAGLE galaxy formation simulation. In the field, galaxies reside in haloes with stellar-to-halo mass ratios of $1.03^{+0.50}_{-0.31}\times10^{-2}$ (68% confidence level); systems like the LMC, which have an SMC-mass satellite, reside in haloes about 1.3 times more massive, which suggests an LMC halo mass at infall, $M_{200}=3.4^{+1.8}_{-1.2}\times10^{11}M_\odot$ (68% confidence level). The colour distribution of dwarfs is bimodal, with the red galaxies ($g-r>0.6$) being mostly satellites. The fraction of red LMC-mass dwarfs is 15% for centrals, and for satellites this fraction increases rapidly with host mass: from 10% for satellites of Milky Way (MW)-mass haloes to nearly 90% for satellites of groups and clusters. The quenching timescale, defined as the time after infall when half of the satellites have acquired red colours, decreases with host mass from ${>}5$ Gyrs for MW-mass hosts to $2.5$ Gyrs for cluster mass hosts. The satellites of MW-mass haloes have higher star formation rates and bluer colours than field galaxies. This is due to enhanced star formation triggered by gas compression shortly after accretion. Both the LMC and M33 have enhanced recent star formation that could be a manifestation of this process. After infall into their MW-mass hosts, the $g-r$ colours of LMC-mass dwarfs become bluer for the first 2 Gyrs, after which they rapidly redden. LMC-mass dwarfs fell into their MW-mass hosts only relatively recently, with more than half having an infall time of less than 3.5 Gyrs.

The Hubble sequence provides a useful classification of galaxy morphology at low redshift. However, morphologies are not static, but rather evolve as the growth of structure proceeds through mergers, accretion and secular processes. We investigate how kinematically defined disc and spheroidal structures form and evolve in the EAGLE hydrodynamic simulation of galaxy formation. At high redshift most galaxies of all masses are asymmetric. By redshift $z\simeq 1.5$ the Hubble sequence is established and after this time most of the stellar mass is in spheroids whose contribution to the stellar mass budget continues to rise to the present day. The stellar mass fraction in discs peaks at $z\simeq 0.5$ but overall remains subdominant at all times although discs contribute most of the stellar mass in systems of mass $M^*\sim 10^{10.5}{\rm M_\odot}$ at $z \leq 1.5$. Star formation occurs predominantly in disc structures throughout most of cosmic time but morphological transformations rearrange stars, thus establishing the low-redshift morphological mix. Morphological transformations are common and we quantify the rates at which they occur. The rate of growth of spheroids decreases at $z < 2$ while the rate of decay of discs remains roughly constant at $z < 1$. Finally, we find that the prograde component of galaxies becomes increasingly dynamically cold with time.

In a previous paper based on dark matter only simulations we show that, in the approximation of an analytic and static potential describing the strongly triaxial and cuspy shape of Milky Way-sized haloes, diffusion due to chaotic mixing in the neighbourhood of the Sun does not efficiently erase phase space signatures of past accretion events. In this second paper we further explore the effect of chaotic mixing using multicomponent Galactic potential models and solar neighbourhood-like volumes extracted from fully cosmological hydrodynamic simulations, thus naturally accounting for the gravitational potential associated with baryonic components, such as the bulge and disc. Despite the strong change in the global Galactic potentials with respect to those obtained in dark matter only simulations, our results confirm that a large fraction of halo particles evolving on chaotic orbits exhibit their chaotic behaviour after periods of time significantly larger than a Hubble time. In addition, significant diffusion in phase space is not observed on those particles that do exhibit chaotic behaviour within a Hubble time.

The massive stars that ionised the Universe have short lifetimes and can only be studied near the time of formation, but any low mass stars that formed contemporaneously might be observable in the local Universe today. We study the abundance pattern and spatial distribution of these `siblings of reionizers' (SoRs) in the EAGLE cosmological hydrodynamical simulation. SoRs tend to be enriched to super-solar levels in $\alpha$-elements compared to iron. In {\sc eagle} galaxies resembling the Milky Way, $\sim 40$ percent of carbon-enhanced metal poor (CEMP) stars are SoRs. Conversely, $\sim 10$ percent of all SoRs are CEMP stars. This fraction increases to $\gtrsim 50$ percent for SoRs of metallicity [Fe/H]$<-4$, and at such low metallicities, most of the CEMP stars are of CEMP-no subtype that are lacking neutron capture elements. Although these numbers may well depend on the details of the physical models implemented in EAGLE, the trends we describe are robust as they result from the strong feedback from star formation in early galaxies, itself a key ingredient of most current models of galaxy formation. We further find that most SoRs today reside in halos with mass $M_h\gtrapprox 10^{12}$ M$_\odot$, and 50 percent of them are in the halo of their central galaxy (distance $>10$ kpc), mainly because they were accreted onto their current host rather than formed in-situ. To a good approximation, the SoRs are CEMP-no stars that reside in the stellar halos of massive galaxies, with nearly half of them contributing to the intracluster light in groups and clusters.

Schwinn et al. (2017) have recently compared the abundance and distribution of massive substructures identified in a gravitational lensing analysis of Abell 2744 by Jauzac et al. (2016) and N-body simulation and found no cluster in {\Lambda}CDM simulation that is similar to Abell 2744. Schwinn et al.(2017) identified the measured projected aperture masses with the actual masses associated with subhaloes in the MXXL N-body simulation. We have used the high resolution Phoenix cluster simulations to show that such an identification is incorrect: the aperture mass is dominated by mass in the body of the cluster that happens to be projected along the line-of-sight to the subhalo. This enhancement varies from factors of a few to factors of more than 100, particularly for subhaloes projected near the centre of the cluster. We calculate aperture masses for subhaloes in our simulation and compare them to the measurements for Abell 2744. We find that the data for Abell 2744 are in excellent agreement with the matched predictions from {\Lambda}CDM. We provide further predictions for aperture mass functions of subhaloes in idealized surveys with varying mass detection thresholds.

We investigate how the nature of the galaxies that reionized the Universe affects the duration of reionization. We contrast two sets of models: one in which galaxies on the faint side of the luminosity function dominate the ionizing emissivity, and a second in which the galaxies on the bright side of the luminosity function dominate. The faint-end of the luminosity function evolves slowly, therefore the transition from mostly neutral to mostly ionized state takes a much longer time in the first set of models compared to the second. Existing observational constraints on the duration of this transition are relatively weak, but taken at face value prefer the model in which galaxies on the bright side play a major role. Measurements of the kinetic Sunyaev Zeldovich effect in the cosmic microwave background from the epoch of reionization also point in the same direction.

We present predictions for the quenching of star formation in satellite galaxies of the Local Group from a suite of 30 cosmological zoom simulations of Milky Way-like host galaxies. The Auriga simulations resolve satellites down to the luminosity of the classical dwarf spheroidal galaxies of the Milky Way. We find strong mass-dependent and distance-dependent quenching signals, where dwarf systems beyond 600 kpc are only strongly quenched below a stellar mass of $10^7$ M$_\odot$. Ram pressure stripping appears to be the dominant quenching mechanism and 50% of quenched systems cease star formation within 1 Gyr of first infall. We demonstrate that systems within a host galaxy's $R_{200}$ radius are comprised of two populations: (i) a first infall population that has entered the host halo within the past few Gyrs and (ii) a population of returning `backsplash' systems that have had a much more extended interaction with the host. Backsplash galaxies that do not return to the host galaxy by redshift zero exhibit quenching properties similar to galaxies within $R_{200}$ and are distinct from other external systems. The simulated quenching trend with stellar mass has some tension with observations, but our simulations are able reproduce the range of quenching times measured from resolved stellar populations of Local Group dwarf galaxies.

The shallow faint-end slope of the galaxy mass function is usually reproduced in $\Lambda$CDM galaxy formation models by assuming that the fraction of baryons that turns into stars drops steeply with decreasing halo mass and essentially vanishes in haloes with maximum circular velocities $V_{\rm max}<20$-$30$ km/s. Dark matter-dominated dwarfs should therefore have characteristic velocities of about that value, unless they are small enough to probe only the rising part of the halo circular velocity curve (i.e., half-mass radii, $r_{1/2}\ll 1$ kpc). Many dwarfs have properties in disagreement with this prediction: they are large enough to probe their halo $V_{\rm max}$ but their characteristic velocities are well below $20$ km/s. These `cold faint giants' (an extreme example is the recently discovered Crater 2 Milky Way satellite) can only be reconciled with our $\Lambda$CDM models if they are the remnants of once massive objects heavily affected by tidal stripping. We examine this possibility using the APOSTLE cosmological hydrodynamical simulations of the Local Group. Assuming that low velocity dispersion satellites have been affected by stripping, we infer their progenitor masses, radii, and velocity dispersions, and find them in remarkable agreement with those of isolated dwarfs. Tidal stripping also explains the large scatter in the mass discrepancy-acceleration relation in the dwarf galaxy regime: tides remove preferentially dark matter from satellite galaxies, lowering their accelerations below the $a_{\rm min}\sim 10^{-11} m/s^2$ minimum expected for isolated dwarfs. In many cases, the resulting velocity dispersions are inconsistent with the predictions from Modified Newtonian Dynamics, a result that poses a possibly insurmountable challenge to that scenario.

The spherical Jeans equation (SJE) is widely used in dynamical modelling of the Milky Way (MW) halo potential. We use haloes and galaxies from the cosmological Millennium-II simulation and hydrodynamical APOSTLE simulations to investigate the performance of the SJE in recovering the underlying mass profiles of MW mass haloes. The best-fitting halo mass and concentration parameters scatter by 25% and 40% around their input values, respectively, when dark matter particles are used as tracers. This scatter becomes as large as a factor of 3 when using star particles instead. This is significantly larger than the estimated statistical uncertainty associated with the use of the SJE. The existence of correlated phase-space structures that violate the steady state assumption of the SJE as well as non-spherical geometries are the principal sources of the scatter. Binary haloes show larger scatter because they are more aspherical in shape and have a more perturbed dynamical state. Our results confirm the previous study of Wang et al. (2017) that the number of independent phase-space structures sets an intrinsic limiting precision on dynamical inferences based on the steady state assumption. Modelling with a radius-independent velocity anisotropy, or using tracers within a limited outer radius, result in significantly larger scatter, but the ensemble-averaged measurement over the whole halo sample is approximately unbiased.

We study the incidence of group and filamentary dwarf galaxy accretion into Milky Way (MW) mass haloes using two types of hydrodynamical simulations: EAGLE, which resolves a large cosmological volume, and the AURIGA suite, which are very high resolution zoom-in simulations of individual MW-sized haloes. The present-day 11 most massive satellites are predominantly (75%) accreted in single events, 14% in pairs and 6% in triplets, with higher group multiplicities being unlikely. Group accretion becomes more common for fainter satellites, with 60% of the top 50 satellites accreted singly, 12% in pairs, and 28% in richer groups. A group similar in stellar mass to the Large Magellanic Cloud (LMC) would bring on average 15 members with stellar mass larger than $10^4{~\rm M_\odot}$. Half of the top 11 satellites are accreted along the two richest filaments. The accretion of dwarf galaxies is highly anisotropic, taking place preferentially perpendicular to the halo minor axis, and, within this plane, preferentially along the halo major axis. The satellite entry points tend to be aligned with the present-day central galaxy disc and satellite plane, but to a lesser extent than with the halo shape. Dwarfs accreted in groups or along the richest filament have entry points that show an even larger degree of alignment with the host halo than the full satellite population. We also find that having most satellites accreted as a single group or along a single filament is unlikely to explain the MW disc of satellites.

We compare the transverse velocity dispersions recently measured within the Sculptor dwarf spheroidal galaxy to the predictions of previously published dynamical models. These provide good fits to the observed number count and velocity dispersion profiles of metal-rich and metal-poor stars both in cored and in cusped potentials. At the projected radius where the proper motions were measured, these models predict transverse dispersions in the range 6 to 9.5 km/s, with the tangential dispersion about 1 km/s larger than the (projected) radial dispersion. Both dispersions are predicted to be about 1 km/s larger for metal-poor than for metal-rich stars. At this projected radius, cored and cusped potentials predict almost identical transverse dispersions. The measured tangential dispersion ($8.5 \pm 3.2$ km/s) agrees remarkably well with these predictions, while the measured radial dispersion ($11.5 \pm 4.3$ km/s) differs only at about the $1\sigma$ level. Thus, the proper motion data are in excellent agreement with previous data but do not help to distinguish between cored and cusped potentials. This will require velocity dispersion data (either from proper motions or from radial velocities) with uncertainties well below 1 km/s over a range of projected radii.

We examine the origin of radial and vertical gradients in the age/metallicity of the stellar component of a galaxy disc formed in the APOSTLE cosmological hydrody- namical simulations. Some of these gradients resemble those in the Milky Way, where they have sometimes been interpreted as due to internal evolution, such as scattering off giant molecular clouds, radial migration driven by spiral patterns, or orbital reso- nances with a bar. Secular processes play a minor role in the simulated galaxy, which lacks strong spiral or bar patterns, and where such gradients arise as a result of the gradual enrichment of a gaseous disc that is born thick but thins as it turns into stars and settles into centrifugal equilibrium. The settling is controlled by the feedback of young stars; which links the star formation, enrichment, and equilibration timescales, inducing radial and vertical gradients in the gaseous disc and its descendent stars. The kinematics of coeval stars evolve little after birth and provide a faithful snapshot of the gaseous disc structure at the time of their formation. In this interpretation, the age-velocity dispersion relation would reflect the gradual thinning of the disc rather than the importance of secular orbit scattering; the outward flaring of stars would result from the gas disc flare rather than from radial migration; and vertical gradients would arise because the gas disc gradually thinned as it enriched. Such radial and vertical trends might just reflect the evolving properties of the parent gaseous disc, and are not necessarily the result of secular evolutionary processes.

We investigate the nature of carbon-enhanced metal poor (CEMP) stars in Milky Way (MW) analogues selected from the EAGLE cosmological hydrodynamical simulation. The stellar evolution model in EAGLE includes the physics of enrichment by asymptotic giant branch (AGB) stars, winds from massive stars, and type I and type II supernovae (SNe). In the simulation, star formation in young MW progenitors is bursty due to efficient stellar feedback, which causes poor metal mixing leading to the formation of CEMP stars with extreme abundance patterns. In this scenario, two classes of CEMP stars emerge: those mostly enriched by low-metallicity type II SNe with low Fe yields that drive galactic outflows, and those mostly enriched by AGB stars when a gas-poor progenitor accretes pristine gas. The first class resembles CEMP-no stars with high [C/Fe] and low [C/O], the second class resembles CEMP-s stars overabundant in s-process elements and high values of [C/O]. This scenario explains several trends seen in data: (i) the increase in the scatter and median of [C/O] at low and decreasing [O/H], (ii) the trend of stars with very low [Fe/H] or [C/H] to be of type CEMP-no, and (iii) the reduction in the scatter of [{\alpha}/Fe] with atomic number in metal poor stars. In this interpretation CEMP-no stars were enriched by the stars that enabled galaxies to reionise the Universe.

The existence of two kinematically and chemically distinct stellar subpopulations in the Sculptor and Fornax dwarf galaxies offers the opportunity to constrain the density profile of their matter haloes by measuring the mass contained within the well-separated half-light radii of the two metallicity subpopulations. Walker and Penarrubia have used this approach to argue that data for these galaxies are consistent with constant-density `cores' in their inner regions and rule out `cuspy' Navarro-Frenk-White (NFW) profiles with high statistical significance, particularly in the case of Sculptor. We test the validity of these claims using dwarf galaxies in the APOSTLE (A Project Of Simulating The Local Environment) Lambda cold dark matter cosmological hydrodynamics simulations of analogues of the Local Group. These galaxies all have NFW dark matter density profiles and a subset of them develop two distinct metallicity subpopulations reminiscent of Sculptor and Fornax. We apply a method analogous to that of Walker and Penarrubia to a sample of 53 simulated dwarfs and find that this procedure often leads to a statistically significant detection of a core in the profile when in reality there is a cusp. Although multiple factors contribute to these failures, the main cause is a violation of the assumption of spherical symmetry upon which the mass estimators are based. The stellar populations of the simulated dwarfs tend to be significantly elongated and, in several cases, the two metallicity populations have different asphericity and are misaligned. As a result, a wide range of slopes of the density profile are inferred depending on the angle from which the galaxy is viewed.

Semi-analytic galaxy formation models are widely used to gain insight into the astrophysics of galaxy formation and in model testing, parameter space searching and mock catalogue building. In this work we present a new model for gas cooling in halos in semi-analytic models, which improves over previous cooling models in several ways. Our new treatment explicitly includes the evolution of the density profile of the hot gas driven by the growth of the dark matter halo and by the dynamical adjustment of the gaseous corona as gas cools down. The effect of the past cooling history on the current mass cooling rate is calculated more accurately, by doing an integral over the past history. The evolution of the hot gas angular momentum profile is explicitly followed, leading to a self-consistent and more detailed calculation of the angular momentum of the cooled down gas. This model predicts higher cooled down masses than the cooling models previously used in GALFORM, closer to the predictions of the cooling models in L-GALAXIES and MORGANA, even though those models are formulated differently. It also predicts cooled down angular momenta that are higher than in previous GALFORM cooling models, but generally lower than the predictions of L-GALAXIES and MORGANA. When used in a full galaxy formation model, this cooling model improves the predictions for early-type galaxy sizes in GALFORM.

We study the shape and kinematics of simulated dwarf galaxy discs in the APOSTLE suite of $\Lambda$CDM cosmological hydrodynamical simulations. We find that a large fraction of these gas-rich, star-forming discs show weak bars in their stellar component, despite being dark matter-dominated systems. The bar pattern shape and orientation reflect the ellipticity of the dark matter potential, and its rotation is locked to the slow figure rotation of the triaxial dark halo. The bar-like nature of the potential induces non-circular motions in the gas component, including strong bisymmetric flows that can be readily seen as m=3 harmonic perturbations in the HI line-of-sight velocity fields. Similar bisymmetric flows are seen in many galaxies of the THINGS and LITTLE THINGS surveys, although on average their amplitudes are a factor of ~2 weaker than in our simulated discs. Our results indicate that bar-like patterns may arise even when baryons are not dominant, and that they are common enough to warrant careful consideration when analyzing the gas kinematics of dwarf galaxy discs.

We present predictions for the outcome of deep galaxy surveys with the $James$ $Webb$ $Space$ $Telescope$ ($JWST$) obtained from a physical model of galaxy formation in $\Lambda$CDM. We use the latest version of the GALFORM model, embedded within a new ($800$ Mpc)$^{3}$ dark matter only simulation with a halo mass resolution of $M_{\rm halo}>2\times10^{9}$ $h^{-1}$ M$_{\odot}$. For computing full UV-to-mm galaxy spectral energy distributions, including the absorption and emission of radiation by dust, we use the spectrophotometric radiative transfer code GRASIL. The model is calibrated to reproduce a broad range of observational data at $z\lesssim6$, and we show here that it can also predict evolution of the rest-frame far-UV luminosity function for $7\lesssim z\lesssim10$ which is in good agreement with observations. We make predictions for the evolution of the luminosity function from $z=16$ to $z=0$ in all broadband filters on the Near InfraRed Camera (NIRCam) and Mid InfraRed Instrument (MIRI) on $JWST$ and present the resulting galaxy number counts and redshift distributions. Our fiducial model predicts that $\sim1$ galaxy per field of view will be observable at $z\sim11$ for a $10^4$ s exposure with NIRCam. A variant model, which produces a higher redshift of reionization in better agreement with $Planck$ data, predicts number densities of observable galaxies $\sim5\times$ greater at this redshift. Similar observations with MIRI are predicted not to detect any galaxies at $z\gtrsim6$. We also make predictions for the effect of different exposure times on the redshift distributions of galaxies observable with $JWST$, and for the angular sizes of galaxies in $JWST$ bands.

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

N-body dark matter simulations of structure formation in the $\Lambda$CDM model predict a population of subhalos within Galactic halos that have higher central densities than inferred for satellites of the Milky Way, a tension known as the `too big to fail' problem. Proposed solutions include baryonic effects, a smaller mass for the Milky Way halo, and warm dark matter. We test these three possibilities using a semi-analytic model of galaxy formation to generate luminosity functions for Milky Way halo-analogue satellite populations, the results of which are then coupled to the Jiang & van den Bosch model of subhalo stripping to predict the subhalo $V_\mathrm{max}$ functions for the 10 brightest satellites. We find that selecting the brightest satellites (as opposed to the most massive) and modelling the expulsion of gas by supernovae at early times increases the likelihood of generating the observed Milky Way satellite $V_\mathrm{max}$ function. The preferred halo mass is $6\times10^{11}M_{\odot}$, which has a 14 percent probability to host a $V_\mathrm{max}$ function like that of the Milky Way satellites. This probability is reduced to 8 percent for a $1.0\times10^{12}M_{\odot}$ halo and to 3 percent for a $1.4\times10^{12}M_{\odot}$ halo. We conclude that the Milky Way satellite $V_\mathrm{max}$ function is compatible with a CDM cosmology, as previously found by Sawala et al. using hydrodynamic simulations. Sterile neutrino-warm dark matter models achieve a higher degree of agreement with the observations, with a maximum 35 percent chance of generating the observed Milky Way satellite $V_\mathrm{max}$ function. However, more work is required to check that the semi-analytic stripping model is calibrated correctly in the sterile neutrino cosmology, and to check if our sterile neutrino models produce sufficient numbers of faint satellites.

We study galaxy formation in sterile neutrino dark matter models that differ signifi- cantly from both cold and from `warm thermal relic' models. We use the EAGLE code to carry out hydrodynamic simulations of the evolution of pairs of galaxies chosen to resemble the Local Group, as part of the APOSTLE simulations project. We compare cold dark matter (CDM) with two sterile neutrino models with 7 keV mass: one, the warmest among all models of this mass (LA120) and the other, a relatively cold case (LA10). We show that the lower concentration of sterile neutrino subhalos compared to their CDM counterparts makes the inferred inner dark matter content of galaxies like Fornax (or Magellanic Clouds) less of an outlier in the sterile neutrino cosmologies. In terms of the galaxy number counts the LA10 simulations are indistinguishable from CDM when one takes into account halo-to-halo (or `simulation-to-simulation') scatter. In order for the LA120 model to match the number of Local Group dwarf galaxies, a higher fraction of low mass haloes is required to form galaxies than is predicted by the EAGLE simulations. As the census of the Local Group galaxies nears completion, this population may provide a strong discriminant between cold and warm dark matter models.

Particle tagging is an efficient, but approximate, technique for using cosmological N-body simulations to model the phase-space evolution of the stellar populations predicted, for example, by a semi-analytic model of galaxy formation. We test the technique developed by Cooper et al. (which we call STINGS here) by comparing particle tags with stars in a smooth particle hydrodynamic (SPH) simulation. We focus on the spherically averaged density profile of stars accreted from satellite galaxies in a Milky Way (MW)-like system. The stellar profile in the SPH simulation can be recovered accurately by tagging dark matter (DM) particles in the same simulation according to a prescription based on the rank order of particle binding energy. Applying the same prescription to an N-body version of this simulation produces a density profile differing from that of the SPH simulation by <10 per cent on average between 1 and 200 kpc. This confirms that particle tagging can provide a faithful and robust approximation to a self-consistent hydrodynamical simulation in this regime (in contradiction to previous claims in the literature). We find only one systematic effect, likely due to the collisionless approximation, namely that massive satellites in the SPH simulation are disrupted somewhat earlier than their collisionless counterparts. In most cases this makes remarkably little difference to the spherically averaged distribution of their stellar debris. We conclude that, for galaxy formation models that do not predict strong baryonic effects on the present-day DM distribution of MW-like galaxies or their satellites, differences in stellar halo predictions associated with the treatment of star formation and feedback are much more important than those associated with the dynamical limitations of collisionless particle tagging.

Dark matter subhalos are the remnants of (incomplete) halo mergers. Identifying them and establishing their evolutionary links in the form of merger trees is one of the most important applications of cosmological simulations. The Hierachical Bound-Tracing (HBT) code identifies halos as they form and tracks their evolution as they merge, simultaneously detecting subhalos and building their merger trees. Here we present a new implementation of this approach, HBT+, that is much faster, more user friendly, and more physically complete than the original code. Applying HBT+ to cosmological simulations we show that both the subhalo mass function and the peak-mass function are well fit by similar double-Schechter functions.The ratio between the two is highest at the high mass end, reflecting the resilience of massive subhalos that experience substantial dynamical friction but limited tidal stripping. The radial distribution of the most massive subhalos is more concentrated than the universal radial distribution of lower mass subhalos. Subhalo finders that work in configuration space tend to underestimate the masses of massive subhalos, an effect that is stronger in the host centre. This may explain, at least in part, the excess of massive subhalos in galaxy cluster centres inferred from recent lensing observations. We demonstrate that the peak-mass function is a powerful diagnostic of merger tree defects, and the merger trees constructed using HBT+ do not suffer from the missing or switched links that tend to afflict merger trees constructed from more conventional halo finders. We make the HBT+ code publicly available.

Many dynamical models of the Milky Way halo require assumptions that the distribution function of a tracer population should be independent of time (i.e., a steady state distribution function) and that the underlying potential is spherical. We study the limitations of such modelling by applying a general dynamical model with minimal assumptions to a large sample of galactic haloes from cosmological $N$-body and hydrodynamical simulations. Using dark matter particles as dynamical tracers, we find that the systematic uncertainties in the measured mass and concentration parameters typically have an amplitude of 25% to 40%. When stars are used as tracers, however, the systematic uncertainties can be as large as a factor of $2-3$. The systematic uncertainties are not reduced by increasing the tracer sample size and vary stochastically from halo to halo. These systematic uncertainties are mostly driven by underestimated statistical noise caused by correlated phase-space structures that violate the steady state assumption. The number of independent phase-space structures inferred from the uncertainty level sets a limiting sample size beyond which a further increase no longer significantly improves the accuracy of dynamical inferences. The systematic uncertainty level is determined by the halo merger history, the shape and environment of the halo. Our conclusions apply generally to any spherical steady-state model.

We study the vertical structure of polytropic, $P\propto \rho^\Gamma$, centrifugally-supported gaseous discs embedded in cold dark matter (CDM) halos. At fixed radius $R$, the shape of the vertical density profile depends only weakly on whether the disc is self-gravitating (SG) or not (NSG). The disc thickness, set by the midplane sound speed and circular velocity, $(c_s/V_c)R$, in the NSG case, and by the sound speed and surface density, $c_s^2/G\Sigma$, in SG discs, is smaller than either of these scales. SG discs are typically Toomre unstable, NSG discs are stable. Exponential discs in CDM halos with roughly flat circular velocity curves generally "flare" outwards. For the polytropic equation of state of the EAGLE simulations, discs whose mass and size match observational constraints are stable (NSG) for $M_d< 3\times 10^9\, M_\odot$ and unstable (SG) at higher masses, if fully gaseous. We test these analytic results using a set of idealized SPH simulations and find excellent agreement. Our results clarify the role of the gravitational softening on the thickness of simulated discs, and on the onset of radial instabilities. EAGLE low-mass discs are non-self-gravitating so the softening plays no role in their vertical structure. High-mass discs, on the other hand, are expected to be self-gravitating and unstable, and may be artificially thickened and stabilized unless gravity is well resolved. Simulations with spatial resolution high enough to not compromise the vertical structure of a disc also resolve the onset of their instabilities, but the converse is not true: resolving instabilities does not guarantee that the vertical structure is resolved.

Small distortions in the images of Einstein rings or giant arcs offer the exciting prospect of detecting dark matter haloes or subhaloes of mass below $10^9$M$_{\odot}$, most of which are too small to have made a visible galaxy. A very large number of such haloes are predicted to exist in the cold dark matter model of cosmogony; in contrast other models, such as warm dark matter, predict no haloes below a mass of this order which depends on the properties of the warm dark matter particle. Attempting to detect these small perturbers could therefore discriminate between different kinds of dark matter particles, and even rule out the cold dark matter model altogether. Globular clusters in the lens galaxy also induce distortions in the image which could, in principle, contaminate the test. Here, we investigate the population of globular clusters in six early type galaxies in the Virgo cluster. We find that the number density of globular clusters of mass $\sim10^6$M$_{\odot}$ is comparable to that of the dark matter perturbers (including subhaloes in the lens and haloes along the line-of-sight). We show that the very different degrees of mass concentration in globular clusters and dark matter haloes result in different lensing distortions. These are detectable with milli-arcsecond resolution imaging which can distinguish between globular cluster and dark matter halo signals.

We examine the origin of the mass discrepancy--radial acceleration relation (MDAR) of disk galaxies. This is a tight empirical correlation between the disk centripetal acceleration and that expected from the baryonic component. The MDAR holds for most radii probed by disk kinematic tracers, regardless of galaxy mass or surface brightness. The relation has two characteristic accelerations; $a_0$, above which all galaxies are baryon-dominated; and $a_{\rm min}$, an effective minimum aceleration probed by kinematic tracers in isolated galaxies. We use a simple model to show that these trends arise naturally in $\Lambda$CDM. This is because: (i) disk galaxies in $\Lambda$CDM form at the centre of dark matter haloes spanning a relatively narrow range of virial mass; (ii) cold dark matter halo acceleration profiles are self-similar and have a broad maximum at the centre, reaching values bracketed precisely by $a_{\rm min}$ and $a_0$ in that mass range; and (iii) halo mass and galaxy size scale relatively tightly with the baryonic mass of a galaxy in any successful $\Lambda$CDM galaxy formation model. Explaining the MDAR in $\Lambda$CDM does not require modifications to the cuspy inner mass profiles of dark haloes, although these may help to understand the detailed rotation curves of some dwarf galaxies and the origin of extreme outliers from the main relation. The MDAR is just a reflection of the self-similar nature of cold dark matter haloes and of the physical scales introduced by the galaxy formation process.

The standard galaxy formation theory assumes that baryons and dark matter are initially well-mixed before becoming segregated due to radiative cooling. We use non-radiative hydrodynamical simulations to explicitly examine this assumption and find that baryons and dark matter can also be segregated because of different physics obeyed by gas and dark matter during the build-up of the halo. As a result, baryons in many haloes do not originate from the same Lagrangian region as the dark matter. When using the fraction of corresponding dark matter and gas particles in the initial conditions (the "paired fraction") as a proxy of the dark matter and gas segregation strength of a halo, on average about $25$ percent of the baryonic and dark matter of the final halo are segregated in the initial conditions. This is at odds with the assumption of the standard galaxy formation model. A consequence of this effect is that the baryons and dark matter of the same halo initially experience different tidal torques and thus their angular momentum vectors are often misaligned. The degree of the misalignment is largely preserved during later halo assembly and can be understood with the tidal torque theory. The result challenges the precision of some semi-analytical approaches which utilize dark matter halo merger trees to infer properties of gas associated to dark matter haloes.

The Gaia astrometric mission may offer an unprecedented opportunity to discover new tidal streams in the Galactic halo. To test this, we apply nGC3, a great-circle-cell count method that combines position and proper motion data to identify streams, to ten mock Gaia catalogues of K giants and RR Lyrae stars constructed from cosmological simulations of Milky Way analogues. We analyse two sets of simulations, one using a combination of $N$-body and semi-analytical methods which has extremely high resolution, the other using hydro-dynamical methods, which captures the dynamics of baryons, including the formation of an in situ halo. These ten realisations of plausible Galactic merger histories allow us to assess the potential for the recovery of tidal streams in different Milky Way formation scenarios. We include the Gaia~selection function and observational errors in these mock catalogues. We find that the nGC3 method has a well-defined detection boundary in the space of stream width and projected overdensity, that can be predicted based on direct observables alone. We predict that about 4-13 dwarf galaxy streams can be detected in a typical Milky Way-mass halo with Gaia+nGC3, with an estimated efficiency of $>$80\% inside the detection boundary. The progenitors of these streams are in the mass range of the classical dwarf galaxies and may have been accreted as early as redshift $\sim3$. Finally, we analyse how different possible extensions of the Gaia mission will improve the detection of tidal streams.

We present mock optical images, broad-band and H$\alpha$ fluxes, and D4000 spectral indices for 30,145 galaxies from the EAGLE hydrodynamical simulation at redshift $z=0.1$, modelling dust with the SKIRT Monte Carlo radiative transfer code. The modelling includes a subgrid prescription for dusty star-forming regions, with both the subgrid obscuration of these regions and the fraction of metals in diffuse interstellar dust calibrated against far-infrared fluxes of local galaxies. The predicted optical colours as a function of stellar mass agree well with observation, with the SKIRT model showing marked improvement over a simple dust screen model. The orientation dependence of attenuation is weaker than observed because EAGLE galaxies are generally puffier than real galaxies, due to the pressure floor imposed on the interstellar medium. The mock H$\alpha$ luminosity function agrees reasonably well with the data, and we quantify the extent to which dust obscuration affects observed H$\alpha$ fluxes. The distribution of D4000 break values is bimodal, as observed. In the simulation, 20$\%$ of galaxies deemed `passive' for the SKIRT model, i.e. exhibiting D4000 $> 1.8$, are classified `active' when ISM dust attenuation is not included. The fraction of galaxies with stellar mass greater than $10^{10}$ M$_\odot$ that are deemed passive is slightly smaller than observed, which is due to low levels of residual star formation in these simulated galaxies. Colour images, fluxes and spectra of EAGLE galaxies are to be made available through the public EAGLE database.

The magnetic fields observed in the Milky~Way and nearby galaxies appear to be in equipartition with the turbulent, thermal, and cosmic ray energy densities, and hence are expected to be dynamically important. However, the origin of these strong magnetic fields is still unclear, and most previous attempts to simulate galaxy formation from cosmological initial conditions have ignored them altogether. Here, we analyse the magnetic fields predicted by the simulations of the Auriga Project, a set of 30 high-resolution cosmological zoom simulations of Milky~Way-like galaxies, carried out with a moving-mesh magneto-hydrodynamics code and a detailed galaxy formation physics model. We find that the magnetic fields grow exponentially at early times owing to a small-scale dynamo with an e-folding time of roughly $100\,\rm{Myr}$ in the center of halos until saturation occurs around $z=2-3$, when the magnetic energy density reaches about $10\%$ of the turbulent energy density with a typical strength of $10-50\,\rm{\mu G}$. In the galactic centers the ratio between magnetic and turbulent energy remains nearly constant until $z=0$. At larger radii, differential rotation in the disks leads to linear amplification that typically saturates around $z=0.5$ to $z=0$. The final radial and vertical variations of the magnetic field strength can be well described by two joint exponential profiles, and are in good agreement with observational constraints. Overall, the magnetic fields have only little effect on the global evolution of the galaxies as it takes too long to reach equipartition. We also demonstrate that our results are well converged with numerical resolution.

We characterize the contribution from accreted material to the galactic discs of the Auriga Project, a set of high resolution magnetohydrodynamic cosmological simulations of late-type galaxies performed with the moving-mesh code AREPO. Our goal is to explore whether a significant accreted (or ex-situ) stellar component in the Milky Way disc could be hidden within the near-circular orbit population, which is strongly dominated by stars born in-situ. One third of our models shows a significant ex-situ disc but this fraction would be larger if constraints on orbital circularity were relaxed. Most of the ex-situ material ($\gtrsim 50\%$) comes from single massive satellites ($> 6 \times 10^{10}~M_{\odot}$). These satellites are accreted with a wide range of infall times and inclination angles (up to $85^{\circ}$). Ex-situ discs are thicker, older and more metal-poor than their in-situ counterparts. They show a flat median age profile, which differs from the negative gradient observed in the in-situ component. As a result, the likelihood of identifying an ex-situ disc in samples of old stars on near-circular orbits increases towards the outskirts of the disc. We show three examples that, in addition to ex-situ discs, have a strongly rotating dark matter component. Interestingly, two of these ex-situ stellar discs show an orbital circularity distribution that is consistent with that of the in-situ disc. Thus, they would not be detected in typical kinematic studies.

The observed stellar kinematics of dispersion-supported galaxies are often used to measure dynamical masses. Recently, several analytical relationships between the stellar line-of-sight velocity dispersion, the projected (2D) or deprojected (3D) half-light radius, and the total mass enclosed within the half-light radius, relying on the spherical Jeans equation, have been proposed. Here, we make use of the APOSTLE cosmological hydrodynamical simulations of the Local Group to test the validity and accuracy of such mass estimators for both dispersion and rotation-supported galaxies, for field and satellite galaxies, and for galaxies of varying masses, shapes, and velocity dispersion anisotropies. We find that the mass estimators of Walker et al. and Wolf et al. are able to recover the masses of dispersion-dominated systems with little systematic bias, but with a 1-sigma scatter of 25 and 23 percent, respectively. The error on the estimated mass is dominated by the impact of the 3D shape of the stellar mass distribution, which is difficult to constrain observationally. This intrinsic scatter becomes the dominant source of uncertainty in the masses estimated for galaxies like the dwarf spheroidal (dSph) satellites of the Milky Way, where the observational errors in their sizes and velocity dispersions are small. Such scatter may also affect the inner density slopes of dSphs derived from multiple stellar populations, relaxing the significance with which Navarro-Frenk-White profiles may be excluded, depending on the degree to which the relevant properties of the different stellar populations are correlated. Finally, we derive a new optimal mass estimator that removes the residual biases and achieves a statistically significant reduction in the scatter to 20 percent overall for dispersion-dominated galaxies, allowing more precise and accurate mass estimates.

We analyze the total and baryonic acceleration profiles of a set of well-resolved galaxies identified in the EAGLE suite of hydrodynamic simulations. Our runs start from the same initial conditions but adopt different prescriptions for unresolved stellar and AGN feedback, resulting in diverse populations of galaxies by the present day. Some of them reproduce observed galaxy scaling relations, while others do not. However, regardless of the feedback implementation, all of our galaxies follow closely a simple relationship between the total and baryonic acceleration profiles, consistent with recent observations of rotationally supported galaxies. The relation has small scatter: different feedback implementations -- which produce different galaxy populations -- mainly shift galaxies along the relation, rather than perpendicular to it. Furthermore, galaxies exhibit a characteristic acceleration, $g_{\dagger}$, above which baryons dominate the mass budget, as observed. These observations, consistent with simple modified Newtonian dynamics, can be accommodated within the standard cold dark matter paradigm.

We use a pair of high resolution N-body simulations implementing two dark matter models, namely the standard cold dark matter (CDM) cosmogony and a warm dark matter (WDM) alternative where the dark matter particle is a 1.5keV thermal relic. We combine these simulations with the GALFORM semi-analytical galaxy formation model in order to explore differences between the resulting galaxy populations. We use GALFORM model variants for CDM and WDM that result in the same z=0 galaxy stellar mass function by construction. We find that most of the studied galaxy properties have the same values in these two models, indicating that both dark matter scenarios match current observational data equally well. Even in under-dense regions, where discrepancies in structure formation between CDM and WDM are expected to be most pronounced, the galaxy properties are only slightly different. The only significant difference in the local universe we find is in the galaxy populations of "Local Volumes", regions of radius 1 to 8Mpc around simulated Milky Way analogues. In such regions our WDM model provides a better match to observed local galaxy number counts and is five times more likely than the CDM model to predict sub-regions within them that are as empty as the observed Local Void. Thus, a highly complete census of the Local Volume and future surveys of void regions could provide constraints on the nature of dark matter.

We investigate the properties of the galaxies that reionized the Universe and the history of cosmic reionization using the "Evolution and Assembly of GaLaxies and their environments" (EAGLE) cosmological hydrodynamical simulations. We obtain the evolution of the escape fraction of ionizing photons in galaxies assuming that galactic winds create channels through which 20~percent of photons escape when the local surface density of star formation is greater than $0.1$ M$_\odot$ yr$^{-1}$ kpc$^{-2}$. Such threshold behaviour for the generation of winds is observed, and the rare local objects which have such high star formation surface densities exhibit high escape fractions of $\sim 10$ percent. In our model the luminosity-weighted mean escape fraction increases with redshift as $\bar f_{\rm esc}=0.045~((1+z)/4)^{1.1}$ at $z>3$, and the galaxy number weighted mean as $\langle f_{\rm esc} \rangle=2.2\times10^{-3}~((1+z)/4)^4$, and becomes constant $\approx0.2$ at redshift $z>10$. The escape fraction evolves as an increasingly large fraction of stars forms above the critical surface density of star formation at earlier times. This evolution of the escape fraction, combined with that of the star formation rate density from EAGLE, reproduces the inferred evolution of the filling factor of ionized regions during the reionization epoch ($6<z<8$), the evolution of the post-reionization ($0\leq z<6$) hydrogen photoionisation rate, and the optical depth due to Thomson scattering of the cosmic microwave background photons measured by the Planck satellite.

We present a comprehensive review of keV-scale sterile neutrino Dark Matter, collecting views and insights from all disciplines involved - cosmology, astrophysics, nuclear, and particle physics - in each case viewed from both theoretical and experimental/observational perspectives. After reviewing the role of active neutrinos in particle physics, astrophysics, and cosmology, we focus on sterile neutrinos in the context of the Dark Matter puzzle. Here, we first review the physics motivation for sterile neutrino Dark Matter, based on challenges and tensions in purely cold Dark Matter scenarios. We then round out the discussion by critically summarizing all known constraints on sterile neutrino Dark Matter arising from astrophysical observations, laboratory experiments, and theoretical considerations. In this context, we provide a balanced discourse on the possibly positive signal from X-ray observations. Another focus of the paper concerns the construction of particle physics models, aiming to explain how sterile neutrinos of keV-scale masses could arise in concrete settings beyond the Standard Model of elementary particle physics. The paper ends with an extensive review of current and future astrophysical and laboratory searches, highlighting new ideas and their experimental challenges, as well as future perspectives for the discovery of sterile neutrinos.

We estimate the systemic orbital kinematics of the Milky Way classical satellites and compare them with predictions from the \Lambda{} cold dark matter (\Lambda{}CDM) model derived from a semi-analytical galaxy formation model applied to high resolution cosmological N-body simulations. We find that the Galactic satellite system is atypical of \Lambda{}CDM systems. The subset of 10 Galactic satellites with proper motion measurements has a velocity anisotropy, \beta=-2.2$\pm$0.4, that lies in the 2.9% tail of the \Lambda{}CDM distribution. Individually, the Milky Way satellites have radial velocities that are lower than expected for their proper motions, with 9 out of the 10 having at most 20% of their orbital kinetic energy invested in radial motion. Such extreme values are expected in only 1.5% of \Lambda{}CDM satellites systems. This tangential motion excess is unrelated to the existence of a Galactic "disc of satellites". We present theoretical predictions for larger satellite samples that may become available as more proper motion measurements are obtained.

We analyse the properties of the HI gas distribution in the Auriga project, a set of magnetohydrodynamic cosmological simulations performed with the moving-mesh code AREPO and a physics model for galaxy formation that succeeds in forming realistic late-type galaxies in the 30 Milky Way-sized haloes simulated in this project. We use a simple approach to estimate the neutral hydrogen fraction in our simulation set, which treats low-density and star-forming gas separately, and we explore two different prescriptions to subtract the contribution of molecular hydrogen from the total HI content. The HI gas in the vast majority of the systems forms extended discs although more disturbed morphologies are present. Notwithstanding the general good agreement with observed HI properties -- such as radial profiles and the mass-diameter relation -- the Auriga galaxies are systematically larger and more gas-rich than typical nearby galaxies. Interestingly, the amount of HI gas outside the disc plane correlates with the star formation rate, consistent with a picture where most of this extra-planar HI gas originates from a fountain-like flow. Our findings are robust with respect to the different assumptions adopted for computing the molecular hydrogen fraction and do not vary significantly over a wide range of numerical resolution. The HI modelling introduced in this paper can be used in future work to build artificial interferometric HI data cubes, allowing an even closer comparison of the gas dynamics in simulated galaxies with observations.

We introduce a suite of thirty cosmological magneto-hydrodynamical zoom simulations of the formation of galaxies in isolated Milky Way mass dark haloes. These were carried out with the moving mesh code \textlcsc{AREPO}, together with a comprehensive model for galaxy formation physics, including AGN feedback and magnetic fields, which produces realistic galaxy populations in large cosmological simulations. We demonstrate that our simulations reproduce a wide range of present-day observables, in particular, two component disc dominated galaxies with appropriate stellar masses, sizes, rotation curves, star formation rates and metallicities. We investigate the driving mechanisms that set present-day disc sizes/scale lengths, and find that they are related to the angular momentum of halo material. We show that the largest discs are produced by quiescent mergers that inspiral into the galaxy and deposit high angular momentum material into the pre-existing disc, simultaneously increasing the spin of dark matter and gas in the halo. More violent mergers and strong AGN feedback play roles in limiting disc size by destroying pre-existing discs and by suppressing gas accretion onto the outer disc, respectively. The most important factor that leads to compact discs, however, is simply a low angular momentum for the halo. In these cases, AGN feedback plays an important role in limiting central star formation and the formation of a massive bulge.

The defining characteristic of the cold dark matter (CDM) hypothesis is the presence of a very large number of low-mass haloes, too small to have made a visible galaxy. Other hypotheses for the nature of the dark matter, such as warm dark matter (WDM), predict a much smaller number of such low-mass haloes. Strong lensing systems offer the possibility of detecting small-mass haloes through the distortions they induce in the lensed image. Here we show that the main contribution to the image distortions comes from haloes along the line of sight rather than subhaloes in the lens as has normally been assumed so far. These interlopers enhance the differences between the predictions of CDM and WDM models. We derive the total perturber mass function, including both subhaloes and interlopers, and show that measurements of approximately 20 strong lens systems with a detection limit of $M_{\rm low}=10^7 h^{-1} M_{\odot}$ would distinguish (at 3 sigma) between CDM and a WDM model consisting of 7 keV sterile neutrinos such as those required to explain the recently detected 3.5 keV X-ray emission line from the centres of galaxies and clusters.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

The Sculptor dwarf spheroidal galaxy contains two distinct stellar populations: one metal-rich and the other metal-poor. Several authors have argued that in order for these two populations to reside in the same gravitational potential, the dark matter halo must have a core similar to that observed in the stellar count profile. This would rule out the cuspy Navarro-Frenk-White (NFW) density profiles predicted for halos and subhalos by dark matter only simulations of the Lambda Cold Dark Matter (Lambda-CDM) cosmological model. We present a new theoretical framework to analyse stellar count and velocity observations in a self-consistent manner based on separable models, f(E,J)=g(J)h(E), for the distribution function of an equilibrium spherical system. We use this machinery to analyse available photometric and kinematic data for the two stellar populations in Sculptor. We find, contrary to some previous claims, that the data are consistent with populations in equilibrium within an NFW dark matter potential with structural parameters in the range expected in Lambda-CDM. Our solution gives a maximum circular velocity for Sculptor between 20 and 35 km/s. We discuss why some previous authors came to a different conclusion.

We examine the baryon content of low-mass {\Lambda}CDM halos $(10^8<M_{200}/{\rm M_\odot}<5\times 10^{9})$ using the APOSTLE cosmological hydrodynamical simulations. Most of these systems are free of stars and have a gaseous content set by the combined effects of cosmic reionization, which imposes a mass-dependent upper limit, and of ram pressure stripping, which reduces it further in high-density regions. Halos mainly affected by reionization (RELHICs; REionization-Limited HI Clouds) inhabit preferentially low-density regions and make up a population where the gas is in hydrostatic equilibrium with the dark matter potential and in thermal equilibrium with the ionizing UV background. Their thermodynamic properties are well specified, and their gas density and temperature profiles may be predicted in detail. Gas in RELHICs is nearly fully ionized but with neutral cores that span a large range of HI masses and column densities and have negligible non-thermal broadening. We present predictions for their characteristic sizes and central column densities: the massive tail of the distribution should be within reach of future blind HI surveys. Local Group RELHICs (LGRs) have some properties consistent with observed Ultra Compact High Velocity Clouds (UCHVCs) but the sheer number of the latter suggests that most UCHVCs are not RELHICs. Our results suggest that LGRs (i) should typically be beyond 500 kpc from the Milky Way or M31; (ii) have positive Galactocentric radial velocities; (iii) HI sizes not exceeding 1 kpc, and (iv) should be nearly round. The detection and characterization of RELHICs would offer a unique probe of the small-scale clustering of cold dark matter.

We use the APOSTLE cosmological hydrodynamic simulations to examine the effects of tidal stripping on cold dark matter (CDM) sub haloes that host three of the most luminous Milky Way (MW) dwarf satellite galaxies: Fornax, Sculptor, and Leo I. We identify simulated satellites that match the observed spatial and kinematic distributions of stars in these galaxies, and track their evolution after infall. We find $\sim$ 30$\%$ of subhaloes hosting satellites with present-day stellar mass $10^6$-$10^8$ $M_{\odot}$ experience $>20\%$ stellar mass loss after infall. Fornax analogues have earlier infall times compared to Sculptor and Leo I analogues. Star formation in Fornax analogues continues for $\sim3$-$6$ Gyr after infall, whereas Sculptor and Leo I analogues stop forming stars $< 2$-$3$ Gyr after infall. Fornax analogues typically show more significant stellar mass loss and exhibit stellar tidal tails, whereas Sculptor and Leo I analogues, which are more deeply embedded in their host DM haloes at infall, do not show substantial mass loss due to tides. When additionally comparing the orbital motion of the host subaloes to the measured proper motion of Fornax we find the matching more difficult; host subhaloes tend to have pericentres smaller than that measured for Fornax itself. From the kinematic and orbital data, we estimate that Fornax has lost $10-20\%$ of its infall stellar mass. Our best estimate for the surface brightness of a stellar tidal stream associated with Fornax is $\Sigma \sim$ 32.6 mag $ {\rm arcsec^{-2}}$, which may be detectable with deep imaging surveys such as DES and LSST.

The sterile neutrino is a viable dark matter candidate that can be produced in the early Universe via non-equilibrium processes, and would therefore possess a highly non-thermal spectrum of primordial velocities. In this paper we analyse the process of structure formation with this class of dark matter particles. To this end we construct primordial dark matter power spectra as a function of the lepton asymmetry, $L_6$, that is present in the primordial plasma and leads to resonant sterile neutrino production. We compare these power spectra with those of thermally produced dark matter particles and show that resonantly produced sterile neutrinos are much colder than their thermal relic counterparts. We also demonstrate that the shape of these power spectra is not determined by the free-streaming scale alone. We then use the power spectra as an input for semi-analytic models of galaxy formation in order to predict the number of luminous satellite galaxies in a Milky Way-like halo. By assuming that the mass of the Milky Way halo must be no more than $2\times10^{12}M_{\odot}$ (the adopted upper bound based on current astronomical observations) we are able to constrain the value of $L_6$ for $M_s\le 8$~keV. We also show that the range of $L_6$ that is in best agreement with the 3.5~keV line (if produced by decays of 7~keV sterile neutrino) requires that the Milky Way halo has a mass no smaller than $1.5\times10^{12}M_{\odot}$. Finally, we compare the power spectra obtained by direct integration of the Boltzmann equations for a non-resonantly produced sterile neutrino with the fitting formula of Viel~et~al. and find that the latter significantly underestimates the power amplitude on scales relevant to satellite galaxies.

Galaxies fall into two clearly distinct types: `blue-sequence' galaxies that are rapidly forming young stars, and `red-sequence' galaxies in which star formation has almost completely ceased. Most galaxies more massive than $3\times10^{10} M_\odot$ follow the red-sequence while less massive central galaxies lie on the blue sequence. We show that these sequences are created by a competition between star formation-driven outflows and gas accretion on to the supermassive black hole at the galaxy's center. We develop a simple analytic model for this interaction. In galaxies less massive than $3\times10^{10} M_\odot$, young stars and supernovae drive a high entropy outflow that is more buoyant that any diffuse corona. The outflow balances the rate of gas inflow, preventing high gas densities building up in the central regions. More massive galaxies, however, are surrounded by a hot corona. We argue that above a halo mass of $\sim 10^{12} M_\odot$, the supernova-driven outflow is no longer buoyant and star formation is unable to prevent the build up of gas in the central regions. This triggers a strongly non-linear response from the black hole. Its accretion rate rises rapidly, heating the galaxy's corona, disrupting the incoming supply of cool gas and starving the galaxy of the fuel for star formation. The host galaxy makes a transition to the red sequence, and further growth predominantly occurs through galaxy mergers. We show that the analytic model provides a good description of galaxy evolution in the EAGLE hydrodynamic simulations, and demonstrate that, so long as star formation-driven outflows are present, the transition mass scale is almost independent of subgrid parameter choice. The transition mass disappears entirely, however, if star formation driven outflows are absent.

Massive early-type galaxies have higher metallicities and higher ratios of $\alpha$ elements to iron than their less massive counterparts. Reproducing these correlations has long been a problem for hierarchical galaxy formation theory, both in semi-analytic models and cosmological hydrodynamic simulations. We show that a simulation in which gas cooling in massive dark haloes is quenched by radio-mode active galactic nuclei (AGNs) feedback naturally reproduces the observed trend between $\alpha$/Fe and the velocity dispersion of galaxies, $\sigma$. The quenching occurs earlier for more massive galaxies. Consequently, these galaxies complete their star formation before $\alpha$/Fe is diluted by the contribution from type Ia supernovae. For galaxies more massive than $\sim 10^{11}~M_\odot$ whose $\alpha$/Fe correlates positively with stellar mass, we find an inversely correlated mass-metallicity relation. This is a common problem in simulations in which star formation in massive galaxies is quenched either by quasar- or radio-mode AGN feedback. The early suppression of gas cooling in progenitors of massive galaxies prevents them from recapturing enriched gas ejected as winds. Simultaneously reproducing the [$\alpha$/Fe]-$\sigma$ relation and the mass-metallicity relation is, thus, difficult in the current framework of galaxy formation.

We use the Copernicus Complexio (COCO) high resolution $N$-body simulations to investigate differences in the properties of small-scale structures in the standard cold dark matter (CDM) model and in a model with a cutoff in the initial power spectrum of density fluctuations consistent with both a thermally produced warm dark matter (WDM) particle or a sterile neutrino with mass 7 keV and leptogenesis parameter $L_6=8.7$. The latter corresponds to the "coldest" model with this sterile neutrino mass compatible with the identification of the recently detected 3.5 keV X-ray line as resulting from particle decay. CDM and WDM predict very different number densities of subhaloes with mass $\leq 10^9\,h^{-1}\,M_\odot$ although they predict similar, nearly universal, normalised subhalo radial density distributions. Haloes and subhaloes in both models have cuspy NFW profiles, but WDM subhaloes below the cutoff scale in the power spectrum (corresponding to maximum circular velocities $V_{\mathrm{max}}^{z=0} \leq50~\mathrm{kms}^{-1}$) are less concentrated than their CDM counterparts. We make predictions for observable properties using the GALFORM semi-analytic model of galaxy formation. Both models predict Milky Way satellite luminosity functions consistent with observations, although the WDM model predicts fewer very faint satellites. This model, however, predicts slightly more UV bright galaxies at redshift $z>7$ than CDM, but both are consistent with observations. Gravitational lensing offers the best prospect of distinguishing between the models.

We analyse the mass assembly of central galaxies in the EAGLE hydrodynamical simulations. We build merger trees to connect galaxies to their progenitors at different redshifts and characterize their assembly histories by focusing on the time when half of the galaxy stellar mass was assembled into the main progenitor. We show that galaxies with stellar mass $M_*<10^{10.5}M_{\odot}$ assemble most of their stellar mass through star formation in the main progenitor (`in-situ' star formation). This can be understood as a consequence of the steep rise in star formation efficiency with halo mass for these galaxies. For more massive galaxies, however, an increasing fraction of their stellar mass is formed outside the main progenitor and subsequently accreted. Consequently, while for low-mass galaxies the assembly time is close to the stellar formation time, the stars in high-mass galaxies typically formed long before half of the present-day stellar mass was assembled into a single object, giving rise to the observed anti-hierarchical downsizing trend. In a typical present-day $M_*\geq10^{11}M_{\odot}$ galaxy, around $20\%$ of the stellar mass has an external origin. This fraction decreases with increasing redshift. Bearing in mind that mergers only make an important contribution to the stellar mass growth of massive galaxies, we find that the dominant contribution comes from mergers with galaxies of mass greater than one tenth of the main progenitor's mass. The galaxy merger fraction derived from our simulations agrees with recent observational estimates.

We examine the properties of barred disc galaxies in a LCDM cosmological hydrodynamical simulation from the EAGLE project. Our study follows the formation of 269 discs identified at z = 0 in the stellar mass range 10.6 < log Mstr /M < 11. These discs show a wide range of bar strengths, from unbarred discs to weak bars to strongly barred systems (= 20%). Bars in these systems develop after redshift = 1.3, on timescales that depend sen- sitively on the strength of the pattern. Strong bars develop relatively quickly (in a few Gyr, = 10 disc rotation periods) in systems that are disc dominated, gas poor, and have declining rotation curves. Weak bars develop more slowly in systems where the disc is less gravitation- ally important, and are still growing at z = 0. Unbarred galaxies are comparatively gas-rich discs whose rotation speeds do not exceed the maximum circular velocity of the halos they inhabit. Bar lengths compare favourably with observations, ranging from 0.2 to 0.8 times the radius containing 90% of the stars. Bars slow down remarkably quickly as they grow, causing the inner regions of the surrounding dark halo to expand. At z = 0 strong bars have corota- tion radii roughly ten times the bar length. Such slow bars are inconsistent with the few cases where pattern speeds have been measured or inferred observationally, a discrepancy that, if confirmed, might prove a challenge for disc galaxy formation in LCDM.

We examine the spatial distribution of the oldest and most metal poor stellar populations of Milky Way-sized galaxies using the APOSTLE cosmological hydrodynamical simulations of the Local Group. In agreement with earlier work, we find strong radial gradients in the fraction of the oldest (tform < 0.8 Gyr) and most metal poor ([Fe/H]< -2.5) stars, both of which increase outwards. The most metal poor stars form over an extended period of time; half of them form after z = 5.3, and the last 10% after z = 2.8. The age of the metal poor stellar population also shows significant variation with environment; a high fraction of them are old in the galaxy's central regions and an even higher fraction in some individual dwarf galaxies, with substantial scatter from dwarf to dwarf. Overall, over half of the stars that belong to both the oldest and most metal-poor population are found outside the solar circle. Somewhat counter-intuitively, we find that dwarf galaxies with a large fraction of metal poor stars that are very old are systems where metal poor stars are relatively rare, but where a substantial old population is present. Our results provide guidance for interpreting the results of surveys designed to hunt for the earliest and most pristine stellar component of our Milky Way.

We derive the expression for the observed redshift in the weak field limit in the observer's past light cone, including all relativistic terms up to second order in velocity. We then apply it to compute the cluster-galaxy cross-correlation functions (CGCF) using N-body simulations. The CGCF is asymmetric along the line of sight (LOS) owing to the presence of the small second order terms such as the gravitational redshift (GRedshift). We identify two systematics in the modelling of the GRedshift signal in stacked clusters. First, it is affected by the morphology of dark matter haloes and the large-scale cosmic-web. The non-spherical distribution of galaxies around the central halo and the presence of neighbouring clusters systematically reduce the GRedshift signal. This bias is approximately 20% for $M_{\rm min}\simeq 10^{14} {\rm M_{\odot}}/h$, and is more than 50% for haloes with $M_{\rm min}\simeq 2\times 10^{13} {\rm M_{\odot}}/h$ at $r>$4 Mpc/$h$. Second, the best-fit gravitational redshift profiles as well as the profiles of all other relativistic terms are found to be significantly different in velocity space compared to their real space versions. We discuss some subtleties relating to these effects in velocity space. We also find that the S/N of the GRedshift signal increases with decreasing halo mass.

The predicted abundance and properties of the low-mass substructures embedded inside larger dark matter haloes differ sharply among alternative dark matter models. Too small to host galaxies themselves, these subhaloes may still be detected via gravitational lensing, or via perturbations of the Milky Way's globular cluster streams and its stellar disk. Here we use the Apostle cosmological simulations to predict the abundance and the spatial and velocity distributions of subhaloes in the range 10^6.5-10^8.5 solar masses inside haloes of mass ~ 10^12 solar masses in LCDM. Although these subhaloes are themselves devoid of baryons, we find that baryonic effects are important. Compared to corresponding dark matter only simulations, the loss of baryons from subhaloes and stronger tidal disruption due to the presence of baryons near the centre of the main halo, reduce the number of subhaloes by ~ 1/4 to 1/2, independently of subhalo mass, but increasingly towards the host halo centre. We also find that subhaloes have non-Maxwellian orbital velocity distributions, with centrally rising velocity anisotropy and positive velocity bias which reduces the number of low-velocity subhaloes, particularly near the halo centre. We parameterise the predicted population of subhaloes in terms of mass, galactocentric distance, and velocities. We discuss implications of our results for the prospects of detecting dark matter substructures and for possible inferences about the nature of dark matter.

We present the public data release of halo and galaxy catalogues extracted from the EAGLE suite of cosmological hydrodynamical simulations of galaxy formation. These simulations were performed with an enhanced version of the GADGET code that includes a modified hydrodynamics solver, time-step limiter and subgrid treatments of baryonic physics, such as stellar mass loss, element-by-element radiative cooling, star formation and feedback from star formation and black hole accretion. The simulation suite includes runs performed in volumes ranging from 25 to 100 comoving megaparsecs per side, with numerical resolution chosen to marginally resolve the Jeans mass of the gas at the star formation threshold. The free parameters of the subgrid models for feedback are calibrated to the redshift z=0 galaxy stellar mass function, galaxy sizes and black hole mass - stellar mass relation. The simulations have been shown to match a wide range of observations for present-day and higher-redshift galaxies. The raw particle data have been used to link galaxies across redshifts by creating merger trees. The indexing of the tree produces a simple way to connect a galaxy at one redshift to its progenitors at higher redshift and to identify its descendants at lower redshift. In this paper we present a relational database which we are making available for general use. A large number of properties of haloes and galaxies and their merger trees are stored in the database, including stellar masses, star formation rates, metallicities, photometric measurements and mock gri images. Complex queries can be created to explore the evolution of more than 10^5 galaxies, examples of which are provided in appendix. (abridged)

Theoretical models of galaxy formation based on the cold dark matter cosmogony typically require strong feedback from supernova (SN) explosions in order to reproduce the Milky Way satellite galaxy luminosity function and the faint end of the field galaxy luminosity function. However, too strong a SN feedback also leads to the universe reionizing too late, and the metallicities of Milky Way satellites being too low. The combination of these four observations therefore places tight constraints on SN feedback. We investigate these constraints using the semi-analytical galaxy formation model galform. We find that these observations favour a SN feedback model in which the feedback strength evolves with redshift. We find that, for our best fit model, half of the ionizing photons are emitted by galaxies with rest-frame far-UV absolute magnitudes $M_{\rm AB}(1500{\rm \AA})<-17.5$, which implies that already observed galaxy populations contribute about half of the photons responsible for reionization. The $z=0$ descendants of these galaxies are mainly galaxies with stellar mass $M_*>10^{10}\,{\rm M}_{\odot}$ and preferentially inhabit halos with mass $M_{\rm halo}>10^{13}\,{\rm M}_{\odot}$.

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

We examine the dark matter content of satellite galaxies in Lambda-CDM cosmological hydrodynamical simulations of the Local Group from the APOSTLE project. We find excellent agreement between simulation results and estimates for the 9 brightest Galactic dwarf spheroidals (dSphs) derived from their stellar velocity dispersions and half-light radii. Tidal stripping plays an important role by gradually removing dark matter from the outside in, affecting in particular fainter satellites and systems of larger-than-average size for their luminosity. Our models suggest that tides have significantly reduced the dark matter content of Can Ven I, Sextans, Carina, and Fornax, a prediction that may be tested by comparing them with field galaxies of matching luminosity and size. Uncertainties in observational estimates of the dark matter content of individual dwarfs have been underestimated in the past, at times substantially. We use our improved estimates to revisit the `too-big-to-fail' problem highlighted in earlier N-body work. We reinforce and extend our previous conclusion that the APOSTLE simulations show no sign of this problem. The resolution does not require `cores' in the dark mass profiles, but, rather, relies on revising assumptions and uncertainties in the interpretation of observational data and accounting for `baryon effects' in the theoretical modelling.

The Tully-Fisher relation (TFR) links the stellar mass of a disk galaxy, $M_{\rm str}$, to its rotation speed: it is well approximated by a power law, shows little scatter, and evolves weakly with redshift. The relation has been interpreted as reflecting the mass-velocity scaling ($M\propto V^3$) of dark matter halos, but this interpretation has been called into question by abundance-matching (AM) models, which predict the galaxy-halo mass relation to be non-monotonic and rapidy evolving. We study the TFR of luminous spirals and its relation to AM using the EAGLE set of $\Lambda$CDM cosmological simulations. Matching both relations requires disk sizes to satisfy constraints given by the concentration of halos and their response to galaxy assembly. EAGLE galaxies approximately match these constraints and show a tight mass-velocity scaling that compares favourably with the observed TFR. The TFR is degenerate to changes in galaxy formation efficiency and the mass-size relation; simulations that fail to match the galaxy stellar mass function may fit the observed TFR if galaxies follow a different mass-size relation. The small scatter in the simulated TFR results because, at fixed halo mass, galaxy mass and rotation speed correlate strongly, scattering galaxies along the main relation. EAGLE galaxies evolve with lookback time following approximately the prescriptions of AM models and the observed mass-size relation of bright spirals, leading to a weak TFR evolution consistent with observation out to $z=1$. $\Lambda$CDM models that match both the abundance and size of galaxies as a function of stellar mass have no difficulty reproducing the observed TFR and its evolution.

We investigate the evolution of supermassive black holes in the `Evolution and Assembly of GaLaxies and their Environments' (EAGLE) cosmological hydrodynamic simulations. The largest of the EAGLE volumes covers a $(100 \,\rm cMpc)^3$ and includes state-of-the-art physical models for star formation and black hole growth that depend only on local gas properties. We focus on the black hole mass function, Eddington ratio distribution and the implied duty cycle of nuclear activity. The simulation is broadly consistent with observational constraints on these quantities. In order to make a more direct comparison with observational data, we calculate the soft and hard X-ray luminosity functions of the active galactic nuclei (AGN). Between redshifts $0$ and $1$, the simulation is in agreement with data. At higher redshifts, the simulation tends to underpredict the luminosities of the brightest observed AGN. This may be due to the limited volume of the simulation, or a fundamental deficiency of the underlying model. It seems unlikely that additional unresolved variability can account for this difference. The simulation shows a similar `downsizing' of the AGN population as seen in observational surveys.

The spatial distribution of the satellite populations of the Milky Way and Andromeda are puzzling in that they are nearly perpendicular to the disks of their central galaxies. To understand the origin of such configurations we study the alignment of the central galaxy, satellite system and dark matter halo in the largest of the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) simulation. We find that centrals and their satellite systems tend to be well aligned with their haloes, with a median misalignment angle of $33^{\circ}$ in both cases. While the centrals are better aligned with the inner $10$ kpc halo, the satellite systems are better aligned with the entire halo indicating that satellites preferentially trace the outer halo. The central - satellite alignment is weak (median misalignment angle of $52^{\circ}$) and we find that around $20\%$ of systems have a misalignment angle larger than $78^{\circ}$, which is the value for the Milky Way. The central - satellite alignment is a consequence of the tendency of both components to align with the dark matter halo. As a consequence, when the central is parallel to the satellite system, it also tends to be parallel to the halo. In contrast, if the central is perpendicular to the satellite system, as in the case of the Milky Way and Andromeda, then the central - halo alignment is much weaker. Dispersion-dominated (spheroidal) centrals have a stronger alignment with both their halo and their satellites than rotation-dominated (disk) centrals. We also found that the halo, the central galaxy and the satellite system tend to be aligned with the surrounding large-scale distribution of matter, with the halo being the better aligned of the three.

We expand our previous study on the relationship between changes in the orientation of the angular momentum vector of dark matter haloes ("spin flips") and changes in their mass (Bett & Frenk 2012), to cover the full range of halo masses in a simulation cube of length 100 $h^{-1}$ Mpc. Since strong disturbances to a halo (such as might be indicated by a large change in the spin direction) are likely also to disturb the galaxy evolving within, spin flips could be a mechanism for galaxy morphological transformation without involving major mergers. We find that 35% of haloes have, at some point in their lifetimes, had a spin flip of at least $45\deg$ that does not coincide with a major merger. Over 75% of large spin flips coincide with non-major mergers; only a quarter coincide with major mergers. We find a similar picture for changes to the inner-halo spin orientation, although here there is an increased likelihood of a flip occurring. Changes in halo angular momentum orientation, and other such measures of halo perturbation, are therefore very important quantities to consider, in addition to halo mergers, when modelling the formation and evolution of galaxies and confronting such models with observations.

We investigate the presence and importance of dark matter discs in a sample of 24 simulated Milky Way galaxies in the APOSTLE project, part of the EAGLE programme of hydrodynamic simulations in Lambda-CDM cosmology. It has been suggested that a dark disc in the Milky Way may boost the dark matter density and modify the velocity modulus relative to a smooth halo at the position of the Sun, with ramifications for direct detection experiments. From a kinematic decomposition of the dark matter and a real space analysis of all 24 halos, we find that only one of the simulated Milky Way analogues has a detectable dark disc component. This unique event was caused by a merger at late time with an LMC-mass satellite at very low grazing angle. Considering that even this rare scenario only enhances the dark matter density at the solar radius by 35% and affects the high energy tail of the dark matter velocity distribution by less than 1%, we conclude that the presence of a dark disc in the Milky Way is unlikely, and is very unlikely to have a significant effect on direct detection experiments.

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

We introduce a new method for stacking voids and deriving their profile that greatly increases the potential of voids as a tool for precision cosmology. Given that voids are highly non-spherical and have most of their mass at their edge, voids are better described relative to their boundary rather than relative to their centre, as in the conventional spherical stacking approach. The boundary profile is obtained by computing the distance of each volume element from the void boundary. Voids can then be stacked and their profiles computed as a function of this boundary distance. This approach enhances the weak lensing signal of voids, both shear and convergence, by a factor of two when compared to the spherical stacking method. It also results in steeper void density profiles that are characterised by a very slow rise inside the void and a pronounced density ridge at the void boundary. The resulting boundary density profile is self-similar when rescaled by the thickness of the density ridge, implying that the average rescaled profile is independent of void size. The boundary velocity profile is characterized by outflows in the inner regions whose amplitude scales with void size, and by a strong inflow into the filaments and walls delimiting the void. This new picture enables a straightforward discrimination between collapsing and expanding voids both for individual objects as well as for stacked samples.

The hierarchical nature of LCDM suggests that the Magellanic Clouds must have been surrounded by a number of satellites before their infall into the Milky Way. Many of those satellites should still be in close proximity to the Clouds, but some could have dispersed ahead/behind the Clouds along their Galactic orbit. Either way, prior association with the Clouds results in strong restrictions on the present-day positions and velocities of candidate Magellanic satellites: they must lie close to the nearly-polar orbital plane of the Magellanic stream, and their distances and radial velocities must follow the latitude dependence expected for a tidal stream with the Clouds at pericenter. We use a cosmological numerical simulation of the disruption of a massive subhalo in a Milky Way-sized LCDM halo to test whether any of the 20 dwarfs recently-discovered in the DES, SMASH, Pan-STARRS, and ATLAS surveys are truly associated with the Clouds. Of the 6 systems with kinematic data, only Hydra II and Hor 1 have distances and radial velocities consistent with a Magellanic origin. Of the remaining dwarfs, six (Hor 2, Eri 3, Ret 3, Tuc 4, Tuc 5, and Phx 2) have positions and distances consistent with a Magellanic origin, but kinematic data are needed to substantiate that possibility. Conclusive evidence for association would require proper motions to constrain the orbital angular momentum direction, which, for true Magellanic satellites, must coincide with that of the Clouds. We use this result to predict radial velocities and proper motions for all new dwarfs. Our results are relatively insensitive to the assumption of first or second pericenter for the Clouds.

We investigate the process of reionisation in a model in which the dark matter is a warm elementary particle such as a sterile neutrino. We focus on models that are consistent with the dark matter decay interpretation of the recently detected line at 3.5 keV in the X-ray spectra of galaxies and clusters. In warm dark matter models the primordial spectrum of density perturbations has a cut-off on the scale of dwarf galaxies. Structure formation therefore begins later than in the standard cold dark matter (CDM) model and very few objects form below the cut-off mass scale. To calculate the number of ionising photons, we use the Durham semi-analytic model of galaxy formation, GALFORM. We find that even the most extreme 7 keV sterile neutrino we consider is able to reionise the Universe early enough to be compatible with the bounds on the epoch of reionisation from Planck. This, perhaps surprising, result arises from the rapid build-up of high redshift galaxies in the sterile neutrino models which is also reflected in a faster evolution of their far-UV luminosity function between $10>z>7$ than in CDM. The dominant sources of ionising photons are systematically more massive in the sterile neutrino models than in CDM. As a consistency check on the models, we calculate the present-day luminosity function of satellites of Milky Way-like galaxies. When the satellites recently discovered in the DES survey are taken into account, strong constraints are placed on viable sterile neutrino models.

We study the implications of galaxy formation on dark matter direct detection using high resolution hydrodynamic simulations of Milky Way-like galaxies simulated within the EAGLE and APOSTLE projects. We identify Milky Way analogues that satisfy observational constraints on the Milky Way rotation curve and total stellar mass. We then extract the dark matter density and velocity distribution in the Solar neighbourhood for this set of Milky Way analogues, and use them to analyse the results of current direct detection experiments. For most Milky Way analogues, the event rates in direct detection experiments obtained from the best fit Maxwellian distribution (with peak speed of 223 - 289 km/s) are similar to those obtained directly from the simulations. As a consequence, the allowed regions and exclusion limits set by direct detection experiments in the dark matter mass and spin-independent cross section plane shift by a few GeV compared to the Standard Halo Model, at low dark matter masses. For each dark matter mass, the halo-to-halo variation of the local dark matter density results in an overall shift of the allowed regions and exclusion limits for the cross section. However, the compatibility of the possible hints for a dark matter signal from DAMA and CDMS-Si and null results from LUX and SuperCDMS is not improved.

We use a suite of cosmological simulations to study the mass-concentration-redshift relation, $c({\rm M},z)$, of dark matter halos. Our simulations include standard $\Lambda$-cold dark matter (CDM) models, and additional runs with truncated power spectra, consistent with a thermal warm dark matter (WDM) scenario. We find that the mass profiles of CDM and WDM halos are self-similar and well approximated by the Einasto profile. The $c({\rm M},z)$ relation of CDM halos is monotonic: concentrations decrease with increasing virial mass at fixed redshift, and decrease with increasing redshift at fixed mass. The mass accretion histories (MAHs) of CDM halos are also scale-free, and can be used to infer concentrations directly. These results do not apply to WDM halos: their MAHs are not scale-free because of the characteristic scale imposed by the power-spectrum suppression. Further, the WDM $c({\rm M},z)$ relation is non-monotonic: concentrations peak at a mass scale dictated by the truncation scale, and decrease at higher and lower masses. We show that the assembly history of a halo can still be used to infer its concentration, provided that the total mass of its progenitors is considered (the "collapsed mass history"; CMH), rather than just that of its main ancestor. This exploits the scale-free nature of CMHs to derive a simple scaling that reproduces the mass-concentration-redshift relation of both CDM and WDM halos over a vast range of halo masses and redshifts. Our model therefore provides a robust account of the mass, redshift, cosmology and power spectrum dependence of dark matter halo concentrations.

We use the EAGLE galaxy formation simulation to study the effects of baryons on the power spectrum of the total matter and dark matter distributions and on the velocity fields of dark matter and galaxies. On scales $k{\stackrel{>}{{}_\sim}} 4{h\,{\rm Mpc}^{-1}}$ the effect of baryons on the amplitude of the total matter power spectrum is greater than $1\%$. The back-reaction of baryons affects the density field of the dark matter at the level of $\sim3\%$ on scales of $1\leq k/({h\,{\rm Mpc}^{-1}})\leq 5$. The dark matter velocity divergence power spectrum at $k{\stackrel{<}{{}_\sim}}0.5{h\,{\rm Mpc}^{-1}}$ is changed by less than $1\%$. The 2D redshift-space power spectrum is affected at the level of $\sim6\%$ at $|\vec{k}|{\stackrel{>}{{}_\sim}} 1{h\,{\rm Mpc}^{-1}}$ (for $\mu>0.5$), but for $|\vec{k}|\leq 0.4{h\,{\rm Mpc}^{-1}}$ it differs by less than $1\%$. We report vanishingly small baryonic velocity bias for haloes: the peculiar velocities of haloes with $M_{200}>3\times10^{11}{{\rm M}_{\odot}}$ (hosting galaxies with $M_{*}>10^9{{\rm M}_{\odot}}$) are affected at the level of at most $1~$km/s, which is negligible for $1\%$-precision cosmology. We caution that since EAGLE overestimates cluster gas fractions it may also underestimate the impact of baryons, particularly for the total matter power spectrum. Nevertheless, our findings suggest that for theoretical modelling of redshift space distortions and galaxy velocity-based statistics, baryons and their back-reaction can be safely ignored at the current level of observational accuracy. However, we confirm that the modelling of the total matter power spectrum in weak lensing studies needs to include realistic galaxy formation physics in order to achieve the accuracy required in the precision cosmology era.

The cold dark matter (CDM) cosmological model unambigously predicts that a large number of haloes should survive as subhaloes when they are accreted into a larger halo. The CDM model would be ruled out if such substructures were shown not to exist. By contrast, if the dark matter consists of warm particles (WDM), then below a threshold mass that depends on the particle mass far fewer substructures would be present. Finding subhaloes below a certain mass would then rule out warm particle masses below some value. Strong gravitational lensing provides a clean method to measure the subhalo mass function through distortions in the structure of Einstein rings and giant arcs.Using mock lensing observations constructed from high-resolution N-body simulations, we show that measurements of approximately 100 strong lens systems with a detection limit of $10^7 h^{-1} M_{\odot}$ would clearly distinguish CDM from WDM in the case where this consists of 7 keV sterile neutrinos such as those that might be responsible for the 3.5 keV X-ray emission line recently detected in galaxies and clusters.

A recent observational study of haloes of nearby Milky Way-like galaxies shows that only half (four out of eight) of the current sample exhibits strong negative metallicity ([Fe/H]) gradients. This is at odds with predictions from hydrodynamical simulations where such gradients are ubiquitous. In this Letter, we use high resolution cosmological hydrodynamical simulations to study the [Fe/H] distribution of galactic haloes. We find that kinematically selected stellar haloes, including both in-situ and accreted particles, have an oblate [Fe/H] distribution. Spherical [Fe/H] radial profiles show strong negative gradients within 100 kpc, in agreement with previous numerical results. However, the projected median [Fe/H] profiles along the galactic disc minor axis, typically obtained in observations, are significantly flatter. The median [Fe/H] values at a given radius are larger for the spherical profiles than for the minor axis profiles by as much as 0.4 dex within the inner 50 kpc. Similar results are obtained if only the accreted stellar component is considered indicating that the differences between spherical and minor axis profiles are not purely driven by heated disc star particles formed in situ. Our study highlights the importance of performing careful comparisons between models and observations of halo [Fe/H] distributions.

We introduce Copernicus Complexio (COCO), a high-resolution cosmological N-body simulation of structure formation in the $\Lambda{\rm CDM}{}$ model. COCO follows an approximately spherical region of radius $\sim 17.4h^{-1}\,{\rm Mpc}$ embedded in a much larger periodic cube that is followed at lower resolution. The high resolution volume has a particle mass of $1.135\times10^5h^{-1}{\rm M}_{\odot}$ (60 times higher than the Millennium-II simulation). COCO gives the dark matter halo mass function over eight orders of magnitude in halo mass; it forms $\sim 60$ haloes of galactic size, each resolved with about 10 million particles. We confirm the power-law character of the subhalo mass function, $\bar{N}(>\mu)\propto\mu^{-s}$, down to a reduced subhalo mass $M_{sub}/M_{200}\equiv\mu=10^{-6}$, with a best-fit power-law index, $s=0.94$, for hosts of mass $\langle M_{200}\rangle=10^{12}h^{-1}{\rm M}_{\odot}$. The concentration-mass relation of COCO haloes deviates from a single power law for masses $M_{200}<\textrm{a few}\times 10^{8}h^{-1}{\rm M}_{\odot}$, where it flattens, in agreement with results by Sanchez-Conde et al. The host mass invariance of the reduced maximum circular velocity function of subhaloes, $\nu\equiv V_{max}/V_{200}$, hinted at in previous simulations, is clearly demonstrated over five orders of magnitude in host mass. Similarly, we find that the average, normalised radial distribution of subhaloes is approximately universal (i.e. independent of subhalo mass), as previously suggested by the Aquarius simulations of individual haloes. Finally, we find that at fixed physical subhalo size, subhaloes in lower mass hosts typically have lower central densities than those in higher mass hosts.

Hydrogen in the Universe was (re)ionised between redshifts $z \approx 10$ and $z \approx 6$. The nature of the sources of the ionising radiation is hotly debated, with faint galaxies below current detection limits regarded as prime candidates. Here we consider a scenario in which ionising photons escape through channels punctured in the interstellar medium by outflows powered by starbursts. We take account of the observation that strong outflows occur only when the star formation density is sufficiently high, and estimate the galaxy-averaged escape fraction as a function of redshift and luminosity from the resolved star formation surface densities in the EAGLE cosmological hydrodynamical simulation. We find that the fraction of ionising photons that escape from galaxies increases rapidly with redshift, reaching values of 5-20 percent at $z > 6$, with the brighter galaxies having higher escape fractions. Combining the dependence of escape fraction on luminosity and redshift with the observed luminosity function, we demonstrate that galaxies emit enough ionising photons to match the existing constraints on reionisation while also matching the observed UV-background post-reionisation. Our findings suggest that galaxies above the current Hubble Space Telescope detection limit emit half of the ionising radiation required to reionise the Universe.

The scaling of disk galaxy rotation velocity with baryonic mass (the "Baryonic Tully-Fisher" relation, BTF) has long confounded galaxy formation models. It is steeper than the M ~ V^3 scaling relating halo virial masses and circular velocities and its zero point implies that galaxies comprise a very small fraction of available baryons. Such low galaxy formation efficiencies may in principle be explained by winds driven by evolving stars, but the tightness of the BTF relation argues against the substantial scatter expected from such vigorous feedback mechanism. We use the APOSTLE/EAGLE simulations to show that the BTF relation is well reproduced in LCDM simulations that match the size and number of galaxies as a function of stellar mass. In such models, galaxy rotation velocities are proportional to halo virial velocity and the steep velocity-mass dependence results from the decline in galaxy formation efficiency with decreasing halo mass needed to reconcile the CDM halo mass function with the galaxy luminosity function. Despite the strong feedback, the scatter in the simulated BTF is smaller than observed, even when considering all simulated galaxies and not just rotationally-supported ones. The simulations predict that the BTF should become increasingly steep at the faint end, although the velocity scatter at fixed mass should remain small. Observed galaxies with rotation speeds below ~40 km/s seem to deviate from this prediction. We discuss observational biases and modeling uncertainties that may help to explain this disagreement in the context of LCDM models of dwarf galaxy formation.

We study high-resolution hydrodynamic simulations of Milky Way type galaxies obtained within the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) project, and identify the those that best satisfy observational constraints on the Milky Way total stellar mass, rotation curve, and galaxy shape. Contrary to mock galaxies selected on the basis of their total virial mass, the Milky Way analogues so identified consistently exhibit very similar dark matter profiles inside the solar circle, therefore enabling more accurate predictions for indirect dark matter searches. We find in particular that high resolution simulated haloes satisfying observational constraints exhibit, within the inner few kiloparsecs, dark matter profiles shallower than those required to explain the so-called Fermi GeV excess via dark matter annihilation.

We explore the co-evolution of the specific angular momentum of dark matter haloes and the cold baryons that comprise the galaxies within. We study over two thousand central galaxies within the reference cosmological hydrodynamical simulation of the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) project. We employ a methodology within which the evolutionary history of a system is specified by the time-evolving properties of the Lagrangian particles that define it at z=0. We find a strong correlation between the evolution of the specific angular momentum of today's stars (cold gas) and that of the inner (whole) dark matter halo they are associated with. This link is particularly strong for the stars formed before the epoch of maximum expansion and subsequent collapse of the central dark matter halo (turnaround). Spheroids are typically assembled primarily from stars formed prior to turnaround, and are therefore destined to suffer a net loss of angular momentum associated with the strong merging activity during the assembly of the inner dark matter halo. Stellar discs retain their specific angular momentum since they are comprised of stars formed mainly after turnaround, from gas that mostly preserves the high specific angular momentum it acquired by tidal torques during the linear growth of the halo. Since the specific angular momentum loss of the stars is tied to the galaxy's morphology today, it may be possible to use our results to predict, statistically, the assembly history of a halo given the morphology of the galaxy it hosts.

N-body simulations suggest that the substructures that survive inside dark matter haloes follow universal distributions in mass and radial number density. We demonstrate that a simple analytical model can explain these subhalo distributions as resulting from tidal stripping which increasingly reduces the mass of subhaloes with decreasing halo-centric distance. As a starting point, the spatial distribution of subhaloes of any given infall mass is shown to be largely indistinguishable from the overall mass distribution of the host halo. Using a physically motivated statistical description of the amount of mass stripped from individual subhaloes, the model fully describes the joint distribution of subhaloes in final mass, infall mass and radius. As a result, it can be used to predict several derived distributions involving combinations of these quantities including, but not limited to, the universal subhalo mass function, the subhalo spatial distribution, the gravitational lensing profile, the dark matter annihilation radiation profile and boost factor. This model clarifies a common confusion when comparing the spatial distributions of galaxies and subhaloes, the so called "anti-bias", as a simple selection effect. We provide a Python code SubGen for populating haloes with subhaloes at http://icc.dur.ac.uk/data/

We develop two general methods to infer the gravitational potential of a system using steady-state tracers, i.e., tracers with a time-independent phase-space distribution. Combined with the phase-space continuity equation, the time independence implies a universal Orbital Probability Density Function (oPDF) $\mathrm{d} P(\lambda|{\rm orbit})\propto \mathrm{d} t$, where $\lambda$ is the coordinate of the particle along the orbit. The oPDF is equivalent to Jeans theorem, and is the key physical ingredient behind most dynamical modelling of steady-state tracers. In the case of a spherical potential, we develop a likelihood estimator that fits analytical potentials to the system, and a non-parametric method ("phase-mark") that reconstructs the potential profile, both assuming only the oPDF. The methods involve no extra assumptions about the tracer distribution function and can be applied to tracers with any arbitrary distribution of orbits, with possible extension to non-spherical potentials. The methods are tested on Monte Carlo samples of steady-state tracers in dark matter haloes to show that they are unbiased as well as efficient. A fully documented \textsc{C/Python} code implementing our method is freely available at a GitHub repository linked from \url{http://icc.dur.ac.uk/data/#oPDF}.

Using realistic cosmological simulations of Milky Way sized haloes, we study their dynamical state and the accuracy of inferring their mass profiles with steady-state models of dynamical tracers. We use a new method that describes the phase-space distribution of a steady-state tracer population in a spherical potential without any assumption regarding the distribution of their orbits. Applying the method to five haloes from the Aquarius $\Lambda$CDM N-body simulation, we find that dark matter particles are an accurate tracer that enables the halo mass and concentration parameters to be recovered with an accuracy of $5\%$. Assuming a potential profile of the NFW form does not significantly affect the fits in most cases, except for halo A whose density profile differs significantly from the NFW form, leading to a $30\%$ bias in the dynamically fitted parameters. The existence of substructures in the dark matter tracers only affects the fits by $\sim 1\%$. Applying the method to mock stellar haloes generated by a particle-tagging technique, we find the stars are farther from equilibrium than dark matter particles, yielding a systematic bias of $\sim 20\%$ in the inferred mass and concentration parameter. The level of systematic biases obtained from a conventional distribution function fit to stars is comparable to ours, while similar fits to DM tracers are significantly biased in contrast to our fits. In line with previous studies, the mass bias is much reduced near the tracer half-mass radius.

We obtain predictions for the properties of cold dark matter annihilation radiation using high resolution hydrodynamic zoom-in cosmological simulations of Milky Way-like galaxies (APOSTLE project) carried out as part of the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) programme. Galactic halos in the simulation have significantly different properties from those assumed in the "standard halo model" often used in dark matter detection studies. The formation of the galaxy causes a contraction of the dark matter halo, whose density profile develops a steeper slope than the Navarro-Frenk-White (NFW) profile between $r\approx1.5$ kpc and $r\approx10$ kpc. At smaller radii, $r\lesssim1.5$ kpc, the halos develop a flatter than NFW slope. This unexpected feature may be specific to our particular choice of subgrid physics model but nevertheless the dark matter density profiles agree within 30% as the mass resolution is increased by a factor 150. The inner regions of the halos are almost perfectly spherical (axis ratios $b/a > 0.97$ within $r=1$ kpc) and there is no offset larger than 45 pc between the centre of the stellar distribution and the centre of the dark halo. The morphology of the predicted dark matter annihilation radiation signal is in broad agreement with $\gamma$-ray observations at large Galactic latitudes ($b\gtrsim3^\circ$). At smaller angles, the inferred signal in one of our four galaxies is similar to that which is observed but it is significantly weaker in the other three.

We present the GAMA Panchromatic Data Release (PDR) constituting over 230deg$^2$ of imaging with photometry in 21 bands extending from the far-UV to the far-IR. These data complement our spectroscopic campaign of over 300k galaxies, and are compiled from observations with a variety of facilities including: GALEX, SDSS, VISTA, WISE, and Herschel, with the GAMA regions currently being surveyed by VST and scheduled for observations by ASKAP. These data are processed to a common astrometric solution, from which photometry is derived for 221,373 galaxies with r<19.8 mag. Online tools are provided to access and download data cutouts, or the full mosaics of the GAMA regions in each band. We focus, in particular, on the reduction and analysis of the VISTA VIKING data, and compare to earlier datasets (i.e., 2MASS and UKIDSS) before combining the data and examining its integrity. Having derived the 21-band photometric catalogue we proceed to fit the data using the energy balance code MAGPHYS. These measurements are then used to obtain the first fully empirical measurement of the 0.1-500$\mu$m energy output of the Universe. Exploring the Cosmic Spectral Energy Distribution (CSED) across three time-intervals (0.3-1.1Gyr, 1.1-1.8~Gyr and 1.8---2.4~Gyr), we find that the Universe is currently generating $(1.5 \pm 0.3) \times 10^{35}$ h$_{70}$ W Mpc$^{-3}$, down from $(2.5 \pm 0.2) \times 10^{35}$ h$_{70}$ W Mpc$^{-3}$ 2.3~Gyr ago. More importantly, we identify significant and smooth evolution in the integrated photon escape fraction at all wavelengths, with the UV escape fraction increasing from 27(18)% at z=0.18 in NUV(FUV) to 34(23)% at z=0.06. The GAMA PDR will allow for detailed studies of the energy production and outputs of individual systems, sub-populations, and representative galaxy samples at $z<0.5$. The GAMA PDR can be found at: http://gama-psi.icrar.org/

The Local Group of galaxies offer some of the most discriminating tests of models of cosmic structure formation. For example, observations of the Milky Way (MW) and Andromeda satellite populations appear to be in disagreement with N-body simulations of the "Lambda Cold Dark Matter" ({\Lambda}CDM) model: there are far fewer satellite galaxies than substructures in cold dark matter halos (the "missing satellites" problem); dwarf galaxies seem to avoid the most massive substructures (the "too-big-to-fail" problem); and the brightest satellites appear to orbit their host galaxies on a thin plane (the "planes of satellites" problem). Here we present results from APOSTLE (A Project Of Simulating The Local Environment), a suite of cosmological hydrodynamic simulations of twelve volumes selected to match the kinematics of the Local Group (LG) members. Applying the Eagle code to the LG environment, we find that our simulations match the observed abundance of LG galaxies, including the satellite galaxies of the MW and Andromeda. Due to changes to the structure of halos and the evolution in the LG environment, the simulations reproduce the observed relation between stellar mass and velocity dispersion of individual dwarf spheroidal galaxies without necessitating the formation of cores in their dark matter profiles. Satellite systems form with a range of spatial anisotropies, including one similar to that of the MW, confirming that such a configuration is not unexpected in {\Lambda}CDM. Finally, based on the observed velocity dispersion, size, and stellar mass, we provide new estimates of the maximum circular velocity for the halos of nine MW dwarf spheroidals.

The mass of the dark matter halo of the Milky Way can be estimated by fitting analytical models to the phase-space distribution of dynamical tracers. We test this approach using realistic mock stellar halos constructed from the Aquarius N-body simulations of dark matter halos in the $\Lambda$CDM cosmology. We extend the standard treatment to include a Navarro-Frenk-White (NFW) potential and use a maximum likelihood method to recover the parameters describing the simulated halos from the positions and velocities of their mock halo stars. We find that the estimate of halo mass is highly correlated with the estimate of halo concentration. The best-fit halo masses within the virial radius, $R_{200}$, are biased, ranging from a 40\% underestimate to a 5\% overestimate in the best case (when the tangential velocities of the tracers are included). There are several sources of bias. Deviations from dynamical equilibrium can potentially cause significant bias; deviations from spherical symmetry are relatively less important. Fits to stars at different galactocentric radii can give different mass estimates. By contrast, the model gives good constraints on the mass within the half-mass radius of tracers even when restricted to tracers within 60kpc. The recovered velocity anisotropies of tracers, $\beta$, are biased systematically, but this does not affect other parameters if tangential velocity data are used as constraints.

We study the formation of stellar haloes in three Milky Way-mass galaxies using cosmological smoothed particle hydrodynamics simulations, focusing on the subset of halo stars that form in situ, as opposed to those accreted from satellites. In situ stars in our simulations dominate the stellar halo out to 20 kpc and account for 30-40 per cent of its total mass. We separate in situ halo stars into three straightforward, physically distinct categories according to their origin: stars scattered from the disc of the main galaxy ('heated disc'), stars formed from gas smoothly accreted on to the halo ('smooth' gas) and stars formed in streams of gas stripped from infalling satellites ('stripped' gas). We find that most belong to the stripped gas category. Those originating in smooth gas outside the disc tend to form at the same time and place as the stripped-gas population, suggesting that their formation is associated with the same gas-rich accretion events. The scattered disc star contribution is negligible overall but significant in the Solar neighbourhood, where ~90 per cent of stars on eccentric orbits once belonged to the disc. However, the distinction between halo and thick disc in this region is highly ambiguous. The chemical and kinematic properties of the different components are very similar at the present day, but the global properties of the in situ halo differ substantially between the three galaxies in our study. In our simulations, the hierarchical buildup of structure is the driving force behind not only the accreted stellar halo, but also those halo stars formed in situ.

We use cosmological simulations to identify dark matter subhalo host candidates of the Fornax dwarf spheroidal galaxy using the stellar kinematic properties of Fornax. We consider cold dark matter (CDM), warm dark matter (WDM), and decaying dark matter (DDM) simulations for our models of structure formation. The subhalo candidates in CDM typically have smaller mass and higher concentrations at z = 0 than the corresponding candidates in WDM and DDM. We examine the formation histories of the ~ 100 Fornax candidate subhalos identified in CDM simulations and, using approximate luminosity-mass relationships for subhalos, we find two of these subhalos that are consistent with both the Fornax luminosity and kinematics. These two subhalos have a peak mass over ten times larger than their z = 0 mass. We suggest that in CDM the dark matter halo hosting Fornax must have been severely stripped of mass and that it had an infall time into the Milky Way of ~ 9 Gyr ago. In WDM, we find that candidate subhalos consistent with the properties of Fornax have a similar infall time and a similar degree of mass loss, while in DDM we find a later infall time of ~ 3 - 4 Gyr ago and significantly less mass loss. We discuss these results in the context of the Fornax star formation history, and show that these predicted subhalo infall times can be linked to different star formation quenching mechanisms. This emphasizes the links between the properties of the dark matter and the mechanisms that drive galaxy evolution.

Feedback from energy liberated by gas accretion onto black holes (BHs) is an attractive mechanism to explain the exponential cut-off at the massive end of the galaxy stellar mass function (SMF). Semi-analytic models of galaxy formation in which this form of feedback is assumed to suppress cooling in haloes where the gas cooling time is large compared to the dynamical time do indeed achieve a good match to the observed SMF. Furthermore, hydrodynamic simulations of individual halos in which gas is assumed to accrete onto the central BH at the Bondi rate have shown that a self-regulating regime is established in which the BH grows just enough to liberate an amount of energy comparable to the thermal energy of the halo. However, this process is efficient at suppressing the growth not only of massive galaxies but also of galaxies like the Milky Way, leading to disagreement with the observed SMF. The Bondi accretion rate, however, is inappropriate when the accreting material has angular momentum. We present an improved accretion model that takes into account the circularisation and subsequent viscous transport of infalling material and include it as a "subgrid" model in hydrodynamic simulations of the evolution of halos with a wide range of masses. The resulting accretion rates are generally low in low mass ($\lsim 10^{11.5} \msun$) halos, but show outbursts of Eddington-limited accretion during galaxy mergers. During outbursts these objects strongly resemble quasars. In higher mass haloes, gas accretion occurs continuously, typically at $~10$ % of the Eddington rate, which is conducive to the formation of radio jets. The resulting dependence of the accretion behaviour on halo mass induces a break in the relation between galaxy stellar mass and halo mass in these simulations that matches observations.

The detection of planar structures within the satellite systems of both the Milky Way (MW) and Andromeda (M31) has been reported as being in stark contradiction to the predictions of the standard cosmological model ($\Lambda$CDM). Given the ambiguity in defining a planar configuration, it is unclear how to interpret the low incidence of the MW and M31 planes in $\Lambda$CDM. We investigate the prevalence of satellite planes around galactic mass haloes identified in high resolution cosmological simulations. We find that planar structures are very common, and that ~10% of $\Lambda$CDM haloes have even more prominent planes than those present in the Local Group. While ubiquitous, the planes of satellite galaxies show a large diversity in their properties. This precludes using one or two systems as small scale probes of cosmology, since a large sample of satellite systems is needed to obtain a good measure of the object-to-object variation. This very diversity has been misinterpreted as a discrepancy between the satellite planes observed in the Local Group and $\Lambda$CDM predictions. In fact, ~10% of $\Lambda$CDM galactic haloes have planes of satellites that are as infrequent as the MW and M31 planes. The look-elsewhere effect plays an important role in assessing the detection significance of satellite planes and accounting for it leads to overestimating the significance level by a factor of 30 and 100 for the MW and M31 systems, respectively.

We report the alignment and shape of dark matter, stellar, and hot gas distributions in the EAGLE and cosmo-OWLS simulations. The combination of these state-of-the-art hydro-cosmological simulations enables us to span four orders of magnitude in halo mass ($11 < log_{10}(M_{200}/ [h^{-1}M_\odot]) < 15$), a wide radial range ($-2.3 < log_{10}(r/[h^{-1}Mpc ]) < 1.3$) and redshifts $0 < z < 1$. The shape parameters of the dark matter, stellar and hot gas distributions follow qualitatively similar trends: they become more aspherical (and triaxial) with increasing halo mass, radius and redshift. We measure the misalignment of the baryonic components (hot gas and stars) of galaxies with their host halo as a function of halo mass, radius, redshift, and galaxy type (centrals vs satellites and early- vs late-type). Overall, galaxies align well with the local distribution of the total (mostly dark) matter. However, the stellar distributions on galactic scales exhibit a median misalignment of about 45-50 degrees with respect to their host haloes. This misalignment is reduced to 25-30 degrees in the most massive haloes ($13 < log_{10}(M_{200}/ [h^{-1}M_\odot ]) < 15$). Half of the disc galaxies in the EAGLE simulations have a misalignment angle with respect to their host haloes larger than 40 degrees. We present fitting functions and tabulated values for the probability distribution of galaxy-halo misalignment to enable a straightforward inclusion of our results into models of galaxy formations based on purely collisionless N-body simulations.

We use a large sample of isolated dark matter halo pairs drawn from cosmological N-body simulations to identify candidate systems whose kinematics match that of the Local Group of Galaxies (LG). We find, in agreement with the "timing argument" and earlier work, that the separation and approach velocity of the Milky Way (MW) and Andromeda (M31) galaxies favour a total mass for the pair of ~ 5*10^12 M_sun. A mass this large, however, is difficult to reconcile with the small relative tangential velocity of the pair, as well as with the small deceleration from the Hubble flow observed for the most distant LG members. Halo pairs that match these three criteria have average masses a factor of ~2 times smaller than suggested by the timing argument, but with large dispersion, spanning more than a decade in mass. Guided by these results, we have selected 12 halo pairs with total mass in the range 1.6-3.6 *10^12 M_sun for the APOSTLE project (A Project Of Simulations of The Local Environment), a suite of resimulations at various numerical resolution levels (reaching up to ~10^4 M_sun per gas particle) that use the hydrodynamical code and subgrid physics developed for the EAGLE project. These simulations reproduce, by construction, the main kinematics of the MW-M31 pair, and produce satellite populations whose overall number, luminosities, and kinematics are in good agreement with observations of the MW and M31 companions. These diagnostics are sensitive to the total mass assumed for the MW-M31 pair; indeed, the LG satellite population would be quite difficult to reproduce for pair masses as high as indicated by the timing argument. The APOSTLE candidate systems thus provide an excellent testbed to confront directly many of the predictions of the Lambda-CDM cosmology with observations of our local Universe.

The recent detection of a 3.5 keV X-ray line from the centres of galaxies and clusters by Bulbul et al. (2014a) and Boyarsky et al. (2014a) has been interpreted as emission from the decay of 7 keV sterile neutrinos which could make up the (warm) dark matter (WDM). As part of the COpernicus COmplexio (COCO) programme, we investigate the properties of dark matter haloes formed in a high-resolution cosmological $N$-body simulation from initial conditions similar to those expected in a universe in which the dark matter consists of 7 keV sterile neutrinos. This simulation and its cold dark matter (CDM) counterpart have $\sim13.4$bn particles, each of mass $\sim 10^5\, h^{-1} M_\odot$, providing detailed information about halo structure and evolution down to dwarf galaxy mass scales. Non-linear structure formation on small scales ($M_{200}\, \leq\, 2 \times 10^9\,h^{-1}\,M_\odot$) begins slightly later in COCO-Warm than in COCO-Cold. The halo mass function at the present day in the WDM model begins to drop below its CDM counterpart at a mass $\sim 2 \times 10^{9}\,h^{-1}\,M_\odot$ and declines very rapidly towards lower masses so that there are five times fewer haloes of mass $M_{200}= 10^{8}\,h^{-1}\,M_\odot$ in COCO-Warm than in COCO-Cold. Halo concentrations on dwarf galaxy scales are correspondingly smaller in COCO-Warm, and we provide a simple functional form that describes its evolution with redshift. The shapes of haloes are similar in the two cases, but the smallest haloes in COCO-Warm rotate slightly more slowly than their CDM counterparts.

We examine the circular velocity profiles of galaxies in {\Lambda}CDM cosmological hydrodynamical simulations from the EAGLE and LOCAL GROUPS projects and compare them with a compilation of observed rotation curves of galaxies spanning a wide range in mass. The shape of the circular velocity profiles of simulated galaxies varies systematically as a function of galaxy mass, but shows remarkably little variation at fixed maximum circular velocity. This is especially true for low-mass dark matter-dominated systems, reflecting the expected similarity of the underlying cold dark matter haloes. This is at odds with observed dwarf galaxies, which show a large diversity of rotation curve shapes, even at fixed maximum rotation speed. Some dwarfs have rotation curves that agree well with simulations, others do not. The latter are systems where the inferred mass enclosed in the inner regions is much lower than expected for cold dark matter haloes and include many galaxies where previous work claims the presence of a constant density "core". The "cusp vs core" issue is thus better characterized as an "inner mass deficit" problem than as a density slope mismatch. For several galaxies the magnitude of this inner mass deficit is well in excess of that reported in recent simulations where cores result from baryon-induced fluctuations in the gravitational potential. We conclude that one or more of the following statements must be true: (i) the dark matter is more complex than envisaged by any current model; (ii) current simulations fail to reproduce the effects of baryons on the inner regions of dwarf galaxies; and/or (iii) the mass profiles of "inner mass deficit" galaxies inferred from kinematic data are incorrect.

We simulate the phase-space distribution of stellar mass in 9 massive Lambda-CDM galaxy clusters by applying the semi-analytic particle tagging method of Cooper et al. to the Phoenix suite of high-resolution N-body simulations (M200 = 7.5 to 33 x 10^14 Msol). The resulting surface brightness (SB) profiles of brightest cluster galaxies (BCGs) match well to observations. On average, stars formed in galaxies accreted by the BCG account for ~90 per cent of its total mass (the remainder is formed in situ). In circular BCG-centred apertures, the superposition of multiple debris clouds (each ~10 per cent of the total BCG mass) from different progenitors can result in an extensive outer diffuse component, qualitatively similar to a 'cD envelope'. These clouds typically originate from tidal stripping at z < 1 and comprise both streams and the extended envelopes of other massive galaxies in the cluster. Stars at very low SB contribute a significant fraction of the total cluster stellar mass budget: in the central 1 Mpc^2 of a z ~ 0.15 cluster imaged at SDSS-like resolution, our fiducial model predicts 80-95 per cent of stellar mass below a SB of mu_V = 26.5 mag arcsec^2 is associated with accreted stars in the envelope of the BCG. The ratio of BCG stellar mass (including this diffuse component) to total cluster stellar mass is ~30 per cent.

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

The Galaxy And Mass Assembly (GAMA) survey is one of the largest contemporary spectroscopic surveys of low-redshift galaxies. Covering an area of ~286 deg^2 (split among five survey regions) down to a limiting magnitude of r < 19.8 mag, we have collected spectra and reliable redshifts for 238,000 objects using the AAOmega spectrograph on the Anglo-Australian Telescope. In addition, we have assembled imaging data from a number of independent surveys in order to generate photometry spanning the wavelength range 1 nm - 1 m. Here we report on the recently completed spectroscopic survey and present a series of diagnostics to assess its final state and the quality of the redshift data. We also describe a number of survey aspects and procedures, or updates thereof, including changes to the input catalogue, redshifting and re-redshifting, and the derivation of ultraviolet, optical and near-infrared photometry. Finally, we present the second public release of GAMA data. In this release we provide input catalogue and targeting information, spectra, redshifts, ultraviolet, optical and near-infrared photometry, single-component S\'ersic fits, stellar masses, H$\alpha$-derived star formation rates, environment information, and group properties for all galaxies with r < 19.0 mag in two of our survey regions, and for all galaxies with r < 19.4 mag in a third region (72,225 objects in total). The database serving these data is available at http://www.gama-survey.org/.

We use the "Evolution and assembly of galaxies and their environments" (EAGLE) cosmological simulation to investigate the effect of baryons on the density profiles of rich galaxy clusters. We focus on EAGLE clusters with $M_{200}>10^{14}~M_\odot$ of which we have six examples. The central brightest cluster galaxies (BCGs) in the simulation have steep stellar density profiles, $\rho_*(r) \propto r^{-3}$. Stars dominate the mass density for $r<10~\rm{kpc}$, and, as a result, the total mass density profiles are steeper than the Navarro-Frenk-White (NFW) profile, in remarkable agreement with observations. The dark matter halo itself closely follows the NFW form at all resolved radii ($r\gtrsim3.0~\rm{kpc}$). The EAGLE BCGs have similar surface brightness and line-of-sight velocity dispersion profiles as the BCGs in the sample of Newman et al., which have the most detailed measurements currently available. After subtracting the contribution of the stars to the central density, Newman et al. infer significantly shallower slopes than the NFW value, in contradiction with the EAGLE results. We discuss possible reasons for this discrepancy, and conclude that an inconsistency between the kinematical model adopted by Newman et al. for their BCGs, which assumes isotropic stellar orbits, and the kinematical structure of the EAGLE BCGs, in which the orbital stellar anisotropy varies with radius and tends to be radially biased, could explain at least part of the discrepancy.

We investigate the internal structure and density profiles of halos of mass $10^{10}-10^{14}~M_\odot$ in the Evolution and Assembly of Galaxies and their Environment (EAGLE) simulations. These follow the formation of galaxies in a $\Lambda$CDM Universe and include a treatment of the baryon physics thought to be relevant. The EAGLE simulations reproduce the observed present-day galaxy stellar mass function, as well as many other properties of the galaxy population as a function of time. We find significant differences between the masses of halos in the EAGLE simulations and in simulations that follow only the dark matter component. Nevertheless, halos are well described by the Navarro-Frenk-White (NFW) density profile at radii larger than ~5% of the virial radius but, closer to the centre, the presence of stars can produce cuspier profiles. Central enhancements in the total mass profile are most important in halos of mass $10^{12}-10^{13}M_\odot$, where the stellar fraction peaks. Over the radial range where they are well resolved, the resulting galaxy rotation curves are in very good agreement with observational data for galaxies with stellar mass $M_*<5\times10^{10}M_\odot$. We present an empirical fitting function that describes the total mass profiles and show that its parameters are strongly correlated with halo mass.

FastSound is a galaxy redshift survey using the near-infrared Fiber Multi-Object Spectrograph (FMOS) mounted on the Subaru Telescope, targeting H$\alpha$ emitters at $z \sim 1.18$--$1.54$ down to the sensitivity limit of H$\alpha$ flux $\sim 2 \times 10^{-16} \ \rm erg \ cm^{-2} s^{-1}$. The primary goal of the survey is to detect redshift space distortions (RSD), to test General Relativity by measuring the growth rate of large scale structure and to constrain modified gravity models for the origin of the accelerated expansion of the universe. The target galaxies were selected based on photometric redshifts and H$\alpha$ flux estimates calculated by fitting spectral energy distribution (SED) models to the five optical magnitudes of the Canada France Hawaii Telescope Legacy Survey (CFHTLS) Wide catalog. The survey started in March 2012, and all the observations were completed in July 2014. In total, we achieved $121$ pointings of FMOS (each pointing has a $30$ arcmin diameter circular footprint) covering $20.6$ deg$^2$ by tiling the four fields of the CFHTLS Wide in a hexagonal pattern. Emission lines were detected from $\sim 4,000$ star forming galaxies by an automatic line detection algorithm applied to 2D spectral images. This is the first in a series of papers based on FastSound data, and we describe the details of the survey design, target selection, observations, data reduction, and emission line detections.

We investigate the evolution of galaxy masses and star formation rates in the Evolution and Assembly of Galaxies and their Environment (EAGLE) simulations. These comprise a suite of hydrodynamical simulations in a $\Lambda$CDM cosmogony with subgrid models for radiative cooling, star formation, stellar mass loss, and feedback from stars and accreting black holes. The subgrid feedback was calibrated to reproduce the observed present-day galaxy stellar mass function and galaxy sizes. Here we demonstrate that the simulations reproduce the observed growth of the stellar mass density to within 20 per cent. The simulation also tracks the observed evolution of the galaxy stellar mass function out to redshift z = 7, with differences comparable to the plausible uncertainties in the interpretation of the data. Just as with observed galaxies, the specific star formation rates of simulated galaxies are bimodal, with distinct star forming and passive sequences. The specific star formation rates of star forming galaxies are typically 0.2 to 0.4 dex lower than observed, but the evolution of the rates track the observations closely. The unprecedented level of agreement between simulation and data makes EAGLE a powerful resource to understand the physical processes that govern galaxy formation.

We calculate the colours and luminosities of redshift z = 0.1 galaxies from the EAGLE simulation suite using the GALAXEV population synthesis models. We take into account obscuration by dust in birth clouds and diffuse ISM using a two-component screen model, following the prescription of Charlot and Fall. We compare models in which the dust optical depth is constant to models where it depends on gas metallicity, gas fraction and orientation. The colours of EAGLE galaxies for the more sophisticated models are in broad agreement with those of observed galaxies. In particular, EAGLE produces a red sequence of passive galaxies and a blue cloud of star forming galaxies, with approximately the correct fraction of galaxies in each population and with g-r colours within 0.1 magnitudes of those observed. Luminosity functions from UV to NIR wavelengths differ from observations at a level comparable to systematic shifts resulting from a choice between Petrosian and Kron photometric apertures. Despite the generally good agreement there are clear discrepancies with observations. The blue cloud of EAGLE galaxies extends to somewhat higher luminosities than in the data, consistent with the modest underestimate of the passive fraction in massive EAGLE galaxies. There is also a moderate excess of bright blue galaxies compared to observations. The overall level of agreement with the observed colour distribution suggests that EAGLE galaxies at z = 0.1 have ages, metallicities and levels of obscuration that are comparable to those of observed galaxies.

We investigate the angular and kinematic distributions of satellite galaxies around a large sample of bright isolated primaries in the spectroscopic and photometric catalogues of the Sloan Digital Sky Survey (SDSS). We detect significant anisotropy in the spatial distribution of satellites. To test whether this anisotropy could be related to the rotating disks of satellites recently found by Ibata et al. in a sample of SDSS galaxies, we repeat and extend their analysis. Ibata et al. found an excess of satellites on opposite sides of their primaries having anticorrelated radial velocities. We find that this excess is sensitive to small changes in the sample selection criteria which can greatly reduce its significance. In addition, we find no evidence for correspondingly correlated velocities for satellites observed on the same side of their primaries, which would be expected for rotating disks of satellites. We conclude that the detection of coherent rotation in the satellite population in current observational samples is not robust. We compare our data to the $\Lambda$CDM Millennium simulations populated with galaxies according to the semi-analytic model of Guo et al. We find excellent agreement with the spatial distribution of satellites in the SDSS data and the lack of a strong signal from coherent rotation.

We present results from thirteen cosmological simulations that explore the parameter space of the "Evolution and Assembly of GaLaxies and their Environments" (EAGLE) simulation project. Four of the simulations follow the evolution of a periodic cube L = 50 cMpc on a side, and each employs a different subgrid model of the energetic feedback associated with star formation. The relevant parameters were adjusted so that the simulations each reproduce the observed galaxy stellar mass function at z = 0.1. Three of the simulations fail to form disc galaxies as extended as observed, and we show analytically that this is a consequence of numerical radiative losses that reduce the efficiency of stellar feedback in high-density gas. Such losses are greatly reduced in the fourth simulation - the EAGLE reference model - by injecting more energy in higher density gas. This model produces galaxies with the observed size distribution, and also reproduces many galaxy scaling relations. In the remaining nine simulations, a single parameter or process of the reference model was varied at a time. We find that the properties of galaxies with stellar mass <~ M* (the "knee" of the galaxy stellar mass function) are largely governed by feedback associated with star formation, while those of more massive galaxies are also controlled by feedback from accretion onto their central black holes. Both processes must be efficient in order to reproduce the observed galaxy population. In general, simulations that have been calibrated to reproduce the low-redshift galaxy stellar mass function will still not form realistic galaxies, but the additional requirement that galaxy sizes be acceptable leads to agreement with a large range of observables.

We apply a semi-analytic galaxy formation model to two high resolution cosmological N-body simulations to investigate analogues of the Milky Way system. We select these according to observed properties of the Milky Way rather than by halo mass as in most previous work. For disk-dominated central galaxies with stellar mass (5--7) x 10d10Msun, the median host halo mass is 1.4 x 10d12Msun, with 1 sigma dispersion in the range [0.86, 3.1] x 10d12Msun, consistent with dynamical measurements of the Milky Way halo mass. For any given halo mass, the probability of hosting a Milky Way system is low, with a maximum of ~20% in haloes of mass ~10d12Msun. The model reproduces the V-band luminosity function and radial profile of the bright (MV < -9) Milky Way satellites. Galaxy formation in low mass haloes is found to be highly stochastic, resulting in an extremely large scatter in the relation between MV (or stellar mass) for satellites and the depth of the subhalo potential well in which they live, as measured by the maximum of the rotation curve, Vmax. We conclude that the "too big to fail" problem is an artifact of selecting satellites in N-body simulations according to subhalo properties: in 10% of cases we find that three or fewer of the brightest (or most massive) satellites have Vmax > 30 km/s. Our model predicts that around half of the dark matter subhaloes with Vmax > 20 km/s host satellites fainter than MV = -9 and so may be missing from existing surveys.

We establish a controlled comparison between the properties of galactic stellar halos obtained with hydrodynamical simulations and with `particle tagging'. Tagging is a fast way to obtain stellar population dynamics: instead of tracking gas and star formation, it `paints' stars directly onto a suitably defined subset of dark matter particles in a collisionless, dark-matter-only simulation.Our study shows that there are conditions under which particle tagging generates good fits to the hydrodynamical stellar density profiles of a central Milky-Way-like galaxy and its most prominent substructure. Phase-space diffusion processes are crucial to reshaping the distribution of stars in infalling spheroidal systems and hence the final stellar halo. We conclude that the success of any particular tagging scheme hinges on this diffusion being taken into account, at a minimum by making use of `live' tagging schemes, in which particles are regularly tagged throughout the evolution of a galaxy.

We present distributions of orbital parameters of infalling satellites of $\Lambda$CDM haloes in the mass range $10^{12}-10^{14}$M$_\odot$, which represent the initial conditions for the subsequent evolution of substructures within the host halo. We use merger trees constructed in a high resolution cosmological N-body simulation to trace satellite haloes, and identify the time of infall. We find signficant trends in the distribution of orbital parameters with both the host halo mass and the ratio of satellite-to-host halo masses. For all host halo masses, satellites whose infall mass is a larger fraction of the host halo mass have more eccentric, radially biased orbits. At fixed satellite-to-host halo mass ratio, high mass haloes are biased towards accreting satellites on slightly more radial orbits. To charactise the orbital distributions fully requires fitting the correlated bivariate distribution of two chosen orbital parameters (e.g. radial and tangential velocity or energy and angular momentum). We provide simple fits to one choice of the bivariate distributions, which when transformed faithfully, captures the behaviour of any of the projected one-dimensional distributions.

Discrepancies between the observed and model-predicted radio flux ratios are seen in a number of quadruply-lensed quasars. The most favored interpretation of these anomalies is that CDM substructures present in lensing galaxies perturb the lens potentials and alter image magnifications and thus flux ratios. So far no consensus has emerged regarding whether or not the predicted CDM substructure abundance fully accounts for the lensing flux anomaly observations. Accurate modeling relies on a realistic lens sample in terms of both the lens environment and internal structures and substructures. In this paper we construct samples of generalised and specific lens potentials, to which we add (rescaled) subhalo populations from the galaxy-scale Aquarius and the cluster-scale Phoenix simulation suites. We further investigate the lensing effects from subhalos of masses several orders of magnitude below the simulation resolution limit. The resulting flux ratio distributions are compared to the currently best available sample of radio lenses. The observed anomalies in B0128+437, B0712+472 and B1555+375 are more likely to be caused by propagation effects or oversimplified lens modeling, signs of which are already seen in the data. Among the quadruple systems that have closely located image triplets/pairs, the anomalous flux ratios of MG0414+0534 can be reproduced by adding CDM subhalos to its macroscopic lens potential, with a probability of 5%-20%; for B0712+472, B1422+231, B1555+375 and B2045+265, these probabilities are only of a few percent. We hence find that CDM substructures are unlikely to be the whole reason for radio flux anomalies. We discuss other possible effects that might also be at work.

We investigate the characteristics and the time evolution of the cosmic web from redshift, z=2, to present time, within the framework of the NEXUS+ algorithm. This necessitates the introduction of new analysis tools optimally suited to describe the very intricate and hierarchical pattern that is the cosmic web. In particular, we characterize filaments (walls) in terms of their linear (surface) mass density. This is very good in capturing the evolution of these structures. At early times the cosmos is dominated by tenuous filaments and sheets, which, during subsequent evolution, merge together, such that the present day web is dominated by fewer, but much more massive, structures. We also show that voids are more naturally described in terms of their boundaries and not their centres. We illustrate this for void density profiles, which, when expressed as a function of the distance from void boundary, show a universal profile in good qualitative agreement with the theoretical shell-crossing framework of expanding underdense regions.

The "Lambda Cold Dark Matter" (LCDM) model of cosmic structure formation is eminently falsifiable: once its parameters are fixed on large scales, it becomes testable in the nearby Universe. Observations within our Local Group of galaxies, including the satellite populations of the Milky Way and Andromeda, appear to contradict LCDM predictions: there are far fewer satellite galaxies than dark matter halos (the "missing satellites" problem), galaxies seem to avoid the largest substructures (the "too big to fail" problem), and the brightest satellites appear to orbit their host galaxies on a thin plane (the "planes of satellites" problem). We present results from the first hydrodynamic simulations of the Local Group that match the observed abundance of galaxies. We find that when baryonic and dark matter are followed simultaneously in the context of a realistic galaxy formation model, all three "problems" are resolved within the LCDM paradigm.

We present a maximum-likelihood weak lensing analysis of the mass distribution in optically selected spectroscopic Galaxy Groups (G$^3$Cv5) in the Galaxy And Mass Assembly (GAMA) survey, using background Sloan Digital Sky Survey (SDSS) photometric galaxies. The scaling of halo mass, $M_h$, with various group observables is investigated. Our main results are: 1) the measured relations of halo mass with group luminosity, virial volume and central galaxy stellar mass, $M_\star$, agree very well with predictions from mock group catalogues constructed from a GALFORM semi-analytical galaxy formation model implemented in the Millennium $\Lambda$CDM N-body simulation; 2) the measured relations of halo mass with velocity dispersion and projected half-abundance radius show weak tension with mock predictions, hinting at problems in the mock galaxy dynamics and their small scale distribution; 3) the median $M_h|M_\star$ measured from weak lensing depends more sensitively on the lognormal dispersion in $M_\star$ at fixed $M_h$ than it does on the median $M_\star|M_h$. Our measurements suggest an intrinsic dispersion of $\sigma_{\log(M_\star)}\sim 0.15$; 4) Comparing our mass estimates with those in the catalogue, we find that the G$^3$Cv5 mass can give biased results when used to select subsets of the group sample. Of the various new halo mass estimators that we calibrate using our weak lensing measurements, group luminosity is the best single-proxy estimator of group mass.

We present a new technique for creating mock catalogues of the individual stars that make up the accreted component of stellar haloes in cosmological simulations and show how the catalogues can be used to test and interpret observational data. The catalogues are constructed from a combination of methods. A semi-analytic galaxy formation model is used to calculate the star formation history in haloes in an N-body simulation and dark matter particles are tagged with this stellar mass. The tags are converted into individual stars using a stellar population synthesis model to obtain the number density and evolutionary stage of the stars, together with a phase-space sampling method that distributes the stars while ensuring that the phase-space structure of the original N-body simulation is maintained. A set of catalogues based on the $\Lambda$CDM Aquarius simulations of Milky Way mass haloes have been created and made publicly available on a website. Two example applications are discussed that demonstrate the power and flexibility of the mock catalogues. We show how the rich stellar substructure that survives in the stellar halo precludes a simple measurement of its density profile and demonstrate explicitly how pencil-beam surveys can return almost any value for the slope of the profile. We also show that localized variations in the abundance of particular types of stars, a signature of differences in the composition of stellar populations, allow streams to be easily identified.

We study the substructure population of Milky Way (MW)-mass halos in the $\Lambda$CDM cosmology using a novel procedure to extrapolate subhalo number statistics beyond the resolution limit of N-body simulations. The technique recovers the mean and the variance of the subhalo abundance, but not its spatial distribution. It extends the dynamic range over which precise statistical predictions can be made by the equivalent of performing a simulation with 50 times higher resolution, at no additional computational cost. We apply this technique to MW-mass halos, but it can easily be applied to halos of any mass. We find up to $20\%$ more substructures in MW-mass halos than found in previous studies. Our analysis lowers the mass of the MW halo required to accommodate the observation that the MW has only three satellites with a maximum circular velocity $V_{max}\ge30 km/s$ in the $\Lambda$CDM cosmology. The probability of having a subhalo population similar to that in the MW is $20\%$ for a virial mass, $M_{200}=1\times10^{12} M_\odot$ and practically zero for halos more massive than $M_{200}=2\times10^{12} M_\odot$.

We use the distribution of maximum circular velocities, $V_{max}$, of satellites in the Milky Way (MW) to constrain the virial mass, $M_{200}$, of the Galactic halo under an assumed prior of a $\Lambda$CDM universe. This is done by analysing the subhalo populations of a large sample of halos found in the Millennium II cosmological simulation. The observation that the MW has at most three subhalos with $V_{max}\ge30 km/s$ requires a halo mass $M_{200}\le1.4\times10^{12} M_\odot$, while the existence of the Magellanic Clouds (assumed to have $V_{max}\ge60 km/s$) requires $M_{200}\ge1.0\times10^{12} M_\odot$. The first of these conditions is necessary to avoid the "too-big-to-fail" problem highlighted by Boylan-Kolchin et al., while the second stems from the observation that massive satellites like the Magellanic Clouds are rare. When combining both requirements, we find that the MW halo mass must lie in the range $0.25 \le M_{200}/(10^{12} M_\odot) \le 1.4$ at $90\%$ confidence. The gap in the abundance of Galactic satellites between $30 km/s\le V_{max} \le 60 km/s$ places our galaxy in the tail of the expected satellite distribution.

Lensing flux-ratio anomalies are most likely caused by gravitational lensing by small-scale dark matter structures. These anomalies offer the prospect of testing a fundamental prediction of the cold dark matter (CDM) cosmological model: the existence of numerous substructures that are too small to host visible galaxies. In two previous studies we found that the number of subhalos in the six high-resolution simulations of CDM galactic halos of the Aquarius project is not sufficient to account for the observed frequency of flux ratio anomalies seen in selected quasars from the CLASS survey. These studies were limited by the small number of halos used, their narrow range of masses (1-2E12 solar masses) and the small range of lens ellipticities considered. We address these shortcomings by investigating the lensing properties of a large sample of halos with a wide range of masses in two sets of high resolution simulations of cosmological volumes and comparing them to a currently best available sample of radio quasars. We find that, as expected, substructures do not change the flux-ratio probability distribution of image pairs and triples with large separations, but they have a significant effect on the distribution at small separations. For such systems, CDM substructures can account for a substantial fraction of the observed flux-ratio anomalies. For large close-pair separation systems, the discrepancies existing between the observed flux ratios and predictions from smooth halo models are attributed to simplifications inherent in these models which do not take account of fine details in the lens mass distributions.

Recent work has suggested that the amplitude of the size mass relation of massive early type galaxies evolves with redshift. Here we use a semi-analytical galaxy formation model to study the size evolution of massive early type galaxies. We find this model is able to reproduce the amplitude of present day amplitude and slope of the relation between size and stellar mass for these galaxies, as well as its evolution. The amplitude of this relation reflects the typical compactness of dark halos at the time when most of the stars are formed. This link between size and star formation epoch is propagated in galaxy mergers. Mergers of high or moderate mass ratio (less than 1:3) become increasingly important with increasing present day stellar mass for galaxies more massive than $10^{11.4}M_{\odot}$. At lower masses, low mass ratio mergers play a more important role. In situ star formation contribute more to the size growth than it does to stellar mass growth. We also find that, for ETGs identified at $z=2$, minor mergers dominate subsequent growth both for stellar mass and in size, consistent with earlier theoretical results.

We introduce the Virgo Consortium's EAGLE project, a suite of hydrodynamical simulations that follow the formation of galaxies and black holes in representative volumes. We discuss the limitations of such simulations in light of their finite resolution and poorly constrained subgrid physics, and how these affect their predictive power. One major improvement is our treatment of feedback from massive stars and AGN in which thermal energy is injected into the gas without the need to turn off cooling or hydrodynamical forces, allowing winds to develop without predetermined speed or mass loading factors. Because the feedback efficiencies cannot be predicted from first principles, we calibrate them to the z~0 galaxy stellar mass function and the amplitude of the galaxy-central black hole mass relation, also taking galaxy sizes into account. The observed galaxy mass function is reproduced to $\lesssim 0.2$ dex over the full mass range, $10^8 < M_*/M_\odot \lesssim 10^{11}$, a level of agreement close to that attained by semi-analytic models, and unprecedented for hydrodynamical simulations. We compare our results to a representative set of low-redshift observables not considered in the calibration, and find good agreement with the observed galaxy specific star formation rates, passive fractions, Tully-Fisher relation, total stellar luminosities of galaxy clusters, and column density distributions of intergalactic CIV and OVI. While the mass-metallicity relations for gas and stars are consistent with observations for $M_* \gtrsim 10^9 M_\odot$, they are insufficiently steep at lower masses. The gas fractions and temperatures are too high for clusters of galaxies, but for groups these discrepancies can be resolved by adopting a higher heating temperature in the subgrid prescription for AGN feedback. EAGLE constitutes a valuable new resource for studies of galaxy formation.

We develop a new method to account for the finite lifetimes of stars and trace individual abundances within a semi-analytic model of galaxy formation. At variance with previous methods, based on the storage of the (binned) past star formation history of model galaxies, our method projects the information about the metals produced by each simple stellar population (SSP) in the future. Using this approach, an accurate accounting of the timings and properties of the individual SSPs composing model galaxies is possible. We analyse the dependence of our chemical model on various ingredients, and apply it to six simulated haloes of roughly Milky Way mass and with no massive close neighbour at z=0. For all models considered, the [Fe/H] distributions of the stars in the disc component are in good agreement with Milky Way data, while for the spheroid component (whose formation we model only through mergers) these are offset low with respect to observational measurements for the Milky Way bulge. This is a consequence of narrow star formation histories, with relatively low rates of star formation. The slow recycling of gas and energy from supernovae in our chemical model has important consequences on the predicted star formation rates, which are systematically lower than the corresponding rates in the same physical model but with an instantaneous recycling approximation. The halo that resembles most our Galaxy in terms of its global properties also reproduces the observed relation between the average metallicity and luminosity of the Milky Way satellites, albeit with a slightly steeper slope.

We use the Herschel ATLAS (H-ATLAS) Phase I data to study the conditional luminosity function of far-IR (250 um) selected galaxies in optically-selected galaxy groups from the GAMA spectroscopic survey, as well as environmental effects on the far-IR-to-optical colour. We applied two methods, which gave consistent results for the far-IR conditional luminosity functions. The direct matching method matches H-ATLAS sources to GAMA/SDSS galaxies, then links the optical counterparts to GAMA groups. The stacking method counts the number of far-IR sources within the projected radii of GAMA groups, subtracting the local background. We investigated the dependence of the far-IR (250 um) luminosity function on group mass in the range 10d12 < M_h < 10d14Msun/h and on redshift in the range 0 < z < 0.4, using a sample of 3000 groups containing H-ATLAS sources with GAMA redshifts over an area of 126 sq. deg. We find that the characteristic 250 um luminosity, L^*(250), increases with group mass up to M_h ~ 10d13Msun/h, but is roughly constant above this, while it increases with redshift at high group masses, but less so at low masses. We also find that the group far-IR luminosity-to-mass ratio L(250)/M_h increases with redshift and is higher in low-mass groups. We estimate that around 70% of the 250 um luminosity density in the local universe is contributed by groups with M_h > 10d12Msun/h. We also find that the far-IR-to-optical colours of H-ATLAS galaxies are independent of group mass over the range 10d12 < M_h <10d14Msun/h in the local universe. We also compare our observational results with recent semi-analytical models, and find that none of these galaxy formation model can reproduce the conditional far-IR luminosity functions of galaxy groups.

Since reionization prevents star formation in most halos below 3 x 10^9 solar masses, dwarf galaxies only populate a fraction of existing dark matter halos. We use hydrodynamic cosmological simulations of the Local Group to study the discriminating factors for galaxy formation in the early Universe and connect them to the present-day properties of galaxies and halos. A combination of selection effects related to reionization, and the subsequent evolution of halos in different environments, introduces strong biases between the population of halos that host dwarf galaxies, and the total halo population. Halos that host galaxies formed earlier and are more concentrated. In addition, halos more affected by tidal stripping are more likely to host a galaxy for a given mass or maximum circular velocity, vmax, today. Consequently, satellite halos are populated more frequently than field halos, and satellite halos of 10^8 - 10^9 solar masses or vmax of 12 - 20 km/s, similar to the Local Group dwarf spheroidals, have experienced a greater than average reduction in both mass and vmax after infall. They are on closer, more radial orbits with higher infall velocities and earlier infall times. Together, these effects make dwarf galaxies highly biased tracers of the underlying dark matter distribution.

The velocity field of dark matter and galaxies reflects the continued action of gravity throughout cosmic history. We show that the low-order moments of the pairwise velocity distribution, $v_{12}$, are a powerful diagnostic of the laws of gravity on cosmological scales. In particular, the projected line-of-sight galaxy pairwise velocity dispersion, $\sigma_{12}(r)$, is very sensitive to the presence of modified gravity. Using a set of high-resolution N-body simulations we compute the pairwise velocity distribution and its projected line-of-sight dispersion for a class of modified gravity theories: the chameleon \fR gravity and Galileon gravity (cubic and quartic). The velocities of dark matter halos with a wide range of masses would exhibit deviations from General Relativity at the $(5-10)\sigma$ level. We examine strategies for detecting these deviations in galaxy redshift and peculiar velocity surveys. If detected, this signature would be a "smoking gun" for modified gravity.

The cosmic web is the largest scale manifestation of the anisotropic gravitational collapse of matter. It represents the transitional stage between linear and non-linear structures and contains easily accessible information about the early phases of structure formation processes. Here we investigate the characteristics and the time evolution of morphological components since. Our analysis involves the application of the NEXUS Multiscale Morphology Filter (MMF) technique, predominantly its NEXUS+ version, to high resolution and large volume cosmological simulations. We quantify the cosmic web components in terms of their mass and volume content, their density distribution and halo populations. We employ new analysis techniques to determine the spatial extent of filaments and sheets, like their total length and local width. This analysis identifies cluster and filaments as the most prominent components of the web. In contrast, while voids and sheets take most of the volume, they correspond to underdense environments and are devoid of group-sized and more massive haloes. At early times the cosmos is dominated by tenuous filaments and sheets, which, during subsequent evolution, merge together, such that the present day web is dominated by fewer, but much more massive, structures. The analysis of the mass transport between environments clearly shows how matter flows from voids into walls, and then via filaments into cluster regions, which form the nodes of the cosmic web. We also study the properties of individual filamentary branches, to find long, almost straight, filaments extending to distances larger than 100Mpc/h. These constitute the bridges between massive clusters, which seem to form along approximatively straight lines.

Well-motivated elementary particle candidates for the dark matter, such as the sterile neutrino, behave as warm dark matter (WDM).For particle masses of order a keV, free streaming produces a cutoff in the linear fluctuation power spectrum at a scale corresponding to dwarf galaxies. We investigate the abundance and structure of WDM haloes and subhaloes on these scales using high resolution cosmological N-body simulations of galactic haloes of mass similar to the Milky Way's. On scales larger than the free-streaming cutoff, the initial conditions have the same power spectrum and phases as one of the cold dark matter (CDM) haloes previously simulated by Springel et al as part of the Virgo consortium Aquarius project. We have simulated four haloes with WDM particle masses in the range 1.4-2.3keV and, for one case, we have carried out further simulations at varying resolution. N-body simulations in which the power spectrum cutoff is resolved are known to undergo artificial fragmentation in filaments producing spurious clumps which, for small masses (<10^7Msun in our case) outnumber genuine haloes. We have developed a robust algorithm to identify these spurious objects and remove them from our halo catalogues. We find that the WDM subhalo mass function is suppressed by well over an order magnitude relative to the CDM case for masses <10^9Msun. Requiring that there should be at least as many subhaloes as there are observed satellites in the Milky Way leads to a conservative lower limit to the (thermal equivalent) WDM particle mass of ~1.5\rmn{keV}. WDM haloes and subhaloes have cuspy density distributions that are well described by NFW or Einasto profiles. Their central densities are lower for lower WDM particle masses and none of the models we have considered suffer from the "too big to fail" problem recently highlighted by Boylan-Kolchin et al.

The relation between galaxies and dark matter halos is of vital importance for evaluating theoretical predictions of structure formation and galaxy formation physics. We show that the widely used method of abundance matching based on dark matter only simulations fails at the low mass end because two of its underlying assumptions are broken: only a small fraction of low mass (below 10^9.5 solar masses) halos host a visible galaxy, and halos grow at a lower rate due to the effect of baryons. In this regime, reliance on dark matter only simulations for abundance matching is neither accurate nor self-consistent. We find that the reported discrepancy between observational estimates of the halo masses of dwarf galaxies and the values predicted by abundance matching does not point to a failure of LCDM, but simply to a failure to account for baryonic effects. Our results also imply that the Local Group contains only a few hundred observable galaxies in contrast with the thousands of faint dwarfs that abundance matching would suggest. We show how relations derived from abundance matching can be corrected, so that they can be used self-consistently to calibrate models of galaxy formation.

We present a new algorithm which groups the subhaloes found in cosmological N- body simulations by structure finders such as SUBFIND into dark matter haloes whose formation histories are strictly hierarchical. One advantage of these `Dhaloes' over the commonly used friends-of-friends (FoF) haloes is that they retain their individual identity in cases when FoF haloes are artificially merged by tenuous bridges of particles or by an overlap of their outer diffuse haloes. Dhaloes are thus well-suited for modelling galaxy formation and their merger trees form the basis of the Durham semi-analytic galaxy formation model, GALFORM. Applying the Dhalo construction to the LambdaCDM Millennium-2 simulation we find that approximately 90% of Dhaloes have a one-to-one, bijective match with a corresponding FoF halo. The remaining 10% are typically secondary components of large FoF haloes. Although the mass functions of both types of haloes are similar, the mass of Dhaloes correlates much more tightly with the virial mass, M200, than FoF haloes. Approximately 80% of FoF and bijective and non-bijective Dhaloes are relaxed according to standard criteria. For these relaxed haloes all three types have similar concentration-M200 relations and, at fixed mass, the concentration distributions are described accurately by log-normal distributions.

We measure the average temperature decrement on the cosmic microwave background (CMB) produced by voids selected in the SDSS DR7 spectroscopic redshift galaxy catalog, spanning redshifts $0<z<0.44$. We find an imprint of amplitude between 2.6 and 2.9$\mu K$ as viewed through a compensated top-hat filter scaled to the radius of each void; we assess the statistical significance of the imprint at ~2$\sigma$. We make crucial use of $N$-body simulations to calibrate our analysis. As expected, we find that large voids produce cold spots on the CMB through the Integrated Sachs-Wolfe (ISW) effect. However, we also find that small voids in the halo density field produce hot spots, because they reside in contracting, larger-scale overdense regions. This is an important effect to consider when stacking CMB imprints from voids of different radius. We have found that the same filter radius that gives the largest ISW signal in simulations also yields close to the largest detected signal in the observations. However, although it is low in significance, our measured signal is much higher-amplitude than expected from ISW in the concordance $\Lambda$CDM universe. The discrepancy is also at the ~2$\sigma$ level. We have demonstrated that our result is robust against the varying of thresholds over a wide range.

We study the late-time Integrated Sachs-Wolfe (ISW) effect in $f(R)$ gravity using N-body simulations. In the $f(R)$ model under study, the linear growth rate is larger than that in general relativity (GR). This slows down the decay of the cosmic potential and induces a smaller ISW effect on large scales. Therefore, the $\dot\Phi$ (time derivative of the potential) power spectrum at $k<0.1h$/Mpc is suppressed relative to that in GR. In the non-linear regime, relatively rapid structure formation in $f(R)$ gravity boosts the non-linear ISW effect relative to GR, and the $\dot\Phi$ power spectrum at $k>0.1h$/Mpc is increased (100$\%$ greater on small scales at $z=0$). We explore the detectability of the ISW signal via stacking supercluster and supervoids. The differences in the corresponding ISW cold or hot spots are $\sim 20\%$ for structures of $\sim 100$Mpc/$h$. Such differences are greater for smaller structures, but the amplitude of the signal is lower. The high amplitude of ISW signal detected by Granett et al. can not explained in the $f(R)$ model. We find relatively big differences between $f(R)$ and GR in the transverse bulk motion of matter, and discuss its detectability via the relative frequency shifts of photons from multiple lensed images.

We use the Millennium Simulation series to study the relation between the accretion history (MAH) and mass profile of cold dark matter halos. We find that the mean density within the scale radius, r_{-2} (where the halo density profile has isothermal slope), is directly proportional to the critical density of the Universe at the time when the main progenitor's virial mass equals the mass enclosed within r_{-2}. Scaled to these characteristic values of mass and density, the mean MAH, expressed in terms of the critical density of the Universe, M(\rho_{crit}(z)), resembles that of the enclosed density profile, M(<\rho >), at z=0. Both follow closely the NFW profile, suggesting that the similarity of halo mass profiles originates from the mass-independence of halo MAHs. Support for this interpretation is provided by outlier halos whose accretion histories deviate from the NFW shape; their mass profiles show correlated deviations from NFW and are better approximated by Einasto profiles. Fitting both M(<\rho >) and M(\rho_{crit}) with either NFW or Einasto profiles yield concentration and shape parameters that are correlated, confirming and extending earlier work linking the concentration of a halo with its accretion history. These correlations also confirm that halo structure is insensitive to initial conditions: only halos whose accretion histories differ greatly from the NFW shape show noticeable deviations from NFW in their mass profiles. As a result, the NFW profile provides acceptable fits to hot dark matter halos, which do not form hierarchically, and for fluctuation power spectra other than CDM. Our findings, however, predict a subtle but systematic dependence of mass profile shape on accretion history which, if confirmed, would provide strong support for the link between accretion history and halo structure we propose here.

We use the Millennium Simulation series to investigate the mass and redshift dependence of the concentration of equilibrium cold dark matter (CDM) halos. We extend earlier work on the relation between halo mass profiles and assembly histories to show how the latter may be used to predict concentrations for halos of all masses and at any redshift. Our results clarify the link between concentration and the ``collapse redshift'' of a halo as well as why concentration depends on mass and redshift solely through the dimensionless ``peak height'' mass parameter, $\nu(M,z)=\delta_{\rm crit}(z)/\sigma(M,z)$. We combine these results with analytic mass accretion histories to extrapolate the $c(M,z)$ relations to mass regimes difficult to reach through direct simulation. Our model predicts that, at given $z$, $c(M)$ should deviate systematically from a simple power law at high masses, where concentrations approach a constant value, and at low masses, where concentrations are substantially lower than expected from extrapolating published empirical fits. This correction may reduce the expected self-annihilation boost factor from substructure by about one order of magnitude. The model also reproduces the $c(M,z)$ dependence on cosmological parameters reported in earlier work, and thus provides a simple and robust account of the relation between cosmology and the mass-concentration-redshift relation of CDM halos.

We use the very high resolution, fully cosmological simulations from the Aquarius project, coupled to a semi-analytical model of galaxy formation, to study the phase-space distribution of halo stars in "solar neighbourhood"-like volumes. We find that this distribution is very rich in substructure in the form of stellar streams for all five stellar haloes we have analysed. These streams can be easily identified in velocity space, as well as in spaces of pseudo-conserved quantities such as E vs. Lz. In our best-resolved local volumes, the number of identified streams ranges from ~ 300 to 600, in very good agreement with previous analytical predictions, even in the presence of chaotic mixing. The fraction of particles linked to (massive) stellar streams in these volumes can be as large as 84%. The number of identified streams is found to decrease as a power-law with galactocentric radius. We show that the strongest limitation to the quantification of substructure in our poorest-resolved local volumes is particle resolution rather than strong diffusion due to chaotic mixing.

Well-motivated particle physics theories predict the existence of particles (such as sterile neutrinos) which acquire non-negligible thermal velocities in the early universe. These particles could behave as warm dark matter (WDM) and generate a small-scale cutoff in the linear density power spectrum which scales approximately inversely with the particle mass. If this mass is of order a keV, the cutoff occurs on the scale of dwarf galaxies. Thus, in WDM models the abundance of small galaxies, such as the satellites that orbit in the halo of the Milky Way, depends on the mass of the warm particle. The abundance also scales with the mass of the host galactic halo. We use the \galform semi-analytic model of galaxy formation to calculate the properties of galaxies in universes in which the dark matter is warm. Using this method, we can compare the predicted satellite luminosity functions to the observed data for the Milky Way dwarf spheroidals, and determine a lower bound on the thermally produced WDM particle mass. This depends strongly on the value of the Milky Way halo mass and, to some extent, on the baryonic physics assumed; we examine both of these dependencies. For our fiducial model we find that for a particle mass of 3.3 keV (the 2$\sigma$ lower limit found by Viel et al. from a recent analysis of the Lyman-$\alpha$ forest) the Milky Way halo mass is required to be $> 1.4 \times 10^{12}$ \msun. For this same fiducial model, we also find that all WDM particle masses are ruled out (at 95% confidence) if the halo of the Milky Way has a mass smaller than $1.1 \times 10^{12}$ \msun, while if the mass of the Galactic halo is greater than 1.8 $\times 10^{12}$ \msun, only WDM particle masses larger than 2 keV are allowed.

We present a scheme to extend the halo mass resolution of N-body simulations of the hierarchical clustering of dark matter. The method uses the density field of the simulation to predict the number of sub-resolution dark matter haloes expected in different regions. The technique requires as input the abundance of haloes of a given mass and their average clustering, as expressed through the linear and higher order bias factors. These quantities can be computed analytically or, more accurately, derived from a higher resolution simulation as done here. Our method can recover the abundance and clustering in real- and redshift-space of haloes with mass below $\sim 7.5 \times 10^{13}h^{-1}M_{\odot}$ at $z=0$ to better than 10%. We demonstrate the technique by applying it to an ensemble of 50 low resolution, large-volume $N$-body simulations to compute the correlation function and covariance matrix of luminous red galaxies (LRGs). The limited resolution of the original simulations results in them resolving just two thirds of the LRG population. We extend the resolution of the simulations by a factor of 30 in halo mass in order to recover all LRGs. With existing simulations it is possible to generate a halo catalogue equivalent to that which would be obtained from a $N$-body simulation using more than 20 trillion particles; a direct simulation of this size is likely to remain unachievable for many years. Using our method it is now feasible to build the large numbers of high-resolution large volume mock galaxy catalogues required to compute the covariance matrices necessary to analyse upcoming galaxy surveys designed to probe dark energy.

We apply a recently developed scaling technique to the Millennium-XXL, one of the largest cosmological N-body simulations carried out to date 3x10^11 particles within a cube of volume ~70Gpc^3). This allows us to investigate the cosmological parameter dependence of the mass and evolution of haloes in the extreme high-mass tail of the z=6 distribution. We assume these objects to be likely hosts for the population of rare but ultraluminous high-redshift quasars discovered by the Sloan Digital Sky Survey. Haloes with a similar abundance to these quasars have a median mass of 9x10^12 Msun in the currently preferred cosmology, but do not evolve into equally extreme objects at z=0. Rather, their descendants span the full range conventionally assigned to present-day clusters, 6x10^13 to 2.5x10^15 Msun for this same cosmology. The masses both at z=6 and at z=0 shift up or down by factors exceeding two if cosmological parameters are pushed to the boundaries of the range discussed in published interpretations of data from the WMAP satellite. The main factor determining the future growth of a high-mass z=6 halo is the mean overdensity of its environment on scales of 7 to 14 Mpc, and descendant masses can be predicted 6 to 8 times more accurately if this density is known than if it is not. All these features are not unique to extreme high-z haloes, but are generic to hierarchical growth. Finally, we find that extreme haloes at z=6 typically acquired about half of their total mass in the preceding 100 Myr, implying very large recent accretion rates which may be related to the large black hole masses and high luminosities of the SDSS quasars.

We present a very large high-resolution cosmological N-body simulation, the Millennium-XXL or MXXL, which uses 303 billion particles to represent the formation of dark matter structures throughout a 4.1Gpc box in a LambdaCDM cosmology. We create sky maps and identify large samples of galaxy clusters using surrogates for four different observables: richness estimated from galaxy surveys, X-ray luminosity, integrated Sunyaev-Zeldovich signal, and lensing mass. The unprecedented combination of volume and resolution allows us to explore in detail how these observables scale with each other and with cluster mass. The scatter correlates between different mass-observable relations because of common sensitivities to the internal structure, orientation and environment of clusters, as well as to line-of-sight superposition of uncorrelated structure. We show that this can account for the apparent discrepancies uncovered recently between the mean thermal SZ signals measured for optically and X-ray selected clusters by stacking data from the Planck satellite. Related systematics can also affect inferences from extreme clusters detected at high redshift. Our results illustrate that cosmological conclusions from galaxy cluster surveys depend critically on proper modelling, not only of the relevant physics, but also of the full distribution of the observables and of the selection biases induced by cluster identification procedures.

We have combined the semi-analytic galaxy formation model of Guo et al. (2011) with the particle-tagging technique of Cooper et al. (2010) to predict galaxy surface brightness profiles in a representative sample of ~1900 massive dark matter haloes (10^12--10^14 M_sol) from the Millennium II Lambda-CDM N-body simulation. Here we present our method and basic results focusing on the outer regions of galaxies, consisting of stars accreted in mergers. These simulations cover scales from the stellar haloes of Milky Way-like galaxies to the 'cD envelopes' of groups and clusters, and resolve low surface brightness substructure such as tidal streams. We find that the surface density of accreted stellar mass around the central galaxies of dark matter haloes is well described by a Sersic profile, the radial scale and amplitude of which vary systematically with halo mass (M_200). The total stellar mass surface density profile breaks at the radius where accreted stars start to dominate over stars formed in the galaxy itself. This break disappears with increasing M_200 because accreted stars contribute more of the total mass of galaxies, and is less distinct when the same galaxies are averaged in bins of stellar mass, because of scatter in the relation between M_star and M_200. To test our model we have derived average stellar mass surface density profiles for massive galaxies at z~0.08 by stacking SDSS images. Our model agrees well with these stacked profiles and with other data from the literature, and makes predictions that can be more rigorously tested by future surveys that extend the analysis of the outer structure of galaxies to fainter isophotes. We conclude that it is likely that the outer structure of the spheroidal components of galaxies is largely determined by collisionless merging during their hierarchical assembly

We constrain the distributions of projected radial alignment angles of satellite galaxy shapes within the Galaxy And Mass Assembly survey group catalogue. We identify the galaxy groups using spectroscopic redshifts and measure galaxy projected ellipticities from Sloan Digital Sky Survey imaging. With a sample of 3,850 groups with 13,655 satellite galaxies with high quality shape measurements, we find a less than 2-sigma signal of radial alignments in the mean projected ellipticity components and the projected position angle when using galaxy shape estimates optimized for weak lensing measurements. Our radial alignment measurement increases to greater than 3-sigma significance relative to the expectation for no alignments if we use 2-D S\'ersic model fits to define galaxy orientations. Our weak measurement of radial alignments is in conflict with predictions from dark matter N-body simulations, which we interpret as evidence for large mis-alignments of baryons and dark matter in group and cluster satellites. Within our uncertainties, that are dominated by our small sample size, we find only weak and marginally significant trends of the radial alignment angle distributions on projected distance from the group centre, host halo mass, and redshift that could be consistent with a tidal torquing mechanism for radial alignments. Using our lensing optimized shape estimators, we estimate that intrinsic alignments of galaxy group members may contribute a systematic error to the mean differential projected surface mass density of groups inferred from weak lensing observations by -1 +/- 20% at scales around 300 kpc/h from the group centre assuming a photometric redshift r.m.s. error of 10%, and given our group sample with median redshift of 0.17 and median virial masses ~10^{13} h^{-1}M_{sun}.

We use a suite of high resolution state-of-the-art N-body Dark Matter simulations of chameleon $f(R)$ gravity to study the higher order volume averaged correlation functions $\bar{\xi_n}$ together with the hierarchical $n$-th order correlation amplitudes $S_n=\bar{\xi}_n/\bar{\xi}_2^{n-1}$ and density distribution functions (PDF). We show that under the non-linear modifications of gravity the hierarchical scaling of the reduced cumulants is preserved. This is however characterised by significant changes of both the $\bar{\xi_n}$ and $S_n$'s values and their scale dependence with respect to General Relativity gravity (GR). In addition, we measure a significant increase of the non linear $\sigma_8$ parameter reaching $14, 5$ and 0.5% in excess of the GR value for the three flavours of our $f(R)$ models. We further note that the values of the reduced cumulants up to order $n=9$ are significantly increased in $f(R)$ gravity for all our models at small scales $R\simlt 30\hmpc$. In contrast the values of the hierarchical amplitudes, $S_n$'s, are smaller in $f(R)$ indicating that the modified gravity density distribution functions are deviating from the GR case. Furthermore we find that the redshift evolution of relative deviations of the $f(R)$ hierarchical correlation amplitudes is fastest at high and moderate redshifts $1\leq z \leq4$. The growth of these deviations significantly slows down in the low redshift universe. We also compute the PDFs and show that for scales below $\sim 20\hmpc$ they are significantly shifted in $f(R)$ gravity towards the low densities. Finally we discuss the implications of our theoretical predictions for measurements of the hierarchical clustering in galaxy redshifts surveys, including the important problems of the galaxy biasing and redshifts space distortions.

We present a clustering analysis of X-ray selected AGN by compiling X-ray samples from the literature and re-estimating the dark matter (DM) halo masses of AGN in a uniform manner. We find that moderate luminosity AGN (Lx(2-10 keV)=10^42-10^44 erg/sec) in the z=0-1.3 Universe are typically found in DM haloes with masses of ~10^13 Msun. We then compare our findings to the theoretical predictions of the coupled galaxy and black hole formation model GALFORM. We find good agreement when our calculation includes the hot-halo mode of accretion onto the central black hole. This type of accretion, which is additional to the common cold accretion during disk instabilities and galaxy mergers, is tightly coupled to the AGN feedback in the model. The hot-halo mode becomes prominent in DM haloes with masses greater than ~10^12.5 Msun, where AGN feedback typically operates, giving rise to a distinct class of moderate luminosity AGN that inhabit rich clusters and superclusters. Cold gas fuelling of the black hole cannot produce the observationally inferred DM halo masses of X-ray AGN. Switching off AGN feedback in the model results in a large population of luminous quasars (Lx(2-10 keV) > 10^44 erg/sec) in DM haloes with masses up to ~10^14 Msun, which is inconsistent with the observed clustering of quasars. The abundance of hot-halo AGN decreases significantly in the z~3-4 universe. At such high redshifts, the cold accretion mode is solely responsible for shaping the environment of moderate luminosity AGN. Our analysis supports two accretion modes (cold and hot) for the fuelling of supermassive black holes and strongly underlines the importance of AGN feedback in cosmological models both of galaxy formation and black hole growth.

Quasars represent the brightest Active Galactic Nuclei (AGN) in the Universe and are thought to indicate the location of prodigiously growing Black Holes (BHs), with luminosities as high as 10^48 erg/sec. It is often expected though that such an extremely energetic process will take place in the most massive bound structures in the dark matter (DM) distribution. We show that in contrast to this expectation, in a galaxy formation model which includes AGN feedback, quasars are predicted to live in average DM halo environments with typical masses of a few times 10^12 Msun. This fundamental prediction arises from the fact that quasar activity (i.e., BH accretion with luminosity greater than 10^46 erg/sec) is inhibited in DM haloes where AGN feedback operates. The galaxy hosts of quasars in our simulations are identified with over massive (in gas and stars) spheroidal galaxies, in which BH accretion is triggered via a galaxy merger or secular processes. We further show that the z=0 descendants of high redshift (z~6) QSOs span a wide range of morphologies, galaxy and halo masses. The z~6 BHs typically grow only by a modest factor by the present day. Remarkably, high redshift QSOs never inhabit the largest DM haloes at that time and their descendants are very seldom found in the most massive haloes at z=0. We also show that observationally it is very likely to find an enhancement in the abundance of galaxies around quasars at z~5. However, these enhancements are considerably weaker compared to the overdensities expected at the extreme peaks of the DM distribution. Thus, it is very unlikely that a quasar detected in the $z\gtrsim5$ Universe pinpoints the location of the progenitors of superclusters in the local Universe.

We study the effect of baryons on the abundance of structures and substructures in a Lambda-CDM cosmology, using a pair of high resolution cosmological simulations from the GIMIC project. Both simulations use identical initial conditions, but while one contains only dark matter, the other also includes baryons. We find that gas pressure, reionisation, supernova feedback, stripping, and truncated accretion systematically reduce the total mass and the abundance of structures below ~10^12 solar masses compared to the pure dark matter simulation. Taking this into account and adopting an appropriate detection threshold lowers the abundance of observed galaxies with maximum circular velocities below 100 km/s, significantly reducing the reported discrepancy between Lambda-CDM and the measured HI velocity function of the ALFALFA survey. We also show that the stellar-to-total mass ratios of galaxies with stellar masses of ~10^5 - 10^7 solar masses inferred from abundance matching of the (sub)halo mass function to the observed galaxy mass function increase by a factor of ~2. In addition, we find that an important fraction of low-mass subhaloes are completely devoid of stars. Accounting for the presence of dark subhaloes below 10^10 solar masses further reduces the abundance of observable objects, and leads to an additional increase in the inferred stellar-to-total mass ratio by factors of 2 - 10 for galaxies in haloes of 10^9 - 10^10 solar masses. This largely reconciles the abundance matching results with the kinematics of individual dwarf galaxies in Lambda-CDM. We propose approximate corrections to the masses of objects derived from pure dark matter calculations to account for baryonic effects.

Models of galaxy formation in a CDM universe predict that massive galaxies are surrounded by a hot, quasi-hydrostatic circumgalactic corona of slowly cooling gas, predominantly accreted from the IGM. This prediction is borne out by the cosmological hydrodynamical simulations of Crain et al., which reproduce scaling relations between the X-ray and optical properties of nearby disc galaxies. Such coronae are metal poor, but observations of the X-ray emitting circumgalactic medium (CGM) of local galaxies typically indicate enrichment to near-solar iron abundance, potentially signalling a shortcoming in galaxy formation models. We show here that, while the hot CGM of galaxies formed in the simulations is metal poor in a mass-weighted sense, its X-ray luminosity-weighted metallicity is often close to solar. This bias arises because the soft X-ray emissivity of a typical 0.1 keV corona is dominated by collisionally-excited metal ions that are synthesised in stars and recycled into the hot CGM. We find that these metals are ejected primarily by stars that form in-situ to the main progenitor of the galaxy, rather than in satellites or external galaxies. The enrichment of the hot CGM therefore proceeds in an inside-out fashion throughout the assembly of the galaxy: metals are transported from the central galaxy by SNe-driven winds and convection over several gigayears, establishing a strong negative radial metallicity gradient. Whilst metal ions synthesised by stars are necessary to produce the X-ray emissivity that enables the hot CGM of isolated galaxies to be detected, the electrons that collisionally excite them are equally important. Since our simulations indicate that the electron density of hot coronae is dominated by the metal-poor gas accreted from the IGM, we infer that the hot CGM observed via X-ray emission is the outcome of both hierarchical accretion and stellar recycling.

The Galaxy And Mass Assembly (GAMA) survey is a multiwavelength photometric and spectroscopic survey, using the AAOmega spectrograph on the Anglo-Australian Telescope to obtain spectra for up to ~300000 galaxies over 280 square degrees, to a limiting magnitude of r_pet < 19.8 mag. The target galaxies are distributed over 0<z<0.5 with a median redshift of z~0.2, although the redshift distribution includes a small number of systems, primarily quasars, at higher redshifts, up to and beyond z=1. The redshift accuracy ranges from sigma_v~50km/s to sigma_v~100km/s depending on the signal-to-noise of the spectrum. Here we describe the GAMA spectroscopic reduction and analysis pipeline. We present the steps involved in taking the raw two-dimensional spectroscopic images through to flux-calibrated one-dimensional spectra. The resulting GAMA spectra cover an observed wavelength range of 3750<lambda<8850 A at a resolution of R~1300. The final flux calibration is typically accurate to 10-20%, although the reliability is worse at the extreme wavelength ends, and poorer in the blue than the red. We present details of the measurement of emission and absorption features in the GAMA spectra. These measurements are characterised through a variety of quality control analyses detailing the robustness and reliability of the measurements. We illustrate the quality of the measurements with a brief exploration of elementary emission line properties of the galaxies in the GAMA sample. We demonstrate the luminosity dependence of the Balmer decrement, consistent with previously published results, and explore further how Balmer decrement varies with galaxy mass and redshift. We also investigate the mass and redshift dependencies of the [NII]/Halpha vs [OIII]/Hbeta spectral diagnostic diagram, commonly used to discriminate between star forming and nuclear activity in galaxies.

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

We have performed a series of numerical experiments to investigate how the primordial thermal velocities of fermionic dark matter particles affect the physical and phase space density profiles of the dark matter haloes into which they collect. The initial particle velocities induce central cores in both profiles, which can be understood in the framework of phase space density theory. We find that the maximum coarse-grained phase space density of the simulated haloes (computed in 6 dimensional phase space using the EnBid code) is very close to the theoretical fine-grained upper bound, while the pseudo phase space density, Q ~ {\rho}/{\sigma}^3, overestimates the maximum phase space density by up to an order of magnitude. The density in the inner regions of the simulated haloes is well described by a 'pseudo-isothermal' profile with a core. We have developed a simple model based on this profile which, given the observed surface brightness profile of a galaxy and its central velocity dispersion, accurately predicts its central phase space density. Applying this model to the dwarf spheroidal satellites of the Milky Way yields values close to 0.5 keV for the mass of a hypothetical thermal warm dark matter particle, assuming the satellite haloes have cores produced by warm dark matter free streaming. Such a small value is in conflict with the lower limit of 1.2 keV set by observations of the Lyman-{\alpha} forest. Thus, if the Milky Way dwarf spheroidal satellites have cores, these are likely due to baryonic processes associated with the forming galaxy, perhaps of the kind proposed by Navarro, Eke and Frenk and seen in recent simulations of galaxy formation in the cold dark matter model.

We combine the six high-resolution Aquarius dark matter simulations with a semi-analytic galaxy formation model to investigate the properties of the satellites of Milky Way-like galaxies. We find good correspondence with the observed luminosity function, luminosity-metallicity relation and radial distribution of the Milky Way satellites. The star formation histories of the dwarf galaxies in our model vary widely, in accordance with what is seen observationally. Ram-pressure stripping of hot gas from the satellites leaves a clear imprint of the environment on the characteristics of a dwarf galaxy. We find that the fraction of satellites dominated by old populations of stars matches observations well. However, the internal metallicity distributions of the model satellites appear to be narrower than observed. This may indicate limitations in our treatment of chemical enrichment, which is based on the instantaneous recycling approximation. Our model works best if the dark matter halo of the Milky Way has a mass of ~8 x 10^11 Msun, in agreement with the lower estimates from observations. The galaxy that resembles the Milky Way the most also has the best matching satellite luminosity function, although it does not contain an object as bright as the SMC or LMC. Compared to other semi-analytic models and abundance matching relations we find that central galaxies reside in less massive haloes, but the halo mass-stellar mass relation in our model is consistent both with hydrodynamical simulations and with recent observations.

We investigate the luminosity functions (LFs) and projected number density profiles of galactic satellites around isolated primaries of different luminosities. We measure these quantities for model satellites placed into the Millennium and Millennium II dark matter simulations by the GALFORM semi-analytic galaxy formation model for different bins of primary galaxy magnitude and we investigate their dependence on satellite luminosity. We compare our model predictions to the data of Guo et al. from the Sloan Digital Sky Survey Data Release 8 (SDSS DR8). First, we use a mock light-cone catalogue to verify that the method we used to count satellites in the SDSS DR8 is unbiased. We find that the radial distributions of model satellites are similar to those around comparable primary galaxies in the SDSS DR8, with only slight differences at low luminosities and small projected radii. However, when splitting the satellites by colour, the model and SDSS satellite systems no longer resemble one another, with many red model satellites, in contrast to the dominant blue fraction at similar luminosity in SDSS. The few model blue satellites are also significantly less centrally concentrated in the halo of their stacked primary than their SDSS counterparts. The implications of this result for the GALFORM model are discussed.

Cold dark matter models predict the existence of a large number of substructures within dark matter halos. If the cold dark matter consists of weakly interacting massive particles, their annihilation within these substructures could lead to diffuse GeV emission that would dominate over the annihilation signal of the host halo. In this work we search for GeV emission from three nearby galaxy clusters: Coma, Virgo and Fornax. We first remove known extragalactic and galactic diffuse gamma-ray backgrounds and point sources from the Fermi 2-year catalog and find a significant residual diffuse emission in all three clusters. We then investigate whether this emission is due to (i) unresolved point sources; (ii) dark matter annihilation; or (iii) cosmic rays (CR). Using 45 months of Fermi-LAT data we detect several new point sources (not present in the Fermi 2-year point source catalogue) which contaminate the signal previously analyzed by Han et al.(arxiv:1201.1003). Including these and accounting for the effects of undetected point sources, we find no significant detection of extended emission from the three clusters studied. Instead, we determine upper limits on emission due to dark matter annihilation and cosmic rays. For Fornax and Virgo the limits on CR emission are consistent with theoretical models, but for Coma the upper limit is a factor of 2 below the theoretical expectation. Allowing for systematic uncertainties associated with the treatment of CR, the upper limits on the cross section for dark matter annihilation from our clusters are more stringent than those from analyses of dwarf galaxies in the Milky Way. We rule out the thermal cross section for supersymmetric dark matter particles for masses as large as 100 GeV (depending on the annihilation channel).

Using GALFORM, a semi-analytical model of galaxy formation in the Lambda cold dark matter cosmology, we study the rest-frame ultraviolet (UV) colours of Lyman-break galaxies (LBGs) in the redshift range 2.5 < z < 10. As the impact of dust on UV luminosity can be dramatic, our model includes a self-consistent computation of dust attenuation based on a radiative transfer model. We find that intrinsically brighter galaxies suffer stronger dust attenuation than fainter ones, though the relation has a large scatter. The model predicts galaxies with UV colours consistent with the colour selection regions designed to select LBGs in observational surveys. We find that the drop-out technique that selects LBGs based on two rest-frame UV colours is robust and effective, selecting more than 70 per cent of UV bright galaxies at a given redshift. We investigate the impact on the predicted UV colours of varying selected model parameters. We find that the UV colours are most sensitive to the modelling of dust attenuation and in particular, to the extinction curve used in the radiative transfer calculation. If we assume a Milky Way dust extinction curve, the predicted UV continuum slopes are, in general, bluer than observed. However, we find that the opposite is true when using the Small Magellanic Cloud dust extinction curve. This demonstrates the strong dependence of UV colours on dust properties and highlights the inadequacy of using the UV continuum slope as a tracer of dust attenuation without any further knowledge of the galaxy inclination or dust characteristics in high redshift galaxies.

One of the important unknowns of current cosmology concerns the effects of the large scale distribution of matter on the formation and evolution of dark matter haloes and galaxies. One main difficulty in answering this question lies in the absence of a robust and natural way of identifying the large scale environments and their characteristics. This work summarizes the NEXUS+ formalism which extends and improves our multiscale scale-space MMF method. The new algorithm is very successful in tracing the Cosmic Web components, mainly due to its novel filtering of the density in logarithmic space. The method, due to its multiscale and hierarchical character, has the advantage of detecting all the cosmic structures, either prominent or tenuous, without preference for a certain size or shape. The resulting filamentary and wall networks can easily be characterized by their direction, thickness, mass density and density profile. These additional environmental properties allows to us to investigate not only the effect of environment on haloes, but also how it correlates with the environment characteristics.

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 describe a methodology to accurately compute halo mass functions, progenitor mass functions, merger rates and merger trees in non-cold dark matter universes using a self-consistent treatment of the generalized extended Press-Schechter formalism. Our approach permits rapid exploration of the subhalo population of galactic halos in dark matter models with a variety of different particle properties or universes with rolling, truncated, or more complicated power spectra. We make detailed comparisons of analytically derived mass functions and merger histories with recent warm dark matter cosmological N-body simulations, and find excellent agreement. We show that, once the accretion of smoothly distributed matter is accounted for, coarse-grained statistics such as the mass accretion history of halos can be almost indistinguishable between cold and warm dark matter cases. However, the halo mass function and progenitor mass functions differ significantly, with the warm dark matter cases being strongly suppressed below the free-streaming scale of the dark matter. We demonstrate the importance of using the correct solution for the excursion set barrier first-crossing distribution in warm dark matter - if the solution for a flat barrier is used instead the truncation of the halo mass function is much slower, leading to an overestimate of the number of low mass halos.

We review the current standard model for the evolution of cosmic structure, tracing its development over the last forty years and focusing specifically on the role played by numerical simulations and on aspects related to the nature of dark matter.

Star-forming galaxies which are too faint to be detected individually produce intensity fluctuations in the cosmic background light. This contribution needs to be taken into account as a foreground when using the primordial signal to constrain cosmological parameters. The extragalactic fluctuations are also interesting in their own right as they depend on the star formation history of the Universe and the way in which this connects with the formation of cosmic structure. We present a new framework which allows us to predict the occupation of dark matter haloes by star-forming galaxies and uses this information, in conjunction with an N-body simulation of structure formation, to predict the power spectrum of intensity fluctuations in the infrared background. We compute the emission from galaxies at far-infrared, millimetre and radio wavelengths. Our method gives accurate predictions for the clustering of galaxies both for the one halo and two halo terms. We illustrate our new framework using a previously published model which reproduces the number counts and redshift distribution of galaxies selected by their emission at $850\,\mu$m. Without adjusting any of the model parameters, the predictions show encouraging agreement at high frequencies and on small angular scales with recent estimates of the extragalactic fluctuations in the background made from early data analysed by the Planck Collaboration. There are, however, substantial discrepancies between the model predictions and observations on large angular scales and at low frequencies, which illustrates the usefulness of the intensity fluctuations as a constraint on galaxy formation models.

Context: The mass-metallicity relationship (MMR) of star-forming galaxies is well-established, however there is still some disagreement with respect to its exact shape and its possible dependence on other observables. Aims: We measure the MMR in the Galaxy And Mass Assembly (GAMA) survey. We compare our measured MMR to that measured in the Sloan Digital Sky Survey (SDSS) and study the dependence of the MMR on various selection criteria to identify potential causes for disparities seen in the literature. Methods: We use strong emission line ratio diagnostics to derive oxygen abundances. We then apply a range of selection criteria for the minimum signal-to-noise in various emission lines, as well as the apparent and absolute magnitude to study variations in the inferred MMR. Results: The shape and position of the MMR can differ significantly depending on the metallicity calibration and selection used. After selecting a robust metallicity calibration amongst those tested, we find that the mass-metallicity relation for redshifts 0.061< z<0.35 in GAMA is in reasonable agreement with that found in the SDSS despite the difference in the luminosity range probed. Conclusions: In view of the significant variations of the MMR brought about by reasonable changes in the sample selection criteria and method, we recommend that care be taken when comparing the MMR from different surveys and studies directly. We also conclude that there could be a modest level of evolution over 0.06<z<0.35 within the GAMA sample.

We study the spatial distribution of satellite galaxies around isolated primaries using the Sloan Digital Sky Survey (SDSS) spectroscopic and photometric galaxy catalogues. We select isolated primaries from the spectroscopic sample and search for potential satellites in the much deeper photometric sample. For specific luminosity primaries we obtain robust statistical results by stacking as many as ~50, 000 galaxy systems. We find no evidence for any anistropy in the satellite galaxy distribution relative to the major axes of the primaries. We derive accurate projected number density profiles of satellites down to 4 magnitudes fainter than their primaries. We find the normalized satellite profiles generally have a universal form and can be well fitted by projected NFW profiles. The NFW concentration parameter increases with decreasing satellite luminosity while being independent of the luminosity of the primary except for very bright primaries. The profiles of the faintest satellites show deviations from the NFW form with an excess at small galactocentric projected distances. In addition, we quantify how the radial distribution of satellites depends on the colour of the satellites and on the colour and concentration of their primaries.

We use the GAMA I dataset combined with GALEX, SDSS and UKIDSS imaging to construct the low-redshift (z<0.1) galaxy luminosity functions in FUV, NUV, ugriz, and YJHK bands from within a single well constrained volume of 3.4 x 10^5 (Mpc/h)^{3}. The derived luminosity distributions are normalised to the SDSS DR7 main survey to reduce the estimated cosmic variance to the 5 per cent level. The data are used to construct the cosmic spectral energy distribution (CSED) from 0.1 to 2.1 \mum free from any wavelength dependent cosmic variance for both the elliptical and non-elliptical populations. The two populations exhibit dramatically different CSEDs as expected for a predominantly old and young population respectively. Using the Driver et al. (2008) prescription for the azimuthally averaged photon escape fraction, the non-ellipticals are corrected for the impact of dust attenuation and the combined CSED constructed. The final results show that the Universe is currently generating (1.8 +/- 0.3) x 10^{35} h W Mpc^{-3} of which (1.2 +/- 0.1) x 10^{35} h W Mpc^{-3} is directly released into the inter-galactic medium and (0.6 +/- 0.1) x 10^{35} h W Mpc^{-3} is reprocessed and reradiated by dust in the far-IR. Using the GAMA data and our dust model we predict the mid and far-IR emission which agrees remarkably well with available data. We therefore provide a robust description of the pre- and post dust attenuated energy output of the nearby Universe from 0.1micron to 0.6mm. The largest uncertainty in this measurement lies in the mid and far-IR bands stemming from the dust attenuation correction and its currently poorly constrained dependence on environment, stellar mass, and morphology.

We examine the rotation rates, sizes, and star formation (SF) efficiencies of a representative population of simulated disc galaxies extracted from the Galaxies-Intergalactic Medium Interaction Calculation (GIMIC) suite of cosmological hydrodynamic simulations. These simulations include efficient, but energetically feasible supernova feedback, but have not been tuned in any way to produce 'realistic' disc galaxies. Yet, they generate a large number of discs, without requiring extremely high resolution. Over the wide galaxy stellar mass range, 9.0 < log10[Mstar (Msun)] < 10.5, the simulations reproduce the observed Tully-Fisher relation, the rotation curves of disc galaxies in bins of stellar mass, the mass-size relation of disc galaxies, the optical rotation to virial circular velocity ratio (Vopt/Vvir), and the SF efficiencies of disc galaxies as inferred from stacked weak lensing and stacked satellite kinematics observations. They also reproduce the specific star formation rates of ~L* galaxies but predict too low levels of star formation for low-mass galaxies, which is plausibly due to the finite resolution of the simulations. At higher stellar masses, log10[Mstar (\Msun)] > 10.6, the simulated galaxies are too concentrated and have too high SF efficiencies. We conjecture that this shortcoming reflects the neglect of feedback from accreting supermassive black holes in these simulations. We conclude that it is possible to generate a representative population of disc galaxies that reproduces many of the observed trends of local disc galaxies using standard numerical hydrodynamic techniques and a plausible implementation of the "subgrid" astrophysical processes thought to be relevant to galaxy formation.

In the simplest scenario, disk galaxies form predominantly in halos with high angular momentum and quiet recent assembly history, whereas spheroids are the slowly-rotating remnants of repeated merging events. We explore these assumptions using one hundred systems with halo masses similar to that of the Milky Way, identified in a series of cosmological gasdynamical simulations GIMIC. At z=0, the simulated galaxies exhibit a wide variety of morphologies, from dispersion-dominated spheroids to pure disk galaxies. Surprisingly, these morphological features are very poorly correlated with their halo properties: disks form in halos with high and low net spin, and mergers play a negligible role in the formation of spheroid stars, most of which form in-situ. More important to morphology is the coherent alignment of the angular momentum of baryons that accrete over time to form a galaxy. Spheroids tend to form when the spin of newly-accreted gas is misaligned with that of the extant galaxy, leading to the episodic formation of stars with different kinematics that cancel out the net rotation of the system. Disks, on the other hand, form out of gas that flows in with similar angular momentum to that of earlier-accreted material. Gas accretion from a hot corona thus favours disk formation, whereas gas that flows "cold", often along separate, misaligned filaments, favours the formation of spheroids. In this scenario, most spheroids consist of superpositions of stellar components with distinct kinematics, age, and metallicity, an arrangement that might survive to the present day given the paucity of major mergers. Since angular momentum is acquired largely at turnaround, morphology is imprinted early by the interplay of the tidal field and the shape of the material destined to form the galaxy.

We report evidence for extended gamma-ray emission from the Virgo, Fornax and Coma clusters based on a maximum-likelihood analysis of the 3-year Fermi-LAT data. For all three clusters, excess emission is observed within three degrees of the center, peaking at the GeV scale. This emission cannot be accounted for by known Fermi sources or by the galactic and extragalactic backgrounds. If interpreted as annihilation emission from supersymmetric dark matter (DM) particles, the data prefer models with a particle mass in the range 20-60 GeV annihilating into the b-bbar channel, or 2-10 GeV and >1 TeV annihilating into mu-mu final states. Our results are consistent with those obtained by Hooper and Linden from a recent analysis of Fermi-LAT data in the region of the Galactic Centre. An extended DM annihilation profile dominated by emission from substructures is preferred over a simple point source model. The significance of DM detection is 4.4 sigma in Virgo and lower in the other two clusters. We also consider the possibility that the excess emission arises from cosmic ray (CR) induced gamma-rays, and infer a CR level within a factor of three of that expected from analytical models. However, the significance of a CR component is lower than the significance of a DM component, and there is no need for such a CR component in the presence of a DM component in the preferred DM mass range. We also set flux and cross-section upper limits for DM annihilation into the b-bbar and mu-mu channels in all three clusters.

The Fermi-LAT collaboration has recently reported the detection of angular power above the photon noise level in the diffuse gamma-ray background between 1 and 50 GeV. Such signal can be used to constrain a possible contribution from Dark-Matter-induced photons. We estimate the intensity and features of the angular power spectrum (APS) of this potential Dark Matter (DM) signal, for both decaying and annihilating DM candidates, by constructing template all-sky gamma-ray maps for the emission produced in the galactic halo and its substructures, as well as in extragalactic (sub)halos. The DM distribution is given by state-of-the-art N-body simulations of cosmic structure formation, namely Millennium-II for extragalactic (sub)halos, and Aquarius for the galactic halo and its subhalos. We use a hybrid method of extrapolation to account for (sub)structures that are below the resolution limit of the simulations, allowing us to estimate the total emission all the way down to the minimal self-bound halo mass. We describe in detail the features appearing in the APS of our template maps and we estimate the effect of various uncertainties such as the value of the minimal halo mass, the fraction of substructures hosted in a halo and the shape of the DM density profile. Our results indicate that the fluctuation APS of the DM-induced emission is of the same order as the Fermi-LAT APS, suggesting that one can constrain this hypothetical emission from the comparison with the measured anisotropy. We also quantify the uncertainties affecting our results, finding "theoretical error bands" spanning more than two orders of magnitude and dominated (for a given particle physics model) by our lack of knowledge of the abundance of low-mass (sub)halos.

[abridged] We introduce the Phoenix Project, a set of $\Lambda$CDM simulations of the dark matter component of nine rich galaxy clusters. Each cluster is simulated at least at two different numerical resolutions. For eight of them, the highest resolution corresponds to $\sim 130$ million particles within the virial radius, while for one this number is over one billion. We study the structure and substructure of these systems and contrast them with six galaxy-sized dark matter haloes from the Aquarius Project, simulated at comparable resolution. This comparison highlights the approximate mass invariance of CDM halo structure and substructure. We find little difference in the spherically-averaged mass, pseudo-phase-space density, and velocity anisotropy profiles of Aquarius and Phoenix haloes. When scaled to the virial properties of the host halo, the abundance and radial distribution of subhaloes are also very similar, despite the fact that Aquarius and Phoenix haloes differ by roughly three decades in virial mass. The most notable difference is that cluster haloes have been assembled more recently and are thus significantly less relaxed than galaxy haloes, which leads to decreased regularity, increased halo-to-halo scatter and sizable deviations from the mean trends. This accentuates the effects of the strong asphericity of individual clusters on surface density profiles, which may vary by up to a factor of three at a given radius, depending on projection. The high apparent concentration reported for some strong-lensing clusters might very well reflect these effects. A more recent assembly also explains why substructure in some Phoenix haloes is slightly more abundant than in Aquarius, especially in the inner regions. Resolved subhaloes nevertheless contribute only $11 \pm 3%$ of the virial mass in Phoenix clusters. .

We investigate the spatial distribution of galactic satellites in high resolution simulations of structure formation in the LCDM model: the Aquarius dark matter simulations of individual halos and the Millennium II simulation of a large cosmological volume. To relate the simulations to observations of the Milky Way we use two alternative models to populate dark halos with "visible" galaxies: a semi-analytic model of galaxy formation and an abundance matching technique. We find that the radial density profile of massive satellites roughly follows that of the dark matter halo (unlike the distribution of dark matter subhalos). Furthermore, our two galaxy formation models give results consistent with the observed profile of the 11 classical satellites of the Milky Way. Our simulations predict that larger, fainter samples of satellites should still retain this profile at least up to samples of 100 satellites. The angular distribution of the classical satellites of the Milky Way is known to be highly anisotropic. Depending on the exact measure of flattening, 5--10 per cent of satellite systems in our simulations are as flat as the Milky Way's and this fraction does not change when we correct for possible obscuration of satellites by the Galactic disk. A moderate flattening of satellite systems is a general property of LCDM, best understood as the consequence of preferential accretion along filaments of the cosmic web. Accretion of a single rich group of satellites can enhance the flattening due to such anisotropic accretion. We verify that a typical Milky Way-mass CDM halo does not acquire its 11 most massive satellites from a small number of rich groups. Single--group accretion becomes more likely for less massive satellites. Our model predictions should be testable with forthcoming studies of satellite systems in other galaxies and surveys of fainter satellites in the Milky Way.

We measure the two-point angular correlation function of a sample of 4,289,223 galaxies with r < 19.4 mag from the Sloan Digital Sky Survey as a function of photometric redshift, absolute magnitude and colour down to M_r - 5log h = -14 mag. Photometric redshifts are estimated from ugriz model magnitudes and two Petrosian radii using the artificial neural network package ANNz, taking advantage of the Galaxy and Mass Assembly (GAMA) spectroscopic sample as our training set. The photometric redshifts are then used to determine absolute magnitudes and colours. For all our samples, we estimate the underlying redshift and absolute magnitude distributions using Monte-Carlo resampling. These redshift distributions are used in Limber's equation to obtain spatial correlation function parameters from power law fits to the angular correlation function. We confirm an increase in clustering strength for sub-L* red galaxies compared with ~L* red galaxies at small scales in all redshift bins, whereas for the blue population the correlation length is almost independent of luminosity for ~L* galaxies and fainter. A linear relation between relative bias and log luminosity is found to hold down to luminosities L~0.03L*. We find that the redshift dependence of the bias of the L* population can be described by the passive evolution model of Tegmark & Peebles (1998). A visual inspection of a random sample of our r < 19.4 sample of SDSS galaxies reveals that about 10 per cent are spurious, with a higher contamination rate towards very faint absolute magnitudes due to over-deblended nearby galaxies. We correct for this contamination in our clustering analysis.

Recent studies suggest that only three of the twelve brightest satellites of the Milky Way (MW) inhabit dark matter halos with maximum circular velocity, V_max, exceeding 30km/s. This is in apparent contradiction with the LCDM simulations of the Aquarius Project, which suggest that MW-sized halos should have at least 8 subhalos with V_max>30km/s. The absence of luminous satellites in such massive subhalos is thus puzzling and may present a challenge to the LCDM paradigm. We note, however, that the number of massive subhalos depends sensitively on the (poorly-known) virial mass of the Milky Way, and that their scarcity makes estimates of their abundance from a small simulation set like Aquarius uncertain. We use the Millennium Simulation series and the invariance of the scaled subhalo velocity function (i.e., the number of subhalos as a function of \nu, the ratio of subhalo V_max to host halo virial velocity, V_200) to secure improved estimates of the abundance of rare massive subsystems. In the range 0.1<\nu<0.5, N_sub(>\nu) is approximately Poisson-distributed about an average given by <N_sub>=10.2x(\nu/0.15)^(-3.11). This is slightly lower than in Aquarius halos, but consistent with recent results from the Phoenix Project. The probability that a LCDM halo has 3 or fewer subhalos with V_max above some threshold value, V_th, is then straightforward to compute. It decreases steeply both with decreasing V_th and with increasing halo mass. For V_th=30km/s, ~40% of M_halo=10^12 M_sun halos pass the test; fewer than 5% do so for M_halo>= 2x10^12 M_sun; and the probability effectively vanishes for M_halo>= 3x 10^12 M_sun. Rather than a failure of LCDM, the absence of massive subhalos might simply indicate that the Milky Way is less massive than is commonly thought.

The formation of galaxies is regulated by a balance between the supply of gas and the rate at which it is ejected. Traditional explanations of gas ejection equate the energy required to escape the galaxy or host halo to an estimate for the energy yield from supernovae. This yield is usually assumed to be a constant fraction of the total available from the supernova, or is derived from the assumption of a consistent momentum yield. By applying these ideas in the context of a cold dark matter cosmogony, we derive a 1st-order analytic connection between these working assumptions and the expected relationship between baryon content and galaxy circular velocity, and find that these quick predictions straddle recent observational estimates. To examine the premises behind these theories in more detail, we then explore their applicability to a set of gasdynamical simulations of idealised galaxies. We show that different premises dominate to differing degrees in the simulated outflow, depending on the mass of the system and the resolution with which it is simulated. Using this study to anticipate the emergent behaviour at arbitrarily high resolution, we motivate more comprehensive analytic model which allows for the range of velocities with which the gas may exit the system, and incorporates both momentum and energy-based constraints on the outflow. Using a trial exit velocity distribution, this is shown to be compatible with the observed baryon fractions in intermediate-mass systems, but implies that current estimates for low-mass systems can not be solely accounted for by supernova winds under commonly-held assumptions.

We present a method for attaining sub-arcsecond pointing stability during sub- orbital balloon flights, as designed for in the High Altitude Lensing Observatory (HALO) concept. The pointing method presented here has the potential to perform near-space quality optical astronomical imaging at 1-2% of the cost of space-based missions. We also discuss an architecture that can achieve sufficient thermomechanical stability to match the pointing stability. This concept is motivated by advances in the development and testing of Ultra Long Duration Balloon (ULDB) flights which promise to allow observation campaigns lasting more than three months. The design incorporates a multi-stage pointing architecture comprising: a gondola coarse azimuth control system, a multi-axis nested gimbal frame structure with arcsecond stability, a telescope de-rotator to eliminate field rotation, and a fine guidance stage consisting of both a telescope mounted angular rate sensor and guide CCDs in the focal plane to drive a fast-steering mirror. We discuss the results of pointing tests together with a preliminary thermo-mechanical analysis required for sub-arcsecond pointing at high altitude. Possible future applications in the areas of wide-field surveys and exoplanet searches are also discussed.

We present measurements of the triaxial dark matter halo shapes and alignment correlation functions in the Millennium and Millennium-2 dark matter N-body simulations. These two simulations allow us to measure the distributions of halo shapes down to 10% of the virial radius over a halo mass range of 6E9 - 2E14 M_sun/h. We largely confirm previous results on the distributions of halo axis ratios as a function of halo mass, but we find that the median angle between halo major axes at different halo radii can vary by a factor of 2 between the Millennium-1 and 2 simulations because of the different mass resolution. Thus, error in the shape determinations from limited resolution is potentially degenerate with the misalignment of halo inner and outer shapes used to constrain Brightest Cluster Galaxy alignments in previous works. We also present simplifying parameterizations for the 3-D halo-mass alignment correlation functions that are necessary ingredients for triaxial halo models of large-scale structure and models of galaxy intrinsic alignments as contaminants for cosmic shear surveys. We measure strong alignments between halos of all masses and the surrounding dark matter overdensities out to several tens of Mpc/h, in agreement with observed shear-galaxy and cluster shape correlations. We use these measurements to forecast the contribution to the weak lensing signal around galaxy clusters from correlated mass along the line-of-sight. For prolate clusters with major axes aligned with the line-of-sight the fraction of the weak lensing signal from mass external to the cluster can be twice that predicted if the excess halo alignment correlation is assumed to be zero.

We compare the results of thirteen cosmological gasdynamical codes used to simulate the formation of a galaxy in the LCDM structure formation paradigm. The various runs differ in their hydrodynamical treatment (SPH, moving-mesh and AMR) but share the same initial conditions and adopt their latest published model of cooling, star formation and feedback. Despite the common halo assembly history, we find large code-to-code variations in the stellar mass, size, morphology and gas content of the galaxy at z=0, due mainly to the different implementations of feedback. Compared with observation, most codes tend to produce an overly massive galaxy, smaller and less gas-rich than typical spirals, with a massive bulge and a declining rotation curve. A stellar disk is discernible in most simulations, though its prominence varies widely from code to code. There is a well-defined trend between the effects of feedback and the severity of the disagreement with observation. Models that are more effective at limiting the baryonic mass of the galaxy come closer to matching observed galaxy scaling laws, but often to the detriment of the disk component. Our conclusions hold at two different numerical resolutions. Some differences can also be traced to the numerical techniques: more gas seems able to cool and become available for star formation in grid-based codes than in SPH. However, this effect is small compared to the variations induced by different feedback prescriptions. We conclude that state-of-the-art simulations cannot yet uniquely predict the properties of the baryonic component of a galaxy, even when the assembly history of its host halo is fully specified. Developing feedback algorithms that can effectively regulate the mass of a galaxy without hindering the formation of high-angular momentum stellar disks remains a challenge.

The flux-ratio anomalies observed in multiply-lensed quasar images are most plausibly explained as the result of perturbing structures superposed on the underlying smooth matter distribution of the primary lens. The cold dark matter cosmological model predicts that a large number of substructures should survive inside larger halos but, surprisingly, this population alone has been shown to be insufficient to explain the observed distribution of the flux ratios of quasar's multiple images. Other halos (and their own subhalos) projected along the line of sight to the primary lens have been considered as additional source of perturbation. In this work, we use ray tracing through the Millennium II simulation to investigate the importance of projection effects due to halos and subhalos of mass m>1E8 Msun/h and extend our analysis to lower masses, m>1E6 Msun/h, using Monte-Carlo halo distributions. We find that the magnitude of the violation depends strongly on the density profile and concentration of the intervening halos, but clustering plays only a minor role. For a typical lensing geometry (lens at redshift 0.6 and source at redshift 2), background haloes (behind the main lens) are more likely to cause a violation than foreground halos. We conclude that line-of-sight structures can be as important as intrinsic substructures in causing flux-ratio anomalies. The combined effect of perturbing structures within the lens and along the line of sight in the LCDM universe results in a cusp-violation probability of 20-30%. This alleviates the discrepancy between models and current data, but a larger observational sample is required for a stronger test of the theory.

High resolution N-body simulations of galactic cold dark matter haloes indicate that we should expect to find a few satellite galaxies around the Milky Way whose haloes have a maximum circular velocity in excess of 40 kms. Yet, with the exception of the Magellanic Clouds and the Sagittarius dwarf, which likely reside in subhaloes with significantly larger velocities than this, the bright satellites of the Milky Way all appear to reside in subhaloes with maximum circular velocities below 40 kms. As recently highlighted by Boylan-Kolchin et al., this discrepancy implies that the majority of the most massive subhaloes within a cold dark matter galactic halo are much too concentrated to be consistent with the kinematic data for the bright Milky Way satellites. Here we show that no such discrepancy exists if haloes are made of warm, rather than cold dark matter because these haloes are less concentrated on account of their typically later formation epochs. Warm dark matter is one of several possible explanations for the observed kinematics of the satellites.

Galaxy and Mass Assembly (GAMA) is a project to study galaxy formation and evolution, combining imaging data from ultraviolet to radio with spectroscopic data from the AAOmega spectrograph on the Anglo-Australian Telescope. Using data from phase 1 of GAMA, taken over three observing seasons, and correcting for various minor sources of incompleteness, we calculate galaxy luminosity functions (LFs) and their evolution in the ugriz passbands. At low redshift, z < 0.1, we find that blue galaxies, defined according to a magnitude-dependent but non-evolving colour cut, are reasonably well fit over a range of more than ten magnitudes by simple Schechter functions in all bands. Red galaxies, and the combined blue-plus-red sample, require double power-law Schechter functions to fit a dip in their LF faintward of the characteristic magnitude M* before a steepening faint end. This upturn is at least partly due to dust-reddened disk galaxies. We measure evolution of the galaxy LF over the redshift range 0.002 < z < 0.5 using both a parametric fit and by measuring binned LFs in redshift slices. The characteristic luminosity L* is found to increase with redshift in all bands, with red galaxies showing stronger luminosity evolution than blue galaxies. The comoving number density of blue galaxies increases with redshift, while that of red galaxies decreases, consistent with prevailing movement from blue cloud to red sequence. As well as being more numerous at higher redshift, blue galaxies also dominate the overall luminosity density beyond redshifts z = 0.2. At lower redshifts, the luminosity density is dominated by red galaxies in the riz bands, by blue galaxies in u and g.

During the growth of a cold dark matter halo, the direction of its spin can undergo rapid changes. These could disrupt or even destroy a stellar disc forming in the halo, possibly resulting in the generation of a bulge or spheroid. We investigate the frequency of significant changes in the orientation of the angular momentum vector of dark matter haloes ("spin flips"), and their degree of correlation with mergers. We focus on haloes of mass similar to that of the Milky Way (MW) halo at redshift z=0 (log M/(Msol/h) = 12.0 to 12.5) and consider flips in the spin of the whole halo or just its inner parts. We find that a greater fraction of major mergers are associated with large spin flips than minor mergers. However, since major mergers are rare, the vast majority (93%) of large whole-halo spin flips (angle change >= 45 degrees) coincide with small mass changes, not major mergers. The spin vector of the inner halo experiences much more frequent flips than the halo as a whole. Over their entire lifetimes (i.e. after a halo acquires half of its final mass), over 10% of halos experience a flip of at least 45 degrees in the spin of the entire halo and nearly 60 percent experience a flip this large in the inner halo. These numbers are reduced to 9 percent for the whole halo and 47 percent for the inner halo when we consider only haloes with no major mergers after formation. Our analysis suggests that spin flips (whose effects are not currently included in galaxy formation models) could be an important factor in the morphological transformation of disc galaxies.

We present new deep observations of 'shell' structures in the halo of the nearby elliptical galaxy NGC 7600, alongside a movie of galaxy formation in a cold dark matter universe (available at http://www.virgo.dur.ac.uk/shell-galaxies). The movie, based on an ab initio cosmological simulation, shows how continuous accretion of clumps of dark matter and stars creates a swath of diffuse circumgalactic structures. The disruption of a massive clump on a near-radial orbit creates a complex system of transient concentric shells which bare a striking resemblance to those of NGC 7600. With the aid of the simulation we interpret NGC 7600 in the context of the CDM model.

We characterize substructure in the simulated stellar halos of Cooper et al. (2010) which were formed by the disruption of satellite galaxies within the cosmological N-body simulations of galactic halos of the Aquarius Project. These stellar halos exhibit a wealth of tidal features: broad overdensities and very narrow faint streams akin to those observed around the Milky Way. The substructures are distributed anisotropically on the sky, a characteristic that should become apparent in the next generation of photometric surveys. The normalized RMS of the density of stars on the sky appears to be systematically larger for our halos compared to the value estimated for the Milky Way from main sequence turn-off stars in the Sloan Digital Sky Survey. We show that this is likely to be due in part to contamination by faint QSOs and redder main sequence stars, and might suggest that ~10% of the Milky Way halo stars have formed in-situ.

If the dark matter consists of supersymmetric particles, Gamma-ray observatories such as the Large Area Telescope aboard the Fermi satellite may detect annihilation radiation from the haloes of galaxies and galaxy clusters. Much recent effort has been devoted to searching for this signal around the Milky Way's dwarf satellites. Using a new suite of high-resolution simulations of galaxy cluster haloes (the Phoenix Project), together with the Aquarius simulations of Milky-Way-like galaxy haloes, we show that higher signal-to-noise and equally clean signals are, in fact, predicted to come from nearby rich galaxy clusters. Most of the cluster emission is produced by small subhaloes with masses less than that of the Sun. The large range of mass scales covered by our two sets of simulations allows us to deduce a physically motivated extrapolation to these small (and unresolved) masses. Since tidal effects destroy subhaloes in the dense inner regions of haloes, most cluster emission is then predicted to come from large radii, implying that the nearest and brightest systems should be much more extended than Fermi's angular resolution limit. The most promising targets for detection are clusters such as Coma and Fornax, but detection algorithms must be tuned to the predicted profile of the emission if they are to maximize the chance of finding this weak signal.

We describe a correlation function statistic that quantifies the amount of spatial and kinematic substructure in the stellar halo. We test this statistic using model stellar halo realizations constructed from the Aquarius suite of six high-resolution N-body simulations in combination with the Galform semi-analytic galaxy formation model. These simulations show considerable scatter in the properties of stellar haloes. We find that our statistic can distinguish between these plausible alternatives for the global structure of the Milky Way stellar halo. We compare with observational data and show that pencil beam surveys of ~100 tracer stars (such as the Spaghetti Survey) are not sufficient to constrain the degree of structure in the Milky Way halo with this statistic. Larger area surveys with >1000 tracer stars (such as BHB stars in the Sloan Digital Sky Survey) provide much tighter constraints on comparisons between models and data. In our simulations, we find examples of haloes with spatial and kinematic substructure consistent with the available Milky Way data.

We use cosmological simulations from the Aquarius Project to study the orbital history of the Large Magellanic Cloud (LMC) and its potential association with other satellites of the Milky Way (MW). We search for dynamical analogs to the LMC and find a subhalo that matches the LMC position and velocity at either of its two most recent pericentric passages. This suggests that the LMC is not necessarily on its first approach to the MW, provided that the virial mass of the Milky Way is M_200 ~1.8e12 M_sun. The simulation results yield specific predictions for the position and velocity of systems associated with the LMC prior to infall. If on first approach, most should lie close to the LMC because the Galactic tidal field has not yet had enough time to disperse them. If on second approach, the list of potential associates increases substantially. Interestingly, our analysis rules out an LMC association for Draco and Ursa Minor, two of the dwarf spheroidals suggested by Lynden-Bell & Lynden-Bell to form part of the "Magellanic Ghostly Stream". Our results also indicate that the direction of the orbital angular momentum is a powerful test of LMC association. This test, however, requires precise proper motions, which are unavailable for most MW satellites. Of the 4 satellites with published proper motions, only the Small Magellanic Cloud is clearly associated with the LMC. Taken at face value, the proper motions of Carina, Fornax and Sculptor rule them out as potential associates, but this conclusion should be revisited when better data become available. The dearth of satellites clearly associated with the Clouds might be solved by wide-field imaging surveys that target its surroundings, a region that may prove a fertile hunting ground for faint, previously unnoticed MW satellites.

This paper describes the first catalogue of photometrically-derived stellar mass estimates for intermediate-redshift (z < 0.65) galaxies in the Galaxy And Mass Assembly (GAMA) spectroscopic redshift survey. These masses, as well as the full set of ancillary stellar population parameters, will be made public as part of GAMA data release 2. Although the GAMA database does include NIR photometry, we show that the quality of our stellar population synthesis fits is significantly poorer when these NIR data are included. Further, for a large fraction of galaxies, the stellar population parameters inferred from the optical-plus-NIR photometry are formally inconsistent with those inferred from the optical data alone. This may indicate problems in our stellar population library, or NIR data issues, or both; these issues will be addressed for future versions of the catalogue. For now, we have chosen to base our stellar mass estimates on optical photometry only. In light of our decision to ignore the available NIR data, we examine how well stellar mass can be constrained based on optical data alone. We use generic properties of stellar population synthesis models to demonstrate that restframe colour alone is in principle a very good estimator of stellar mass-to-light ratio, M*/Li. Further, we use the observed relation between restframe (g-i) and M*/Li for real GAMA galaxies to argue that, modulo uncertainties in the stellar evolution models themselves, (g-i) colour can in practice be used to estimate M*/Li to an accuracy of < ~0.1 dex. This 'empirically calibrated' (g-i)-M*/Li relation offers a simple and transparent means for estimating galaxies' stellar masses based on minimal data, and so provides a solid basis for other surveys to compare their results to z < ~0.4 measurements from GAMA.

We use the high resolution cosmological N-body simulations from the Aquarius project to investigate in detail the mechanisms that determine the shape of Milky Way-type dark matter haloes. We find that, when measured at the instantaneous virial radius, the shape of individual haloes changes with time, evolving from a typically prolate configuration at early stages to a more triaxial/oblate geometry at the present day. This evolution in halo shape correlates well with the distribution of the infalling material: prolate configurations arise when haloes are fed through narrow filaments, which characterizes the early epochs of halo assembly, whereas triaxial/oblate configurations result as the accretion turns more isotropic at later times. Interestingly, at redshift z=0, clear imprints of the past history of each halo are recorded in their shapes at different radii, which also exhibit a variation from prolate in the inner regions to triaxial/oblate in the outskirts. Provided that the Aquarius haloes are fair representatives of Milky Way-like 10^12 Msun objects, we conclude that the shape of such dark matter haloes is a complex, time-dependent property, with each radial shell retaining memory of the conditions at the time of collapse.

We present a model for the satellites of the Milky Way in which galaxy formation is followed using semi-analytic techniques applied to the six high-resolution N-body simulations of galactic halos of the Aquarius project. The model, calculated using the Galform code, incorporates improved treatments of the relevant physics in the LambdaCDM cosmogony, particularly a self-consistent calculation of reionization by UV photons emitted by the forming galaxy population, including the progenitors of the central galaxy. Along the merger tree of each halo, the model calculates gas cooling (by Compton scattering off cosmic microwave background photons, molecular hydrogen and atomic processes), gas heating (from hydrogen photoionization and supernova energy), star formation and evolution. The evolution of the intergalactic medium is followed simultaneously with that of the galaxies. Star formation in the more massive progenitor subhalos is suppressed primarily by supernova feedback, while for smaller subhalos it is suppressed primarily by photoionization due to external and internal sources. The model is constrained to match a wide range of properties of the present day galaxy population as a whole, but at high redshift it requires an escape fraction of UV photons near unity in order completely to reionize the universe by redshift z ~ 8. In the most successful model the local sources photoionize the pre-galactic region completely by z ~ 10. In addition to the luminosity function of Milky Way satellites, the model matches their observed luminosity-metallicity relation, their radial distribution and the inferred values of the mass within 300 pc, which in the models increase slowly but significantly with luminosity. There is a large variation in satellite properties from halo to halo, with the luminosity function, for example, varying by a factor of ~ 2 among the six simulations.

We apply a basis function expansion method to create a time-evolving density/potential approximation of the late growth of simulated N-body dark matter haloes. We demonstrate how the potential of a halo from the Aquarius Project can be accurately represented by a small number of basis functions, and show that the halo expansion (HEX) method provides a way to replay simulations. We explore the level of accuracy of the technique as well as some of its limitations. We find that the number of terms included in the expansion must be large enough to resolve the large-scale distribution and shape of the halo but, beyond this, additional terms result in little further improvement. Particle and subhalo orbits can be integrated in this realistic, time-varying halo potential approximation, at much lower cost than the original simulation, with high fidelity for many individual orbits, and a good match to the distributions of orbital energy and angular momentum. Statistically, the evolution of structural subhalo properties, such as mass, half-mass radius and characteristic circular velocity, are very well reproduced in the halo expansion approximation over several gigayears. We demonstrate an application of the technique by following the evolution of an orbiting subhalo at much higher resolution than can be achieved in the original simulation. Our method represents a significant improvement over commonly used techniques based on static analytical descriptions of the halo potential.

We investigate the properties of satellite galaxies that surround isolated hosts within the redshift range 0.01 < z < 0.15, using data taken as part of the Galaxy And Mass Assembly survey. Making use of isolation and satellite criteria that take into account stellar mass estimates, we find 3514 isolated galaxies of which 1426 host a total of 2998 satellites. Separating the red and blue populations of satellites and hosts, using colour-mass diagrams, we investigate the radial distribution of satellite galaxies and determine how the red fraction of satellites varies as a function of satellite mass, host mass and the projected distance from their host. Comparing the red fraction of satellites to a control sample of small neighbours at greater projected radii, we show that the increase in red fraction is primarily a function of host mass. The satellite red fraction is about 0.2 higher than the control sample for hosts with 11.0 < log M_* < 11.5, while the red fractions show no difference for hosts with 10.0 < log M_* < 10.5. For the satellites of more massive hosts the red fraction also increases as a function of decreasing projected distance. Our results suggest that the likely main mechanism for the quenching of star formation in satellites hosted by isolated galaxies is strangulation.

We investigate the formation and evolution of satellite galaxies using smoothed particle hydrodynamics (SPH) simulations of a Milky Way(MW)-like system, focussing on the best resolved examples, analogous to the classical MW satellites. Comparing with a pure dark matter simulation, we find that the condensation of baryons has had a relatively minor effect on the structure of the satellites' dark matter halos. The stellar mass that forms in each satellite agrees relatively well over three levels of resolution (a factor of ~64 in particle mass) and scales with (sub)halo mass in a similar way in an independent semi-analytical model. Our model provides a relatively good match to the average luminosity function of the MW and M31. To establish whether the potential wells of our satellites are realistic, we measure their masses within observationally determined half-light radii, finding that the most massive examples have somewhat higher mass-to-light ratios than those derived for the MW dSphs from stellar kinematic data. A statistical test yields a ~9 percent probability that the simulated and observationally derived distributions of masses are consistent. Our results may suggest that either the MW halo is less massive than assumed in our simulations (~1.4e12 M_sun) or that there is substantial scatter in the satellite luminosity function or distribution of mass-to-light ratios at fixed host halo mass. Alternatively, feedback processes not properly captured by our simulations may have reduced the central densities of (sub)halos, or the subhalos may have initially formed with lower concentrations as would be the case, for example, if the dark matter were made of warm, rather than cold particles.

Using the complete GAMA-I survey covering ~142 sq. deg. to r=19.4, of which ~47 sq. deg. is to r=19.8, we create the GAMA-I galaxy group catalogue (G3Cv1), generated using a friends-of-friends (FoF) based grouping algorithm. Our algorithm has been tested extensively on one family of mock GAMA lightcones, constructed from Lambda-CDM N-body simulations populated with semi-analytic galaxies. Recovered group properties are robust to the effects of interlopers and are median unbiased in the most important respects. G3Cv1 contains 14,388 galaxy groups (with multiplicity >= 2$), including 44,186 galaxies out of a possible 110,192 galaxies, implying ~40% of all galaxies are assigned to a group. The similarities of the mock group catalogues and G3Cv1 are multiple: global characteristics are in general well recovered. However, we do find a noticeable deficit in the number of high multiplicity groups in GAMA compared to the mocks. Additionally, despite exceptionally good local spatial completeness, G3Cv1 contains significantly fewer compact groups with 5 or more members, this effect becoming most evident for high multiplicity systems. These two differences are most likely due to limitations in the physics included of the current GAMA lightcone mock. Further studies using a variety of galaxy formation models are required to confirm their exact origin.

(Abridged) We study the phase-space distribution of satellite galaxies associated with late-type galaxies in the GIMIC suite of simulations. GIMIC consists of re-simulations of 5 cosmologically representative regions from the Millennium simulation, which have higher resolution and incorporate baryonic physics. Whilst the disc of the galaxy is well aligned with the inner regions (r ~ 0.1r_200) of the dark matter halo, both in shape and angular momentum, there can be substantial misalignments at larger radii (r ~ r_200). Misalignments of > 45 deg are seen in ~ 30% of our sample. We find that the satellite population aligns with the shape (and angular momentum) of the outer dark matter halo. However, the alignment with the galaxy is weak owing to the mismatch between the disc and dark matter halo. Roughly 20% of the satellite systems with ten bright galaxies within r_200 exhibit a polar spatial alignment with respect to the galaxy --- an orientation reminiscent of the classical satellites of the Milky Way. We find that a small fraction (~10%) of satellite systems show evidence for rotational support which we attribute to group infall. There is a bias towards satellites on prograde orbits relative to the spin of the dark matter halo (and to a lesser extent with the angular momentum of the disc). This preference towards co-rotation is stronger in the inner regions of the halo where the most massive satellites accreted at relatively early times are located. We attribute the anisotropic spatial distribution and angular momentum bias of the satellites at z=0 to their directional accretion along the major axes of the dark matter halo. The satellite galaxies have been accreted relatively recently compared to the dark matter mass and have experienced less phase-mixing and relaxation --- the memory of their accretion history can remain intact to z=0.

BigBOSS is a Stage IV ground-based dark energy experiment to study baryon acoustic oscillations (BAO) and the growth of structure with a wide-area galaxy and quasar redshift survey over 14,000 square degrees. It has been conditionally accepted by NOAO in response to a call for major new instrumentation and a high-impact science program for the 4-m Mayall telescope at Kitt Peak. The BigBOSS instrument is a robotically-actuated, fiber-fed spectrograph capable of taking 5000 simultaneous spectra over a wavelength range from 340 nm to 1060 nm, with a resolution R = 3000-4800. Using data from imaging surveys that are already underway, spectroscopic targets are selected that trace the underlying dark matter distribution. In particular, targets include luminous red galaxies (LRGs) up to z = 1.0, extending the BOSS LRG survey in both redshift and survey area. To probe the universe out to even higher redshift, BigBOSS will target bright [OII] emission line galaxies (ELGs) up to z = 1.7. In total, 20 million galaxy redshifts are obtained to measure the BAO feature, trace the matter power spectrum at smaller scales, and detect redshift space distortions. BigBOSS will provide additional constraints on early dark energy and on the curvature of the universe by measuring the Ly-alpha forest in the spectra of over 600,000 2.2 < z < 3.5 quasars. BigBOSS galaxy BAO measurements combined with an analysis of the broadband power, including the Ly-alpha forest in BigBOSS quasar spectra, achieves a FOM of 395 with Planck plus Stage III priors. This FOM is based on conservative assumptions for the analysis of broad band power (kmax = 0.15), and could grow to over 600 if current work allows us to push the analysis to higher wave numbers (kmax = 0.3). BigBOSS will also place constraints on theories of modified gravity and inflation, and will measure the sum of neutrino masses to 0.024 eV accuracy.

We study the luminosity function of satellite galaxies around isolated primaries using the Sloan Digital Sky Survey (SDSS) spectroscopic and photometric galaxy samples. We select isolated primaries from the spectroscopic sample and search for potential satellites in the much deeper photometric sample. For primaries of similar luminosity to the Milky Way and M31, we are able to stack as many as $\sim\negthinspace 20,000$ galaxy systems to obtain robust statistical results. We derive the satellite luminosity function extending almost 8 magnitudes fainter than the primary galaxy. We also determine how the satellite luminosity function varies with the luminosity, colour and concentration of the primary. We find that, in the mean, isolated primaries of comparable luminosity to the Milky Way and M31 contain about a factor of two fewer satellites brighter than $M_V=-14$ than the average of the Milky Way and M31.

We present an algorithm for quickly generating multiple realizations of N-body simulations to be used, for example, for cosmological parameter estimation from surveys of large-scale structure. Our algorithm uses a new method to resample the large-scale (Gaussian-distributed) Fourier modes in a periodic N-body simulation box in a manner that properly accounts for the nonlinear mode-coupling between large and small scales. We find that our method for adding new large-scale mode realizations recovers the nonlinear power spectrum to sub-percent accuracy on scales larger than about half the Nyquist frequency of the simulation box. Using 20 N-body simulations, we obtain a power spectrum covariance matrix estimate that matches the estimator in Takahashi et al. (2009) (from 5000 simulations) with < 20% errors in all matrix elements. Comparing the rates of convergence, we determine that our algorithm requires ~ 8 times fewer simulations to achieve a given error tolerance in estimates of the power spectrum covariance matrix. The degree of success of our algorithm indicates that we understand the main physical processes that give rise to the correlations in the matter power spectrum. Namely, the large-scale Fourier modes modulate both the degree of structure growth through the variation in the effective local matter density and also the spatial frequency of small-scale perturbations through large-scale displacements. We expect our algorithm to be useful for noise modeling when constraining cosmological parameters from weak lensing (cosmic shear) and galaxy surveys, rescaling summary statistics of N-body simulations for new cosmological parameter values, and any applications where the influence of Fourier modes larger than the simulation size must be accounted for.

We predict the formation histories, properties and descendants of Lyman-break galaxies (LBGs) in the Lambda-CDM cosmology. In our model, which incorporates a top-heavy IMF in starbursts, we find that most LBGs are starbursts triggered by minor mergers of galaxies. The duration of the LBG phase is predicted to be quite short, ~20-60 Myr. We investigate the distributions of stellar and halo masses and morphologies for bright (L_UV > L*_UV) and faint (L_UV > 0.1 L*_UV) LBGs at z=3, z=6 and z=10 (where we classify LBGs according to their rest-frame UV luminosities relative the observed characteristic luminosity L*_UV at z \approx 3). Bright LBGs at z=3 are predicted to have median stellar masses ~ 1x10^9 Msun/h and host halo masses ~ 3x10^{11} Msun/h, and to be typically mildly disk-dominated in stellar mass. On the other hand, faint LBGs at z=10 are predicted to have median stellar masses of only ~ 1x10^7 Msun/h and host halo masses 2x10^{10} Msun/h, and to be generally bulge-dominated. Bright LBGs at z=3 evolve to present-day galaxies with median stellar mass ~ 5x^{10} Msun/h (comparable to the Milky Way), consisting of roughly equal numbers of disk- and bulge-dominated systems, and hosted by halos with median mass ~2x10^{13} Msun/h (corresponding to medium-size galaxy groups). The model predicts that 40% of Milky Way mass galaxies at the present-day have a bright LBG progenitor in the redshift range 3<z<4, while 95% have a faint LBG progenitor in the same redshift range, and 7% have a faint LBG progenitor at 10<z<11. With our multiwavelength model, we also investigate the overlap between the LBG population and that of submillimetre selected galaxies (SMGs); at z=3, only ~1% of bright LBGs are also predicted to also be bright SMGs (with an 850 mum flux in excess of 5 mJy).

We calculate the orbital angular momentum of dark matter subhaloes in the Aquarius simulations of cold dark matter galactic haloes. We calculate the orientation of their angular momentum relative to that of the spin vector of their host halo and find a variety of different configurations. All six Aquarius haloes contain statistically significant populations of subhalo orbits that are aligned with the main halo spin. All haloes posses a population of subhaloes that rotates in the same direction as the main halo and three of them possess, in addition, a population that rotates in the opposite direction. These configurations arise from the filamentary accretion of subhaloes. Quasi-planar distributions of coherently rotating satellites, such as those inferred in the Milky Way and other galaxies, arise naturally in simulations of a $\Lambda$CDM universe.

The stellar initial mass function (IMF) describes the distribution in stellar masses produced from a burst of star formation. For more than fifty years, the implicit assumption underpinning most areas of research involving the IMF has been that it is universal, regardless of time and environment. We measure the high-mass IMF slope for a sample of low-to-moderate redshift galaxies from the Galaxy And Mass Assembly survey. The large range in luminosities and galaxy masses of the sample permits the exploration of underlying IMF dependencies. A strong IMF-star formation rate dependency is discovered, which shows that highly star forming galaxies form proportionally more massive stars (they have IMFs with flatter power-law slopes) than galaxies with low star formation rates. This has a significant impact on a wide variety of galaxy evolution studies, all of which rely on assumptions about the slope of the IMF. Our result is supported by, and provides an explanation for, the results of numerous recent explorations suggesting a variation of or evolution in the IMF.

We use multiwavelength data from the Galaxy And Mass Assembly (GAMA) and Herschel ATLAS (H-ATLAS) surveys to compare the relationship between various dust obscuration measures in galaxies. We explore the connections between the ultraviolet (UV) spectral slope, $\beta$, the Balmer decrement, and the far infrared (IR) to $150\,$nm far ultraviolet (FUV) luminosity ratio. We explore trends with galaxy mass, star formation rate (SFR) and redshift in order to identify possible systematics in these various measures. We reiterate the finding of other authors that there is a large scatter between the Balmer decrement and the $\beta$ parameter, and that $\beta$ may be poorly constrained when derived from only two broad passbands in the UV. We also emphasise that FUV derived SFRs, corrected for dust obscuration using $\beta$, will be overestimated unless a modified relation between $\beta$ and the attenuation factor is used. Even in the optimum case, the resulting SFRs have a significant scatter, well over an order of magnitude. While there is a stronger correlation between the IR to FUV luminosity ratio and $\beta$ parameter than with the Balmer decrement, neither of these correlations are particularly tight, and dust corrections based on $\beta$ for high redshift galaxy SFRs must be treated with caution. We conclude with a description of the extent to which the different obscuration measures are consistent with each other as well as the effects of including other galactic properties on these correlations.

Cosmological N-body simulations indicate that the spherically-averaged density profiles of cold dark matter halos are accurately described by Einasto profiles, where the logarithmic slope is a power-law of adjustable exponent, \gamma =dln\rho /dlnr ~ r^\alpha $. The pseudo-phase-space density (PPSD) profiles of CDM halos also show remarkable regularity, and are well approximated by simple power laws, Q(r)=\rho /\sigma ^3 ~ r^-\chi . We show that this is expected from dynamical equilibrium considerations, since Jeans' equations predict that the pseudo-phase-space density profiles of Einasto halos should resemble power laws over a wide range of radii. For the values of \alpha typical of CDM halos, the inner Q profiles of equilibrium halos deviate significantly from a power law only very close to the center, and simulations of extremely high-resolution would be needed to detect such deviations unambiguously. We use an ensemble of halos drawn from the Millennium-II simulation to study which of these two alternatives describe best the mass profile of CDM halos. Our analysis indicates that at the resolution of the best available simulations, both Einasto and power-law PPSD profiles (with adjustable exponents \alpha and \chi, respectively) provide equally acceptable fits to the simulations. A full account of the structure of CDM halos requires understanding how the shape parameters that characterize departures from self-similarity, like \alpha or \chi, are determined by evolutionary history, environment or initial conditions.

We describe and apply a simple prescription for defining connected structures in galaxy redshift surveys. The method is based upon two passes with a friends-of-friends groupfinder. The first pass uses a cylindrical linking volume to find galaxy groups and clusters, in order to suppress the line-of-sight smearing introduced by the large random velocities of galaxies within these deep potential wells. The second pass, performed with a spherical linking volume, identifies the connected components. This algorithm has been applied to the 2dFGRS, within which it picks out a total of 7,603 systems containing at least two galaxies and having a mean redshift less than 0.12. Connected systems with many members appear filamentary in nature, and the algorithm recovers two particularly large filaments within the 2dFGRS. For comparison, the algorithm is has also been applied to LambdaCDM mock galaxy surveys. While the model population of such systems is broadly similar to that in the 2dFGRS, it does not generally contain such extremely large structures.

We present an analysis of the properties of the lowest Halpha-luminosity galaxies (L_Halpha<4x10^32 W; SFR<0.02 Msun/yr) in the Galaxy And Mass Assembly (GAMA) survey. These galaxies make up the the rise above a Schechter function in the number density of systems seen at the faint end of the Halpha luminosity function. Above our flux limit we find that these galaxies are principally composed of intrinsically low stellar mass systems (median stellar mass =2.5x10^8 Msun) with only 5/90 having stellar masses M>10^10 Msun. The low SFR systems are found to exist predominantly in the lowest density environments (median density ~0.02 galaxy Mpc^-2 with none in environments more dense than ~1.5 galaxy Mpc^-2). Their current specific star formation rates (SSFR; -8.5 < log(SSFR[yr^-1])<-12.) are consistent with their having had a variety of star formation histories. The low density environments of these galaxies demonstrates that such low-mass, star-forming systems can only remain as low-mass and forming stars if they reside sufficiently far from other galaxies to avoid being accreted, dispersed through tidal effects or having their gas reservoirs rendered ineffective through external processes.

We study the role of submillimetre galaxies (SMGs) in the galaxy formation process in the Lambda Cold Dark Matter cosmology. We use the Baugh et al. (2005) semi-analytical model, which matches the observed SMG number counts and redshift distribution by assuming a top-heavy initial mass function (IMF) in bursts triggered by galaxy mergers. We build galaxy merger trees and follow the evolution and properties of SMGs and their descendants. Our primary sample of model SMGs consists of galaxies which had 850 mu fluxes brighter than 5 mJy at some redshift z>1. Our model predicts that the present-day descendants of such SMGs cover a wide range of stellar masses ~ 10^{10} - 10^{12} Msun/h, with a median ~ 10^{11} Msun/h, and that more than 70% of these descendants are bulge-dominated. More than 50% of present day galaxies with stellar masses larger than 7 x 10^{11} Msun/h are predicted to be descendants of such SMGs. We find that although SMGs make an important contribution to the total star formation rate at z~2, the final stellar mass produced in the submillimetre phase contributes only 0.2% of the total present-day stellar mass, and 2% of the stellar mass of SMG descendants, in stark contrast to the popular picture in which the SMG phase marks the production of the bulk of the mass of present day massive ellipticals.

Star formation rate and accummulated stellar mass are two fundamental physical quantities that describe the evolutionary state of a forming galaxy. Two recent attempts to determine the relationship between these quantities, by interpreting a sample of star-forming galaxies at redshift of z~4, have led to opposite conclusions. We use a model galaxy population to investigate possible causes for this discrepancy and conclude that minor errors in the conversion from observables to physical quantities can lead to major misrepresentation when applied without awareness of sample selection. We also investigate, in a general way, the physical origin of the correlation between star formation rate and stellar mass within hierarchical galaxy formation theory.

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

We make a detailed investigation of the properties of Lyman-break galaxies (LBGs) in the LambdaCDM model. We present predictions for two published variants of the GALFORM semi-analytical model: the Baugh et al. (2005) model, which has star formation at high redshifts dominated by merger-driven starbursts with a top-heavy IMF, and the Bower et al. (2006) model, which has AGN feedback and a standard Solar neighbourhood IMF throughout. We show predictions for the evolution of the rest-frame far-UV luminosity function in the redshift range z=3-20, and compare with the observed luminosity functions of LBGs at z=3-10. We find that the Baugh et al. model is in excellent agreement with these observations, while the Bower et al. model predicts too many high-luminosity LBGs. Dust extinction, which is predicted self-consistently based on galaxy gas contents, metallicities and sizes, is found to have a large effect on LBG luminosities. We compare predictions for the size evolution of LBGs at different luminosities with observational data for 2<z<7, and find the Baugh et al. model to be in good agreement. We present predictions for stellar, halo and gas masses, star formation rates, circular velocities, bulge-to-disk ratios, gas and stellar metallicities and clustering bias, as functions of far-UV luminosity and redshift. We find broad consistency with current observational constraints. We then present predictions for the abundance and angular sizes of LBGs out to very high redshift (z<20), finding that planned deep surveys with JWST should detect objects out to z<15. The typical UV luminosities of galaxies are predicted to be very low at high redshifts, which has implications for detecting the galaxies responsible for reionizing the IGM; for example, at z=10, 50% of the ionizing photons are expected to be produced by galaxies fainter than M_AB(1500A)-5logh ~ -15.

We measure the projected cross-correlation between low redshift (z < 0.5) far-IR selected galaxies in the SDP field of the Herschel-ATLAS (H-ATLAS) survey and optically selected galaxies from the Galaxy and Mass Assembly (GAMA) redshift survey. In order to obtain robust correlation functions, we restrict the analysis to a subset of 969 out of 6900 H-ATLAS galaxies, which have reliable optical counterparts with r<19.4 mag and well-determined spectroscopic redshifts. The overlap region between the two surveys is 12.6 sq. deg; the matched sample has a median redshift of z ~ 0.2. The cross-correlation of GAMA and H-ATLAS galaxies within this region can be fitted by a power law, with correlation length r_0 ~ 4.63 +/- 0.51 Mpc. Comparing with the corresponding auto-correlation function of GAMA galaxies within the SDP field yields a relative bias (averaged over 2-8 Mpc) of H-ATLAS and GAMA galaxies of b_H/b_G ~ 0.6. Combined with clustering measurements from previous optical studies, this indicates that most of the low redshift H-ATLAS sources are hosted by halos with masses comparable to that of the Milky Way. The correlation function appears to depend on the 250 um luminosity, L_250, with bright (median luminosity \nu L_250 ~ 1.6 x 10^10 L_sun) objects being somewhat more strongly clustered than faint (\nu L_250 ~ 4.0 x 10^9 L_sun) objects. This implies that galaxies with higher dust-obscured star formation rates are hosted by more massive halos.

The existence of hot, X-ray luminous gaseous coronae surrounding present day L* galaxies is a generic prediction of galaxy formation theory in the cold dark matter cosmogony. While extended X-ray emission has been known to exist around elliptical galaxies for a long time, diffuse extra-planar emission has only recently been detected around disc galaxies. We compile samples of elliptical and disc galaxies that have Chandra and XMM-Newton measurements, and compare the scaling of the coronal X-ray luminosity (L_X) with both the K-band luminosity (L_K) and the coronal X-ray temperature (T_X). The X-ray flux measurements are corrected for non-thermal point source contamination by spatial excision and spectral subtraction for resolved and unresolved sources respectively. We find that the properties of the extended X-ray emission from galaxies of different morphological types are similar: for both elliptical and disc galaxies, the L_X - L_K and L_X - T_X relations have similar slope, normalisation and scatter. The observed universality of coronal X-ray properties suggests that the bulk of this emission originates from gas that has been accreted, shock-heated and compressed during the assembly of the galaxy and that outflows triggered by stellar processes make only a minor contribution to the X-ray emission. This reservoir of cooling gas is a potential source of morphological transformation; it provides a fresh supply of material for discs to grow around galaxies of all morphological types.

The spectral energy distribution of galaxies is a complex function of the star formation history and geometrical arrangement of stars and gas in galaxies. The computation of the radiative transfer of stellar radiation through the dust distribution is time-consuming. This aspect becomes unacceptable in particular when dealing with the predictions by semi-analytical galaxy formation models populating cosmological volumes, to be then compared with multi-wavelength surveys. Mainly for this aim, we have implemented an artificial neural network algorithm into the spectro-photometric and radiative transfer code GRASIL in order to compute the spectral energy distribution of galaxies in a short computing time. This allows to avoid the adoption of empirical templates that may have nothing to do with the mock galaxies output by models. The ANN has been implemented to compute the dust emission spectrum (the bottleneck of the computation), and separately for the star-forming molecular clouds and the diffuse dust (due to their different properties and dependencies). We have defined the input neurons effectively determining their emission, which means this implementation has a general applicability and is not linked to a particular galaxy formation model. We have trained the net for the disc and spherical geometries, and tested its performance to reproduce the SED of disc and starburst galaxies, as well as for a semi-analytical model for spheroidal galaxies. We have checked that for this model both the SEDs and the galaxy counts in the Herschel bands obtained with the ANN approximation are almost superimposed to the same quantities obtained with the full GRASIL. We conclude that this method appears robust and advantageous, and will present the application to a more complex SAM in another paper.

In order to generate credible 0.1-2 {\mu}m SEDs, the GAMA project requires many Gigabytes of imaging data from a number of instruments to be re-processed into a standard format. In this paper we discuss the software infrastructure we use, and create self-consistent ugrizYJHK photometry for all sources within the GAMA sample. Using UKIDSS and SDSS archive data, we outline the pre-processing necessary to standardise all images to a common zeropoint, the steps taken to correct for seeing bias across the dataset, and the creation of Gigapixel-scale mosaics of the three 4x12 deg GAMA regions in each filter. From these mosaics, we extract source catalogues for the GAMA regions using elliptical Kron and Petrosian matched apertures. We also calculate S\'ersic magnitudes for all galaxies within the GAMA sample using SIGMA, a galaxy component modelling wrapper for GALFIT 3. We compare the resultant photometry directly, and also calculate the r band galaxy LF for all photometric datasets to highlight the uncertainty introduced by the photometric method. We find that (1) Changing the object detection threshold has a minor effect on the best-fitting Schechter parameters of the overall population (M* +/- 0.055mag, {\alpha} +/- 0.014, {\Phi}* +/- 0.0005 h^3 Mpc^{-3}). (2) An offset between datasets that use Kron or Petrosian photometry regardless of the filter. (3) The decision to use circular or elliptical apertures causes an offset in M* of 0.20mag. (4) The best-fitting Schechter parameters from total-magnitude photometric systems (such as SDSS modelmag or S\'ersic magnitudes) have a steeper faint-end slope than photometry dependent on Kron or Petrosian magnitudes. (5) Our Universe's total luminosity density, when calculated using Kron or Petrosian r-band photometry, is underestimated by at least 15%.

We present self-consistent star formation rates derived through pan-spectral analysis of galaxies drawn from the Galaxy and Mass Assembly (GAMA) survey. We determine the most appropriate form of dust obscuration correction via application of a range of extinction laws drawn from the literature as applied to Halpha, [O{II}] and UV luminosities. These corrections are applied to a sample of 31,508 galaxies from the GAMA survey at z < 0.35. We consider several different obscuration curves, including those of Milky Way, Calzetti (2001) and Fischera and Dopita (2005) curves and their effects on the observed luminosities. At the core of this technique is the observed Balmer decrement, and we provide a prescription to apply optimal obscuration corrections using the Balmer decrement. We carry out an analysis of the star formation history (SFH) using stellar population synthesis tools to investigate the evolutionary history of our sample of galaxies as well as to understand the effects of variation in the Initial Mass Function (IMF) and the effects this has on the evolutionary history of galaxies. We find that the Fischera and Dopita (2005) obscuration curve with an R_{v} value of 4.5 gives the best agreement between the different SFR indicators. The 2200A feature needed to be removed from this curve to obtain complete consistency between all SFR indicators suggesting that this feature may not be common in the average integrated attenuation of galaxy emission. We also find that the UV dust obscuration is strongly dependent on the SFR.

The Galaxy And Mass Assembly (GAMA) survey has been operating since February 2008 on the 3.9-m Anglo-Australian Telescope using the AAOmega fibre-fed spectrograph facility to acquire spectra with a resolution of R~1300 for 120,862 SDSS selected galaxies. The target catalogue constitutes three contiguous equatorial regions centred at 9h (G09), 12h (G12) and 14.5h (G15) each of 12 x 4 sq.deg to limiting fluxes of r < 19.4, r < 19.8, and r < 19.4 mag respectively (and additional limits at other wavelengths). Spectra and reliable redshifts have been acquired for over 98 per cent of the galaxies within these limits. Here we present the survey footprint, progression, data reduction, redshifting, re-redshifting, an assessment of data quality after 3 years, additional image analysis products (including ugrizYJHK photometry, Sersic profiles and photometric redshifts), observing mask, and construction of our core survey catalogue (GamaCore). From this we create three science ready catalogues: GamaCoreDR1 for public release, which includes data acquired during year 1 of operations within specified magnitude limits (February 2008 to April 2008); GamaCoreMainSurvey containing all data above our survey limits for use by the GAMA team and collaborators; and GamaCoreAtlasSv containing year 1, 2 and 3 data matched to Herschel-ATLAS Science Demonstration data. These catalogues along with the associated spectra, stamps and profiles can be accessed via the GAMA website: http://www.gama-survey.org/

We use the Aquarius simulation series to study the imprint of assembly history on the structure of Galaxy-mass cold dark matter halos. Our results confirm earlier work regarding the influence of mergers on the mass density profile and the inside-out growth of halos. The inner regions that contain the visible galaxies are stable since early times and are significantly affected only by major mergers. Particles accreted diffusely or in minor mergers are found predominantly in the outskirts of halos. Our analysis reveals trends that run counter to current perceptions of hierarchical halo assembly. For example, major mergers (i.e. those with progenitor mass ratios greater than 1:10) contribute little to the total mass growth of a halo, on average less than 20 per cent for our six Aquarius halos. The bulk is contributed roughly equally by minor mergers and by "diffuse" material which is not resolved into individual objects. This is consistent with modeling based on excursion-set theory which suggests that about half of this diffuse material should not be part of a halo of any scale. Interestingly, the simulations themselves suggest that a significantly fraction is not truly diffuse, since it was ejected from earlier halos by mergers prior to their joining the main system. The Aquarius simulations resolve halos to much lower mass scales than are expected to retain gas or form stars. These results thus confirm that most of the baryons from which visible galaxies form are accreted diffusely, rather than through mergers, and they suggest that only relatively rare major mergers will affect galaxy structure at later times.

We study the population statistics of the surviving subhaloes of LCDM dark matter haloes using a set of very high resolution N-body simulations. These include both simulations of representative regions of the Universe and ultra-high resolution resimulations of individual dark matter haloes. We find that more massive haloes tend to have a larger mass fraction in subhaloes. For example, cluster size haloes typically have 7.5 percent of their mass within R200 in substructures of fractional mass larger than 1e-5, which is 25 percent higher than galactic haloes. There is, however, a large variance in the subhalo mass fraction from halo to halo, whereas the subhalo abundance shows much higher regularity. For dark matter haloes of fixed mass, the subhalo abundance decreases by 30 percent between redshift 2 and 0. The subhalo abundance function correlates with the host halo concentration parameter and formation redshift. However, the intrinsic scatter is not significantly reduced for narrow ranges of concentration parameter or formation redshift, showing that they are not the dominant parameters that determine the subhalo abundance in a halo.

We track the co-evolution of supermassive black holes (SMBHs) and their host galaxies. The calculation is embedded in the GALFORM semi-analytical model which simulates the formation and evolution of galaxies in a cold dark matter (CDM) universe. During the evolution of the host galaxy, hot and cold gas are added to the SMBH by flows triggered by halo gas cooling, disc instabilities and galaxy mergers. This builds up the mass and spin of the BH, and the resulting accretion power regulates the gas cooling and subsequent star formation. The accretion flow is assumed to form a geometrically thin cool disc when the accretion rate exceeds 0.01\dot{M}_Edd, and a geometrically thick, radiatively inefficient hot flow when the accretion rate falls below this value. The resulting quasar optical luminosity function matches observations very well, and the mass of the SMBH correlates with the mass of the galaxy bulge as observed. The BH spin distribution depends strongly on whether the gas in any given accretion episode remains in the same plane (prolonged accretion) or whether, due to self-gravity, it fragments into multiple, randomly aligned accretion episodes (chaotic accretion). In the chaotic accretion model there is a clear correlation of spin with SMBH mass. Massive BHs (M>5\times10^8\Msun) are hosted by giant elliptical galaxies and are rapidly spinning, while lower mass BHs are hosted in spiral galaxies and have much lower spin. Using the Blandford-Znajek mechanism for jet production to calculate the jet power, our model is able to reproduce the radio loudness of radio galaxies, LINERS and Seyferts. This is the first confirmation that a CDM galaxy formation model can reproduce the observed phenomenology of AGN.

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

Lensing flux-ratio anomalies have been frequently observed and taken as evidence for the presence of abundant dark matter substructures in lensing galaxies, as predicted by the cold dark matter (CDM) model of cosmogony. In previous work, we examined the cusp-caustic relations of the multiple images of background quasars lensed by galaxy-scale dark matter haloes, using a suite of high-resolution N-body simulations (the Aquarius simulations). In this work, we extend our previous calculations to incorporate both the baryonic and diffuse dark components in lensing haloes. We include in each lensing simulation: (1) a satellite galaxy population derived from a semi-analytic model applied to the Aquarius haloes, (2) an empirical Milky-Way globular cluster population and (3) satellite streams (diffuse dark component) identified in the simulations. Accounting for these extra components, we confirm our earlier conclusion that the abundance of intrinsic substructures (dark or bright, bound or diffuse) is not sufficient to explain the observed frequency of cusp-caustic violations in the CLASS survey. We conclude that the observed effect could be the result of the small number statistics of CLASS, or intergalactic haloes along the line of sight acting as additional sources of lensing flux anomalies. Another possibility is that this discrepancy signals a failure of the CDM model.

Semi-analytic models are a powerful tool for studying the formation of galaxies. However, these models inevitably involve a significant number of poorly constrained parameters that must be adjusted to provide an acceptable match to the observed universe. In this paper, we set out to quantify the degree to which observational data-sets can constrain the model parameters. By revealing degeneracies in the parameter space we can hope to better understand the key physical processes probed by the data. We use novel mathematical techniques to explore the parameter space of the GALFORM semi-analytic model. We base our investigation on the Bower et al. 2006 version of GALFORM, adopting the same methodology of selecting model parameters based on an acceptable match to the local bJ and K luminosity functions. The model contains 16 parameters that are poorly constrained, and we investigate this parameter space using the Model Emulator technique, constructing a Bayesian approximation to the GALFORM model that can be rapidly evaluated at any point in parameter space. By combining successive waves of emulation, we show that only 0.26% of the initial volume is of interest for further exploration. However, within this region we show that the Bower et al. 2006 model is only one choice from an extended sub-space of model parameters that can provide equally acceptable fits. We explore the geometry of this region and begin to explore the physical connections between parameters that are exposed by this analysis. We also consider the impact of adding additional observational data to further constrain the parameter space.

The existence of X-ray luminous gaseous coronae around massive disc galaxies is a long-standing prediction of galaxy formation theory in the cold dark matter cosmogony. This prediction has garnered little observational support, with non-detections commonplace and detections for only a relatively small number of galaxies which are much less luminous than expected. We investigate the coronal properties of a large sample of bright, disc-dominated galaxies extracted from the GIMIC suite of cosmological hydrodynamic simulations recently presented by Crain et al. Remarkably, the simulations reproduce the observed scalings of X-ray luminosity with K-band luminosity and star formation rate and, when account is taken of the density structure of the halo, with disc rotation velocity as well. Most of the star formation in the simulated galaxies (which have realistic stellar mass fractions) is fuelled by gas cooling from a quasi-hydrostatic hot corona. However, these coronae are more diffuse, and of a lower luminosity, than predicted by the analytic models of White & Frenk because of a substantial increase in entropy at z ~ 1-3. Both the removal of low entropy gas by star formation and energy injection from supernovae contribute to this increase in entropy, but the latter is dominant for halo masses M_200 <~ 10^(12.5) Msun. Only a small fraction of the mass of the hot gas is outflowing as a wind but, because of its high density and metallicity, it contributes disproportionally to the X-ray emission. The bulk of the X-ray emission, however, comes from the diffuse quasi-hydrostatic corona which supplies the fuel for ongoing star formation in discs today. Future deep X-ray observations with high spectral resolution (e.g. with NeXT/ASTRO-H or IXO) should be able to map the velocity structure of the hot gas and test this fundamental prediction of current galaxy formation theory.

We present six simulations of galactic stellar haloes formed by the tidal disruption of accreted dwarf galaxies in a fully cosmological setting. Our model is based on the Aquarius project, a suite of high resolution N-body simulations of individual dark matter haloes. We tag subsets of particles in these simulations with stellar populations predicted by the Galform semi-analytic model. Our method self-consistently tracks the dynamical evolution and disruption of satellites from high redshift. The luminosity function and structural properties of surviving satellites, which agree well with observations, suggest that this technique is appropriate. We find that accreted stellar haloes are assembled between 1<z<7 from less than 5 significant progenitors. These progenitors are old, metal-rich satellites with stellar masses similar to the brightest Milky Way dwarf spheroidals. In contrast to previous stellar halo simulations, we find that several of these major contributors survive as self-bound systems to the present day. Both the number of these significant progenitors and their infall times are inherently stochastic. This results in great diversity among our stellar haloes, which amplifies small differences between the formation histories of their dark halo hosts. The masses and density/surface-brightness profiles of the stellar haloes are consistent with expectations from the Milky Way and M31. Each halo has a complex structure, consisting of well-mixed components, tidal streams, shells and other subcomponents. This structure is not adequately described by smooth models. We find one example of an accreted thick disk. Contrasts in age and metallicity between halo stars and those in surviving satellites are in broad agreement with recent observations. [Abridged]

We investigate the properties of satellite galaxies in cosmological N-body/SPH simulations of galaxy formation in Milky Way-sized haloes. Because of their shallow potential wells, satellite galaxies are very sensitive to heating processes which affect their gas content. Their properties can therefore be used to constrain the nature of feedback processes that regulate galaxy formation. In our simulations, we assume that all the energy produced by supernovae is used as kinetic energy to drive galactic winds. Several of our simulations produce bright, disc-dominated galaxies. We find that wind models in which the wind speed, v_w, is proportional to local velocity dispersion of dark matter, sigma, (and thus the wind mass-loading, eta_w \propto sigma^{-2}) have episodic star formation histories, reproduce the observed satellite luminosity function reasonably well (down to M_v=-7) and match the luminosity-metallicity relation observed in the Local Group satellites. By contrast, models that assume a constant wind speed overproduce faint satellites and predict an incorrect luminosity-metallicity relation. Our simulations therefore suggest that the feedback processes that operate on the scale of satellite galaxies should generate galactic outflows whose intensity varies inversely with the depth of the potential.

We investigate whether the subhalos of Lambda-CDM galaxy halos have potentials consistent with the observed properties of Milky Way satellites, particularly those with high-quality photometric and kinematic data: Fornax, Leo I, Sculptor, Sextans, and Carina. We compare spherical models with isotropic velocity dispersion tensors to the observed, circularly averaged star counts, line-of-sight velocity dispersion profiles and line-of-sight velocity distributions. We identify subhalos within the six high-resolution dark matter halos of the Aquarius Project for which the spherically averaged potentials result in excellent fits to each of the five galaxies. In particular, our simple one-integral models reproduce the observations in the inner regions, proving that these data are fully consistent with Lambda-CDM expectations and do not require cored dark matter distributions. For four of the five satellites the fits require moderately cusped {\it stellar} density profiles. The star count data for Leo I, however, do require a cored distribution of star counts. Current data suggest that these five satellites may be hosted by Lambda-CDM subhalos with maximum circular velocities in the range 10 to 30 km/s.

We present a new method for constructing maps of the secondary temperature fluctuations imprinted on the cosmic microwave background (CMB) radiation by photons propagating through the evolving cosmic gravitational potential. Large cosmological N-body simulations are used to calculate the complete non-linear evolution of the peculiar gravitational potential. Tracing light rays back through the past lightcone of a chosen observer accurately captures the temperature perturbations generated by linear (the integrated Sachs-Wolfe or ISW effect) and non-linear (the Rees-Sciama or RS effect) evolution. These effects give rise to three kinds of non-linear features in the temperature maps. (a) In overdense regions, converging flows of matter induce cold spots of order 100 Mpc in extent which can dominate over the ISW effect at high redshift, and are surrounded by hot rings. (b) In underdense regions, the RS effect enhances ISW cold spots which can be surrounded by weak hot rings. (c) Transverse motions of large lumps of matter produce characteristic dipole features, consisting of adjacent hot and cold spots separated by a few tens of Megaparsecs. These non-linear features are not easily detectable; they modulate the ISW sky maps at about the 10 percent level. The RS effect causes the angular power spectrum to deviate from linear theory at l~50 and generates non-Gaussianity, skewing the one-point distribution function to negative values. Cold spots of similar angular size, but much smaller amplitude than the CMB cold spot reported by Cruz et al. are produced. Joint analysis of our maps and the corresponding galaxy distribution may enable techniques to be developed to detect these non-linear, non-Gaussian features. Our maps are available at http://astro.dur.ac.uk/~cai/ISW

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.

We use a model for the evolution of galaxies in the far-IR based on the LambdaCDM cosmology to make detailed predictions for upcoming cosmological surveys with the Herschel Space Observatory. We use the combined GALFORM semi-analytical galaxy formation model and GRASIL spectrophotometric code to compute galaxy SEDs including the reprocessing of radiation by dust. The model, which is the same as that in Baugh et al. (2005), assumes two different IMFs: a normal solar neighbourhood IMF for quiescent star formation in disks, and a very top-heavy IMF in starbursts triggered by galaxy mergers. We have shown previously that the top-heavy IMF appears necessary to explain the number counts and redshifts of faint sub-mm galaxies. In this paper, we present predictions for galaxy luminosity functions, number counts and redshift distributions in the Herschel imaging bands. We find that source confusion will be a serious problem in the deepest planned surveys. We also show predictions for physical properties such as star formation rates and stellar, gas and halo masses, together with fluxes at other wavelengths (from the far-UV to the radio) relevant for multi-wavelength follow-up observations. We investigate what fraction of the total IR emission from dust and of the high-mass star formation over the history of the Universe should be resolved by planned surveys with Herschel, and find a fraction ~30-50%, depending on confusion. Finally, we show that galaxies in Herschel surveys should be significantly clustered.

We investigate the magnitude and internal alignment of the angular momentum of cold dark matter haloes in simulations with and without baryons. We analyse the cumulative angular momentum profiles of hundreds of thousands of haloes in the Millennium simulation and in a smaller, but higher resolution, simulation, in total spanning 5 orders of magnitude in mass. For haloes of a given mass, the median specific angular momentum increases as j(<r) proportional to r. The direction of the vector varies considerably with radius: the median angle between the inner (< 0.25 Rvir) and total (< Rvir) angular momentum vectors is about 25degr. To investigate how baryons affect halo spin, we use another high resolution simulation, which includes gas cooling, star formation and feedback. This simulation produces a sample of galaxies with a realistic distribution of disc-to-total ratios. The formation of the galaxy spins up the dark matter within 0.1 Rvir such that the specific halo angular momentum increases by approx 50% in the median. The dark matter angular momentum becomes better aligned, but there remains a broad distribution of (mis-)alignments between the halo and the central galaxy, with a median angle between their angular momenta of ~ 30degr. Galaxies have a range of orientations relative to the shape of the halo: half of them have their minor axes misaligned by more than 45degr. Finally, we align a sample of haloes according to the orientation of their galaxies and stack the projected mass distributions. Although the individual haloes are aspherical, galaxy--halo misalignments produce a stacked mass distribution that is indistinguishable from circular. If the misalignment found in our simulations is realistic, it will be extremely difficult for weak lensing studies to measure the shape of dark haloes using this technique.

We use N-body simulations to investigate the radial dependence of the density and velocity dispersion in cold dark matter (CDM) halos. In particular, we explore how closely Q rho/sigma^3, a surrogate measure of the phase-space density, follows a power-law in radius. Our study extends earlier work by considering, in addition to spherically-averaged profiles, local Q-estimates for individual particles, Q_i; profiles based on the ellipsoidal radius dictated by the triaxial structure of the halo, Q_i(r'); and by carefully removing substructures in order to focus on the profile of the smooth halo, Q^s. The resulting Q_i^s(r') profiles follow closely a power law near the center, but show a clear upturn from this trend near the virial radius, r_{200}. The location and magnitude of the deviations are in excellent agreement with the predictions from Bertschinger's spherical secondary-infall similarity solution. In this model, Q \propto r^{-1.875} in the inner, virialized regions, but departures from a power-law occur near r_{200} because of the proximity of this radius to the location of the first shell crossing - the shock radius in the case of a collisional fluid. Particles there have not yet fully virialized, and so Q departs from the inner power-law profile. Our results imply that the power-law nature of $Q$ profiles only applies to the inner regions and cannot be used to predict accurately the structure of CDM halos beyond their characteristic scale radius.

We have implemented a simple model to identify the likely sites of the first stars and galaxies in the high-resolution simulations of the formation of galactic dark matter halos of the Aquarius Project. The first star in a galaxy like the Milky Way formed around redshift z=35; by z=10, the young galaxy contained up to 30000 dark matter haloes capable of forming stars by molecular hydrogen cooling. These minihaloes were strongly clustered and feedback may have severely limited the actual number of Population III stars that formed. By the present day, the remnants of the first stars would be strongly concentrated to the centre of the main halo. If a second generation of long-lived stars formed near the first (the first star relics), we would expect to find half of them within 30kpc/h of the galactic centre and a significant fraction in satellites where they may be singled out by their anomalous metallicity patterns. The first halo in which gas could cool by atomic hydrogen line radiation formed at z=25; by z=10, the number of such first galaxies had increased to ~300. Feedback might have decreased the number of first galaxies at the time they undergo significant star formation, but not the number that survive to the present because near neighbours merge. Half of all the "first galaxies" that form before z=10 merge with the main halo before z~3 and most lose a significant fraction of their mass. However, today there should still be more than 20 remnants orbiting within the central ~30kpc/h of the Milky Way. These satellites have circular velocities of a few kilometers per second or more, comparable to those of known Milky Way dwarfs. They are a promising hunting ground for the remnants of the earliest epoch of star formation.

We describe the spectroscopic target selection for the Galaxy And Mass Assembly (GAMA) survey. The input catalogue is drawn from the Sloan Digital Sky Survey (SDSS) and UKIRT Infrared Deep Sky Survey (UKIDSS). The aim is to measure redshifts for galaxies in three 4x12 degree regions at 9h, 12h and 14.5h, on the celestial equator, with magnitude selections r < 19.4, z < 18.2 and K(AB) < 17.6 over all three regions, and r < 19.8 in the 12-h region. The target density is 1080 per sq. deg. in the 12-h region and 720 per sq. deg. in the other regions. The average GAMA target density and area are compared with completed and ongoing galaxy redshift surveys. The GAMA survey implements a highly complete star-galaxy separation that jointly uses an intensity-profile separator (delta(sg) = r-band psf mag - model mag) as per the SDSS) and a colour separator. The colour separator is defined as delta(sg,jk) = J - K - f(g-i), where f(g-i) is a quadratic fit to the J-K colour of the stellar locus over the range 0.3 < g-i < 2.3. All galaxy populations investigated are well separated with delta(sg,jk) > 0.2. From two years out of a three-year AAOmega program on the Anglo-Australian Telescope, we have obtained 79599 unique galaxy redshifts. Previously known redshifts in the GAMA region bring the total up to 98497. The median galaxy redshift is 0.2 with 99% at z < 0.5. We present some of the global statistical properties of the survey, including K-band galaxy counts, colour-redshift relations and preliminary n(z).

Recent analysis of the Milky Way's satellite galaxies reveals that these objects share a common central mass density, even though their luminosities range over five orders of magnitude. This observation can be understood in the context of galaxy formation theory by quantifying the factors which restrict the central mass density to a small range. One limit is set by the maximum mass that can collapse into a given region by the hierarchical growth of structure in the standard cold dark matter cosmology. Another limit comes from the natural thresholds which exist for gas to be able to cool and form a galaxy. The wide range of luminosities in these satellites reflect the effects of supernova feedback on the fraction of cooled baryons which are retained.

The Galaxy And Mass Assembly (GAMA) project is the latest in a tradition of large galaxy redshift surveys, and is now underway on the 3.9m Anglo-Australian Telescope at Siding Spring Observatory. GAMA is designed to map extragalactic structures on scales of 1kpc - 1Mpc in complete detail to a redshift of z~0.2, and to trace the distribution of luminous galaxies out to z~0.5. The principal science aim is to test the standard hierarchical structure formation paradigm of Cold Dark Matter (CDM) on scales of galaxy groups, pairs, discs, bulges and bars. We will measure (1) the Dark Matter Halo Mass Function (as inferred from galaxy group velocity dispersions); (2) baryonic processes, such as star formation and galaxy formation efficiency (as derived from Galaxy Stellar Mass Functions); and (3) the evolution of galaxy merger rates (via galaxy close pairs and galaxy asymmetries). Additionally, GAMA will form the central part of a new galaxy database, which aims to contain 275,000 galaxies with multi-wavelength coverage from coordinated observations with the latest international ground- and space-based facilities: GALEX, VST, VISTA, WISE, HERSCHEL, GMRT and ASKAP. Together, these data will provide increased depth (over 2 magnitudes), doubled spatial resolution (0.7"), and significantly extended wavelength coverage (UV through Far-IR to radio) over the main SDSS spectroscopic survey for five regions, each of around 50 deg^2. This database will permit detailed investigations of the structural, chemical, and dynamical properties of all galaxy types, across all environments, and over a 5 billion year timeline.

A heuristic greedy algorithm is developed for efficiently tiling spatially dense redshift surveys. In its first application to the Galaxy and Mass Assembly (GAMA) redshift survey we find it rapidly improves the spatial uniformity of our data, and naturally corrects for any spatial bias introduced by the 2dF multi object spectrograph. We make conservative predictions for the final state of the GAMA redshift survey after our final allocation of time, and can be confident that even if worse than typical weather affects our observations, all of our main survey requirements will be met.

The Herschel ATLAS is the largest open-time key project that will be carried out on the Herschel Space Observatory. It will survey 510 square degrees of the extragalactic sky, four times larger than all the other Herschel surveys combined, in five far-infrared and submillimetre bands. We describe the survey, the complementary multi-wavelength datasets that will be combined with the Herschel data, and the six major science programmes we are undertaking. Using new models based on a previous submillimetre survey of galaxies, we present predictions of the properties of the ATLAS sources in other wavebands.

We study the formation histories and present-day structure of satellite galaxies formed in a high resolution hydrodynamic simulation of a Milky Way-like galaxy. The simulated satellites span nearly 4 orders of magnitude in luminosity but have a very similar mass within their inner 600 pc, ~ 3 10^7 M_solar, with very little scatter. This result is in agreement with the recent measurements for dwarf spheroidal galaxies (dSphs) in the Milky Way by Strigari et al. In our simulations a preferred mass scale arises naturally from the effects of the early reionisation of gas. These impose a sharp threshold of ~ 12 km/s on the circular velocity of haloes which can cool gas and make stars. At the present day, subhaloes that host satellites as luminous as the classical Milky Way dwarfs (L_V > 2.6 10^5 L_solar), have typically grown to have circular velocities of $> 20 km/s. There are, however, subhaloes with similar circular velocities today which were, nevertheless, below threshold at reionisation and thus remain dark. Star formation in above-threshold haloes is truncated when the halo is accreted into the main galaxy progenitor. Thus, many properties of today's dwarf satellites such as their luminosity and star formation rate are set by their accretion time.

Virtually all massive galaxies, including our own, host central black holes ranging in mass from millions to billions of solar masses. The growth of these black holes releases vast amounts of energy that powers quasars and other weaker active galactic nuclei. A tiny fraction of this energy, if absorbed by the host galaxy, could halt star formation by heating and ejecting ambient gas. A central question in galaxy evolution is the degree to which this process has caused the decline of star formation in large elliptical galaxies, which typically have little cold gas and few young stars, unlike spiral galaxies.

[Abridged] We present the first results of hydrodynamical simulations that follow the formation of galaxies to z=0 in spherical regions of radius ~20 Mpc/h drawn from the Millennium Simulation. The regions have overdensities that deviate by (-2, -1, 0, +1, +2)sigma from the cosmic mean, where sigma is the rms mass fluctuation on a scale of ~20Mpc/h at z=1.5. The simulations have mass resolution of up to 10^6 Msun/h, cover the entire range of large-scale environments and allow extrapolation of statistics to the entire 500 (Mpc/h)^3 Millennium volume. They include gas cooling, photoheating from an ionising background, SNe feedback and winds, but no AGN. We find that the specific SFR density at z <~ 10 varies systematically from region to region by up to an order of magnitude, but the global value, averaged over all volumes, reproduces observational data. Massive, compact galaxies, similar to those observed in the GOODS fields, form in the overdense regions as early as z=6, but do not appear in the underdense regions until z~3. These environmental variations are not caused by a dependence of the star formation properties on environment, but rather by a strong variation of the halo mass function from one environment to another, with more massive haloes forming preferentially in the denser regions. At all epochs, stars form most efficiently in haloes of circular velocity ~ 250 km/s. However, the star formation history exhibits a form of "downsizing" (even in the absence of AGN): the stars comprising massive galaxies at z=0 have mostly formed by z=1-2, whilst those comprising smaller galaxies typically form at later times. However, additional feedback is required to limit star formation in massive galaxies at late times.

We introduce a new technique based on artificial neural networks which allows us to make accurate predictions for the spectral energy distributions (SEDs) of large samples of galaxies, at wavelengths ranging from the far-ultra-violet to the sub-millimetre and radio. The neural net is trained to reproduce the SEDs predicted by a hybrid code comprised of the GALFORM semi-analytical model of galaxy formation, which predicts the full star formation and galaxy merger histories, and the GRASIL spectro-photometric code, which carries out a self-consistent calculation of the SED, including absorption and emission of radiation by dust. Using a small number of galaxy properties predicted by GALFORM, the method reproduces the luminosities of galaxies in the majority of cases to within 10% of those computed directly using GRASIL. The method performs best in the sub-mm and reasonably well in the mid-infrared and the far-ultra-violet. The luminosity error introduced by the method has negligible impact on predicted statistical distributions, such as luminosity functions or colour distributions of galaxies. We use the neural net to predict the overlap between galaxies selected in the rest-frame UV and in the observer-frame sub-mm at z=2. We find that around half of the galaxies with a 850um flux above 5 mJy should have optical magnitudes brighter than R_AB < 25 mag. However, only 1% of the galaxies selected in the rest-frame UV down to R_AB < 25 mag should have 850um fluxes brighter than 5 mJy. Our technique will allow the generation of wide-angle mock catalogues of galaxies selected at rest-frame UV or mid- and far-infrared wavelengths.

We study the mass, velocity dispersion, and anisotropy profiles of $\Lambda$CDM halos using a suite of N-body simulations of unprecedented numerical resolution (the {\it Aquarius Project}). Our analysis confirms a number of results claimed by earlier work, and clarifies a few issues where conflicting claims may be found in the recent literature. The spherically-averaged density profile becomes progressively shallower inwards and, at the innermost resolved radius, the logarithmic slope is $\gamma \equiv -$d$\ln\rho/$d$\ln r \simlt 1$. Asymptotic inner slopes as steep as the recently claimed $\rho \propto r^{-1.2}$ are clearly ruled out. The radial dependence of $\gamma$ is well approximated by a power-law, $\gamma \propto r^{\alpha}$ (the Einasto profile). The shape parameter, $\alpha$, varies slightly but significantly from halo to halo, implying that the mass profiles of $\Lambda$CDM halos are not strictly universal: different halos cannot, in general, be rescaled to look identical. Departures from similarity are also seen in velocity dispersion profiles and correlate with those in density profiles so as to preserve a power-law form for the spherically averaged pseudo-phase-space density, $\rho/\sigma^3\propto r^{-1.875}$. Our conclusions are reliable down to radii below 0.4% of the virial radius, providing well-defined predictions for halo structure when baryonic effects are neglected, and thus an instructive theoretical template against which the modifications induced by the baryonic components of real galaxies can be judged.

We use high-resolution Aquarius simulations of Milky Way-sized haloes in the LCDM cosmology to study the effects of dark matter substructures on gravitational lensing. Each halo is resolved with ~ 10^8 particles (at a mass resolution ~ 10^3-4 M_sun/h) within its virial radius. Subhaloes with masses larger than 10^5 M_sun/h are well resolved, an improvement of at least two orders of magnitude over previous lensing studies. We incorporate a baryonic component modelled as a Hernquist profile and account for the response of the dark matter via adiabatic contraction. We focus on the "anomalous" flux ratio problem, in particular on the violation of the cusp-caustic relation due to substructures. We find that subhaloes with masses less than ~ 10^8 M_sun/h play an important role in causing flux anomalies; such low mass subhaloes have been unresolved in previous studies. There is large scatter in the predicted flux ratios between different haloes and between different projections of the same halo. In some cases, the frequency of predicted anomalous flux ratios is comparable to that observed for the radio lenses, although in most cases it is not. The probability for the simulations to reproduce the observed violations of the cusp lenses is about 0.001. We therefore conclude that the amount of substructure in the central regions of the Aquarius haloes is insufficient to explain the observed frequency of violations of the cusp-caustic relation. These conclusions are based purely on our dark matter simulations which ignore the effect of baryons on subhalo survivability.

We propose an X-ray mission called Xenia to search for decaying superweakly interacting Dark Matter particles (super-WIMP) with a mass in the keV range. The mission and its observation plan are capable of providing a major break through in our understanding of the nature of Dark Matter (DM). It will confirm, or reject, predictions of a number of particle physics models by increasing the sensitivity of the search for decaying DM by about two orders of magnitude through a wide-field imaging X-ray spectrometer in combination with a dedicated observation program. The proposed mission will provide unique limits on the mixing angle and mass of neutral leptons, right handed partners of neutrinos, which are important Dark Matter candidates. The existence of these particles is strongly motivated by observed neutrino flavor oscillations and the problem of baryon asymmetry of the Universe. In super-WIMP models, the details of the formation of the cosmic web are different from those of LambdaCDM. The proposed mission will, in addition to the search for decaying Dark Matter, provide crucial insight into the nature of DM by studying the structure of the "cosmic web". This will be done by searching for missing baryons in emission, and by using gamma-ray bursts as backlight to observe the warm-hot intergalactic media in absorption.

The huge size and uniformity of the Sloan Digital Sky Survey makes possible an exacting test of current models of galaxy formation. We compare the predictions of the GALFORM semi-analytical galaxy formation model for the luminosities, morphologies, colours and scale-lengths of local galaxies. GALFORM models the luminosity and size of the disk and bulge components of a galaxy, and so we can compute quantities which can be compared directly with SDSS observations, such as the Petrosian magnitude and the Sersic index. We test the predictions of two published models set in the cold dark matter cosmology: the Baugh et al. (2005) model, which assumes a top-heavy initial mass function (IMF) in starbursts and superwind feedback, and the Bower et al. (2006) model, which uses AGN feedback and a standard IMF. The Bower et al model better reproduces the overall shape of the luminosity function, the morphology-luminosity relation and the colour bimodality observed in the SDSS data, but gives a poor match to the size-luminosity relation. The \Baugh et al. model successfully predicts the size-luminosity relation for late-type galaxies. Both models fail to reproduce the sizes of bright early-type galaxies. These problems highlight the need to understand better both the role of feedback processes in determining galaxy sizes, in particular the treatment of the angular momentum of gas reheated by supernovae, and the sizes of the stellar spheroids formed by galaxy mergers and disk instabilities.