Source author record

Stuart McAlpine

Stuart McAlpine appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

11works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2022arXiv

SIBELIUS-DARK: a galaxy catalogue of the Local Volume from a constrained realisation simulation

We present SIBELIUS-DARK, a constrained realisation simulation of the local volume to a distance of 200~Mpc from the Milky Way. SIBELIUS-DARK is the first study of the \textit{Simulations Beyond The Local Universe} (SIBELIUS) project, which has the goal of embedding a model Local Group-like system within the correct cosmic environment. The simulation is dark-matter-only, with the galaxy population calculated using the semi-analytic model of galaxy formation, GALFORM. We demonstrate that the large-scale structure that emerges from the SIBELIUS constrained initial conditions matches well the observational data. The inferred galaxy population of SIBELIUS-DARK also match well the observational data, both statistically for the whole volume and on an object-by-object basis for the most massive clusters. For example, the $K$-band number counts across the whole sky, and when divided between the northern and southern Galactic hemispheres, are well reproduced by SIBELIUS-DARK. We find that the local volume is somewhat unusual in the wider context of $Λ$CDM: it contains an abnormally high number of supermassive clusters, as well as an overall large-scale underdensity at the level of $\approx 5$\% relative to the cosmic mean. However, whilst rare, the extent of these peculiarities does not significantly challenge the $Λ$CDM model. SIBELIUS-DARK is the most comprehensive constrained realisation simulation of the local volume to date, and with this paper we publicly release the halo and galaxy catalogues at $z=0$, which we hope will be useful to the wider astronomy community.

preprint2022arXiv

The Milky Way's plane of satellites: consistent with $Λ$CDM

The "plane of satellites problem" describes the arrangement of the Milky Way's 11 brightest satellite galaxies in a remarkably thin plane, possibly supported by rotation. This is in apparent contradiction to the standard cosmological model, wherein the Galaxy is surrounded by a dispersion-supported dark matter halo. Here, we show that the reported exceptional anisotropy of the satellite system is strongly contingent on a lopsided radial distribution, which earlier simulations have failed to reproduce, combined with the close but fleeting conjunction of the two most distant satellites, Leo I and Leo II. Using Gaia proper motions, we show that the orbital pole alignment is much more common than previously reported, and reveal the plane of satellites to be transient rather than rotationally supported. Comparing to new simulations, where such short-lived planes are common, we find the Milky Way satellites to be compatible with standard model expectations.

preprint2020arXiv

Galaxy mergers in EAGLE do not induce a significant amount of black hole growth yet do increase the rate of luminous AGN

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 fraction 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

preprint2020arXiv

Interacting galaxies in the IllustrisTNG simulations -- I: Triggered star formation in a cosmological context

We use the IllustrisTNG cosmological hydrodynamical simulations to investigate how the specific star formation rates (sSFRs) of massive galaxies $(M_* > 10^{10} M_\odot)$ depend on the distance to their closest companions. We estimate sSFR enhancements by comparing with control samples that are matched in redshift, stellar mass, local density and isolation, and we restrict our analysis to pairs with stellar mass ratios of 0.1 to 10. At small separations (~15 kpc), the mean sSFR is enhanced by a factor of 2.0 $\pm$ 0.1 in the flagship (110.7 Mpc)$^3$ simulation (TNG100-1). Statistically significant enhancements extend out to 3D separations of 280 kpc in the (302.6 Mpc)$^3$ simulation (TNG300-1). We find similar trends in the EAGLE and Illustris simulations, although their sSFR enhancements are lower than those in TNG100-1 by about a factor of two. Enhancements in IllustrisTNG galaxies are seen throughout the redshift range explored ($0 \leq z < 1$), with the strength of the enhancements decreasing with increasing redshift for galaxies with close companions. In order to more closely compare with observational results, we separately consider 2D projected distances between galaxies in IllustrisTNG. We detect significant sSFR enhancements out to projected separations of 260 kpc in TNG300-1, with projection effects diluting the size of the enhancements by about 20 per cent below 50 kpc. We find similar sSFR enhancements in TNG100-1 and Sloan Digital Sky Survey galaxies, with enhancements extending out to projected separations of about 150 kpc for star-forming galaxies at $z < 0.2$. Finally, by summing over all separations, we estimate that the presence of closest companions boosts the average sSFR of massive galaxies in TNG100-1 by 14.5 per cent.

preprint2020arXiv

Setting the Stage: Structures from Gaussian Random Fields

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

preprint2020arXiv

The Star-Formation Properties of the Observed and Simulated AGN Universe: BAT vs EAGLE

In this paper we present data from 72 low redshift, hard X-ray selected AGN taken from the {\it Swift}-BAT 58 month catalogue. We utilise spectral energy distribution fitting to the optical to IR photometry in order to estimate host galaxy properties. We compare this observational sample to a volume and flux matched sample of AGN from the EAGLE hydrodynamical simulations in order to verify how accurately the simulations can reproduce observed AGN host galaxy properties. After correcting for the known +0.2 dex offset in the SFRs between EAGLE and previous observations, we find agreement in the SFR and X-ray luminosity distributions; however we find that the stellar masses in EAGLE are $0.2 - 0.4$ dex greater than the observational sample, which consequently leads to lower sSFRs. We compare these results to our previous study at high redshift, finding agreement in both the observations and simulations, whereby the widths of sSFR distributions are similar ($\sim0.4-0.6$ dex) and the median of the SFR distributions lie below the star forming main sequence by $\sim0.3-0.5$ dex across all samples. We also use EAGLE to select a sample of AGN host galaxies at high and low redshift and follow their characteristic evolution from $z=8$ to $z=0$. We find similar behaviour between these two samples, whereby star formation is quenched when the black hole goes through its phase of most rapid growth. Utilising EAGLE we find that 23\% of AGN selected at $z\sim0$ are also AGN at high redshift, and that their host galaxies are among the most massive objects in the simulation. Overall we find EAGLE reproduces the observations well, with some minor inconsistencies ($\sim$ 0.2 dex in stellar masses and $\sim$ 0.4 dex in sSFRs).

preprint2016arXiv

A chronicle of galaxy mass assembly in the EAGLE simulation

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

preprint2016arXiv

Music from the heavens - Gravitational waves from supermassive black hole mergers in the EAGLE simulations

We estimate the expected event rate of gravitational wave signals from mergers of supermassive black holes that could be resolved by a space-based interferometer, such as the Evolved Laser Interferometer Space Antenna (eLISA), utilising the reference cosmological hydrodynamical simulation from the EAGLE suite. These simulations assume a $Λ$CDM cosmogony with state-of-the-art subgrid models for radiative cooling, star formation, stellar mass loss, and feedback from stars and accreting black holes. They have been shown to reproduce the observed galaxy population with unprecedented fidelity. We combine the merger rates of supermassive black holes in EAGLE with the latest phenomenological waveform models to calculate the gravitational waves signals from the intrinsic parameters of the merging black holes. The EAGLE models predict $\sim2$ detections per year by a gravitational wave detector such as eLISA. We find that these signals are largely dominated by mergers between seed mass black holes merging at redshifts between $z\sim2$ and $z\sim1$. In order to investigate the dependence on the assumed black hole seed mass, we introduce an additional model with a black hole seed mass an order of magnitude smaller than in our reference model. We also consider a variation of the reference model where a prescription for the expected delays in the black hole merger timescale has been included after their host galaxies merge. We find that the merger rate is similar in all models, but that the initial black hole seed mass could be distinguished through their detected gravitational waveforms. Hence, the characteristic gravitational wave signals detected by eLISA will provide profound insight into the origin of supermassive black holes and the initial mass distribution of black hole seeds.

preprint2016arXiv

Supermassive black holes in the EAGLE Universe. Revealing the observables of their growth

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

preprint2016arXiv

The dark nemesis of galaxy formation: why hot haloes trigger black hole growth and bring star formation to an end

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

preprint2016arXiv

The EAGLE simulations of galaxy formation: public release of halo and galaxy catalogues

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)