Showing 1–50 of 173 results for author: Haardt, F

Cosmic Reionization in the JWST Era: Back to AGNs?

Authors: Piero Madau, Emanuele Giallongo, Andrea Grazian, Francesco Haardt

Abstract: Deep surveys with the James Webb Space Telescope (JWST) have revealed an emergent population of moderate-luminosity, broad-line active galactic nuclei (AGNs) at 4< z< 14 powered by accretion onto early massive black holes. The high number densities reported, together with the large Lyman-continuum (LyC) production efficiency and leakiness into the intergalactic medium (IGM) that are typical of UV-… ▽ More Deep surveys with the James Webb Space Telescope (JWST) have revealed an emergent population of moderate-luminosity, broad-line active galactic nuclei (AGNs) at 4< z< 14 powered by accretion onto early massive black holes. The high number densities reported, together with the large Lyman-continuum (LyC) production efficiency and leakiness into the intergalactic medium (IGM) that are typical of UV-selected AGNs, lead us to reassess a scenario where AGNs are the sole drivers of the cosmic hydrogen/helium reionization process. Our approach is based on the assumptions, grounded in recent observations, that: (a) the fraction of broad-line AGNs among galaxies is around 10-15%; (b) the mean escape fraction of hydrogen LyC radiation is high, >80%, in AGN hosts and is negligible otherwise; and (c) internal absorption at 4 ryd or a steep ionizing EUV spectrum delay full reionization of HeII until z~2.8-3.0, in agreement with observations of the HeII Lyman-alpha forest. In our fiducial models: 1) hydrogen reionization is 99% completed by redshift z~5.3-5.5, and reaches its midpoint at z~6.5-6.7; (2) the integrated Thomson scattering optical depth to reionization is ~0.05, consistent with constraints from cosmic microwave background (CMB) anisotropy data; and (3) the abundant AGN population detected by JWST does not violate constraints on the unresolved X-ray background. △ Less

Submitted 26 June, 2024; originally announced June 2024.

Comments: 13 pages, 5 figures, accepted for publication in The Astrophysical Journal

On the universal validity of Case B recombination theory

Authors: C. Scarlata, M. Hayes, N. Panagia, V. Mehta, F. Haardt, M. Bagley

Abstract: In an ongoing search for low-mass extreme emission line galaxies, we identified a galaxy with a Ha/Hb Balmer line ratio of 2.620 +- 0.078. Ha/Hb Balmer ratios lower than the dust-free Case~B value appear relatively frequently in extreme emission line galaxies. These low values suggest that the Case~B assumption may not be valid in these objects. After ruling out the possibility that the low Ha/Hb… ▽ More In an ongoing search for low-mass extreme emission line galaxies, we identified a galaxy with a Ha/Hb Balmer line ratio of 2.620 +- 0.078. Ha/Hb Balmer ratios lower than the dust-free Case~B value appear relatively frequently in extreme emission line galaxies. These low values suggest that the Case~B assumption may not be valid in these objects. After ruling out the possibility that the low Ha/Hb ratio is due to systematic errors introduced by observational effects, we use constraints from the total Hb luminosity, the [OIII]/[OII] line ratio and the Balmer line equivalent widths, to suggest that the gas is optically thick to both Ha and Lya photons, and the geometry and orientation of the scattering gas causes Ha photons to be preferentially removed from the line of sight with respect to higher order Balmer series photons. Finally, we use data from the SDSS survey to show that Balmer self-absorption may be more important than previously assumed in high excitation emission line galaxies, where Lya pumping of the hydrogen excited state can be effective. If not recognized, Balmer self-absorption could lead to inaccurate estimates of galaxy physical properties. As an example, the effect of dust extinction could be over-estimated, for spherically symmetric scattering medium, or under-estimated, for a not spherically-symmetric distribution. △ Less

Submitted 29 April, 2024; v1 submitted 13 April, 2024; originally announced April 2024.

BANG-MaNGA: A census of kinematic discs and bulges across mass and star formation in the local Universe

Authors: Fabio Rigamonti, Luca Cortese, Francesco Bollati, Stefano Covino, Massimo Dotti, A. Fraser-McKelvie, Francesco Haardt

Abstract: We investigate the relevance of kinematically identified bulges, discs and their role relative to galaxy quenching. We utilize an analysis of the SDSS-MaNGA survey conducted with the GPU-based code BANG which simultaneously models galaxy photometry and kinematics to decompose galaxies into their structural components. Below M~1011 Msun, galaxies exhibit a wide range of dynamical properties, determ… ▽ More We investigate the relevance of kinematically identified bulges, discs and their role relative to galaxy quenching. We utilize an analysis of the SDSS-MaNGA survey conducted with the GPU-based code BANG which simultaneously models galaxy photometry and kinematics to decompose galaxies into their structural components. Below M~1011 Msun, galaxies exhibit a wide range of dynamical properties, determined by the relative prominence of a dispersion-supported inner region and a rotationally-supported disc. Our analysis reveals a natural separation between these classes, with only a minor fraction of stellar mass retained by structures exhibiting intermediate dynamical support. When examining galaxies in terms of their star formation activity, an apparent decrease in rotational support is observed as they move below the star-forming main sequence. This behaviour is evident with luminosity-weighted tracers of kinematics, while it almost vanishes with mass-weighted tracers. Luminosity-weighted quantities not only capture differences in kinematics but also in the stellar population, potentially leading to biased interpretations of galaxy dynamical properties and quenching. Our findings suggest that quenching does not imply almost any structural transformation in galaxies below M~10^11 Msun. Processes as disc fading more likely account for observed differences in mass-weighted and luminosity-weighted galaxy properties; when the galactic disc ceases star formation, its mass-to-light ratio grows without any significant morphological transformation. The picture is remarkably different above M~10^11 Msun. Regardless of the tracer used, a substantial increase in galaxy dispersion support is observed along with a significant structural change. A different quenching mechanism, most likely associated with mergers, dominates. Notably, this mechanism is confined to a very limited range of high masses. △ Less

Submitted 14 March, 2024; originally announced March 2024.

Comments: Accepted for publication in A&A. 16 pages, 10 figures

Evaporation of Close-in Sub-Neptunes by Cooling White Dwarfs

Authors: Elena Gallo, Andrea Caldiroli, Riccardo Spinelli, Federico Biassoni, Francesco Haardt, Mary Anne Limbach, Juliette Becker, Fred Adams

Abstract: Motivated by the recent surge in interest concerning white dwarf (WD) planets, this work presents the first numerical exploration of WD-driven atmospheric escape, whereby the high-energy radiation from a hot/young WD can trigger the outflow of the hydrogen-helium envelope for close-in planets. As a pilot investigation, we focus on two specific cases: a gas giant and a sub-Neptune-sized planet, bot… ▽ More Motivated by the recent surge in interest concerning white dwarf (WD) planets, this work presents the first numerical exploration of WD-driven atmospheric escape, whereby the high-energy radiation from a hot/young WD can trigger the outflow of the hydrogen-helium envelope for close-in planets. As a pilot investigation, we focus on two specific cases: a gas giant and a sub-Neptune-sized planet, both orbiting a rapidly cooling WD with mass $M_\ast$ = 0.6 \msun\ and separation $a$ = 0.02 AU. In both cases, the ensuing mass outflow rates exceed $10^{14}$ g sec$^{-1}$ for WD temperatures greater than $T_{\rm WD} \simeq$ 50,000 K. At $T_{\rm WD} \simeq$ 18,000 K [/22,000 K], the sub-Neptune [/gas giant] mass outflow rate approaches $10^{12}$ g sec$^{-1}$, i.e., comparable to the strongest outflows expected from close-in planets around late main-sequence stars. Whereas the gas giant remains virtually unaffected from an evolutionary standpoint, atmospheric escape may have sizable effects for the sub-Neptune, depending on its dynamical history, e.g., assuming that the hydrogen-helium envelope makes up 1 [/4] per cent of the planet mass, the entire envelope would be evaporated away so long as the planet reaches 0.02 AU within the first 230 [/130] Myr of the WD formation. We discuss how these results can be generalized to eccentric orbits with effective semi-major axis $a'=a/(1-e^2)^{1/4}$, which receive the same orbit-averaged irradiation. Extended to a much broader parameter space, this approach can be exploited to model the expected demographics of WD planets as a function of their initial mass, composition and migration history, as well as their potential for habitability. △ Less

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

Comments: ApJ Letters, in press

A weakly-parametric approach to stochastic background inference in LISA

Authors: Federico Pozzoli, Riccardo Buscicchio, Christopher J. Moore, Francesco Haardt, Alberto Sesana

Abstract: Detecting stochastic gravitational wave backgrounds (SGWBs) with The Laser Interferometer Space Antenna (LISA) is among the mission science objectives. Disentangling SGWBs of astrophysical and cosmological origin is a challenging task, further complicated by the noise level uncertainties. In this study, we introduce a Bayesian methodology to infer upon SGWBs, taking inspiration from Gaussian stoch… ▽ More Detecting stochastic gravitational wave backgrounds (SGWBs) with The Laser Interferometer Space Antenna (LISA) is among the mission science objectives. Disentangling SGWBs of astrophysical and cosmological origin is a challenging task, further complicated by the noise level uncertainties. In this study, we introduce a Bayesian methodology to infer upon SGWBs, taking inspiration from Gaussian stochastic processes. We investigate the suitability of the approach for signal of unknown spectral shape. We do by discretely exploring the model hyperparameters, a first step towards a more efficient transdimensional exploration. We apply the proposed method to a representative astrophysical scenario: the inference on the astrophysical foreground of Extreme Mass Ratio Inspirals, recently estimated in~\cite{Pozzoli2023}. We find the algorithm to be capable of recovering the injected signal even with large priors, while simultaneously providing estimate of the noise level. △ Less

Submitted 20 November, 2023; originally announced November 2023.

Comments: 16 pages, 10 figures, 0 tables

On the connection between AGN radiative feedback and massive black hole spin

Authors: Francesco Bollati, Alessandro Lupi, Massimo Dotti, Francesco Haardt

Abstract: We present a novel implementation for active galactic nucleus (AGN) feedback through ultra-fast winds in the code gizmo. Our feedback recipe accounts for the angular dependence of radiative feedback upon black hole spin. We self-consistently evolve in time i) the gas accretion process from resolved scales to an unresolved AGN disc, ii) the evolution of the spin of the massive black hole (MBH), iii… ▽ More We present a novel implementation for active galactic nucleus (AGN) feedback through ultra-fast winds in the code gizmo. Our feedback recipe accounts for the angular dependence of radiative feedback upon black hole spin. We self-consistently evolve in time i) the gas accretion process from resolved scales to an unresolved AGN disc, ii) the evolution of the spin of the massive black hole (MBH), iii) the injection of AGN-driven winds into the resolved scales, and iv) the spin-induced anisotropy of the overall feedback process. We test our implementation by following the propagation of the wind-driven outflow into a homogeneous medium, and we compare the results against simple analytical models. Then, we consider an isolated galaxy setup and there we study the impact of the AGN feedback on the evolution of the MBH and the of the host galaxy. We find that: i) AGN feedback limits the gas inflow that powers the MBH, with a consequent weak impact on the host galaxy characterized by a star formation (SF) suppression of about a factor of two in the nuclear region; ii) the impact of AGN feedback on the host galaxy and on MBH growth is primarily determined by the AGN luminosity, rather than by its angular pattern set by the MBH spin; iii) the imprint of the angular pattern of the AGN radiation emission manifest in a more clear way at high accretion rates. At such high rates the more isotropic angular patterns, proper to higher spin values, sweep away gas in the nuclear region more easily, hence causing a slower MBH mass and spin growths and a higher quenching of the SF. We argue that the influence of spin-dependent anisotropy of AGN feedback on MBH and galaxy evolution is likely to be relevant in those scenarios characterized by high and prolonged MBH accretion episodes and by high AGN wind-galaxy coupling. Such conditions are more frequently met in galaxy mergers and/or high redshift galaxies. △ Less

Submitted 13 November, 2023; originally announced November 2023.

Self-Consistent Modeling of Metastable Helium Exoplanet Transits

Authors: Federico Biassoni, Andrea Caldiroli, Elena Gallo, Francesco Haardt, Riccardo Spinelli, Francesco Borsa

Abstract: Absorption of stellar X-ray and Extreme Ultraviolet radiation in the upper atmosphere of close-in exoplanets can give rise to hydrodynamic outflows, which may lead to the gradual shedding of their primordial, light element envelopes. Excess absorption by neutral helium atoms in the metastable state has recently emerged as a viable diagnostic of atmospheric escape. Here we present a public module t… ▽ More Absorption of stellar X-ray and Extreme Ultraviolet radiation in the upper atmosphere of close-in exoplanets can give rise to hydrodynamic outflows, which may lead to the gradual shedding of their primordial, light element envelopes. Excess absorption by neutral helium atoms in the metastable state has recently emerged as a viable diagnostic of atmospheric escape. Here we present a public module to the 1D photo-ionization hydrodynamic code ATES, designed to calculate the HeI triplet transmission probability for a broad range of planetary parameters. By relaxing the isothermal outflow assumption, the code enables a self-consistent assessment of the HeI triplet absorption depth along with the atmospheric mass loss rate and the outflow temperature profile, which strongly affects the recombination rate of HeII into HeI triplet. We investigate how the transit signal can be expected to depend upon known system parameters, including host spectral type, orbital distance, as well as planet gravity. At variance with previous studies, which identified K-type stars as favorable hosts, we conclude that late M-dwarfs with Neptune-sized planets orbiting at ~0.05-0.1 AU can be expected to yield the strongest transit signal well in excess of 30% for near-cosmological He/H abundances. More generally, we show that the physics which regulates the population and depletion of the metastable state, combined with geometrical effects, can yield somewhat counter-intuitive results, such as a non-monotonic dependence of the transit depth on orbital distance. These are compounded by a strong degeneracy between the stellar EUV flux intensity and the atmospheric He/H abundance, both of which are highly uncertain. Compared against spectroscopy data our modelling suggests that either a large fraction of the targets have helium depleted envelopes, or, that the input stellar EUV spectra are systematically overestimated. △ Less

Submitted 19 October, 2023; originally announced October 2023.

Comments: A&A accepted. The ATES code and the Transmission Probability Module (TPM) can be found at https://github.com/AndreaCaldiroli/ATES-Code

Iron K-alpha echoes from the circumnuclear disk orbiting Sgr A*

Authors: G. Stel, G. Ponti, F. Haardt

Abstract: Molecular clouds in the Galactic center (GC) reprocess radiation from past outbursts of nearby high-energy sources, generating a bright Fe K-alpha fluorescence at 6.4 keV. The closest clouds to the GC are only $\simeq 1.5$ pc from Sgr A*, forming a torus-like structure known as the circumnuclear disk (CND). The study of fluorescence emission can lead to a characterization of the illuminating sourc… ▽ More Molecular clouds in the Galactic center (GC) reprocess radiation from past outbursts of nearby high-energy sources, generating a bright Fe K-alpha fluorescence at 6.4 keV. The closest clouds to the GC are only $\simeq 1.5$ pc from Sgr A*, forming a torus-like structure known as the circumnuclear disk (CND). The study of fluorescence emission can lead to a characterization of the illuminating source(s), the reflecting clouds, and the global geometry of such a system lying in the GC. The primary purpose of our study is to analyze possible fluorescence signals arising in the CND. This signal would allow us to constrain the CND's physical properties and the source-reflector system's geometry. By exploiting the last $\simeq 20$ years of XMM-Newton observations of the GC, we studied the variability of the Fe K-alpha line in the region around Sgr A*. We identified regions with a flux excess and computed the spectrum therein. We then derived the hydrogen column density of the CND after relating the intensity of the 6.4 keV line to the total energy emitted by known transient sources in the region. Starting from data collected in 2019, we find significant line excesses in a region compatible with the eastern portion of the CND. The echo radiation can be linked to the 2013 outburst of the magnetar SGR J1745-2900. We derive a mean effective hydrogen column density of the CND in the eastern region of $\simeq 10^{23}$ cm$^{-2}$. The scenario depicted is physically plausible, given the luminosity, the position of the illuminating source, and the expected density of the CND. △ Less

Submitted 15 September, 2023; originally announced September 2023.

Comments: 15 pages, 18 figures, to be published in Astronomy & Astrophysics

Journal ref: A&A 679, A44 (2023)

GRMHD simulations of accretion flows onto unequal-mass, precessing massive binary black hole mergers

Authors: Federico Cattorini, Bruno Giacomazzo, Monica Colpi, Francesco Haardt

Abstract: In this work, we use general relativistic magnetohydrodynamics simulations to explore the effect of spin orientation on the dynamics of gas in the vicinity of merging black holes. We present a suite of eight simulations of unequal-mass, spinning black hole binaries embedded in magnetized clouds of matter. Each binary evolution covers approximately 15 orbits before the coalescence. The geometry of… ▽ More In this work, we use general relativistic magnetohydrodynamics simulations to explore the effect of spin orientation on the dynamics of gas in the vicinity of merging black holes. We present a suite of eight simulations of unequal-mass, spinning black hole binaries embedded in magnetized clouds of matter. Each binary evolution covers approximately 15 orbits before the coalescence. The geometry of the accretion flows in the vicinity of the black holes is significantly altered by the orientation of the individual spins with respect to the orbital angular momentum, with the primary black hole dominating the mass accretion rate $\dot{M}$. We observe quasiperiodic modulations of $\dot{M}$ in most of the configurations, whose amplitude is dependent on the orientation of the black hole spins. We find the presence of a relation between the average amplitude of $\dot{M}$ and the spin precession parameter $χ_{\mathrm{p}}$ showing that spin misalignment systematically leads to stronger modulation, whereas configurations with spins aligned to the orbital angular momentum damp out the quasiperiodicity. This finding suggests a possible signature imprinted in the accretion luminosity of precessing binaries approaching merger and has possible consequences on future multimessenger observations of massive binary black hole systems. △ Less

Submitted 11 September, 2023; originally announced September 2023.

Comments: 12 pages, 8 figures. Submitted to PRD

The MUSE Ultra Deep Field (MUDF). IV. A pair of X-ray weak quasars at the heart of two extended Lyα nebulae

Authors: Elisabeta Lusso, Emanuele Nardini, Michele Fumagalli, Matteo Fossati, Fabrizio Arrigoni Battaia, Mitchell Revalski, Marc Rafelski, Valentina D'Odorico, Celine Peroux, Stefano Cristiani, Pratika Dayal, Francesco Haardt, Emma K. Lofthouse

Abstract: We present the results obtained from follow-up observations of the MUSE Ultra Deep Field (MUDF) at X-ray energies with XMM-Newton. The MUDF is centred on a unique field with two bright, physically associated quasars at $z\simeq3.23$, separated by $\sim$500 kpc in projection. Both quasars are embedded within extended Ly$α$ nebulae ($\gtrsim 100~\rm kpc$ at a surface brightness flux level of… ▽ More We present the results obtained from follow-up observations of the MUSE Ultra Deep Field (MUDF) at X-ray energies with XMM-Newton. The MUDF is centred on a unique field with two bright, physically associated quasars at $z\simeq3.23$, separated by $\sim$500 kpc in projection. Both quasars are embedded within extended Ly$α$ nebulae ($\gtrsim 100~\rm kpc$ at a surface brightness flux level of $\approx 6\times 10^{-19} \rm erg~s^{-1}~cm^{-2}~arcsec^{-2}$), whose elongated morphology is suggestive of an extended filament connecting the quasar haloes. The new X-ray observations presented here allow us to characterise the physical properties (e.g. X-ray slope, luminosities, gas column densities) in the innermost region of the MUDF quasars. We find that both quasars are X-ray underluminous compared to objects at similar ultraviolet luminosities. Based on our X-ray spectral analysis, absorbing columns of $N_H(z)\gtrsim$ 10$^{23}$ cm$^{-2}$ appear unlikely, therefore such a weakness is possibly intrinsic. When also including literature data, we do not observe any detectable trend between the area of the nebulae and nuclear luminosities at both the rest-frame 2 keV and 2500 $\rm \mathring{A}$. The area is also not correlated with the X-ray photon index nor with the integrated band flux in the hard band (2$-$10 keV). We also do not find any trend between the extended Ly$α$ emission of the nebulae and the nuclear X-ray luminosity. Finally, the properties of the MUDF quasars' nebulae are consistent with the observed relation between the Ly$α$ integrated luminosity of the nebulae and their area. Our results suggest that the quasar ionization power is not a strong driver of the morphology and size of the nebulae. △ Less

Submitted 4 September, 2023; v1 submitted 24 August, 2023; originally announced August 2023.

Comments: 15 pages, 9 figures, reference added, published in MNRAS

The GAPS programme at TNG. XLV. HI Balmer lines transmission spectroscopy and NLTE atmospheric modelling of the ultra-hot Jupiter KELT-20b/MASCARA-2b

Authors: L. Fossati, F. Biassoni, G. M. Cappello, F. Borsa, D. Shulyak, A. S. Bonomo, D. Gandolfi, F. Haardt, T. Koskinen, A. F. Lanza, V. Nascimbeni, D. Sicilia, M. Young, G. Aresu, A. Bignamini, M. Brogi, I. Carleo, R. Claudi, R. Cosentino, G. Guilluy, C. Knapic, L. Malavolta, L. Mancini, D. Nardiello, M. Pinamonti , et al. (5 additional authors not shown)

Abstract: We aim at extracting the transmission spectrum of the HI Balmer lines of the ultra-hot Jupiter (UHJ) KELT-20b/MASCARA-2b from observations and to further compare the results with what obtained through forward modelling accounting for non-local thermodynamic equilibrium (NLTE) effects. We extract the line profiles from six transits obtained with the HARPS-N high-resolution spectrograph attached to… ▽ More We aim at extracting the transmission spectrum of the HI Balmer lines of the ultra-hot Jupiter (UHJ) KELT-20b/MASCARA-2b from observations and to further compare the results with what obtained through forward modelling accounting for non-local thermodynamic equilibrium (NLTE) effects. We extract the line profiles from six transits obtained with the HARPS-N high-resolution spectrograph attached to the Telescopio Nazionale Galileo telescope. We compute the temperature-pressure (TP) profile employing the helios code in the lower atmosphere and the Cloudy NLTE code in the middle and upper atmosphere. We further use Cloudy to compute the theoretical planetary transmission spectrum in LTE and NLTE for comparison with observations. We detected the Halpha (0.79+/-0.03%; 1.25 Rp), Hbeta (0.52+/-0.03%; 1.17 Rp), and Hgamma (0.39+/-0.06%; 1.13 Rp) lines, while we detected the Hdelta line at almost 4 sigma (0.27+/-0.07%; 1.09 Rp). The models predict an isothermal temperature of about2200 K at pressures >10^-2 bar and of about 7700 K at pressures <10^-8 bar, with a roughly linear temperature rise in between. In the middle and upper atmosphere, the NLTE TP profile is up to about 3000 K hotter than in LTE. The synthetic transmission spectrum derived from the NLTE TP profile is in good agreement with the observed HI Balmer line profiles, validating our obtained atmospheric structure. Instead, the synthetic transmission spectrum derived from the LTE TP profile leads to significantly weaker absorption compared to the observations. Metals appear to be the primary agents leading to the temperature inversion in UHJs and the impact of NLTE effects on them increases the magnitude of the inversion. We find that the impact of NLTE effects on the TP profile of KELT-20b/MASCARA-2b is larger than for the hotter UHJ KELT-9b, and thus NLTE effects might be relevant also for planets cooler than KELT-20b/MASCARA-2b. △ Less

Submitted 27 June, 2023; originally announced June 2023.

Comments: Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 676, A99 (2023)

The Impact of Cosmic Variance on Inferences of Global Neutral Fraction Derived from Ly$α$ Luminosity Functions During Reionization

Authors: Sean Bruton, Claudia Scarlata, Francesco Haardt, Matthew J. Hayes, Charlotte Mason, Alexa M. Morales, Andrei Mesinger

Abstract: We investigate the impact of field-to-field variation, deriving from cosmic variance, in measured Lyman-$α$ emitter (LAE) luminosity functions (LFs) and this variation's impact on inferences of the neutral fraction of the intergalactic medium (IGM) during reionization. We post-process a z=7 IGM simulation to populate the dark matter halos with LAEs. These LAEs have realistic UV magnitudes, Ly$α$ f… ▽ More We investigate the impact of field-to-field variation, deriving from cosmic variance, in measured Lyman-$α$ emitter (LAE) luminosity functions (LFs) and this variation's impact on inferences of the neutral fraction of the intergalactic medium (IGM) during reionization. We post-process a z=7 IGM simulation to populate the dark matter halos with LAEs. These LAEs have realistic UV magnitudes, Ly$α$ fluxes, and Ly$α$ line profiles. We calculate the attenuation of Ly$α$ emission in universes with varying IGM neutral fraction, $\bar{\rm{x}}_{\rm{HI}}$. In a $\bar{\rm{x}}_{\rm{HI}}=0.3$ simulation, we perform 100 realizations of a mock 2 square degree survey with a redshift window $Δz = 0.5$ and flux limit $\rm{f}_{Lyα}>1\times10^{-17}\:\rm{ergs}\:\: \rm{s}^{-1} \: \rm{cm}^{-2}$; such a survey is typical in depth and volume of the largest LAE surveys conducted today. For each realization, we compute the LAE LF and use it to recover the input $\bar{\rm{x}}_{\rm{HI}}$. Comparing the inferred values of $\bar{\rm{x}}_{\rm{HI}}$ across the ensemble of the surveys, we find that cosmic variance, deriving from large-scale structure and variation in the neutral gas along the sightline, imposes a floor in the uncertainty of $Δ\bar{\rm{x}}_{\rm{HI}} \sim 0.2$ when $\bar{\rm{x}}_{\rm{HI}}$ $=0.3$. We explore mitigation strategies to decrease this uncertainty, such as increasing the volume, decreasing the flux limit, or probing the volume with many independent fields. Increasing the area and/or depth of the survey does not mitigate the uncertainty, but composing a survey with many independent fields is effective. This finding highlights the best strategy for LAE surveys aiming at constraining $\bar{\rm{x}}_{\rm{HI}}$ of the universe during reionization. △ Less

Submitted 8 May, 2023; originally announced May 2023.

Comments: 17 pages, 13 figures

Decomposing galaxies with BANG: an automated morpho-kinematical decomposition of the SDSS-DR17 MaNGA survey

Authors: Fabio Rigamonti, Massimo Dotti, Stefano Covino, Francesco Haardt, Luca Cortese, Marco Landoni, Ludovica Varisco

Abstract: From a purely photometric perspective galaxies are generally decomposed into a bulge+disc system, with bulges being dispersion-dominated and discs rotationally-supported. However, recent observations have demonstrated that such a framework oversimplifies complexity, especially if one considers galaxy kinematics. To address this issue we introduced with the GPU-based code \textsc{bang} a novel appr… ▽ More From a purely photometric perspective galaxies are generally decomposed into a bulge+disc system, with bulges being dispersion-dominated and discs rotationally-supported. However, recent observations have demonstrated that such a framework oversimplifies complexity, especially if one considers galaxy kinematics. To address this issue we introduced with the GPU-based code \textsc{bang} a novel approach that employs analytical potential-density pairs as galactic components, allowing for a computationally fast, still reliable fit of the morphological and kinematic properties of galaxies. Here we apply \textsc{bang} to the SDSS-MaNGA survey, estimating key parameters such as mass, radial extensions, and dynamics, for both bulges and discs of +10,000 objects. We test our methodology against a smaller subsample of galaxies independently analysed with an orbit-based algorithm, finding agreement in the recovered total stellar mass. We also manage to reproduce well-established scaling relations, demonstrating how proper dynamical modelling can result in tighter correlations and provide corrections to standard approaches. Finally, we propose a more general way of decomposing galaxies into "hot" and "cold" components, showing a correlation with orbit-based approaches and visually determined morphological type. Unexpected tails in the "hot-to-total" mass-ratio distribution are present for galaxies of all morphologies, possibly due to visual morphology misclassifications. △ Less

Submitted 9 January, 2024; v1 submitted 5 May, 2023; originally announced May 2023.

Comments: 13 pages, 10 figures, accepted for pubblication on MNRAS

HYPerluminous quasars at the Epoch of ReionizatION (HYPERION). A new regime for the X-ray nuclear properties of the first quasars

Authors: L. Zappacosta, E. Piconcelli, F. Fiore, I. Saccheo, R. Valiante, C. Vignali, F. Vito, M. Volonteri, M. Bischetti, A. Comastri, C. Done, M. Elvis, E. Giallongo, F. La Franca, G. Lanzuisi, M. Laurenti, G. Miniutti, A. Bongiorno, M. Brusa, F. Civano, S. Carniani, V. D'Odorico, C. Feruglio, S. Gallerani, R. Gilli , et al. (18 additional authors not shown)

Abstract: The existence of luminous quasars (QSO) at the Epoch of Reionization (EoR; i.e. z>6) powered by supermassive black holes (SMBH) with masses $\gtrsim10^9~M_\odot$ challenges models of early SMBH formation. To shed light on the nature of these sources we started a multiwavelength programme based on a sample of 18 HYPerluminous quasars at the Epoch of ReionizatION (HYPERION). These are the luminous Q… ▽ More The existence of luminous quasars (QSO) at the Epoch of Reionization (EoR; i.e. z>6) powered by supermassive black holes (SMBH) with masses $\gtrsim10^9~M_\odot$ challenges models of early SMBH formation. To shed light on the nature of these sources we started a multiwavelength programme based on a sample of 18 HYPerluminous quasars at the Epoch of ReionizatION (HYPERION). These are the luminous QSOs whose SMBH must have had the fastest mass growth during the Universe first Gyr. In this paper we present the HYPERION sample and report on the first of the 3 years planned observations of the 2.4 Ms XMM-Newton Multi-Year Heritage program on which HYPERION is based. The goal of this program is to accurately characterize the X-ray nuclear properties of QSOs at the EoR. Through a joint X-ray spectral analysis of 10 sources, in the rest-frame $\sim2-50$ keV range, we report a steep average photon index ($Γ\sim2.4\pm0.1$). Absorption is not required. The average $Γ$ is inconsistent at $\geq4σ$ level with the canonical 1.8-2 value measured in QSO at z<6. This spectral slope is also much steeper than that reported in lower-z QSOs with similar luminosity or accretion rate, thus suggesting a genuine redshift evolution. Alternatively, we can interpret this result as the presence of an unusually low-energy cutoff $E_{cut}\sim20$ keV on a standard $Γ=1.9$ power-law. We also report on mild indications that HYPERION QSOs show higher soft X-ray emission at 2 keV compared to the UV one at 2500A than expected by lower-z luminous AGN. We speculate that a redshift-dependent coupling between the corona and accretion disc or intrinsically different coronal properties may account for the steep spectral slopes, especially in the presence of powerful winds. The reported slopes, if confirmed at lower luminosities, may have an important impact on future X-ray AGN studies in the early Universe. △ Less

Submitted 19 July, 2023; v1 submitted 3 May, 2023; originally announced May 2023.

Comments: 21 pages (including appendix), 12 figures, 4 tables. Accepted for pubblication in A&A

Journal ref: A&A 678, A201 (2023)

The ultraviolet habitable zone of exoplanets

Authors: Riccardo Spinelli, Francesco Borsa, Giancarlo Ghirlanda, Gabriele Ghisellini, Francesco Haardt

Abstract: The dozens of rocky exoplanets discovered in the Circumstellar Habitable Zone (CHZ) currently represent the most suitable places to host life as we know it outside the Solar System. However, the presumed presence of liquid water on the CHZ planets does not guarantee suitable environments for the emergence of life. According to experimental studies, the building blocks of life are most likely produ… ▽ More The dozens of rocky exoplanets discovered in the Circumstellar Habitable Zone (CHZ) currently represent the most suitable places to host life as we know it outside the Solar System. However, the presumed presence of liquid water on the CHZ planets does not guarantee suitable environments for the emergence of life. According to experimental studies, the building blocks of life are most likely produced photochemically in presence of a minimum ultraviolet (UV) flux. On the other hand, high UV flux can be life-threatening, leading to atmospheric erosion and damaging biomolecules essential to life. These arguments raise questions about the actual habitability of CHZ planets around stars other than Solar-type ones, with different UV to bolometric luminosity ratios. By combining the "principle of mediocricy" and recent experimental studies, we define UV boundary conditions (UV-habitable Zone, UHZ) within which life can possibly emerge and evolve. We investigate whether exoplanets discovered in CHZs do indeed experience such conditions. By analysing Swift-UV/Optical Telescope data, we measure the near ultraviolet (NUV) luminosities of 17 stars harbouring 23 planets in their CHZ. We derive an empirical relation between NUV luminosity and stellar effective temperature. We find that eighteen of the CHZ exoplanets actually orbit outside the UHZ, i.e., the NUV luminosity of their M-dwarf hosts is decisively too low to trigger abiogenesis - through cyanosulfidic chemistry - on them. Only stars with effective temperature >3900 K illuminate their CHZ planets with enough NUV radiation to trigger abiogenesis. Alternatively, colder stars would require a high-energy flaring activity. △ Less

Submitted 28 March, 2023; originally announced March 2023.

Comments: 8 pages, 3 figures

Gravitational waves from an eccentric population of primordial black holes orbiting Sgr A$^{\star}$

Authors: Stefano Bondani, Matteo Bonetti, Luca Broggi, Francesco Haardt, Alberto Sesana, Massimo Dotti

Abstract: Primordial black holes (PBH), supposedly formed in the very early Universe, have been proposed as a possible viable dark matter candidate. In this work we characterize the expected gravitational wave (GW) losses from a population of PBHs orbiting Sgr A$^{\star}$, the super-massive black hole at the Galactic center (GC), and assess the signal detectability by the planned space-borne interferometer… ▽ More Primordial black holes (PBH), supposedly formed in the very early Universe, have been proposed as a possible viable dark matter candidate. In this work we characterize the expected gravitational wave (GW) losses from a population of PBHs orbiting Sgr A$^{\star}$, the super-massive black hole at the Galactic center (GC), and assess the signal detectability by the planned space-borne interferometer LISA and by the proposed next generation space-borne interferometer $μ$Ares. Assuming that PBHs indeed form the entire diffuse mass allowed to reside within the orbit of the S2 star, we compute an upper limit to the expected GW signal both from resolved and non-resolved sources, under the further assumptions of monochromatic mass function and thermally distributed eccentricities. By comparing with our previous work where PBHs on circular orbits were assumed, we show for 1 M$_{\odot}$ PBHs how the GW signal from high harmonics over a 10 year data stream increases by a factor of six the chances of LISA detectability, from the $\approx 10\%$ of the circular case, to $\approx 60\%$, whereas multiple sources can be identified in $20\%$ of our mock populations. The background signal, made by summing up all non resolved sources, should be certainly detectable thanks to the PBHs with higher eccentricity evolving under two body relaxation. In the case of $μ$Ares, because of its improved sensitivity in the $μ$Hz band, one third of the entire population of PBHs orbiting Sgr A$^{\star}$ would be resolved. The background noise from the remaining non resolved sources should be detectable as well. Finally we present the results for different PBH masses. △ Less

Submitted 19 February, 2024; v1 submitted 22 March, 2023; originally announced March 2023.

Comments: 11 pages, 9 figures

Journal ref: Phys. Rev. D 109, 043005 (2024)

Eddington accreting Black Holes in the Epoch of Reionization

Authors: Fabio Fontanot, Stefano Cristiani, Andrea Grazian, Francesco Haardt, Valentina D'Odorico, Konstantina Boutsia, Giorgio Calderone, Guido Cupani, Francesco Guarneri, Chiara Fiorin, Giulia Rodighiero

Abstract: The evolution of the luminosity function (LF) of Active Galactic Nuclei (AGNs) at $z \gtrsim 5$ represents a key constraint to understand their contribution to the ionizing photon budget necessary to trigger the last phase transition in the Universe, i.e. the epoch of Reionization. Recent searches for bright high-z AGNs suggest that the space densities of this population at $z>4$ has to be revised… ▽ More The evolution of the luminosity function (LF) of Active Galactic Nuclei (AGNs) at $z \gtrsim 5$ represents a key constraint to understand their contribution to the ionizing photon budget necessary to trigger the last phase transition in the Universe, i.e. the epoch of Reionization. Recent searches for bright high-z AGNs suggest that the space densities of this population at $z>4$ has to be revised upwards, and sparks new questions about their evolutionary paths. Gas accretion is the key physical mechanism to understand both the distribution of luminous sources and the growth of central Super-Massive Black Holes (SMBHs). In this work, we model the high-z AGN-LF assuming that high-z luminous AGN shine at their Eddington limit: we derive the expected evolution as a function of the ``duty-cycle'' ($f_{\rm dc}$), i.e. the fraction of life-time that a given SMBH spends accreting at the Eddington rate. Our results show that intermediate values ($f_{\rm dc} \simeq 0.1$) predict the best agreement with the ionizing background and photoionization rate, but do not provide enough ionizing photons to account for the observed evolution of the hydrogen neutral fraction. Smaller values ($f_{\rm dc} \lesssim 0.05$) are required for AGNs to be the dominant population responsible for Hydrogen reionization in the Early Universe. We then show that this low-$f_{\rm dc}$ evolution can be reconciled with the current constraints on Helium reionization, although it implies a relatively large number of inactive SMBHs at $z\gtrsim5$, in tension with SMBH growth models based on heavy seeding. △ Less

Submitted 17 January, 2023; originally announced January 2023.

Comments: 11 pages, 5 figures, MNRAS accepted

Dynamical evolution of massive black hole pairs in the presence of spin-dependent radiative feedback

Authors: Francesco Bollati, Alessandro Lupi, Massimo Dotti, Francesco Haardt

Abstract: The putative ubiquity of massive black holes (MBH) at the center of galaxies, and the hierarchical progress of structure formation along the cosmic history, together necessarily imply the existence of a large population of cosmic MBH binaries. Such systems are understood to be the loudest sources of gravitational waves at mHz frequencies, the regime that will be probed by the next Laser Interferom… ▽ More The putative ubiquity of massive black holes (MBH) at the center of galaxies, and the hierarchical progress of structure formation along the cosmic history, together necessarily imply the existence of a large population of cosmic MBH binaries. Such systems are understood to be the loudest sources of gravitational waves at mHz frequencies, the regime that will be probed by the next Laser Interferometer Space Antenna (LISA). It has been proposed that the rate at which MBHs pair and then bind to form binaries is critically dependent upon the feedback exerted by the MBHs on the surrounding gaseous environment. Using the publicly available code GIZMO, we perform a suite of simulations aimed at studying the dynamics of a MBH pair embedded in a gaseous disk on 100 pc scale. By means of dedicated modules, we follow the dynamics of MBHs in the presence of different spin-dependent radiative feedback models, and compare the results to a benchmark case with no feedback at all. Our main finding is that feedback causes the secondary MBH to shrink its orbit at a reduced pace, when compared to models where feedback is absent. Moreover, such slower inspiral occurs on eccentric orbits, as feedback has the net effect of hampering the circularization process. Though idealized in many aspects, our study highlights and quantifies the importance of including spin-dependent feedback recipes in hydrodynamic simulations of MBH pairs, and ultimately in assessing the cosmological coalescence rate of such systems in view of their detection through gravitational waves. △ Less

Submitted 16 December, 2022; originally announced December 2022.

Planetary parameters, XUV environments and mass-loss rates for nearby gaseous planets with X-ray detected host-stars

Authors: R. Spinelli, E. Gallo, F. Haardt, A. Caldiroli, F. Biassoni, F. Borsa, E. Rauscher

Abstract: We leverage Gaia DR2 parallactic distances to deliver new or revised estimates of planetary parameters and X-ray irradiation for a distance-limited ($\lesssim 100$ pc) sample of 27 gaseous planets (from super-Earths to hot Jupiters) with publicly available Chandra and/or XMM observations, for which we carry out a homogeneous data reduction. For 20 planets with X-ray detected host stars we make use… ▽ More We leverage Gaia DR2 parallactic distances to deliver new or revised estimates of planetary parameters and X-ray irradiation for a distance-limited ($\lesssim 100$ pc) sample of 27 gaseous planets (from super-Earths to hot Jupiters) with publicly available Chandra and/or XMM observations, for which we carry out a homogeneous data reduction. For 20 planets with X-ray detected host stars we make use of the photoionization hydrodynamics code ATES to derive updated atmospheric mass outflow rates. The newly derived masses/radii are not consistent with the exoplanet.eu values for five systems; HD 149026b and WASP-38, for mass; and Au Mic b, HAT-P-20 and HAT-P-2 for radii. Notably, the lower mass implies a (Saturn-like) density of $0.86\pm 0.09$ g cm$^{-3}$ for HD 149026 b. This independent estimate is consistent with the lowest values reported in the literature. Separately, we report on the X-ray detection of GJ 9827, HD 219134 and LHS 1140 for the first time. The inferred stellar X-ray luminosity of LHS 1140 ($1.34^{+0.19}_{-0.21} \times 10^{26}$ erg sec$^{-1}$) implies that LHS 1140 b is the least irradiated transiting super-Earth known to orbit within the habitable zone of a nearby M-dwarf. △ Less

Submitted 7 April, 2023; v1 submitted 2 August, 2022; originally announced August 2022.

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

Subarcsecond view on the high-redshift blazar GB 1508+5714 by the International LOFAR Telescope

Authors: A. Kappes, P. R. Burd, M. Kadler, G. Ghisellini, E. Bonnassieux, M. Perucho, M. Brüggen, C. C. Cheung, B. Ciardi, E. Gallo, F. Haardt, L. K. Morabito, T. Sbarrato, A. Drabent, J. Harwood, N. Jackson, J. Moldon

Abstract: Studies of the most distant AGNs allow us to test our current understanding of the physics present in radio-jetted AGNs across a range of environments. The decrease in apparent luminosity with distance is the primary difficulty to overcome in the study of these distant AGNs, which requires highly sensitive instruments. Our goal is to employ new long wavelength radio data to better parametrise the… ▽ More Studies of the most distant AGNs allow us to test our current understanding of the physics present in radio-jetted AGNs across a range of environments. The decrease in apparent luminosity with distance is the primary difficulty to overcome in the study of these distant AGNs, which requires highly sensitive instruments. Our goal is to employ new long wavelength radio data to better parametrise the broad-band SED of GB 1508+5714, a high-redshift (z=4.30) AGN. Its high redshift, high intrinsic luminosity, and classification as a blazar allow us to test emission models that consider the efficient cooling of jet electrons via inverse Compton losses in interactions with the dense CMB photon field at high redshifts. A significant detection of this effect in GB 1508+5714 may partly explain the apparent sparsity of high-redshift radio galaxies in wide-field surveys; detections of this kind are only becoming possible with the current generation of SKA precursors. We used international LOFAR telescope to image the long wavelength radio emission around the high-redshift blazar GB 1508+5714 on arcsecond scales at frequencies between 128 MHz and 160 MHz. This allowed us to compare the spatially resolved structure with higher frequency observations, and to construct spectral index maps. The LOFAR image shows a compact unresolved core and two resolved emission regions around 2 arcsec to the east and to the west of the radio core. We find structure consistent with previous VLA observations, as well as a previously unreported emission region to the east. We interpret the arcsecond-scale radio structure of GB 1508+5714 as a FR II-like radio galaxy at a small viewing angle. Our SED modelling shows that a scenario featuring significant quenching effects caused by interaction with the CMB provides a good description of the data, and notably explains the suppressed radio emission. △ Less

Submitted 23 May, 2022; originally announced May 2022.

Comments: 11 pages, 10 figures, 2 tables

MSC Class: 85-02

Journal ref: A&A 663, A44 (2022)

Maximally informed Bayesian modelling of disc galaxies

Authors: Fabio Rigamonti, Massimo Dotti, Stefano Covino, Francesco Haardt, Marco Landoni, Walter Del Pozzo, Alessandro Lupi, Stefano Zibetti

Abstract: Dissecting the underlying structure of galaxies is of main importance in the framework of galaxy formation and evolution theories. While a classical bulge+disc decomposition of disc galaxies is usually taken as granted, this is only rarely solidly founded upon the full exploitation of the richness of data arising from spectroscopic studies with integral field units. In this work we describe a full… ▽ More Dissecting the underlying structure of galaxies is of main importance in the framework of galaxy formation and evolution theories. While a classical bulge+disc decomposition of disc galaxies is usually taken as granted, this is only rarely solidly founded upon the full exploitation of the richness of data arising from spectroscopic studies with integral field units. In this work we describe a fully Bayesian estimation method of the global structure of disc galaxies which makes use of the wealth of photometric, kinematic, and mass-to-light ratio data, and that can be seen as a first step towards a machine-learning approach, certainly needed when dealing with larger samples of galaxies. Ours is a novel, hybrid line of action in tackling the problem of galactic parameter estimation, neither purely photometric nor orbit-based. Being rooted on a nested sampler, our code, which is available publicly as an online repository, allows for a statistical assessment of the need for multiple components in the dissecting process. As a first case-study the GPU-optimized code is applied to the S0 galaxy NGC-7683, finding that in this galaxy a pseudo-bulge, possibly the remnant of a bar-like structure, does exist in the center of the system. These results are then tested against the publicly available, orbit-based code DYNAMITE, finding substantial agreement. △ Less

Submitted 17 May, 2022; v1 submitted 13 May, 2022; originally announced May 2022.

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

Re-weighting dark matter in disc galaxies: a new general relativistic observational test

Authors: Davide Astesiano, Sergio L. Cacciatori, Massimo Dotti, Francesco Haardt, Federico Re

Abstract: A recent analysis of data from the ESA Gaia mission demonstrated that the kinematics of stars in the Milky Way can be modelled without invoking the presence of dark matter whatsoever. Indeed, the higher-than-Keplerian velocities observed in outer stars can be ascribed to the properties of the general relativistic (GR) metric assumed to describe the Galaxy. Here we generalize the concept, and deriv… ▽ More A recent analysis of data from the ESA Gaia mission demonstrated that the kinematics of stars in the Milky Way can be modelled without invoking the presence of dark matter whatsoever. Indeed, the higher-than-Keplerian velocities observed in outer stars can be ascribed to the properties of the general relativistic (GR) metric assumed to describe the Galaxy. Here we generalize the concept, and derive the most general exact GR model for a stationary axi-symmetric dynamically cold dust structure. We explicitly show how deviations from the commonly adopted Newtonian dynamics are indeed manifest even at low velocities and low densities. We provide for the first time a detailed description of the frequency shift experienced by photons travelling from any emission site within an external disc galaxy to the observer, relating the outcome of the shift measurement to the gravitational properties of the galaxy. Finally, we devise a novel, groundbreaking observational test potentially able to fully characterize the GR metric under the minimal set of assumptions mentioned above. The proposed experiment exploits the effects non-diagonal GR terms have on the frequency shift of photons, ultimately providing a test to evaluate whether dark matter is actually required by disc galaxy kinematics. △ Less

Submitted 11 April, 2022; originally announced April 2022.

Comments: 7 pages, 0 figures

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)

On the detectability of gravitational waves from primordial black holes orbiting Sgr A*

Authors: Stefano Bondani, Francesco Haardt, Alberto Sesana, Enrico Barausse, Massimo Dotti

Abstract: In this work we characterize the expected gravitational wave signal detectable by the planned space-borne interferometer LISA and the proposed next generation space-borne interferometer $μ$Ares arising from a population of primordial black holes orbiting Sgr A*, the super-massive black hole at the Galactic center. Assuming that such objects indeed form the entire diffuse mass allowed by the observ… ▽ More In this work we characterize the expected gravitational wave signal detectable by the planned space-borne interferometer LISA and the proposed next generation space-borne interferometer $μ$Ares arising from a population of primordial black holes orbiting Sgr A*, the super-massive black hole at the Galactic center. Assuming that such objects indeed form the entire diffuse mass allowed by the observed orbit of S2 in the Galactic center, under the simplified assumption of circular orbits and monochromatic mass function, we assess the expected signal in gravitational waves, either from resolved and non-resolved sources. We estimate a small but non negligible chance of $\simeq$ 10% of detecting one single 1 M$_{\odot}$ primordial black hole with LISA in a 10-year-long data stream, while the background signal due to unresolved sources would essentially elude any reasonable chance of detection. On the contrary, $μ$Ares, with a $\simeq$ 3 orders-of-magnitude better sensitivity at $\simeq$ 10$^{-5}$ Hz, would be able to resolve $\simeq$ 140 solar mass primordial black holes in the same amount of time, while the unresolved background should be observable with an integrated signal-to-noise ratio $\gtrsim$ 100. Allowing the typical PBH mass to be in the range 0.01-10 M$_{\odot}$ would increase LISA chance of detection to $\simeq$ 40% towards the lower limit of the mass spectrum. In the case of $μ$Ares, instead, we find a "sweet spot" just about 1 M$_{\odot}$, a mass for which the number of resolvable events is indeed maximized. △ Less

Submitted 2 August, 2022; v1 submitted 10 March, 2022; originally announced March 2022.

Comments: 10 pages, 11 figures; Referee reviewed

Journal ref: Phys. Rev. D 106, 043015, 2022

Misaligned Spinning Binary Black Hole Mergers in Hot Magnetized Plasma

Authors: Federico Cattorini, Sofia Maggioni, Bruno Giacomazzo, Francesco Haardt, Monica Colpi, Stefano Covino

Abstract: We present general relativistic magneto-hydrodynamical simulations of equal-mass spinning black hole binary mergers embedded in a magnetized gas cloud. We focus on the effect of the spin orientation relative to the orbital angular momentum on the flow dynamics, mass accretion rate and Poynting luminosity. We find that, across the inspiral, the gas accreting onto the individual black holes concentr… ▽ More We present general relativistic magneto-hydrodynamical simulations of equal-mass spinning black hole binary mergers embedded in a magnetized gas cloud. We focus on the effect of the spin orientation relative to the orbital angular momentum on the flow dynamics, mass accretion rate and Poynting luminosity. We find that, across the inspiral, the gas accreting onto the individual black holes concentrates into disk-like overdensities, whose angular momenta are oriented towards the spin axes and which persist until merger. We identify quasi-periodic modulations occurring in the mass accretion rate at the level of 1-20%, evolving in parallel with the gravitational wave chirp. The similarity between the accretion rate time-series and the gravitational strain is a consequence of the interplay between strong, dynamical gravitational fields and magnetic fields in the vicinity of the inspiralling black holes. This result suggests that quasi-periodicity in the pre-merger accretion rate of massive binaries is not exclusive of environments in which the black holes are embedded in a circumbinary accretion disk, and could provide an additional useful signature of electromagnetic emission concurrent to low-frequency gravitational wave detection. △ Less

Submitted 4 May, 2022; v1 submitted 16 February, 2022; originally announced February 2022.

Comments: 9 pages, 4 figures. Accepted for publication in The Astrophysical Journal Letters

Irradiation-driven escape of primordial planetary atmospheres II. Evaporation efficiency of sub-Neptunes through hot Jupiters

Authors: Andrea Caldiroli, Francesco Haardt, Elena Gallo, Riccardo Spinelli, Isaac Malsky, Emily Rauscher

Abstract: Making use of the publicly available 1D photoionization hydrodynamics code ATES we set out to investigate the combined effects of planetary gravitational potential energy ($φ_p\equiv GM_p/R_p$) and stellar X-ray and Extreme Ultraviolet (XUV) irradiation ($F_{\rm XUV}$) on the evaporation efficiency ($η$) of moderately-to-highly irradiated gaseous planets, from sub-Neptunes through hot Jupiters. We… ▽ More Making use of the publicly available 1D photoionization hydrodynamics code ATES we set out to investigate the combined effects of planetary gravitational potential energy ($φ_p\equiv GM_p/R_p$) and stellar X-ray and Extreme Ultraviolet (XUV) irradiation ($F_{\rm XUV}$) on the evaporation efficiency ($η$) of moderately-to-highly irradiated gaseous planets, from sub-Neptunes through hot Jupiters. We show that the (known) existence of a threshold potential above which energy-limited escape (i.e., $η\simeq 1$) is unattainable can be inferred analytically. For $\log φ_p\gtrsim \log φ_p^{\rm thr}\approx [12.9-13.2]$ (in cgs units), most of the energy absorption occurs where the average kinetic energy acquired by the ions through photo-electron collisions is insufficient for escape. This causes the evaporation efficiency to plummet with increasing $φ_p$,. Whether or not planets with $φ_p\lesssim φ_p^{\rm thr}$ exhibit energy-limited outflows is regulated primarily by the stellar irradiation level. Specifically, for low-gravity planets, above $F_{\rm XUV}\simeq 10^{4-5}$ erg cm$^{-2}$s$^{-1}$ Ly$α$ losses overtake adiabatic and advective cooling and the evaporation efficiency of low-gravity planets drops below the energy-limited approximation, albeit remaining largely independent of $φ_p$Further, we show that whereas $η$ increases as $F_{\rm XUV}$ increases for planets above $φ^{\rm thr}_p$, the opposite is true for low-gravity planets. This behavior can be understood by examining the relative fractional contributions of advective and radiative losses as a function of atmospheric temperature. This novel framework enables a reliable, physically motivated prediction of the expected evaporation efficiency for a given planetary system; an analytical approximation of the best-fitting $η$ is given in the appendix. △ Less

Submitted 6 April, 2022; v1 submitted 1 December, 2021; originally announced December 2021.

Comments: 12 pages, 8 figures. Accepted for publication in A&A. In the new version, section 4 and 5 have been completely revised after the referee revision. Comments are welcome

Journal ref: A&A 663, A122 (2022)

Irradiation-driven escape of primordial planetary atmospheres I. The ATES photoionization hydrodynamics code

Authors: Andrea Caldiroli, Francesco Haardt, Elena Gallo, Riccardo Spinelli, Isaac Malsky, Emily Rauscher

Abstract: Intense X-ray and ultraviolet stellar irradiation can heat and inflate the atmospheres of closely orbiting exoplanets, driving mass outflows that may be significant enough to evaporate a sizable fraction of the planet atmosphere over the system lifetime. The recent surge in the number of known exoplanets, together with the imminent deployment of new ground and space-based facilities for exoplanet… ▽ More Intense X-ray and ultraviolet stellar irradiation can heat and inflate the atmospheres of closely orbiting exoplanets, driving mass outflows that may be significant enough to evaporate a sizable fraction of the planet atmosphere over the system lifetime. The recent surge in the number of known exoplanets, together with the imminent deployment of new ground and space-based facilities for exoplanet discovery and characterization, requires a prompt and efficient assessment of the most promising targets for intensive spectroscopic follow-ups. To this purpose, we developed ATES (ATmospheric EScape); a new hydrodynamics code that is specifically designed to compute the temperature, density, velocity and ionization fraction profiles of highly irradiated planetary atmospheres, along with the current, steady-state mass loss rate. ATES solves the one-dimensional Euler, mass and energy conservation equations in radial coordinates through a finite-volume scheme. The hydrodynamics module is paired with a photoionization equilibrium solver that includes cooling via bremsstrahlung, recombination and collisional excitation/ionization for the case of a primordial atmosphere entirely composed of atomic hydrogen and helium, whilst also accounting for advection of the different ion species. Compared against the results of 14 moderately-to-highly irradiated planets simulated with The PLUTO-CLOUDY Interface (TPCI), ATES yields remarkably good agreement at a significantly smaller fraction of the computational time. A convergence study shows that ATES recovers stable, steady-state hydrodynamic solutions for systems with $\log(-φ_p) \lesssim 12.9 + 0.17\log F_{\rm XUV}$. Incidentally, atmospheres of systems above this threshold are generally thought to be undergoing Jeans escape. The code, which also features a user-friendly graphic interface, is available publicly as an online repository. △ Less

Submitted 16 August, 2021; v1 submitted 18 June, 2021; originally announced June 2021.

Comments: 17 pages, 11 figures. Revised version. Accepted by A&A. Minor changes only; this new version includes a convergence study. Comments are welcome

Journal ref: A&A 655, A30 (2021)

Proof of CMB-driven X-ray brightening of high-z radio galaxies

Authors: Edmund Hodges-Kluck, Elena Gallo, Gabriele Ghisellini, Francesco Haardt, Jianfeng Wu, Benedetta Ciardi

Abstract: We present a definitive assessment of the role of Inverse Compton scattering of Cosmic Microwave Background photons (IC/CMB) in the context of radio galaxies. Owing to the steep increase of the CMB radiation energy density, IC/CMB is supposed to become progressively more important with respect to radio synchrotron cooling as the redshift increases. For typical energies at play, this process will u… ▽ More We present a definitive assessment of the role of Inverse Compton scattering of Cosmic Microwave Background photons (IC/CMB) in the context of radio galaxies. Owing to the steep increase of the CMB radiation energy density, IC/CMB is supposed to become progressively more important with respect to radio synchrotron cooling as the redshift increases. For typical energies at play, this process will up-scatter the CMB photons into the X-ray band, and is thus expected to yield a redshift-dependent, concurrent X-ray brightening and radio dimming of the jet-powered structures. Here we show how a conclusive proof of this effect hinges on high-resolution imaging data in which the extended lobes can be distinguished from the compact hot spots where synchrotron-self-Compton dominates the X-ray emission regardless of redshift. We analyze Chandra and Very Large Array data of 11 radio galaxies between 1.3<z<4.3, and demonstrate that the emission from their lobes is fully consistent with the expectations from IC/CMB in equipartition. Once the dependence on size and radio luminosity are properly accounted for, the measured lobe X-ray luminosities bear the characteristic (1+z)**4 proportionality expected of a CMB seed radiation field. Whereas this effect can effectively quench the (rest-frame) GHz radio emission from z>3 radio galaxies below <1 mJy, IC/CMB alone can not be responsible for a deficit in high-z, radio-loud AGN if--as we argue--such AGN typically have bright, compact hot spots. △ Less

Submitted 7 May, 2021; originally announced May 2021.

Comments: 16 pages, 5 figures, 5 tables. Accepted for publication in MNRAS

Intermediate mass black holes from stellar mergers in young star clusters

Authors: Ugo N. Di Carlo, Michela Mapelli, Mario Pasquato, Sara Rastello, Alessandro Ballone, Marco Dall'Amico, Nicola Giacobbo, Giuliano Iorio, Mario Spera, Stefano Torniamenti, Francesco Haardt

Abstract: Intermediate mass black holes (IMBHs) in the mass range $10^2-10^5\,\mathrm{M_{\odot}}$ bridge the gap between stellar black holes (BHs) and supermassive BHs. Here, we investigate the possibility that IMBHs form in young star clusters via runaway collisions and BH mergers. We analyze $10^4$ simulations of dense young star clusters, featuring up-to-date stellar wind models and prescriptions for cor… ▽ More Intermediate mass black holes (IMBHs) in the mass range $10^2-10^5\,\mathrm{M_{\odot}}$ bridge the gap between stellar black holes (BHs) and supermassive BHs. Here, we investigate the possibility that IMBHs form in young star clusters via runaway collisions and BH mergers. We analyze $10^4$ simulations of dense young star clusters, featuring up-to-date stellar wind models and prescriptions for core collapse and (pulsational) pair instability. In our simulations, only 9 IMBHs out of 218 form via binary BH mergers, with a mass $\sim{}100-140$ M$_\odot$. This channel is strongly suppressed by the low escape velocity of our star clusters. In contrast, IMBHs with masses up to $\sim{}438$ M$_{\odot}$ efficiently form via runaway stellar collisions, especially at low metallicity. Up to $\sim{}0.2$~% of all the simulated BHs are IMBHs, depending on progenitor's metallicity. The runaway formation channel is strongly suppressed in metal-rich ($Z=0.02$) star clusters, because of stellar winds. IMBHs are extremely efficient in pairing with other BHs: $\sim{}70$% of them are members of a binary BH at the end of the simulations. However, we do not find any IMBH-BH merger. More massive star clusters are more efficient in forming IMBHs: $\sim{}8$% ($\sim{}1$%) of the simulated clusters with initial mass $10^4-3\times{}10^4$ M$_\odot$ ($10^3-5\times{}10^3$ M$_\odot$) host at least one IMBH. △ Less

Submitted 3 May, 2021; originally announced May 2021.

Comments: 12 pages, 16 figures, 2 tables. Comments welcome

Fully General Relativistic Magnetohydrodynamic Simulations of Accretion Flows onto Spinning Massive Black Hole Binary Mergers

Authors: Federico Cattorini, Bruno Giacomazzo, Francesco Haardt, Monica Colpi

Abstract: We perform the first suite of fully general relativistic magnetohydrodynamic simulations of spinning massive black hole binary mergers. We consider binary black holes with spins of different magnitudes aligned to the orbital angular momentum, which are immersed in a hot, magnetized gas cloud. We investigate the effect of the spin and degree of magnetization (defined through the fluid parameter… ▽ More We perform the first suite of fully general relativistic magnetohydrodynamic simulations of spinning massive black hole binary mergers. We consider binary black holes with spins of different magnitudes aligned to the orbital angular momentum, which are immersed in a hot, magnetized gas cloud. We investigate the effect of the spin and degree of magnetization (defined through the fluid parameter $β^{-1}\equiv p_{\mathrm{mag}}/p_{\mathrm{fluid}}$) on the properties of the accretion flow. We find that magnetized accretion flows are characterized by more turbulent dynamics, as the magnetic field lines are twisted and compressed during the late inspiral. Post-merger, the polar regions around the spin axis of the remnant Kerr black hole are magnetically dominated, and the magnetic field strength is increased by a factor $\sim$10$^2$ (independently from the initial value of $β^{-1}$). The magnetized gas in the equatorial plane acquires higher angular momentum, and settles in a thin circular structure around the black hole. We find that mass accretion rates of magnetized configurations are generally smaller than in the unmagnetized cases by up to a factor $\sim$3. Black hole spins have also a suppressing effect on the accretion rate, as large as $\sim$48\%. As a potential driver for electromagnetic emission we follow the evolution of the Poynting luminosity, which increases after merger up to a factor $\sim2$ with increasing spin, regardless of the initial level of magnetization of the fluid. Our results stress the importance of taking into account both spins and magnetic fields when studying accretion processes onto merging massive black holes. △ Less

Submitted 3 June, 2021; v1 submitted 25 February, 2021; originally announced February 2021.

Comments: 15 pages, 12 figures

Journal ref: Phys. Rev. D 103, 103022 (2021)

The Evolution of the Lyman-Alpha Luminosity Function During Reionization

Authors: Alexa Morales, Charlotte Mason, Sean Bruton, Max Gronke, Francesco Haardt, Claudia Scarlata

Abstract: The time frame in which hydrogen reionization occurred is highly uncertain, but can be constrained by observations of Lyman-alpha (Ly$α$) emission from distant sources. Neutral hydrogen in the intergalactic medium (IGM) attenuates Ly$α$~photons emitted by galaxies. As reionization progressed the IGM opacity decreased, increasing Ly$α$~visibility. The galaxy Ly$α$~luminosity function (LF) is thus a… ▽ More The time frame in which hydrogen reionization occurred is highly uncertain, but can be constrained by observations of Lyman-alpha (Ly$α$) emission from distant sources. Neutral hydrogen in the intergalactic medium (IGM) attenuates Ly$α$~photons emitted by galaxies. As reionization progressed the IGM opacity decreased, increasing Ly$α$~visibility. The galaxy Ly$α$~luminosity function (LF) is thus a useful tool to constrain the timeline of reionization. In this work, we model the Ly$α$~LF as a function of redshift, $z=5-10$, and average IGM neutral hydrogen fraction, $\overline{x}_\textsc{hi}$. We combine the Ly$α$~luminosity probability distribution obtained from inhomogeneous reionization simulations with a model for the UV LF to model the Ly$α$~LF. As the neutral fraction increases, the average number density of Ly$α$~emitting galaxies decreases, and are less luminous, though for $\overline{x}_\textsc{hi} \lesssim 0.4$ there is only a small decrease of the Ly$α$~LF. We use our model to infer the IGM neutral fraction at $z=6.6, 7.0, 7.3$ from observed Ly$α$~LFs. We conclude that there is a significant increase in the neutral fraction with increasing redshift: $\overline{x}_\textsc{hi}(z=6.6)=0.08^{+ 0.08}_{- 0.05}, \, \overline{x}_\textsc{hi}(z=7.0)=0.28 \pm 0.05$ and $\overline{x}_\textsc{hi}(z=7.3)=0.83^{+ 0.06}_{- 0.07}$. We predict trends in the Ly$α$~luminosity density and Schechter parameters as a function of redshift and the neutral fraction. We find that the Ly$α$~luminosity density decreases as the universe becomes more neutral. Furthermore, as the neutral fraction increases, the faint-end slope of the Ly$α$~LF steepens, and the characteristic Ly$α$~luminosity shifts to lower values, concluding that the evolving shape of the Ly$α$~LF -- not just its integral -- is an important tool to study reionization. △ Less

Submitted 2 July, 2021; v1 submitted 4 January, 2021; originally announced January 2021.

Comments: 20 pages, 10 figures, submitted to ApJ

Unveiling early black hole growth with multi-frequency gravitational wave observations

Authors: Rosa Valiante, Monica Colpi, Raffaella Schneider, Alberto Mangiagli, Matteo Bonetti, Giulia Cerini, Stephen Fairhurst, Francesco Haardt, Cameron Mills, Alberto Sesana

Abstract: Third Generation ground based Gravitational Wave Interferometers, like the Einstein Telescope (ET), Cosmic Explorer (CE), and the Laser Interferometer Space Antenna (LISA) will detectcoalescing binary black holes over a wide mass spectrum and across all cosmic epochs. We track the cosmological growth of the earliest light and heavy seeds that swiftly transit into the supermassive domain using a se… ▽ More Third Generation ground based Gravitational Wave Interferometers, like the Einstein Telescope (ET), Cosmic Explorer (CE), and the Laser Interferometer Space Antenna (LISA) will detectcoalescing binary black holes over a wide mass spectrum and across all cosmic epochs. We track the cosmological growth of the earliest light and heavy seeds that swiftly transit into the supermassive domain using a semi analytical model for the formation of quasars at $z=6.4$, 2 and $0.2$, in which we follow black hole coalescences driven by triple interactions. We find that light seed binaries of several $10^2$ M$_\odot$ are accessible to ET with a signal-to-noise ratio ($S/N$) of $10-20$ at $6<z<15$. They then enter the LISA domain with larger $S/N$ as they grow toa few $10^4$ M$_\odot$. Detecting their gravitational signal would provide first time evidence that light seeds form, grow and dynamically pair during galaxy mergers. The electromagnetic emission of accreting black holes of similar mass and redshift is too faint to be detected even for the deepest future facilities. ET will be our only chance to discover light seeds forming at cosmicdawn. At $2<z<8$, we predict a population of "starved binaries", long-lived marginally-growing light seed pairs, to be loud sources in the ET bandwidth ($S/N>20$). Mergers involving heavy seeds ($\sim 10^5 M_\odot - 10^6 M_\odot$) would be within reach up to $z=20$ in the LISA frequency domain. The lower-z model predicts $11.25(18.7)$ ET(LISA) events per year, overall. △ Less

Submitted 28 October, 2020; originally announced October 2020.

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

The best place and time to live in the Milky Way

Authors: R. Spinelli, G. Ghirlanda, F. Haardt, G. Ghisellini, G. Scuderi

Abstract: Counted among the most powerful cosmic events, supernovae (SNe) and gamma-ray bursts (GRBs) can be highly disruptive for life: Their radiation can be harmful for biota or induce extinction by removing most of the protective atmospheric ozone layer from terrestrial planets (TPs). Nearby high-energy transient astrophysical events have been proposed as possible triggers of mass extinctions on Earth.… ▽ More Counted among the most powerful cosmic events, supernovae (SNe) and gamma-ray bursts (GRBs) can be highly disruptive for life: Their radiation can be harmful for biota or induce extinction by removing most of the protective atmospheric ozone layer from terrestrial planets (TPs). Nearby high-energy transient astrophysical events have been proposed as possible triggers of mass extinctions on Earth. We assess the habitability of the Milky Way (MW) throughout its cosmic history against potentially disruptive astrophysical transients with the aim of identifying the safest places and epochs within our Galaxy. We also test the hypothesis that one long GRB played a leading role in the late Ordovician mass-extinction event (~445 Myr ago). We characterised the habitability of the MW throughout its cosmic history as a function of galactocentric distance of TPs. We estimated the dangerous effects of transient astrophysical events (long and short GRBs and SNe) with a model that connects their rate to the specific star formation and metallicity evolution within the Galaxy throughout its cosmic history. Our model also accounts for the probability that TPs form around FGK and M stars. Until about six billion years ago, the outskirts of the Galaxy were the safest places to live, despite the relatively low density of TPs. In the last about four billion years, regions between 2 and 8 kpc from the center, which had a higher density of TPs, became the best places for a relatively safer biotic life growth. We confirm the hypothesis that one long GRB played a leading role in the late Ordovician mass-extinction event. In the last 500 Myr, the safest neighborhood in the Galaxy was a region at a distance of 2 to 8 kpc from the Galactic center, whereas the MW outskirts were sterilized by two to five long GRBs. △ Less

Submitted 26 January, 2021; v1 submitted 28 September, 2020; originally announced September 2020.

Comments: 11 pages, 10 figures, 4 tables. Accepted for pubblication in Astronomy & Astrophysics

Journal ref: A&A 647, A41 (2021)

Extended X-ray emission from the z=4.26 radio galaxy 4C 63.20

Authors: Kate Napier, Adi Foord, Elena Gallo, Gabriele Ghisellini, Edmund Hodges-Kluck, Jianfeng Wu, Francesco Haardt, Benedetta Ciardi

Abstract: We report on deep Chandra X-ray Telescope imaging observations of 4C 63.20, one of the few known radio galaxies at z>3.5. The X-ray counterpart is resolved into a core plus two off-nuclear sources that (combined) account for close to 30% of the total X-ray flux. Their morphology and orientation are consistent with a diffuse, lobe-like nature, albeit compact hotspots cannot be ruled out. The broadb… ▽ More We report on deep Chandra X-ray Telescope imaging observations of 4C 63.20, one of the few known radio galaxies at z>3.5. The X-ray counterpart is resolved into a core plus two off-nuclear sources that (combined) account for close to 30% of the total X-ray flux. Their morphology and orientation are consistent with a diffuse, lobe-like nature, albeit compact hotspots cannot be ruled out. The broadband spectral energy distribution of 4C 63.20 can be reproduced with a jet model where the majority of the radio flux can be ascribed to synchrotron emission from the hotspots, whereas the (non-nuclear) X-ray emission is produced via Inverse Compton (IC) off of Cosmic Microwave Background (CMB) photons within the extended lobes. This scenario is broadly consistent with the expectation from highly magnetized lobes in a hotter CMB, and supports the view that IC/CMB may quench less extreme radio lobes at high redshifts. △ Less

Submitted 20 July, 2020; originally announced July 2020.

Comments: 10 pages, 8 figures, accepted to MNRAS

Lyman Continuum Escape Fraction from Low-mass Starbursts at z=1.3

Authors: Anahita Alavi, James Colbert, Harry I. Teplitz, Brian Siana, Claudia Scarlata, Michael Rutkowski, Vihang Mehta, Alaina Henry, Y. Sophia Dai, Francesco Haardt, Micaela Bagley

Abstract: We present a new constraint on the Lyman Continuum (LyC) escape fraction at z~1.3. We obtain deep, high sensitivity far-UV imaging with the Advanced Camera for Surveys (ACS) Solar Blind Channel (SBC) on the Hubble Space Telescope (HST), targeting 11 star-forming galaxies at 1.2<z<1.4. The galaxies are selected from the 3D-HST survey to have high H$α$ equivalent width (EW) with EW > 190 Å, low stel… ▽ More We present a new constraint on the Lyman Continuum (LyC) escape fraction at z~1.3. We obtain deep, high sensitivity far-UV imaging with the Advanced Camera for Surveys (ACS) Solar Blind Channel (SBC) on the Hubble Space Telescope (HST), targeting 11 star-forming galaxies at 1.2<z<1.4. The galaxies are selected from the 3D-HST survey to have high H$α$ equivalent width (EW) with EW > 190 Å, low stellar mass (M* < 10^10 M_sun) and U-band magnitude of U<24.2. These criteria identify young, low metallicity star bursting populations similar to the primordial star-forming galaxies believed to have reionized the universe. We do not detect any LyC signal (with S/N >3) in the individual galaxies or in the stack in the far-UV images. We place $3σ$ limits on the relative escape fraction of individual galaxies to be f_{esc,rel}<[0.10-0.22] and a stacked $3σ$ limit of f_{esc,rel}<0.07. Comparing to the confirmed LyC emitters from the literature, the galaxies in our sample span similar ranges of various galaxy properties including stellar mass, dust attenuation, and star formation rate (SFR). In particular, we compare the distribution of H$α$ and [OIII] EWs of confirmed LyC emitters and non-detections including the galaxies in this study. Finally, we discuss if a dichotomy seen in the distribution of H$α$ EWs can perhaps distinguish the LyC emitters from the non-detections. △ Less

Submitted 10 July, 2020; originally announced July 2020.

Comments: 16 pages, 10 Figures, 2 Tables

Binary black holes in young star clusters: the impact of metallicity

Authors: Ugo N. Di Carlo, Michela Mapelli, Nicola Giacobbo, Mario Spera, Yann Bouffanais, Sara Rastello, Filippo Santoliquido, Mario Pasquato, Alessandro Ballone, Alessandro A. Trani, Stefano Torniamenti, Francesco Haardt

Abstract: Young star clusters are the most common birth-place of massive stars and are dynamically active environments. Here, we study the formation of black holes (BHs) and binary black holes (BBHs) in young star clusters, by means of 6000 N-body simulations coupled with binary population synthesis. We probe three different stellar metallicities (Z=0.02, 0.002 and 0.0002) and two initial density regimes (d… ▽ More Young star clusters are the most common birth-place of massive stars and are dynamically active environments. Here, we study the formation of black holes (BHs) and binary black holes (BBHs) in young star clusters, by means of 6000 N-body simulations coupled with binary population synthesis. We probe three different stellar metallicities (Z=0.02, 0.002 and 0.0002) and two initial density regimes (density at the half-mass radius $ρ_{\rm h}\ge{}3.4\times10^4$ and $\ge{1.5\times10^2}$ M$_\odot$ pc$^{-3}$ in dense and loose star clusters, respectively). Metal-poor clusters tend to form more massive BHs than metal-rich ones. We find $\sim{}6$, $\sim{}2$, and $<1$% of BHs with mass $m_{\rm BH}>60$ M$_\odot$ at Z=0.0002, 0.002 and 0.02, respectively. In metal-poor clusters, we form intermediate-mass BHs with mass up to $\sim{}320$ M$_\odot$. BBH mergers born via dynamical exchanges (exchanged BBHs) can be more massive than BBH mergers formed from binary evolution: the former (latter) reach total mass up to $\sim{}140$ M$_\odot$ ($\sim{}80$ M$_\odot$). The most massive BBH merger in our simulations has primary mass $\sim{}88$ M$_\odot$, inside the pair-instability mass gap, and a mass ratio of $\sim{}0.5$. Only BBHs born in young star clusters from metal-poor progenitors can match the masses of GW170729, the most massive event in O1 and O2, and those of GW190412, the first unequal-mass merger. We estimate a local BBH merger rate density $\sim{}110$ and $\sim{}55$ Gpc$^{-3}$ yr$^{-1}$, if we assume that all stars form in loose and dense star clusters, respectively. △ Less

Submitted 26 July, 2020; v1 submitted 20 April, 2020; originally announced April 2020.

Comments: 14 pages, 10 figures, 5 tables, accepted for publication in MNRAS

On the eccentricity evolution of massive black hole binaries in stellar backgrounds

Authors: Matteo Bonetti, Alexander Rasskazov, Alberto Sesana, Massimo Dotti, Francesco Haardt, Nathan W. C. Leigh, Manuel Arca Sedda, Giacomo Fragione, Elena Rossi

Abstract: We study the dynamical evolution of eccentric massive black hole binaries (MBHBs) interacting with unbound stars by means of an extensive set of three body scattering experiments. Compared to previous studies, we extend the investigation down to a MBHB mass ratio of $q=m_2/m_1=10^{-4}$, where $m_1$ and $m_2$ are the masses of the primary and secondary hole respectively. Contrary to a simple extrap… ▽ More We study the dynamical evolution of eccentric massive black hole binaries (MBHBs) interacting with unbound stars by means of an extensive set of three body scattering experiments. Compared to previous studies, we extend the investigation down to a MBHB mass ratio of $q=m_2/m_1=10^{-4}$, where $m_1$ and $m_2$ are the masses of the primary and secondary hole respectively. Contrary to a simple extrapolation from higher mass ratios, we find that for $q\lesssim 10^{-3}$ the eccentricity growth rate becomes negative, i.e., the binary {\it circularises} as it shrinks. This behaviour is due to the subset of interacting stars captured in metastable counter-rotating orbits; those stars tend to extract angular momentum from the binary, promoting eccentricity growth for $q>10^{-3}$, but tend to inject angular momentum into the binary driving it towards circularisation for $q<10^{-3}$. The physical origin of this behaviour requires a detailed study of the orbits of this subset of stars and is currently under investigation. Our findings might have important consequences for intermediate MBHs (IMBHs) inspiralling onto MBHs (e.g. a putative $10^3\rm M_{\odot}$ black hole inspiralling onto SgrA$^*$). △ Less

Submitted 4 February, 2020; v1 submitted 7 January, 2020; originally announced January 2020.

Comments: 6 pages, 5 figures; accepted for publication in MNRAS

Binary black holes in the pair-instability mass gap

Authors: Ugo N. Di Carlo, Michela Mapelli, Yann Bouffanais, Nicola Giacobbo, Filippo Santoliquido, Sandro Bressan, Mario Spera, Francesco Haardt

Abstract: Pair instability (PI) and pulsational PI prevent the formation of black holes (BHs) with mass $\gtrsim{}60$ M$_\odot$ from single star evolution. Here, we investigate the possibility that BHs with mass in the PI gap form via stellar mergers and multiple stellar mergers, facilitated by dynamical encounters in young star clusters. We analyze $10^4$ simulations, run with the direct N-body code nbody6… ▽ More Pair instability (PI) and pulsational PI prevent the formation of black holes (BHs) with mass $\gtrsim{}60$ M$_\odot$ from single star evolution. Here, we investigate the possibility that BHs with mass in the PI gap form via stellar mergers and multiple stellar mergers, facilitated by dynamical encounters in young star clusters. We analyze $10^4$ simulations, run with the direct N-body code nbody6++gpu coupled with the population synthesis code MOBSE. We find that up to $\sim{}6$~% of all simulated BHs have mass in the PI gap, depending on progenitor's metallicity. This formation channel is strongly suppressed in metal-rich (Z = 0.02) star clusters, because of stellar winds. BHs with mass in the PI gap are initially single BHs but can efficiently acquire companions through dynamical exchanges. We find that $\sim{}$21%, 10% and 0.5% of all binary BHs have at least one component in the PI mass gap at metallicity Z = 0.0002, 0.002 and 0.02, respectively. Based on the evolution of the cosmic star formation rate and metallicity, and under the assumption that all stars form in young star clusters, we predict that $\sim{}5$~% of all binary BH mergers detectable by advanced LIGO and Virgo at their design sensitivity have at least one component in the PI mass gap. △ Less

Submitted 6 July, 2020; v1 submitted 4 November, 2019; originally announced November 2019.

Comments: 7 pages, 3 figures, 2 tables, MNRAS, in press

The MUSE Ultra Deep Field (MUDF). II. Survey design and the gaseous properties of galaxy groups at 0.5 < z < 1.5

Authors: M. Fossati, M. Fumagalli, E. K. Lofthouse, V. D'Odorico, E. Lusso, S. Cantalupo, R. J. Cooke, S. Cristiani, F. Haardt, S. L. Morris, C. Peroux, L. J. Prichard, M. Rafelski, I. Smail, T. Theuns

Abstract: We present the goals, design, and first results of the MUSE Ultra Deep Field (MUDF) survey, a large programme using the Multi Unit Spectroscopic Explorer (MUSE) instrument at the ESO Very Large Telescope. The MUDF survey is collecting ~ 150 hours on-source of integral field optical spectroscopy in a 1.5 x 1.2 square arcmin region which hosts several astrophysical structures along the line of sight… ▽ More We present the goals, design, and first results of the MUSE Ultra Deep Field (MUDF) survey, a large programme using the Multi Unit Spectroscopic Explorer (MUSE) instrument at the ESO Very Large Telescope. The MUDF survey is collecting ~ 150 hours on-source of integral field optical spectroscopy in a 1.5 x 1.2 square arcmin region which hosts several astrophysical structures along the line of sight, including two bright z ~ 3.2 quasars with close separation (~ 500 kpc). Following the description of the data reduction procedures, we present the analysis of the galaxy environment and gaseous properties of seven groups detected at redshifts 0.5 < z < 1.5, spanning a large dynamic range in halo mass, log(Mh/Msun) ~ 11 - 13.5. For four of the groups, we find associated MgII absorbers tracing cool gas in high-resolution spectroscopy of the two quasars, including one case of correlated absorption in both sightlines at distance ~ 480 kpc. The absorption strength associated with the groups is higher than what has been reported for more isolated galaxies of comparable mass and impact parameters. We do not find evidence for widespread cool gas giving rise to strong absorption within these groups. Combining these results with the distribution of neutral and ionised gas seen in emission in lower-redshift groups, we conclude that gravitational interactions in the group environment strip gas from the galaxy haloes into the intragroup medium, boosting the cross section of cool gas and leading to the high fraction of strong MgII absorbers that we detect. △ Less

Submitted 21 September, 2019; v1 submitted 10 September, 2019; originally announced September 2019.

Comments: 20 pages, 12 figures, 4 tables, Accepted for publication in MNRAS

The high-energy radiation environment of the habitable-zone super-Earth LHS 1140b

Authors: R. Spinelli, F. Borsa, G. Ghirlanda, G. Ghisellini, S. Campana, F. Haardt, E. Poretti

Abstract: In the last few years many exoplanets in the habitable zone (HZ) of M-dwarfs have been discovered, but the X-ray/UV activity of cool stars is very different from that of our Sun. The high-energy radiation environment influences the habitability, plays a crucial role for abiogenesis, and impacts planetary atmospheres. LHS 1140b is a super-Earth-size planet orbiting in the HZ of LHS 1140, an M4.5 dw… ▽ More In the last few years many exoplanets in the habitable zone (HZ) of M-dwarfs have been discovered, but the X-ray/UV activity of cool stars is very different from that of our Sun. The high-energy radiation environment influences the habitability, plays a crucial role for abiogenesis, and impacts planetary atmospheres. LHS 1140b is a super-Earth-size planet orbiting in the HZ of LHS 1140, an M4.5 dwarf at ~15 parsecs. We present the results of a Swift X-ray/UV observing campaign. We characterize for the first time the X-ray/UV radiation environment of LHS 1140b. We measure the variability of the near ultraviolet (NUV) flux and estimate the far ultraviolet (FUV) flux with a correlation between FUV and NUV flux of a sample of low-mass stars in the GALEX archive. We highlight the presence of a dominating X-ray source close to the J2000 coordinates of LHS 1140, characterize its spectrum, and derive an X-ray flux upper limit for LHS 1140. We find that this contaminant source could have influenced the previously estimated spectral energy distribution. No significant variation of the NUV flux of LHS 1140 is found over 3 months, and we do not observe any flare during the 38 ks on the target. LHS 1140 is in the 25th percentile of least variable M4-M5 dwarfs of the GALEX sample. Analyzing the UV flux experienced by the HZ planet LHS 1140b, we find that outside the atmosphere it receives a NUV flux <2% with respect to that of the present-day Earth, while the FUV/NUV ratio is ~100-200 times higher. This represents a lower limit to the true FUV/NUV ratio since the GALEX FUV band does not include Lyman-alpha, which dominates the FUV output of low-mass stars. This is a warning for future searches for biomarkers, which must take into account this high ratio. The relatively low level and stability of UV flux experienced by LHS 1140b should be favorable for its present-day habitability. △ Less

Submitted 20 June, 2019; originally announced June 2019.

Comments: 6 pages, 4 figures, 2 tables. Accepted for pubblication in Astronomy & Astrophysics

Journal ref: A&A 627, A144 (2019)

The MUSE Ultra Deep Field (MUDF). I. Discovery of a group of Ly$α$ nebulae associated with a bright $z\approx 3.23$ quasar pair

Authors: E. Lusso, M. Fumagalli, M. Fossati, R. Mackenzie, R. M. Bielby, F. Arrigoni Battaia, S. Cantalupo, R. Cooke, S. Cristiani, P. Dayal, V. D'Odorico, F. Haardt, E. Lofthouse, S. Morris, C. Peroux, L. Prichard, M. Rafelski, R. Simcoe, A. M. Swinbank, T. Theuns

Abstract: We present first results from Multi Unit Spectroscopic Explorer (MUSE) observations at the Very Large Telescope in the MUSE Ultra Deep Field (MUDF), a $\approx 1.2\times 1.4$ arcmin$^2$ region for which we are collecting $\approx$200 hours of integral field spectroscopy. The $\approx 40$-hour observation completed to date reveals the presence of a group of three Ly$α$ nebulae associated with a bri… ▽ More We present first results from Multi Unit Spectroscopic Explorer (MUSE) observations at the Very Large Telescope in the MUSE Ultra Deep Field (MUDF), a $\approx 1.2\times 1.4$ arcmin$^2$ region for which we are collecting $\approx$200 hours of integral field spectroscopy. The $\approx 40$-hour observation completed to date reveals the presence of a group of three Ly$α$ nebulae associated with a bright quasar pair at $z\simeq3.23$ with projected separation of $\approx 500\rm~kpc$. Two of the nebulae are physically associated with the quasars which are likely powering the Ly$α$ emission, and extend for $\gtrsim 100~\rm kpc$ at a surface brightness level of $\approx 6\times 10^{-19}~\rm erg~s^{-1}~cm^{-2}~arcsec^{-2}$. A third smaller ($\approx$35 kpc) nebula lies at a velocity offset of $\approx 1550$ km s$^{-1}$. Despite their clustered nature, the two large nebulae have properties similar to those observed in isolated quasars and exhibit no sharp decline in flux at the current depth, suggesting an even more extended distribution of gas around the quasars. We interpret the shape and the alignment of the two brighter nebulae as suggestive of the presence of an extended structure connecting the two quasar host galaxies, as seen for massive galaxies forming within gas-rich filaments in cosmological simulations. △ Less

Submitted 12 March, 2019; v1 submitted 1 March, 2019; originally announced March 2019.

Comments: MNRAS Letters in press, 3 figures, revised version: minor changes following proofs

Merging black holes in young star clusters

Authors: Ugo N. Di Carlo, Nicola Giacobbo, Michela Mapelli, Mario Pasquato, Mario Spera, Long Wang, Francesco Haardt

Abstract: Searching for distinctive signatures, which characterize different formation channels of binary black holes (BBHs), is a crucial step towards the interpretation of current and future gravitational wave detections. Here, we investigate the demography of merging BBHs in young star clusters (SCs), which are the nursery of massive stars. We performed $4\times{} 10^3$ N-body simulations of SCs with met… ▽ More Searching for distinctive signatures, which characterize different formation channels of binary black holes (BBHs), is a crucial step towards the interpretation of current and future gravitational wave detections. Here, we investigate the demography of merging BBHs in young star clusters (SCs), which are the nursery of massive stars. We performed $4\times{} 10^3$ N-body simulations of SCs with metallicity $Z=0.002$, initial binary fraction $0.4$ and fractal initial conditions, to mimic the clumpiness of star forming regions. Our simulations include a novel population-synthesis approach based on the code MOBSE. We find that SC dynamics does not affect the merger rate significantly, but leaves a strong fingerprint on the properties of merging BBHs. More than 50 % of merging BBHs in young SCs form by dynamical exchanges in the first few Myr. Dynamically formed merging BBHs are significantly heavier than merging BBHs in isolated binaries: merging BBHs with total mass up to $\sim{}120$ M$_\odot$ form in young SCs, while the maximum total mass of merging BBHs in isolated binaries with the same metallicity is only $\sim{}70$ M$_\odot$. Merging BBHs born via dynamical exchanges tend to have smaller mass ratios than BBHs in isolated binaries. Furthermore, SC dynamics speeds up the merger: the delay time between star formation and coalescence is significantly shorter in young SCs. In our simulations, massive systems such as GW170729 form only via dynamical exchanges. Finally $\sim{}2$ % of merging BBHs in young SCs have mass in the pair-instability mass gap ($\sim{}60-120$ M$_\odot$). This represents a unique fingerprint of merging BBHs in SCs. △ Less

Submitted 3 January, 2019; originally announced January 2019.

Comments: 15 pages, 11 figures, 5 tables

Post-Newtonian Evolution of Massive Black Hole Triplets in Galactic Nuclei: IV. Implications for LISA

Authors: Matteo Bonetti, Alberto Sesana, Francesco Haardt, Enrico Barausse, Monica Colpi

Abstract: Coalescing massive black hole binaries (MBHBs) of $10^{4-7} \rm M_{\odot}$, forming in the aftermath of galaxy mergers, are primary targets of the space mission LISA, the {\it Laser Interferometer Space Antenna}. An assessment of LISA detection prospects requires an estimate of the abundance and properties of MBHBs that form and evolve during the assembly of cosmic structures. To this aim, we empl… ▽ More Coalescing massive black hole binaries (MBHBs) of $10^{4-7} \rm M_{\odot}$, forming in the aftermath of galaxy mergers, are primary targets of the space mission LISA, the {\it Laser Interferometer Space Antenna}. An assessment of LISA detection prospects requires an estimate of the abundance and properties of MBHBs that form and evolve during the assembly of cosmic structures. To this aim, we employ a semi-analytic model to follow the co-evolution of MBHBs within their host galaxies. We identify three major evolutionary channels driving the binaries to coalescence: two standard paths along which the binary evolution is driven by interactions with the stellar and/or gaseous environment, and a novel channel where MBHB coalescence occurs during the interaction with a third black hole. For each channel, we follow the orbital evolution of MBHBs with physically motivated models that include a self-consistent treatment of the orbital eccentricity. We find that LISA will detect between $\approx 25$ and $\approx 75$ events per year depending on the seed model. We show that triple-induced coalescences can range from a few detected events up to $\sim 30\%$ of the total detected mergers. Moreover, even if the standard gas/stars-driven evolutionary channels should fail and MBHBs were to stall, triple interactions would still occur as a result of the hierarchical nature of galaxy formation, resulting in about $\approx 10$ to $\approx 20$ LISA detections per year. Remarkably, triple interactions among the black holes can produce coalescing binaries with large eccentricities ($\gtrsim 0.9$) upon entrance into the LISA band. This eccentricity will remain significant ($\sim 0.1$) also at merger, requiring suitable templates for parameter estimation. △ Less

Submitted 18 May, 2019; v1 submitted 3 December, 2018; originally announced December 2018.

Comments: 17 pages, 10 figures, Accepted MNRAS

Barred galaxies in cosmological zoom-in simulations: the importance of feedback

Authors: Tommaso Zana, Pedro R. Capelo, Massimo Dotti, Lucio Mayer, Alessandro Lupi, Francesco Haardt, Silvia Bonoli, Sijing Shen

Abstract: Bars are a key factor in the long-term evolution of spiral galaxies, in their unique role in redistributing angular momentum and transporting gas and stars on large scales. The Eris-suite simulations are cosmological zoom-in, N-body, smoothed-particle hydrodynamic simulations built to follow the formation and evolution of a Milky Way-sized galaxy across the build-up of the large scale structure. H… ▽ More Bars are a key factor in the long-term evolution of spiral galaxies, in their unique role in redistributing angular momentum and transporting gas and stars on large scales. The Eris-suite simulations are cosmological zoom-in, N-body, smoothed-particle hydrodynamic simulations built to follow the formation and evolution of a Milky Way-sized galaxy across the build-up of the large scale structure. Here we analyse and describe the outcome of two particular simulations taken from the Eris suite - ErisBH and Eris2k - which mainly differ in the prescriptions employed for gas cooling, star formation, and feedback from supernovae and black holes. Our study shows that the enhanced effective feedback in Eris2k, due to the collective effect of the different micro-physics implementations, results in a galaxy which is less massive than its ErisBH counterpart till z~1. However, when the stellar content is large enough so that global dynamical instabilities can be triggered, the galaxy in Eris2k develops a stronger and more extended bar with respect to ErisBH. We demonstrate that he structural properties and time evolution of the two bars are very different. Our results highlight the importance of accurate sub-grid prescriptions in cosmological zoom-in simulations of the process of galaxy formation and evolution, and the possible use of a statistical sample of barred galaxies to assess the strength of the stellar feedback. △ Less

Submitted 2 July, 2019; v1 submitted 17 October, 2018; originally announced October 2018.

Comments: 16 pages, 14 Figure; accepted for publication in MNRAS; added references and figures, more details, results unchanged

MCMC determination of the cosmic UV background at $z\simeq 0$ from H$α$ fluorescence

Authors: Davide Caruso, Francesco Haardt, Michele Fumagalli, Sebastiano Cantalupo

Abstract: In a recent paper (Fumagalli et al. 2017) we reported on the detection of a diffuse H$α$ glow in the outskirts of the nearby, edge-on disc galaxy UGC 7321 observed with the Multi Unit Spectroscopic Explorer (MUSE) at the ESO Very Large Telescope. By interpreting the H$α$ emission as fluorescence arising from hydrogen ionised by an external (i.e., extragalactic) radiation field, we estimated the UV… ▽ More In a recent paper (Fumagalli et al. 2017) we reported on the detection of a diffuse H$α$ glow in the outskirts of the nearby, edge-on disc galaxy UGC 7321 observed with the Multi Unit Spectroscopic Explorer (MUSE) at the ESO Very Large Telescope. By interpreting the H$α$ emission as fluorescence arising from hydrogen ionised by an external (i.e., extragalactic) radiation field, we estimated the UV background (UVB) intensity in terms of HI ionisation rate (per ion) at $z\simeq 0$ to be in the range $Γ_{\rm HI}\sim 6-8\times 10^{-14}$ s$^{-1}$. In the present work, by performing radiative transfer calculations over a large set of models of the gaseous disc of UGC 7321, we refine our estimate and through an MCMC analysis derive a value for the photoionisaton rate of $Γ_{\rm HI}=7.27^{+2.93}_{-2.90}\times 10^{-14}$ s$^{-1}$. In particular, our analysis demonstrates that this value is robust against large variations in the galaxy model and that the uncertainties are mainly driven by the errors associated with the observed H$α$ surface brightness. Our measurement is consistent with several recent determinations of the same quantity by a completely independent technique (i.e., flux decrement analysis of the Ly$α$-forest), and support the notion that the low redshift UVB is largely dominated by active galactic nuclei (AGNs), possibly with no need of further contribution from star forming galaxies. △ Less

Submitted 29 August, 2018; originally announced August 2018.

Comments: 6 pages, 2 figures, submitted to MNRAS

Black holes, gravitational waves and fundamental physics: a roadmap

Authors: Leor Barack, Vitor Cardoso, Samaya Nissanke, Thomas P. Sotiriou, Abbas Askar, Krzysztof Belczynski, Gianfranco Bertone, Edi Bon, Diego Blas, Richard Brito, Tomasz Bulik, Clare Burrage, Christian T. Byrnes, Chiara Caprini, Masha Chernyakova, Piotr Chrusciel, Monica Colpi, Valeria Ferrari, Daniele Gaggero, Jonathan Gair, Juan Garcia-Bellido, S. F. Hassan, Lavinia Heisenberg, Martin Hendry, Ik Siong Heng , et al. (181 additional authors not shown)

Abstract: The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horiz… ▽ More The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress. △ Less

Submitted 1 February, 2019; v1 submitted 13 June, 2018; originally announced June 2018.

Comments: White Paper for the COST action "Gravitational Waves, Black Holes, and Fundamental Physics", 272 pages, 12 figures; v4: updated references and author list. Overall improvements and corrections. To appear in Classical and Quantum Gravity

Bar resilience to flybys in a cosmological framework

Authors: Tommaso Zana, Massimo Dotti, Pedro R. Capelo, Lucio Mayer, Francesco Haardt, Sijing Shen, Silvia Bonoli

Abstract: It has been proposed that close interactions with satellite galaxies can significantly perturb the morphology of the main galaxy. However, the dynamics of an already formed bar following the interaction with the external environment has not been studied in detail in a fully cosmological context. In this work, analysing the cosmological zoom-in simulation Eris2k, we study the effects that a very un… ▽ More It has been proposed that close interactions with satellite galaxies can significantly perturb the morphology of the main galaxy. However, the dynamics of an already formed bar following the interaction with the external environment has not been studied in detail in a fully cosmological context. In this work, analysing the cosmological zoom-in simulation Eris2k, we study the effects that a very unequal-mass flyby crossing the stellar disc has on the stability of the pre-existing bar. We characterize the evolution of the bar strength and length showing that the perturbation exerted by the flyby shuffles the orbits of stars for less than one Gyr. After this time, the bar shows a remarkable resilience, reforming with properties comparable to those it had before the interaction. Our work shows that close unequal-mass encounters, the most frequent interactions occurring during the evolution of cosmic structures, have (i) an overall minor impact on the global evolution of the bar in the long term, still (ii) the effect is destructive and (iii) a very weak interaction is sufficient to dismantle a strong bar leading to its "apparent death". As a consequence, due to the non-negligible duration of the bar-less period, a fraction of observed spiral galaxies classified as non-barred could be prone to bar formation. △ Less

Submitted 20 July, 2018; v1 submitted 9 May, 2018; originally announced May 2018.

Comments: 6 pages, 3 Figure; Accepted for publication in MNRAS

Consistent modelling of the meta-galactic UV background and the thermal/ionization history of the intergalactic medium

Authors: Ewald Puchwein, Francesco Haardt, Martin G. Haehnelt, Piero Madau

Abstract: Recent observations suggest that hydrogen reionization ends late ($z \simeq 6$) and proceeds quickly. We present here a new model of the meta-galactic UV/X-ray background (UVB) that is consistent with this. It adopts the most recent determinations of the ionizing emissivity due to stars and AGN, as well as of the HI absorber column density distribution. Another major improvement is a new treatment… ▽ More Recent observations suggest that hydrogen reionization ends late ($z \simeq 6$) and proceeds quickly. We present here a new model of the meta-galactic UV/X-ray background (UVB) that is consistent with this. It adopts the most recent determinations of the ionizing emissivity due to stars and AGN, as well as of the HI absorber column density distribution. Another major improvement is a new treatment of the intergalactic medium (IGM) opacity for ionizing photons that is able to consistently capture the transition from a neutral to an ionized IGM. Previous synthesis models of the UVB, when used in simulations, yield reionization and thermal histories that are inconsistent with the assumed ionizing emissivities. With our new treatment, this discrepancy is fully resolved. In our fiducial model, galaxies leaking $\lesssim 18\%$ of their Lyman continuum emission drive HI reionization, while AGN drive HeII reionization (completing at $z \simeq 6.2$ and $2.8$ respectively). Due to the limited time available for cooling between HI and HeII reionization, higher IGM temperatures are predicted for late reionization scenarios. In our fiducial model, the predicted temperatures agree well with observational constraints at $z \lesssim 4$, while being slightly high compared to (somewhat uncertain) data above that. Models with a larger contribution of AGN are instead disfavoured by the temperature data, as well as by measurements of the HI and HeII Lyman-$α$ forest opacities. We also present "equivalent-equilibrium" ionization/heating rates that mimic our fiducial UVB model for use in simulation codes that assume ionization equilibrium. △ Less

Submitted 15 January, 2018; originally announced January 2018.

Comments: 23 pages, 11 figures, submitted to MNRAS, see ancillary files for machine-readable rate tables

Post-Newtonian evolution of massive black hole triplets in galactic nuclei -- III. A robust lower limit to the nHz stochastic background of gravitational waves

Authors: Matteo Bonetti, Alberto Sesana, Enrico Barausse, Francesco Haardt

Abstract: Inspiraling massive black-hole binaries (MBHBs) forming in the aftermath of galaxy mergers are expected to be the loudest gravitational-wave (GW) sources relevant for pulsar-timing arrays (PTAs) at nHz frequencies. The incoherent overlap of signals from a cosmic population of MBHBs gives rise to a stochastic GW background (GWB) with characteristic strain around $h_c\sim10^{-15}$ at a reference fre… ▽ More Inspiraling massive black-hole binaries (MBHBs) forming in the aftermath of galaxy mergers are expected to be the loudest gravitational-wave (GW) sources relevant for pulsar-timing arrays (PTAs) at nHz frequencies. The incoherent overlap of signals from a cosmic population of MBHBs gives rise to a stochastic GW background (GWB) with characteristic strain around $h_c\sim10^{-15}$ at a reference frequency of 1 yr$^{-1}$, although uncertainties around this value are large. Current PTAs are piercing into the GW amplitude range predicted by MBHB-population models, but no detection has been reported so far. To assess the future success prospects of PTA experiments, it is therefore important to estimate the minimum GWB level consistent with our current understanding of the formation and evolution of galaxies and massive black holes (MBHs). To this purpose, we couple a semianalytic model of galaxy evolution and an extensive study of the statistical outcome of triple MBH interactions. We show that even in the most pessimistic scenario where all MBHBs stall before entering the GW-dominated regime, triple interactions resulting from subsequent galaxy mergers inevitably drive a considerable fraction of the MBHB population to coalescence. At frequencies relevant for PTA, the resulting GWB is only a factor of 2-to-3 suppressed compared to a fiducial model where binaries are allowed to merge over Gyr timescales. Coupled with current estimates of the expected GWB amplitude range, our findings suggest that the minimum GWB from cosmic MBHBs is unlikely to be lower than $h_c\sim 10^{-16}$ (at $f = 1$ yr$^{-1}$), well within the expected sensitivity of projected PTAs based on future observations with FAST, MeerKAT and SKA. △ Less

Submitted 9 April, 2018; v1 submitted 18 September, 2017; originally announced September 2017.

Comments: 16 pages, 9 figures, accepted for publication in MNRAS

Journal ref: MNRAS, 2018, 477, 2599

Post-Newtonian evolution of massive black hole triplets in galactic nuclei -- II. Survey of the parameter space

Authors: Matteo Bonetti, Francesco Haardt, Alberto Sesana, Enrico Barausse

Abstract: Massive black hole binaries (MBHBs) are expected to form at the centre of merging galaxies during the hierarchical assembly of the cosmic structure, and are expected to be the loudest sources of gravitational waves (GWs) in the low frequency domain. However, because of the dearth of energy exchanges with background stars and gas, many of these MBHBs may stall at separations too large for GW emissi… ▽ More Massive black hole binaries (MBHBs) are expected to form at the centre of merging galaxies during the hierarchical assembly of the cosmic structure, and are expected to be the loudest sources of gravitational waves (GWs) in the low frequency domain. However, because of the dearth of energy exchanges with background stars and gas, many of these MBHBs may stall at separations too large for GW emission to drive them to coalescence in less than a Hubble time. Triple MBH systems are then bound to form after a further galaxy merger, triggering a complex and rich dynamics that can eventually lead to MBH coalescence. Here we report on the results of a large set of numerical simulations, where MBH triplets are set in spherical stellar potentials and MBH dynamics is followed through 2.5 post-Newtonian order in the equations of motion. From our full suite of simulated systems we find that a fraction $\simeq 20-30$% of the MBH binaries that would otherwise stall are led to coalesce within a Hubble time. The corresponding coalescence timescale peaks around 300 Myr, while the eccentricity close to the plunge, albeit small, is non-negligible ($\lesssim 0.1$). We construct and discuss marginalised probability distributions of the main parameters involved and, in a companion paper of the series, we will use the results presented here to forecast the contribution of MBH triplets to the GW signal in the nHz regime probed by Pulsar Timing Array experiments. △ Less

Submitted 9 April, 2018; v1 submitted 18 September, 2017; originally announced September 2017.

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

Journal ref: MNRAS, 2018, 477, 3910