Showing 1–40 of 40 results for author: Habouzit, M

The First Billion Years, According to JWST

Authors: Angela Adamo, Hakim Atek, Micaela B. Bagley, Eduardo Bañados, Kirk S. S. Barrow, Danielle A. Berg, Rachel Bezanson, Maruša Bradač, Gabriel Brammer, Adam C. Carnall, John Chisholm, Dan Coe, Pratika Dayal, Daniel J. Eisenstein, Jan J. Eldridge, Andrea Ferrara, Seiji Fujimoto, Anna de Graaff, Melanie Habouzit, Taylor A. Hutchison, Jeyhan S. Kartaltepe, Susan A. Kassin, Mariska Kriek, Ivo Labbé, Roberto Maiolino , et al. (24 additional authors not shown)

Abstract: With stunning clarity, JWST has revealed the Universe's first billion years. The scientific community is analyzing a wealth of JWST imaging and spectroscopic data from that era, and is in the process of rewriting the astronomy textbooks. Here, 1.5 years into the JWST science mission, we provide a snapshot of the great progress made towards understanding the initial chapters of our cosmic history.… ▽ More With stunning clarity, JWST has revealed the Universe's first billion years. The scientific community is analyzing a wealth of JWST imaging and spectroscopic data from that era, and is in the process of rewriting the astronomy textbooks. Here, 1.5 years into the JWST science mission, we provide a snapshot of the great progress made towards understanding the initial chapters of our cosmic history. We highlight discoveries and breakthroughs, topics and issues that are not yet understood, and questions that will be addressed in the coming years, as JWST continues its revolutionary observations of the Early Universe. While this compendium is written by a small number of authors, invited to ISSI Bern in March 2024 as part of the 2024 ISSI Breakthrough Workshop, we acknowledge the work of a large community that is advancing our collective understanding of the evolution of the Early Universe. △ Less

Submitted 31 May, 2024; originally announced May 2024.

Comments: review article written by the attendees of the 2024 ISSI breakthrough workshop "The first billion year of the Universe", submitted. Comments welcome

Is the James Webb Space Telescope detecting too many AGN candidates?

Authors: Melanie Habouzit

Abstract: In less than two years of operation, the James Webb Space Telescope (JWST) has already accelerated significantly our quest to identify active massive black holes (BHs) in the first billion years of the Universe's history. At the time of writing, about 50 AGN detections and candidates have been identified through spectroscopy, photometry, and/or morphology. Broad-line AGN are about a hundred times… ▽ More In less than two years of operation, the James Webb Space Telescope (JWST) has already accelerated significantly our quest to identify active massive black holes (BHs) in the first billion years of the Universe's history. At the time of writing, about 50 AGN detections and candidates have been identified through spectroscopy, photometry, and/or morphology. Broad-line AGN are about a hundred times more numerous than the faint end of the UV-bright quasar population at z~5-6. In this paper, we compare the observational constraints on the abundance of these AGN at z~5 to the populations of AGN produced in large-scale cosmological simulations. Assuming a null fraction of obscured simulated AGN, we find that while some simulations produce more AGN than discovered so far, some others produce a similar abundance or even fewer AGN in the bolometric luminosity range probed by JWST. Keeping in mind the large uncertainty on the constraints, we discuss the implications for the theoretical modeling of BH formation and evolution in case similar constraints continue to accumulate. At the redshift of interest, the simulated AGN populations diverge the most at Lbol~1e44 erg/s (by more than a dex in the bolometric luminosity function). This regime is most affected by incompleteness in JWST surveys. However, it holds significant potential for constraining the physical processes determining the assembly of BHs (e.g., seeding, feedback from supernova and AGN) and the current abundance of broad-line AGN with >1e44.5 erg/s. △ Less

Submitted 8 May, 2024; originally announced May 2024.

Comments: 11 pages, submitted to the journal

A SPectroscopic survey of biased halos In the Reionization Era (ASPIRE): Impact of Galaxies on the CGM Metal Enrichment at z > 6 Using the JWST and VLT

Authors: Siwei Zou, Zheng Cai, Feige Wang, Xiaohui Fan, Jaclyn B. Champagne, Joseph F. Hennawi, Jan-Torge Schindler, Emanuele P. Farina, Jinyi Yang, Kohei Inayoshi, Eduardo Banados, Sarah E. I. Bosman, Zihao Li, Xiaojing Lin, Yunjing Wu, Fengwu Sun, Zi-Yi Guo, Girish Kulkarni, Melanie Habouzit, Stephane Charlot, Jacopo Chevallard, Thomas Connor, Anna-Christina Eilers, Linhua Jiang, Xiangyu Jin , et al. (5 additional authors not shown)

Abstract: We characterize the multiphase circumgalactic medium and galaxy properties at z = 6.0-6.5 in four quasar fields from the James Webb Space Telescope A SPectroscopic survey of biased halos In the Reionization Era (ASPIRE) program. We use the Very Large Telescope/X-shooter spectra of quasar J0305-3150 to identify one new metal absorber at z = 6.2713 with multiple transitions (OI, MgI, FeII and CII).… ▽ More We characterize the multiphase circumgalactic medium and galaxy properties at z = 6.0-6.5 in four quasar fields from the James Webb Space Telescope A SPectroscopic survey of biased halos In the Reionization Era (ASPIRE) program. We use the Very Large Telescope/X-shooter spectra of quasar J0305-3150 to identify one new metal absorber at z = 6.2713 with multiple transitions (OI, MgI, FeII and CII). They are combined with the published absorbing systems in Davies et al. (2023a) at the same redshift range to form of a sample of nine metal absorbers at z = 6.03 to 6.49. We identify eight galaxies within 1000 km s$^{-1}$ and 350 kpc around the absorbing gas from the ASPIRE spectroscopic data, with their redshifts secured by [OIII]($λλ$4959, 5007) doublets and H$β$ emission lines. Our spectral energy distribution fitting indicates that the absorbing galaxies have stellar mass ranging from 10$^{7.2}$ to 10$^{8.8}M_{\odot}$ and metallicity between 0.02 and 0.4 solar. Notably, the z = 6.2713 system in the J0305-3150 field resides in a galaxy overdensity region, which contains two (tentatively) merging galaxies within 350 kpc and seven galaxies within 1 Mpc. We measure the relative abundances of $α$ elements to iron ([$α$/Fe]) and find that the CGM gas in the most overdense region exhibits a lower [$α$/Fe] ratio. Our modeling of the galaxy's chemical abundance favors a top-heavy stellar initial mass function, and hints that we may be witnessing the contribution of the first generation Population III stars to the CGM at the end of reionization epoch. △ Less

Submitted 31 January, 2024; originally announced February 2024.

Comments: 21 pages, 4 figures in the main text. Accepted for publication in ApJL

AGN X-ray luminosity function and absorption function in the Early Universe ($3\leq z \leq 6$)

Authors: E. Pouliasis, A. Ruiz, I. Georgantopoulos, F. Vito, R. Gilli, C. Vignali, Y. Ueda, E. Koulouridis, M. Akiyama, S. Marchesi, B. Laloux, T. Nagao, S. Paltani, M. Pierre, Y. Toba, M. Habouzit, B. Vijarnwannaluk, C. Garrel

Abstract: The XLF of AGN offers a robust tool to study the evolution and the growth of SMBHs over cosmic time. Owing to the limited area probed by X-ray surveys, optical surveys are routinely used to probe the accretion in the high redshift Universe $z\geq 3$. However, optical surveys may be incomplete because they are strongly affected by dust redenning. In this work, we derive the XLF and its evolution at… ▽ More The XLF of AGN offers a robust tool to study the evolution and the growth of SMBHs over cosmic time. Owing to the limited area probed by X-ray surveys, optical surveys are routinely used to probe the accretion in the high redshift Universe $z\geq 3$. However, optical surveys may be incomplete because they are strongly affected by dust redenning. In this work, we derive the XLF and its evolution at high redshifts using a large sample of AGNs selected in different fields with various areas and depths covering a wide range of luminosities. Additionally, we put the tightest yet constraints on the absorption function in this redshift regime. In particular, we use more than 600 soft X-ray selected high-z sources in the Chandra Deep fields, the Chandra COSMOS Legacy survey and the XMM-XXL northern field. We derive the X-ray spectral properties for all sources via spectral fitting, using a consistent technique and model. For modeling the parametric form of the XLF and the absorption function, we use a Bayesian methodology allowing us to correctly propagate the uncertainties for the observed X-ray properties of our sources and also the absorption effects. The evolution of XLF is in agreement with a pure density evolution model similar to what is witnessed at optical wavelengths, although a luminosity dependent density evolution model cannot be securely ruled out. A large fraction ($60\%)$ of our sources are absorbed by column densities of $\rm N_H \geq 10^{23} cm^{-2} $, while $17$\% of the sources are CTK. Our results favor a scenario where both the ISM of the host and the AGN torus contribute to the obscuration. The derived BHAD is in agreement with the simulations, if one takes into account that the X-ray AGN are hosted by massive galaxies, while it differs from the one derived using JWST data. The latter could be due to the differences in the AGN and host-galaxy properties. △ Less

Submitted 24 January, 2024; originally announced January 2024.

Comments: Accepted for publication in A&A. 25 pages, 15 figures (+4 in Appendix), 2 tables (+1 in Appendix)

Journal ref: A&A 685, A97 (2024)

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

Surveying the onset and evolution of supermassive black holes at high-z with AXIS

Authors: Nico Cappelluti, Adi Foord, Stefano Marchesi, Fabio Pacucci, Angelo Ricarte, Melanie Habouzit, Fabio Vito, Meredith Powell, Michael Koss, Richard Mushotzky, the AXIS AGN-SWG

Abstract: The nature and origin of supermassive black holes (SMBHs) remain an open matter of debate within the scientific community. While various theoretical scenarios have been proposed, each with specific observational signatures, the lack of sufficiently sensitive X-ray observations hinders the progress of observational tests. In this white paper, we present how AXIS will contribute to solving this issu… ▽ More The nature and origin of supermassive black holes (SMBHs) remain an open matter of debate within the scientific community. While various theoretical scenarios have been proposed, each with specific observational signatures, the lack of sufficiently sensitive X-ray observations hinders the progress of observational tests. In this white paper, we present how AXIS will contribute to solving this issue. With an angular resolution of 1.5$^{\prime\prime}$ on-axis and minimal off-axis degradation, we have designed a deep survey capable of reaching flux limits in the [0.5-2] keV range of approximately 2$\times$10$^{-18}$ \fcgs~ over an area of 0.13 deg$^2$ in approximately 7 million seconds (7 Ms). Furthermore, we have planned an intermediate depth survey covering approximately 2 deg$^2$ and reaching flux limits of about 2$\times$10$^{-17}$ \fcgs ~ in order to detect a significant number of SMBHs with X-ray luminosities (L$_X$) of approximately 10$^{42}$ \lx up to z$\sim$10. These observations will enable AXIS to detect SMBHs with masses smaller than 10$^5$ \ms, assuming Eddington-limited accretion and a typical bolometric correction for Type II AGN. AXIS will provide valuable information on the seeding and population synthesis models of SMBH, allowing for more accurate constraints on their initial mass function (IMF) and accretion history from z$\sim$0-10. To accomplish this, AXIS will leverage the unique synergy of survey telescopes such as JWST, Roman, Euclid, LSST, and the new generation of 30m class telescopes. These instruments will provide optical identification and redshift measurements, while AXIS will discover the smoking gun of nuclear activity, particularly in the case of highly obscured AGN or peculiar UV spectra as predicted and recently observed in the early Universe. △ Less

Submitted 13 November, 2023; originally announced November 2023.

Comments: This White Paper is part of a series commissioned for the AXIS Probe Concept Mission. This White Paper is part of a series commissioned for the AXIS Probe Concept Mission; additional AXIS White Papers can be found at http://axis.astro.umd.edu with a mission overview at arXiv:2311.00780

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)

Overview of the Advanced X-ray Imaging Satellite (AXIS)

Authors: Christopher S. Reynolds, Erin A. Kara, Richard F. Mushotzky, Andrew Ptak, Michael J. Koss, Brian J. Williams, Steven W. Allen, Franz E. Bauer, Marshall Bautz, Arash Bodaghee, Kevin B. Burdge, Nico Cappelluti, Brad Cenko, George Chartas, Kai-Wing Chan, Lía Corrales, Tansu Daylan, Abraham D. Falcone, Adi Foord, Catherine E. Grant, Mélanie Habouzit, Daryl Haggard, Sven Herrmann, Edmund Hodges-Kluck, Oleg Kargaltsev , et al. (18 additional authors not shown)

Abstract: The Advanced X-ray Imaging Satellite (AXIS) is a Probe-class concept that will build on the legacy of the Chandra X-ray Observatory by providing low-background, arcsecond-resolution imaging in the 0.3-10 keV band across a 450 arcminute$^2$ field of view, with an order of magnitude improvement in sensitivity. AXIS utilizes breakthroughs in the construction of lightweight segmented X-ray optics usin… ▽ More The Advanced X-ray Imaging Satellite (AXIS) is a Probe-class concept that will build on the legacy of the Chandra X-ray Observatory by providing low-background, arcsecond-resolution imaging in the 0.3-10 keV band across a 450 arcminute$^2$ field of view, with an order of magnitude improvement in sensitivity. AXIS utilizes breakthroughs in the construction of lightweight segmented X-ray optics using single-crystal silicon, and developments in the fabrication of large-format, small-pixel, high readout rate CCD detectors with good spectral resolution, allowing a robust and cost-effective design. Further, AXIS will be responsive to target-of-opportunity alerts and, with onboard transient detection, will be a powerful facility for studying the time-varying X-ray universe, following on from the legacy of the Neil Gehrels (Swift) X-ray observatory that revolutionized studies of the transient X-ray Universe. In this paper, we present an overview of AXIS, highlighting the prime science objectives driving the AXIS concept and how the observatory design will achieve these objectives. △ Less

Submitted 1 November, 2023; originally announced November 2023.

Comments: Published in Proceedings of SPIE Optics & Photonics 2023, San Diego

Discovery and Characterization of Galactic-scale Dual Supermassive Black Holes Across Cosmic Time

Authors: Yue Shen, J. Andrew Casey-Clyde, Yu-Ching Chen, Arran Gross, Melanie Habouzit, Hsiang-Chih Hwang, Yuzo Ishikawa, Jun-Yao Li, Xin Liu, Chiara M. F. Mingarelli, D. Porquet, Aaron Stemo, Ming-Yang Zhuang

Abstract: The hierarchical structure formation paradigm predicts the formation of pairs of supermassive black holes in merging galaxies. When both (or one) members of the SMBH pair are unobscured AGNs, the system can be identified as a dual (or offset) AGN. Quantifying the abundance of these AGN pairs as functions of separation, redshift and host properties is crucial to understanding SMBH formation and AGN… ▽ More The hierarchical structure formation paradigm predicts the formation of pairs of supermassive black holes in merging galaxies. When both (or one) members of the SMBH pair are unobscured AGNs, the system can be identified as a dual (or offset) AGN. Quantifying the abundance of these AGN pairs as functions of separation, redshift and host properties is crucial to understanding SMBH formation and AGN fueling in the broad context of galaxy formation. The High Latitude Wide Area Survey with Roman, with its unprecedented combination of sensitivity, spatial resolution, area and NIR wavelength coverage, will revolutionize the study of galactic-scale environments of SMBH pairs. This white paper summarizes the science opportunities and technical requirements on the discovery and characterization of SMBH pairs down to galactic scales (i.e., less than tens of kpc) over broad ranges of redshift (1<z<7) and luminosity (Lbol>1E42 erg/s). △ Less

Submitted 27 June, 2023; originally announced June 2023.

Comments: Roman Core Community Survey White Paper, focusing on the High Latitude Wide Area Survey

Massive Black Hole Binaries as LISA Precursors in the Roman High Latitude Time Domain Survey

Authors: Zoltán Haiman, Chengcheng Xin, Tamara Bogdanović, Pau Amaro Seoane, Matteo Bonetti, J. Andrew Casey-Clyde, Maria Charisi, Monica Colpi, Jordy Davelaar, Alessandra De Rosa, Daniel J. D'Orazio, Kate Futrowsky, Poshak Gandhi, Alister W. Graham, Jenny E. Greene, Melanie Habouzit, Daryl Haggard, Kelly Holley-Bockelmann, Xin Liu, Alberto Mangiagli, Alessandra Mastrobuono-Battisti, Sean McGee, Chiara M. F. Mingarelli, Rodrigo Nemmen, Antonella Palmese , et al. (5 additional authors not shown)

Abstract: With its capacity to observe $\sim 10^{5-6}$ faint active galactic nuclei (AGN) out to redshift $z\approx 6$, Roman is poised to reveal a population of $10^{4-6}\, {\rm M_\odot}$ black holes during an epoch of vigorous galaxy assembly. By measuring the light curves of a subset of these AGN and looking for periodicity, Roman can identify several hundred massive black hole binaries (MBHBs) with 5-12… ▽ More With its capacity to observe $\sim 10^{5-6}$ faint active galactic nuclei (AGN) out to redshift $z\approx 6$, Roman is poised to reveal a population of $10^{4-6}\, {\rm M_\odot}$ black holes during an epoch of vigorous galaxy assembly. By measuring the light curves of a subset of these AGN and looking for periodicity, Roman can identify several hundred massive black hole binaries (MBHBs) with 5-12 day orbital periods, which emit copious gravitational radiation and will inevitably merge on timescales of $10^{3-5}$ years. During the last few months of their merger, such binaries are observable with the Laser Interferometer Space Antenna (LISA), a joint ESA/NASA gravitational wave mission set to launch in the mid-2030s. Roman can thus find LISA precursors, provide uniquely robust constraints on the LISA source population, help identify the host galaxies of LISA mergers, and unlock the potential of multi-messenger astrophysics with massive black hole binaries. △ Less

Submitted 26 June, 2023; originally announced June 2023.

Comments: White Paper for the Nancy Grace Roman Space Telescope's Core Community Surveys (https://roman.gsfc.nasa.gov/science/ccs_white_papers.html)

Mass-redshift dependency of Supermassive Black Hole Binaries for the Gravitational Wave Background

Authors: Musfar Muhamed Kozhikkal, Siyuan Chen, Gilles Theureau, Melanie Habouzit, Alberto Sesana

Abstract: Studying how the black hole (BH) - (galaxy) bulge mass relation evolves with redshift provides valuable insights into the co-evolution of supermassive black holes and their host galaxies. However, obtaining accurate measurement of BH masses is challenging due to the bias towards the most massive and luminous galaxies. Instead we focus on the BH and bulge masses as they vary with redshift using the… ▽ More Studying how the black hole (BH) - (galaxy) bulge mass relation evolves with redshift provides valuable insights into the co-evolution of supermassive black holes and their host galaxies. However, obtaining accurate measurement of BH masses is challenging due to the bias towards the most massive and luminous galaxies. Instead we focus on the BH and bulge masses as they vary with redshift using the EAGLE, Illustris, TNG100, TNG300, Horizon-AGN and SIMBA large-scale cosmological simulations. We use an analytical astrophysical model with galaxy stellar mass function, pair fraction, merger timescale and BH-bulge mass relation extended to include redshift evolution. The model can predict the intensity of the gravitational wave background (GWB) produced by a population of supermassive black hole binary (SMBHB) as a function of the frequency. This allows us to compare the predictions of this model with the constraints of Pulsar Timing Array observations. Here, we employ Bayesian analysis for the parameter inference. We find that all six simulations are consistent $\leq 3.5σ$ with a range of simulated GWB spectra. By fixing the BH-bulge mass parameters to the simulations we analyze the changes in the constraints on the other astrophysical parameters. Furthermore, we also examine the variation in SMBHB merger rate with mass and redshift between these large-scale simulations. △ Less

Submitted 12 March, 2024; v1 submitted 29 May, 2023; originally announced May 2023.

Comments: 20 pages, 16 figures, 3 tables, 1 appendix figure and 1 appendix table

A SPectroscopic survey of biased halos In the Reionization Era (ASPIRE): JWST Reveals a Filamentary Structure around a z=6.61 Quasar

Authors: Feige Wang, Jinyi Yang, Joseph F. Hennawi, Xiaohui Fan, Fengwu Sun, Jaclyn B. Champagne, Tiago Costa, Melanie Habouzit, Ryan Endsley, Zihao Li, Xiaojing Lin, Romain A. Meyer, Jan-Torge Schindler, Yunjing Wu, Eduardo Bañados, Aaron J. Barth, Aklant K. Bhowmick, Rebekka Bieri, Laura Blecha, Sarah Bosman, Zheng Cai, Luis Colina, Thomas Connor, Frederick B. Davies, Roberto Decarli , et al. (34 additional authors not shown)

Abstract: We present the first results from the JWST ASPIRE program (A SPectroscopic survey of biased halos In the Reionization Era). This program represents an imaging and spectroscopic survey of 25 reionization-era quasars and their environments by utilizing the unprecedented capabilities of NIRCam Wide Field Slitless Spectroscopy (WFSS) mode. ASPIRE will deliver the largest ($\sim280~{\rm arcmin}^2$) gal… ▽ More We present the first results from the JWST ASPIRE program (A SPectroscopic survey of biased halos In the Reionization Era). This program represents an imaging and spectroscopic survey of 25 reionization-era quasars and their environments by utilizing the unprecedented capabilities of NIRCam Wide Field Slitless Spectroscopy (WFSS) mode. ASPIRE will deliver the largest ($\sim280~{\rm arcmin}^2$) galaxy redshift survey at 3-4 $μ$m among JWST Cycle-1 programs and provide extensive legacy values for studying the formation of the earliest supermassive black holes (SMBHs), the assembly of galaxies, early metal enrichment, and cosmic reionization. In this first ASPIRE paper, we report the discovery of a filamentary structure traced by the luminous quasar J0305-3150 and ten [OIII] emitters at $z=6.6$. This structure has a 3D galaxy overdensity of $δ_{\rm gal}=12.6$ over 637 cMpc$^3$, one of the most overdense structures known in the early universe, and could eventually evolve into a massive galaxy cluster. Together with existing VLT/MUSE and ALMA observations of this field, our JWST observations reveal that J0305-3150 traces a complex environment where both UV-bright and dusty galaxies are present, and indicate that the early evolution of galaxies around the quasar is not simultaneous. In addition, we discovered 31 [OIII] emitters in this field at other redshifts, $5.3<z<6.7$, with half of them situated at $z\sim5.4$ and $z\sim6.2$. This indicates that star-forming galaxies, such as [OIII] emitters, are generally clustered at high redshifts. These discoveries demonstrate the unparalleled redshift survey capabilities of NIRCam WFSS and the potential of the full ASPIRE survey dataset. △ Less

Submitted 19 April, 2023; originally announced April 2023.

Comments: accepted for publication in ApJL

A SPectroscopic survey of biased halos In the Reionization Era (ASPIRE): A First Look at the Rest-frame Optical Spectra of $z > 6.5$ Quasars Using JWST

Authors: Jinyi Yang, Feige Wang, Xiaohui Fan, Joseph F. Hennawi, Aaron J. Barth, Eduardo Bañados, Fengwu Sun, Weizhe Liu, Zheng Cai, Linhua Jiang, Zihao Li, Masafusa Onoue, Jan-Torge Schindler, Yue Shen, Yunjing Wu, Aklant K. Bhowmick, Rebekka Bieri, Laura Blecha, Sarah Bosman, Jaclyn B. Champagne, Luis Colina, Thomas Connor, Tiago Costa, Frederick B. Davies, Roberto Decarli , et al. (31 additional authors not shown)

Abstract: Studies of rest-frame optical emission in quasars at $z>6$ have historically been limited by the wavelengths accessible by ground-based telescopes. The James Webb Space Telescope (JWST) now offers the opportunity to probe this emission deep into the reionization epoch. We report the observations of eight quasars at $z>6.5$ using the JWST/NIRCam Wide Field Slitless Spectroscopy, as a part of the ''… ▽ More Studies of rest-frame optical emission in quasars at $z>6$ have historically been limited by the wavelengths accessible by ground-based telescopes. The James Webb Space Telescope (JWST) now offers the opportunity to probe this emission deep into the reionization epoch. We report the observations of eight quasars at $z>6.5$ using the JWST/NIRCam Wide Field Slitless Spectroscopy, as a part of the ''A SPectroscopic survey of biased halos In the Reionization Era (ASPIRE)" program. Our JWST spectra cover the quasars' emission between rest frame $\sim$ 4100 and 5100 Å. The profiles of these quasars' broad H$β$ emission lines span a FWHM from 3000 to 6000 $\rm{km~s^{-1}}$. The H$β$-based virial black hole (BH) masses, ranging from 0.6 to 2.1 billion solar masses, are generally consistent with their MgII-based BH masses. The new measurements based on the more reliable H$β$ tracer thus confirm the existence of billion solar-mass BHs in the reionization epoch. In the observed [OIII] $λλ$4960,5008 doublets of these luminous quasars, broad components are more common than narrow core components ($\le~1200~\rm{km~s^{-1}}$), and only one quasar shows stronger narrow components than broad. Two quasars exhibit significantly broad and blueshifted [OIII] emission, thought to trace galactic-scale outflows, with median velocities of $-610~\rm{km~s^{-1}}$ and $-1430~\rm{km~s^{-1}}$ relative to the [CII] $158\,μ$m line. All eight quasars show strong optical FeII emission, and follow the Eigenvector 1 relations defined by low-redshift quasars. The entire ASPIRE program will eventually cover 25 quasars and provide a statistical sample for the studies of the BHs and quasar spectral properties. △ Less

Submitted 19 April, 2023; originally announced April 2023.

Comments: 13 pages, 4 figures, accepted for publication in ApJL

What if young z>9 JWST galaxies hosted massive black holes?

Authors: Marta Volonteri, Melanie Habouzit, Monica Colpi

Abstract: JWST is discovering star forming `candidate' galaxies with photometric redshifts $z>9$ and little attenuation. We model presumptive massive black holes (MBHs) in such galaxies and find that their unobscured emission is fainter than the galaxy starlight in JWST filters, and difficult to be detected via color-color selection, and X-ray and radio observations. Only MBHs overmassive relative to expect… ▽ More JWST is discovering star forming `candidate' galaxies with photometric redshifts $z>9$ and little attenuation. We model presumptive massive black holes (MBHs) in such galaxies and find that their unobscured emission is fainter than the galaxy starlight in JWST filters, and difficult to be detected via color-color selection, and X-ray and radio observations. Only MBHs overmassive relative to expected galaxy scaling relations, accreting at high Eddington rates, would be detectable. Their discovery would point to the presence of heavy MBH seeds, but care is needed to exclude the existence of lighter seeds as only overmassive MBHs are detectable in this type of galaxies. Conversely, if no overmassive MBHs are hosted in these galaxies, either there are no heavy seeds or they are rare. The most massive/highest redshift candidate galaxies can attain stellar masses in excess of 5e10 Msun by z~6 if they grow along the SFR-mass sequence, and can nurse a MBH growing from ~1e5 Msun up to >3e7 Msun by z~6, to become hosts of some z>6 quasars. Candidate galaxies of log(M_{gal}/Msun)~8 can not grow their putative seeds fast, unless seeds are >1e6 Msun. The number density of the JWST candidate galaxies far outnumbers that of the highest-z quasar hosts and this allows for about only 1 bright z~6-7 quasar every 1000 of these galaxies. △ Less

Submitted 3 March, 2023; v1 submitted 9 December, 2022; originally announced December 2022.

Comments: Published in MNRAS

Detection of stellar light from quasar host galaxies at redshifts above 6

Authors: Xuheng Ding, Masafusa Onoue, John D. Silverman, Yoshiki Matsuoka, Takuma Izumi, Michael A. Strauss, Knud Jahnke, Camryn L. Phillips, Junyao Li, Marta Volonteri, Zoltan Haiman, Irham Taufik Andika, Kentaro Aoki, Shunsuke Baba, Rebekka Bieri, Sarah E. I. Bosman, Connor Bottrell, Anna-Christina Eilers, Seiji Fujimoto, Melanie Habouzit, Masatoshi Imanishi, Kohei Inayoshi, Kazushi Iwasawa, Nobunari Kashikawa, Toshihiro Kawaguchi , et al. (19 additional authors not shown)

Abstract: The detection of starlight from the host galaxies of quasars during the reionization epoch ($z>6$) has been elusive, even with deep HST observations. The current highest redshift quasar host detected, at $z=4.5$, required the magnifying effect of a foreground lensing galaxy. Low-luminosity quasars from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) mitigate the challenge of detecting the… ▽ More The detection of starlight from the host galaxies of quasars during the reionization epoch ($z>6$) has been elusive, even with deep HST observations. The current highest redshift quasar host detected, at $z=4.5$, required the magnifying effect of a foreground lensing galaxy. Low-luminosity quasars from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) mitigate the challenge of detecting their underlying, previously-undetected host galaxies. Here we report rest-frame optical images and spectroscopy of two HSC-SSP quasars at $z>6$ with JWST. Using NIRCam imaging at 3.6$μ$m and 1.5$μ$m and subtracting the light from the unresolved quasars, we find that the host galaxies are massive (stellar masses of $13\times$ and $3.4\times$ $10^{10}$ M$_{\odot}$, respectively), compact, and disk-like. NIRSpec medium-resolution spectroscopy shows stellar absorption lines in the more massive quasar, confirming the detection of the host. Velocity-broadened gas in the vicinity of these quasars enables measurements of their black hole masses ($1.4\times 10^9$ and $2.0\times$ $10^{8}$ M$_{\odot}$, respectively). Their location in the black hole mass - stellar mass plane is consistent with the distribution at low redshift, suggesting that the relation between black holes and their host galaxies was already in place less than a billion years after the Big Bang. △ Less

Submitted 23 June, 2023; v1 submitted 25 November, 2022; originally announced November 2022.

Comments: Matched to the published Nature version of the article. 27 pages, 4 main figures, 1 table, 6 supplementary figures, 2 supplementary table

Astrophysics with the Laser Interferometer Space Antenna

Authors: Pau Amaro Seoane, Jeff Andrews, Manuel Arca Sedda, Abbas Askar, Quentin Baghi, Razvan Balasov, Imre Bartos, Simone S. Bavera, Jillian Bellovary, Christopher P. L. Berry, Emanuele Berti, Stefano Bianchi, Laura Blecha, Stephane Blondin, Tamara Bogdanović, Samuel Boissier, Matteo Bonetti, Silvia Bonoli, Elisa Bortolas, Katelyn Breivik, Pedro R. Capelo, Laurentiu Caramete, Federico Cattorini, Maria Charisi, Sylvain Chaty , et al. (134 additional authors not shown)

Abstract: The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery… ▽ More The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery potential of LISA. The next decade is crucial to prepare the astrophysical community for LISA's first observations. This review outlines the extensive landscape of astrophysical theory, numerical simulations, and astronomical observations that are instrumental for modeling and interpreting the upcoming LISA datastream. To this aim, the current knowledge in three main source classes for LISA is reviewed; ultracompact stellar-mass binaries, massive black hole binaries, and extreme or intermediate mass ratio inspirals. The relevant astrophysical processes and the established modeling techniques are summarized. Likewise, open issues and gaps in our understanding of these sources are highlighted, along with an indication of how LISA could help making progress in the different areas. New research avenues that LISA itself, or its joint exploitation with upcoming studies in the electromagnetic domain, will enable, are also illustrated. Improvements in modeling and analysis approaches, such as the combination of numerical simulations and modern data science techniques, are discussed. This review is intended to be a starting point for using LISA as a new discovery tool for understanding our Universe. △ Less

Submitted 25 May, 2023; v1 submitted 11 March, 2022; originally announced March 2022.

Journal ref: Living Reviews in Relativity, Volume 26, Article number: 2 (2023)

Co-evolution of massive black holes and their host galaxies at high redshift: discrepancies from six cosmological simulations and the key role of JWST

Authors: Melanie Habouzit, Masafusa Onoue, Eduardo Banados, Marcel Neeleman, Daniel Angles-Alcazar, Fabian Walter, Annalisa Pillepich, Romeel Dave, Knud Jahnke, Yohan Dubois

Abstract: The James Webb Space Telescope will have the power to characterize high-redshift quasars at z>6 with an unprecedented depth and spatial resolution. While the brightest quasars at such redshift (i.e., with bolometric luminosity L_bol> 10^46 erg/s) provide us with key information on the most extreme objects in the Universe, measuring the black hole (BH) mass and Eddington ratios of fainter quasars w… ▽ More The James Webb Space Telescope will have the power to characterize high-redshift quasars at z>6 with an unprecedented depth and spatial resolution. While the brightest quasars at such redshift (i.e., with bolometric luminosity L_bol> 10^46 erg/s) provide us with key information on the most extreme objects in the Universe, measuring the black hole (BH) mass and Eddington ratios of fainter quasars with L_bol= 10^45-10^46 erg/s opens a path to understand the build-up of more normal BHs at z>6. In this paper, we show that the Illustris, TNG100, TNG300, Horizon-AGN, EAGLE, and SIMBA large-scale cosmological simulations do not agree on whether BHs at z>4 are overmassive or undermassive at fixed galaxy stellar mass with respect to the M_BH-M_star scaling relation at z=0 (BH mass offsets). Our conclusions are unchanged when using the local scaling relation produced by each simulation or empirical relations. We find that the BH mass offsets of the simulated faint quasar population at z>4, unlike those of bright quasars, represent the BH mass offsets of the entire BH population, for all the simulations. Thus, a population of faint quasars with L_bol= 10^45-10^46 erg/s observed by JWST can provide key constraints on the assembly of BHs at high redshift. Moreover, this will help constraining the high-redshift regime of cosmological simulations, including BH seeding, early growth, and co-evolution with the host galaxies. Our results also motivate the need for simulations of larger cosmological volumes down to z=6, with the same diversity of sub-grid physics, in order to gain statistics on the most extreme objects at high redshift. △ Less

Submitted 24 January, 2022; originally announced January 2022.

Comments: published in MNRAS, 19 pages, 8 figures, key figures: Fig. 3, Fig.5, and Fig. 7

The black hole population in low-mass galaxies in large-scale cosmological simulations

Authors: Houda Haidar, Melanie Habouzit, Marta Volonteri, Mar Mezcua, Jenny Greene, Nadine Neumayer, Daniel Angles-Alcazar, Ignacio Martin-Navarro, Nils Hoyer, Yohan Dubois, Romeel Dave

Abstract: Recent systematic searches for massive black holes (BHs) in local dwarf galaxies led to the discovery of a population of faint Active Galactic Nuclei (AGN). We investigate the agreement of the BH and AGN populations in the Illustris, TNG, Horizon-AGN, EAGLE, and SIMBA simulations with current observational constraints in low-mass galaxies. We find that some of these simulations produce BHs that ar… ▽ More Recent systematic searches for massive black holes (BHs) in local dwarf galaxies led to the discovery of a population of faint Active Galactic Nuclei (AGN). We investigate the agreement of the BH and AGN populations in the Illustris, TNG, Horizon-AGN, EAGLE, and SIMBA simulations with current observational constraints in low-mass galaxies. We find that some of these simulations produce BHs that are too massive, and that the BH occupation fraction at z=0 is not inherited from the simulation seeding modeling. The ability of BHs and their host galaxies to power an AGN depends on BH and galaxy subgrid modeling. The fraction of AGN in low-mass galaxies is not used to calibrate the simulations, and thus can be used to differentiate galaxy formation models. AGN fractions at z=0 span two orders of magnitude at fixed galaxy stellar mass in simulations, similarly to observational constraints, but uncertainties and degeneracies affect both observations and simulations. The agreement is difficult to interpret due to differences in the masses of simulated and observed BHs, BH occupation fraction affected by numerical choices, and an unknown fraction of obscured AGN. Our work advocates for more thorough comparisons with observations to improve the modeling of cosmological simulations, and our understanding of BH and galaxy physics in the low-mass regime. The mass of BHs, their ability to efficiently accrete gas, and the AGN fraction in low-mass galaxies have important implications for the build-up of the entire BH and galaxy populations with time. △ Less

Submitted 14 June, 2022; v1 submitted 24 January, 2022; originally announced January 2022.

Comments: Accepted in MNRAS, 21 pages, 11 figures, 1 table

Supermassive black holes in cosmological simulations II: the AGN population and predictions for upcoming X-ray missions

Authors: Melanie Habouzit, Rachel S. Somerville, Yuan Li, Shy Genel, James Aird, Daniel Anglés-Alcázar, Romeel Davé, Iskren Y. Georgiev, Stuart McAlpine, Yetli Rosas-Guevara, Yohan Dubois, Dylan Nelson, Eduardo Bañados, Lars Hernquist, Sébastien Peirani, Mark Vogelsberger

Abstract: In large-scale hydrodynamical cosmological simulations, the fate of massive galaxies is mainly dictated by the modeling of feedback from active galactic nuclei (AGN). The amount of energy released by AGN feedback is proportional to the mass that has been accreted onto the BHs, but the exact sub-grid modeling of AGN feedback differs in all simulations. Whilst modern simulations reliably produce pop… ▽ More In large-scale hydrodynamical cosmological simulations, the fate of massive galaxies is mainly dictated by the modeling of feedback from active galactic nuclei (AGN). The amount of energy released by AGN feedback is proportional to the mass that has been accreted onto the BHs, but the exact sub-grid modeling of AGN feedback differs in all simulations. Whilst modern simulations reliably produce populations of quiescent massive galaxies at z<2, it is also crucial to assess the similarities and differences of the responsible AGN populations. Here, we compare the AGN population of the Illustris, TNG100, TNG300, Horizon-AGN, EAGLE, and SIMBA simulations. The AGN luminosity function (LF) varies significantly between simulations. Although in agreement with current observational constraints at z=0, at higher redshift the agreement of the LFs deteriorates with most simulations producing too many AGN of L_{x, 2-10 keV}~10^43-10^44 erg/s. AGN feedback in some simulations prevents the existence of any bright AGN with L_{x, 2-10 keV}>=10^45 erg/s (although this is sensitive to AGN variability), and leads to smaller fractions of AGN in massive galaxies than in the observations at z<=2. We find that all the simulations fail at producing a number density of AGN in good agreement with observational constraints for both luminous (L_{x, 2-10 keV}~10^43-10^45 erg/s) and fainter (L_{x, 2-10 keV}~10^42-10^43 erg/s) AGN, and at both low and high redshift. These differences can aid us in improving future BH and galaxy subgrid modeling in simulations. Upcoming X-ray missions (e.g., Athena, AXIS, and LynX) will bring faint AGN to light and new powerful constraints. After accounting for AGN obscuration, we find that the predicted number density of detectable AGN in future surveys spans at least one order of magnitude across the simulations, at any redshift. △ Less

Submitted 2 November, 2021; originally announced November 2021.

Comments: Published in MNRAS

The origins of massive black holes

Authors: Marta Volonteri, Melanie Habouzit, Monica Colpi

Abstract: Massive black holes (MBHs) inhabit galaxy centers, power luminous quasars and Active Galactic Nuclei (AGN) and shape their cosmic environment with the energy they produce. The origins of MBHs remain a mystery and the recent detection by LIGO/Virgo of an almost 150 solar mass black hole has revitalized the question of whether there is a continuum between "stellar" and "massive" black holes and what… ▽ More Massive black holes (MBHs) inhabit galaxy centers, power luminous quasars and Active Galactic Nuclei (AGN) and shape their cosmic environment with the energy they produce. The origins of MBHs remain a mystery and the recent detection by LIGO/Virgo of an almost 150 solar mass black hole has revitalized the question of whether there is a continuum between "stellar" and "massive" black holes and what the seeds of MBHs are. Seeds could have formed in the first galaxies, or could be also related to the collapse of horizon-sized regions in the early Universe. Understanding the origins of MBHs straddles fundamental physics, cosmology and astrophysics and it bridges the fields of gravitational wave physics and traditional astronomy. With several facilities in the next 10-15 years we foresee the possibility of discovering MBHs' avenues of formation. In this article we link three main topics: the channels of black hole seed formation, the journey from seeds to massive black holes, the diagnostics on the origins of MBHs. We highlight and critically discuss current unsolved problems and touch on recent developments that stirred the community. △ Less

Submitted 15 April, 2022; v1 submitted 19 October, 2021; originally announced October 2021.

Comments: This is a post-peer-review, pre-copyedit version of an article published online in Nature Reviews Physics on September 20 2021. The final authenticated version is available online at https://doi.org/10.1038/s42254-021-00364-9. SharedIt link: https://rdcu.be/cx5Dq

Quasar clustering at redshift 6

Authors: J. Greiner, J. Bolmer, R. M. Yates, M. Habouzit, E. Banados, P. M. J. Afonso, P. Schady

Abstract: Large-scale surveys over the last years have revealed about 300 QSOs at redshift above 6. Follow-up observations identified surprising properties, such as the very high black hole (BH) masses, spatial correlations with surrounding cold gas of the host galaxy, or high CIV-MgII velocity shifts. In particular, the discovery of luminous high-redshift quasars suggests that at least some black holes lik… ▽ More Large-scale surveys over the last years have revealed about 300 QSOs at redshift above 6. Follow-up observations identified surprising properties, such as the very high black hole (BH) masses, spatial correlations with surrounding cold gas of the host galaxy, or high CIV-MgII velocity shifts. In particular, the discovery of luminous high-redshift quasars suggests that at least some black holes likely have large masses at birth and grow efficiently. We aim at quantifying quasar pairs at high redshift for a large sample of objects. This provides a new key constraint on a combination of parameters related to the origin and assembly for the most massive black holes: BH formation efficiency and clustering, growth efficiency and relative contribution of BH mergers. We observed 116 spectroscopically confirmed QSOs around redshift 6 with the simultaneous 7-channel imager GROND in order to search for companions. Applying identical colour-colour cuts as for those which led to the spectroscopically confirmed QSO, we perform LePHARE fits to the 26 best QSO pair candidates, and obtained spectroscopic observations for 11 of those. e do not find any QSO pair with a companion brighter than M1450(AB) < -26 mag within our 0.1-3.3 h^-1 cMpc search radius, in contrast to the serendipitous findings in the redshift range 4--5. However, a low fraction of such pairs at this luminosity and redshift is consistent with indications from present-day cosmological-scale galaxy evolution models. In turn, the incidence of L- and T-type brown dwarfs which occupy a similar colour space as z ~ 6 QSOs, is higher than expected, by a factor of 5 and 20, respectively. △ Less

Submitted 20 July, 2021; originally announced July 2021.

Comments: 20 pages, 18 figures, A&A (accepted)

Journal ref: A&A 654, A79 (2021)

Bringing faint active galactic nuclei (AGNs) to light: a view from large-scale cosmological simulations

Authors: Adrian P. Schirra, Melanie Habouzit, Ralf S. Klessen, Francesca Fornasini, Dylan Nelson, Annalisa Pillepich, Daniel Anglés-Alcázar, Romeel Davé, Francesca Civano

Abstract: The sensitivity of X-ray facilities and our ability to detect fainter active galactic nuclei (AGNs) will increase with the upcoming Athena mission and the AXIS and Lynx concept missions, thus improving our understanding of supermassive black holes (BHs) in a luminosity regime that can be dominated by X-ray binaries. We analyze the population of faint AGN (L_x (2-10 keV) < 10^42 erg/s) in the Illus… ▽ More The sensitivity of X-ray facilities and our ability to detect fainter active galactic nuclei (AGNs) will increase with the upcoming Athena mission and the AXIS and Lynx concept missions, thus improving our understanding of supermassive black holes (BHs) in a luminosity regime that can be dominated by X-ray binaries. We analyze the population of faint AGN (L_x (2-10 keV) < 10^42 erg/s) in the Illustris, TNG100, EAGLE, and SIMBA cosmological simulations, and find that the properties of their host galaxies vary from one simulation to another. In Illustris and EAGLE, faint AGN are powered by low-mass BHs located in low-mass star-forming galaxies. In TNG100 and SIMBA, they are mostly associated with more massive BHs in quenched massive galaxies. We model the X-ray binary populations (XRB) of the simulated galaxies, and find that AGN often dominate the galaxy AGN+XRB hard X-ray luminosity at z>2, while XRBs dominate in some simulations at z<2. Whether the AGN or XRB emission dominates in star-forming and quenched galaxies depends on the simulations. These differences in simulations can be used to discriminate between galaxy formation models with future high-resolution X-ray observations. We compare the luminosity of simulated faint AGN host galaxies to observations of stacked galaxies from Chandra. Our comparison indicates that the simulations post-processed with our X-ray modeling tend to overestimate the AGN+XRB X-ray luminosity; luminosity that can be strongly affected by AGN obscuration. Some simulations reveal clear AGN trends as a function of stellar mass (e.g., galaxy luminosity drop in massive galaxies), which are not apparent in the observations. △ Less

Submitted 20 November, 2021; v1 submitted 4 November, 2020; originally announced November 2020.

Comments: Published version - MNRAS

Journal ref: MNRAS, Volume 508, Issue 4, December 2021, Pages 4816-4843

Supermassive black holes in cosmological simulations I: M_BH-M_star relation and black hole mass function

Authors: Melanie Habouzit, Yuan Li, Rachel S. Somerville, Shy Genel, Annalisa Pillepich, Marta Volonteri, Romeel Davé, Yetli Rosas-Guevara, Stuart McAlpine, Sébastien Peirani, Lars Hernquist, Daniel Anglés-Alcázar, Amy Reines, Richard Bower, Yohan Dubois, Dylan Nelson, Christophe Pichon, Mark Vogelsberger

Abstract: The past decade has seen significant progress in understanding galaxy formation and evolution using large-scale cosmological simulations. While these simulations produce galaxies in overall good agreement with observations, they employ different sub-grid models for galaxies and supermassive black holes (BHs). We investigate the impact of the sub-grid models on the BH mass properties of the Illustr… ▽ More The past decade has seen significant progress in understanding galaxy formation and evolution using large-scale cosmological simulations. While these simulations produce galaxies in overall good agreement with observations, they employ different sub-grid models for galaxies and supermassive black holes (BHs). We investigate the impact of the sub-grid models on the BH mass properties of the Illustris, TNG100, TNG300, Horizon-AGN, EAGLE, and SIMBA simulations, focusing on the M_BH-M_star relation and the BH mass function. All simulations predict tight M_BH-M_star relations, and struggle to produce the lowest (M_BH< 10^7.5 Msun) in galaxies of M_star~10^10.5-10^11.5 Msun. While the time evolution of the mean M_BH-M_star relation is mild (<1 dex in BH mass for 0<z<5) for all the simulations, its linearity (shape) and normalization varies from simulation to simulation. The strength of SN feedback has a large impact on the linearity and time evolution for M_star<10^10.5 Msun. We find that the low-mass end is a good discriminant of the simulation models, and highlights the need for new observational constraints. At the high-mass end, strong AGN feedback can suppress the time evolution of the relation normalization. Compared with the observations of the local universe, we find an excess of BHs with M_BH>10^9 Msun in most of the simulations. The BH mass function is dominated by efficiently accreting BHs (log10 f_Edd >-2$) at high redshifts, and transitions progressively from the high-mass to the low-mass end to be governed by inactive BHs. The transition time and the contribution of active BHs are different among the simulations, and can be used to evaluate models against observations. △ Less

Submitted 1 November, 2021; v1 submitted 17 June, 2020; originally announced June 2020.

Comments: Published version - MNRAS

Properties of simulated galaxies and supermassive black holes in cosmic voids

Authors: Melanie Habouzit, Alice Pisani, Andy Goulding, Yohan Dubois, Rachel S. Somerville, Jenny E. Greene

Abstract: Cosmic voids, the under-dense regions of the cosmic web, are widely used to constrain cosmology. Voids contain few, isolated galaxies, presumably expected to be less evolved and preserving memory of the pristine Universe. We use the cosmological hydrodynamical simulation Horizon-AGN coupled to the void finder {\sc \texttt{VIDE}} to investigate properties of galaxies in voids at z=0. We find that,… ▽ More Cosmic voids, the under-dense regions of the cosmic web, are widely used to constrain cosmology. Voids contain few, isolated galaxies, presumably expected to be less evolved and preserving memory of the pristine Universe. We use the cosmological hydrodynamical simulation Horizon-AGN coupled to the void finder {\sc \texttt{VIDE}} to investigate properties of galaxies in voids at z=0. We find that, closer to void centers, low-mass galaxies are more common than their massive counterparts. At fixed dark matter halo mass, they have smaller stellar masses than in denser regions. The star formation rate of void galaxies diminishes when approaching void centers, but their sSFR slightly increases, suggesting that void galaxies form stars more efficiently with respect to their stellar mass. We find that this can not only be attributed to the prevalence of low-mass galaxies. The inner region of voids also predominantly host low-mass BHs. However, the BH mass to galaxy mass ratios resemble those of the whole simulation at z=0. Our results suggest that even if the growth channels in cosmic voids are different than in denser environments, voids grow their galaxies and BHs in a similar way. While a large fraction of the BHs have low Eddington ratios, we find that 20\% could be observed as AGN with log10 L=41.5-42.5 erg/s in hard X-ray (2-10 keV). These results pave the way to future work with larger next-generation hydro simulations, aiming to confirm our findings and prepare the application on data from upcoming large surveys such as PFS, Euclid and WFIRST. △ Less

Submitted 21 January, 2020; v1 submitted 13 December, 2019; originally announced December 2019.

Comments: accepted by MNRAS without modification, 24 pages, 17 figures

Decoupling the rotation of stars and gas -- II: the link between black hole activity and MaNGA kinematics in TNG

Authors: Christopher Duckworth, Tjitske K. Starkenburg, Shy Genel, Timothy Davis, Melanie Habouzit, Katarina Kraljic, Rita Tojeiro

Abstract: We study the relationship between supermassive black hole (BH) feedback, BH luminosity and the kinematics of stars and gas for galaxies in IllustrisTNG. We use a sample of galaxies with mock MaNGA observations to identify kinematic misalignment at $z=0$ (difference in rotation of stars and gas), for which we follow the evolutionary history of BH activity and gas properties over the last 8 Gyrs. Mi… ▽ More We study the relationship between supermassive black hole (BH) feedback, BH luminosity and the kinematics of stars and gas for galaxies in IllustrisTNG. We use a sample of galaxies with mock MaNGA observations to identify kinematic misalignment at $z=0$ (difference in rotation of stars and gas), for which we follow the evolutionary history of BH activity and gas properties over the last 8 Gyrs. Misaligned low mass galaxies typically have boosted BH luminosity, BH growth and have had more energy injected into the gas over the last 8 Gyr in comparison to aligned galaxies. These properties likely lead to outflows and gas loss, in agreement with active low mass galaxies in observations. We show that splitting on BH luminosity at $z=0$ produces statistically consistent distributions of kinematic misalignment at $z=0$, however, splitting on the maximum BH luminosity over the last 8 Gyrs does not. While instantaneous correlation at $z=0$ is difficult due to misalignment persisting on longer timescales, the relationship between BH activity and misalignment is clear. High mass quenched galaxies with misalignment typically have similar BH luminosities, show no overall gas loss, and have typically lower gas phase metallicity over the last 8 Gyrs in comparison to those aligned; suggesting external origin. △ Less

Submitted 12 November, 2019; originally announced November 2019.

Comments: 5 pages, 3 Figures, Submitted to MNRAS letters. For first paper in series, see arXiv:1910.10744

Correlations Between Black Holes and Host Galaxies in the Illustris and IllustrisTNG Simulations

Authors: Yuan Li, Melanie Habouzit, Shy Genel, Rachel Somerville, Bryan A. Terrazas, Eric F. Bell, Annalisa Pillepich, Dylan Nelson, Rainer Weinberger, Vicente Rodriguez-Gomez, Chung-Pei Ma, Ruediger Pakmor, Lars Hernquist, Mark Vogelsberger

Abstract: We study black hole - host galaxy correlations, and the relation between the over-massiveness (the distance from the average $M_{BH}-σ$ relation) of super-massive black holes (SMBHs) and star formation histories of their host galaxies in the Illustris and TNG100 simulations. We find that both simulations are able to produce black hole scaling relations in general agreement with observations at… ▽ More We study black hole - host galaxy correlations, and the relation between the over-massiveness (the distance from the average $M_{BH}-σ$ relation) of super-massive black holes (SMBHs) and star formation histories of their host galaxies in the Illustris and TNG100 simulations. We find that both simulations are able to produce black hole scaling relations in general agreement with observations at $z=0$, but with noticeable discrepancies. Both simulations show an offset from the observations for the $M_{BH}-σ$ relation, and the relation between $M_{BH}$ and the Sersic index. The relation between $M_{BH}$ and stellar mass $M_*$ is tighter than the observations, especially for TNG100. For massive galaxies in both simulations, the hosts of over-massive SMBHs (those above the mean $M_{BH}-σ$ relation) tend to have larger Sersic indices and lower baryon conversion efficiency, suggesting a multidimensional link between SMBHs and properties of their hosts. In Illustris, the hosts of over-massive SMBHs have formed earlier and have lower present-day star formation rates, in qualitative agreement with the observations for massive galaxies with $σ>100 \rm km/s$. For low-mass galaxies, such a correlation still holds in Illustris but does not exist in the observed data. For TNG100, the correlation between SMBH over-massiveness and star formation history is much weaker. The hosts of over-massive SMBHs generally have consistently larger star formation rates throughout history. These galaxies have higher stellar mass as well, due to the strong $M_{BH}-M_*$ correlation. Our findings show that simulated SMBH scaling relations and correlations are sensitive to features in the modeling of SMBHs. △ Less

Submitted 30 September, 2019; originally announced October 2019.

Comments: 18 pages, 11 figures, submitted to ApJ

The relationship between black hole mass and galaxy properties: Examining the black hole feedback model in IllustrisTNG

Authors: Bryan A. Terrazas, Eric F. Bell, Annalisa Pillepich, Dylan Nelson, Rachel S. Somerville, Shy Genel, Rainer Weinberger, Mélanie Habouzit, Yuan Li, Lars Hernquist, Mark Vogelsberger

Abstract: Supermassive black hole feedback is thought to be responsible for the lack of star formation, or quiescence, in a significant fraction of galaxies. We explore how observable correlations between the specific star formation rate (sSFR), stellar mass (M$_{\rm{star}}$), and black hole mass (M$_{\rm{BH}}$) are sensitive to the physics of black hole feedback in a galaxy formation model. We use the Illu… ▽ More Supermassive black hole feedback is thought to be responsible for the lack of star formation, or quiescence, in a significant fraction of galaxies. We explore how observable correlations between the specific star formation rate (sSFR), stellar mass (M$_{\rm{star}}$), and black hole mass (M$_{\rm{BH}}$) are sensitive to the physics of black hole feedback in a galaxy formation model. We use the IllustrisTNG simulation suite, specifically the TNG100 simulation and ten model variations that alter the parameters of the black hole model. Focusing on central galaxies at $z = 0$ with M$_{\rm{star}} > 10^{10}$ M$_{\odot}$, we find that the sSFR of galaxies in IllustrisTNG decreases once the energy from black hole kinetic winds at low accretion rates becomes larger than the gravitational binding energy of gas within the galaxy stellar radius. This occurs at a particular M$_{\rm{BH}}$ threshold above which galaxies are found to sharply transition from being mostly star-forming to mostly quiescent. As a result of this behavior, the fraction of quiescent galaxies as a function of M$_{\rm{star}}$ is sensitive to both the normalization of the M$_{\rm{BH}}$-M$_{\rm{star}}$ relation and the M$_{\rm{BH}}$ threshold for quiescence in IllustrisTNG. Finally, we compare these model results to observations of 91 central galaxies with dynamical M$_{\rm{BH}}$ measurements with the caveat that this sample is not representative of the whole galaxy population. While IllustrisTNG reproduces the observed trend that quiescent galaxies host more massive black holes, the observations exhibit a broader scatter in M$_{\rm{BH}}$ at a given M$_{\rm{star}}$ and show a smoother decline in sSFR with M$_{\rm{BH}}$. △ Less

Submitted 6 June, 2019; originally announced June 2019.

Comments: 20 pages, submitted to MNRAS

Cosmic voids: a novel probe to shed light on our Universe

Authors: Alice Pisani, Elena Massara, David N. Spergel, David Alonso, Tessa Baker, Yan-Chuan Cai, Marius Cautun, Christopher Davies, Vasiliy Demchenko, Olivier Doré, Andy Goulding, Mélanie Habouzit, Nico Hamaus, Adam Hawken, Christopher M. Hirata, Shirley Ho, Bhuvnesh Jain, Christina D. Kreisch, Federico Marulli, Nelson Padilla, Giorgia Pollina, Martin Sahlén, Ravi K. Sheth, Rachel Somerville, Istvan Szapudi , et al. (4 additional authors not shown)

Abstract: Cosmic voids, the less dense patches of the Universe, are promising laboratories to extract cosmological information. Thanks to their unique low density character, voids are extremely sensitive to diffuse components such as neutrinos and dark energy, and represent ideal environments to study modifications of gravity, where the effects of such modifications are expected to be more prominent. Robust… ▽ More Cosmic voids, the less dense patches of the Universe, are promising laboratories to extract cosmological information. Thanks to their unique low density character, voids are extremely sensitive to diffuse components such as neutrinos and dark energy, and represent ideal environments to study modifications of gravity, where the effects of such modifications are expected to be more prominent. Robust void-related observables, including for example redshift-space distortions (RSD) and weak lensing around voids, are a promising way to chase and test new physics. Cosmological analysis of the large-scale structure of the Universe predominantly relies on the high density regions. Current and upcoming surveys are designed to optimize the extraction of cosmological information from these zones, but leave voids under-exploited. A dense, large area spectroscopic survey with imaging capabilities is ideal to exploit the power of voids fully. Besides helping illuminate the nature of dark energy, modified gravity, and neutrinos, this survey will give access to a detailed map of under-dense regions, providing an unprecedented opportunity to observe and study a so far under-explored galaxy population. △ Less

Submitted 12 March, 2019; originally announced March 2019.

Comments: Science White Paper for the Astro2020 Decadal Survey, 5 pages text, 14 pages with references and affiliations

The First Luminous Quasars and Their Host Galaxies

Authors: Xiaohui Fan, Aaron Barth, Eduardo Banados, Gisella De Rosa, Roberto Decarli, Anna-Christina Eilers, Emanuele Paolo Farina, Jenny Greene, Melanie Habouzit, Linhua Jiang, Hyunsung D. Jun, Anton Koekemoer, Sangeeta Malhotra, Chiara Mazzucchelli, Fabio Pacucci, James Rhoads, Dominik Riechers, Jane Rigby, Yue Shen, Robert A. Simcoe, Dan Stern, Michael A. Strauss, Tommaso Treu, Bram Venemans, Marianne Vestergaard , et al. (4 additional authors not shown)

Abstract: The discovery of luminous quasars at redshifts up to 7.5 demonstrates the existence of several billion M_sun supermassive black holes (SMBHs) less than a billion years after the Big Bang. They are accompanied by intense star formation in their host galaxies, pinpointing sites of massive galaxy assembly in the early universe, while their absorption spectra reveal an increasing neutral intergalactic… ▽ More The discovery of luminous quasars at redshifts up to 7.5 demonstrates the existence of several billion M_sun supermassive black holes (SMBHs) less than a billion years after the Big Bang. They are accompanied by intense star formation in their host galaxies, pinpointing sites of massive galaxy assembly in the early universe, while their absorption spectra reveal an increasing neutral intergalactic medium (IGM) at the epoch of reionization. Extrapolating from the rapid evolution of the quasar density at z=5-7, we expect that there is only one luminous quasar powered by a billion M_sun SMBH in the entire observable universe at z~9. In the next decade, new wide-field, deep near-infrared (NIR) sky surveys will push the redshift frontier to the first luminous quasars at z~9-10; the combination with new deep X-ray surveys will probe fainter quasar populations that trace earlier phases of SMBH growth. The identification of these record-breaking quasars, and the measurements of their BH masses and accretion properties require sensitive spectroscopic observations with next generation of ground-based and space telescopes at NIR wavelengths. High-resolution integral-field spectroscopy at NIR, and observations at millimeter and radio wavelengths, will together provide a panchromatic view of the quasar host galaxies and their galactic environment at cosmic dawn, connecting SMBH growth with the rise of the earliest massive galaxies. Systematic surveys and multiwavelength follow-up observations of the earliest luminous quasars will strongly constrain the seeding and growth of the first SMBHs in the universe, and provide the best lines of sight to study the history of reionization. △ Less

Submitted 10 March, 2019; originally announced March 2019.

Comments: 7 pages, 3 figures, Science white paper submitted to the Astro2020 Decadal Survey

The diverse galaxy counts in the environment of high-redshift massive black holes in Horizon-AGN

Authors: Mélanie Habouzit, Marta Volonteri, Rachel S. Somerville, Yohan Dubois, Sébastien Peirani, Christophe Pichon, Julien Devriendt

Abstract: High-redshift quasars are believed to reside in highly biased regions of the Universe, where black hole (BH) growth is sustained by an enhanced number of mergers and by being at the intersection of filaments bringing fresh gas. This assumption should be supported by an enhancement of the number counts of galaxies in the field of view of quasars. While the current observations of quasar environment… ▽ More High-redshift quasars are believed to reside in highly biased regions of the Universe, where black hole (BH) growth is sustained by an enhanced number of mergers and by being at the intersection of filaments bringing fresh gas. This assumption should be supported by an enhancement of the number counts of galaxies in the field of view of quasars. While the current observations of quasar environments do not lead to a consensus on a possible excess of galaxies, the future missions JWST, WFIRST, and Euclid will provide new insights on quasar environments, and will substantially increase the number of study-cases. We are in a crucial period, where we need to both understand the current observations and predict how upcoming missions will improve our understanding of BH environments. Using the large-scale simulation Horizon-AGN, we find that statistically the most massive BHs reside in environments with the largest number counts of galaxies. However, we find a large variance in galaxy number counts, and some massive BHs do not show enhanced counts in their neighborhood. Interestingly, some massive BHs have a very close galaxy companion but no further enhancement of the galaxy number counts at larger scales, in agreement with recent observations. We find that AGN feedback in the surrounding galaxies is able to decrease their luminosity and stellar mass, and therefore to make them un-observable when using restrictive galaxy selection criteria. Radiation from the quasars can spread over large distances, which could affect the formation history of surrounding galaxies, but a careful analysis of these processes requires radiative transfer simulations. △ Less

Submitted 13 December, 2019; v1 submitted 26 October, 2018; originally announced October 2018.

Comments: published version, 26 pages, 15 figures

Linking galaxy structural properties and star formation activity to black hole activity with IllustrisTNG

Authors: Melanie Habouzit, Shy Genel, Rachel S. Somerville, Dale Kocevski, Michaela Hirschmann, Avishai Dekel, Ena Choi, Dylan Nelson, Annalisa Pillepich, Paul Torrey, Lars Hernquist, Mark Vogelsberger, Rainer Weinberger, Volker Springel

Abstract: We study the connection between active galactic nuclei (AGN) and their host galaxies through cosmic time in the large-scale cosmological IllustrisTNG simulations. We first compare BH properties, i.e. the hard X-ray BH luminosity function, AGN galaxy occupation fraction, and distribution of Eddington ratios, to available observational constraints. The simulations produce a population of BHs in good… ▽ More We study the connection between active galactic nuclei (AGN) and their host galaxies through cosmic time in the large-scale cosmological IllustrisTNG simulations. We first compare BH properties, i.e. the hard X-ray BH luminosity function, AGN galaxy occupation fraction, and distribution of Eddington ratios, to available observational constraints. The simulations produce a population of BHs in good agreement with observations, but we note an excess of faint AGN in hard X-ray (L_x ~ 10^{43-44} erg/s), and a lower number of bright AGN (L_x>10^{44} erg/s), a conclusion that varies quantitatively but not qualitatively with BH luminosity estimation method. The lower Eddington ratios of the 10^{9} Msun BHs compared to observations suggest that AGN feedback may be too efficient in this regime. We study galaxy star formation activity and structural properties, and design sample-dependent criteria to identify different galaxy types (star-forming/quiescent, extended/compact) that we apply both to the simulations and observations from the candels fields. We analyze how the simulated and observed galaxies populate the specific star formation rate - stellar mass surface density diagram. A large fraction of the z=0 M_{star}>10^{11} Msun quiescent galaxies first experienced a compaction phase (i.e. reduction of galaxy size) while still forming stars, and then a quenching event. We measure the dependence of AGN fraction on galaxies' locations in this diagram. After correcting the simulations with a redshift and AGN luminosity-dependent model for AGN obscuration, we find good qualitative and quantitative agreement with observations. The AGN fraction is the highest among compact star-forming galaxies (16-20% at z~1.5-2), and the lowest among compact quiescent galaxies (6-10% at z~1.5-2). △ Less

Submitted 8 January, 2019; v1 submitted 14 September, 2018; originally announced September 2018.

Comments: 35 pages, 22 figures, accepted for publication in MNRAS

On the formation of the first quasars

Authors: Rosa Valiante, Bhaskar Agarwal, Melanie Habouzit, Edwige Pezzulli

Abstract: Observations of the most luminous quasars at redshift z>6 reveal the existence of numerous supermasssive black holes (>10^9 Msun) already in place about twelve billion years ago. In addition, the interstellar medium of the galaxies hosting these black holes are observed to be chemically mature systems, with metallicities (Z>Zsun) and dust masses (>10^8 Msun) similar to that of more evolved, local… ▽ More Observations of the most luminous quasars at redshift z>6 reveal the existence of numerous supermasssive black holes (>10^9 Msun) already in place about twelve billion years ago. In addition, the interstellar medium of the galaxies hosting these black holes are observed to be chemically mature systems, with metallicities (Z>Zsun) and dust masses (>10^8 Msun) similar to that of more evolved, local galaxies. The connection between the rapid growth of the first supermassive black holes and the fast chemical evolution of the host galaxy is one of the most puzzling issues for theoretical models. Here we review state-of-the-art theoretical models that focus on this problem with particular emphasis on the conditions that lead to the formation of quasar seeds and their subsequent evolution at z>6. △ Less

Submitted 10 March, 2017; originally announced March 2017.

Comments: To be published in PASA. Comments are welcome

Density profile of dark matter haloes and galaxies in the Horizon-AGN simulation: the impact of AGN feedback

Authors: S. Peirani, Y. Dubois, M. Volonteri, J. Devriendt, K. Bundy, J. Silk, C. Pichon, S. Kaviraj, R. Gavazzi, M. Habouzit

Abstract: Using a suite of three large cosmological hydrodynamical simulations, Horizon-AGN, Horizon-noAGN (no AGN feedback) and Horizon-DM (no baryons), we investigate how a typical sub-grid model for AGN feedback affects the evolution of the inner density profiles of massive dark matter haloes and galaxies. Based on direct object-to-object comparisons, we find that the integrated inner mass and density sl… ▽ More Using a suite of three large cosmological hydrodynamical simulations, Horizon-AGN, Horizon-noAGN (no AGN feedback) and Horizon-DM (no baryons), we investigate how a typical sub-grid model for AGN feedback affects the evolution of the inner density profiles of massive dark matter haloes and galaxies. Based on direct object-to-object comparisons, we find that the integrated inner mass and density slope differences between objects formed in these three simulations (hereafter, H_AGN, H_noAGN and H_DM) significantly evolve with time. More specifically, at high redshift (z~5), the mean central density profiles of H_AGN and H_noAGN dark matter haloes tend to be much steeper than their H_DM counterparts owing to the rapidly growing baryonic component and ensuing adiabatic contraction. By z~1.5, these mean halo density profiles in H_AGN have flattened, pummelled by powerful AGN activity ("quasar mode"): the integrated inner mass difference gaps with H_noAGN haloes have widened, and those with H_DM haloes have narrowed. Fast forward 9.5 billion years, down to z=0, and the trend reverses: H_AGN halo mean density profiles drift back to a more cusped shape as AGN feedback efficiency dwindles ("radio mode"), and the gaps in integrated central mass difference with H_noAGN and H_DM close and broaden respectively. On the galaxy side, the story differs noticeably. Averaged stellar profile central densities and inner slopes are monotonically reduced by AGN activity as a function of cosmic time, resulting in better agreement with local observations. As both dark matter and stellar inner density profiles respond quite sensitively to the presence of a central AGN, there is hope that future observational determinations of these quantities can be used constrain AGN feedback models. △ Less

Submitted 16 August, 2017; v1 submitted 29 November, 2016; originally announced November 2016.

Comments: 18 pages, 15 figures, accepted for publication in MNRAS. New plots added

Blossoms from black hole seeds: properties and early growth regulated by supernova feedback

Authors: Melanie Habouzit, Marta Volonteri, Yohan Dubois

Abstract: Massive black holes (BHs) inhabit local galaxies, including the Milky Way and some dwarf galaxies. BH formation, occurring at early cosmic times, must account for the properties of BHs in today's galaxies, notably why some galaxies host a BH, and others do not. We investigate the formation, distribution and growth of BH `seeds' by using the adaptive mesh refinement code Ramses. We develop an imple… ▽ More Massive black holes (BHs) inhabit local galaxies, including the Milky Way and some dwarf galaxies. BH formation, occurring at early cosmic times, must account for the properties of BHs in today's galaxies, notably why some galaxies host a BH, and others do not. We investigate the formation, distribution and growth of BH `seeds' by using the adaptive mesh refinement code Ramses. We develop an implementation of BH formation in dense, low-metallicity environments, as advocated by models invoking the collapse of the first generation of stars, or of dense nuclear star clusters. The seed masses are computed one-by-one on-the-fly, based on the star formation rate and the stellar initial mass function. This self-consistent method to seed BHs allows us to study the distribution of BHs in a cosmological context and their evolution over cosmic time. We find that all high-mass galaxies tend to a host a BH, whereas low-mass counterparts have a lower probability of hosting a BH. After the end of the epoch of BH formation, this probability is modulated by the growth of the galaxy. The simulated BHs connect to low-redshift observational samples, and span a similar range in accretion properties as Lyman-Break Analogs. The growth of BHs in low-mass galaxies is stunted by strong supernova feedback. The properties of BHs in dwarf galaxies thus remain a testbed for BH formation. Simulations with strong supernova feedback, which is able to quench BH accretion in shallow potential wells, produce galaxies and BHs in better agreement with observational constraints. △ Less

Submitted 26 July, 2017; v1 submitted 30 May, 2016; originally announced May 2016.

Comments: 16 pages, 16 figures, matched the published version

On the number density of "direct collapse" black hole seeds

Authors: Melanie Habouzit, Marta Volonteri, Muhammad Latif, Yohan Dubois, Sebastien Peirani

Abstract: Supermassive black holes (BHs) of millions solar masses and above reside in the center of most local galaxies, but they also power active galactic nuclei and quasars, detected up to z=7. This observational evidence puts strong constraints on the BH growth and the mass of the first BH seeds. The scenario of "direct collapse" is very appealing as it leads to the formation of large mass BH seeds in t… ▽ More Supermassive black holes (BHs) of millions solar masses and above reside in the center of most local galaxies, but they also power active galactic nuclei and quasars, detected up to z=7. This observational evidence puts strong constraints on the BH growth and the mass of the first BH seeds. The scenario of "direct collapse" is very appealing as it leads to the formation of large mass BH seeds in the range 10^4-10^6 Msun, which eases explaining how quasars at z=6-7 are powered by BHs with masses >10^9 Msun. Direct collapse, however, appears to be rare, as the conditions required by the scenario are that gas is metal-free, the presence of a strong photo-dissociating Lyman-Werner flux, and large inflows of gas at the center of the halo, sustained for 10-100 Myr. We performed several cosmological hydrodynamical simulations that cover a large range of box sizes and resolutions, thus allowing us to understand the impact of several physical processes on the distribution of direct collapse BHs. We identify halos where direct collapse can happen, and derive the number density of BHs. We also investigate the discrepancies between hydrodynamical simulations, direct or post-processed, and semi-analytical studies. We find that for direct collapse to account for BHs in normal galaxies, the critical Lyman-Werner flux required for direct collapse must be much less than predicted by 3D simulations that include detailed chemical models. However, when supernova feedback is relatively weak, enough direct collapse BHs to explain z=6-7 quasars can be obtained for more realistic values of the critical Lyman Werner flux. △ Less

Submitted 13 December, 2019; v1 submitted 4 January, 2016; originally announced January 2016.

Comments: matched the version published in 2016, 13 pages, 7 figures

The evolution of high-redshift massive black holes

Authors: Marta Volonteri, Melanie Habouzit, Fabio Pacucci, Michael Tremmel

Abstract: Massive black holes (MBHs) are nowadays recognized as integral parts of galaxy evolution. Both the approximate proportionality between MBH and galaxy mass, and the expected importance of feedback from active MBHs in regulating star formation in their host galaxies point to a strong interplay between MBHs and galaxies. MBHs must form in the first galaxies and be fed by gas in these galaxies, with c… ▽ More Massive black holes (MBHs) are nowadays recognized as integral parts of galaxy evolution. Both the approximate proportionality between MBH and galaxy mass, and the expected importance of feedback from active MBHs in regulating star formation in their host galaxies point to a strong interplay between MBHs and galaxies. MBHs must form in the first galaxies and be fed by gas in these galaxies, with continuous or intermittent inflows that, at times, can be larger than the Eddington rate. Feedback from supernovae and from the MBHs themselves modulates the growth of the first MBHs. While current observational data only probe the most massive and luminous MBHs, the tip of the iceberg, we will soon be able to test theoretical models of MBH evolution on more "normal" MBHs: the MBHs that are indeed relevant in building the population that we observe in local galaxies, including our own Milky Way. △ Less

Submitted 9 November, 2015; originally announced November 2015.

Comments: 8 pages. Brief review to appear in Galaxies at High Redshift and Their Evolution over Cosmic Time, Proceedings IAU Symposium No. 319, 2015

Impact of dust cooling on direct collapse black hole formation

Authors: M. A. Latif, K. Omukai, M. Habouzit, D. R. G. Schleicher, M. Volonteri

Abstract: Observations of quasars at $ z > 6$ suggest the presence of black holes with a few times $\rm 10^9 ~M_{\odot}$. Numerous models have been proposed to explain their existence including the direct collapse which provides massive seeds of $\rm 10^5~M_{\odot}$. The isothermal direct collapse requires a strong Lyman-Werner flux to quench $\rm H_2$ formation in massive primordial halos. In this study, w… ▽ More Observations of quasars at $ z > 6$ suggest the presence of black holes with a few times $\rm 10^9 ~M_{\odot}$. Numerous models have been proposed to explain their existence including the direct collapse which provides massive seeds of $\rm 10^5~M_{\odot}$. The isothermal direct collapse requires a strong Lyman-Werner flux to quench $\rm H_2$ formation in massive primordial halos. In this study, we explore the impact of trace amounts of metals and dust enrichment. We perform three dimensional cosmological simulations for two halos of $\rm > 10^7~M_{\odot}$ with $\rm Z/Z_{\odot}= 10^{-4}-10^{-6}$ illuminated by an intense Lyman Werner flux of $\rm J_{21}=10^5$. Our results show that initially the collapse proceeds isothermally with $\rm T \sim 8000$ K but dust cooling becomes effective at densities of $\rm 10^{8}-10^{12} ~cm^{-3}$ and brings the gas temperature down to a few 100-1000 K for $\rm Z/Z_{\odot} \geq 10^{-6}$. No gravitationally bound clumps are found in $\rm Z/Z_{\odot} \leq 10^{-5}$ cases by the end of our simulations in contrast to the case with $\rm Z/Z_{\odot} = 10^{-4}$. Large inflow rates of $\rm \geq 0.1~M_{\odot}/yr$ are observed for $\rm Z/Z_{\odot} \leq 10^{-5}$ similar to a zero-metallicity case while for $\rm Z/Z_{\odot} = 10^{-4}$ the inflow rate starts to decline earlier due to the dust cooling and fragmentation. For given large inflow rates a central star of $\rm \sim 10^4~M_{\odot}$ may form for $\rm Z/Z_{\odot} \leq 10^{-5}$. △ Less

Submitted 24 March, 2016; v1 submitted 23 September, 2015; originally announced September 2015.

Comments: Accepted for publication in ApJ, comments are still welcome

Black hole formation and growth with non-Gaussian primordial density perturbations

Authors: Mélanie Habouzit, Marta Volonteri, Muhammad Latif, Takahiro Nishimichi, Sébastien Peirani, Yohan Dubois, Gary A. Mamon, Joseph Silk, Jacopo Chevallard

Abstract: Quasars powered by massive black holes (BHs) with mass estimates above a billion solar masses have been identified at redshift 6 and beyond. The existence of such BHs requires almost continuous growth at the Eddington limit for their whole lifetime, of order of one billion years. In this paper, we explore the possibility that positively skewed scale-dependent non-Gaussian primordial fluctuations m… ▽ More Quasars powered by massive black holes (BHs) with mass estimates above a billion solar masses have been identified at redshift 6 and beyond. The existence of such BHs requires almost continuous growth at the Eddington limit for their whole lifetime, of order of one billion years. In this paper, we explore the possibility that positively skewed scale-dependent non-Gaussian primordial fluctuations may ease the assembly of massive BHs. In particular, they produce more low-mass halos at high redshift, thus altering the production of metals and ultra-violet flux, believed to be important factors in BH formation. Additionally, a higher number of progenitors and of nearly equal-mass halo mergers would boost the mass increase provided by BH-BH mergers and merger-driven accretion. We use a set of two cosmological simulations, with either Gaussian or scale-dependent non-Gaussian primordial fluctuations to perform a proof-of-concept experiment to estimate how BH formation and growth are altered. We estimate the BH number density and the fraction of halos where BHs form, for both simulations and for two popular scenarios of BH formation (remnants of the first generation of stars and direct collapse in the absence of metals and molecular hydrogen). We find that the fractions of halos where BHs form are almost identical, but that non-Gaussian primordial perturbations increase the total number density of BHs for the both BH formation scenarios. We also evolve BHs using merger trees extracted from the simulations and find that non-Gaussianities increase both the BH masses and the number of the most massive BHs. △ Less

Submitted 13 December, 2015; v1 submitted 21 July, 2015; originally announced July 2015.

Comments: 11 pages, 10 figures, MNRAS accepted

Report number: MN-15-2256-MJ.R2

Effect of primordial non-Gaussianities on the far-UV luminosity function of high-redshift galaxies: implications for cosmic reionization

Authors: Jacopo Chevallard, Joseph Silk, Takahiro Nishimichi, Melanie Habouzit, Gary A. Mamon, Sébastien Peirani

Abstract: [Abridged] Understanding how the intergalactic medium (IGM) was reionized at z > 6 is one of the big challenges of current high redshift astronomy. It requires modelling the collapse of the first astrophysical objects (Pop III stars, first galaxies) and their interaction with the IGM, while at the same time pushing current observational facilities to their limits. The observational and theoretical… ▽ More [Abridged] Understanding how the intergalactic medium (IGM) was reionized at z > 6 is one of the big challenges of current high redshift astronomy. It requires modelling the collapse of the first astrophysical objects (Pop III stars, first galaxies) and their interaction with the IGM, while at the same time pushing current observational facilities to their limits. The observational and theoretical progress of the last few years have led to the emergence of a coherent picture in which the budget of hydrogen-ionizing photons is dominated by low-mass star-forming galaxies, with little contribution from Pop III stars and quasars. The reionization history of the Universe therefore critically depends on the number density of low-mass galaxies at high redshift. In this work, we explore how changes in the statistical properties of initial density fluctuations affect the formation of early galaxies. Following Habouzit et al. (2014), we run 5 N-body simulations with Gaussian and (scale-dependent) non-Gaussian initial conditions, all consistent with Planck constraints. By appealing to a galaxy formation model and to a population synthesis code, we compute the far-UV galaxy luminosity function down to M_UV = -14 at redshift 7 < z < 15. We find that models with strong primordial non-Gaussianities on < Mpc scales show a far-UV luminosity function significantly enhanced in low-mass galaxies. We adopt a reionization model calibrated from state-of-the-art hydrodynamical simulations and show that such non-Gaussianities leave a clear imprint on the Universe reionization history and electron Thomson scattering optical depth tau_E. Although current uncertainties in the physics of reionization and on the determination of tau_E still dominate the signatures of non-Gaussianities, our results suggest that tau_E could ultimately be used to constrain the statistical properties of initial density fluctuations. △ Less

Submitted 28 October, 2014; originally announced October 2014.

Comments: 18 pages, 12 figures, accepted for publication in MNRAS

Testing primordial non-Gaussianities on galactic scales at high redshift

Authors: Mélanie Habouzit, Takahiro Nishimichi, Sébastien Peirani, Gary A. Mamon, Joseph Silk, Jacopo Chevallard

Abstract: Primordial non-Gaussianities provide an important test of inflationary models. Although the Planck CMB experiment has produced strong limits on non-Gaussianity on scales of clusters, there is still room for considerable non-Gaussianity on galactic scales. We have tested the effect of local non-Gaussianity on the high redshift galaxy population by running five cosmological N-body simulations down t… ▽ More Primordial non-Gaussianities provide an important test of inflationary models. Although the Planck CMB experiment has produced strong limits on non-Gaussianity on scales of clusters, there is still room for considerable non-Gaussianity on galactic scales. We have tested the effect of local non-Gaussianity on the high redshift galaxy population by running five cosmological N-body simulations down to z=6.5. For these simulations, we adopt the same initial phases, and either Gaussian or scale-dependent non-Gaussian primordial fluctuations, all consistent with the constraints set by Planck on clusters scales. We then assign stellar masses to each halo using the halo - stellar mass empirical relation of Behroozi et al. (2013). Our simulations with non-Gaussian initial conditions produce halo mass functions that show clear departures from those obtained from the analogous simulations with Gaussian initial conditions at z>~10. We observe a >0.3 dex enhancement of the low-end of the halo mass function, which leads to a similar effect on the galaxy stellar mass function, which should be testable with future galaxy surveys at z>10. As cosmic reionization is thought to be driven by dwarf galaxies at high redshift, our findings may have implications for the reionization history of the Universe. △ Less

Submitted 9 September, 2014; v1 submitted 30 July, 2014; originally announced July 2014.

Comments: 6 pages, 3 figures, 1 table, MNRAS (Letters) in press