We examine the influence of dark matter halo assembly on the evolution of a simulated $\sim L^\star$ galaxy. Starting from a zoom-in simulation of a star-forming galaxy evolved with the EAGLE galaxy formation model, we use the genetic modification technique to create a pair of complementary assembly histories: one in which the halo assembles later than in the unmodified case, and one in which it assembles earlier. Delayed assembly leads to the galaxy exhibiting a greater present-day star formation rate than its unmodified counterpart, whilst in the accelerated case the galaxy quenches at $z\simeq 1$, and becomes spheroidal. We simulate each assembly history nine times, adopting different seeds for the random number generator used by EAGLE's stochastic subgrid implementations of star formation and feedback. The systematic changes driven by differences in assembly history are significantly stronger than the random scatter induced by this stochasticity. The sensitivity of $\sim L^\star$ galaxy evolution to dark matter halo assembly follows from the close coupling of the growth histories of the central black hole (BH) and the halo, such that earlier assembly fosters the formation of a more massive BH, and more efficient expulsion of circumgalactic gas. In response to this expulsion, the circumgalactic medium reconfigures at a lower density, extending its cooling time and thus inhibiting the replenishment of the interstellar medium. Our results indicate that halo assembly history significantly influences the evolution of $\sim L^\star$ central galaxies, and that the expulsion of circumgalactic gas is a crucial step in quenching them.

We present genetIC, a new code for generating initial conditions for cosmological N-body simulations. The code allows precise, user-specified alterations to be made to arbitrary regions of the simulation (while maintaining consistency with the statistical ensemble). These "genetic modifications" allow, for example, the history, mass, or environment of a target halo to be altered in order to study the effect on their evolution. The code natively supports initial conditions with nested zoom regions at progressively increasing resolution. Modifications in the high-resolution region must propagate self-consistently onto the lower resolution grids; to enable this while maintaining a small memory footprint, we introduce a Fourier-space filtering approach to generating fields at variable resolution. Due to a close correspondence with modifications, constrained initial conditions can also be produced by genetIC (for example with the aim of matching structures in the local Universe). We test the accuracy of modifications performed within zoom initial conditions. The code achieves sub-percent precision, which is easily sufficient for current applications in galaxy formation.

We study how star formation is regulated in low-mass field dwarf galaxies ($10^5 \leq M_{\star} \leq 10^6 \, \text{M}_{\odot}$), using cosmological high-resolution ($3 \, \text{pc}$) hydrodynamical simulations. Cosmic reionization quenches star formation in all our simulated dwarfs, but three galaxies with final dynamical masses of $3 \times 10^{9} \,\text{M}_{\odot}$ are subsequently able to replenish their interstellar medium by slowly accreting gas. Two of these galaxies re-ignite and sustain star formation until the present day at an average rate of $10^{-5} \, \text{M}_{\odot} \, \text{yr}^{-1}$, highly reminiscent of observed low-mass star-forming dwarf irregulars such as Leo T. The resumption of star formation is delayed by several billion years due to residual feedback from stellar winds and Type Ia supernovae; even at $z=0$, the third galaxy remains in a temporary equilibrium with a large gas content but without any ongoing star formation. Using the "genetic modification'' approach, we create an alternative mass growth history for this gas-rich quiescent dwarf and show how a small $(0.2\,\mathrm{dex})$ increase in dynamical mass can overcome residual stellar feedback, re-igniting star formation. The interaction between feedback and mass build-up produces a diversity in the stellar ages and gas content of low-mass dwarfs, which will be probed by combining next-generation HI and imaging surveys.

Using cosmological hydrodynamic simulations with physically motivated models of super-massive black hole (SMBH) formation and growth, we compare the assembly of Milky Way-mass ( $M_{\mathrm{vir}} \approx 7 \times 10^{11}$ $M_{\odot}$ at $z = 0$) galaxies in cold dark matter (CDM) and self-interacting dark matter (SIDM) models. Our SIDM model adopts a constant cross-section of 1 cm$^2$/g. We find that SMBH formation is suppressed in the early universe due to SIDM coring. SMBH-SMBH mergers are also suppressed in SIDM as a consequence of the lower number of SMBHs formed. Lack of initial merger-driven SMBH growth in turn delays SMBH growth by billions of years in SIDM compared to CDM. Further, we find that this delayed growth suppresses SMBH accretion in the largest progenitors of the main SIDM galaxies during the first 5 Gyrs of their evolution. Nonetheless, by $z = 0.8$ the CDM and SIDM SMBH masses differ only by around 0.2 dex, so that both remain compatible with the $M_{BH}-M_{*}$ relation. We show that the reduced accretion causes the SIDM SMBHs to less aggressively regulate star formation in their host galaxies than their CDM counterparts, resulting in a factor of 3 or more stars being produced over the lifetime of the SIDM galaxies compared to the CDM galaxies. Our results highlight a new way in which SIDM can affect the growth and merger history of SMBHs and ultimately give rise to very different galaxy evolution compared to the classic CDM model.

Using a suite of hydrodynamical simulations with cold dark matter, baryons, and neutrinos, we present a detailed study of the effect of massive neutrinos on the 1-D and 3-D flux power spectra of the Lyman-$\alpha$ (Ly$\alpha$) forest. The presence of massive neutrinos in cosmology induces a scale- and time-dependent suppression of structure formation that is strongest on small scales. Measuring this suppression is a key method for inferring neutrino masses from cosmological data, and is one of the main goals of ongoing and future surveys like eBOSS, DES, LSST, Euclid or DESI. The clustering in the Ly$\alpha$ forest traces the quasi-linear power at late times and on small scales. In combination with observations of the cosmic microwave background, the forest therefore provides some of the tightest constraints on the sum of the neutrino masses. However there is a well-known degeneracy between $\Sigma m_{\nu}$ and the amplitude of perturbations in the linear matter power spectrum. We study the corresponding degeneracy in the 1-D flux power spectrum of the Ly$\alpha$ forest, and for the first time also study this degeneracy in the 3-D flux power spectrum. We show that the non-linear effects of massive neutrinos on the Ly$\alpha$ forest, beyond the effect of linear power amplitude suppression, are negligible, and this degeneracy persists in the Ly$\alpha$ forest observables to a high precision. We discuss the implications of this degeneracy for choosing parametrisations of the Ly$\alpha$ forest for cosmological analysis.

We determine the cosmic abundance of molecular hydrogen (H2) in the local universe from the xCOLD GASS survey. To constrain the H2 mass function at low masses and correct for the effect of the lower stellar mass limit of 10^9 Msun in the xCOLD GASS survey, we use an empirical approach based on an observed scaling relation between star formation rate and gas mass. We also constrain the HI and HI+H2 mass functions using the xGASS survey, and compare it to the HI mass function from the ALFALFA survey. We find the cosmic abundance of molecular gas in the local Universe to be Omega_H2=(5.34+/-0.47)x10^-5 h^-1. Molecular gas accounts for 19.6 +/- 3.9% of the total abundance of cold gas, Omega_HI+H2=(4.66+/-0.70)x10^-4 h^-1. Galaxies with stellar masses in excess of 10^9 Msun account for 89% of the molecular gas in the local Universe, while in comparison such galaxies only contain 73% of the cold atomic gas as traced by the HI 21cm line. The xCOLD GASS CO, molecular gas and cold gas mass functions and Omega_H2 measurements provide constraints for models of galaxy evolution and help to anchor blind ALMA and NOEMA surveys attempting to determine the abundance of molecular gas at high redshifts.

We introduce the "Engineering Dwarfs at Galaxy formation's Edge" (EDGE) project to study the cosmological formation and evolution of the smallest galaxies in the Universe. In this first paper, we explore the effects of resolution and sub-grid physics on a single low mass halo ($M_{\rm halo}=10^{9}~M_\odot$), simulated to redshift $z=0$ at a mass and spatial resolution of $\sim 20~M_\odot$ and $\sim 3$ pc. We consider different star formation prescriptions, supernova feedback strengths and on-the-fly radiative transfer (RT). We show that RT changes the mode of galactic self-regulation at this halo mass, suppressing star formation by causing the interstellar and circumgalactic gas to remain predominantly warm ($\sim 10^4$ K) even before cosmic reionisation. By contrast, without RT, star formation regulation occurs only through starbursts and their associated vigorous galactic outflows. In spite of this difference, the entire simulation suite (with the exception of models without any feedback) matches observed dwarf galaxy sizes, velocity dispersions, $V$-band magnitudes and dynamical mass-to-light-ratios. This is because such structural scaling relations are predominantly set by the host dark matter halo, with the remaining model-to-model variation being smaller than the observational scatter. We find that only the stellar mass-metallicity relation differentiates the galaxy formation models. Explosive feedback ejects more metals from the dwarf, leading to a lower metallicity at a fixed stellar mass. We conclude that the stellar mass-metallicity relation of the very smallest galaxies provides a unique constraint on galaxy formation physics.

The statistical power of Lyman-${\alpha}$ forest Baryon Acoustic Oscillation (BAO) measurements is set to increase significantly in the coming years as new instruments such as the Dark Energy Spectroscopic Instrument deliver progressively more constraining data. Generating mock datasets for such measurements will be important for validating analysis pipelines and evaluating the effects of systematics. With such studies in mind, we present LyaCoLoRe: a package for producing synthetic Lyman-${\alpha}$ forest survey datasets for BAO analyses. LyaCoLoRe transforms initial Gaussian random field skewers into skewers of transmitted flux fraction via a number of fast approximations. In this work we explain the methods of producing mock datasets used in LyaCoLoRe, and then measure correlation functions on a suite of realisations of such data. We demonstrate that we are able to recover the correct BAO signal, as well as large-scale bias parameters similar to literature values. Finally, we briefly describe methods to add further astrophysical effects to our skewers - high column density systems and metal absorbers - which act as potential complications for BAO analyses.

We demonstrate how the least luminous galaxies in the Universe, ultra-faint dwarf galaxies, are sensitive to their dynamical mass at the time of cosmic reionization. We select a low-mass ($\sim \text{1.5} \times 10^{9} \, \text{M}_{\odot}$) dark matter halo from a cosmological volume, and perform zoom hydrodynamical simulations with multiple alternative histories using "genetically modified" initial conditions. Earlier forming ultra-faints have higher stellar mass today, due to a longer period of star formation before their quenching by reionization. Our histories all converge to the same final dynamical mass, demonstrating the existence of extended scatter ($\geq$ 1 dex) in stellar masses at fixed halo mass due to the diversity of possible histories. One of our variants builds less than 2 % of its final dynamical mass before reionization, rapidly quenching in-situ star formation. The bulk of its final stellar mass is later grown by dry mergers, depositing stars in the galaxy's outskirts and hence expanding its effective radius. This mechanism constitutes a new formation scenario for highly diffuse ($\text{r}_{1 /2} \sim 820 \, \text{pc}$, $\sim 32 \, \text{mag arcsec}^2$), metal-poor ($\big[ \mathrm{Fe}\, / \mathrm{H} \big]= -2.9$), ultra-faint ($\mathcal{M}_V= -5.7$) dwarf galaxies within the reach of next-generation low surface brightness surveys.

We investigate the nonlinear dynamics of cold atom systems that can in principle serve as quantum simulators of false vacuum decay. The analog false vacuum manifests as a metastable vacuum state for the relative phase in a two-species Bose-Einstein condensate (BEC), induced by a driven periodic coupling between the two species. In the appropriate low energy limit, the evolution of the relative phase is approximately governed by a relativistic wave equation exhibiting true and false vacuum configurations. In previous work, a linear stability analysis identified exponentially growing short-wavelength modes driven by the time-dependent coupling. These modes threaten to destabilize the analog false vacuum. Here, we employ numerical simulations of the coupled Gross-Pitaevski equations (GPEs) to determine the non-linear evolution of these linearly unstable modes. We find that unless a physical mechanism modifies the GPE on short length scales, the analog false vacuum is indeed destabilized. We briefly discuss various physically expected corrections to the GPEs that may act to remove the exponentially unstable modes. To investigate the resulting dynamics in cases where such a removal mechanism exists, we implement a hard UV cutoff that excludes the unstable modes as a simple model for these corrections. We use this to study the range of phenomena arising from such a system. In particular, we show that by modulating the strength of the time-dependent coupling, it is possible to observe the crossover between a second and first order phase transition out of the false vacuum.

We present a generalization of our recently proposed machine learning framework, aiming to provide new physical insights into dark matter halo formation. We investigate the impact of the initial density and tidal shear fields on the formation of haloes over the mass range $11.4 \leq \log(M/M_{\odot}) \leq 13.4$. The algorithm is trained on an N-body simulation to infer the final mass of the halo to which each dark matter particle will later belong. We then quantify the difference in the predictive accuracy between machine learning models using a metric based on the Kullback-Leibler divergence. We first train the algorithm with information about the density contrast in the particles' local environment. The addition of tidal shear information does not yield an improved halo collapse model over one based on density information alone; the difference in their predictive performance is consistent with the statistical uncertainty of the density-only based model. This implies that our machine learning setup does not identify any significant role for the tidal shear in determining halo masses. This result is confirmed as we verify the ability of the initial conditions-to-halo mass mapping learnt from one simulation to generalize to independent simulations. Our work illustrates the broader potential of developing interpretable machine learning frameworks to gain physical understanding of non-linear large-scale structure formation.

The modeling of cosmological observables is based on the statistics of the matter density, velocity and gravitational fields in the Universe as a function of time. Typically, calculations are restricted to equal time correlations, where any given fields are evaluated at the same redshift. For some applications, it is necessary to make accurate predictions of unequal time correlators, where the fields considered are evaluated at different redshifts. In this work, we show that the Zeldovich approximation provides an accurate (<10%) analytical prescription to model unequal time correlators, which we validate against numerical N-body simulations. The Zeldovich approximation introduces a scale-dependent exponential suppression of unequal time correlators, which depends on cosmology and the redshifts of the fields considered. Comparing the Zeldovich case to previous approximations, we show that it can yield accurate predictions for wavenumbers that extend well into the nonlinear regime. However, we also show that correlations over such scales are typically suppressed by the geometry of the lightcone, and thus should normally be negligible for cosmology with galaxy surveys. We discuss potential exceptions, such as intrinsic galaxy alignments, where unequal time correlators could play a role in the modeling of the observables.

We introduce a new picture of vacuum decay which, in contrast to existing semiclassical techniques, provides a real-time description and does not rely on classically-forbidden tunneling paths. Using lattice simulations, we observe vacuum decay via bubble formation by generating realizations of vacuum fluctuations and evolving with the classical equations of motion. The decay rate obtained from an ensemble of simulations is in excellent agreement with existing techniques. Future applications include bubble correlation functions, fast decay rates, and decay of non-vacuum states.

We apply our recently proposed "quadratic genetic modification" approach to generating and testing the effects of alternative mass accretion histories for a single $\Lambda$CDM halo. The goal of the technique is to construct different formation histories, varying the overall contribution of mergers to the fixed final mass. This enables targeted studies of galaxy and dark matter halo formation's sensitivity to the smoothness of mass accretion. Here, we focus on two dark matter haloes, each with four different mass accretion histories. We find that the concentration of both haloes systematically decreases as their merger history becomes smoother. This causal trend tracks the known correlation between formation time and concentration parameters in the overall halo population. At fixed formation time, we further establish that halo concentrations are sensitive to the order in which mergers happen. This ability to study an individual halo's response to variations in its history is highly complementary to traditional methods based on emergent correlations from an extended halo population.

We present methods for interpolating between the 1-D flux power spectrum of the Lyman-$\alpha$ forest, as output by cosmological hydrodynamic simulations. Interpolation is necessary for cosmological parameter estimation due to the limited number of simulations possible. We construct an emulator for the Lyman-$\alpha$ forest flux power spectrum from $21$ small simulations using Latin hypercube sampling and Gaussian process interpolation. We show that this emulator has a typical accuracy of 1.5% and a worst-case accuracy of 4%, which compares well to the current statistical error of 3 - 5% at $z < 3$ from BOSS DR9. We compare to the previous state of the art, quadratic polynomial interpolation. The Latin hypercube samples the entire volume of parameter space, while quadratic polynomial emulation samples only lower-dimensional subspaces. The Gaussian process provides an estimate of the emulation error and we show using test simulations that this estimate is reasonable. We construct a likelihood function and use it to show that the posterior constraints generated using the emulator are unbiased. We show that our Gaussian process emulator has lower emulation error than quadratic polynomial interpolation and thus produces tighter posterior confidence intervals, which will be essential for future Lyman-$\alpha$ surveys such as DESI.

The Lyman-alpha forest provides strong constraints on both cosmological parameters and intergalactic medium astrophysics, which are forecast to improve further with the next generation of surveys including eBOSS and DESI. As is generic in cosmological inference, extracting this information requires a likelihood to be computed throughout a high-dimensional parameter space. Evaluating the likelihood requires a robust and accurate mapping between the parameters and observables, in this case the 1D flux power spectrum. Cosmological simulations enable such a mapping, but due to computational time constraints can only be evaluated at a handful of sample points; "emulators" are designed to interpolate between these. The problem then reduces to placing the sample points such that an accurate mapping is obtained while minimising the number of expensive simulations required. To address this, we introduce an emulation procedure that employs Bayesian optimisation of the training set for a Gaussian process interpolation scheme. Starting with a Latin hypercube sampling (other schemes with good space-filling properties can be used), we iteratively augment the training set with extra simulations at new parameter positions which balance the need to reduce interpolation error while focussing on regions of high likelihood. We show that smaller emulator error from the Bayesian optimisation propagates to smaller widths on the posterior distribution. Even with fewer simulations than a Latin hypercube, Bayesian optimisation shrinks the 95% credible volume by 90% and, e.g., the 1 sigma error on the amplitude of small-scale primordial fluctuations by 38%. This is the first demonstration of Bayesian optimisation applied to large-scale structure emulation, and we anticipate the technique will generalise to many other probes such as galaxy clustering, weak lensing and 21cm.

We present the first results from RomulusC, the highest resolution cosmological hydrodynamic simulation of a galaxy cluster run to date. RomulusC, a zoom-in simulation of a halo with $z=0$ mass $10^{14}$ M$_{\odot}$, is run with the same sub-grid physics and resolution as Romulus25 (Tremmel et al. 2017). With unprecedented mass and spatial resolution, RomulusC represents a unique opportunity to study the evolution of galaxies in dense environments down to dwarf masses. We demonstrate that RomulusC results in an intracluster medium (ICM) consistent with observations. The star formation history and stellar mass of the brightest cluster galaxy (BCG) is consistent with observations and abundance matching results, indicating that our sub-grid models, optimized only to reproduce observations of field dwarf and Milky Way mass galaxies, are able to produce reasonable galaxy masses and star formation histories in much higher mass systems. Feedback from supermassive black holes (SMBHs) regulates star formation by driving large-scale, collimated outflows that coexist with a low entropy core. We find that non-BCG cluster member galaxies are substantially quenched compared to the field down to dwarf galaxy masses and, at low masses, quenching is seen to have no dependence on mass or distance from the cluster center. This enhanced quenched population extends beyond $R_{200}$ and is in place at high redshift. Similarly, we predict that SMBH activity is significantly suppressed within clusters outside of the BCG, but show how the effect could be lost when only focusing on the brightest AGN in the most massive galaxies.

We test a method to reduce unwanted sample variance when predicting Lyman-$\alpha$ (ly$\alpha$) forest power spectra from cosmological hydrodynamical simulations. Sample variance arises due to sparse sampling of modes on large scales and propagates to small scales through non-linear gravitational evolution. To tackle this, we generate initial conditions in which the density perturbation amplitudes are {\it fixed} to the ensemble average power spectrum -- and are generated in {\it pairs} with exactly opposite phases. We run $50$ such simulations ($25$ pairs) and compare their performance against $50$ standard simulations by measuring the ly$\alpha$ 1D and 3D power spectra at redshifts $z=2$, 3, and 4. Both ensembles use periodic boxes of $40$ Mpc/h containing $512^3$ particles each of dark matter and gas. As a typical example of improvement, for wavenumbers $k=0.25$ h/Mpc at $z=3$, we find estimates of the 1D and 3D power spectra converge $34$ and $12$ times faster in a paired-fixed ensemble compared with a standard ensemble. We conclude that, by reducing the computational time required to achieve fixed accuracy on predicted power spectra, the method frees up resources for exploration of varying thermal and cosmological parameters -- ultimately allowing the improved precision and accuracy of statistical inference.

We examine the effects of SMBH feedback on the CGM using a cosmological hydrodynamic simulation \citep[{\sc Romulus25};][]{Tremmel2017} and a set of four zoom-in `genetically modified' Milky Way-mass galaxies sampling different evolutionary paths. By tracing the distribution of metals in the circumgalactic medium (CGM), we show that \ion{O}{6} is a sensitive indicator of supermassive black hole (SMBH) feedback. First, we calculate the column densities of \ion{O}{6} in simulated Milky Way-mass galaxies and compare them with observations from the COS-Halos Survey. Our simulations show column densities of \ion{O}{6} in the CGM consistent with those of COS-Halos star forming and quenched galaxies. These results contrast with those from previous simulation studies which typically underproduce CGM column densities of \ion{O}{6}. We determine that a galaxy's star formation history and assembly record have little effect on the amount of \ion{O}{6} in its CGM. Instead, column densities of \ion{O}{6} are closely tied to galaxy halo mass and BH growth history. The set of zoom-in, genetically modified Milky Way-mass galaxies indicates that the SMBH drives highly metal-enriched material out into its host galaxy's halo which in turn elevates the column densities of \ion{O}{6} in the CGM.

We train a machine learning algorithm to learn cosmological structure formation from N-body simulations. The algorithm infers the relationship between the initial conditions and the final dark matter haloes, without the need to introduce approximate halo collapse models. We gain insights into the physics driving halo formation by evaluating the predictive performance of the algorithm when provided with different types of information about the local environment around dark matter particles. The algorithm learns to predict whether or not dark matter particles will end up in haloes of a given mass range, based on spherical overdensities. We show that the resulting predictions match those of spherical collapse approximations such as extended Press-Schechter theory. Additional information on the shape of the local gravitational potential is not able to improve halo collapse predictions; the linear density field contains sufficient information for the algorithm to also reproduce ellipsoidal collapse predictions based on the Sheth-Tormen model. We investigate the algorithm's performance in terms of halo mass and radial position and perform blind analyses on independent initial conditions realisations to demonstrate the generality of our results.

We present Tangos, a Python framework and web interface for database-driven analysis of numerical structure formation simulations. To understand the role that such a tool can play, consider constructing a history for the absolute magnitude of each galaxy within a simulation. The magnitudes must first be calculated for all halos at all timesteps and then linked using a merger tree; folding the required information into a final analysis can entail significant effort. Tangos is a generic solution to this information organization problem, aiming to free users from the details of data management. At the querying stage, our example of gathering properties over history is reduced to a few clicks or a simple, single-line Python command. The framework is highly extensible; in particular, users are expected to define their own properties which tangos will write into the database. A variety of parallelization options are available and the raw simulation data can be read using existing libraries such as pynbody or yt. Finally, tangos-based databases and analysis pipelines can easily be shared with collaborators or the broader community to ensure reproducibility. User documentation is provided separately.

The initial conditions of cosmological simulations are commonly drawn from a Gaussian ensemble. The limited number of modes inside a simulation volume gives rise to statistical fluctuations known as \textit{sample variance}, limiting the accuracy of simulation predictions. Fixed fields offer an alternative initialization strategy; they have the same power spectrum as standard Gaussian fields but without intrinsic amplitude scatter at linear order. Paired fixed fields consists of two fixed fields with opposite phases that cancel phase correlations which otherwise induce second-order scatter in the non-linear power spectrum. We study the statistical properties of those fields for 19 different quantities at different redshifts through a large set of 600 N-body and 506 state-of-the-art magneto-hydrodynamic simulations covering a wide range of scales, mass and spatial resolutions. We find that paired fixed simulations do not introduce a bias on any of the examined quantities. We quantify the statistical improvement brought by these simulations, over standard ones, on different power spectra such as matter, halos, CDM, gas, stars, black-holes and magnetic fields, finding that they can reduce their variance by factors as large as $10^6$. We quantify the improvement achieved by fixing and by pairing, showing that sample variance in some quantities can be highly suppressed by pairing after fixing. Paired fixed simulations do not change the scatter in quantities such as the probability distribution function of matter density, or the halo, void or stellar mass functions. We argue that procedures aiming at reducing the sample variance of those quantities are unlikely to work. Our results show that paired fixed simulations do not affect either mean relations or scatter of galaxy properties, and suggest that the information embedded in 1-pt statistics is highly complementary to that in clustering.

We examine the scalings of galactic outflows with halo mass across a suite of twenty high-resolution cosmological zoom galaxy simulations covering halo masses from 10^9.5 - 10^12 M_sun. These simulations self-consistently generate outflows from the available supernova energy in a manner that successfully reproduces key galaxy observables including the stellar mass-halo mass, Tully-Fisher, and mass-metallicity relations. We quantify the importance of ejective feedback to setting the stellar mass relative to the efficiency of gas accretion and star formation. Ejective feedback is increasingly important as galaxy mass decreases; we find an effective mass loading factor that scales as v_circ^(-2.2), with an amplitude and shape that is invariant with redshift. These scalings are consistent with analytic models for energy-driven wind, based solely on the halo potential. Recycling is common: about half the outflow mass across all galaxy masses is later re-accreted. The recycling timescale is typically about 1 Gyr, virtually independent of halo mass. Recycled material is re-accreted farther out in the disk and with typically about 2-3 times more angular momentum. These results elucidate and quantify how the baryon cycle plausibly regulates star formation and alters the angular momentum distribution of disk material across the halo mass range where most of cosmic star formation occurs.

We present a self-consistent prediction from a large-scale cosmological simulation for the population of `wandering' supermassive black holes (SMBHs) of mass greater than $10^6$ M$_{\odot}$ on long-lived, kpc-scale orbits within Milky Way (MW)-mass galaxies. We extract a sample of MW-mass halos from the Romulus25 cosmological simulation (Tremmel et al. 2017), which is uniquely able to capture the orbital evolution of SMBHs during and following galaxy mergers. We predict that such halos, regardless of recent merger history or morphology, host an average of $5.1 \pm 3.3$ SMBHs, including their central black hole, within 10 kpc from the galactic center and an average of $12.2 \pm 8.4$ SMBHs total within their virial radius, not counting those in satellite halos. Wandering SMBHs exist within their host galaxies for several Gyrs, often accreted by their host halo in the early Universe. We find, with $>4\sigma$ significance, that wandering SMBHs are preferentially found outside of galactic disks.

Correlations measured in three dimensions in the Lyman-alpha forest are contaminated by the presence of the damping wings of high column density (HCD) absorbing systems of neutral hydrogen (HI; having column densities $N(\mathrm{HI}) > 1.6\times10^{17}\,\mathrm{atoms}\,\mathrm{cm}^{-2}$), which extend significantly beyond the redshift-space location of the absorber. We measure this effect as a function of the column density of the HCD absorbers and redshift by measuring 3D flux power spectra in cosmological hydrodynamical simulations from the Illustris project. Survey pipelines exclude regions containing the largest damping wings. We find that, even after this procedure, there is a scale-dependent correction to the 3D Lyman-alpha forest flux power spectrum from residual contamination. We model this residual using a simple physical model of the HCD absorbers as linearly biased tracers of the matter density distribution, convolved with their Voigt profiles and integrated over the column density distribution function. We recommend the use of this model over existing models used in data analysis, which approximate the damping wings as top-hats and so miss shape information in the extended wings. The simple 'linear Voigt model' is statistically consistent with our simulation results for a mock residual contamination up to small scales ($|k| < 1\,h\,\mathrm{Mpc}^{-1}$). It does not account for the effect of the highest column density absorbers on the smallest scales (e.g., $|k| > 0.4\,h\,\mathrm{Mpc}^{-1}$ for small damped Lyman-alpha absorbers; HCD absorbers with $N(\mathrm{HI}) \sim 10^{21}\,\mathrm{atoms}\,\mathrm{cm}^{-2}$). However, these systems are in any case preferentially removed from survey data. Our model is appropriate for an accurate analysis of the baryon acoustic oscillations feature. It is additionally essential for reconstructing the full shape of the 3D flux power spectrum.

We present the first self-consistent prediction for the distribution of formation timescales for close Supermassive Black Hole (SMBH) pairs following galaxy mergers. Using ROMULUS25, the first large-scale cosmological simulation to accurately track the orbital evolution of SMBHs within their host galaxies down to sub-kpc scales, we predict an average formation rate density of close SMBH pairs of 0.013 cMpc^-3 Gyr^-1. We find that it is relatively rare for galaxy mergers to result in the formation of close SMBH pairs with sub-kpc separation and those that do form are often the result of Gyrs of orbital evolution following the galaxy merger. The likelihood and timescale to form a close SMBH pair depends strongly on the mass ratio of the merging galaxies, as well as the presence of dense stellar cores. Low stellar mass ratio mergers with galaxies that lack a dense stellar core are more likely to become tidally disrupted and deposit their SMBH at large radii without any stellar core to aid in their orbital decay, resulting in a population of long-lived 'wandering' SMBHs. Conversely, SMBHs in galaxies that remain embedded within a stellar core form close pairs in much shorter timescales on average. This timescale is a crucial, though often ignored or very simplified, ingredient to models predicting SMBH mergers rates and the connection between SMBH and star formation activity.

We measure the effect of high column density absorbing systems of neutral hydrogen (HI) on the one-dimensional (1D) Lyman-alpha forest flux power spectrum using cosmological hydrodynamical simulations from the Illustris project. High column density absorbers (which we define to be those with HI column densities $N(\mathrm{HI}) > 1.6 \times 10^{17}\,\mathrm{atoms}\,\mathrm{cm}^{-2}$) cause broadened absorption lines with characteristic damping wings. These damping wings bias the 1D Lyman-alpha forest flux power spectrum by causing absorption in quasar spectra away from the location of the absorber itself. We investigate the effect of high column density absorbers on the Lyman-alpha forest using hydrodynamical simulations for the first time. We provide templates as a function of column density and redshift, allowing the flexibility to accurately model residual contamination, i.e., if an analysis selectively clips out the largest damping wings. This flexibility will improve cosmological parameter estimation, e.g., allowing more accurate measurement of the shape of the power spectrum, with implications for cosmological models containing massive neutrinos or a running of the spectral index. We provide fitting functions to reproduce these results so that they can be incorporated straightforwardly into a data analysis pipeline.

Recent work has studied the interplay between a galaxy's history and its observable properties using "genetically modified" cosmological zoom simulations. The approach systematically generates alternative histories for a halo, while keeping its cosmological environment fixed. Applications to date altered linear properties of the initial conditions such as the mean overdensity of specified regions; we extend the formulation to include quadratic features such as local variance, which determines the overall importance of smooth accretion relative to mergers in a galaxy's history. We introduce an efficient algorithm for this new class of modification and demonstrate its ability to control the variance of a region in a one-dimensional toy model. Outcomes of this work are two-fold: (i) a clarification of the formulation of genetic modifications and (ii) a proof of concept for quadratic modifications leading the way to a forthcoming implementation in cosmological simulations.

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.

When modelling the ionisation of gas in the intergalactic medium after reionisation, it is standard practice to assume a uniform radiation background. This assumption is not always appropriate; models with radiative transfer show that large-scale ionisation rate fluctuations can have an observable impact on statistics of the Lyman-alpha forest. We extend such calculations to include beaming of sources, which has previously been neglected but which is expected to be important if quasars dominate the ionising photon budget. Beaming has two effects: first, the physical number density of ionising sources is enhanced relative to that directly observed; and second, the radiative transfer itself is altered. We calculate both effects in a hard-edged beaming model where each source has a random orientation, using an equilibrium Boltzmann hierarchy in terms of spherical harmonics. By studying the statistical properties of the resulting ionisation rate and HI density fields at redshift $z\sim 2.3$, we find that the two effects partially cancel each other; combined, they constitute a maximum $5\%$ correction to the power spectrum $P_{\mathrm{HI}}(k)$ at $k=0.04 \, h/\mathrm{Mpc}$. On very large scales ($k<0.01\, h/\mathrm{Mpc}$) the source density renormalisation dominates; it can reduce, by an order of magnitude, the contribution of ionising shot-noise to the intergalactic HI power spectrum. The effects of beaming should be considered when interpreting future observational datasets.

We present a novel implementation of supermassive black hole (SMBH) formation, dynamics, and accretion in the massively parallel tree+SPH code, ChaNGa. This approach improves the modeling of SMBHs in fully cosmological simulations, allowing for a more de- tailed analysis of SMBH-galaxy co-evolution throughout cosmic time. Our scheme includes novel, physically motivated models for SMBH formation, dynamics and sinking timescales within galaxies, and SMBH accretion of rotationally supported gas. The sub-grid parameters that regulate star formation (SF) and feedback from SMBHs and SNe are optimized against a comprehensive set of z = 0 galaxy scaling relations using a novel, multi-dimensional parameter search. We have incorporated our new SMBH implementation and parameter optimization into a new set of high resolution, large-scale cosmological simulations called Romulus. We present initial results from our flagship simulation, Romulus25, showing that our SMBH model results in SF efficiency, SMBH masses, and global SF and SMBH accretion histories at high redshift that are consistent with observations. We discuss the importance of SMBH physics in shaping the evolution of massive galaxies and show how SMBH feedback is much more effective at regulating star formation compared to SNe feedback in this regime. Further, we show how each aspect of our SMBH model impacts this evolution compared to more common approaches. Finally, we present a science application of this scheme studying the properties and time evolution of an example dual AGN system, highlighting how our approach allows simulations to better study galaxy interactions and SMBH mergers in the context of galaxy-BH co-evolution.

We explore for the first time the effect of self-interacting dark matter (SIDM) on the dark matter (DM) and baryonic distribution in massive galaxies formed in hydrodynamical cosmological simulations, including explicit baryonic physics treatment. A novel implementation of Super-Massive Black Hole (SMBH) formation and evolution is used, as in Tremmel et al.(2015, 2016), allowing to explicitly follow SMBH dynamics at the center of galaxies. A high SIDM constant cross-section is chosen, $\sigma$=10 $\rm cm^2/gr$, to amplify differences from CDM models. Milky Way-like galaxies form a shallower DM density profile in SIDM than they do in CDM, with differences already at 20 kpc scales. This demonstrates that even for the most massive spirals the effect of SIDM dominates over the adiabatic contraction due to baryons. Strikingly, the dynamics of SMBHs differs in the SIDM and reference CDM case. SMBHs in massive spirals have sunk to the centre of their host galaxy in both the SIDM and CDM run, while in less massive galaxies about 80$\%$ of the SMBH population is off-centered in the SIDM case, as opposed to the CDM case in which $\sim$90$\%$ of SMBHs have reached their host's centre. SMBHs are found as far as $\sim$9 kpc away from the centre of their host SIDM galaxy. This difference is due to the increased dynamical friction timescale caused by the lower DM density in SIDM galaxies compared to CDM, resulting in 'core stalling'. This pilot work highlights the importance of simulating in a full hydrodynamical context different DM models combined to SMBH physics to study their influence on galaxy formation.

We show how the interplay between active galactic nuclei (AGN) and merger history determines whether a galaxy quenches star formation at high redshift. We first simulate, in a full cosmological context, a galaxy of total dynamical mass $10^{12}\,M_{\odot}$ at $z=2$. Then we systematically alter the accretion history of the galaxy by minimally changing the linear overdensity in the initial conditions. This "genetic modification" approach allows the generation of three sets of $\Lambda$CDM initial conditions leading to maximum merger ratios of 1:10, 1:5 and 2:3 respectively. The changes leave the final halo mass, large scale structure and local environment unchanged, providing a controlled numerical experiment. Interaction between the AGN physics and mergers in the three cases lead respectively to a star-forming, temporarily-quenched and permanently-quenched galaxy. However the differences do not primarily lie in the black hole accretion rates, but in the kinetic effects of the merger: the galaxy is resilient against AGN feedback unless its gaseous disk is first disrupted. Typical accretion rates are comparable in the three cases, falling below $0.1\,M_{\odot}$ yr$^{-1}$, equivalent to around $2\%$ of the Eddington rate or $10^{-3}$ times the pre-quenching star formation rate, in agreement with observations. This low level of black hole accretion can be sustained even when there is insufficient dense cold gas for star formation. Conversely, supernova feedback is too distributed to generate outflows in high-mass systems, and cannot maintain quenching over periods longer than the halo gas cooling time.

A fundamental assumption in the standard model of cosmology is that the Universe is isotropic on large scales. Breaking this assumption leads to a set of solutions to Einstein's field equations, known as Bianchi cosmologies, only a subset of which have ever been tested against data. For the first time, we consider all degrees of freedom in these solutions to conduct a general test of isotropy using cosmic microwave background temperature and polarization data from Planck. For the vector mode (associated with vorticity), we obtain a limit on the anisotropic expansion of $(\sigma_V/H)_0 < 4.7 \times 10^{-11}$ (95% CI), which is an order of magnitude tighter than previous Planck results that used CMB temperature only. We also place upper limits on other modes of anisotropic expansion, with the weakest limit arising from the regular tensor mode, $(\sigma_{T,\rm reg}/H)_0<1.0 \times 10^{-6}$ (95% CI). Including all degrees of freedom simultaneously for the first time, anisotropic expansion of the Universe is strongly disfavoured, with odds of 121,000:1 against.

We present a new framework for testing the isotropy of the Universe using cosmic microwave background data, building on the nested-sampling ANICOSMO code. Uniquely, we are able to constrain the scalar, vector and tensor degrees of freedom alike; previous studies only considered the vector mode (linked to vorticity). We employ Bianchi type VII$_h$ cosmologies to model the anisotropic Universe, from which other types may be obtained by taking suitable limits. In a separate development, we improve the statistical analysis by including the effect of Bianchi power in the high-$\ell$, as well as the low-$\ell$, likelihood. To understand the effect of all these changes, we apply our new techniques to WMAP data. We find no evidence for anisotropy, constraining shear in the vector mode to $(\sigma_V/H)_0 < 1.7 \times 10^{-10}$ (95% CL). For the first time, we place limits on the tensor mode; unlike other modes, the tensor shear can grow from a near-isotropic early Universe. The limit on this type of shear is $(\sigma_{T,\rm reg}/H)_0 < 2.4 \times 10^{-7}$ (95% CL).

We present Spin-SILC, a new foreground component separation method that accurately extracts the cosmic microwave background (CMB) polarisation $E$ and $B$ modes from raw multifrequency Stokes $Q$ and $U$ measurements of the microwave sky. Spin-SILC is an internal linear combination method that uses spin wavelets to analyse the spin-2 polarisation signal $P = Q + iU$. The wavelets are additionally directional (non-axisymmetric). This allows different morphologies of signals to be separated and therefore the cleaning algorithm is localised using an additional domain of information. The advantage of spin wavelets over standard scalar wavelets is to simultaneously and self-consistently probe scales and directions in the polarisation signal $P = Q + iU$ and in the underlying $E$ and $B$ modes, therefore providing the ability to perform component separation and $E$-$B$ decomposition concurrently for the first time. We test Spin-SILC on full-mission Planck simulations and data and show the capacity to correctly recover the underlying cosmological $E$ and $B$ modes. We also demonstrate a strong consistency of our CMB maps with those derived from existing component separation methods. Spin-SILC can be combined with the pseudo- and pure $E$-$B$ spin wavelet estimators presented in a companion paper to reliably extract the cosmological signal in the presence of complicated sky cuts and noise. Therefore, it will provide a computationally-efficient method to accurately extract the CMB $E$ and $B$ modes for future polarisation experiments.

We study how outflows of gas launched from a central galaxy undergoing repeated starbursts propagate through the circumgalactic medium (CGM), using the simulation code RAMSES. We assume that the outflow from the disk can be modelled as a rapidly moving bubble of hot gas at $\mathrm{\sim1\;kpc}$ above disk, then ask what happens as it moves out further into the halo around the galaxy on $\mathrm{\sim 100\;kpc}$ scales. To do this we run 60 two-dimensional simulations scanning over parameters of the outflow. Each of these is repeated with and without radiative cooling, assuming a primordial gas composition to give a lower bound on the importance of cooling. In a large fraction of radiative-cooling cases we are able to form rapidly outflowing cool gas from in situ cooling of the flow. We show that the amount of cool gas formed depends strongly on the 'burstiness' of energy injection; sharper, stronger bursts typically lead to a larger fraction of cool gas forming in the outflow. The abundance ratio of ions in the CGM may therefore change in response to the detailed historical pattern of star formation. For instance, outflows generated by star formation with short, intense bursts contain up to 60 per cent of their gas mass at temperatures $<5 \times 10^4\,\mathrm{K}$; for near-continuous star formation the figure is $\lesssim$ 5 per cent. Further study of cosmological simulations, and of idealised simulations with e.g., metal-cooling, magnetic fields and/or thermal conduction, will help to understand the precise signature of bursty outflows on observed ion abundances.

We present and test a method that dramatically reduces variance arising from the sparse sampling of wavemodes in cosmological simulations. The method uses two simulations which are fixed (the initial Fourier mode amplitudes are fixed to the ensemble average power spectrum) and paired (with initial modes exactly out of phase). We measure the power spectrum, monopole and quadrupole redshift-space correlation functions, halo mass function and reduced bispectrum at $z=1$. By these measures, predictions from a fixed pair can be as precise on non-linear scales as an average over 50 traditional simulations. The fixing procedure introduces a non-Gaussian correction to the initial conditions; we give an analytic argument showing why the simulations are still able to predict the mean properties of the Gaussian ensemble. We anticipate that the method will drive down the computational time requirements for accurate large-scale explorations of galaxy bias and clustering statistics, enabling more precise comparisons with theoretical models, and facilitating the use of numerical simulations in cosmological data interpretation.

We present new clean maps of the CMB temperature anisotropies (as measured by Planck) constructed with a novel internal linear combination (ILC) algorithm using directional, scale-discretised wavelets --- Scale-discretised, directional wavelet ILC or SILC. Directional wavelets, when convolved with signals on the sphere, can separate the anisotropic filamentary structures which are characteristic of both the CMB and foregrounds. Extending previous component separation methods, which use the frequency, spatial and harmonic signatures of foregrounds to separate them from the cosmological background signal, SILC can additionally use morphological information in the foregrounds and CMB to better localise the cleaning algorithm. We test the method on Planck data and simulations, demonstrating consistency with existing component separation algorithms, and discuss how to optimise the use of morphological information by varying the number of directional wavelets as a function of spatial scale. We find that combining the use of directional and axisymmetric wavelets depending on scale could yield higher quality CMB temperature maps. Our results set the stage for the application of SILC to polarisation anisotropies through an extension to spin wavelets.

We introduce and explore "paired" cosmological simulations. A pair consists of an A and B simulation with initial conditions related by the inversion $\delta_A(x, t_{initial})=-\delta_B(x,t_{initial})$ (underdensities substituted for overdensities and vice versa). We argue that the technique is valuable for improving our understanding of cosmic structure formation. The A and B fields are by definition equally likely draws from {\Lambda}CDM initial conditions, and in the linear regime evolve identically up to the overall sign. As non-linear evolution takes hold, a region that collapses to form a halo in simulation A will tend to expand to create a void in simulation B. Applications include (i) contrasting the growth of A-halos and B-voids to test excursion-set theories of structure formation; (ii) cross-correlating the density field of the A and B universes as a novel test for perturbation theory; and (iii) canceling error terms by averaging power spectra between the two boxes. Generalizations of the method to more elaborate field transformations are suggested.

We propose a method to generate `genetically-modified' (GM) initial conditions for high-resolution simulations of galaxy formation in a cosmological context. Building on the Hoffman-Ribak algorithm, we start from a reference simulation with fully random initial conditions, then make controlled changes to specific properties of a single halo (such as its mass and merger history). The algorithm demonstrably makes minimal changes to other properties of the halo and its environment, allowing us to isolate the impact of a given modification. As a significant improvement over previous work, we are able to calculate the abundance of the resulting objects relative to the reference simulation. Our approach can be applied to a wide range of cosmic structures and epochs; here we study two problems as a proof-of-concept. First, we investigate the change in density profile and concentration as the collapse time of three individual halos are varied at fixed final mass, showing good agreement with previous statistical studies using large simulation suites. Second, we modify the $z=0$ mass of halos to show that our theoretical abundance calculations correctly recover the halo mass function. The results demonstrate that the technique is robust, opening the way to controlled experiments in galaxy formation using hydrodynamic zoom simulations.

We use cosmological hydrodynamic simulations to consistently compare the assembly of dwarf galaxies in both $\Lambda$ dominated, Cold (CDM) and Self--Interacting (SIDM) dark matter models. The SIDM model adopts a constant cross section of 2 $cm^{2}/g$, a relatively large value to maximize its effects. These are the first SIDM simulations that are combined with a description of stellar feedback that naturally drives potential fluctuations able to create dark matter cores. Remarkably, SIDM fails to significantly lower the central dark matter density at halo peak velocities V$_{max}$ $<$ 30 Km/s. This is due to the fact that the central regions of very low--mass field halos have relatively low central velocity dispersion and densities, leading to time scales for SIDM collisions greater than a Hubble time. CDM halos with V$_{max}$ $<$ 30 km/s have inefficient star formation, and hence weak supernova feedback. At a fixed 2 cm2/g SIDM cross section, the DM content of very low mass CDM and SIDM halos differs by no more than a factor of two within 100-200pc. At larger halo masses ($\sim$ 10$^{10}$ solar masses), the introduction of baryonic processes creates field dwarf galaxies with dark matter cores and central DM$+$baryon distributions that are effectively indistinguishable between CDM and SIDM. Both models are in broad agreement with observed Local Group field galaxies across the range of masses explored. To significantly differentiate SIDM from CDM at the scale of faint dwarf galaxies, a velocity dependent cross section that rapidly increases to values larger than 2 $cm^{2}/g$ for halos with V$_{max}$ < 25-30 Km/s needs to be introduced.

Galaxies and the dark matter halos that host them are not spherically symmetric, yet spherical symmetry is a helpful simplifying approximation for idealised calculations and analysis of observational data. The assumption leads to an exact conservation of angular momentum for every particle, making the dynamics unrealistic. But how much does that inaccuracy matter in practice for analyses of stellar distribution functions, collisionless relaxation, or dark matter core-creation? We provide a general answer to this question for a wide class of aspherical systems; specifically, we consider distribution functions that are "maximally stable", i.e. that do not evolve at first order when external potentials (which arise from baryons, large scale tidal fields or infalling substructure) are applied. We show that a spherically-symmetric analysis of such systems gives rise to the false conclusion that the density of particles in phase space is ergodic (a function of energy alone). Using this idea we are able to demonstrate that: (a) observational analyses that falsely assume spherical symmetry are made more accurate by imposing a strong prior preference for near-isotropic velocity dispersions in the centre of spheroids; (b) numerical simulations that use an idealised spherically-symmetric setup can yield misleading results and should be avoided where possible; and (c) triaxial dark matter halos (formed in collisionless cosmological simulations) nearly attain our maximally-stable limit, but their evolution freezes out before reaching it.

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.

The archetypal theory of dark energy is quintessence: a minimally coupled scalar field with a canonical kinetic energy and potential. By studying random potentials we show that quintessence imposes a restricted set of priors on the equation of state of dark energy. Focusing on the commonly-used parametrisation, $w(a)\approx w_0+w_a(1-a)$, we show that there is a natural scale and direction in the $(w_0, w_a)$ plane that distinguishes quintessence as a general framework. We calculate the expected information gain for a given survey and show that, because of the non-trivial prior information, it is a function of more than just the figure of merit. This allows us to make a quantitative case for novel survey strategies. We show that the scale of the prior sets target observational requirements for gaining significant information. This corresponds to a figure of merit FOM$\gtrsim 200$, a requirement that future galaxy redshift surveys will meet.

We use high resolution Hydro$+$N-Body cosmological simulations to compare the assembly and evolution of a small field dwarf (stellar mass ~ 10$^{6-7}$ M$\odot$, total mass 10$^{10}$ M$\odot$ in $\Lambda$ dominated CDM and 2keV WDM cosmologies. We find that star formation (SF) in the WDM model is reduced and delayed by 1-2 Gyr relative to the CDM model, independently of the details of SF and feedback. Independent of the DM model, but proportionally to the SF efficiency, gas outflows lower the central mass density through `dynamical heating', such that all realizations have circular velocities $<$ 20kms at 500$~$pc, in agreement with local kinematic constraints. As a result of dynamical heating, older stars are less centrally concentrated than younger stars, similar to stellar population gradients observed in nearby dwarf galaxies. Introducing an important diagnostic of SF and feedback models, we translate our simulations into artificial color-magnitude diagrams and star formation histories in order to directly compare to available observations. The simulated galaxies formed most of their stars in many $\sim$10 Myr long bursts. The CDM galaxy has a global SFH, HI abundance and Fe/H and alpha-elements distribution well matched to current observations of dwarf galaxies. These results highlight the importance of directly including `baryon physics' in simulations when 1) comparing predictions of galaxy formation models with the kinematics and number density of local dwarf galaxies and 2) differentiating between CDM and non-standard models with different DM or power spectra.

Recent work has shown that the z~2.5 Lyman-alpha forest on large scales encodes information about the galaxy and quasar populations that keep the intergalactic medium photoionized. We present the first forecasts for constraining the populations with data from current and next-generation surveys. At a minimum the forest should tell us whether galaxies or, conversely, quasars dominate the photon production. The number density and clustering strength of the ionising sources might be estimated to sub-10% precision with a DESI-like survey if degeneracies (e.g., with the photon mean-free-path, small-scale clustering power normalization and potentially other astrophysical effects) can be broken by prior information. We demonstrate that, when inhomogeneous ionisation is correctly handled, constraints on dark energy do not degrade.

I discuss fluctuations in the neutral hydrogen (HI) density of the z~2.3 intergalactic medium and show that their relation to cosmic overdensity is strongly scale-dependent. This behaviour arises from a linearized version of the well-known "proximity effect", in which bright sources suppress atomic hydrogen density. Using a novel, systematic and detailed linear-theory radiative transfer calculation, I demonstrate how HI density consequently anti-correlates with total matter density when averaged on scales exceeding the Lyman-limit mean-free-path. The radiative transfer thumbprint is highly distinctive and should be measurable in the Lyman-alpha forest. Effects extend to sufficiently small scales to generate significant distortion of the correlation function shape around the baryon acoustic oscillation peak, although the peak location shifts only by 1.2 percent for a mean source bias of b_j=3. The distortion changes significantly with b_j and other astrophysical parameters; measuring it should provide a helpful observational constraint on the nature of ionizing photon sources in the near future.

One of the principal discoveries in modern cosmology is that standard model particles (including baryons, leptons and photons) together comprise only 5% of the mass-energy budget of the Universe. The remaining 95% consists of dark energy and dark matter (DM). Consequently our picture of the universe is known as {\Lambda}CDM, with {\Lambda} denoting dark energy and CDM cold dark matter. {\Lambda}CDM is being challenged by its apparent inability to explain the low density of DM measured at the centre of cosmological systems, ranging from faint dwarf galaxies to massive clusters containing tens of galaxies the size of the Milky Way. But before making conclusions one should carefully include the effect of gas and stars, which were historically seen as merely a passive component during the assembly of galaxies. We now understand that these can in fact significantly alter the DM component, through a coupling based on rapid gravitational potential fluctuations.

The angular power spectrum is a powerful statistic for analysing cosmological signals imprinted in the clustering of matter. However, current galaxy and quasar surveys cover limited portions of the sky, and are contaminated by systematics that can mimic cosmological signatures and jeopardise the interpretation of the measured power spectra. We provide a framework for obtaining unbiased estimates of the angular power spectra of large-scale structure surveys at the largest scales using quadratic estimators. The method is tested by analysing the 600 CMASS mock catalogues constructed by Manera et al. (2013) for the Baryon Oscillation Spectroscopic Survey (BOSS). We then consider the Richards et al. (2009) catalogue of photometric quasars from the Sixth Data Release (DR6) of the Sloan Digital Sky Survey (SDSS), which is known to include significant stellar contamination and systematic uncertainties. Focusing on the sample of ultraviolet-excess (UVX) sources, we show that the excess clustering power present on the largest-scales can be largely mitigated by making use of improved sky masks and projecting out the modes corresponding to the principal systematics. In particular, we find that the sample of objects with photometric redshift $1.3 < z_p < 2.2$ exhibits no evidence of contamination when using our most conservative mask and mode projection. This indicates that any residual systematics are well within the statistical uncertainties. We conclude that, using our approach, this sample can be used for cosmological studies.

We present the first public version (v0.2) of the open-source and community-developed Python package, Astropy. This package provides core astronomy-related functionality to the community, including support for domain-specific file formats such as Flexible Image Transport System (FITS) files, Virtual Observatory (VO) tables, and common ASCII table formats, unit and physical quantity conversions, physical constants specific to astronomy, celestial coordinate and time transformations, world coordinate system (WCS) support, generalized containers for representing gridded as well as tabular data, and a framework for cosmological transformations and conversions. Significant functionality is under active development, such as a model fitting framework, VO client and server tools, and aperture and point spread function (PSF) photometry tools. The core development team is actively making additions and enhancements to the current code base, and we encourage anyone interested to participate in the development of future Astropy versions.

Chronos is our response to ESA's call for white papers to define the science for the future L2, L3 missions. Chronos targets the formation and evolution of galaxies, by collecting the deepest NIR spectroscopic data, from the formation of the first galaxies at z~10 to the peak of formation activity at z~1-3. The strong emission from the atmospheric background makes this type of survey impossible from a ground-based observatory. The spectra of galaxies represent the equivalent of a DNA fingerprint, containing information about the past history of star formation and chemical enrichment. The proposed survey will allow us to dissect the formation process of galaxies including the timescales of quenching triggered by star formation or AGN activity, the effect of environment, the role of infall/outflow processes, or the connection between the galaxies and their underlying dark matter haloes. To provide these data, the mission requires a 2.5m space telescope optimised for a campaign of very deep NIR spectroscopy. A combination of a high multiplex and very long integration times will result in the deepest, largest, high-quality spectroscopic dataset of galaxies from z=1 to 12, spanning the history of the Universe, from 400 million to 6 billion years after the big bang, i.e. covering the most active half of cosmic history.

The presence of a dark matter core in the central kiloparsec of many dwarf galaxies has been a long standing problem in galaxy formation theories based on the standard cold dark matter paradigm. Recent cosmological simulations, based on Smooth Particle Hydrodynamics and rather strong feedback recipes have shown that it was indeed possible to form extended dark matter cores using baryonic processes related to a more realistic treatment of the interstellar medium. Using adaptive mesh refinement, together with a new, stronger supernovae feedback scheme that we have recently implemented in the RAMSES code, we show that it is also possible to form a prominent dark matter core within the well-controlled framework of an isolated, initially cuspy, 10 billion solar masses dark matter halo. Although our numerical experiment is idealized, it allows a clean and unambiguous identification of the dark matter core formation process. Our dark matter inner profile is well fitted by a pseudo-isothermal profile with a core radius of 800 pc. The core formation mechanism is consistent with the one proposed recently by Pontzen & Governato. We highlight two key observational predictions of all simulations that find cusp-core transformations: (i) a bursty star formation history with peak to trough ratio of 5 to 10 and a duty cycle comparable to the local dynamical time; and (ii) a stellar distribution that is hot with v/sigma=1. We compare the observational properties of our model galaxy with recent measurements of the isolated dwarf WLM. We show that the spatial and kinematical distribution of stars and HI gas are in striking agreement with observations, supporting the fundamental role played by stellar feedback in shaping both the stellar and dark matter distribution.

Using high resolution cosmological hydrodynamical simulations of Milky Way-massed disk galaxies, we demonstrate that supernovae feedback and tidal stripping lower the central masses of bright (-15 < M_V < -8) satellite galaxies. These simulations resolve high density regions, comparable to giant molecular clouds, where stars form. This resolution allows us to adopt a prescription for H_2 formation and destruction that ties star formation to the presence of shielded, molecular gas. Before infall, supernova feedback from the clumpy, bursty star formation captured by this physically motivated model leads to reduced dark matter (DM) densities and shallower inner density profiles in the massive satellite progenitors (Mvir > 10^9 Msun, Mstar > 10^7 Msun) compared to DM-only simulations. The progenitors of the lower mass satellites are unable to maintain bursty star formation histories, due to both heating at reionization and gas loss from initial star forming events, preserving the steep inner density profile predicted by DM-only simulations. After infall, tidal stripping acts to further reduce the central densities of the luminous satellites, particularly those that enter with cored dark matter halos, increasing the discrepancy in the central masses predicted by baryon+DM and DM-only simulations. We show that DM-only simulations, which neglect the baryonic effects described in this work, produce denser satellites with larger central velocities. We provide a simple correction to the central DM mass predicted for satellites by DM-only simulations. We conclude that DM-only simulations should be used with great caution when interpreting kinematic observations of the Milky Way's dwarf satellites.

We use maximum entropy arguments to derive the phase space distribution of a virialized dark matter halo. Our distribution function gives an improved representation of the end product of violent relaxation. This is achieved by incorporating physically motivated dynamical constraints (specifically on orbital actions) which prevent arbitrary redistribution of energy. We compare the predictions with three high-resolution dark matter simulations of widely varying mass. The numerical distribution function is accurately predicted by our argument, producing an excellent match for the vast majority of particles. The remaining particles constitute the central cusp of the halo (<4% of the dark matter). They can be accounted for within the presented framework once the short dynamical timescales of the centre are taken into account.

We explore the circumgalactic medium (CGM) of two simulated star-forming galaxies with luminosities L ~ 0.1 and 1 L* generated using the smooth particle hydrodynamic code GASOLINE. These simulations are part of the Making Galaxies In a Cosmological Context (MAGICC) program in which the stellar feedback is tuned to match the stellar mass-halo mass relationship. For comparison, each galaxy was also simulated using a 'lower feedback' (LF) model which has strength comparable to other implementations in the literature. The 'MAGICC feedback' (MF) model has a higher incidence of massive stars and an approximately two times higher energy input per supernova. Apart from the low-mass halo using LF, each galaxy exhibits a metal-enriched CGM that extends to approximately the virial radius. A significant fraction of this gas has been heated in supernova explosions in the disc and subsequently ejected into the CGM where it is predicted to give rise to substantial O VI absorption. The simulations do not yet address the question of what happens to the O VI when the galaxies stop forming stars. Our models also predict a reservoir of cool H I clouds that show strong Ly\alpha absorption to several hundred kpc. Comparing these models to recent surveys with the Hubble Space Telescope, we find that only the MF models have sufficient O VI and H I gas in the CGM to reproduce the observed distributions. In separate analyses, these same MF models also show better agreement with other galaxy observables (e.g. rotation curves, surface brightness profiles and H I gas distribution). We infer that the CGM is the dominant reservoir of baryons for galaxy haloes.

We calculate the energy that baryons must inject in cold dark matter (CDM) haloes in order to remove centrally-divergent DM cusps on scales relevant to observations of dwarf spheroidal galaxies (dSphs). We estimate that the CDM haloes often associated with the Milky Way's dSphs (M_vir/M_\odot \sim 10^{9-10}) require \Delta E/erg \sim 10^{53-55} in order to form cores on scales comparable to the luminous size of these galaxies. While supernova type II (SNeII) explosions can in principle generate this energy, the actual contribution is limited by the low star formation efficiency implied by the abundance of luminous satellites. Considering that CDM's well-known `core/cusp' and `missing satellite' problems place opposing demands on star formation efficiencies, existing observational evidences for large cores in the most luminous dSphs require that CDM models invoke some combination of the following: (i) efficient (of order unity) coupling of SNeII energy into dark matter particles, (ii) star formation histories peaking at unexpectedly high redshifts (z>6), (iii) a top-heavy stellar IMF, and/or (iv) substantial satellite disruption or other stochastic effects to ease the substructure abundance constraints. Our models show that the tension between CDM problems on small scales would increase if cored DM profiles were to be found in fainter dwarves.

We examine the evolution of the inner dark matter (DM) and baryonic density profile of a new sample of simulated field galaxies using fully cosmological, Lambda CDM, high resolution SPH + N-Body simulations. These simulations include explicit H2 and metal cooling, star formation (SF) and supernovae (SNe) driven gas outflows. Starting at high redshift, rapid, repeated gas outflows following bursty SF transfer energy to the DM component and significantly flatten the originally `cuspy' central DM mass profile of galaxies with present day stellar masses in the 10^4.5 -- 10^9.8 Msolar range. At z=0, the central slope of the DM density profile of our galaxies (measured between 0.3 and 0.7 kpc from their centre) is well fitted by rhoDM propto r^alpha with alpha \simeq -0.5 + 0.35 log_10(Mstar/10^8Msolar) where Mstar is the stellar mass of the galaxy and 4 < log_10 Mstar < 9.4. These values imply DM profiles flatter than those obtained in DM--only simulations and in close agreement with those inferred in galaxies from the THINGS and LITTLE THINGS survey. Only in very small halos, where by z=0 star formation has converted less than ~ 0.03% of the original baryon abundance into stars, outflows do not flatten the original cuspy DM profile out to radii resolved by our simulations. The mass (DM and baryonic) measured within the inner 500 pc of each simulated galaxy remains nearly constant over four orders of magnitudes in stellar mass for Mstar 10^9 Msolar. This finding is consistent with estimates for faint Local Group dwarfs and field galaxies. These results address one of the outstanding problems faced by the CDM model, namely the strong discrepancy between the original predictions of cuspy DM profiles and the shallower central DM distribution observed in galaxies.

We propose and successfully test against new cosmological simulations a novel analytical description of the physical processes associated with the origin of cored dark matter density profiles. In the simulations, the potential in the central kiloparsec changes on sub-dynamical timescales over the redshift interval 4 > z > 2 as repeated, energetic feedback generates large underdense bubbles of expanding gas from centrally-concentrated bursts of star formation. The model demonstrates how fluctuations in the central potential irreversibly transfer energy into collisionless particles, thus generating a dark matter core. A supply of gas undergoing collapse and rapid expansion is therefore the essential ingredient. The framework, based on a novel impulsive approximation, breaks with the reliance on adiabatic approximations which are inappropriate in the rapidly-changing limit. It shows that both outflows and galactic fountains can give rise to cusp-flattening, even when only a few per cent of the baryons form stars. Dwarf galaxies maintain their core to the present time. The model suggests that constant density dark matter cores will be generated in systems of a wide mass range if central starbursts or AGN phases are sufficiently frequent and energetic.

Obscuration due to Galactic emission complicates the extraction of information from cosmological surveys, and requires some combination of the (typically imperfect) modeling and subtraction of foregrounds, or the removal of part of the sky. This particularly affects the extraction of information from the largest observable scales. Maximum-likelihood estimators for reconstructing the full-sky spherical harmonic coefficients from partial-sky maps have recently been shown to be susceptible to contamination from within the sky cut, arising due to the necessity to band-limit the data by smoothing prior to reconstruction. Using the WMAP 7-year data, we investigate modified implementations of such estimators which are robust to the leakage of contaminants from within masked regions. We provide a measure, based on the expected amplitude of residual foregrounds, for selecting the most appropriate estimator for the task at hand. We explain why the related quadratic maximum-likelihood estimator of the angular power spectrum does not suffer from smoothing-induced bias.

Using high resolution, fully cosmological smoothed particle hydro-dynamical simulations of dwarf galaxies in a Lambda cold dark matter Universe, we show how baryons attain a final angular momentum distribution which allows pure disc galaxies to form. Blowing out substantial amounts of gas through supernovae and stellar winds, which is well supported observationally, is a key ingredient in forming bulgeless discs. We outline why galactic outflows preferentially remove low angular momentum material, and show that this is a natural result when structure forms in a cold dark matter cosmology. The driving factors are a) the mean angular momentum of accreted material increases with time, b) lower potentials at early times, c) the existence of an extended reservoir of high angular momentum gas which is not within star forming regions, meaning that only gas from the inner region (low angular momentum gas) is expelled and d) the tendency for outflows to follow the path of least resistance which is perpendicular to the disc. We also show that outflows are enhanced during mergers, thus expelling much of the gas which has lost its angular momentum during these events, and preventing the formation of "classical", merger driven bulges in low mass systems. Stars formed prior to such mergers form a diffuse, extended stellar halo component.

Homogeneous, nearly-isotropic Bianchi cosmological models are considered. Their time evolution is expressed as a complete set of non-interacting linear modes on top of a Friedmann-Robertson-Walker background model. This connects the extensive literature on Bianchi models with the more commonly-adopted perturbation approach to general relativistic cosmological evolution. Expressions for the relevant metric perturbations in familiar coordinate systems can be extracted straightforwardly. Amongst other possibilities, this allows for future analysis of anisotropic matter sources in a more general geometry than usually attempted. We discuss the geometric mechanisms by which maximal symmetry is broken in the context of these models, shedding light on the origin of different Bianchi types. When all relevant length-scales are super-horizon, the simplest Bianchi I models emerge (in which anisotropic quantities appear parallel transported). Finally we highlight the existence of arbitrarily long near-isotropic epochs in models of general Bianchi type (including those without an exact isotropic limit).

The observed angular correlation function of the cosmic microwave background has previously been reported to be anomalous, particularly when measured in regions of the sky uncontaminated by Galactic emission. Recent work by Efstathiou et al. presents a Bayesian comparison of isotropic theories, casting doubt on the significance of the purported anomaly. We extend this analysis to all anisotropic Gaussian theories with vanishing mean (<delta T> = 0), using the much wider class of models to confirm that the anomaly is not likely to point to new physics. On the other hand if there is any new physics to be gleaned, it results from low-l alignments which will be better quantified by a full-sky statistic. We also consider quadratic maximum likelihood power spectrum estimators that are constructed assuming isotropy. The underlying assumptions are therefore false if the ensemble is anisotropic. Nonetheless we demonstrate that, for theories compatible with the observed sky, these estimators (while no longer optimal) remain statistically superior to pseudo-C_l power spectrum estimators.

In recent work, we have shown that it is possible to link quantitatively many aspects of damped Lyman alpha (DLA) absorbers in the spectra of quasars to high resolution simulations of galaxy formation. Using runs from the same series of hydrodynamic numerical studies, we consider the expected properties of Lyman alpha absorbers seen in the spectra of high redshift (z>2) gamma ray burst afterglows (GRB-DLAs). If GRBs are associated with the death of massive stars, their afterglows provide insights into otherwise unprobed regions of protogalactic objects, but detailed physical interpretations are currently embryonic. We find that median impact parameters (measured from the potential minimum) are approximately 1 kpc for GRBs compared with 4 kpc for QSO-DLAs. However, an equally important difference is that GRB-DLAs are predominantly associated with halos of mass 10^10<M_vir/M_sol<10^12, an order of magnitude larger than the hosts of QSO-DLAs. Accordingly, there are differences in the stellar properties of hosts. Our simulations accurately predict the form of the GRB-DLA HI column density distribution, producing quantitative agreement for N_HI>10^19 cm^-2, but they somewhat underpredict the incidence of low column densities N_HI<10^19 cm^-2. Line-of-sight neutral gas metallicities predicted by our simulations (10^-2 < Z/Z_sol < 1) are consistent with the modest observational constraints. Because of large internal dispersions in gas metallicities, this agreement is not significantly compromised by imposing a cut-off on the metallicity of stars able to launch GRBs (Z_star<Z_sol/3), confounding claims that the observed metallicity of GRB-DLAs poses a challenge to current GRB models. (Abridged.)

Combining a recent derivation of the CMB evolution equations for homogeneous but anisotropic (Bianchi) cosmologies with an account of the full linearized dynamical freedoms available in such models, I calculate and discuss the various temperature and polarisation anisotropy patterns which may be formed. Certain anisotropies can be hidden in superhorizon modes at early times, thus avoiding any constraints from nucleosynthesis while nevertheless producing non-trivial redshift-zero temperature patterns in flat and open universes. The results are likely to be more of pedagogical than observational interest, but future work will assess whether such patterns can be matched to anomalies in WMAP results.

If damped Lyman alpha systems (DLAs) contain even modest amounts of dust, the ultraviolet luminosity of the background quasar can be severely diminished. When the spectrum is redshifted, this leads to a bias in optical surveys for DLAs. Previous estimates of the magnitude of this effect are in some tension; in particular, the distribution of DLAs in the column-density -- metallicity plane has led to claims that we may be missing a considerable fraction of metal rich, high column density DLAs, whereas radio surveys do not unveil a substantial population of otherwise hidden systems. Motivated by this tension, we perform a Bayesian parameter estimation analysis of a simple dust obscuration model. We include radio and optical observations of DLAs in our overall likelihood analysis and show that these do not, in fact, constitute conflicting constraints. Our model gives statistical limits on the biasing effects of dust, predicting that only 7% of DLAs are missing from optical samples due to dust obscuration; at 2 sigma confidence, this figure takes a maximum value of 17%. This contrasts with recent claims that DLA incidence rates are underestimated by 30 -- 50%. Optical measures of the mean metallicities of DLAs are found to underestimate the true value by just 0.1 dex (or at most 0.4 dex, 2 sigma confidence limit), in agreement with the radio survey results of Akerman et al. The quantity most affected by dust biasing is the total cosmic density of metals in DLAs, Omega_{Z,DLA}, which is underestimated in optical surveys by a factor of approximately two. (Abridged.)

We investigate the population of z=3 damped Lyman alpha absorbers (DLAs) in a recent series of high resolution galaxy formation simulations. The simulations are of interest because they form at z=0 some of the most realistic disk galaxies to date. No free parameters are available in our study: these have been fixed by physical and z=0 observational constraints, and thus our study provides a genuine consistency test. The precise role of DLAs in galaxy formation remains in debate, but they provide a number of strong constraints on the nature of our simulated bound systems at z=3 because of their coupled information on neutral H I densities, kinematics, metallicity and estimates of star formation activity. Our results, without any parameter-tuning, closely match the observed incidence rate and column density distributions of DLAs. Our simulations are the first to reproduce the distribution of metallicities (with a median of Z_{DLA} = Z_{solar}/20) without invoking observationally unsupported mechanisms such as dust biasing. This is especially encouraging given that these simulations have previously been shown to have a realistic 0<z<2 stellar mass-metallicity relation. Additionally, we see a strong positive correlation between sightline metallicity and low-ion velocity width, the normalization and slope of which comes close to matching recent observational results. However, we somewhat underestimate the number of observed high velocity width systems; the severity of this disagreement is comparable to other recent DLA focused studies. By z=0 the majority of the z=3 neutral gas forming the DLAs has been converted into stars, in agreement with rough physical expectations. [Abridged]

We derive the CMB radiative transfer equation in the form of a multipole hierarchy in the nearly-Friedmann-Robertson-Walker limit of homogeneous, but anisotropic, universes classified via their Bianchi type. Compared with previous calculations, this allows a more sophisticated treatment of recombination, produces predictions for the polarization of the radiation, and allows for reionization. Our derivation is independent of any assumptions about the dynamical behaviour of the field equations, except that it requires anisotropies to be small back to recombination; this is already demanded by observations. We calculate the polarization signal in the Bianchi VIIh case, with the parameters recently advocated to mimic the several large-angle anomalous features observed in the CMB. We find that the peak polarization signal is ~ 1.2 micro K for the best-fit model to the temperature anisotropies, and is mostly confined to multipoles l<10. Remarkably, the predicted large-angle EE and TE power spectra in the Bianchi model are consistent with WMAP observations that are usually interpreted as evidence of early reionization. However, the power in B-mode polarisation is predicted to be similar to the E-mode power and parity-violating correlations are also predicted by the model; the WMAP non-detection of either of these signals casts further strong doubts on the veracity of attempts to explain the large-angle anomalies with global anisotropy. On the other hand, given that there exist further dynamical degrees of freedom in the VIIh universes that are yet to be compared with CMB observations, we cannot at this time definitively reject the anisotropy explanation.

Recent observations suggest the incidence of strong intervening MgII absorption systems along the line-of-sight to gamma ray burst (GRB) afterglows is significantly higher than expected from analogous quasar sightlines. One possible explanation is a geometric effect, arising because MgII absorbers only partially cover the quasar continuum regions, in which case MgII absorbers must be considerably smaller than previous estimates. We investigate the production of abnormal absorption profiles by partial coverage and conclude that the lack of any known anomalous profiles in observed systems, whilst constraining, cannot on its own rule out patchy MgII absorbers. In a separate test, we look for differences in the distribution function of MgII equivalent widths over quasar continuum regions and CIII] emission lines. We show that these anomalies should be observable in any scenario where MgII absorbers are very small, but they are not present in the data. We conclude that models invoking small MgII cloudlets to explain the excess of absorbers seen towards GRBs are ruled out.

Several recent studies suggest a correlation between dark matter halo mass and the shape of the density profile. We re-analyze simulations from Ricotti (2003) in which such a correlation was proposed. We use a standard analysis of the halo density profiles and compare the old simulations to new ones performed with Gadget2, including higher resolution runs. We confirm Ricotti's result that, at virialization, the central log slopes alpha, at 5%-10% of the virial radius are correlated with the halo mass and that the halo concentration is a universal constant. Our results do not contradict the majority of published papers: when using a split power law to fit the density profiles, due to the alpha-concentration degeneracy, the fits are consistent with halos having a universal shape with alpha=1 or 1.5 and concentrations that depend on the mass, in agreement with results published elsewhere. Recently, several groups have found no evidence for convergence of the inner halo profile to a constant power law. The choice of a split power law parameterization used in this letter is motivated by the need to compare our results to previous ones and is formally valid because we are not able to resolve regions where the slope of the fitting function reaches its asymptotic constant value. Using a non-parameterized technique, we also show that the density profiles of dwarf galaxies at z ~ 10 have a log slope shallower than 0.5 within 5% of the virial radius.

Using analytic calculations and N-body simulations we show that in constant density (harmonic) cores, sinking satellites undergo an initial phase of very rapid (super-Chandrasekhar) dynamical friction, after which they experience no dynamical friction at all. For density profiles with a central power law profile of log-slope, $-\alpha$, the infalling satellite heats the background and causes $\alpha$ to decrease. For $\alpha < 0.5$ initially, the satellite generates a small central constant density core and stalls as in the $\alpha = 0$ case. We discuss some astrophysical applications of our results to decaying satellite orbits, galactic bars and mergers of supermassive black hole binaries. In a companion paper we show that a central constant density core can provide a natural solution to the timing problem for Fornax's globular clusters.

Using analytic arguments and a suite of very high resolution (10^3 Msun per particle) cosmological hydro-dynamical simulations, we argue that high redshift, z ~ 10, M ~ 10^8 Msun halos, form the smallest `baryonic building block' (BBB) for galaxy formation. These halos are just massive enough to efficiently form stars through atomic line cooling and to hold onto their gas in the presence of supernovae winds and reionisation. These combined effects, in particular that of the supernovae feedback, create a sharp transition: over the mass range 3-10x10^7 Msun, the BBBs drop two orders ofmagnitude in stellar mass. Below ~2x10^7 Msun, galaxies will be dark with almost no stars and no gas. Above this scale is the smallest unit of galaxy formation: the BBB. A small fraction (~100) of these gas rich BBBs fall in to a galaxy the size of the Milky Way. Ten percent of these survive to become the observed LG dwarf galaxies at the present epoch. Those in-falling halos on benign orbits which keep them far away from the Milky Way or Andromeda manage to retain their gas and slowly form stars - these become the smallest dwarf irregular galax ies; those on more severe orbits lose their gas faster than they can form stars and become the dwarf spheroidals. The remaining 90% of the BBBs will be accreted. We show that this gives a metallicity and total stellar mass consistent with the Milky Way old stellar halo (abridged).