-
SWIFT: A modern highly-parallel gravity and smoothed particle hydrodynamics solver for astrophysical and cosmological applications
Authors:
Matthieu Schaller,
Josh Borrow,
Peter W. Draper,
Mladen Ivkovic,
Stuart McAlpine,
Bert Vandenbroucke,
Yannick Bahé,
Evgenii Chaikin,
Aidan B. G. Chalk,
Tsang Keung Chan,
Camila Correa,
Marcel van Daalen,
Willem Elbers,
Pedro Gonnet,
Loïc Hausammann,
John Helly,
Filip Huško,
Jacob A. Kegerreis,
Folkert S. J. Nobels,
Sylvia Ploeckinger,
Yves Revaz,
William J. Roper,
Sergio Ruiz-Bonilla,
Thomas D. Sandnes,
Yolan Uyttenhove
, et al. (2 additional authors not shown)
Abstract:
Numerical simulations have become one of the key tools used by theorists in all the fields of astrophysics and cosmology. The development of modern tools that target the largest existing computing systems and exploit state-of-the-art numerical methods and algorithms is thus crucial. In this paper, we introduce the fully open-source highly-parallel, versatile, and modular coupled hydrodynamics, gra…
▽ More
Numerical simulations have become one of the key tools used by theorists in all the fields of astrophysics and cosmology. The development of modern tools that target the largest existing computing systems and exploit state-of-the-art numerical methods and algorithms is thus crucial. In this paper, we introduce the fully open-source highly-parallel, versatile, and modular coupled hydrodynamics, gravity, cosmology, and galaxy-formation code SWIFT. The software package exploits hybrid shared- and distributed-memory task-based parallelism, asynchronous communications, and domain-decomposition algorithms based on balancing the workload, rather than the data, to efficiently exploit modern high-performance computing cluster architectures. Gravity is solved for using a fast-multipole-method, optionally coupled to a particle mesh solver in Fourier space to handle periodic volumes. For gas evolution, multiple modern flavours of Smoothed Particle Hydrodynamics are implemented. SWIFT also evolves neutrinos using a state-of-the-art particle-based method. Two complementary networks of sub-grid models for galaxy formation as well as extensions to simulate planetary physics are also released as part of the code. An extensive set of output options, including snapshots, light-cones, power spectra, and a coupling to structure finders are also included. We describe the overall code architecture, summarise the consistency and accuracy tests that were performed, and demonstrate the excellent weak-scaling performance of the code using a representative cosmological hydrodynamical problem with $\approx$$300$ billion particles. The code is released to the community alongside extensive documentation for both users and developers, a large selection of example test problems, and a suite of tools to aid in the analysis of large simulations run with SWIFT.
△ Less
Submitted 29 March, 2024; v1 submitted 22 May, 2023;
originally announced May 2023.
-
Modelling the accretion and feedback of supermassive black hole binaries in gas-rich galaxy mergers
Authors:
Shihong Liao,
Peter H. Johansson,
Matias Mannerkoski,
Dimitrios Irodotou,
Francesco Paolo Rizzuto,
Stuart McAlpine,
Antti Rantala,
Alexander Rawlings,
Till Sawala
Abstract:
We introduce a new model for the accretion and feedback of supermassive black hole (SMBH) binaries to the KETJU code, which enables us to resolve the evolution of SMBH binaries down to separations of tens of Schwarzschild radii in gas-rich galaxy mergers. Our subgrid binary accretion model extends the widely used Bondi--Hoyle--Lyttleton accretion into the binary phase and incorporates preferential…
▽ More
We introduce a new model for the accretion and feedback of supermassive black hole (SMBH) binaries to the KETJU code, which enables us to resolve the evolution of SMBH binaries down to separations of tens of Schwarzschild radii in gas-rich galaxy mergers. Our subgrid binary accretion model extends the widely used Bondi--Hoyle--Lyttleton accretion into the binary phase and incorporates preferential mass accretion onto the secondary SMBH, which is motivated by results from small-scale hydrodynamical circumbinary disc simulations. We perform idealised gas-rich disc galaxy merger simulations using pure thermal or pure kinetic active galactic nuclei (AGN) feedback. Our binary accretion model provides more physically motivated SMBH mass ratios, which are one of the key parameters for computing gravitational wave (GW) induced recoil velocities. The merger time-scales of our simulated SMBH binaries are in the range $t_{\rm merge}{\sim} 10$--$400$ Myr. Prograde in-plane equal-mass galaxy mergers lead to the shortest merger time-scales, as they experience the strongest starbursts, with the ensuing high stellar density resulting in a rapid SMBH coalescence. Compared to the thermal AGN feedback, the kinetic AGN feedback predicts longer merger time-scales and results in more core-like stellar profiles, as it is more effective in removing gas from the galaxy centre and quenching star formation. This suggests that the AGN feedback implementation plays a critical role in modelling SMBH coalescences. Our model will be useful for improving the modelling of SMBH mergers in gas-rich galaxies, the prime targets for the upcoming LISA GW observatory.
△ Less
Submitted 18 February, 2023; v1 submitted 21 November, 2022;
originally announced November 2022.
-
The Milky Way's plane of satellites: consistent with $Λ$CDM
Authors:
Till Sawala,
Marius Cautun,
Carlos S. Frenk,
John Helly,
Jens Jasche,
Adrian Jenkins,
Peter H. Johansson,
Guilhem Lavaux,
Stuart McAlpine,
Matthieu Schaller
Abstract:
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…
▽ More
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.
△ Less
Submitted 5 May, 2022;
originally announced May 2022.
-
SIBELIUS-DARK: a galaxy catalogue of the Local Volume from a constrained realisation simulation
Authors:
Stuart McAlpine,
John C. Helly,
Matthieu Schaller,
Till Sawala,
Guilhem Lavaux,
Jens Jasche,
Carlos S. Frenk,
Adrian Jenkins,
John R. Lucey,
Peter H. Johansson
Abstract:
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…
▽ More
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.
△ Less
Submitted 8 February, 2022;
originally announced February 2022.
-
Supermassive black holes in cosmological simulations II: the AGN population and predictions for upcoming X-ray missions
Authors:
Melanie Habouzit,
Rachel S. Somerville,
Yuan Li,
Shy Genel,
James Aird,
Daniel Anglés-Alcázar,
Romeel Davé,
Iskren Y. Georgiev,
Stuart McAlpine,
Yetli Rosas-Guevara,
Yohan Dubois,
Dylan Nelson,
Eduardo Bañados,
Lars Hernquist,
Sébastien Peirani,
Mark Vogelsberger
Abstract:
In large-scale hydrodynamical cosmological simulations, the fate of massive galaxies is mainly dictated by the modeling of feedback from active galactic nuclei (AGN). The amount of energy released by AGN feedback is proportional to the mass that has been accreted onto the BHs, but the exact sub-grid modeling of AGN feedback differs in all simulations. Whilst modern simulations reliably produce pop…
▽ More
In large-scale hydrodynamical cosmological simulations, the fate of massive galaxies is mainly dictated by the modeling of feedback from active galactic nuclei (AGN). The amount of energy released by AGN feedback is proportional to the mass that has been accreted onto the BHs, but the exact sub-grid modeling of AGN feedback differs in all simulations. Whilst modern simulations reliably produce populations of quiescent massive galaxies at z<2, it is also crucial to assess the similarities and differences of the responsible AGN populations. Here, we compare the AGN population of the Illustris, TNG100, TNG300, Horizon-AGN, EAGLE, and SIMBA simulations. The AGN luminosity function (LF) varies significantly between simulations. Although in agreement with current observational constraints at z=0, at higher redshift the agreement of the LFs deteriorates with most simulations producing too many AGN of L_{x, 2-10 keV}~10^43-10^44 erg/s. AGN feedback in some simulations prevents the existence of any bright AGN with L_{x, 2-10 keV}>=10^45 erg/s (although this is sensitive to AGN variability), and leads to smaller fractions of AGN in massive galaxies than in the observations at z<=2. We find that all the simulations fail at producing a number density of AGN in good agreement with observational constraints for both luminous (L_{x, 2-10 keV}~10^43-10^45 erg/s) and fainter (L_{x, 2-10 keV}~10^42-10^43 erg/s) AGN, and at both low and high redshift. These differences can aid us in improving future BH and galaxy subgrid modeling in simulations. Upcoming X-ray missions (e.g., Athena, AXIS, and LynX) will bring faint AGN to light and new powerful constraints. After accounting for AGN obscuration, we find that the predicted number density of detectable AGN in future surveys spans at least one order of magnitude across the simulations, at any redshift.
△ Less
Submitted 2 November, 2021;
originally announced November 2021.
-
The SIBELIUS Project: E Pluribus Unum
Authors:
Till Sawala,
Stuart McAlpine,
Jens Jasche,
Guilhem Lavaux,
Adrian Jenkins,
Peter H. Johansson,
Carlos S. Frenk
Abstract:
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 provides 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…
▽ More
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 provides 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 $λ$ = 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.
△ Less
Submitted 22 March, 2021;
originally announced March 2021.
-
The Star-Formation Properties of the Observed and Simulated AGN Universe: BAT vs EAGLE
Authors:
Thomas M. Jackson,
D. J. Rosario,
D. M. Alexander,
J. Scholtz,
Stuart McAlpine,
R. G. Bower
Abstract:
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 acc…
▽ More
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).
△ Less
Submitted 12 August, 2020;
originally announced August 2020.
-
Supermassive black holes in cosmological simulations I: M_BH-M_star relation and black hole mass function
Authors:
Melanie Habouzit,
Yuan Li,
Rachel S. Somerville,
Shy Genel,
Annalisa Pillepich,
Marta Volonteri,
Romeel Davé,
Yetli Rosas-Guevara,
Stuart McAlpine,
Sébastien Peirani,
Lars Hernquist,
Daniel Anglés-Alcázar,
Amy Reines,
Richard Bower,
Yohan Dubois,
Dylan Nelson,
Christophe Pichon,
Mark Vogelsberger
Abstract:
The past decade has seen significant progress in understanding galaxy formation and evolution using large-scale cosmological simulations. While these simulations produce galaxies in overall good agreement with observations, they employ different sub-grid models for galaxies and supermassive black holes (BHs). We investigate the impact of the sub-grid models on the BH mass properties of the Illustr…
▽ More
The past decade has seen significant progress in understanding galaxy formation and evolution using large-scale cosmological simulations. While these simulations produce galaxies in overall good agreement with observations, they employ different sub-grid models for galaxies and supermassive black holes (BHs). We investigate the impact of the sub-grid models on the BH mass properties of the Illustris, TNG100, TNG300, Horizon-AGN, EAGLE, and SIMBA simulations, focusing on the M_BH-M_star relation and the BH mass function. All simulations predict tight M_BH-M_star relations, and struggle to produce the lowest (M_BH< 10^7.5 Msun) in galaxies of M_star~10^10.5-10^11.5 Msun. While the time evolution of the mean M_BH-M_star relation is mild (<1 dex in BH mass for 0<z<5) for all the simulations, its linearity (shape) and normalization varies from simulation to simulation. The strength of SN feedback has a large impact on the linearity and time evolution for M_star<10^10.5 Msun. We find that the low-mass end is a good discriminant of the simulation models, and highlights the need for new observational constraints. At the high-mass end, strong AGN feedback can suppress the time evolution of the relation normalization. Compared with the observations of the local universe, we find an excess of BHs with M_BH>10^9 Msun in most of the simulations. The BH mass function is dominated by efficiently accreting BHs (log10 f_Edd >-2$) at high redshifts, and transitions progressively from the high-mass to the low-mass end to be governed by inactive BHs. The transition time and the contribution of active BHs are different among the simulations, and can be used to evaluate models against observations.
△ Less
Submitted 1 November, 2021; v1 submitted 17 June, 2020;
originally announced June 2020.
-
Setting the Stage: Structures from Gaussian Random Fields
Authors:
Till Sawala,
Adrian Jenkins,
Stuart McAlpine,
Jens Jasche,
Guilhem Lavaux,
Peter H. Johansson,
Carlos S. Frenk
Abstract:
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 o…
▽ More
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.
△ Less
Submitted 9 March, 2020;
originally announced March 2020.
-
Interacting galaxies in the IllustrisTNG simulations -- I: Triggered star formation in a cosmological context
Authors:
David R. Patton,
Kieran D. Wilson,
Colin J. Metrow,
Sara L. Ellison,
Paul Torrey,
Westley Brown,
Maan H. Hani,
Stuart McAlpine,
Jorge Moreno,
Joanna Woo
Abstract:
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 wi…
▽ More
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.
△ Less
Submitted 30 March, 2020; v1 submitted 29 February, 2020;
originally announced March 2020.
-
Galaxy mergers in EAGLE do not induce a significant amount of black hole growth yet do increase the rate of luminous AGN
Authors:
Stuart McAlpine,
Chris M. Harrison,
David J. Rosario,
David M. Alexander,
Sara L. Ellison,
Peter H. Johansson,
David R. Patton
Abstract:
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 sys…
▽ More
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
△ Less
Submitted 7 May, 2020; v1 submitted 3 February, 2020;
originally announced February 2020.
-
The diverse evolutionary pathways of post-starburst galaxies
Authors:
M. M. Pawlik,
S. McAlpine,
J. W. Trayford,
V. Wild,
R. Bower,
R. A. Crain,
M. Schaller,
J. Schaye
Abstract:
About 35 years ago a class of galaxies with unusually strong Balmer absorption lines and weak emission lines was discovered in distant galaxy clusters. These objects, alternatively referred to as post-starburst, E+A or k+a galaxies, are now known to occur in all environments and at all redshifts, with many exhibiting compact morphologies and low-surface brightness features indicative of past galax…
▽ More
About 35 years ago a class of galaxies with unusually strong Balmer absorption lines and weak emission lines was discovered in distant galaxy clusters. These objects, alternatively referred to as post-starburst, E+A or k+a galaxies, are now known to occur in all environments and at all redshifts, with many exhibiting compact morphologies and low-surface brightness features indicative of past galaxy mergers. They are commonly thought to represent galaxies that are transitioning from blue to red sequence, making them critical to our understanding of the origins of galaxy bimodality. However, recent observational studies have questioned this simple interpretation. From observations alone, it is challenging to disentangle the different mechanisms that lead to the quenching of star formation in galaxies. Here we present examples of three different evolutionary pathways that lead to galaxies with strong Balmer absorption lines in the EAGLE simulation: classical blue-to-red quenching, blue-to-blue cycle and red-to-red rejuvenation. The first two are found in both post-starburst galaxies and galaxies with truncated star formation. Each pathway is consistent with scenarios hypothesised for observational samples. The fact that 'post-starburst' signatures can be attained via various evolutionary channels explains the diversity of observed properties, and lends support to the idea that slower quenching channels are important at low redshift.
△ Less
Submitted 26 March, 2019;
originally announced March 2019.
-
The nature of sub-millimeter and highly star-forming galaxies in the EAGLE simulation
Authors:
Stuart McAlpine,
Ian Smail,
Richard G. Bower,
Mark A. Swinbank,
James W. Trayford,
Tom Theuns,
Maarten Baes,
Peter Camps,
Robert A. Crain,
Joop Schaye
Abstract:
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 reasonable agreement between the galaxies within this sample and the properties of the observed submm population, such as their star formation rates (SFRs…
▽ More
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 reasonable 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.
△ Less
Submitted 16 January, 2019;
originally announced January 2019.
-
The gas fractions of dark matter haloes hosting simulated $\sim L^\star$ galaxies are governed by the feedback history of their black holes
Authors:
Jonathan J. Davies,
Robert A. Crain,
Ian G. McCarthy,
Benjamin D. Oppenheimer,
Joop Schaye,
Matthieu Schaller,
Stuart McAlpine
Abstract:
We examine the origin of scatter in the relationship between the gas fraction and mass of dark matter haloes hosting present-day $\sim L^\star$ central galaxies in the EAGLE simulations. The scatter is uncorrelated with the accretion rate of the central galaxy's black hole (BH), but correlates strongly and negatively with the BH's mass, implicating differences in the expulsion of gas by active gal…
▽ More
We examine the origin of scatter in the relationship between the gas fraction and mass of dark matter haloes hosting present-day $\sim L^\star$ central galaxies in the EAGLE simulations. The scatter is uncorrelated with the accretion rate of the central galaxy's black hole (BH), but correlates strongly and negatively with the BH's mass, implicating differences in the expulsion of gas by active galactic nucleus feedback, throughout the assembly of the halo, as the main cause of scatter. Haloes whose central galaxies host undermassive BHs also tend to retain a higher gas fraction, and exhibit elevated star formation rates (SFRs). Diversity in the mass of central BHs stems primarily from diversity in the dark matter halo binding energy, as these quantities are strongly and positively correlated at fixed halo mass, such that $\sim L^\star$ galaxies hosted by haloes that are more (less) tightly-bound develop central BHs that are more (less) massive than is typical for their halo mass. Variations in the halo gas fraction at fixed halo mass are reflected in both the soft X-ray luminosity and thermal Sunyaev-Zel'dovich flux, suggesting that the prediction of a strong coupling between the properties of galaxies and their halo gas fractions can be tested with measurements of these diagnostics for galaxies with diverse SFRs but similar halo masses.
△ Less
Submitted 4 June, 2019; v1 submitted 17 October, 2018;
originally announced October 2018.
-
Evolution of the cold gas properties of simulated post-starburst galaxies
Authors:
Timothy A. Davis,
Freeke van de Voort,
Kate Rowlands,
Stuart McAlpine,
Vivienne Wild,
Robert A. Crain
Abstract:
Post-starburst galaxies are typically considered to be a transition population, en route to the red sequence after a recent quenching event. Despite this, recent observations have shown that these objects typically have large reservoirs of cold molecular gas. In this paper we study the star-forming gas properties of a large sample of post-starburst galaxies selected from the cosmological, hydrodyn…
▽ More
Post-starburst galaxies are typically considered to be a transition population, en route to the red sequence after a recent quenching event. Despite this, recent observations have shown that these objects typically have large reservoirs of cold molecular gas. In this paper we study the star-forming gas properties of a large sample of post-starburst galaxies selected from the cosmological, hydrodynamical EAGLE simulations. These objects resemble observed high-mass post-starburst galaxies both spectroscopically and in terms of their space density, stellar mass distribution and sizes. We find that the vast majority of simulated post-starburst galaxies have significant gas reservoirs, with star-forming gas masses of ~10$^9$ M$_{\odot}$, in good agreement with those seen in observational samples. The simulation reproduces the observed time evolution of the gas fraction of the post-starburst galaxy population, with the average galaxy losing ~90 per cent of its star-forming interstellar medium in only ~600 Myr. A variety of gas consumption/loss processes are responsible for this rapid evolution, including mergers and environmental effects, while active galactic nuclei play only a secondary role. The fast evolution in the gas fraction of post-starburst galaxies is accompanied by a clear decrease in the efficiency of star formation, due to a decrease in the dense gas fraction. We predict that forthcoming ALMA observations of the gas reservoirs of low-redshift post-starburst galaxies will show that the molecular gas is typically compact and has disturbed kinematics, reflecting the disruptive nature of many of the evolutionary pathways that build up the post-starburst galaxy population.
△ Less
Submitted 15 January, 2019; v1 submitted 1 October, 2018;
originally announced October 2018.
-
The aftermath of the Great Collision between our Galaxy and the Large Magellanic Cloud
Authors:
Marius Cautun,
Alis J. Deason,
Carlos S. Frenk,
Stuart McAlpine
Abstract:
The Milky Way (MW) offers a uniquely detailed view of galactic structure and is often regarded as a prototypical spiral galaxy. But recent observations indicate that the MW is atypical: it has an undersized supermassive black hole at its centre; it is surrounded by a very low mass, excessively metal-poor stellar halo; and it has an unusually large nearby satellite galaxy, the Large Magellanic Clou…
▽ More
The Milky Way (MW) offers a uniquely detailed view of galactic structure and is often regarded as a prototypical spiral galaxy. But recent observations indicate that the MW is atypical: it has an undersized supermassive black hole at its centre; it is surrounded by a very low mass, excessively metal-poor stellar halo; and it has an unusually large nearby satellite galaxy, the Large Magellanic Cloud (LMC). Here we show that the LMC is on a collision course with the MW with which it will merge in $2.4^{+1.2}_{-0.8}$ Gyrs (68% confidence level). This catastrophic and long-overdue event will restore the MW to normality. Using the EAGLE galaxy formation simulation, we show that, as a result of the merger, the central supermassive black hole will increase in mass by up to a factor of 8. The Galactic stellar halo will undergo an equally impressive transformation, becoming 5 times more massive. The additional stars will come predominantly from the disrupted LMC, but a sizeable number will be ejected onto the halo from the stellar disc. The post-merger stellar halo will have the median metallicity of the LMC, [Fe/H]=-0.5 dex, which is typical of other galaxies of similar mass to the MW. At the end of this exceptional event, the MW will become a true benchmark for spiral galaxies, at least temporarily.
△ Less
Submitted 17 January, 2019; v1 submitted 24 September, 2018;
originally announced September 2018.
-
The rapid growth phase of supermassive black holes
Authors:
Stuart McAlpine,
Richard G. Bower,
David J. Rosario,
Robert A. Crain,
Joop Schaye,
Tom Theuns
Abstract:
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 galaxy host: in sub $L_{*}$ galaxies efficient stellar feedback regulates the gas inflow onto the galaxy and significantly reduces the growth of the ce…
▽ More
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 galaxy 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 ($μ\geq \frac{1}{4}$, where $μ$ 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 μ< \frac{1}{4}$) have a substantially lower influence in triggering the rapid growth phase at all epochs.
△ Less
Submitted 7 September, 2018; v1 submitted 21 May, 2018;
originally announced May 2018.
-
The abundances and properties of Dual AGN and their host galaxies in the EAGLE simulations
Authors:
Yetli Rosas-Guevara,
Richard Bower,
Stuart McAlpine,
Silvia Bonoli,
Patricia Tissera
Abstract:
We look into the abundance of Dual AGN in the largest hydrodynamical simulation from the EAGLE project. We define a Dual AGN as two active black holes (BHs) with a separation below 30 kpc. We find that only 1 per cent of AGN with $\LhX\geq 10^{42}\ergs$ are part of a Dual AGN system at $z=0.8-1$. During the evolution of a typical binary BH system, the rapid variability of the hard X-ray luminosity…
▽ More
We look into the abundance of Dual AGN in the largest hydrodynamical simulation from the EAGLE project. We define a Dual AGN as two active black holes (BHs) with a separation below 30 kpc. We find that only 1 per cent of AGN with $\LhX\geq 10^{42}\ergs$ are part of a Dual AGN system at $z=0.8-1$. During the evolution of a typical binary BH system, the rapid variability of the hard X-ray luminosity on Myr time-scales severely limits the detectability of Dual AGN. To quantify this effect, we calculate a probability of detection, $t_{\rm on}/t_{\rm 30}$, where $t_{\rm 30}$ is the time in which the two black holes are separated at distances below 30 kpc and $t_{\rm on}$, the time that both AGN are visible (e.g. when both AGN have $\LhX\geq 10^{42}\ergs$) in this period. We find that the average fraction of visible Dual systems is 3 per cent. The visible Dual AGN distribution as a function of BH separation presents a pronounced peak at $\sim 20$ kpc that can be understood as a result of the rapid orbital decay of the host galaxies after their first encounter. We also find that $75$ per cent of the host galaxies have recently undergone or are undergoing a merger with stellar mass ratio $\geq 0.1$. Finally, we find that the fraction of visible Dual AGN increases with redshift as found in observations.
△ Less
Submitted 27 November, 2018; v1 submitted 3 May, 2018;
originally announced May 2018.
-
Data release of UV to submm broadband fluxes for simulated galaxies from the EAGLE project
Authors:
Peter Camps,
Ana Trčka,
James Trayford,
Maarten Baes,
Tom Theuns,
Robert A. Crain,
Stuart McAlpine,
Matthieu Schaller,
Joop Schaye
Abstract:
We present dust-attenuated and dust emission fluxes for sufficiently resolved galaxies in the EAGLE suite of cosmological hydrodynamical simulations, calculated with the SKIRT radiative transfer code. The post-processing procedure includes specific components for star formation regions, stellar sources, and diffuse dust, and takes into account stochastic heating of dust grains to obtain realistic…
▽ More
We present dust-attenuated and dust emission fluxes for sufficiently resolved galaxies in the EAGLE suite of cosmological hydrodynamical simulations, calculated with the SKIRT radiative transfer code. The post-processing procedure includes specific components for star formation regions, stellar sources, and diffuse dust, and takes into account stochastic heating of dust grains to obtain realistic broad-band fluxes in the wavelength range from ultraviolet to sub-millimeter. The mock survey includes nearly half a million simulated galaxies with stellar masses above 10^8.5 solar masses across six EAGLE models. About two thirds of these galaxies, residing in 23 redshift bins up to z=6, have a sufficiently resolved metallic gas distribution to derive meaningful dust attenuation and emission, with the important caveat that the same dust properties were used at all redshifts. These newly released data complement the already publicly available information about the EAGLE galaxies, which includes intrinsic properties derived by aggregating the properties of the smoothed particles representing matter in the simulation. We further provide an open source framework of Python procedures for post-processing simulated galaxies with the radiative transfer code SKIRT. The framework allows any third party to calculate synthetic images, SEDs, and broadband fluxes for EAGLE galaxies, taking into account the effects of dust attenuation and emission.
△ Less
Submitted 15 December, 2017;
originally announced December 2017.
-
Identifying the subtle signatures of feedback from distant AGN using ALMA observations and the EAGLE hydrodynamical simulations
Authors:
J. Scholtz,
D. M. Alexander,
C. M. Harrison,
D. J. Rosario,
S. McAlpine,
J. R Mullaney,
F. Stanley,
J. Simpson,
T. Theuns,
R. G. Bower,
R. C. Hickox,
P. Santini,
A. M. Swinbank
Abstract:
We present sensitive 870$μ$m continuum measurements from our ALMA programmes of 114 X-ray selected AGN in the CDF-S and COSMOS fields. We use these observations in combination with data from Spitzer and Herschel to construct a sample of 86 X-ray selected AGN, 63 with ALMA constraints at $z=1.5-3.2$ with stellar mass $>2\times10^{10}M_{\odot}$. We constructed broad-band spectral energy distribution…
▽ More
We present sensitive 870$μ$m continuum measurements from our ALMA programmes of 114 X-ray selected AGN in the CDF-S and COSMOS fields. We use these observations in combination with data from Spitzer and Herschel to construct a sample of 86 X-ray selected AGN, 63 with ALMA constraints at $z=1.5-3.2$ with stellar mass $>2\times10^{10}M_{\odot}$. We constructed broad-band spectral energy distributions in the infrared band (8 - 1000$μ$m) and constrain star-formation rates (SFRs) uncontaminated by the AGN. Using a hierarchical Bayesian method that takes into account the information from upper limits, we fit SFR and specific SFR (sSFR) distributions. We explore these distributions as a function of both X-ray luminosity and stellar mass. We compare our measurements to two versions of the EAGLE hydrodynamical simulations: the reference model with AGN feedback and the model without AGN. We find good agreement between the observations and that predicted by the EAGLE reference model for the modes and widths of the sSFR distributions as a function of both X-ray luminosity and stellar mass; however, we found that the EAGLE model without AGN feedback predicts a significantly narrower width when compared to the data. Overall, from the combination of the observations with the model predictions, we conclude that (1) even with AGN feedback, we expect no strong relationship between the sSFR distribution parameters and instantaneous AGN luminosity and (2) a signature of AGN feedback is a broad distribution of sSFRs for all galaxies (not just those hosting an AGN) with stellar masses above $\approx 10^{10}$M$_{\odot}$.
△ Less
Submitted 7 December, 2017;
originally announced December 2017.
-
The link between galaxy and black hole growth in the EAGLE simulation
Authors:
Stuart McAlpine,
Richard G. Bower,
Chris M. Harrison,
Robert A. Crain,
Matthieu Schaller,
Joop Schaye,
Tom Theuns
Abstract:
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 studies, that the <SFR>--BHAR relation for an AGN selected sample produces a relatively flat trend, whilst the <BHAR>--SFR relation for a SFR selected sam…
▽ More
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 studies, 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.
△ Less
Submitted 17 March, 2017; v1 submitted 4 January, 2017;
originally announced January 2017.
-
A chronicle of galaxy mass assembly in the EAGLE simulation
Authors:
Yan Qu,
John C. Helly,
Richard G. Bower,
Tom Theuns,
Robert A. Crain,
Carlos S. Frenk,
Michelle Furlong,
Stuart McAlpine,
Matthieu Schaller,
Joop Schaye,
Simon D. M. White
Abstract:
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…
▽ More
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.
△ Less
Submitted 23 September, 2016;
originally announced September 2016.
-
The dark nemesis of galaxy formation: why hot haloes trigger black hole growth and bring star formation to an end
Authors:
Richard Bower,
Joop Schaye,
Carlos S. Frenk,
Tom Theuns,
Matthieu Schaller,
Robert A. Crain,
Stuart McAlpine
Abstract:
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…
▽ More
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.
△ Less
Submitted 4 November, 2016; v1 submitted 25 July, 2016;
originally announced July 2016.
-
Supermassive black holes in the EAGLE Universe. Revealing the observables of their growth
Authors:
Yetli Rosas-Guevara,
Richard G. Bower,
Joop Schaye,
Stuart McAlpine,
Claudio Dalla-Vecchia,
Carlos S. Frenk,
Matthieu Schaller,
Tom Theuns
Abstract:
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 funct…
▽ More
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.
△ Less
Submitted 8 July, 2016; v1 submitted 31 March, 2016;
originally announced April 2016.
-
Music from the heavens - Gravitational waves from supermassive black hole mergers in the EAGLE simulations
Authors:
Jaime Salcido,
Richard G. Bower,
Tom Theuns,
Stuart McAlpine,
Matthieu Schaller,
Robert A. Crain,
Joop Schaye,
John Regan
Abstract:
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 rad…
▽ More
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.
△ Less
Submitted 29 July, 2016; v1 submitted 21 January, 2016;
originally announced January 2016.
-
The EAGLE simulations of galaxy formation: public release of halo and galaxy catalogues
Authors:
Stuart McAlpine,
John C. Helly,
Matthieu Schaller,
James W. Trayford,
Yan Qu,
Michelle Furlong,
Richard G. Bower,
Robert A. Crain,
Joop Schaye,
Tom Theuns,
Claudio Dalla Vecchia,
Carlos S. Frenk,
Ian G. McCarthy,
Adrian Jenkins,
Yetli Rosas-Guevara,
Simon D. M. White,
Maarten Baes,
Peter Camps,
Gerard Lemson
Abstract:
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 radiativ…
▽ More
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)
△ Less
Submitted 25 August, 2016; v1 submitted 5 October, 2015;
originally announced October 2015.