We introduce Simulations Beyond The Local Universe (SIBELIUS) that connect the Local Group to its cosmic environment. We show that introducing hierarchical small-scale perturbations to a density field constrained on large scales by observations provi
des an efficient way to explore the sample space of Local Group analogues. From more than 60 000 simulations, we identify a hierarchy of Local Group characteristics emanating from different scales: the total mass, orientation, orbital energy and the angular momentum are largely determined by modes above $lambda$ = 1.6 comoving Mpc (cMpc) in the primordial density field. Smaller scale variations are mostly manifest as perturbations to the MW-M31 orbit, and we find that the observables commonly used to describe the Local Group -- the MW-M31 separation and radial velocity -- are transient and depend on specifying scales down to 0.2 cMpc in the primordial density field. We further find that the presence of M33/LMC analogues significantly affects the MW-M31 orbit and its sensitivity to small-scale perturbations. We construct initial conditions that lead to the formation of a Local Group whose primary observables precisely match the current observations.
We investigate the connection between galaxy--galaxy mergers and enhanced black hole (BH) growth using the cosmological hydrodynamical EAGLE simulation. We do this via three methods of analysis, investigating: the merger fraction of AGN, the AGN frac
tion of merging systems and the AGN fraction of galaxies with close companions. In each case, we find an increased abundance of AGN within merging systems relative to control samples of inactive or isolated galaxies (by up to a factor of $approx 3$ depending on the analysis method used), confirming that mergers are enhancing BH accretion rates for at least a subset of the galaxy population. The greatest excess of AGN triggered via a merger are found in lower mass ($M_* sim 10^{10}$~Msol) gas rich ($f_{mathrm{gas}} > 0.2$) central galaxies with lower mass BHs ($M_{mathrm{BH}} sim 10^{7}$~Msol) at lower redshifts ($z<1$). We find no enhancement of AGN triggered via mergers in more massive galaxies ($M_* gtrsim 10^{11}$~Msol). The enhancement of AGN is not uniform throughout the phases of a merger, and instead peaks within the early emph{remnants} of merging systems (typically lagging $approx 300$~Myr post-coalescence of the two galaxies at $z=0.5$). We argue that neither major ($M_{mathrm{*,1}} / M_{mathrm{*,2}} geq frac{1}{4}$) nor minor mergers ($frac{1}{10} < M_{mathrm{*,1}} / M_{mathrm{*,2}} < frac{1}{4}$) are statistically relevant for enhancing BH masses globally. Whilst at all redshifts the galaxies experiencing a merger have accretion rates that are on average 2--3 times that of isolated galaxies, the majority of mass that is accreted onto BHs occurs outside the periods of a merger. We compute that on average no more than 15% of a BHs final day mass comes from the enhanced accretion rates triggered via a merger
We exploit EAGLE, a cosmological hydrodynamical simulation, to reproduce the selection of the observed sub-millimeter (submm) galaxy population by selecting the model galaxies at $z geq 1$ with mock submm fluxes $S_{850} geq 1$ mJy. There is a reason
able agreement between the galaxies within this sample and the properties of the observed submm population, such as their star formation rates (SFRs) at $z<3$, redshift distribution and many integrated galaxy properties. We find that the bulk of the $S_{850} geq 1$ mJy model population is at $z = 2.5$, and that they are massive galaxies ($M_* sim 10^{11}$ Msol) with high dust masses ($M_{mathrm{dust}} sim 10^{8}$ Msol), gas fractions ($f_{mathrm{gas}} approx 50$%) and SFRs ($dot M_* approx 100$ Msol/yr). They have major and minor merger fractions similar to the general population, suggesting that mergers are not the primary driver of the model submm galaxies. Instead, the $S_{850} geq 1$ mJy model galaxies yield high SFRs primarily because they maintain a significant gas reservoir as a result of hosting an undermassive black hole. In addition, we find that not all highly star-forming EAGLE galaxies have submm fluxes $S_{850} > 1$ mJy. Thus, we investigate the nature of $z geq 1$ highly star-forming Submm-Faint galaxies (i.e., $dot M_* geq 80$ Msol/yr but $S_{850}< 1$ mJy). We find they are similar to the model submm galaxies; being gas rich and hosting undermassive black holes, however they are typically lower mass ($M_* sim 10^{10}$ Msol) and are at higher redshifts ($z>4$). These typically higher-$z$ galaxies show stronger evidence for having been triggered by major mergers, and critically, they are likely missed by current submm surveys due to their higher dust temperatures. This suggests a potentially even larger contribution to the SFR density at $z > 3$ from dust-obscured systems than implied by current observations.
We investigate the rapid growth phase of supermassive black holes (BHs) within the hydrodynamical cosmological eagle simulation. This non-linear phase of BH growth occurs within $sim$$L_{*}$ galaxies, embedded between two regulatory states of the gal
axy host: in sub $L_{*}$ galaxies efficient stellar feedback regulates the gas inflow onto the galaxy and significantly reduces the growth of the central BH, while in galaxies more massive than $L_{*}$ efficient AGN feedback regulates the gas inflow onto the galaxy and curbs further non-linear BH growth. We find evolving critical galaxy and halo mass scales at which rapid BH growth begins. Galaxies in the low-redshift Universe transition into the rapid BH growth phase in haloes that are approximately an order of magnitude more massive than their high-redshift counterparts (M{200} $approx 10^{12.4}$~Msol at $z approx 0$ decreasing to M{200} $approx 10^{11.2}$~Msol at $z approx 6$). Instead, BHs enter the rapid growth phase at a fixed critical halo virial temperature ($T_{mathrm{vir}} approx 10^{5.6}$~K). We additionally show that major galaxy--galaxy interactions ($mu geq frac{1}{4}$, where $mu$ is the stellar mass ratio) play a substantial role in triggering the rapid growth phase of BHs in the low-redshift Universe, whilst potentially having a lower influence at high redshift. Approximately 40% of BHs that initiate the rapid BH growth phase at $z approx 0$ do so within $pm 0.5$ dynamical times of a major galaxy--galaxy merger, a fourfold increase above what is expected from the background merger rate. We find that minor mergers ($frac{1}{10} leq mu < frac{1}{4}$) have a substantially lower influence in triggering the rapid growth phase at all epochs.
We investigate the connection between the star formation rate (SFR) of galaxies and their central black hole accretion rate (BHAR) using the EAGLE cosmological hydrodynamical simulation. We find, in striking concurrence with recent observational stud
ies, that the <SFR>--BHAR relation for an AGN selected sample produces a relatively flat trend, whilst the <BHAR>--SFR relation for a SFR selected sample yields an approximately linear trend. These trends remain consistent with their instantaneous equivalents even when both SFR and BHAR are time-averaged over a period of 100~Myr. There is no universal relationship between the two growth rates. Instead, SFR and BHAR evolve through distinct paths that depend strongly on the mass of the host dark matter halo. The galaxies hosted by haloes of mass M200 $lesssim 10^{11.5}$Msol grow steadily, yet black holes (BHs) in these systems hardly grow, yielding a lack of correlation between SFR and BHAR. As haloes grow through the mass range $10^{11.5} lesssim$ M200 $lesssim 10^{12.5 }$Msol BHs undergo a rapid phase of non-linear growth. These systems yield a highly non-linear correlation between the SFR and BHAR, which are non-causally connected via the mass of the host halo. In massive haloes (M200 $gtrsim 10^{12.5}$Msol) both SFR and BHAR decline on average with a roughly constant scaling of SFR/BHAR $sim 10^{3}$. Given the complexity of the full SFR--BHAR plane built from multiple behaviours, and from the large dynamic range of BHARs, we find the primary driver of the different observed trends in the <SFR>--BHAR and <BHAR>--SFR relationships are due to sampling considerably different regions of this plane.
We present the public data release of halo and galaxy catalogues extracted from the EAGLE suite of cosmological hydrodynamical simulations of galaxy formation. These simulations were performed with an enhanced version of the GADGET code that includes
a modified hydrodynamics solver, time-step limiter and subgrid treatments of baryonic physics, such as stellar mass loss, element-by-element radiative cooling, star formation and feedback from star formation and black hole accretion. The simulation suite includes runs performed in volumes ranging from 25 to 100 comoving megaparsecs per side, with numerical resolution chosen to marginally resolve the Jeans mass of the gas at the star formation threshold. The free parameters of the subgrid models for feedback are calibrated to the redshift z=0 galaxy stellar mass function, galaxy sizes and black hole mass - stellar mass relation. The simulations have been shown to match a wide range of observations for present-day and higher-redshift galaxies. The raw particle data have been used to link galaxies across redshifts by creating merger trees. The indexing of the tree produces a simple way to connect a galaxy at one redshift to its progenitors at higher redshift and to identify its descendants at lower redshift. In this paper we present a relational database which we are making available for general use. A large number of properties of haloes and galaxies and their merger trees are stored in the database, including stellar masses, star formation rates, metallicities, photometric measurements and mock gri images. Complex queries can be created to explore the evolution of more than 10^5 galaxies, examples of which are provided in appendix. (abridged)
