We study halo mass functions with the Press-Schechter formalism for interacting dark matter models, where matter power spectra are damped due to dark acoustic oscillations in the early universe. After adopting a smooth window function, we calibrate the analytical model with numerical simulations from the "effective theory of structure formation" (ETHOS) project and fix the model parameters in the high mass regime, $M_{\rm h}\gtrsim3\times10^{10}\;{\rm M}_{\odot}$. We also perform high-resolution cosmological simulations with halo masses down to $M_{\rm h}\sim10^8\;{\rm M}_{\odot}$ to cover a wide mass range for comparison. Although the model is calibrated with ETHOS1 and CDM simulations for high halo masses at redshift $z=0$, it successfully reproduces simulations for two other ETHOS models in the low mass regime at low and high redshifts. As an application, we compare the cumulative number density of haloes to that of observed galaxies at $z=6$, and find the interacting dark matter models with a kinetic decoupling temperature below $0.5\ \rm{keV}$ is disfavored. We also perform the abundance-matching analysis and derive the stellar-halo mass relation for these models at $z=4$. Suppression in halo abundance leads to less massive haloes that host observed galaxies in the stellar mass range $M_*\simeq 10^5-10^7\ {\rm M}_{\odot}$.

We explore a scenario where metal poor globular clusters (GCs) form at the centres of their own dark matter halos in the early universe before reionization. This hypothesis leads to predictions about the abundance, distribution and kinematics of GCs today that we explore using cosmological N-body simulations and analytical modelling. We find that selecting the massive tail of collapsed objects at $z\gtrsim 9$ as GC formation sites leads to four main predictions: i) a highly clustered population of GCs around galaxies today, ii) a natural scaling between number of GCs and halo virial mass that follows roughly the observed trend, iii) a very low number of free floating GCs outside massive halos and iv) GCs should be embedded within massive and extended dark matter (sub)halos. We find that the strongest constraint to the model is given by the combination of (i) and (ii): a mass cut to tagged GCs halos which accounts for the number density of metal poor GCs today predicts a radial distribution that is too extended compared to recent observations. On the other hand, a mass cut sufficient to match the observed half number radius could only explain 60% of the metal poor population. In all cases, observations favour early redshifts for GC formation ($z\geq 15$).

I describe a fast algorithm for the identification of connected sets of points where the point-wise connections are determined by a fixed spatial distance - a task commonly referred to in the cosmological simulation community as Friends-of-Friends (FOF) group finding. This technique sorts particles into fine cells sufficiently compact to guarantee their cohabitants are linked, and uses locality sensitive hashing to search for neighbouring (blocks of) cells. Tests on N-body simulations of up to a billion particles exhibit speed increases of factors up to 20x compared with FOF via trees (a factor around 8 is typical), and is consistently complete in less than the time of a k-d tree construction, giving it an intrinsic advantage over tree-based methods. The code is open-source and available online at https://github.com/pec27/hfof .

We present an update to the chemical enrichment component of the smoothed-particle hydrodynamics model for galaxy formation presented in Scannapieco et al. (2005) in order to address the needs of modelling galactic chemical evolution in realistic cosmological environments. Attribution of the galaxy-scale abundance patterns to individual enrichment mechanisms such as the winds from asymptotic giant branch (AGB) stars or the presence of a prompt fraction of Type Ia supernovae is complicated by the interaction between them and gas cooling, subsequent star formation and gas ejection. In this work we address the resulting degeneracies by extending our implementation to a suite of mechanisms that encompasses different IMFs, models for yields from the aforementioned stars, models for the prompt component of the delay-time-distribution (DTDs) for Type Ia SNe and metallicity-dependent gas cooling rates, and then applying these to both isolated initial conditions and cosmological hydrodynamical zoom simulations. We find DTDs with a large prompt fraction (such as the bimodal and power-law models) have, at z=0, similar abundance patterns compared to the low-prompt component time distributions (uniform or wide Gaussian models). However, some differences appear, such as the former having systematically higher [X/Fe] ratios and narrower [O/Fe] distributions compared to the latter, and a distinct evolution of the [Fe/H] abundance.

We employ isolated N-body simulations to study the response of self-interacting dark matter (SIDM) halos in the presence of the baryonic potentials. Dark matter self-interactions lead to kinematic thermalization in the inner halo, resulting in a tight correlation between the dark matter and baryon distributions. A deep baryonic potential shortens the phase of SIDM core expansion and triggers core contraction. This effect can be further enhanced by a large self-scattering cross section. We find the final SIDM density profile is sensitive to the baryonic concentration and the strength of dark matter self-interactions. Assuming a spherical initial halo, we also study evolution of the SIDM halo shape together with the density profile. The halo shape at later epochs deviates from spherical symmetry due to the influence of the non-spherical disc potential, and its significance depends on the baryonic contribution to the total gravitational potential, relative to the dark matter one. In addition, we construct a multi-component model for the Milky Way, including an SIDM halo, a stellar disc and a bulge, and show it is consistent with observations from stellar kinematics and streams.

The efficient simulation of isotropic Gaussian random fields on the unit sphere is a task encountered frequently in numerical applications. A fast algorithm based on Markov properties and fast Fourier Transforms in 1d is presented that generates samples on an n x n grid in O(n^2 log n). Furthermore, an efficient method to set up the necessary conditional covariance matrices is derived and simulations demonstrate the performance of the algorithm. An open source implementation of the code has been made available at https://github.com/pec27/smerfs .

The existence of stellar halos around galaxies is a natural prediction of the hierarchical nature of the LambdaCDM model. Recent observations of Milky Way-like galaxies have revealed a wide range in stellar halo mass, including cases with no statistically significant detection of a stellar halo, as in the case of M101, NGC3351 and NGC1042. We use the Illustris simulation to investigate the scatter in stellar halo content and, in particular, to study the formation of galaxies with the smallest fraction of this diffuse component. Stellar halos are far from spherical, which diminishes the surface brightness of the stellar halo for face-on disks. Once accounting for projection effects, we find that the stellar halo fraction f_SH correlates strongly with galaxy morphology and star formation rate, but not with environment, in agreement with observations. Galaxies with the lowest stellar halo fractions are disk-dominated, star-forming and assemble their dark matter halos on average earlier than galaxies with similar stellar masses. Accreted satellites are also lower in stellar mass and have earlier infall times than centrals with high f_SH. In situ rather than accreted stars dominate the stellar halos of galaxies with the lowest stellar halo fractions, with a transition radius from in situ to accretion-dominated r~45 kpc. Our results extrapolated to real galaxies such as M101 may indicate that these galaxies inhabit old halos which endured mergers only at higher redshifts and evolved relatively unperturbed in the last ~10 Gyrs.

Galactic rotation curves are a fundamental constraint for any cosmological model. We use controlled N-body simulations of galaxies to study the gravitational effect of baryons in a scenario with collisionless cold dark matter (CDM) versus one with a self-interacting dark matter (SIDM) component. In particular, we examine the inner profiles of the rotation curves in the velocity range Vmax = 30-250 km/s, whose diversity has been found to be greater than predicted by the Lambda-CDM scenario. We find that the scatter in the observed rotation curves exceeds that predicted by dark matter only mass-concentration relations in either the CDM nor SIDM models. Allowing for realistic baryonic content and spatial distributions, however, helps create a large variety of rotation curve shapes; which is in better agreement with observations in the case of self-interactions due to the characteristic cored profiles being more accommodating to the slowly rising rotation curves than CDM. We find individual fits to model two of the most remarkable outliers of similar Vmax, UGC 5721 and IC 2574, the former a cusp-like rotation curve and the latter a seemingly 8 kpc cored profile. This diversity in SIDM arises as permutations of overly concentrated halos with compact baryonic distributions versus underdense halos with extended baryonic disks. The SIDM solution is promising and its feasibility ultimately depends on the sampling of the halo mass-concentration relation and its interplay with the baryonic profiles, emphasising the need for a better understanding of the frequency of extreme outliers present in current observational samples.

Intermediate surface brightness (ISB) galaxies are less numerous than their counterparts at high and low surface brightness (HSB and LSB). Investigating ISB characteristics from a sample from the S4G survey, complete down to M_B=-16, we find that they have intermediate stellar, gas and baryonic masses and on average as much gas as stars. They lie on the (baryonic) Tully-Fisher relation between HSBs and LSBs, although they present a higher scatter than the latter. Their stellar to baryonic mass ratios have intermediate values unlike their condensed baryonic fractions. By comparing their environments, as classified by the eigenvalues of the velocity shear tensor of local constrained simulations, ISBs have a 5-10% probability higher (smaller) to be in sheets (filaments) with respect to HSBs and LSBs. Additionally, for galaxies in filaments (with close neighbors), the mass and mu_0 are correlated at 2.5 (2) sigma more than for those in sheets. ISBs live in regions where the divergence of the velocity field is smaller than where HSBs and LSBs live, a result at more than 50% significance. ISBs may exist as an unstable transition state between LSBs and HSBs, the low flow activity environment maximally encouraging their formation. Interaction events altering the central baryon fraction could happen at a lower rate in these less dense environment, whilst in the higher density environments the LSBs are primarily satellite galaxies, whose accretion is sufficiently constrained that it fails to promote them to HSBs.

We study the two main constituent galaxies of a constrained simulation of the Local Group as candidates for the Milky Way (MW) and Andromeda (M31). We focus on the formation of the stellar discs and its relation to the formation of the group as a rich system with two massive galaxies, and investigate the effects of mergers and accretion as drivers of morphological transformations. We use a state-of-the-art hydrodynamical code which includes star formation, feedback and chemical enrichment to carry out our study. We run two simulations, where we include or neglect the effects of radiation pressure from stars, to investigate the impact of this process on the morphologies and star formation rates of the simulated galaxies. We find that the simulated M31 and MW have different formation histories, even though both inhabit, at z=0, the same environment. These differences directly translate into and explain variations in their star formation rates, in-situ fractions and final morphologies. The M31 candidate has an active merger history, as a result of which its stellar disc is unable to survive unaffected until the present time. In contrast, the MW candidate has a smoother history with no major mergers at late times, and forms a disc that grows steadily; at z=0 the simulated MW has an extended, rotationally-supported disc which is dominant over the bulge. Our two feedback implementations predict similar evolution of the galaxies and their discs, although some variations are detected, the most important of which is the formation time of the discs: in the model with weaker/stronger feedback the discs form earlier/later. In summary, by comparing the formation histories of the two galaxies, we conclude that the particular merger/accretion history of a galaxy rather than its environment at the LG-scales is the main driver of the formation and subsequent growth or destruction of galaxy discs.

In this letter we present, for the first time, a study of star formation rate, gas fraction and galaxy morphology of a constrained simulation of the Milky Way (MW) and Andromeda (M31) galaxies, compared to other MW-mass galaxies. By combining with unconstrained simulations we cover a sufficient volume to compare these galaxies environmental densities ranging from the field to that of the Local Group (LG). This is particularly relevant as it has been shown that, quite generally, galaxy properties depend intimately upon their environment, most prominently when galaxies in clusters are compared to those in the field. For galaxies in loose groups such as the LG, however, environmental effects have been less clear. We consider the galaxy's environmental density in spheres of 1200 kpc (comoving) and find that whilst environment does not appear to directly affect morphology, there is a positive trend with star formation rates. This enhancement in star formation occurs systematically for galaxies in higher density environments, regardless whether they are part of the LG or in filaments. Our simulations suggest that the richer environment at Mpc-scales may help replenish the star-forming gas, allowing higher specific star formation rates in galaxies such as the MW.

The abundance evolution of galaxies depends critically on the balance between the mixing of metals in their interstellar medium, the inflow of new gas and the outflow of enriched gas. We study these processes in gas columns perpendicular to a galactic disk using sub-parsec resolution simulations that track stellar ejecta with the Flash code. We model a simplified interstellar medium stirred and enriched by supernovae and their progenitors. We vary the density distribution of the gas column and integrate our results over an exponential disk to predict wind and ISM enrichment properties for disk galaxies. We find that winds from more massive galaxies are hotter and more highly enriched, in stark contrast to that which is often assumed in galaxy formation models. We use these findings in a simple model of galactic enrichment evolution, in which the metallicity of forming galaxies is the result of accretion of nearly pristine gas and outflow of enriched gas along an equilibrium sequence. We compare these predictions to the observed mass-metallicity relation, and demonstrate how the galaxy's gas fraction is a key controlling parameter. This explains the observed flattening of the mass-metallicity relation at higher stellar masses.

Feedback from supernovae is an essential aspect of galaxy formation. In order to improve subgrid models of feedback we perform a series of numerical experiments to investigate how supernova explosions power galactic winds. We use the Flash hydrodynamic code to model a simplified ISM, including gravity, hydrodynamics, radiative cooling above 10,000 K, and star formation that reproduces the Kennicutt-Schmidt relation. By simulating a small patch of the ISM in a tall box perpendicular to the disk, we obtain sub-parsec resolution allowing us to resolve individual supernova events and we investigate how the wind properties depend on those of the ISM and the galaxy. We find that outflows are more efficient in disks with lower surface densities or gas fractions. A simple model in which the warm cloudy medium is the barrier that limits the expansion of blast waves reproduces the scaling of outflow properties with disk parameters at high star formation rates. The scaling we find sets the investigation of galaxy winds on a new footing, providing a physically motivated sub-grid description of winds that can be implemented in cosmological hydrodynamic simulations and phenomenological models. [Abridged]

We present a study of cooling in radiative shocks simulated with smoothed particle hydrodynamics (SPH) and adaptive mesh refinement codes. We obtain a similarity solution for a shock-tube problem in the presence of radiative cooling, and test how well the solution is reproduced in Gadget and Flash. Shock broadening governed by the details of the numerical scheme (artificial viscosity or Riemann solvers) leads to potentially significant overcooling in both codes. We interpret our findings in terms of a resolution criterion, and apply it to realistic simulations of cosmological accretion shocks onto galaxy haloes, cold accretion and thermal feedback from supernovae or active galactic nuclei. To avoid numerical overcooling of accretion shocks onto haloes that should develop a hot corona requires a particle or cell mass resolution of 10^6 M_sun, which is within reach of current state-of-the-art simulations. At this mass resolution, thermal feedback in the interstellar medium of a galaxy requires temperatures of supernova or AGN driven bubbles to be in excess of 10^7 K at densities of n_H=1.0 cm^-3, in order to avoid spurious suppression of the feedback by numerical overcooling.