Source author record

Francesco Haardt

Francesco Haardt appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

31works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

31 published item(s)

preprint2022arXiv

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

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.

preprint2022arXiv

Maximally informed Bayesian modelling of disc galaxies

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.

preprint2022arXiv

Misaligned Spinning Binary Black Hole Mergers in Hot Magnetized Plasma

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.

preprint2022arXiv

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

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.

preprint2022arXiv

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

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.

preprint2020arXiv

Binary black holes in the pair-instability mass gap

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.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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$^*$).

preprint2016arXiv

Post-Newtonian Evolution of Massive Black Hole Triplets in Galactic Nuclei: I. Numerical Implementation and Tests

Massive black-hole binaries (MBHBs) are thought to be the main source of gravitational waves (GWs) in the low-frequency domain surveyed by ongoing and forthcoming Pulsar Timing Array campaigns and future space-borne missions, such as {\it eLISA}. However, many low-redshift MBHBs in realistic astrophysical environments may not reach separations small enough to allow significant GW emission, but rather stall on (sub)pc-scale orbits. This "last-parsec problem" can be eased by the appearance of a third massive black hole (MBH) -- the "intruder" -- whose action can force, under certain conditions, the inner MBHB on a very eccentric orbit, hence allowing intense GW emission eventually leading to coalescence. A detailed assessment of the process, ultimately driven by the induced Kozai-Lidov oscillations of the MBHB orbit, requires a general relativistic treatment and the inclusion of external factors, such as the Newtonian precession of the intruder orbit in the galactic potential and its hardening by scattering off background stars. In order to tackle this problem, we developed a three-body Post-Newtonian (PN) code framed in a realistic galactic potential, including both non-dissipative 1PN and 2PN terms, and dissipative terms such as 2.5PN effects, orbital hardening of the outer binary, and the effect of the dynamical friction on the early stages of the intruder dynamics. In this first paper of a series devoted at studing the dynamics of MBH triplets from a cosmological perspective, we describe, test and validate our code.

preprint2016arXiv

The Early Growth of the First Black Holes

With detections of quasars powered by increasingly massive black holes (BHs) at increasingly early times in cosmic history over the past decade, there has been correspondingly rapid progress made on the theory of early BH formation and growth. Here we review the emerging picture of how the first massive BHs formed from the primordial gas and then grew to supermassive scales. We discuss the initial conditions for the formation of the progenitors of these seed BHs, the factors dictating the initial masses with which they form, and their initial stages of growth via accretion, which may occur at super-Eddington rates. Finally, we briefly discuss how these results connect to large-scale simulations of the growth of supermassive BHs over the course of the first billion years following the Big Bang.

preprint2015arXiv

Cosmic Reionization after Planck: Could Quasars Do It All?

We assess a model of late cosmic reionization in which the ionizing background radiation arises entirely from high redshift quasars and other active galactic nuclei (AGNs). The low optical depth to Thomson scattering reported by the Planck Collaboration pushes the redshift of instantaneous reionization down to z=8.8^{+1.7}_{-1.4} and greatly reduces the need for significant Lyman-continuum emission at very early times. We show that, if recent claims of a numerous population of faint AGNs at z=4-6 are upheld, and the high inferred AGN comoving emissivity at these epochs persists to higher, z~10, redshifts, then active galaxies may drive the reionization of hydrogen and helium with little contribution from normal star-forming galaxies. We discuss an AGN-dominated scenario that satisfies a number of observational constraints: the HI photoionization rate is relatively flat over the range 2<z<5, hydrogen gets fully reionized by z=5.7, and the integrated Thomson scattering optical depth is tau=0.056, in agreement with measurements based on the Lya opacity of the intergalactic medium (IGM) and cosmic microwave background (CMB) polarization. It is a prediction of the model that helium gets doubly reionized before redshift 4, the heat input from helium reionization dominates the thermal balance of the IGM after hydrogen reionization, and z>5 AGNs provide a significant fraction of the unresolved X-ray background at 2 keV. Singly- and doubly-ionized helium contribute about 13% to tau, and the HeIII volume fraction is already 50% when hydrogen becomes fully reionized.

preprint2015arXiv

High redshift AGNs and HI reionisation: limits from the unresolved X-ray background

The rapidly declining population of bright quasars at z~3 appears to make an increasingly small contribution to the ionising background at the HI Lyman limit. It is then generally though that massive stars in (pre-)galactic systems may provide the additional ionising flux needed to complete HI reionisation by z>6. A galaxy dominated background, however, may require that the escape fraction of Lyman continuum radiation from high redshift galaxies is as high as 10%, a value somewhat at odds with (admittedly scarce) observational constraints. High escape fractions from dwarf galaxies have been advocated, or, alternatively, a so-far undetected (or barely detected) population of unobscured, high-redshift faint AGNs. Here we question the latter hypothesis, and show that such sources, to be consistent with the measured level of the unresolved X-ray background at z=0, can provide a fraction of the HII filling factor not larger than 13% by z=6. The fraction rises to <27% in the somewhat extreme case of a constant comoving redshift evolution of the AGN emissivity. This still calls for a mean escape fraction of ionising photons from high-z galaxies >10%.

preprint2015arXiv

Massive black hole and gas dynamics in mergers of galaxy nuclei - II. Black hole sinking in star-forming nuclear discs

Mergers of gas-rich galaxies are key events in the hierarchical built-up of cosmic structures, and can lead to the formation of massive black hole binaries. By means of high-resolution hydrodynamical simulations we consider the late stages of a gas-rich major merger, detailing the dynamics of two circumnuclear discs, and of the hosted massive black holes during their pairing phase. During the merger gas clumps with masses of a fraction of the black hole mass form because of fragmentation. Such high-density gas is very effective in forming stars, and the most massive clumps can substantially perturb the black hole orbits. After $\sim 10$ Myr from the start of the merger a gravitationally bound black hole binary forms at a separation of a few parsecs, and soon after, the separation falls below our resolution limit of $0.39$ pc. At the time of binary formation the original discs are almost completely disrupted because of SNa feedback, while on pc scales the residual gas settles in a circumbinary disc with mass $\sim 10^5 M_\odot$. We also test that binary dynamics is robust against the details of the SNa feedback employed in the simulations, while gas dynamics is not. We finally highlight the importance of the SNa time-scale on our results.

preprint2015arXiv

The Lyman Continuum Escape Fraction of Low-Mass Star-Forming Galaxies at z~1

To date no direct detection of Lyman continuum emission has been measured for intermediate--redshift z~1 star-forming galaxies . We combine HST grism spectroscopy with GALEX UV and ground--based optical imaging to extend the search for escaping Lyman continuum to a large (~600) sample of z~1 low-mass, moderately star-forming galaxies selected initially on H$α$ emission. The characteristic escape fraction of LyC from SFGs that populate this parameter space remains weakly constrained by previous surveys, but these faint SFGs are assumed to play a significant role in the reionization of neutral hydrogen in the intergalactic medium (IGM) at high redshift (z>6). We do not make an unambiguous detection of escaping LyC radiation from this $z\sim1$ sample, individual non--detections to constrain the absolute Lyman continuum escape fraction, $f_{esc}$<2.1% (3$σ$). We measure upper limits of $f_{esc}$<9.6% from a sample of SFGs selected on high H$α$ equivalent width (EW>200Å), which are thought to be close analogs of high redshift sources of reionization. For reference, we also present an emissivity--weighted escape fraction which is useful as a measurement of the general contribution of the SFGs to the z~1 ionizing UV background. In the discussion, we consider the implications of these intermediate redshift constraints for the reionization of hydrogen in the intergalactic medium at high (z > 6) redshift. If we assume our $z\sim1$ SFGs, for which we measure this emissivity-weighted $f_{esc}$, are analogs to the high redshift sources of reionization, we find is difficult reconcile reionization by faint (M$_{UV}<-13$) SFGs with a low escape fraction ($f_{esc}<$3%), with constraints from independent high redshift observations. If $f_{esc}$ evolves with redshift, reionization by SFGs may be consistent with observations from Planck.

preprint2015arXiv

The photoheating of the intergalactic medium in synthesis models of the UV background

We compare cosmological hydrodynamical simulations combined with the homogeneous metagalactic UV background (UVB) of Haardt & Madau (2012) (HM2012) to observations of the Lyman-alpha forest that are sensitive to the thermal and ionization state of the intergalactic medium (IGM). The transition from optically thick to thin photoheating predicted by the simple one-zone, radiative transfer model implemented by HM2012 predicts a thermal history that is in remarkably good agreement with the observed rise of the IGM temperature at z~3 if we account for the expected evolution of the volume filling factor of HeIII. Our simulations indicate that there may be, however, some tension between the observed peak in the temperature evolution and the rather slow evolution of the HeII opacities suggested by recent Hubble Space Telescope/COS measurements. The HM2012 UVB also underpredicts the metagalactic hydrogen photoionization rate required by our simulations to match the observed opacity of the forest at z>4 and z<2.

preprint2014arXiv

Massive black hole and gas dynamics in galaxy nuclei mergers. I. Numerical implementation

Numerical effects are known to plague adaptive mesh refinement (AMR) codes when treating massive particles, e.g. representing massive black holes (MBHs). In an evolving background, they can experience strong, spurious perturbations and then follow unphysical orbits. We study by means of numerical simulations the dynamical evolution of a pair MBHs in the rapidly and violently evolving gaseous and stellar background that follows a galaxy major merger. We confirm that spurious numerical effects alter the MBH orbits in AMR simulations, and show that numerical issues are ultimately due to a drop in the spatial resolution during the simulation, drastically reducing the accuracy in the gravitational force computation. We therefore propose a new refinement criterion suited for massive particles, able to solve in a fast and precise way for their orbits in highly dynamical backgrounds. The new refinement criterion we designed enforces the region around each massive particle to remain at the maximum resolution allowed, independently upon the local gas density. Such maximally-resolved regions then follow the MBHs along their orbits, and effectively avoids all spurious effects caused by resolution changes. Our suite of high resolution, adaptive mesh-refinement hydrodynamic simulations, including different prescriptions for the sub-grid gas physics, shows that the new refinement implementation has the advantage of not altering the physical evolution of the MBHs, accounting for all the non trivial physical processes taking place in violent dynamical scenarios, such as the final stages of a galaxy major merger.

preprint2014arXiv

Super-Critical Growth of Massive Black Holes from Stellar-Mass Seeds

We consider super-critical accretion with angular momentum onto stellar-mass black holes as a possible mechanism for growing billion-solar-mass holes from light seeds at early times. We use the radiatively-inefficient "slim disk" solution -- advective, optically thick flows that generalize the standard geometrically thin disk model -- to show how mildly super-Eddington intermittent accretion may significantly ease the problem of assembling the first massive black holes when the Universe was less than 0.8 Gyr old. Because of the low radiative efficiencies of slim disks around non-rotating as well as rapidly rotating holes, the mass e-folding timescale in this regime is nearly independent of the spin parameter. The conditions that may lead to super-critical growth in the early Universe are briefly discussed.

preprint2014arXiv

The Photon Underproduction Crisis

We examine the statistics of the low-redshift Lyman-alpha forest from smoothed particle hydrodynamic simulations in light of recent improvements in the estimated evolution of the cosmic ultraviolet background (UVB) and recent observations from the Cosmic Origins Spectrograph (COS). We find that the value of the metagalactic photoionization rate required by our simulations to match the observed properties of the low-redshift Lyman-alpha forest is a factor of 5 larger than the value predicted by state-of-the art models for the evolution of this quantity. This mismatch results in the mean flux decrement of the Lyman-alpha forest being underpredicted by at least a factor of 2 (a 10-sigma discrepancy with observations) and a column density distribution of Lyman-alpha forest absorbers systematically and significantly elevated compared to observations over nearly two decades in column density. We examine potential resolutions to this mismatch and find that either conventional sources of ionizing photons (galaxies and quasars) must be significantly elevated relative to current observational estimates or our theoretical understanding of the low-redshift universe is in need of substantial revision.

preprint2013arXiv

Black-hole mass estimates for a homogeneous sample of bright flat-spectrum radio quasars

We have selected a complete sample of flat-spectrum radio quasars (FSRQs) from the WMAP 7-yr catalog within the SDSS area, all with measured redshift, and have compared the black hole mass estimates based on fitting a standard accretion disk model to the `blue bump' with those obtained from the commonly used single epoch virial method. The sample comprises 79 objects with a flux density limit of 1 Jy at 23 GHz, 54 of which (68%) have a clearly detected `blue bump'. Thirty-four of the latter have, in the literature, black hole mass estimates obtained with the virial method. The mass estimates obtained from the two methods are well correlated. If the calibration factor of the virial relation is set to $f=4.5$, well within the range of recent estimates, the mean logarithmic ratio of the two mass estimates is equal to zero with a dispersion close to the estimated uncertainty of the virial method. The fact that the two independent methods agree so closely in spite of the potentially large uncertainties associated with each lends strong support to both of them. The distribution of black-hole masses for the 54 FSRQs in our sample with a well detected blue bump has a median value of $7.4\times 10^{8}\,M_\odot$. It declines at the low mass end, consistent with other indications that radio loud AGNs are generally associated with the most massive black holes, although the decline may be, at least partly, due to the source selection. The distribution drops above $\log(M_\bullet/M_\odot) = 9.4$, implying that ultra-massive black holes associated with FSRQs must be rare.

preprint2012arXiv

Evolution of binary black holes in self gravitating discs: dissecting the torques

We study the interplay between gas accretion and gravity torques in changing a binary elements and its total angular momentum (L) budget. Especially, we analyse the physical origin of the gravity torques (T_g) and their location within the disc. We analyse 3D SPH simulations of the evolution of initially quasi-circular massive black hole binaries (BHBs) residing in the central hollow of massive self-gravitating circumbinary discs. We use different thermodynamics within the cavity and for the numerical size of the black holes to show that (i) the BHB eccentricity growth found previously is a general result, independent of the accretion and the adopted thermodynamics; (ii) the semi-major axis decay depends both on the T_g and on the interplay with the disc-binary L-transfer due to accretion; (iii) the spectral structure of the T_g is predominately caused by disc edge overdensities and spiral arms developing in the body of the disc and, in general, does not reflect directly the period of the binary; (iv) the net T_g changes sign across the BHB corotation radius. We quantify the relative importance of the two, which appear to depend on the thermodynamical properties of the instreaming gas, and which is crucial in assessing the disc-binary L-transfer; (v) the net torque manifests as a purely kinematic (non-resonant) effect as it stems from the cavity, where the material flows in and out in highly eccentric orbits. Both accretion onto the black holes and the interaction with gas streams inside the cavity must be taken into account to assess the fate of the BHB. Moreover, the total torque exerted by the disc affects L(BHB) by changing all the elements (mass, mass ratio, eccentricity, semimajor axis) of the BHB. Common prescriptions equating tidal torque to semi-major axis shrinking might therefore be poor approximations for real astrophysical systems.

preprint2011arXiv

A new model for the extragalactic gamma-ray background

We present a two-parameter model of the extragalactic gamma-ray background (EGB) in the 0.1-100 GeV range as measured by the Large Area Telescope (LAT) onboard the Fermi satellite. The EGB can be fully explained as the sum of three distinct components, namely blazars, non-beamed AGNs (Seyfert galaxies and QSOs), and cosmic rays from star-forming galaxies. The contribution to the background from beamed sources is obtained by fitting the Fermi-LAT blazar differential number counts assuming that the gamma-ray luminosity function is directly proportional to the radio luminosity function of FRI and FRII galaxies. The high energy emission from non-beamed AGNs is instead determined by popular synthesis models of the observed X-ray background. Finally, the EGB is fit by adding a third component arising from pion decay in cosmic rays, assuming that such component is closely linked to the cosmic star formation history. We find that blazars dominate at energies \geq 10 GeV, for E \leq 0.2 GeV the main contribution is from non-beamed AGNs, while cosmic rays are required in between. Because of absorption due to interaction of gamma-rays with the extra-galactic background light, our model falls short at the highest energies probed by LAT, (\geq 70 GeV), leaving room to a possible contribution from dark matter particle annihilation. As an example, a particle of mass ~ 0.5 TeV and cross section <σv>~ 5*10^{-26} cm^3 s^-1 can accomodate the data.

preprint2011arXiv

A new model for the Fermi-LAT Extragalactic Gamma-Ray Background

We compute the contribution of blazars (Flat Spectrum Radio Quasars and BL-Lacs) to the Fermi-LAT Extragalactic gamma-ray Background, assuming a gamma-ray luminosity function proportional to the luminosity function of radio galaxies and the Spectral Energy Distribution of blazars given by the so-called blazar sequence. We find that blazars account for the 45% of the total Fermi-LAT EGB. However, at low (<10 GeV) and high (>50 GeV) energies the blazar contribution falls short to account for the observed background light. At E<10 GeV, the Extragalactic gamma-ray Background can be fully explained by considering the additional contribution of star-forming galaxies, whereas at very high energies, where the absorption due to Extragalactic Background Light dominates, the observed background calls for an additional unknown component.

preprint2011arXiv

Blazars in the early Universe

We investigate the relative occurrence of radio--loud and radio-quiet quasars in the first billion years of the Universe, powered by black holes heavier than one billion solar masses. We consider the sample of high-redshfit blazars detected in the hard X-ray band in the 3-years all sky survey performed by the Burst Alert Telescope (BAT) onboard the Swift satellite. All the black holes powering these blazars exceed a billion solar mass, with accretion luminosities close to the Eddington limit. For each blazar pointing at us, there must be hundreds of similar sources (having black holes of similar masses) pointing elsewhere. This puts constraints on the density of billion solar masses black holes at high redshift (z>4), and on the relative importance of (jetted) radio-loud vs radio-quiet sources. We compare the expected number of high redshift radio--loud sources with the high luminosity radio-loud quasars detected in the Sloan Digital Sky Survey (SDSS), finding agreement up to z~3, but a serious deficit at z>3 of SDSS radio-loud quasars with respect to the expectations. We suggest that the most likely explanations for this disagreement are: i) the ratio of blazar to misaligned radio-sources decreases by an order of magnitude above z=3, possibly as a result of a decrease of the average bulk Lorentz factor; ii) the SDSS misses a large fraction of radio-loud sources at high redshifts, iii) the SDSS misses both radio-loud and radio-quiet quasars at high redshift, possibly because of obscuration or because of collimation of the optical-UV continuum in systems accreting near Eddington. These explanations imply very different number density of heavy black holes at high redshifts, that we discuss in the framework of the current ideas about the relations of dark matter haloes at high redshifts and the black hole they host.

preprint2011arXiv

New Population Synthesis Models of the Cosmic UV Background

We present improved synthesis models of the evolving spectrum of the UV/X-ray diffuse background, updating and extending our previous results. A "minimal cosmic reionization model" is presented in which the galaxy UV emissivity traces recent determinations of the cosmic history of star formation, the luminosity-weighted escape fraction of hydrogen-ionizing radiation increases rapidly with redshift, from 1% at z=3 to 50% at z=8.5, the clumping factor of the intergalactic medium is C=3, and Population III stars make a negligible contribution to the metagalactic flux. The model provides a good fit to the hydrogen-ionization rates inferred from flux decrement measurements, and yields an optical depth to Thomson scattering, tau_e=0.085, that is in agreement with WMAP results.

preprint2011arXiv

Radiative transfer in a clumpy universe: IV. New synthesis models of the cosmic UV/X-ray background

We present improved synthesis models of the evolving spectrum of the UV/X-ray diffuse background, updating and extending our previous results. Five new main components are added to our radiative transfer code CUBA: (1) the sawtooth modulation of the background intensity from resonant line absorption in the Lyman series of cosmic hydrogen and helium; (2) the X-ray emission from obscured and unobscured quasars; (3) a piecewise parameterization of the distribution in redshift and column density of intergalactic absorbers that fits recent measurements of the mean free path of 1 ryd photons; (4) an accurate treatment of the photoionization structure of absorbers; and (5) the UV emission from star-forming galaxies at all redshifts. We provide tables of the predicted HI and HeII photoionization and photoheating rates for use, e.g., in cosmological hydrodynamics simulations of the Lya forest, and a new metallicity-dependent calibration to the UV luminosity density-star formation rate density relation. A "minimal cosmic reionization model" is also presented in which the galaxy UV emissivity traces recent determinations of the cosmic history of star formation, the luminosity-weighted escape fraction of hydrogen-ionizing radiation increases rapidly with lookback time, the clumping factor of the high-redshift intergalactic medium evolves following the results of recent hydrodynamic simulations, and Population III stars and miniquasars make a negligible contribution to the metagalactic flux. The model provides a good fit to the hydrogen-ionization rates inferred from flux decrement and proximity effect measurements, predicts that cosmological HII (HeIII) regions overlap at redshift 6.7 (2.8), and yields an optical depth to Thomson scattering that is in agreement with WMAP results. (Abridged)

preprint2009arXiv

GeV Gamma-Ray Attenuation and the High-Redshift UV Background

We present new calculations of the evolving UV background out to the epoch of cosmological reionization and make predictions for the amount of GeV gamma-ray attenuation by electron-positron pair production. Our results are based on recent semi-analytic models of galaxy formation, which provide predictions of the dust-extinguished UV radiation field due to starlight, and empirical estimates of the contribution due to quasars. We account for the reprocessing of ionizing photons by the intergalactic medium. We test whether our models can reproduce estimates of the ionizing background at high redshift from flux decrement analysis and proximity effect measurements from quasar spectra, and identify a range of models that can satisfy these constraints. Pair-production against soft diffuse photons leads to a spectral cutoff feature for gamma rays observed between 10 and 100 GeV. This cutoff varies with redshift and the assumed star formation and quasar evolution models. We find only negligible amounts of absorption for gamma rays observed below 10 GeV for any emission redshift. With observations of high-redshift sources in sufficient numbers by the Fermi Gamma-ray Space Telescope and new ground-based instruments it should be possible to constrain the extragalactic background light in the UV and optical portion of the spectrum.

preprint2009arXiv

He II absorption and the sawtooth spectrum of the cosmic far-UV background

Cosmic ultraviolet background radiation between 3 and 4 Ryd is reprocessed by resonant line absorption in the Lyman series of intergalactic He II. This process results in a sawtooth modulation of the radiation spectrum from the He II Lya frequency to the Lyman limit. The size of this modulation is a sensitive probe of the epoch of helium reionization and of the sources that keep the intergalactic medium (IGM) highly ionized. For large absorption opacities, the background intensity will peak at frequencies just above each resonance, go to zero at resonance, and fluctuate greatly just below resonance. The He II sawtooth modulation may be one of the missing ingredients needed in the modelling of the abundances of metal ions like C III and Si IV observed in the IGM at redshift 3.

preprint1997arXiv

Simultaneous X-ray and gamma-ray observations of Cyg X-1 in the hard state by Ginga and OSSE

We present four X-ray/gamma-ray spectra of Cyg X-1 observed in the hard ('low') state simultaneously by Ginga and GRO/OSSE on 1991 July 6. The 3-30 keV Ginga spectra are well represented by power laws with an energy spectral index of alpha~0.6 and a Compton reflection component including a fluorescent Fe K-alpha corresponding to the solid angle of the reflector of ~0.3 times 2 Pi. The overall Ginga/OSSE spectra can be modelled by repeated Compton scattering in a mildly-relativistic, tau ~1, plasma. However, the high-energy cutoff is steeper than that due to single-temperature thermal Comptonisation. It can be described by a superposition of dominant optically-thin, thermal emission at kT~140 keV and a Wien-like component from an optically-thick plasma at kT~50 keV. The X-ray spectra do not show the presence of an anisotropy break required if thermal Compton scattering takes place in a corona above a cold disc. Also, the flat spectral index shows that the plasma is soft-photon starved, i.e., the luminosity in incident soft X-ray seed photons is very much less than that in the hard X-rays. Furthermore, the observed solid angle of the reflector is significantly less than 2 pi. These facts taken together strongly rule out a disc-corona geometry. Rather, the observed spectra are consistent with a geometry in which the cold accretion disc (which both supplies the seed soft X-rays and reflects hard X-rays) only exists at large radii, while the Comptonising hot plasma is located in an inner region with no cold disc. This hot plasma consists of either pure electron-positron pairs if the source size is ~5 Schwarzschild radii or it contains also protons if the size is larger.

preprint1996arXiv

The Synchrotron Boiler: a Thermalizer in Seyfert Galaxies

There are difficulties in understanding what keeps the plasma thermalized in compact sources, especially during rapid variations of the emitted flux. Particle-particle collisions are too inefficient in hot rarefied plasmas, and a faster process is called for. Synchrotron absorption is such a process. We show that relativistic electrons can thermalize in a few synchrotron cooling times by emitting and absorbing cyclo-synchrotron photons. The resulting equilibrium distribution is a Maxwellian at low energies, with a high energy power law tail when Compton cooling is important. Assuming that the particles emit completely self absorbed synchrotron radiation while they at the same time Compton scatter ambient UV photons, we calculate the time dependent behavior of the distribution function, and the final high energy spectra.

preprint1994arXiv

The Contribution of the Obscuring Torus to the X--Ray Spectrum of Seyfert Galaxies: A Test for the Unification Model

The presence of an obscuring torus around the nucleus of Seyfert galaxies, as supposed in the popular unification scheme, can strongly modify the X--ray spectrum for both type 1 and 2 Seyfert galaxies. We find that if the torus is Compton thick, it can scatter at small angles a significant fraction of the nuclear radiation, and contribute to the continuum of Seyfert 1 galaxies above $\sim$10 keV, and to the fluorescence iron line at 6.4 keV. At large inclination angles and for large torus column densities, the spectrum is attenuated by photoabsorption and Compton scattering, while the iron fluorescence line produced by the torus can have large equivalent widths. Even after dilution by the continuum scattered by the warm material outside the torus, this iron line could be strong enough to explain the ``cold" component in the spectrum of NGC 1068. We stress that the complex pattern of the predicted variability can be a powerful tool for constraining the parameters of the model, such as the column density of the torus, its inclination and the amount of warm scattering material.