Search for: All records

We examine the dual [both black hole (BH) active] and offset (one BH active and in distinct galaxies) active galactic nucleus (AGN) population (comprising ∼ 2000 pairs at $0.5\, \text{kpc}\lesssim \Delta r\lt 30\, \text{kpc}$) at z = 2 ∼ 3 in the ASTRID simulation covering (360 cMpc)3. The dual (offset) AGN make up 3.0(0.5) per cent of all AGN at z = 2. The dual fraction is roughly constant while the offset fraction increases by a factor of 10 from z = 4 ∼ 2. Compared with the full AGN population, duals are characterized by low MBH/M* ratios, high specific star formation rates (sSFR) of $\sim 1\, \text{Gyr}^{-1}$, and high Eddington ratios (∼0.05, double that of single AGN). Dual AGNs are formed in major galaxy mergers (typically involving $M_\text{halo}\lt 10^{13}\, M_\odot$), with simular-mass BHs. At small separations (when host galaxies are in the late phase of the merger), duals become 2 ∼ 8 times brighter (albeit more obscured) than at larger separations. 80  per cent of the bright, close duals would merge within $\sim 500\, \text{Myr}$. Notably, the initially less-massive BHs in duals frequently become the brighter AGN during galaxy mergers. In offset AGN, the active BH is typically ≳ 10 times more massive than its non-activemore »counterpart and than most BHs in duals. Offsets are predominantly formed in minor galaxy mergers with the active BH residing in the centre of massive haloes ($M_\text{ halo}\sim 10^{13-14}\, \mathrm{M}_\odot$). In these deep potentials, gas stripping is common and the secondary quickly deactivates. The stripping also leads to inefficient orbital decay amongst offsets, which stall at $\Delta r\sim 5\, \text{kpc}$ for a few hundred Myrs.

« less

We introduce the Astrid  simulation, a large-scale cosmological hydrodynamic simulation in a $250 \, h^{-1}\mathrm{Mpc}$ box with 2 × 55003 particles. Astrid contains a large number of high redshift galaxies, which can be compared to future survey data, and resolves galaxies in haloes more massive than $2\times 10^9 \, \mathrm{M}_{\odot }$. Astrid  has been run from z = 99 to 3. As a particular focus is modelling the high redshift Universe, it contains models for inhomogeneous hydrogen and helium reionization, baryon relative velocities and massive neutrinos, as well as supernova and AGN feedback. The black hole model includes mergers driven by dynamical friction rather than repositioning. We briefly summarize the implemented models, and the technical choices we took when developing the simulation code. We validate the model, showing good agreement with observed ultraviolet luminosity functions, galaxy stellar mass functions and specific star formation rates (SFRs). We show that the redshift at which a given galaxy underwent hydrogen reionization has a large effect on the halo gas fraction. Finally, at z = 6, haloes with $M \sim 2\times 10^9 \, \mathrm{M}_{\odot }$ which have been reionized have an SFR 1.5 times greater than those which have not yet been reionized.

The launch of space-based gravitational-wave (GW) detectors (e.g. Laser Interferometry Space Antenna; LISA) and current and upcoming Pulsar Timing Arrays will extend the GW window to low frequencies, opening new investigations into dynamical processes involving massive black hole binaries (MBHBs) and their mergers across cosmic time. MBHBs are expected to be among the primary sources for the upcoming low-frequency (10−4–10−1 Hz) window probed by LISA. It is important to investigate the expected supermassive BH merger rates and associated signals, to determine how potential LISA events are affected by physics included in current models. To study this, we post-process the large population of MBHBs in the Illustris simulation to account for dynamical friction time delays associated with BH infall/inspiral. We show that merger delays associated with binary evolution have the potential to decrease the expected merger rates, with $M_{\rm {BH}}\ \gt\ 10^6\ \mathrm{M}_\odot$ MBHBs (the lowest mass in Illustris) decreasing from ∼3 to ∼0.1 yr−1, and shifting the merger peak from z ∼2 to ∼1.25. During this time, we estimate that accretion grows the total merging mass by as much as 7x the original mass. Importantly, however, dynamical friction-associated delays (which shift the mergers toward lower redshift and higher masses) leadmore »to a stronger signal/strain for the emitted GWs in the LISA band, increasing mean frequency from 10−3.1 to 10−3.4–10−4.0 Hz, and mean strain from 10−17.2 to 10−16.3–10−15.3. Finally, we show that after including a merger delay and associated MBH growth, mergers still tend to lie on the typical MBH–M* relation, but with an increased likelihood of an undermassive BH.

« less

A population of supermassive black hole (SMBH) binaries is expected to generate a stochastic gravitational wave background (SGWB) in the pulsar timing array (PTA) frequency range of 10−9 to $10^{-7}\, {\rm Hz}$. Detection of this signal is a current observational goal and so predictions of its characteristics are of significant interest. In this work, we use SMBH binary mergers from the MassiveBlackII simulation to estimate the characteristic strain of the stochastic background. We examine both a gravitational wave (GW) driven model of binary evolution and a model which also includes the effects of stellar scattering and a circumbinary gas disc. Results are consistent with PTA upper limits and similar to estimates in the literature. The characteristic strain at a reference frequency of $1\, {\rm yr}^{-1}$ is found to be $A_{\rm {yr}^{-1}}= 6.9 \times 10^{-16}$ and $A_{\rm {yr}^{-1}}= 6.4 \times 10^{-16}$ in the GW-driven and stellar scattering/gas disc cases, respectively. Using the latter approach, our models show that the SGWB is mildly suppressed compared to the purely GW-driven model as frequency decreases inside the PTA frequency band.

In this work, we establish and test methods for implementing dynamical friction (DF) for massive black hole pairs that form in large volume cosmological hydrodynamical simulations that include galaxy formation and black hole growth. We verify our models and parameters both for individual black hole dynamics and for the black hole population in cosmological volumes. Using our model of DF from collisionless particles, black holes can effectively sink close to the galaxy centre, provided that the black hole’s dynamical mass is at least twice that of the lowest mass resolution particles in the simulation. Gas drag also plays a role in assisting the black holes’ orbital decay, but it is typically less effective than that from collisionless particles, especially after the first billion years of the black hole’s evolution. DF from gas becomes less than $1{{\ \rm per\ cent}}$ of DF from collisionless particles for BH masses >107 M⊙. Using our best DF model, we calculate the merger rate down to z = 1.1 using an Lbox = 35 Mpc h−1 simulation box. We predict ∼2 mergers per year for z > 1.1 peaking at z ∼ 2. These merger rates are within the range obtained in previous work using similar resolution hydrodynamical simulations.more »We show that the rate is enhanced by factor of ∼2 when DF is taken into account in the simulations compared to the no-DF run. This is due to ${\gt}40{{\ \rm per\ cent}}$ more black holes reaching the centre of their host halo when DF is added.

« less

We study the sizes of galaxies in the Epoch of Reionization using a sample of ${\sim 100\, 000}$ galaxies from the BlueTides cosmological hydrodynamical simulation from z = 7 to 11. We measure the galaxy sizes from stellar mass and luminosity maps, defining the effective radius as the minimum radius that could enclose the pixels containing 50 per cent of the total mass/light in the image. We find an inverse relationship between stellar mass and effective half-mass radius, suggesting that the most massive galaxies are more compact and dense than lower mass galaxies, which have flatter mass distributions. We find a mildly negative relation between intrinsic far-ultraviolet luminosity and size, while we find a positive size–luminosity relation when measured from dust-attenuated images. This suggests that dust is the predominant cause of the observed positive size–luminosity relation, with dust preferentially attenuating bright sightlines resulting in a flatter emission profile and thus larger measured effective radii. We study the size–luminosity relation across the rest-frame ultraviolet and optical, and find that the slope decreases at longer wavelengths; this is a consequence of the relation being caused by dust, which produces less attenuation at longer wavelengths. We find that the far-ultraviolet size–luminosity relation shows mild evolutionmore »from z = 7 to 11, and galaxy size evolves with redshift as R ∝ (1 + z)−m, where m = 0.662 ± 0.009. Finally, we investigate the sizes of z = 7 quasar host galaxies, and find that while the intrinsic sizes of quasar hosts are small relative to the overall galaxy sample, they have comparable sizes when measured from dust-attenuated images.

« less

ABSTRACT In this work, we expand and test the capabilities of our recently developed superresolution (SR) model to generate high-resolution (HR) realizations of the full phase-space matter distribution, including both displacement and velocity, from computationally cheap low-resolution (LR) cosmological N-body simulations. The SR model enhances the simulation resolution by generating 512 times more tracer particles, extending into the deeply nonlinear regime where complex structure formation processes take place. We validate the SR model by deploying the model in 10 test simulations of box size 100 h−1 Mpc, and examine the matter power spectra, bispectra, and two-dimensional power spectra in redshift space. We find the generated SR field matches the true HR result at per cent level down to scales of k ∼ 10 h  Mpc−1. We also identify and inspect dark matter haloes and their substructures. Our SR model generates visually authentic small-scale structures that cannot be resolved by the LR input, and are in good statistical agreement with the real HR results. The SR model performs satisfactorily on the halo occupation distribution, halo correlations in both real and redshift space, and the pairwise velocity distribution, matching the HR results with comparable scatter, thus demonstrating its potential in making mock halo catalogues. The SR technique can be amore »powerful and promising tool for modelling small-scale galaxy formation physics in large cosmological volumes.« less

ABSTRACT We study the alignments of satellite galaxies, and their anisotropic distribution, with respect to location and orientation of their host central galaxy in MassiveBlack-II (MB-II) and IllustrisTNG simulations. We find that: the shape of the satellite system in haloes of mass ($\gt 10^{13}\, h^{-1}\, \mathrm{M}_{\odot }$) is well aligned with the shape of the central galaxy at z = 0.06 with the mean alignment between the major axes being ∼Δθ = 12° when compared to a uniform random distribution; that satellite galaxies tend to be anisotropically distributed along the major axis of the central galaxy with a stronger alignment in haloes of higher mass or luminosity; and that the satellite distribution is more anisotropic for central galaxies with lower star formation rate, which are spheroidal, and for red central galaxies. Radially, we find that satellites tend to be distributed along the major axis of the shape of the stellar component of central galaxies at smaller scales and the dark matter component on larger scales. We find that the dependence of satellite anisotropy on central galaxy properties and the radial distance is similar in both the simulations with a larger amplitude in MB-II. The orientation of satellite galaxies tends tomore »point toward the location of the central galaxy at small scales and this correlation decreases with increasing distance, and the amplitude of satellite alignment is higher in high-mass haloes. However, the projected ellipticities do not exhibit a scale-dependent radial alignment, as has been seen in some observational measurements.« less

Cosmological simulations of galaxy formation are limited by finite computational resources. We draw from the ongoing rapid advances in artificial intelligence (AI; specifically deep learning) to address this problem. Neural networks have been developed to learn from high-resolution (HR) image data and then make accurate superresolution (SR) versions of different low-resolution (LR) images. We apply such techniques to LR cosmological N-body simulations, generating SR versions. Specifically, we are able to enhance the simulation resolution by generating 512 times more particles and predicting their displacements from the initial positions. Therefore, our results can be viewed as simulation realizations themselves, rather than projections, e.g., to their density fields. Furthermore, the generation process is stochastic, enabling us to sample the small-scale modes conditioning on the large-scale environment. Our model learns from only 16 pairs of small-volume LR-HR simulations and is then able to generate SR simulations that successfully reproduce the HR matter power spectrum to percent level up to$16\hspace{0.17em}{h}^{-1}\mathrm{\text{Mpc}}$and the HR halo mass function to within$10\%$down to$10^{11}\hspace{0.17em}{M}_{\odot }$. We successfully deploy the model in a box 1,000 times larger than the training simulation box, showing that high-resolution mock surveys can be generated rapidly. We conclude that AI assistance has themore »potential to revolutionize modeling of small-scale galaxy-formation physics in large cosmological volumes.

« less