Showing 1–50 of 175 results for author: Di Matteo, T

Active Galactic Nuclei and Host Galaxies in COSMOS-Web. II. First Look at the Kpc-scale Dual and Offset AGN Population

Authors: Junyao Li, Ming-Yang Zhuang, Yue Shen, Marta Volonteri, Nianyi Chen, Tiziana Di Matteo

Abstract: 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 systematicall… ▽ More 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. △ Less

Submitted 23 May, 2024; originally announced May 2024.

Comments: 20 pages, 8 figures, submitted to ApJ

Astrometric Jitter as a Detection Diagnostic for Recoiling and Slingshot Supermassive Black Hole Candidates

Authors: Anavi Uppal, Charlotte Ward, Suvi Gezari, Priyamvada Natarajan, Nianyi Chen, Patrick LaChance, Tiziana Di Matteo

Abstract: 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… ▽ More 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 3pi 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. △ Less

Submitted 17 May, 2024; originally announced May 2024.

The evolution of galaxy morphology from redshift z=6 to 3: Mock JWST observations of galaxies in the ASTRID simulation

Authors: Patrick LaChance, Rupert Croft, Yueying Ni, Nianyi Chen, Tiziana Di Matteo, Simeon Bird

Abstract: We present mock JWST observations for more than 215,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… ▽ More We present mock JWST observations for more than 215,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. △ Less

Submitted 29 January, 2024; originally announced January 2024.

Comments: 17 pages, 14 figures

Simulated host galaxy analogs of high-z quasars observed with JWST

Authors: Sabrina Berger, Madeline A. Marshall, J. Stuart B. Wyithe, Tiziana di Matteo, Yueying Ni, Stephen M. Wilkins

Abstract: The hosts of two low-luminosity high-z quasars, J2255+0251 and J2236+0032, were recently detected using JWST's NIRCam instrument. These represent the first high-z quasar host galaxy stellar detections and open a new window into studying high-z quasars. We examine the implications of the measured properties of J2255+0251 and J2236+0032 within the context of the hydrodynamic simulation BlueTides at… ▽ More The hosts of two low-luminosity high-z quasars, J2255+0251 and J2236+0032, were recently detected using JWST's NIRCam instrument. These represent the first high-z quasar host galaxy stellar detections and open a new window into studying high-z quasars. We examine the implications of the measured properties of J2255+0251 and J2236+0032 within the context of the hydrodynamic simulation BlueTides at z = 6.5. We find that these observed quasars fall on the BlueTides stellar to black hole mass relation and have similar luminosities to the brightest simulated quasars. We predict their star formation rates, estimating approximately $10^{2-3}$ $M_{\odot}/ \rm yr$ for both quasar hosts. J2255+0251 and J2236+0032's host galaxy radii also fall within estimates of the radii of the simulated host galaxies of similar luminosity quasars. We generate mock JWST NIRCam images of analogs to the observed quasars within BlueTides and perform a point source removal to illustrate both a qualitative and quantitative comparison of the measured and simulated radii and magnitudes. The quasar subtraction works well for similar luminosity quasars, and the recovered host images are consistent with what was observed for J2255+0251 and J2236+0032, further supporting the success of those observations. We also use our mock imaging pipeline to make predictions for the detection of J2255+0251 and J2236+0032's hosts in upcoming JWST observations. We anticipate that the simulation analogs of future high-z quasar host discoveries will allow us to make accurate predictions of their properties beyond the capabilities of JWST. △ Less

Submitted 18 April, 2024; v1 submitted 16 January, 2024; originally announced January 2024.

Comments: Accepted to MNRAS. 15 pages, 11 figures

z~2 dual AGN host galaxies are disky: stellar kinematics in the ASTRID Simulation

Authors: Ekaterine Dadiani, Tiziana Di Matteo, Nianyi Chen, Patrick Lachance, Yue Shen, Yu-Ching Chen, Rupert Croft, Yueying Ni, Simeon Bird

Abstract: We study dual AGN host galaxy morphologies at $z=2$ using the ASTRID simulation, selecting black hole (BH) pairs with small separation ($Δ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 'b… ▽ More We study dual AGN host galaxy morphologies at $z=2$ using the ASTRID simulation, selecting black hole (BH) pairs with small separation ($Δ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ŕsic 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ŕsic 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$. △ Less

Submitted 21 December, 2023; originally announced December 2023.

Comments: 15 pages, 12 figures, submitted to the Open Journal of Astrophysics

MAGICS I. The First Few Orbits Encode the Fate of Seed Massive Black Hole Pairs

Authors: Nianyi Chen, Diptajyoti Mukherjee, Tiziana Di Matteo, Yueying Ni, Simeon Bird, Rupert Croft

Abstract: 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 fro… ▽ More 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 ($ρ_\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. △ Less

Submitted 25 April, 2024; v1 submitted 14 December, 2023; originally announced December 2023.

Comments: 21 pages, 20 Figures. Submitted to The Open Journal of Astrophysics. Comments welcome!

Searching for the Highest-z Dual AGN in the Deepest Chandra Surveys

Authors: Brandon Sandoval, Adi Foord, Steven W. Allen, Marta Volonteri, Aaron Stemo, Nianyi Chen, Tiziana Di Matteo, Kayhan Gultekin, Melanie Habouzit, Clara Puerto-Sanchez, Edmund Hodges-Kluck, Yohan Dubois

Abstract: 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 to… ▽ More 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<z<3.5$ of 4.8\%. Expanding substantially on these results at X-ray wavelengths will require future surveys spanning larger sky areas and extending to fainter fluxes than has been possible with Chandra. We illustrate the potential of the AXIS mission concept in this regard. △ Less

Submitted 4 December, 2023; originally announced December 2023.

Comments: 19 pages, 5 figures, 2 tales. Submitted to ApJ

Tracking Supermassive Black Hole Mergers from kpc to sub-pc Scales with AXIS

Authors: Adi Foord, Nico Cappelluti, Tingting Liu, Marta Volonteri, Melanie Habouzit, Fabio Pacucci, Stefano Marchesi, Nianyi Chen, Tiziana Di Matteo, Labani Mallick, Michael Koss

Abstract: 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… ▽ More 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. △ Less

Submitted 27 May, 2024; v1 submitted 13 November, 2023; originally announced November 2023.

Comments: 22 pages, 5 figures. Accepted to Universe (05/24/2024)

Modeling Quasar Proximity Zones in a Realistic Cosmological Environment with a Self-consistent Light Curve

Authors: Yihao Zhou, Huanqing Chen, Tiziana Di Matteo, Yueying Ni, Rupert A. C. Croft, Simeon Bird

Abstract: We study quasar proximity zones in a simulation that includes a self-consistent quasar formation model and realistic IGM environments. The quasar host halo is $10^{13}\ M_{\mathrm{\odot}}$ at $z=6$, more massive than typical halos studied in previous work. Between $6<z<7.5$, the quasar luminosity varies rapidly, with a mean magnitude of $M_{UV,mean}=-24.8$ and the fluctuation reaching up to two or… ▽ More We study quasar proximity zones in a simulation that includes a self-consistent quasar formation model and realistic IGM environments. The quasar host halo is $10^{13}\ M_{\mathrm{\odot}}$ at $z=6$, more massive than typical halos studied in previous work. Between $6<z<7.5$, the quasar luminosity varies rapidly, with a mean magnitude of $M_{UV,mean}=-24.8$ and the fluctuation reaching up to two orders of magnitude. Using this light curve to post-process the dense environment around the quasar, we find that the proximity zone size ($R_{p}$) ranges between $0.5-5$ pMpc. We show that the light curve variability causes a similar degree of scatter in $R_{p}$ as does the density fluctuation, both of which result in a standard deviation of $\sim 0.3$ pMpc). The $R_{p}$ traces the light curve fluctuations closely but with a time delay of $\sim 10^4\ \mathrm{yr}$, breaking the correspondence between the $R_{p}$ and the contemporaneous $M_{UV}$. This also indicates that we can only infer quasar activity within the past $\sim 10^4$ years instead of the integrated lifetime from $R_{p}$ in the later part of cosmic reionization. Compared with the variable light curve, a constant light curve underestimates the $R_{p}$ by 13% at the dim end ($M_{UV}\sim -23.5$), and overestimates the $R_{p}$ by 30% at the bright end ($M_{UV}\sim -26$). By calculating the $R_{p}$ generated by a number of quasars, we show that variable light curves predict a wider $R_{p}$ distribution than lightbulb models, and readily explain the extremely small $R_{p}$ values that have been observed. △ Less

Submitted 20 September, 2023; originally announced September 2023.

Comments: 15 pages, 12 figures, comments welcome

Fly-by galaxy encounters with multiple black holes produce star-forming linear wakes

Authors: Nianyi Chen, Patrick LaChance, Yueying Ni, Tiziana Di Matteo, Rupert Croft, Priyamvada Natarajan, Simeon Bird

Abstract: 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 blac… ▽ More 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. △ Less

Submitted 3 July, 2023; originally announced July 2023.

Comments: 8 pages, 5 figures, comments welcome

Efficient Reionization in a Large Hydrodynamic Galaxy Formation Simulation

Authors: James E. Davies, Simeon Bird, Simon Mutch, Yueying Ni, Yu Feng, Rupert Croft, Tiziana Di Matteo, J. Stuart B. Wyithe

Abstract: Accuracy in the topology and statistics of a simulated Epoch of Reionization (EoR) are vital to draw connections between observations and physical processes. While full radiative transfer models produce the most accurate reionization models, they are highly computationally expensive, and are infeasible for the largest cosmological simulations. Instead, large simulations often include EoR models th… ▽ More Accuracy in the topology and statistics of a simulated Epoch of Reionization (EoR) are vital to draw connections between observations and physical processes. While full radiative transfer models produce the most accurate reionization models, they are highly computationally expensive, and are infeasible for the largest cosmological simulations. Instead, large simulations often include EoR models that are pre-computed via the initial density field, or post-processed where feedback effects are ignored. We introduce Astrid-ES, a resimulation of the Astrid epoch of reionisation $20 > z > 5.5$ which includes an on-the-fly excursion-set reionization algorithm. Astrid-ES produces more accurate reionization histories without significantly impacting the computational time. This model directly utilises the star particles produced in the simulation to calculate the EoR history and includes a UV background which heats the gas particles after their reionization. We contrast the reionization topology and statistics in Astrid-ES with the previously employed parametric reionisation model, finding that in Astrid-ES, ionised regions are more correlated with galaxies, and the 21cm power-spectrum shows an increase in large scale power. We calculate the relation between the size of HII regions and the UV luminosity of the brightest galaxy within them. Prior to the overlap phase, we find a power-law fit of $\mathrm{log} (R) = -0.314 M_\mathrm{UV} - 2.550 \mathrm{log}(1+z) + 7.408$ with a standard deviation $σ_R < 0.15 \mathrm{dex}$ across all mass bins. We also examine the properties of halos throughout reionization, finding that while the properties of halos in the simulation are correlated with the redshift of reionisation, they are not greatly affected by reionisation itself. △ Less

Submitted 8 August, 2023; v1 submitted 13 June, 2023; originally announced June 2023.

Comments: 12 pages, 11 figures. Accepted for publication in Monthly Notices of the Royal Astronomical Society

PRIYA: A New Suite of Lyman-alpha Forest Simulations for Cosmology

Authors: Simeon Bird, Martin Fernandez, Ming-Feng Ho, Mahdi Qezlou, Reza Monadi, Yueying Ni, Nianyi Chen, Rupert Croft, Tiziana Di Matteo

Abstract: 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-$α$ 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… ▽ More 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-$α$ 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-$α$ forest 1D flux power spectra from SDSS and future DESI data releases. △ Less

Submitted 2 January, 2024; v1 submitted 8 June, 2023; originally announced June 2023.

Comments: 25 pages, 11 figures, matches version published in JCAP

Journal ref: JCAP 10 (2023) 037

AI-assisted super-resolution cosmological simulations III: Time evolution

Authors: Xiaowen Zhang, Patrick Lachance, Yueying Ni, Yin Li, Rupert A. C. Croft, Tiziana Di Matteo, Simeon Bird, Yu Feng

Abstract: In this work, we extend our recently developed super-resolution (SR) model for cosmological simulations to produce fully time consistent evolving representations of the particle phase-space distribution. We employ a style-based constrained generative adversarial network (Style-GAN) where the changing cosmic time is an input style parameter to the network. The matter power spectrum and halo mass fu… ▽ More In this work, we extend our recently developed super-resolution (SR) model for cosmological simulations to produce fully time consistent evolving representations of the particle phase-space distribution. We employ a style-based constrained generative adversarial network (Style-GAN) where the changing cosmic time is an input style parameter to the network. The matter power spectrum and halo mass function agree well with results from high-resolution N-body simulations over the full trained redshift range ($10 \le z \le 0$). Furthermore, we assess the temporal consistency of our SR model by constructing halo merger trees. We examine progenitors, descendants and mass growth along the tree branches. All statistical indicators demonstrate the ability of our SR model to generate satisfactory high-resolution simulations based on low-resolution inputs. △ Less

Submitted 20 May, 2023; originally announced May 2023.

Comments: 12 pages, 11 figures, code and movie available in https://github.com/sagasv5-xw/map2map on styled srsgan branch

Over-massive Central Black Holes in the Cosmological Simulations ASTRID and Illustris TNG50

Authors: Emma Jane Weller, Fabio Pacucci, Priyamvada Natarajan, Tiziana Di Matteo

Abstract: Recent dynamical measurements indicate the presence of a central SMBH with mass $\sim 3\times 10^6 \, \rm M_\odot$ in the dwarf galaxy Leo I, placing the system $\sim50$ times above the standard, local $M_{BH} - M_\star$ relation. While a few over-massive central SMBHs are reported in nearby isolated galaxies, this is the first detected in a Milky Way satellite. We used the ASTRID and Illustris TN… ▽ More Recent dynamical measurements indicate the presence of a central SMBH with mass $\sim 3\times 10^6 \, \rm M_\odot$ in the dwarf galaxy Leo I, placing the system $\sim50$ times above the standard, local $M_{BH} - M_\star$ relation. While a few over-massive central SMBHs are reported in nearby isolated galaxies, this is the first detected in a Milky Way satellite. We used the ASTRID and Illustris TNG50 LCDM cosmological simulations to investigate the assembly history of galaxies hosting over-massive SMBHs. We estimate that, at the stellar mass of Leo I, $\sim15\%$ of galaxies above the $M_{BH} - M_\star$ relation lie $>10$ times above it. Leo I-like systems are rare but exist in LCDM simulations: they occur in $\sim0.005\%$ of all over-massive systems. Examining the properties of simulated galaxies harboring over-massive central SMBHs, we find that: (i) stars assemble more slowly in galaxies above the $M_{BH} - M_\star$ relation; (ii) the gas fraction in these galaxies experiences a significantly steeper decline over time; and (iii) $>95\%$ of satellite host galaxies in over-dense regions are located above the $M_{BH} - M_\star$ relation. This suggests that massive satellite infall and consequent tidal stripping in a group/dense environment can drive systems away from the $M_{BH} - M_\star$ relation, causing them to become over-massive. As the merging histories of over-massive and under-massive systems do not differ, we conclude that additional environmental effects, such as being in overdense regions, must play a crucial role. In the high-$z$ Universe, central over-massive SMBHs are a signature of heavy black hole seeds; we demonstrate, in contrast, that low-$z$ over-massive systems result from complex environmental interactions. △ Less

Submitted 3 May, 2023; originally announced May 2023.

Comments: Accepted for publication in MNRAS. 9 pages, 8 figures

Journal ref: Monthly Notices of the Royal Astronomical Society, 2023, Volume 522, Issue 4, pp.4963-4971

Massive black holes in galactic nuclei: Theory and Simulations

Authors: Tiziana Di Matteo, Daniel Angles-Alcazar, Francesco Shankar

Abstract: Massive black holes are fundamental constituents of our cosmos, from the Big Bang to today. Understanding their formation from cosmic dawn, their growth, and the emergence of the first, rare quasars in the early Universe remains one of our greatest theoretical and observational challenges. Hydrodynamic cosmological simulations self-consistently combine the processes of structure formation at cosmo… ▽ More Massive black holes are fundamental constituents of our cosmos, from the Big Bang to today. Understanding their formation from cosmic dawn, their growth, and the emergence of the first, rare quasars in the early Universe remains one of our greatest theoretical and observational challenges. Hydrodynamic cosmological simulations self-consistently combine the processes of structure formation at cosmological scales with the physics of smaller, galaxy scales. They capture our most realistic understanding of massive black holes and their connection to galaxy formation and have become the primary avenue for theoretical research in this field. The space-based gravitational wave interferometer, LISA, will open up new investigations into the dynamical processes involving massive black holes. Multi-messenger astrophysics brings new exciting prospects for tracing the origin, growth and merger history of massive black holes across cosmic ages. △ Less

Submitted 23 April, 2023; originally announced April 2023.

Comments: 77 pages, 20 figures, This chapter is the preprint of the version currently in production. Please cite this chapter as the following: T.DiMatteo, D. Angles-Alcazar, and F. Shankar. Massive black holes in galactic nuclei: Theory and simulations, in The Encyclopedia of Cosmology (Set 2): Black Holes, edited by Z. Haiman (World Scientific, New Jersey, 2023)

MSC Class: none

The CAMELS project: Expanding the galaxy formation model space with new ASTRID and 28-parameter TNG and SIMBA suites

Authors: Yueying Ni, Shy Genel, Daniel Anglés-Alcázar, Francisco Villaescusa-Navarro, Yongseok Jo, Simeon Bird, Tiziana Di Matteo, Rupert Croft, Nianyi Chen, Natalí S. M. de Santi, Matthew Gebhardt, Helen Shao, Shivam Pandey, Lars Hernquist, Romeel Dave

Abstract: 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.… ▽ More 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 ($Ω_m$, $σ_8$, $Ω_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. △ Less

Submitted 4 April, 2023; originally announced April 2023.

Direct formation of massive black holes via dynamical collapse in metal-enriched merging galaxies at $z \sim 10$: fully cosmological simulations

Authors: Lucio Mayer, Pedro R. Capelo, Lorenz Zwick, Tiziana Di Matteo

Abstract: We present the results of the first fully cosmological hydrodynamical simulations studying the merger-driven model for massive black hole (BH) seed formation via direct collapse. Using the zoom-in technique as well as particle splitting, we achieve a final spatial resolution of $2$ pc. We show that the major merger of two massive galaxies at redshift $z \sim 8$ results in the formation of a nuclea… ▽ More We present the results of the first fully cosmological hydrodynamical simulations studying the merger-driven model for massive black hole (BH) seed formation via direct collapse. Using the zoom-in technique as well as particle splitting, we achieve a final spatial resolution of $2$ pc. We show that the major merger of two massive galaxies at redshift $z \sim 8$ results in the formation of a nuclear supermassive disk (SMD) of only $4$ pc in radius, owing to a prodigious gas inflow sustained at $100$-$1000$ $M_{\odot}$ yr$^{-1}$. The core of the merger remnant is metal-rich, well above solar abundance, and the SMD reaches a gaseous mass of $3 \times 10^8$ $M_{\odot}$ in less than a million years after the merger, despite a concurrent prominent nuclear starburst. Dynamical heating as gas falls into the deepest part of the potential well, and heating and stirring by supernova blastwaves, generate a turbulent multi-phase interstellar medium, with a gas velocity dispersion exceeding 100 km s$^{-1}$. As a result, only moderate fragmentation occurs in the inner $10$-$20$ pc despite the temperature falls below $1000$ K. The SMD is Jeans-unstable as well as bar-unstable and will collapse further adiabatically, becoming warm and ionized. We show that the SMD, following inevitable contraction, will become general relativistic unstable and directly form a supermassive BH of mass in the range $10^6$-$10^8$ $M_{\odot}$, essentially skipping the stage of BH seed formation. These results confirm that mergers between the most massive galaxies at $z \sim 8$-$10$ can naturally explain the rapid emergence of bright high-redshift quasars. △ Less

Submitted 4 April, 2023; originally announced April 2023.

Comments: 14 pages, 8 figures, submitted to ApJ

Boosting Line Intensity Map Signal-to-Noise with the Ly-$α$ Forest Cross-Correlation

Authors: Mahdi Qezlou, Simeon Bird, Adam Lidz, Guochao Sun, Andrew B. Newman, Gwen C. Rudie, Yueying Ni, Rupert Croft, Tiziana Di Matteo

Abstract: We forecast the prospects for cross-correlating future line intensity mapping (LIM) surveys with the current and future Ly-$α$ forest data. We use large cosmological hydrodynamic simulations to model the expected emission signal for the CO rotational transition in the COMAP LIM experiment at the 5-year benchmark and the Ly-$α$ forest absorption signal for various surveys, including eBOSS, DESI, an… ▽ More We forecast the prospects for cross-correlating future line intensity mapping (LIM) surveys with the current and future Ly-$α$ forest data. We use large cosmological hydrodynamic simulations to model the expected emission signal for the CO rotational transition in the COMAP LIM experiment at the 5-year benchmark and the Ly-$α$ forest absorption signal for various surveys, including eBOSS, DESI, and PFS. We show that CO$\times$Ly-$α$ forest can significantly enhance the detection signal-to-noise ratio of CO, with a $200$ to $300 \%$ improvement when cross-correlated with the forest observed in the Prime Focus Spectrograph (PFS) survey and a $50$ to $75\%$ enhancement for the currently available eBOSS or the upcoming DESI observations. We compare to the signal-to-noise improvements expected for a galaxy survey and show that CO$\times$Ly-$α$ is competitive with even a spectroscopic galaxy survey in raw signal-to-noise. Furthermore, our study suggests that the clustering of CO emission is tightly constrained by CO$\times$Ly-$α$ forest, due to the increased signal-to-noise ratio and the simplicity of Ly-$α$ absorption power spectrum modeling. Any foreground contamination or systematics are expected not to be shared between LIM surveys and Ly-$α$ forest observations; this provides an unbiased inference. Our findings highlight the potential benefits of utilizing the Ly-$α$ forest to aid in the initial detection of signals in line intensity experiments. For example, we also estimate that [CII]$\times$Ly-$α$ forest measurements from EXCLAIM and DESI/eBOSS, respectively, should have a larger signal-to-noise ratio than planned [CII]$\times$quasar observations by about an order of magnitude. Our results can be readily applied to actual data thanks to the observed quasar spectra in eBOSS Stripe 82, which overlaps with several LIM surveys. △ Less

Submitted 1 November, 2023; v1 submitted 30 March, 2023; originally announced March 2023.

Comments: Codes and the produced data are available at https://github.com/qezlou/lali

Triple and Quadruple Black Holes in the ASTRID Simulation at $z \sim 2$

Authors: Calvin Hoffman, Nianyi Chen, Tiziana Di Matteo, Yueying Ni, Simeon Bird, Rupert Croft, Abraham Loeb

Abstract: 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 $Δr_{\rm max} < 200\,{\rm kpc}$. The tuple systems span a range of separations with the… ▽ More 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 $Δr_{\rm max} < 200\,{\rm kpc}$. The tuple systems span a range of separations with the majority of the observable AGN systems at $Δ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\%$. △ Less

Submitted 8 March, 2023; originally announced March 2023.

Comments: 10 pages, 9 figures; comments welcome

High-redshift supermassive black hole mergers in simulations with dynamical friction modelling

Authors: Colin DeGraf, Nianyi Chen, Yueying Ni, Tiziana Di Matteo, Simeon Bird, Michael Tremmel, Rupert Croft

Abstract: 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 si… ▽ More 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. △ Less

Submitted 1 February, 2023; originally announced February 2023.

Comments: 11 pages, 10 figures, submitted to MNRAS

Inferring the impact of feedback on the matter distribution using the Sunyaev Zel'dovich effect: Insights from CAMELS simulations and ACT+DES data

Authors: Shivam Pandey, Kai Lehman, Eric J. Baxter, Yueying Ni, Daniel Anglés-Alcázar, Shy Genel, Francisco Villaescusa-Navarro, Ana Maria Delgado, Tiziana di Matteo

Abstract: Feedback from active galactic nuclei and stellar processes changes the matter distribution on small scales, leading to significant systematic uncertainty in weak lensing constraints on cosmology. We investigate how the observable properties of group-scale halos can constrain feedback's impact on the matter distribution using Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS). Ex… ▽ More Feedback from active galactic nuclei and stellar processes changes the matter distribution on small scales, leading to significant systematic uncertainty in weak lensing constraints on cosmology. We investigate how the observable properties of group-scale halos can constrain feedback's impact on the matter distribution using Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS). Extending the results of previous work to smaller halo masses and higher wavenumber, $k$, we find that the baryon fraction in halos contains significant information about the impact of feedback on the matter power spectrum. We explore how the thermal Sunyaev Zel'dovich (tSZ) signal from group-scale halos contains similar information. Using recent Dark Energy Survey (DES) weak lensing and Atacama Cosmology Telescope (ACT) tSZ cross-correlation measurements and models trained on CAMELS, we obtain $10\%$ constraints on feedback effects on the power spectrum at $k \sim 5\, h/{\rm Mpc}$. We show that with future surveys, it will be possible to constrain baryonic effects on the power spectrum to $\mathcal{O}(<1\%)$ at $k = 1\, h/{\rm Mpc}$ and $\mathcal{O}(3\%)$ at $k = 5\, h/{\rm Mpc}$ using the methods that we introduce here. Finally, we investigate the impact of feedback on the matter bispectrum, finding that tSZ observables are highly informative in this case. △ Less

Submitted 5 January, 2023; originally announced January 2023.

Comments: 18 pages, 15 figures. Comments are welcome

Orbital and Radiative Properties of Wandering Intermediate-Mass Black Holes in the ASTRID Simulation

Authors: Emma Jane Weller, Fabio Pacucci, Yueying Ni, Nianyi Chen, Tiziana Di Matteo, Magdalena Siwek, Lars Hernquist

Abstract: 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 populat… ▽ More 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 \, μ\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. △ Less

Submitted 2 April, 2023; v1 submitted 28 October, 2022; originally announced October 2022.

Comments: Accepted for publication in MNRAS. This is the final version of the manuscript. 9 pages, 7 figures

A vast population of wandering and merging IMBHs at cosmic noon

Authors: Tiziana Di Matteo, Yueying Ni, Nianyi Chen, Rupert Croft, Simeon Bird, Fabio Pacucci, Angelo Ricarte, Michael Tremmel

Abstract: 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… ▽ More 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. △ Less

Submitted 26 October, 2022; originally announced October 2022.

Super-resolution simulation of the Fuzzy Dark Matter cosmological model

Authors: Meris Sipp, Patrick LaChance, Rupert Croft, Yueying Ni, Tiziana Di Matteo

Abstract: AI super-resolution, combining deep learning and N-body simulations has been shown to successfully reproduce the large scale structure and halo abundances in the Lambda Cold Dark Matter cosmological model. Here, we extend its use to models with a different dark matter content, in this case Fuzzy Dark Matter (FDM), in the approximation that the difference is encoded in the initial power spectrum. W… ▽ More AI super-resolution, combining deep learning and N-body simulations has been shown to successfully reproduce the large scale structure and halo abundances in the Lambda Cold Dark Matter cosmological model. Here, we extend its use to models with a different dark matter content, in this case Fuzzy Dark Matter (FDM), in the approximation that the difference is encoded in the initial power spectrum. We focus on redshift z = 2, with simulations that model smaller scales and lower masses, the latter by two orders of magnitude, than has been done in previous AI super-resolution work. We find that the super-resolution technique can reproduce the power spectrum and halo mass function to within a few percent of full high resolution calculations. We also find that halo artifacts, caused by spurious numerical fragmentation of filaments, are equally present in the super-resolution outputs. Although we have not trained the super-resolution algorithm using full quantum pressure FDM simulations, the fact that it performs well at the relevant length and mass scales means that it has promise as technique which could avoid the very high computational cost of the latter, in some contexts. We conclude that AI super-resolution can become a useful tool to extend the range of dark matter models covered in mock catalogs. △ Less

Submitted 23 October, 2022; originally announced October 2022.

Comments: 7 pages, 4 figures

A close quasar pair in a disk-disk galaxy merger at z = 2.17

Authors: Yu-Ching Chen, Xin Liu, Adi Foord, Yue Shen, Masamune Oguri, Nianyi Chen, Tiziana Di Matteo, Miguel Holgado, Hsiang-Chih Hwang, Nadia Zakamska

Abstract: 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 nu… ▽ More 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. △ Less

Submitted 5 April, 2023; v1 submitted 22 September, 2022; originally announced September 2022.

Comments: 88 pages, 18 figures, 8 tables; published

Journal ref: Nature, 616, 45 (2023)

Statistics of Galactic-Scale Quasar Pairs at Cosmic Noon

Authors: Yue Shen, Hsiang-Chih Hwang, Masamune Oguri, Nianyi Chen, Tiziana Di Matteo, Yueying Ni, Simeon Bird, Nadia Zakamska, Xin Liu, Yu-Ching Chen, Kaitlin M. Kratter

Abstract: 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… ▽ More 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<z<3.5$, using Gaia EDR3-resolved pairs around SDSS DR16 quasars. The measurement was based on a sample of 60 Gaia-resolved double quasars (out of 487 Gaia pairs dominated by quasar+star superpositions) at these separations, corrected for pair completeness in Gaia, which we quantify as functions of pair separation, magnitude of the primary, and magnitude contrast. The double quasar fraction increases towards smaller separations by a factor of $\sim 5$ over these scales. The division between physical quasar pairs and lensed quasars in our sample is currently unknown, requiring dedicated follow-up observations (in particular, deep, sub-arcsec-resolution IR imaging for the closest pairs). Intriguingly, at this point the observed pair statistics are in rough agreement with theoretical predictions both for the lensed quasar population in mock catalogs and for dual quasars in cosmological hydrodynamic simulations. Upcoming wide-field imaging/spectroscopic space missions such as Euclid, CSST and Roman, combined with targeted follow-up observations, will conclusively measure the abundances and host galaxy properties of galactic-scale quasar pairs, offset AGNs, and sub-arcsec lensed quasars across cosmic time. △ Less

Submitted 9 August, 2022; originally announced August 2022.

Comments: 19 pages, 9 figures; submitted to ApJ

Journal ref: ApJ, 943, 38 (2023)

Properties and Evolution of Dual and Offset AGN in the ASTRID Simulation at $z \sim 2$

Authors: Nianyi Chen, Tiziana Di Matteo, Yueying Ni, Michael Tremmel, Colin DeGraf, Yue Shen, A. Miguel Holgado, Simeon Bird, Rupert Croft, Yu Feng

Abstract: We examine the dual (both BHs active) and offset (one BH active) AGN population (comprising $\sim$ 2000 pairs at $0.5\,\text{kpc}\lesssimΔ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$. Com… ▽ More We examine the dual (both BHs active) and offset (one BH active) AGN population (comprising $\sim$ 2000 pairs at $0.5\,\text{kpc}\lesssimΔ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 $Δr\sim5\,\text{kpc}$ for a few hundred Myrs. △ Less

Submitted 9 August, 2022; originally announced August 2022.

Comments: 17 Pages, 18 Figures; comments welcome!

The BlueTides Mock Image Catalogue: Simulated observations of high-redshift galaxies and predictions for JWST imaging surveys

Authors: Madeline A. Marshall, Katelyn Watts, Stephen Wilkins, Tiziana Di Matteo, Jussi K. Kuusisto, William J. Roper, Aswin P. Vijayan, Yueying Ni, Yu Feng, Rupert A. C. Croft

Abstract: We present a mock image catalogue of ~100,000 MUV=-22.5 to -19.6 mag galaxies at z=7-12 from the BlueTides cosmological simulation. We create mock images of each galaxy with the James Webb (JWST), Hubble, Roman, and Euclid Space Telescopes, as well as Subaru, and VISTA, with a range of near- and mid-infrared filters. We perform photometry on the mock images to estimate the success of these instrum… ▽ More We present a mock image catalogue of ~100,000 MUV=-22.5 to -19.6 mag galaxies at z=7-12 from the BlueTides cosmological simulation. We create mock images of each galaxy with the James Webb (JWST), Hubble, Roman, and Euclid Space Telescopes, as well as Subaru, and VISTA, with a range of near- and mid-infrared filters. We perform photometry on the mock images to estimate the success of these instruments for detecting high-z galaxies. We predict that JWST will have unprecedented power in detecting high-z galaxies, with a 95% completeness limit at least 2.5 magnitudes fainter than VISTA and Subaru, 1.1 magnitudes fainter than Hubble, and 0.9 magnitudes fainter than Roman, for the same wavelength and exposure time. Focusing on JWST, we consider a range of exposure times and filters, and find that the NIRCam F356W and F277W filters will detect the faintest galaxies, with 95% completeness at m=27.4 mag in 10ks exposures. We also predict the number of high-z galaxies that will be discovered by upcoming JWST imaging surveys. We predict that the COSMOS-Web survey will detect ~1000 MUV<-20.1 mag galaxies at 6.5<z<7.5, by virtue of its large survey area. JADES-Medium will detect almost 100% of MUV<-20 mag galaxies at z<8.5 due to its significant depth, however with its smaller survey area it will detect only ~100 of these galaxies at 6.5<z<7.5. Cosmic variance results in a large range in the number of predicted galaxies each survey will detect, which is more evident in smaller surveys such as CEERS and the PEARLS NEP and GOODS-S fields. △ Less

Submitted 25 August, 2022; v1 submitted 17 June, 2022; originally announced June 2022.

Comments: Published in MNRAS. The BlueTides Mock Image Catalogue is publicly available as a High Level Science Product (HLSP) via the Mikulski Archive for Space Telescopes (MAST), at https://doi.org/10.17909/er09-4527. Example codes for using the catalogue data are available at https://github.com/madelinemarshall/BlueTidesMockImageCatalogue

Probing the $z\gtrsim6$ quasars in a universe with IllustrisTNG physics: Impact of gas-based black hole seeding models

Authors: Aklant Kumar Bhowmick, Laura Blecha, Yueying Ni, Tiziana Di Matteo, Paul Torrey, Luke Zoltan Kelley, Mark Vogelsberger, Rainer Weinberger, Lars Hernquist

Abstract: We explore implications of a range of black hole (BH) seeding prescriptions on the formation of the brightest $z\gtrsim6$ quasars in cosmological hydrodynamic simulations. The underlying galaxy formation model is the same as in IllustrisTNG. Using constrained initial conditions, we study the growth of BHs in rare overdense regions (forming $\gtrsim10^{12}M_{\odot}/h$ halos by $z=7$) using a… ▽ More We explore implications of a range of black hole (BH) seeding prescriptions on the formation of the brightest $z\gtrsim6$ quasars in cosmological hydrodynamic simulations. The underlying galaxy formation model is the same as in IllustrisTNG. Using constrained initial conditions, we study the growth of BHs in rare overdense regions (forming $\gtrsim10^{12}M_{\odot}/h$ halos by $z=7$) using a $(9~\mathrm{Mpc}/h)^3$ simulated volume. BH growth is maximal within halos that are compact and have a low tidal field. For these halos, we consider an array of gas-based seeding prescriptions wherein $M_{\mathrm{seed}}=10^4-10^6~M_{\odot}/h$ seeds are inserted in halos above critical thresholds for halo mass and dense, metal-poor gas mass (defined as $\tilde{M}_{\mathrm{h}}$ and $\tilde{M}_{\mathrm{sf,mp}}$, respectively, in units of $M_{\mathrm{seed}}$). We find that a seed model with $\tilde{M}_{\mathrm{sf,mp}}=5$ and $\tilde{M}_{\mathrm{h}}=3000$ successfully produces a $z\sim6$ quasar with $\sim10^9~M_{\odot}$ mass and $\sim10^{47}~\mathrm{ergs~s^ {-1}}$ luminosity. BH mergers play a crucial role at $z\gtrsim9$, causing an early boost in BH mass at a time when accretion-driven BH growth is negligible. When more stringent seeding conditions are applied (for e.g., $\tilde{M}_{\mathrm{sf,mp}}=1000$), the relative paucity of BH seeds results in a much lower merger rate. In this case, $z\gtrsim6$ quasars can only be formed if we enhance the maximum allowed BH accretion rates (by factors $\gtrsim10$) compared to the accretion model used in IllustrisTNG. This can be achieved either by allowing for super-Eddington accretion, or by reducing the radiative efficiency. Our results show that progenitors of $z\sim6$ quasars have distinct BH merger histories for different seeding models, which will be distinguishable with LISA observations. △ Less

Submitted 11 May, 2022; originally announced May 2022.

Comments: 21 pages, 15 figures

Concordance between observations and simulations in the evolution of the mass relation between supermassive black holes and their host galaxies

Authors: Xuheng Ding, John D. Silverman, Tommaso Treu, Junyao Li, Aklant K. Bhowmick, Nicola Menci, Marta Volonteri, Laura Blecha, Tiziana Di Matteo, Yohan Dubois

Abstract: We carry out a comparative analysis of the relation between the mass of supermassive black holes (BHs) and the stellar mass of their host galaxies at $0.2<z<1.7$ using well-matched observations and multiple state-of-the-art simulations (e.g., Massive Black II, Horizon-AGN, Illustris, TNG and a semi-analytic model). The observed sample consists of 646 uniformly-selected SDSS quasars (… ▽ More We carry out a comparative analysis of the relation between the mass of supermassive black holes (BHs) and the stellar mass of their host galaxies at $0.2<z<1.7$ using well-matched observations and multiple state-of-the-art simulations (e.g., Massive Black II, Horizon-AGN, Illustris, TNG and a semi-analytic model). The observed sample consists of 646 uniformly-selected SDSS quasars ($0.2 < z < 0.8$) and 32 broad-line active galactic nuclei (AGNs; $1.2<z<1.7$) with imaging from Hyper Suprime-Cam (HSC) for the former and Hubble Space Telescope (HST) for the latter. We first add realistic observational uncertainties to the simulation data and then construct a simulated sample in the same manner as the observations. Over the full redshift range, our analysis demonstrates that all simulations predict a level of intrinsic scatter of the scaling relations comparable to the observations which appear to agree with the dispersion of the local relation. Regarding the mean relation, Horizon-AGN and TNG are in closest agreement with the observations at low and high redshift ($z\sim$ 0.2 and 1.5, respectively) while the other simulations show subtle differences within the uncertainties. For insight into the physics involved, the scatter of the scaling relation, seen in the SAM, is reduced by a factor of two and closer to the observations after adopting a new feedback model that considers the geometry of the AGN outflow. The consistency in the dispersion with redshift in our analysis supports the importance of both quasar- and radio-mode feedback prescriptions in the simulations. Finally, we highlight the importance of increasing the sensitivity (e.g., using the James Webb Space Telescope), thereby pushing to lower masses and minimizing biases due to selection effects. △ Less

Submitted 21 July, 2022; v1 submitted 9 May, 2022; originally announced May 2022.

Comments: 17 pages, 9 figures, Published by ApJ

Journal ref: 2022, ApJ, 933, 132D

An estimate of the stochastic gravitational wave background from the MassiveBlackII simulation

Authors: Bailey Sykes, Hannah Middleton, Andrew Melatos, Tiziana Di Matteo, Colin DeGraf, Aklant Bhowmick

Abstract: A population of super-massive black hole binaries is expected to generate a stochastic gravitational wave background (SGWB) in the pulsar timing array (PTA) frequency range of $10^{-9}$--$10^{-7}$ 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 super-massive black hole binary mergers from the Ma… ▽ More A population of super-massive black hole binaries is expected to generate a stochastic gravitational wave background (SGWB) in the pulsar timing array (PTA) frequency range of $10^{-9}$--$10^{-7}$ 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 super-massive black hole binary mergers from the MassiveBlackII simulation to estimate the characteristic strain of the stochastic background. We examine both a gravitational wave driven model of binary evolution and a model which also includes the effects of stellar scattering and a circumbinary gas disk. Results are consistent with PTA upper limits and similar to estimates in the literature. The characteristic strain at a reference frequency of $1 yr^{-1}$ is found to be $A_{yr^{-1}} = 6.9 \times 10^{-16}$ and $A_{yr^{-1}} = 6.4 \times 10^{-16}$ in the gravitational-wave driven and stellar scattering/gas disk cases, respectively. Using the latter approach, our models show that the SGWB is mildly suppressed compared to the purely gravitational wave driven model as frequency decreases inside the PTA frequency band. △ Less

Submitted 10 February, 2022; originally announced February 2022.

Comments: 10 pages, 3 figures, 1 table, accepted for publication in MNRAS

Massive Black Hole Mergers with Orbital Information: Predictions from the ASTRID Simulation

Authors: Nianyi Chen, Yueying Ni, A. Miguel Holgado, Tiziana Di Matteo, Michael Tremmel, Colin DeGraf, Simeon Bird, Rupert Croft, Yu Feng

Abstract: 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… ▽ More 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. △ Less

Submitted 15 December, 2021; originally announced December 2021.

Comments: 17 pages, 13 Figures; comments are welcome

The ASTRID Simulation: Galaxy Formation and Reionization

Authors: Simeon Bird, Yueying Ni, Tiziana Di Matteo, Rupert Croft, Yu Feng, Nianyi Chen

Abstract: 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 hi… ▽ More 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. △ Less

Submitted 23 March, 2022; v1 submitted 1 November, 2021; originally announced November 2021.

Comments: 15 pages, 12 figures. Accepted by MNRAS, added star formation rate density

The ASTRID simulation: the evolution of Supermassive Black Holes

Authors: Yueying Ni, Tiziana Di Matteo, Simeon Bird, Rupert Croft, Yu Feng, Nianyi Chen, Michael Tremmel, Colin DeGraf, Yin Li

Abstract: 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… ▽ More 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. △ Less

Submitted 20 March, 2022; v1 submitted 27 October, 2021; originally announced October 2021.

Comments: Match the published version

The Impact of Dust on the Sizes of Galaxies in the Epoch of Reionization

Authors: Madeline A. Marshall, Stephen Wilkins, Tiziana Di Matteo, William J. Roper, Aswin P. Vijayan, Yueying Ni, Yu Feng, Rupert A. C. Croft

Abstract: We study the sizes of galaxies in the Epoch of Reionization using a sample of ~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 which could enclose the pixels containing 50% of the total mass/light in the image. We find an inverse relation… ▽ More We study the sizes of galaxies in the Epoch of Reionization using a sample of ~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 which could enclose the pixels containing 50% 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 sight lines 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 evolution from z=7 to 11, and galaxy size evolves with redshift as $R\propto(1+z)^{-m}$, where $m=0.662\pm0.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

Submitted 9 February, 2022; v1 submitted 22 October, 2021; originally announced October 2021.

Comments: 12 pages, 10 figures. Accepted for publication in MNRAS

On the detectability of massive black hole merger eventsby LISA

Authors: Samuel Banks, Katharine Lee, Nazanin Azimi, Kendall Scarborough, Nikolai Stefanov, Indra Periwal, Colin DeGraf, Tiziana Di Matteo

Abstract: The launch of space based gravitational wave (GW) detectors (e.g. Laser Interferometry Space Antenna; LISA) and current and upcoming Pulsar Timing Arrays (PTAs) 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 th… ▽ More The launch of space based gravitational wave (GW) detectors (e.g. Laser Interferometry Space Antenna; LISA) and current and upcoming Pulsar Timing Arrays (PTAs) 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 MBH 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}} > 10^6 M_\odot$ MBHBs (the lowest mass in Illustris) decreasing from $\sim 3$ yr$^{-1}$ to $\sim 0.1 $yr$^{-1}$, and shifting the merger peak from z $\sim 2$ to $\sim 1.25$. During this time, we estimate that accretion grows the total merging mass by as much as 7x from the original mass. Importantly, however, dynamical friction associated delays (which shift the mergers toward lower-redshift and higher-masses) lead 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 $M_{\rm{BH}}$ growth, mergers still tend to lie on the typical $M_{\rm{BH}}-M_*$ relation, but with an increased likelihood of an undermassive black hole. △ Less

Submitted 19 July, 2021; originally announced July 2021.

Comments: 13 pages, 10 figures. Submitted to MNRAS

The Diversity of Environments Around Luminous Quasars at Redshift $z \sim 6$

Authors: Keven Ren, Michele Trenti, Madeline A. Marshall, Tiziana Di Matteo, Yueying Ni

Abstract: Significant clustering around the rarest luminous quasars is a feature predicted by dark matter theory combined with number density matching arguments. However, this expectation is not reflected by observations of quasars residing in a diverse range of environments. Here, we assess the tension in the diverse clustering of visible $i$-band dropout galaxies around luminous $z\sim6$ quasars. Our appr… ▽ More Significant clustering around the rarest luminous quasars is a feature predicted by dark matter theory combined with number density matching arguments. However, this expectation is not reflected by observations of quasars residing in a diverse range of environments. Here, we assess the tension in the diverse clustering of visible $i$-band dropout galaxies around luminous $z\sim6$ quasars. Our approach uses a simple empirical method to derive the median luminosity to halo mass relation, $L_{c}(M_{h})$ for both quasars and galaxies under the assumption of log-normal luminosity scatter, $Σ_{Q}$ and $Σ_{G}$. We show that higher $Σ_{Q}$ reduces the average halo mass hosting a quasar of a given luminosity, thus introducing at least a partial reversion to the mean in the number count distribution of nearby Lyman-Break galaxies. We generate a large sample of mock Hubble Space Telescope fields-of-view centred across rare $z\sim6$ quasars by resampling pencil beams traced through the dark matter component of the BlueTides cosmological simulation. We find that diverse quasar environments are expected for $Σ_{Q}>0.4$, consistent with numerous observations and theoretical studies. However, we note that the average number of galaxies around the central quasar is primarily driven by galaxy evolutionary processes in neighbouring halos, as embodied by our parameter $Σ_{G}$, instead of a difference in the large scale structure around the central quasar host, embodied by $Σ_{Q}$. We conclude that models with $Σ_{G}>0.3$ are consistent with current observational constraints on high-z quasars, and that such a value is comparable to the scatter estimated from hydrodynamical simulations of galaxy formation. △ Less

Submitted 13 June, 2021; originally announced June 2021.

Comments: 19 pages, 8 figures; accepted for publication in ApJ

AI-assisted super-resolution cosmological simulations II: Halo substructures, velocities and higher order statistics

Authors: Yueying Ni, Yin Li, Patrick Lachance, Rupert A. C. Croft, Tiziana Di Matteo, Simeon Bird, Yu Feng

Abstract: In this work, we expand and test the capabilities of our recently developed super-resolution (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… ▽ More In this work, we expand and test the capabilities of our recently developed super-resolution (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 non-linear 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 Mpc/h, and examine the matter power spectra, bispectra and 2D power spectra in redshift space. We find the generated SR field matches the true HR result at percent level down to scales of k ~ 10 h/Mpc. We also identify and inspect dark matter halos and their substructures. Our SR model generate 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 catalogs. The SR technique can be a powerful and promising tool for modelling small-scale galaxy formation physics in large cosmological volumes. △ Less

Submitted 17 September, 2021; v1 submitted 3 May, 2021; originally announced May 2021.

Comments: 13 pages, 11 figures, published version

Dynamical Friction Modeling of Massive Black Holes in Cosmological Simulations and Effects on Merger Rate Predictions

Authors: Nianyi Chen, Yueying Ni, Michael Tremmel, Tiziana Di Matteo, Simeon Bird, Colin DeGraf, Yu Feng

Abstract: 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… ▽ More 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. △ Less

Submitted 2 April, 2021; v1 submitted 31 March, 2021; originally announced April 2021.

Comments: 17 pages, 12 figures; Comments welcome!

Observing the host galaxies of high-redshift quasars with JWST: predictions from the BlueTides simulation

Authors: Madeline A. Marshall, J. Stuart B. Wyithe, Rogier A. Windhorst, Tiziana Di Matteo, Yueying Ni, Stephen Wilkins, Rupert A. C. Croft, Mira Mechtley

Abstract: The bright emission from high-redshift quasars completely conceals their host galaxies in the rest-frame ultraviolet/optical, with detection of the hosts in these wavelengths eluding even the Hubble Space Telescope (HST) using detailed point spread function (PSF) modelling techniques. In this study we produce mock images of a sample of z=7 quasars extracted from the BlueTides simulation, and apply… ▽ More The bright emission from high-redshift quasars completely conceals their host galaxies in the rest-frame ultraviolet/optical, with detection of the hosts in these wavelengths eluding even the Hubble Space Telescope (HST) using detailed point spread function (PSF) modelling techniques. In this study we produce mock images of a sample of z=7 quasars extracted from the BlueTides simulation, and apply Markov Chain Monte Carlo-based PSF modelling to determine the detectability of their host galaxies with the James Webb Space Telescope (JWST). While no statistically significant detections are made with HST, we predict that at the same wavelengths and exposure times JWST NIRCam imaging will detect ~50% of quasar host galaxies. We investigate various observational strategies, and find that NIRCam wide-band imaging in the long-wavelength filters results in the highest fraction of successful quasar host detections, detecting >80% of the hosts of bright quasars in exposure times of 5 ks. Exposure times of ~5 ks are required to detect the majority of host galaxies in the NIRCam wide-band filters, however even 10 ks exposures with MIRI result in <30% successful host detections. We find no significant trends between galaxy properties and their detectability. The PSF modelling can accurately recover the host magnitudes, radii, and spatial distribution of the larger-scale emission, when accounting for the central core being contaminated by residual quasar flux. Care should be made when interpreting the host properties measured using PSF modelling. △ Less

Submitted 9 June, 2021; v1 submitted 4 January, 2021; originally announced January 2021.

Comments: 21 pages, 14 figures. Resubmitted to MNRAS following referee feedback, with most notable improvements in Section 5

Not all peaks are created equal: the early growth of Supermassive Black Holes

Authors: Yueying Ni, Tiziana Di Matteo, Yu Feng

Abstract: In this work, we use the constrained Gaussian realization technique to study the early growth of supermassive black holes (SMBHs) in cosmological hydrodynamic simulations, exploring its relationship with features of the initial density peaks on large scales, ~1 Mpc/h. Our constrained simulations of volume (20 Mpc/h)^3 successfully reconstruct the large-scale structure as well as the black hole gro… ▽ More In this work, we use the constrained Gaussian realization technique to study the early growth of supermassive black holes (SMBHs) in cosmological hydrodynamic simulations, exploring its relationship with features of the initial density peaks on large scales, ~1 Mpc/h. Our constrained simulations of volume (20 Mpc/h)^3 successfully reconstruct the large-scale structure as well as the black hole growth for the hosts of the rare 10^9 Msun SMBHs found in the BlueTides simulation at z~7. We run a set of simulations with constrained initial conditions by imposing a 5 σ_0(R_G) peak on scale of R_G = 1 Mpc/h varying different peak features, such as the shape and compactness as well as the tidal field surrounding the peak. We find that initial density peaks with high compactness and low tidal field induce the most rapid BH growth at early epochs. This is because compact density peaks with a more spherical large scale matter distribution lead to the formation of high density gas clumps in the centers of halos, and thus boost early BH accretion. Moreover, such initially compact density peaks in low tidal field regions also lead to a more compact BH host galaxy morphology. This can explain the tight correlation between BH growth and host galaxy compactness seen in observations. △ Less

Submitted 6 November, 2021; v1 submitted 8 December, 2020; originally announced December 2020.

Comments: matches MNRAS accepted version

Morphological evolution of supermassive black hole merger hosts and multimessenger signatures

Authors: Colin DeGraf, Debora Sijacki, Tiziana Di Matteo, Kelly Holley-Bockelmann, Greg Snyder, Volker Springel

Abstract: With projects such as Laser Interferometer Space Antenna (LISA) and Pulsar Timing Arrays expected to detect gravitational waves from supermassive black hole mergers in the near future, it is key that we understand what we expect those detections to be, and maximize what we can learn from them. To address this, we study the mergers of supermassive black holes in the Illustris simulation, the overal… ▽ More With projects such as Laser Interferometer Space Antenna (LISA) and Pulsar Timing Arrays expected to detect gravitational waves from supermassive black hole mergers in the near future, it is key that we understand what we expect those detections to be, and maximize what we can learn from them. To address this, we study the mergers of supermassive black holes in the Illustris simulation, the overall rate of mergers, and the correlation between merging black holes and their host galaxies. We find that these mergers occur in typical galaxies along the $M_{\rm{BH}}-M_*$ relation, and that between LISA and PTAs we expect to probe the full range of galaxy masses. As galaxy mergers can trigger increased star formation, we find that galaxies hosting low-mass black hole mergers tend to show a slight increase in star formation rates compared to a mass-matched sample. However, high-mass merger hosts have typical star formation rates, due to a combination of low gas fractions and powerful AGN feedback. Although minor black hole mergers do not correlate with disturbed morphologies, major mergers (especially at high-masses) tend to show morphological evidence of recent galaxy mergers which survives for ~500 Myr. This is on the same scale as the infall/hardening time of the merging black holes, suggesting that electromagnetic followups to gravitational wave signals may not be able to observe this correlation. We further find that incorporating a realistic timescale delay for the black hole mergers could shift the distribution of merger masses toward higher-masses, decreasing the rate of LISA detections while increasing the rate of PTA detections. △ Less

Submitted 28 July, 2021; v1 submitted 1 December, 2020; originally announced December 2020.

Comments: 14 pages, 17 figures. Published in MNRAS

Journal ref: MNRAS (2021), 503, 3629D

Simulated X-ray Emission in Galaxy Clusters with Feedback from Active Galactic Nuclei

Authors: Rudrani Kar Chowdhury, Soumya Roy, Suchetana Chatterjee, Nishikanta Khandai, Craig L. Sarazin, Tiziana Di Matteo

Abstract: To investigate the effect of feedback from active galactic nuclei (AGN) on their surrounding medium, we study the diffuse X-ray emission from galaxy groups and clusters by coupling the Astrophysical Plasma Emission Code (APEC) with the cosmological hydrodynamic simulation involving AGN feedback. We construct a statistical sample of synthetic Chandra X-ray photon maps to observationally characteriz… ▽ More To investigate the effect of feedback from active galactic nuclei (AGN) on their surrounding medium, we study the diffuse X-ray emission from galaxy groups and clusters by coupling the Astrophysical Plasma Emission Code (APEC) with the cosmological hydrodynamic simulation involving AGN feedback. We construct a statistical sample of synthetic Chandra X-ray photon maps to observationally characterize the effect of AGN on the ambient medium. We show that AGN are effective in displacing the hot X-ray emitting gas from the centers of groups and clusters, and that these signatures remain evident in observations of the X-ray surface brightness profiles. △ Less

Submitted 11 November, 2020; originally announced November 2020.

Comments: Accepted for publication in Astronomische Nachrichten, proceedings of the 9th International Workshop on Astronomy and Relativistic Astrophysics

Mocking Faint Black Holes during Reionization

Authors: Marius B. Eide, Benedetta Ciardi, Yu Feng, Tiziana Di Matteo

Abstract: To investigate the potential abundance and impact of nuclear black holes (BHs) during reionization, we generate a neural network that estimates their masses and accretion rates by training it on 23 properties of galaxies harbouring them at $z=6$ in the cosmological hydrodynamical simulation Massive-Black II. We then populate all galaxies in the simulation from $z=18$ to $z=5$ with BHs from this ne… ▽ More To investigate the potential abundance and impact of nuclear black holes (BHs) during reionization, we generate a neural network that estimates their masses and accretion rates by training it on 23 properties of galaxies harbouring them at $z=6$ in the cosmological hydrodynamical simulation Massive-Black II. We then populate all galaxies in the simulation from $z=18$ to $z=5$ with BHs from this network. As the network allows to robustly extrapolate to BH masses below those of the BH seeds, we predict a population of faint BHs with a turnover-free luminosity function, while retaining the bright (and observed) BHs, and together they predict a Universe in which intergalactic hydrogen is $15\%$ ionized at $z=6$ for a clumping factor of 5. Faint BHs may play a stronger role in H reionization without violating any observational constraints. This is expected to have an impact also on pre-heating and -ionization, which is relevant to observations of the 21 cm line from neutral H. We also find that BHs grow more efficiently at higher $z$, but mainly follow a redshift-independent galaxy-BH relation. We provide a power law parametrisation of the hydrogen ionizing emissivity of BHs. △ Less

Submitted 19 October, 2020; originally announced October 2020.

Comments: 8 pages, 7 figures, accepted for publication in MNRAS

AI-assisted super-resolution cosmological simulations

Authors: Yin Li, Yueying Ni, Rupert A. C. Croft, Tiziana Di Matteo, Simeon Bird, Yu Feng

Abstract: Cosmological simulations of galaxy formation are limited by finite computational resources. We draw from the ongoing rapid advances in Artificial Intelligence (specifically Deep Learning) to address this problem. Neural networks have been developed to learn from high-resolution (HR) image data, and then make accurate super-resolution (SR) versions of different low-resolution (LR) images. We apply… ▽ More Cosmological simulations of galaxy formation are limited by finite computational resources. We draw from the ongoing rapid advances in Artificial Intelligence (specifically Deep Learning) to address this problem. Neural networks have been developed to learn from high-resolution (HR) image data, and then make accurate super-resolution (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 new 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\,h^{-1}\mathrm{Mpc}$, and the HR halo mass function to within $10 \%$ down to $10^{11} \, M_\odot$. We successfully deploy the model in a box 1000 times larger than the training simulation box, showing that high-resolution mock surveys can be generated rapidly. We conclude that AI assistance has the potential to revolutionize modeling of small-scale galaxy formation physics in large cosmological volumes. △ Less

Submitted 4 May, 2021; v1 submitted 13 October, 2020; originally announced October 2020.

Comments: 10 pages, 8 figures, matches PNAS accepted version

Journal ref: PNAS May 11, 2021 118 (19) e2022038118

Large scale simulations of H and He reionization and heating driven by stars and more energetic sources

Authors: Marius B. Eide, Benedetta Ciardi, Luca Graziani, Philipp Busch, Yu Feng, Tiziana Di Matteo

Abstract: We present simulations of cosmic reionization and reheating from $z=18$ to $z=5$, investigating the role of stars (emitting soft UV-photons), nuclear black holes (BHs, with power-law spectra), X-ray binaries (XRBs, with hard X-ray dominated spectra), and the supernova-associated thermal bremsstrahlung of the diffuse interstellar medium (ISM, with soft X-ray spectra). We post-process the hydrodynam… ▽ More We present simulations of cosmic reionization and reheating from $z=18$ to $z=5$, investigating the role of stars (emitting soft UV-photons), nuclear black holes (BHs, with power-law spectra), X-ray binaries (XRBs, with hard X-ray dominated spectra), and the supernova-associated thermal bremsstrahlung of the diffuse interstellar medium (ISM, with soft X-ray spectra). We post-process the hydrodynamical simulation Massive-Black II (MBII) with multifrequency ionizing radiative transfer. The source properties are directly derived from the physical environment of MBII, and our only real free parameter is the ionizing escape fraction $f_{\rm esc}$. We find that, among the models explored here, the one with an escape fraction that decreases with decreasing redshift yields results most in line with observations, such as of the neutral hydrogen fraction and the Thomson scattering optical depth. Stars are the main driver of hydrogen reionization and consequently of the thermal history of the intergalactic medium (IGM). We obtain $\langle x_{\rm HII} \rangle = 0.99998$ at $z=6$ for all source types, with volume averaged temperatures $\langle T \rangle \sim 20,000~{\rm K}$. BHs are rare and negligible to hydrogen reionization, but conversely they are the only sources which can fully ionize helium, increasing local temperatures by $\sim 10^4~{\rm K}$. The thermal and ionization state of the neutral and lowly ionized hydrogen differs significantly with different source combinations, with ISM and (to a lesser extent) XRBs, playing a significant role and, as a consequence, determining the transition from absorption to emission of the 21 cm signal from neutral hydrogen. △ Less

Submitted 14 September, 2020; originally announced September 2020.

Comments: 17 pages, 19 figures, accepted for publication in MNRAS

Stochastic Processes as the Origin of the Double-Power Law Shape of the Quasar Luminosity Function

Authors: Keven Ren, Michele Trenti, Tiziana Di Matteo

Abstract: The Quasar Luminosity Function (QLF) offers insight into the early co-evolution of black holes and galaxies. It has been characterized observationally up to redshift $z\sim6$ with clear evidence of a double power-law shape, in contrast to the Schechter-like form of the underlying dark-matter halo mass function. We investigate a physical origin for the difference in these distributions by consideri… ▽ More The Quasar Luminosity Function (QLF) offers insight into the early co-evolution of black holes and galaxies. It has been characterized observationally up to redshift $z\sim6$ with clear evidence of a double power-law shape, in contrast to the Schechter-like form of the underlying dark-matter halo mass function. We investigate a physical origin for the difference in these distributions by considering the impact of stochasticity induced by the processes that determine the quasar luminosity for a given host halo and redshift. We employ a conditional luminosity function and construct the relation between median quasar magnitude versus halo mass $M_{UV,\rm{c}}(M_{\rm{h}})$ with log-normal in luminosity scatter $Σ$, and duty-cycle $ε_{\rm{DC}}$, and focus on high redshift $z\gtrsim4$. We show that, in order to reproduce the observed QLF, the $Σ=0$ abundance matching requires all of the brightest quasars to be hosted in the rarest most massive dark-matter halos (with an increasing $M_{UV,\rm{c}}/M_{\rm{h}}$ in halo mass). Conversely, for $Σ>0$ the brightest quasars can be over-luminous outliers hosted in relatively common dark-matter halos. In this case, the median quasar magnitude versus halo mass relation, $M_{UV,\rm{c}}$, flattens at the high-end, as expected in self-regulated growth due to feedback. We sample the parameter space of $Σ$ and $ε_{\rm{DC}}$ and show that $M_{UV,\rm{c}}$ flattens above $M_{\rm{h}}\sim 10^{12}M_{\odot}$ for $ε_{\rm{DC}}<10^{-2}$. Models with $ε_{\rm{DC}}\sim1$ instead require a high mass threshold close to $M_{\rm{h}}\gtrsim10^{13}M_{\odot}$. We investigate the impact of $ε_{\rm{DC}}$ and $Σ$ on measurements of clustering and find there is no luminosity dependence on clustering for $Σ>0.3$, consistent with recent observations from Subaru HSC. △ Less

Submitted 15 April, 2020; originally announced April 2020.

Comments: 12 pages, 6 figures; accepted for publication in ApJ

TXS 2116$-$077: A Gamma-ray Emitting Relativistic Jet Hosted in a Galaxy Merger

Authors: Vaidehi S. Paliya, Enrique Pérez, Rubén García-Benito, Marco Ajello, Francisco Prada, Antxon Alberdi, Hyewon Suh, C. H. Ishwara Chandra, Alberto Domínguez, Stefano Marchesi, Tiziana Di Matteo, Dieter Hartmann, Marco Chiaberge

Abstract: What triggers collimated relativistic outflows or jets, from the centers of galaxies remains a fundamental question in astrophysics. The merging of two galaxies has been proposed to realize the conditions to successfully launch and drive such jets into the intergalactic medium. However, evidences for the operation of this mechanism are scarce. Here we report the first unambiguous detection of an o… ▽ More What triggers collimated relativistic outflows or jets, from the centers of galaxies remains a fundamental question in astrophysics. The merging of two galaxies has been proposed to realize the conditions to successfully launch and drive such jets into the intergalactic medium. However, evidences for the operation of this mechanism are scarce. Here we report the first unambiguous detection of an ongoing merger of a narrow-line Seyfert 1 galaxy, TXS 2116$-$077, hosting a closely aligned, $γ$-ray emitting relativistic jet with a Seyfert 2 galaxy at a separation of $\sim$12 kpc, using the observations taken with 8.2 m Subaru telescope. Our subsequent followup observations with 10.4 m Gran Telescopio Canarias, 4.2 m William Herschel Telescope, and Chandra X-ray observatory have provided what is likely to be the first glimpse of the merging environment hosting a closely aligned relativistic jet. Our finding that the jet is considerably younger than the merger demonstrates that jet activity can be triggered by galaxy mergers and that $γ$-ray detected narrow-line Seyfert 1 galaxies represent the beginning phase of that activity. These results also highlight the crucial role of mergers in shaping the fate of galaxies in their cosmological evolution and are consistent with recent studies focused on the host galaxy imaging of this enigmatic class of active galactic nuclei. △ Less

Submitted 6 April, 2020; originally announced April 2020.

Comments: 13 pages, 7 figures, 1 table, published in the Astrophysical Journal

Group-scale intrinsic galaxy alignments in the Illustris-TNG and MassiveBlack-II simulations

Authors: Ananth Tenneti, Thomas D. Kitching, Benjamin Joachimi, Tiziana Di Matteo

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 and IllustrisTNG simulations. We find that: the shape of the satellite system in halos of mass ($> 10^{13}h^{-1}M_{\odot}$) is well aligned with the shape of the central galaxy at $z=0.06$ with the mean alignment between the maj… ▽ More 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 and IllustrisTNG simulations. We find that: the shape of the satellite system in halos of mass ($> 10^{13}h^{-1}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 $\sim Δθ= 12^{\circ}$ 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 halos 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 MassiveBlack-II. The orientation of satellite galaxies tends to 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 halos. However, the projected ellipticities do not exhibit a scale-dependent radial alignment, as has been seen in some observational measurements. △ Less

Submitted 22 January, 2021; v1 submitted 27 February, 2020; originally announced February 2020.

Comments: 15 pages, 9 figures, revised after referee comments, accepted for publication in MNRAS

Testing the fidelity of simulations of black hole - galaxy co-evolution at z ~ 1.5 with observations

Authors: Xuheng Ding, Tommaso Treu, John D. Silverman, Aklant K. Bhowmick, N. Menci, Tiziana Di Matteo

Abstract: We examine the scaling relations between the mass of a supermassive black hole (SMBH) and its host galaxy properties at $1.2<z<1.7$ using both observational data and simulations. Recent measurements of 32 X-ray-selected broad-line Active Galactic Nucleus (AGNs) are compared with two independent state-of-the-art efforts, including the hydrodynamic simulation MassiveBlackII (MBII) and a semi-analyti… ▽ More We examine the scaling relations between the mass of a supermassive black hole (SMBH) and its host galaxy properties at $1.2<z<1.7$ using both observational data and simulations. Recent measurements of 32 X-ray-selected broad-line Active Galactic Nucleus (AGNs) are compared with two independent state-of-the-art efforts, including the hydrodynamic simulation MassiveBlackII (MBII) and a semi-analytic model (SAM). After applying an observational selection function to the simulations, we find that both MBII and SAM agree well with the data, in terms of the central distribution. However, the dispersion in the mass ratio between black hole mass and stellar mass is significantly more consistent with the MBII prediction ($\sim0.3~$dex), than with the SAM ($\sim0.7~$dex), even when accounting for observational uncertainties. Hence, our observations can distinguish between the different recipes adopted in the models. The mass relations in the MBII are highly dependent on AGN feedback while the relations in the SAM are more sensitive to galaxy merger events triggering nuclear activity. Moreover, the intrinsic scatter in the mass ratio of our high-$z$ sample is comparable to that observed in the local sample, all but ruling out the proposed scenario the correlations are purely stochastic in nature arising from some sort of cosmic central limit theorem. Our results support the hypothesis of AGN feedback being responsible for a causal link between the SMBH and its host galaxy, resulting in a tight correlation between their respective masses. △ Less

Submitted 21 May, 2020; v1 submitted 18 February, 2020; originally announced February 2020.

Comments: 11 pages, 5 figures. ApJ in press. Comments are welcome

Journal ref: ApJ, (2020), 896, 159