Impact of dust cooling on direct collapse black hole formation

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}$.Comment: Accepted for publication in ApJ, comments are still welcom

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

This article has been accepted for publication in Monthly Notices of the Royal Astronomical Society. ©: 2017 The Author(s). Published by Oxford University Press on behalf of the Royal Astronomical Society.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.Peer reviewedFinal Published versio

Active galactic nucleus X-ray luminosity function and absorption function in the Early Universe (3 ≤ z ≤ 6)

The X-ray luminosity function (XLF) of active galactic nuclei (AGN) offers a robust tool to study the evolution and the growth of the supermassive black-hole population 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 ≥ 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 (z ≥ 3) using a large sample of AGN 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 used more than 600 soft X-ray selected (0.5 − 2 keV) high-z sources in the Chandra deep fields, the Chandra COSMOS Legacy survey, and the XMM-XXL northern field. We derived the X-ray spectral properties for all sources via spectral fitting, using a consistent technique and model. To model the parametric form of the XLF and the absorption function, we used 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 NH ≥ 1023 cm−2, while ∼17% of the sources are Compton-Thick. Our results favour a scenario where both the interstellar medium of the host and the AGN torus contribute to the obscuration. The derived black hole accretion rate density is roughly in agreement with the large-scale cosmological hydrodynamical simulations, if one takes into account the results that the X-ray AGN are hosted by massive galaxies, while it differs from that derived using JWST data. The latter could be due to the differences in the AGN and host-galaxy properties

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

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).Comment: Roman Core Community Survey White Paper, focusing on the High Latitude Wide Area Surve

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

In this paper we present the case for void science, arguing that cosmic voids are a novel probe to constrain modified gravity, dark energy, the sum of neutrino masses and galaxy evolution. Voids will answer some of the most relevant questions in cosmology and astrophysics over the next decade. <p/

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

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

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

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.Comment: White Paper for the Nancy Grace Roman Space Telescope's Core Community Surveys (https://roman.gsfc.nasa.gov/science/ccs_white_papers.html

Astrophysics with the Laser Interferometer Space Antenna

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; ultra-compact stellar-mass binaries, massive black hole binaries, and extreme or interme-diate 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

Astrophysics with the Laser Interferometer Space Antenna

Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy as it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and other space-based 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: ultra-compact 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 make progress in the different areas. New research avenues that LISA itself, or its joint exploitation with 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