ABSTRACT
We study how statistical properties of supermassive black holes depend on the frequency and conditions for massive seed formation in cosmological simulations of structure formation. We ...develop a novel method to recalculate detailed growth histories and merger trees of black holes within the framework of the Illustris simulation for several seed formation models, including a physically motivated model where black hole seeds only form in progenitor galaxies that conform to the conditions for direct collapse black hole formation. While all seed models considered here are in a broad agreement with present observational constraints on black hole populations from optical, UV, and X-ray studies, we find that they lead to widely different black hole number densities and halo occupation fractions, which are currently observationally unconstrained. In terms of future electromagnetic spectrum observations, the faint-end quasar luminosity function and the low-mass-end black hole–host galaxy scaling relations are very sensitive to the specific massive seed prescription. Specifically, the direct collapse model exhibits a seeding efficiency that decreases rapidly with cosmic time and produces much fewer black holes in low-mass galaxies, in contrast to the original Illustris simulation. We further find that the total black hole merger rate varies by more than one order of magnitude for different seed models, with the redshift evolution of the chirp mass changing as well. Supermassive black hole merger detections with LISA and International Pulsar Timing Array may hence provide the most direct means of constraining massive black hole seed formation in the early Universe.
Abstract
We use the high-resolution cosmological simulation Illustris to investigate the clustering of supermassive black holes across cosmic time, the link between black hole clustering and host ...halo masses, and the implications for black hole duty cycles. Our predicted black hole correlation length and bias match the observational data very well across the full redshift range probed. Black hole clustering is strongly luminosity dependent on small, 1-halo scales, with some moderate dependence on larger scales of a few Mpc at intermediate redshifts. We find black hole clustering to evolve only weakly with redshift, initially following the behaviour of their hosts. However, below z ∼ 2 black hole clustering increases faster than that of their hosts, which leads to a significant overestimate of the clustering-predicted host halo mass. The full distribution of host halo masses is very wide, including a low-mass tail extending up to an order of magnitude below the naive prediction for minimum host mass. Our black hole duty cycles, f
duty, follow a power-law dependence on black hole mass and decrease with redshift, and we provide accurate analytic fits to these. The increase in clustering amplitude at late times, however, means that duty cycle estimates based on black hole clustering can overestimate f
duty substantially, by more than two orders of magnitude. We find the best agreement when the minimum host mass is assumed to be 1011.2 M⊙, which provides an accurate measure across all redshifts and luminosity ranges probed by our simulation.
Previous simulations of the growth of cosmic structures have broadly reproduced the 'cosmic web' of galaxies that we see in the Universe, but failed to create a mixed population of elliptical and ...spiral galaxies, because of numerical inaccuracies and incomplete physical models. Moreover, they were unable to track the small-scale evolution of gas and stars to the present epoch within a representative portion of the Universe. Here we report a simulation that starts 12 million years after the Big Bang, and traces 13 billion years of cosmic evolution with 12 billion resolution elements in a cube of 106.5 megaparsecs a side. It yields a reasonable population of ellipticals and spirals, reproduces the observed distribution of galaxies in clusters and characteristics of hydrogen on large scales, and at the same time matches the 'metal' and hydrogen content of galaxies on small scales.
Celotno besedilo
Dostopno za:
DOBA, IJS, IZUM, KILJ, KISLJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
We present the full public release of all data from the Illustris simulation project. Illustris is a suite of large volume, cosmological hydrodynamical simulations run with the moving-mesh code Arepo ...and including a comprehensive set of physical models critical for following the formation and evolution of galaxies across cosmic time. Each simulates a volume of (106.5 Mpc)3 and self-consistently evolves five different types of resolution elements from a starting redshift of z=127 to the present day, z=0. These components are: dark matter particles, gas cells, passive gas tracers, stars and stellar wind particles, and supermassive black holes. This data release includes the snapshots at all 136 available redshifts, halo and subhalo catalogs at each snapshot, and two distinct merger trees. Six primary realizations of the Illustris volume are released, including the flagship Illustris-1 run. These include three resolution levels with the fiducial “full” baryonic physics model, and a dark matter only analog for each. In addition, we provide four distinct, high time resolution, smaller volume “subboxes”. The total data volume is ∼265 TB, including ∼800 full volume snapshots and ∼30,000 subbox snapshots. We describe the released data products as well as tools we have developed for their analysis. All data may be directly downloaded in its native HDF5 format. Additionally, we release a comprehensive, web-based API which allows programmatic access to search and data processing tasks. In both cases we provide example scripts and a getting-started guide in several languages: currently, IDL, Python, and Matlab. This paper addresses scientific issues relevant for the interpretation of the simulations, serves as a pointer to published and on-line documentation of the project, describes planned future additional data releases, and discusses technical aspects of the release.
We present relations between X-ray luminosity and velocity dispersion (L − σ), X-ray luminosity and gas mass (L − Mgas), and cluster radius and velocity dispersion (r500 − σ) for 62 galaxy clusters ...in the HIFLUGCS, an X-ray flux-limited sample minimizing bias toward any cluster morphology. Our analysis in total is based on ~1.3 Ms of clean X-ray XMM-Newton data and 13439 cluster member galaxies with redshifts. Cool cores are among the major contributors to the scatter in the L − σ relation. When the cool-core-corrected X-ray luminosity is used the intrinsic scatter decreases to 0.27 dex. Even after the X-ray luminosity is corrected for the cool core, the scatter caused by the presence of cool cores dominates for the low-mass systems. The scatter caused by the non-cool-core clusters does not strongly depend on the mass range, and becomes dominant in the high-mass regime. The observed L − σ relation agrees with the self-similar prediction, matches that of a simulated sample with AGN feedback disregarding six clusters with <45 cluster members with spectroscopic redshifts, and shows a common trend of increasing scatter toward the low-mass end, i.e., systems with σ ≤ 500 kms-1. A comparison of observations with simulations indicates an AGN-feedback-driven impact in the low-mass regime. The best fits to the L − Mgas relations for the disturbed clusters and undisturbed clusters in the observational sample closely match those of the simulated samples with and without AGN feedback, respectively. This suggests that one main cause of the scatter is AGN activity providing feedback in different phases, e.g. during a feedback cycle. The slope and scatter in the observed r500 − σ relation is similar to that of the simulated sample with AGN feedback except for a small offset but still within the scatter.
We explore how radiative cooling, supernova feedback, cosmic rays, and a new model of the energetic feedback from active galactic nuclei (AGNs) affect the thermal and kinetic Sunyaev-Zel'dovich (SZ) ...power spectra. To do this, we use a suite of hydrodynamical TreePM-SPH simulations of the cosmic web in large periodic boxes and tailored higher resolution simulations of individual galaxy clusters. Our AGN feedback simulations match the recent universal pressure profile and cluster mass scaling relations of the REXCESS X-ray cluster sample better than previous analytical or numerical approaches. For multipoles l 2000, our power spectra with and without enhanced feedback are similar, suggesting that theoretical uncertainties over that range are relatively small, although current analytic and semi-analytic approaches overestimate this SZ power. We find the power at high 2000-1,000 multipoles in which the Atacama Cosmology Telescope (ACT) and South Pole Telescope (SPT) probe is sensitive to the feedback prescription, and hence can constrain the theory of intracluster gas, in particular for the highly uncertain redshifts >0.8. The apparent tension between Delta *s8 from primary cosmic microwave background power and from analytic SZ spectra inferred using ACT and SPT data is lessened with our AGN feedback spectra.
ABSTRACT
We analyse a hydrodynamical simulation model for the recurrent heating of the central intra‐cluster medium (ICM) by active galactic nuclei (AGN). Besides the self‐gravity of the dark matter ...and gas components, our approach includes the radiative cooling and photoheating of the gas, as well as a subresolution multiphase model for star formation and supernova feedback. Additionally, we incorporate a periodic heating mechanism in the form of hot, buoyant bubbles, injected into the intragalactic medium (IGM) during the active phases of the accreting central AGN. We use simulations of isolated cluster haloes of different masses to study the bubble dynamics and heat transport into the IGM. We also apply our model to self‐consistent cosmological simulations of the formation of galaxy clusters with a range of masses. Our numerical schemes explore a variety of different assumptions for the spatial configuration of AGN‐driven bubbles, for their duty cycles and for the energy injection mechanism, in order to obtain better constraints on the underlying physical picture. We argue that AGN heating can substantially affect the properties of both the stellar and gaseous components of clusters of galaxies. Most importantly, it alters the properties of the central dominant (cD) galaxy by reducing the mass deposition rate of freshly cooled gas out of the ICM, thereby offering an energetically plausible solution to the cooling‐flow problem. At the same time, this leads to reduced or eliminated star formation in the central cD galaxy, giving it red stellar colours as observed.
We present a ‘two-fluid’ implementation of dust in smoothed particle hydrodynamics (SPH) in the test particle limit. The scheme is able to handle both short and long stopping times and reproduces the ...short friction time limit, which is not properly handled in other implementations. We apply novel tests to verify its accuracy and limitations, including multidimensional tests that have not been previously applied to the drag-coupled dust problem and which are particularly relevant to self-gravitating protoplanetary discs. Our tests demonstrate several key requirements for accurate simulations of gas–dust mixtures. First, in standard SPH particle jitter can degrade the dust solution, even when the gas density is well reproduced. The use of integral gradients, a Wendland kernel and a large number of neighbours can control this, albeit at a greater computational cost. Secondly, when it is necessary to limit the artificial viscosity we recommend using the Cullen & Dehnen switch, since the alternative, using α ∼ 0.1, can generate a large velocity noise up to σ
v
≲ 0.3c
s in the dust particles. Thirdly, we find that an accurate dust density estimate requires >400 neighbours, since, unlike the gas, the dust particles do not feel regularization forces. This density noise applies to all particle-based two-fluid implementations of dust, irrespective of the hydro solver and could lead to numerically induced fragmentation. Although our tests show accurate dusty gas simulations are possible, care must be taken to minimize the contribution from numerical noise.