Search

Detecting dual active galactic nuclei (DAGN) in observations and understanding theoretically which massive black holes (MBHs) compose them and in which galactic and large-scale environment they reside are becoming increasingly important questions as we enter the multi-messenger era of MBH astronomy. This paper presents the abundance and properties of DAGN produced in nine large-scale cosmological hydrodynamical simulations. We focus on DAGN powered by AGN with Lbol > 1e43 erg/s and belonging to distinct galaxies, i.e. pairs that can be characterised with current and near-future electromagnetic observations. We find that the number density of DAGN separated by a few to 30 proper kpc varies from 1e-8 (or none) to 1e-3 comoving Mpc^3 in the redshift range z=0-7. At a given redshift, the densities of the DAGN numbers vary by up to two orders of magnitude from one simulation to another. However, for all simulations, the DAGN peak is in the range z=1-3, right before the peak of cosmic star formation or cosmic AGN activity. The corresponding fractions of DAGN (with respect to the total number of AGN) range from 0 to 6 percent. We find that simulations could produce too few DAGN at z=0 (or merge pairs too quickly) compared to current observational constraints while being consistent with preliminary constraints at high redshift (z = 3). Next-generation observatories (e.g., AXIS) will be of paramount importance to detect DAGN across cosmic times. We predict the detectability of DAGN with future X-ray telescopes and discuss DAGN as progenitors for future LISA gravitational wave detections., Comment: Accepted for publication in MNRAS

We study the co-evolution of black holes (BHs) and their host galaxies in the ASTRID and Illustris-TNG300 cosmological simulations and the Dark Sage Semi-Analytic Model (SAM), focusing on the evolution of the BH mass - stellar mass ($M_{\rm BH}-M_*$) relation. Due to differences in the adopted sub-grid modeling of BH seeding, dynamics, and feedback, the models differ in their predicted redshift evolution of the $M_{\rm BH}-M_*$ relation. We find that it is the interplay between the star formation rate (SFR) and the black hole accretion rate (BHAR) which drives the evolution of the mean relation. We define a quantity $\mathcal{R}$, the ratio between the specific BHAR and SFR (i.e. $\mathcal{R} \equiv\ $sBHAR/sSFR), and demonstrate that it is $\mathcal{R}$ that governs the evolution of individual sources in the $M_{\rm BH}-M_*$ plane. The efficiency of BH growth versus stellar mass growth in the sSFR-sBHAR plane reflects the partitioning of gas between fueling star formation versus BH accretion. This partitioning depends on the implementation of BH dynamics and the nature of how AGN feedback quenches galaxies. In the cosmological simulations (ASTRID and Illustris-TNG300), the BHAR and SFR are intrinsically linked, resulting in a tight $M_{\rm BH}-M_*$ correlation, while the Dark Sage SAM produces a significantly larger scatter. We discuss these results in the context of recently discovered over-massive BHs and massive quenched galaxies at high redshift by the James Webb Space Telescope., Comment: 17 pages, 8 figures. Submitted to ApJ. Comments welcome

MBH seed mergers are expected to be among the loudest sources of gravitational waves detected by the Laser Interferometer Space Antenna (LISA), providing a unique window into the birth and early growth of SMBH. We present the MAGICS-II simulation suite, consisting of 6 galaxy mergers that result in MBH seeds mergers identified in the cosmological simulation ASTRID. With the enhanced resolution (mass resolution: 500 $\mathrm{M}_{\odot}$; softening length: $5$ pc), improved subgrid models for the MBH dynamics and accretion, and the accurate regularized gravity integrator included in KETJU, we trace MBH seeds dynamics down to 0.1 pc. After evolving all the systems for $\approx$ 1.2 Gyr in three stages (MAGICS-2000, MAGICS-500, and MAGICS-K), we find in 4 of the 6 systems the MBHs stall at separations $ \Delta r \gtrsim 200$ pc. Only in 2 systems, the MBHs manage to sink further, and only in one of them a bound binary forms. In the sinking systems, the MBH retains a population of bound stars. The final separation between the MBH is related to the surrounding unstripped stellar (and/or dark matter) mass: if more than 90\% of the surrounding stellar system is stripped away, the MBHs stall. Besides the unstripped stars from the original host galaxy, we find that newly formed stars bound to the MBH significantly contribute to its sinking. Resolving the stellar system around MBH seeds, and its induced tidal interactions and dynamical friction is key for accurately capturing MBH dynamics. For this, high resolution simulations are required. In a companion paper (MAGICS-III), we resimulate the central regions of these systems with increased resolution to model directly the effects of actual star clusters around MBHs., Comment: 26 pages, 20 figures. Submitted to ApJ. Comments welcome!

Merger rate predictions of Massive Black Hole (MBH) seeds from large-scale cosmological simulations differ widely, with recent studies highlighting the challenge of low-mass MBH seeds failing to reach the galactic center, a phenomenon known as the seed sinking problem. In this work, we tackle this issue by integrating cosmological simulations and galaxy merger simulations from the MAGICS-I and MAGICS-II resimulation suites with high-resolution $N$-body simulations. Building on the findings of MAGICS-II, which showed that only MBH seeds embedded in stellar systems are able to sink to the center, we extend the investigation by incorporating nuclear star clusters (NSCs) into our models. Utilizing $N$-body resimulations with up to $10^7$ particles, we demonstrate that interactions between NSCs and their surrounding galactic environment, particularly tidal forces triggered by cluster interactions, significantly accelerate the sinking of MBHs to the galactic center. This process leads to the formation of a hard binary in $\lesssim 500$ Myr after the onset of a galaxy merger. Our results show that in 8 out of 12 models, the high stellar density of the surrounding NSCs enhances MBH hardening, facilitating gravitational wave (GW) mergers by redshift $z = 4$. We conclude that at $z > 4$, dense NSCs serve as the dominant channel for MBH seed mergers, producing a merger rate of $0.3$--$0.6\, \mathrm{yr}^{-1}$ at $z = 4$, which is approximately 300--600 times higher than in non-NSC environments. In contrast, in environments without NSCs, surrounding dark matter plays a more significant role in loss-cone scattering., Comment: 26 pages, 13 figures. Submitted to ApJ. Comments welcome!

We present new results from the ASTRID simulation from $z=3$ to $z=0.5$, covering the epoch of cosmic noon. The galaxy stellar mass function, as well as the black hole mass and luminosity functions in ASTRID, exhibit good agreement with recent observational constraints. We study the $M_{\rm BH}$-$M_*$ scaling relation and its connections to AGN luminosity, galaxy color, and star formation rate, demonstrating that AGN feedback plays a crucial role in the quenching of massive galaxies ($M_*>10^{10.5} M_{\odot}$). Although AGN feedback suppresses star formation through quenching, AGN-host galaxies still exhibit statistically higher levels of star formation compared to inactive ones, due to the positive correlation between AGN activity and star formation, both fueled by a shared gas reservoir. The fraction of quiescent galaxies in ASTRID increases with both galaxy mass and redshift evolution, aligning well with observational trends. We find that different quenching mechanisms can leave distinct morphological imprints on quenched galaxies. Massive, compact quiescent galaxies typically experience shorter quenching timescales, have younger central regions, and host overmassive black holes. This is usually due to a compaction-like quenching mechanism that funnels gas into the galaxy center, leading to starbursts and triggering AGN kinetic feedback. In contrast, quiescent galaxies with more diffuse morphologies generally experience `inside-out' quenching, which is characterized by older central regions compared to the outskirts. These galaxies typically experience longer quenching timescales due to quenching processes operating on a larger halo scale, which gradually deplete the galactic star-forming gas. Data of the \astrid simulation down to $z=0.5$ is available at \url{https://astrid.psc.edu}., Comment: Summited to APJ, comments are welcome

Kpc-scale dual and offset Active Galactic Nuclei (AGNs) are signposts of accreting supermassive black holes (SMBHs) triggered during late-stage galaxy mergers, offering crucial insights into the coevolution of SMBHs and galaxies. However, robustly confirmed systems at high redshift (e.g., $z>1$) are scarce and biased towards the most luminous and unobscured systems. In this study, we systematically search for kpc-scale (projected separation $<15$ kpc) dual and offset AGNs around 571 moderate-luminosity, X-ray-selected AGNs including the obscured population, utilizing deep HST ACS/F814W and multiband JWST NIRCam imaging from the COSMOS-Web survey. We identify 59 dual and 30 offset AGN candidates in late stage major mergers based on spatially-resolved spectral energy distribution analyses. This translates to $\sim28$ and $\sim10$ bona-fide dual and offset AGNs using a probabilistic pair counting scheme to minimize chance superpositions. Notably, the fraction of dual and offset AGNs among moderate-luminosity ($L_{\rm bol}\sim10^{43}-10^{46}\ \rm erg\ s^{-1}$), obscured AGNs is nearly two orders of magnitude higher than that of the most luminous, unobscured quasar pairs. We find tentative evidence for an increasing pair fraction among AGNs with redshift (from a few percent at $z\sim 0.5$ to $\sim22.9_{-17.7}^{+27.5}\%$ at $z\sim4.5$) and a higher occurrence rate of dual over offset AGNs. There is no pileup of dual/offset AGNs below $\sim 2~{\rm kpc}$ separations. These results generally align with predictions from the ASTRID and Horizon-AGN cosmological simulations when matching sample selection criteria, implying a high probability of both BHs being active simultaneously in late-stage major mergers., Comment: 20 pages, 8 figures, submitted to ApJ

Supermassive black holes (SMBHs) can be ejected from their galactic centers due to gravitational wave recoil or the slingshot mechanism following a galaxy merger. If an ejected SMBH retains its inner accretion disk, it may be visible as an off-nuclear active galactic nucleus (AGN). At present, only a handful of offset AGNs that are recoil or slingshot candidates have been found, and none have been robustly confirmed. Compiling a large sample of runaway SMBHs would enable us to constrain the mass and spin evolution of binary SMBHs and study feedback effects of displaced AGNs. We adapt the method of varstrometry -- which was developed for Gaia observations to identify off-center, dual, and lensed AGNs -- in order to quickly identify off-nuclear AGNs in optical survey data by looking for an excess of blue versus red astrometric jitter. We apply this to the Pan-STARRS1 3$\pi$ Survey and report on five new runaway AGN candidates. We focus on ZTF18aajyzfv: a luminous quasar offset by 6.7 $\pm$ 0.2 kpc from an adjacent galaxy at $z$=0.224, and conclude after Keck LRIS spectroscopy and comparison to ASTRID simulation analogs that it is likely a dual AGN. This selection method can be easily adapted to work with data from the soon-to-be commissioned Vera C. Rubin Telescope Legacy Survey of Space and Time (LSST). LSST will have a higher cadence and deeper magnitude limit than Pan-STARRS1, and should permit detection of many more runaway SMBH candidates., Comment: Accepted for publication by ApJ. 13 pages, 12 figures, 3 tables

We present mock JWST observations for more than 250,000 different galaxies from the Astrid simulation with $3 \leq z \leq 6$. The mock observations are made using the BPASS stellar SED model, and a simple dust model. They are then viewed through NIRCam filters, convolved with a PSF, have noise added, and are drizzled together to emulate the Cosmic Evolution Early Release Science (CEERS) survey. We analyse this dataset by computing a number of morphological measures and find our catalog to have comparable statistics to similar mock catalogs, and the first release of CEERS data. We find that most of the Sersic indices of galaxies in our redshift range are lower than observed, with most having n less than one. Additionally, we observe the sizes of galaxies of all masses to increase from redshift z=6 to redshift z=3 consistent with other results. The number of galaxies in our catalog allows us to examine how relationships like the mass-size relation evolve with redshift, and compare the accuracy of a variety of traditional galaxy classification techniques (Sersic fit, Asymmetry-Concentration, and Gini-$M_{20}$) within our redshift range. We find the mass-size relation to be nearly flat at redshift z=6, and consistently increases as redshift decreases, and find the galaxy classification methods have minimal correlation with each other in our redshift range. We also investigate the impact that different stages of our imaging pipeline have on these morphological measures to determine how robust mock catalogs are to different choices at each step. Finally, we test the addition of incorporating light from AGNs into our pipeline and find that while the population of galaxies that have significant AGN luminosity is low, those galaxies do tend to have higher Sersic indices once the AGN luminosity is added, rectifying some of the systematic bias towards lower Sersic indices present in our dataset., Comment: 17 pages, 14 figures

We study dual AGN host galaxy morphologies at $z=2$ using the ASTRID simulation, selecting black hole (BH) pairs with small separation ($\Delta r<30\rm{kpc}$), high mass ($M_{\text{BH,12}}>10^7M_\odot$), and luminosity ($L_{\text{bol,12}}>10^{43}\rm{erg/s}$). We kinematically decompose (using MORDOR) $\sim1000$ dual AGN hosts into standard components - a `disk' (thin and thick disk, pseudo-bulge) and 'bulge' (bulge and halo) and define disk-dominated galaxies by the disk-to-total $D/T\geq0.5$. In ASTRID, $60.9\pm2.1\%$ of dual AGN hosts (independent of separation) are disk-dominated, with the $D/T$ distribution peaking at $\sim0.7$. Notably, hosts of BH pairs have similar morphologies (most either both disk or bulge-dominated). In dual-AGN hosts, the $D/T$ increases from $\sim17\% $ at $M_{\rm *}\sim 10^{9} M_{\odot}$ to $ 64\% $ for $M_{\rm *} \sim 10^{11.5} M_{\odot}$, and the pseudo-bulge is the dominant component of the disk fraction at the high mass end. Moreover, dual AGN hosts exhibit a higher fraction of disk/large pseudo-bulge than single-AGN hosts. The Disk-to-Total ratio is approximately constant with BH mass or AGN luminosity. We also create mock images of dual AGN host galaxies, employing morphological fitting software Statmorph to calculate morphological parameters and compare them with our kinematic decomposition results. Around $83.3\pm2.4\%$ of galaxies display disk-like profiles, of which $\sim60.7\pm2.2\%$ are kinematically confirmed as disks. Se\'rsic indices and half-mass radii of dual AGN host galaxies align with observational measurements from HST at $z\sim2$. Around $34\%$ are identified as mergers from the $\text{Gini}-M_{20}$ relation. We find two dual AGN hosted by galaxies that exhibit disk-like se\'rsic index $n_{12}<1$ and $(D/T)_{12}>0.5$, which are in remarkable agreement with properties of recently discovered dual quasars in disk galaxies at $z\sim 2$., Comment: 17 pages, 14 figures, submitted to the Open Journal of Astrophysics

The elusive massive black hole (MBH) seeds stand to be revealed by the Laser Space Antenna Interferometer through mergers. As an aftermath of galaxy mergers, MBH coalescence is a vastly multi-scale process connected to galaxy formation. We introduce the "Massive black hole Assembly in Galaxies Informed by Cosmological Simulations" (MAGICS) suite, with galaxy/MBH properties and orbits recovered from large-volume cosmological simulation ASTRID. The simulations include subgrid star formation, supernovae feedback, and MBH accretion/feedback. In this first suite, we extract fifteen representative galaxy mergers with seed MBHs to examine their dynamics at an improved mass and spatial resolution (by $\sim2000$ and $\sim20$) and follow MBH orbits down to $\sim10\,\text{pc}$. We find that the seed MBH energy loss and orbital decay are largely governed by global torques induced by the galaxy merger process on scales resolvable by cosmological simulations. Specifically, pairs sink quickly if their orbits shrink rapidly below $1\,\text{kpc}$ during the first $\sim200\,\text{Myr}$ of pairing due to effective energy loss in major galaxy mergers, whereas MBHs gaining energy in minor galaxy mergers with head-on collisions are likely to stall. High initial eccentricities ($e_\text{init}>0.5$) and high stellar densities at kpc scales ($\rho_\text{star}>0.05\,M_\odot/\text{pc}^3$) also lead to most efficient decays. $\sim50\%$ high-redshift seed MBH pairs experience consecutive galaxy mergers and are more likely to stall at $\sim1\,\text{kpc}$. For a subset of systems, we carry out N-Body re-simulations until binary formation and find that some stalled systems merge at high-z when embedded in sufficient nuclear star clusters., Comment: 21 pages, 20 Figures. Submitted to The Open Journal of Astrophysics. Comments welcome!

We present an analysis searching for dual AGN among 62 high-redshift ($2.5 < z < 3.5$) X-ray sources selected from publicly available deep Chandra fields. We aim to quantify the frequency of dual AGN in the high-redshift Universe, which holds implications for black hole merger timescales and low-frequency gravitational wave detection rates. We analyze each X-ray source using BAYMAX, an analysis tool that calculates the Bayes factor for whether a given archival Chandra AGN is more likely a single or dual point source. We find no strong evidence for dual AGN in any individual source in our sample. We then increase our sensitivity to search for dual AGN across the sample by comparing our measured distribution of Bayes factors to that expected from a sample composed entirely of single point sources, and again find no evidence for dual AGN in the observed sample distribution. Although our analysis utilizes one of the largest Chandra catalogs of high-$z$ X-ray point sources available to study, the findings remain limited by the modest number of sources observed at the highest spatial resolution with Chandra and the typical count rates of the detected sources. Our non-detection allows us to place an upper-limit on the X-ray dual AGN fraction between $2.5

We present an analysis showcasing how the Advanced X-ray Imaging Satellite (AXIS), a proposed NASA Probe-class mission, will significantly increase our understanding of supermassive black holes undergoing mergers -- from kpc to sub-pc scales. In particular, the AXIS point spread function, field of view, and effective area are expected to result in (1) the detection of hundreds to thousands of new dual AGNs across the redshift range 0 < z < 5 and (2) blind searches for binary AGNs that are exhibiting merger signatures in their light curves and spectra. AXIS will detect some of the highest-redshift dual AGNs to date, over a large range of physical separations. The large sample of AGN pairs detected by AXIS (over a magnitude more than currently known) will result in the first X-ray study that quantifies the frequency of dual AGNs as a function of redshift up to z = 4., Comment: 22 pages, 5 figures. Accepted to Universe (05/24/2024)

We look for simulated star-forming linear wakes such as the one recently discovered by van Dokkum et al. (2023) in the cosmological hydrodynamical simulation ASTRID. Amongst the runaway black holes in ASTRID, none are able to produce clear star-forming wakes. Meanwhile, fly-by encounters, typically involving a compact galaxy (with a central black hole) and a star-forming galaxy (with a duo of black holes) reproduce remarkably well many of the key properties (its length and linearity; recent star formation, etc.) of the observed star-forming linear feature. We predict the feature to persist for approximately 100 Myr in such a system and hence constitute a rare event. The feature contains a partly stripped galaxy (with $M_{\rm gal}=10^9 \sim 10^{10}M_\odot$) and a dual BH system ($M_{\rm BH}=10^5 \sim 10^7\,M_\odot$) in its brightest knot. X-ray emission from AGN in the knot should be detectable in such systems. After $100\sim 200\,{\rm Myrs}$ from the first fly-by, the galaxies merge leaving behind a triple black hole system in a (still) actively star-forming early-type remnant of mass $\sim 5\times 10^{10}\,M_\odot$. Follow-up JWST observations may be key for revealing the nature of these linear features by potentially detecting the older stellar populations constituting the bright knot. Confirmation of such detections may therefore help discriminate a fly-by encounter from a massive BH wake to reveal the origin of such features., Comment: 8 pages, 5 figures, comments welcome

We present the PRIYA suite of cosmological simulations, based on the code and hydrodynamic model of the ASTRID simulation, and designed for cosmological analyses of the Lyman-$\alpha$ forest. Our simulation suite spans a $9$-dimensional parameter space, including $4$ cosmological parameters and $5$ astrophysical/thermal parameters. We have run $48$ low fidelity simulations with $1536^3$ particles in a $120$ Mpc/h box and $3$ high fidelity simulations with $3072^3$ particles in a $120$ Mpc/h box. All our simulations include a full physics model for galaxy formation, including supernova and AGN feedback, and thus also contain a realistic population of DLAs. We advance on earlier simulations suites by larger particle loads, by incorporating new physical models for patchy hydrogen and helium reionization, and by self-consistently incorporating a model for AGN feedback. We show that patchy helium reionization imprints an excess in the 1D flux power spectrum on large scales, which may allow future measurements of helium reionization bubble sizes. Simulation parameters are chosen based on a Latin hypercube design and a Gaussian process is used to interpolate to arbitrary parameter combinations. We build a multi-fidelity emulator for the 1D flux power spectrum and the mean IGM temperature. We show that our final interpolation error is $< 1\%$ and that our simulations produce a flux power spectrum converged at the percent level for $z=5.4$ - $2.2$. Our simulation suite will be used to interpret Lyman-$\alpha$ forest 1D flux power spectra from SDSS and future DESI data releases., Comment: 25 pages, 11 figures, matches version published in JCAP

We present CAMELS-ASTRID, the third suite of hydrodynamical simulations in the Cosmology and Astrophysics with MachinE Learning (CAMELS) project, along with new simulation sets that extend the model parameter space based on the previous frameworks of CAMELS-TNG and CAMELS-SIMBA, to provide broader training sets and testing grounds for machine-learning algorithms designed for cosmological studies. CAMELS-ASTRID employs the galaxy formation model following the ASTRID simulation and contains 2,124 hydrodynamic simulation runs that vary 3 cosmological parameters ($\Omega_m$, $\sigma_8$, $\Omega_b$) and 4 parameters controlling stellar and AGN feedback. Compared to the existing TNG and SIMBA simulation suites in CAMELS, the fiducial model of ASTRID features the mildest AGN feedback and predicts the least baryonic effect on the matter power spectrum. The training set of ASTRID covers a broader variation in the galaxy populations and the baryonic impact on the matter power spectrum compared to its TNG and SIMBA counterparts, which can make machine-learning models trained on the ASTRID suite exhibit better extrapolation performance when tested on other hydrodynamic simulation sets. We also introduce extension simulation sets in CAMELS that widely explore 28 parameters in the TNG and SIMBA models, demonstrating the enormity of the overall galaxy formation model parameter space and the complex non-linear interplay between cosmology and astrophysical processes. With the new simulation suites, we show that building robust machine-learning models favors training and testing on the largest possible diversity of galaxy formation models. We also demonstrate that it is possible to train accurate neural networks to infer cosmological parameters using the high-dimensional TNG-SB28 simulation set.

We use the ASTRID cosmological hydrodynamic simulation to investigate the properties and evolution of triple and quadruple Massive Black Hole (MBH) systems at $z = 2-3$. Only a handful of MBH tuple systems have been detected to date. In ASTRID, we find $4\%$ of the $M_{\rm BH}>10^7\,M_\odot$ are in tuples with $\Delta r_{\rm max} < 200\,{\rm kpc}$. The tuple systems span a range of separations with the majority of the observable AGN systems at $\Delta r \sim 50-100$ kpc. They include some of the most massive BHs (up to $10^{10} \,M_\odot$) but with at least one of the components of $M_{\rm BH} \sim 10^7 \,M_\odot$. Tuples' host galaxies are typically massive with $M_* \sim 10^{10-11} \,M_\odot$. We find that $>10\%$ massive halos with $M_{\rm halo} > 10^{13} M_\odot$ host MBH tuples. Following the subsequent interactions between MBHs in tuples, we found that in $\sim 5\%$ of the triplets all three MBHs merge within a Gyr, and $15\%$ go through one merger. As a by-product of the complex multi-galaxy interaction of these systems, we also find that up to $\sim 5\%$ of tuples lead to runaway MBHs. In ASTRID, virtually all of the ultramassive black holes ($>10^{10} \,M_\odot $) have undergone a triple quasar phase while for BHs with $M_{\rm BH} \sim 10^9 \,M_\odot$ this fraction drops to $50\%$., Comment: 10 pages, 9 figures; comments welcome

In the near future, projects like LISA and Pulsar Timing Arrays are expected to detect gravitational waves from mergers between supermassive black holes, and it is crucial to precisely model the underlying merger populations now to maximize what we can learn from this new data. Here we characterize expected high-redshift (z > 2) black hole mergers using the very large volume Astrid cosmological simulation, which uses a range of seed masses to probe down to low-mass BHs, and directly incorporates dynamical friction so as to accurately model the dynamical processes which bring black holes to the galaxy center where binary formation and coalescence will occur. The black hole populations in Astrid include black holes down to 10$^{4.5}$ M$_\odot$, and remain broadly consistent with the TNG simulations at scales > 10$^6$ M$_\odot$ (the seed mass used in TNG). By resolving lower-mass black holes, the overall merger rate is ~5x higher than in TNG. However, incorporating dynamical friction delays mergers compared to a recentering scheme, reducing the high-z merger rate mass-matched mergers by a factor of ~2x. We also calculate the expected LISA Signal-to-Noise values, and show that the distribution peaks at high SNR (>100), emphasizing the importance of implementing a seed mass well below LISA's peak sensitivity (10$^6$ M$_\odot$) to resolve the majority of LISA's GW detections., Comment: 11 pages, 10 figures, submitted to MNRAS

Detection of massive binary black hole (BBH) mergers at high redshifts is a target for LISA space mission. While the individual masses of a BBH merger are redshifted, the mass ratio of BBH mergers is independent of their redshift. Therefore, if there is an independent correlation between the mass ratio and redshift, such a relationship can be used to i) infer the redshift of the merging binaries, and together with the luminosity distance measurement ($D_L$), constrain the expansion rate of the universe at high redshifts $H(z)$, and ii) constrain models of supermassive black hole seed formation in the universe assuming a fixed cosmology. We discuss why there is an expected relation between the mass ratio of the massive BBHs with their redshift and show the clues for this relation by analyzing cosmological hydrodynamical simulations of BBH mergers. This approach opens up the possibility of directly measuring the expansion rate at redshift $z \approx 2$ with LISA for the first time. Moreover, we discover a trend between seed mass and mass ratio of massive BBHs which by itself is a major result that could be exploited to constrain the formation scenarios of supermassive BH seeds., Comment: to be submitted to ApJL. Comments are welcome

Intermediate-Mass Black Holes (IMBHs) of $10^3-10^6 \, M_\odot$ are commonly found at the center of dwarf galaxies. Simulations and observations convincingly show that a sizable population of IMBHs could wander off-center in galaxies. We use the cosmological simulation ASTRID to study the orbital and radiative properties of wandering IMBHs in massive galaxies at $z\sim3$. We find that this population of black holes has large orbital inclinations ($60^\circ\pm22^\circ$) with respect to the principal plane of the host. The eccentricity of their orbits is also significant ($0.6\pm0.2$) and decreases with time. Wandering IMBHs undergo spikes of accretion activity around the pericenter of their orbits, with rates $10^{-3}-10^{-5}$ times the Eddington rate and a median accretion duty cycle of $\sim 12\%$. Their typical spectral energy distribution peaks in the infrared at $\sim 11 \, \mu \rm m$ rest-frame. Assuming a standard value of $10\%$ for the matter-to-energy radiative efficiency, IMBHs reach $2-10$ keV X-ray luminosities $>10^{37} \, \mathrm{erg\,s^{-1}}$ for $\sim10\%$ of the time. This luminosity corresponds to fluxes $>10^{-15} \, \mathrm{erg \, s^{-1} \, cm^{-2}}$ within $10$ Mpc. They could be challenging to detect because of competing emissions from X-ray binaries and the interstellar medium. X-ray luminosities $> 10^{41} \, \mathrm{erg \, s^{-1}}$, in the hyper-luminous X-ray sources (HLXs) regime, are reached by $\sim 7\%$ of the IMBHs. These findings suggest that HLXs are a small subset of the wandering IMBH population, which is characterized by luminosities $10^3-10^4$ times fainter. Dedicated surveys are needed to assess the demographics of this missing population of black holes., Comment: Accepted for publication in MNRAS. This is the final version of the manuscript. 9 pages, 7 figures

Massive black holes in the centers of galaxies today must have grown by several orders of magnitude from seed black holes formed at early times. Detecting a population of intermediate mass black holes (IMBHs) can provide constraints on these elusive BH seeds. Here we use the large volume, cosmological hydrodynamical simulation Astrid, which includes IMBH seeds and dynamical friction to investigate the population of IMBH seeds. Dynamical friction is largely inefficient at sinking and merging seed IMBHs at high-z. This leads to an extensive population (several hundred per galaxy) of wandering IMBHs in large halos at z~2. A small fraction of these IMBHs are detectable as HLXs, Hyper Luminous X-ray sources. Importantly, at z ~ 2, IMBHs mergers produce the peak of GW events. We find close to a million GW events in Astrid between z=2-3 involving seed IMBH mergers. These GW events (almost all detectable by LISA) at cosmic noon should provide strong constraints on IMBH seed models and their formation mechanisms. At the center of massive galaxies, where the number of IMBHs can be as high as 10-100, SMBH-IMBH pairs can form. These Intermediate mass ratio inspirals (IMRIs) and extreme mass ratio inspirals (EMRIs), will require the next generation of milli-muHz space-based GW interferometers to be detected. Large populations of IMBHs around massive black holes will probe their environments and MBH causal structure.

Galaxy mergers produce pairs of supermassive black holes (SMBHs), which may be witnessed as dual quasars if both SMBHs are rapidly accreting. The kiloparsec (kpc)-scale separation represents a physical regime sufficiently close for merger-induced effects to be important yet wide enough to be directly resolvable with the facilities currently available. Whereas many kpc-scale dual active galactic nuclei--the low-luminosity counterparts of quasars--have been observed in low-redshift mergers, no unambiguous dual quasar is known at cosmic noon ($z\approx2$), the peak of global star formation and quasar activity. Here we report multiwavelength observations of SDSS J0749+2255 as a kpc-scale, dual-quasar system hosted by a galaxy merger at cosmic noon ($z=2.17$). We discover extended host galaxies associated with the much brighter compact quasar nuclei (separated by 0.46" or 3.8 kpc) and low-surface-brightness tidal features as evidence for galactic interactions. Unlike its low-redshift and low-luminosity counterparts, SDSS J0749+2255 is hosted by massive compact disk-dominated galaxies. The apparent lack of stellar bulges and the fact that SDSS J0749+2255 already follows the local SMBH mass--host stellar mass relation, suggest that at least some SMBHs may have formed before their host stellar bulges. While still at kpc-scale separations where the host-galaxy gravitational potential dominates, the two SMBHs may evolve into a gravitationally bound binary system in around 0.22 Gyr., Comment: 88 pages, 18 figures, 8 tables; published

The origin of rare and elusive ultramassive black holes (UMBH, with MBH > 1e10 Msun) is an open question. Using the large volume cosmological hydrodynamic simulation ASTRID, we report on the formation of an extremely massive UMBH with MBH ~ 1e11 Msun at z~2. The UMBH is assembled as a result of two successive mergers of massive galaxies each with stellar mass M* > 3e11 Msun that also produces a bright, rare triple quasar system powered by three ~10^9 Msun black holes. The second merger of supermassive black holes (SMBHs) follows the first after 150 Myrs. The merger events lead to sustained Eddington accretion onto the central SMBH, forming an UMBH in the center of a massive compact stellar core with M* > 2e12 Msun. The strong feedback of the UMBH quenches the surrounding star formation to < 10 Msun/yr in the inner 50 kpc/h region. There are two more UMBHs with MBH > 5e10 Msun at z>2 in ASTRID which are also produced by major mergers of galaxies, and their progenitors can be observed as quasar triplets of lower luminosity. The rarely observed quasar multiples can be the cradle of UMBHs at high redshift, and likely end up in the center of the most massive clusters., Comment: 7 pages, 4 figures, comments welcome

The statistics of galactic-scale quasar pairs can elucidate our understanding of the dynamical evolution of supermassive black hole (SMBH) pairs, the duty cycles of quasar activity in mergers, or even the nature of dark matter, but have been challenging to measure at cosmic noon, the prime epoch of massive galaxy and SMBH formation. Here we measure a double quasar fraction of $\sim 6.2\pm0.5\times 10^{-4}$ integrated over $\sim 0.3-3$ arcsec separations (projected physical separations of $\sim 3-30\,{\rm kpc}$ at $z\sim 2$) in luminous ($L_{\rm bol}>10^{45.8}\,{\rm erg\,s^{-1}}$) unobscured quasars at $1.5

We examine the dual (both BHs active) and offset (one BH active) AGN population (comprising $\sim$ 2000 pairs at $0.5\,\text{kpc}\lesssim\Delta r<30\,\text{kpc}$) at $z=2\sim3$ in the ASTRID simulation covering (360 cMpc)${^3}$. The dual (offset) AGN make up $3.0(2.2)\%$ of all AGN at $z=2$. The dual fraction is roughly constant while the offset fraction increases by a factor of ten from $z=4\sim2$. Compared with the full AGN population, duals are characterized by a low $M_\text{BH}/M_*$ ratio, a high specific star-formation rate (sSFR) of $\sim 1\,\text{Gyr}^{-1}$, and a high Eddington ratio ($\sim 0.05$, double that of single AGN). The dual AGN are formed in major galaxy mergers (typically involving $M_\text{halo}<10^{13}\,M_\odot$), with BHs that have similar masses. At small separations (when their host galaxies are in the late phase of the merger) duals become $2\sim8$ times brighter (albeit more obscured) than at larger separations. $80\%$ of these bright, close duals merge in the simulation within $\sim500\,\text{Myrs}$. Notably, the initially less-massive BH in duals frequently becomes the brighter AGN during the galaxy merger. In offset AGN, the active BH is typically $\gtrsim 10$ times more massive than its non-active counterpart and than most BHs in duals. Offsets are predominantly formed in minor galaxy mergers with the active BH residing in the center of massive halos ($ M_\text{ halo}\sim 10^{13-14}\,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\sim5\,\text{kpc}$ for a few hundred Myrs., Comment: 17 Pages, 18 Figures; comments welcome!

Using the novel semi-numerical code for reionization AMBER, we model the patchy kinetic Sunyaev-Zel'dovich (kSZ) effect by directly specifying the reionization history with the redshift midpoint $z_\mathrm{mid}$, duration $\Delta_\mathrm{z}$, and asymmetry $A_\mathrm{z}$. We further control the ionizing sources and radiation through the minimum halo mass $M_\mathrm{h}$ and the radiation mean free path $\lambda_\mathrm{mfp}$. AMBER reproduces the free electron number density and the patchy kSZ power spectrum of radiation-hydrodynamic simulations at the target resolution ($1\,{\rm Mpc}/h$) with matched reionization parameters. With a suite of $(2\,{\rm Gpc}/h)^3$ simulations using AMBER, we first constrain the redshift midpoint $6.08$ is needed to ensure an end of reionization before $z=5.5$). We show that constraints on the asymmetry require $\sim 0.1\,\mu k^2$ measurement accuracy at multipoles other than $\ell=3000$. Finally, we find that the amplitude and shape of the kSZ spectrum are only weakly sensitive to $M_h$ under a fixed reionization history and radiation mean-free path., Comment: 20 pages, 17 Figures; comments welcome

We examine massive black hole (MBH) mergers and their associated gravitational wave signals from the large-volume cosmological simulation Astrid. Astrid includes galaxy formation and black hole models recently updated with a MBH seed population between $3\times 10^4M_{\odot}/h$ and $3\times 10^5M_{\odot}/h$ and a sub-grid dynamical friction (DF) model to follow the MBH dynamics down to $1.5\;\text{ckpc}/h$. We calculate initial eccentricities of MBH orbits directly from the simulation at kpc-scales, and find orbital eccentricities above $0.7$ for most MBH pairs before the numerical merger. After approximating unresolved evolution on scales below ${\sim 200\,\text{pc}}$, we find that the in-simulation DF on large scales accounts for more than half of the total orbital decay time ($\sim 500\,\text{Myrs}$) due to DF. The binary hardening time is an order of magnitude longer than the DF time, especially for the seed-mass binaries ($M_\text{BH}<2M_\text{seed}$). As a result, only $\lesssim20\%$ of seed MBH pairs merge at $z>3$ after considering both unresolved DF evolution and binary hardening. These $z>3$ seed-mass mergers are hosted in a biased population of galaxies with the highest stellar masses of $>10^9\,M_\odot$. With the higher initial eccentricity prediction from Astrid, we estimate an expected merger rate of $0.3-0.7$ per year from the $z>3$ MBH population. This is a factor of $\sim 7$ higher than the prediction using the circular orbit assumption. The LISA events are expected at a similar rate, and comprise $\gtrsim 60\%$ seed-seed mergers, $\sim 30\%$ involving only one seed-mass MBH, and $\sim 10\%$ mergers of non-seed MBHs., Comment: 17 pages, 13 Figures; comments are welcome

We introduce the ASTRID simulation, a large-scale cosmological hydrodynamic simulation in a $250$ Mpc/h box with $2\times 5500^3$ particles. ASTRID contains a large number of high redshift galaxies, which can be compared to future survey data, and resolves galaxies in halos more massive than $2\times 10^9 M_\odot$. ASTRID has been run from $z=99$ to $z=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 summarise the implemented models, and the technical choices we took when developing the simulation code. We validate the model, showing good agreement with observed UV luminosity functions, galaxy stellar mass functions and specific star formation rates. 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$, halos with $M \sim 2\times 10^9 M_\odot$ which have been reionized have a star formation rate $1.5$ times greater than those which have not yet been reionized., Comment: 15 pages, 12 figures. Accepted by MNRAS, added star formation rate density

We present the evolution of black holes (BHs) and their relationship with their host galaxies in Astrid, a large-volume cosmological hydrodynamical simulation with box size 250 $h^{-1} \rm Mpc$ containing $2\times5500^3$ particles evolved to z=3. Astrid statistically models BH gas accretion and AGN feedback to their environments, applies a power-law distribution for BH seed mass $M_{\rm sd}$, uses a dynamical friction model for BH dynamics and executes a physical treatment of BH mergers. The BH population is broadly consistent with empirical constraints on the BH mass function, the bright end of the luminosity functions, and the time evolution of BH mass and accretion rate density. The BH mass and accretion exhibit a tight correlation with host stellar mass and star formation rate. We trace BHs seeded before z>10 down to z=3, finding that BHs carry virtually no imprint of the initial $M_{\rm sd}$ except those with the smallest $M_{\rm sd}$, where less than 50\% of them have doubled in mass. Gas accretion is the dominant channel for BH growth compared to BH mergers. With dynamical friction, Astrid predicts a significant delay for BH mergers after the first encounter of a BH pair, with a typical elapse time of about 200 Myrs. There are in total $4.5 \times 10^5$ BH mergers in Astrid at z>3, $\sim 10^3$ of which have X-ray detectable EM counterparts: a bright kpc scale dual AGN with $L_X>10^{43}$ erg/s. BHs with $M_{\rm BH} \sim 10^{7-8} M_{\odot}$ experience the most frequent mergers. Galaxies that host BH mergers are unbiased tracers of the overall $M_{\rm BH} - M_{*}$ relation. Massive ($>10^{11} M_{\odot}$) galaxies have a high occupation number (>10) of BHs, and hence host the majority of BH mergers., Comment: Match the published version

The Abundance Matching Box for the Epoch of Reionization (AMBER) is a semi-numerical code for modeling the cosmic dawn. The new algorithm is not based on the excursion set formalism for reionization, but takes the novel approach of calculating the reionization-redshift field $z_\mathrm{re}(\boldsymbol{x})$ assuming that hydrogen gas encountering higher radiation intensity are photoionized earlier. Redshift values are assigned while matching the abundance of ionized mass according to a given mass-weighted ionization fraction $\bar{x}_\mathrm{i}(z)$. The code has the unique advantage of allowing users to directly specify the reionization history through the redshift midpoint $z_\mathrm{mid}$, duration $\Delta_\mathrm{z}$, and asymmetry $A_\mathrm{z}$ input parameters. The reionization process is further controlled through the minimum halo mass $M_\mathrm{min}$ for galaxy formation and the radiation mean free path $l_\mathrm{mfp}$ for radiative transfer. We implement improved methods for constructing density, velocity, halo, and radiation fields, which are essential components for modeling reionization observables. We compare AMBER with two other semi-numerical methods and find that our code more accurately reproduces the results from radiation-hydrodynamic simulations. The parallelized code is over four orders of magnitude faster than radiative transfer simulations and will efficiently enable large-volume models, full-sky mock observations, and parameter-space studies. AMBER will be made publicly available to facilitate and transform studies of the EoR., Comment: AMBER is publicly available at https://github.com/hytrac/amber

In this work we establish and test methods for implementing dynamical friction for massive black hole pairs that form in large volume cosmological hydrodynamical simulations which 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 dynamical friction (DF) from collisionless particles, black holes can effectively sink close to the galaxy center, 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\%$ of DF from collisionless particles for BH masses $> 10^{7}$ M$_{\odot}$. Using our best DF model, we calculate the merger rate down to $z=1.1$ using an $L_{\rm box}=35$ Mpc$/h$ simulation box. We predict $\sim 2$ mergers per year for $z>1.1$ peaking at $z\sim 2$. These merger rates are within the range obtained in previous work using similar-resolution hydro-dynamical simulations. We show that the rate is enhanced by factor of $\sim 2$ when DF is taken into account in the simulations compared to the no-DF run. This is due to $>40\%$ more black holes reaching the center of their host halo when DF is added., Comment: 17 pages, 12 figures; Comments welcome!

Covariance matrices are among the most difficult pieces of end-to-end cosmological analyses. In principle, for two-point functions, each component involves a four-point function, and the resulting covariance often has hundreds of thousands of elements. We investigate various compression mechanisms capable of vastly reducing the size of the covariance matrix in the context of cosmic shear statistics. This helps identify which of its parts are most crucial to parameter estimation. We start with simple compression methods, by isolating and "removing" 200 modes associated with the lowest eigenvalues, then those with the lowest signal-to-noise ratio, before moving on to more sophisticated schemes like compression at the tomographic level and, finally, with the Massively Optimized Parameter Estimation and Data compression (MOPED). We find that, while most of these approaches prove useful for a few parameters of interest, like $\Omega_m$, the simplest yield a loss of constraining power on the intrinsic alignment (IA) parameters as well as $S_8$. For the case considered -- cosmic shear from the first year of data from the Dark Energy Survey -- only MOPED was able to replicate the original constraints in the 16-parameter space. Finally, we apply a tolerance test to the elements of the compressed covariance matrix obtained with MOPED and confirm that the IA parameter $A_{\mathrm{IA}}$ is the most susceptible to inaccuracies in the covariance matrix., Comment: 13 pages, 12 figures, 1 table. V2 reflects version accepted for publication in PRD

In the Simulations and Constructions of the Reionization of Cosmic Hydrogen (SCORCH) project, we compare analytical models of the hydrogen ionization fraction with radiation-hydrodynamic simulations. We derive analytical models of the mass-weighted hydrogen ionization fraction from the local ionization balance equations as a more accurate alternative to the widely adopted model based on the volume filling factor. In particular, our model has a recombination term quadratic in the ionization fraction, which is consistent with the two-body interaction nature of recombination. Then, we use the radiation-hydrodynamic simulations to study the clumping factors needed to solve the analytical equations, and provide accurate fitting functions. We find that the ionized hydrogen clumping factors from our radiative transfer simulations are significantly different than those from other simulations that use a uniform photoionization background. In addition to redshift dependence, we also see the dependence of ionized hydrogen clumping factor on ionization fraction, and we incorporate this into our fits. We calculate the reionization histories using our analytical models and clumping factors and compare with widely adopted models, and all of our models achieve $<7\%$ difference from simulation results while the other models have $>20\%$ deviations. The Thomson optical depths from reionization calculated from our analytical models result in $<5\%$ deviation from simulations, while the previous analytical models have $>20\%$ difference in and could result in biased conclusions of the IGM reionization., Comment: 14 pages, 8 figures