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$\alpha$ equivalent width (EW) with EW > 190 \AA, 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\sigma$ limits on the relative escape fraction of individual galaxies to be f_{esc,rel}<[0.10-0.22] and a stacked $3\sigma$ 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$\alpha$ 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$\alpha$ EWs can perhaps distinguish the LyC emitters from the non-detections.

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.

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 $\rho_{\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.

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

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.

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.

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.

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.

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$\alpha$ 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$\alpha$ 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.

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.

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.

In a recent paper (Fumagalli et al. 2017) we reported on the detection of a diffuse H$\alpha$ 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$\alpha$ 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 $\Gamma_{\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 $\Gamma_{\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$\alpha$ 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$\alpha$-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.

The emergence of a large-scale stellar bar is one of the most striking features in disc galaxies. By means of state-of-the-art cosmological zoom-in simulations, we study the formation and evolution of bars in Milky Way-like galaxies in a fully cosmological context, including the physics of gas dissipation, star formation, and supernova feedback. Our goal is to characterise the actual trigger of the non-axisymmetric perturbation that leads to the strong bar observable in the simulations at z=0, discriminating between an internal/secular versus an external/tidal origin. To this aim, we run a suite of cosmological zoom-in simulations altering the original history of galaxy-satellite interactions at a time when the main galaxy, though already bar-unstable, does not feature any non-axisymmetric structure yet. We find that the main effect of a late minor merger and of a close fly-by is to delay the time of bar formation and those two dynamical events are not directly responsible for the development of the bar and do not alter significantly its global properties (e.g. its final extension). We conclude that, once the disc has grown to a mass large enough to sustain global non-axisymmetric modes, then bar formation is inevitable.

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.

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.

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.

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-$\alpha$ 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.

We highlight some subtleties that affect naive implementations of quadrupolar and octupolar gravitational waveforms from numerically-integrated trajectories of three-body systems. Some of those subtleties arise from the requirement that the source be contained in its "coordinate near zone" when applying the standard PN formulae for gravitational-wave emission, and from the need to use the non-linear Einstein equations to correctly derive the quadrupole emission formula. We show that some of these subtleties were occasionally overlooked in the literature, with consequences for published results. We also provide prescriptions that lead to correct and robust predictions for the waveforms computed from numerically-integrated orbits.

Recent work suggests that strong emission line, star-forming galaxies may be significant Lyman Continuum leakers. We combine archival HST broadband ultraviolet and optical imaging (F275W and F606W, respectively) with emission line catalogs derived from WFC3 IR G141 grism spectroscopy to search for escaping Lyman Continuum (LyC) emission from homogeneously selected $z\simeq$2.5 SFGs. We detect no escaping Lyman Continuum from SFGs selected on [OII] nebular emission (N=208) and, within a narrow redshift range, on [OIII]/[OII]. We measure 1$\sigma$ upper limits to the LyC escape fraction relative to the non-ionizing UV continuum from [OII] emitters, $f_{esc}<$5.6%, and strong [OIII]/[OII]$>$5 ELGs, $f_{esc}<$14.0%. Our observations are not deep enough to detect $f_{esc}\lesssim$10% typical of the low redshift Lyman continuum emitters. However, we find that this population represents a small fraction of the star-forming galaxy population at $z\simeq$2. Thus, unless the number of extreme emission line galaxies grows substantially to z$>$6, such galaxies may be insufficient for reionization. Deeper survey data in the rest-frame ionizing UV will be necessary to determine whether strong line ratios could be useful for pre-selecting LyC leakers at high redshift.

We present results from the "Ponos" simulation suite on the early evolution of a massive, $M_{\rm vir}(z=0)=1.2\times 10^{13}$ M$_{\odot}$ galaxy. At $z\gtrsim6$, before feedback from a central supermassive black hole becomes dominant, the main galaxy has a stellar mass $\sim 2\times 10^{9}$ M$_{\odot}$ and a star formation rate $\sim 20$ M$_{\odot}$ yr$^{-1}$. The galaxy sits near the expected main sequence of star-forming galaxies at those redshifts, and resembles moderately star-forming systems observed at $z>5$. The high specific star formation rate results in vigorous heating and stirring of the gas by supernovae feedback, and the galaxy develops a thick and turbulent disc, with gas velocity dispersion $\sim 40$ km s$^{-1}$, rotation to dispersion ratio $\sim 2$, and with a significant amount of gas at $\sim 10^5$ K. The Toomre parameter always exceeds the critical value for gravito-turbulence, $Q\sim 1.5-2$, mainly due to the contribution of warm/hot gas inside the disc. Without feedback, a nearly gravito-turbulent regime establishes with similar gas velocity dispersion and lower $Q$. We propose that the "hot and turbulent" disc regime seen in our simulations, unlike the "cold and turbulent" gravito-turbulent regime of massive clumpy disc galaxies at $z\sim 1-2$, is a fundamental characterisation of main sequence galaxies at $z\gtrsim 6$, as they can sustain star formation rates comparable to those of low-mass starbursts at $z=0$. This results in no sustained coherent gas inflows through the disc, and in fluctuating and anisotropic mass transport, possibly postponing the assembly of the bulge and causing the initial feeding of the central black hole to be highly intermittent.

We report the detection of extended Halpha emission from the tip of the HI disk of the nearby edge-on galaxy UGC 7321, observed with the Multi Unit Spectroscopic Explorer (MUSE) instrument at the Very Large Telescope. The Halpha surface brightness fades rapidly where the HI column density drops below N(HI) = 10^19 cm^-2 , consistent with fluorescence arising at the ionisation front from gas that is photoionized by the extragalactic ultraviolet background (UVB). The surface brightness measured at this location is (1.2 +/- 0.5)x10^-19 erg/s/cm^2/arcsec^2, where the error is mostly systematic and results from the proximity of the signal to the edge of the MUSE field of view, and from the presence of a sky line next to the redshifted Halpha wavelength. By combining the Halpha and the HI 21 cm maps with a radiative transfer calculation of an exponential disk illuminated by the UVB, we derive a value for the HI photoionization rate of Gamma ~ (6-8)x10^-14 1/s . This value is consistent with transmission statistics of the Lyalpha forest and with recent models of a UVB which is dominated by quasars.

In an effort to understand the cause of the apparent depletion in the number density of radio-loud AGNs at $z>3$, this work investigates the viability of the so-called Cosmic Microwave Background (CMB) quenching mechanism of intrinsically jetted, high-z AGNs, whereby Inverse Compton scattering of CMB photons off electrons within the extended lobes results in a substantial dimming of the lobe synchrotron emission at GHz frequencies, while simultaneously boosting their diffuse X-ray signal. We focus on five $z>3.5$ radio galaxies that have sufficiently deep Chandra exposure (> 50 ks) to warrant a meaningful investigation of any extended X-ray emission. For those objects with evidence for statistically significant extended X-ray lobes (4C 41.17 and 4C 03.24), we combine the Chandra measurements with literature data at lower frequencies to assemble the systems' Spectral Energy Distributions (SEDs), and utilize state-of-the-art SED modelling (Ghisellini et al. 2015) -- including emission from the disk, torus, jet, hotspots, and lobes -- to infer their physical parameters. For both radio galaxies, the magnetic energy density in the hotspots is found to exceed the energy density in CMB photons, wheres the opposite is true for the lobes. This implies that any extended synchrotron emission likely originates from the hotspots themselves, rather than the lobes. Conversely, Inverse Compton scattering of CMB photons dominates the extended X-ray emission from the lobes, which are effectively "radio-quenched". As a result, CMB quenching is effective in these systems in spite of the fact that the observed X-ray to radio luminosity ratio does not bear the signature $(1+z)^4$ dependence of the CMB energy density.

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.

We present new Chandra X-ray observations of seven low-mass black holes (~1e6 Msun) accreting at low Eddington ratios between -2.0<log L/Ledd<-1.5. We compare the X-ray properties of these seven low-mass active galactic nuclei (AGN) to a total of 73 other low-mass AGN in the literature with published Chandra observations (with Eddington ratios extending from -2.0<log L/Ledd<-0.1). We do not find any statistical differences between low- and high-Eddington ratio low-mass AGN in the distributions of their X-ray to ultraviolet luminosity ratios (aox), or in their X-ray spectral shapes. Furthermore, the aox distribution of low-L/Ledd AGN displays an X-ray weak tail that is also observed within high-L/Ledd objects. Our results indicate that between -2<log L/Ledd<-0.1, there is no systematic change in the structure of the accretion flow for active galaxies hosting 1e6 Msun black holes. We examine the accuracy of current bolometric luminosity estimates for our low-L/Ledd objects with new Chandra observations, and it is plausible that their Eddington ratios could be underestimated by up to an order of magnitude. If so, then in analogy with weak emission line quasars, we suggest that accretion from a geometrically thick, radiatively inefficient `slim disk' could explain their diverse properties in aox. Alternatively, if current Eddington ratios are in fact correct (or overestimated), then the X-ray weak tail would imply that there is diversity in disk/corona couplings among individual low-mass objects. Finally, we conclude by noting that the aox distribution for low-mass black holes may have favorable consequences for the epoch of cosmic reionization being driven by AGN.

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.

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$\alpha$ 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$\sigma$). We measure upper limits of $f_{esc}$<9.6% from a sample of SFGs selected on high H$\alpha$ equivalent width (EW>200\AA), 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.

The rapid assembly of the massive black holes that power the luminous quasars observed at $z \sim 6-7$ remains a puzzle. Various direct collapse models have been proposed to head-start black hole growth from initial seeds with masses $\sim 10^5\,\rm M_\odot$, which can then reach a billion solar mass while accreting at the Eddington limit. Here we propose an alternative scenario based on radiatively inefficient super-critical accretion of stellar-mass holes embedded in the gaseous circum-nuclear discs (CNDs) expected to exist in the cores of high redshift galaxies. Our sub-pc resolution hydrodynamical simulations show that stellar-mass holes orbiting within the central 100 pc of the CND bind to very high density gas clumps that arise from the fragmentation of the surrounding gas. Owing to the large reservoir of dense cold gas available, a stellar-mass black hole allowed to grow at super-Eddington rates according to the "slim disc" solution can increase its mass by 3 orders of magnitudes within a few million years. These findings are supported by simulations run with two different hydro codes, RAMSES based on the Adaptive Mesh Refinement technique and GIZMO based on a new Lagrangian Godunov-type method, and with similar, but not identical, sub-grid recipes for star formation, supernova feedback, black hole accretion and feedback. The low radiative efficiency of super-critical accretion flows are instrumental to the rapid mass growth of our black holes, as they imply modest radiative heating of the surrounding nuclear environment.

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.

Stellar bars are a common feature in massive disc galaxies. On a theoretical ground, the response of gas to a bar is generally thought to cause nuclear starbursts and, possibly, AGN activity once the perturbed gas reaches the central super-massive black hole. By means of high resolution numerical simulations we detail the purely dynamical effects that a forming bar exerts on the gas of an isolated disc galaxy. The galaxy is initially unstable to the formation of non-axisymmetric structures, and within 1 Gyr it develops spiral arms that eventually evolve into a central stellar bar on kpc scale. A first major episode of gas inflow occurs during the formation of the spiral arms while at later times, when the stellar bar is establishing, a low density region is carved between the bar co-rotational and inner Lindblad resonance radii. The development of such "dead zone" inhibits further massive gas inflows. Indeed, the gas inflow reaches its maximum during the relatively fast bar formation phase and not, as often assumed, when the bar is fully formed. We conclude that the low efficiency of long-lived, evolved bars in driving gas toward galactic nuclei is the reason why observational studies have failed to establish an indisputable link between bars and AGNs. On the other hand, the high efficiency in driving strong gas inflows of the intrinsically transient process of bar formation suggests that the importance of bars as drivers of AGN activity in disc galaxies has been overlooked so far. We finally prove that our conclusions are robust against different numerical implementations of the hydrodynamics routinely used in galaxy evolution studies.

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.

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.

The very existence of more than a dozen of high-redshift (z>4) blazars indicates that a much larger population of misaligned powerful jetted AGN was already in place when the Universe was <1.5 Gyr old. Such parent population proved to be very elusive, and escaped direct detection in radio surveys so far. High redshift blazars themselves seem to be failing in producing extended radio-lobes, raising questions about the connection between such class and the vaster population of radio-galaxies. We show that the interaction of the jet electrons with the intense cosmic microwave background (CMB) radiation explains the lack of extended radio emission in high redshift blazars and in their parent population, helping to explain the apparently missing misaligned counterparts of high redshift blazars. On the other hand, the emission from the more compact and more magnetised hot spots are less affected by the enhanced CMB energy density. By modelling the spectral energy distribution of blazar lobes and hot spots we find that most of them should be detectable by low frequency deep radio observations, e.g., by LOw-Frequency ARray for radio astronomy (LOFAR) and by relatively deep X-ray observations with good angular resolution, e.g., by the Chandra satellite. At high redshifts, the emission of a misaligned relativistic jet, being de-beamed, is missed by current large sky area surveys. The isotropic flux produced in the hot spots can be below ~1 mJy and the isotropic lobe radio emission is quenched by the CMB cooling. Consequently, even sources with very powerful jets can go undetected in current radio surveys, and misclassified as radio-quiet AGNs.

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%.

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.

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.

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.

We discuss how the interaction between the electrons in a relativistic jet and the Cosmic Microwave Background (CMB) affects the observable properties of radio-loud AGN at early epochs. At high z the magnetic energy density in the radio lobes of powerful radio-loud quasars can be exceeded by the energy density of the CMB (because of its (1+z)^4 dependance). In this case, relativistic electrons cool preferentially by scattering off CMB photons, rather than by synchrotron. Thus, sources sharing the same intrinsic properties have different extended radio and X-ray luminosities when located at different z: more distant sources are less luminous in radio and more luminous in X-rays than their closer counterparts. Instead, in compact regions where the local magnetic field still exceeds the CMB in terms of energy density, synchrotron radiation would be unaffected by the presence of the CMB. Such regions include the compact inner jet and the so-called hot spots in the radio lobes. The decrease in radio luminosity is larger in misaligned sources, whose radio flux is dominated by the extended isotropic component. These sources can fail detection in current flux limited radio surveys, and therefore they are possibly under-represented in the associated samples. As the cooling time is longer for lower energy electrons, the radio luminosity deficit due to the CMB photons is less important at low radio frequencies. Therefore objects not detected so far in current surveys at a few GHz could be picked up by low frequency deep surveys, such as LOFAR and SKA. Until then, we can estimate the number of high redshift radio-loud AGNs through the census of their aligned proxies, i.e., blazars. Indeed, their observed radio emission arises in the inner and strongly magnetized compact core of the relativistic jet, and not affected by inverse Compton scattering off CMB photons.

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.

It has been recently suggested that supermassive black holes at z = 5-6 might form from super-fast (\dot M > 10^4 Msun/yr) accretion occurring in unstable, massive nuclear gas disks produced by mergers of Milky-Way size galaxies. Interestingly, such mechanism is claimed to work also for gas enriched to solar metallicity. These results are based on an idealized polytropic equation of state assumption, essentially preventing the gas from cooling. We show that under more realistic conditions, the disk rapidly (< 1 yr) cools, the accretion rate drops, and the central core can grow only to \approx 100 Msun. In addition, most of the disk becomes gravitationally unstable in about 100 yr, further quenching the accretion. We conclude that this scenario encounters a number of difficulties that possibly make it untenable.

In powerful radio-quiet active galactic nuclei (AGN), black holes heavier than one billion solar masses form at a redshift ~1.5-2. Supermassive black holes in jetted radio-loud AGN seems to form earlier, at a redshift close to 4. The ratio of active radio-loud to radio-quiet AGN hosting heavy black holes is therefore a rather a strong function of redshift. We report on some recent evidence supporting this conclusion, gathered from the Burst Alert Telescope (BAT, onboard Swift) and by the Large Area Telescope (LAT, onboard Fermi). We suggest that the more frequent occurrence of relativistic jets in the most massive black holes at high redshifts, compared to later times, could be due to the average black hole spin being greater in the distant past, or else to the jet helping a fast accretion rate (or some combination of the two scenarios). We emphasize that the large total accretion efficiency of rapidly spinning black holes inhibits a fast growth, unless a large fraction of the available gravitational energy of the accreted mass is not converted into radiation, but used to form and maintain a powerful jet.

We place firm upper limits on the global accretion history of massive black holes at z>5 from the recently measured unresolved fraction of the cosmic X-ray background. The maximum allowed unresolved intensity observed at 1.5 keV implies a maximum accreted-mass density onto massive black holes of rho_acc < 1.4E4 M_sun Mpc^{-3} for z>5. Considering the contribution of lower-z AGNs, the value reduces to rho_acc < 0.66E4 M_sun Mpc^{-3}. The tension between the need for the efficient and rapid accretion required by the observation of massive black holes already in place at z>7 and the strict upper limit on the accreted mass derived from the X-ray background may indicate that black holes are rare in high redshift galaxies, or that accretion is efficient only for black holes hosted by rare galaxies.

We study the spectral properties of the unresolved cosmic X-ray background (CXRB) in the 1.5-7.0 keV energy band with the aim of providing an observational constraint on the statistical properties of those sources that are too faint to be individually probed. We made use of the Swift X-ray observation of the Chandra Deep Field South complemented by the Chandra data. Exploiting the lowest instrument background (Swift) together with the deepest observation ever performed (Chandra) we measured the unresolved emission at the deepest level and with the best accuracy available today. We find that the unresolved CXRB emission can be modeled by a single power law with a very hard photon index Gamma=0.1+/-0.7 and a flux of 5(+/-3)E-12 cgs in the 2.0-10 keV energy band (1 sigma error). Thanks to the low instrument background of the Swift-XRT, we significantly improved the accuracy with respect to previous measurements. These results point towards a novel ingredient in AGN population synthesis models, namely a positive evolution of the Compton-thick AGN population from local Universe to high redshift.

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.

Massive black hole binaries at sub-parsec separations may display in their spectra anomalously small flux ratios between the MgII and CIV broad emission lines, i.e. F_MgII/F_CIV <~ 0.1, due to the erosion of the broad line region around the active, secondary black hole, by the tidal field of the primary. In Paper I by Montuori et al. (2011), we focussed on broad lines emitted by gas bound to the lighter accreting member of a binary when the binary is at the center of a hollow density region (the gap) inside a circum-binary disc. The main aim of this new study is at exploring the potential contribution to the broad line emission by the circum-binary disc and by gaseous streams flowing toward the black hole through the gap. We carry out a post-process analysis of data extracted from a SPH simulation of a circum-binary disc around a black hole binary. Our main result is that the MgII to CIV flux ratio can be reduced to ~ 0.1 within an interval of sub-pc binary separations of the order of a ~ (0.01-0.2)(f_Edd/0.1)^(1/2) pc corresponding to orbital periods of ~ (20-200) (f_Edd/0.1)^(3/4) years for a secondary BH mass in the range M_2 ~ 10^7-10^9 M_sun and a binary mass ratio of 0.3. At even closer separations this ratio returns to increase to values that are indistinguishable from the case of a single AGN (typically F_MgII/F_CIV ~ 0.3-0.4) because of the contribution to the MgII line from gas in the circum-binary disc.

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)

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.

Narrow Line Seyfert 1 galaxies (NLS1) are generally believed to host "under-massive" black holes with respect to the predictions from the host galaxy -- black hole mass scale relations. Black hole masses in NLS1 are estimated from the continuum luminosity and the width of broad emission lines. Here we show that the "mass deficit" can be canceled out if we assume that the broad line region (BLR) in type-1 AGN has a flat geometry, which is seen face-on in NLS1. The detection of relativistic jets aligned along the line of sight in a number of NLS1 supports this picture. Moreover, a flat geometry of the BLR is also suggested by a general trend of the mass deficit as a function of the line width observed in other type-1 AGN, from quasars to BL Lac objects, and is consistent with a simple extension of the Unified Model of AGN to the BLR geometry.

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.

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 <{\sigma}v>~ 5*10^{-26} cm^3 s^-1 can accomodate the data.

ORIGIN is a proposal for the M3 mission call of ESA aimed at the study of metal creation from the epoch of cosmic dawn. Using high-spectral resolution in the soft X-ray band, ORIGIN will be able to identify the physical conditions of all abundant elements between C and Ni to red-shifts of z=10, and beyond. The mission will answer questions such as: When were the first metals created? How does the cosmic metal content evolve? Where do most of the metals reside in the Universe? What is the role of metals in structure formation and evolution? To reach out to the early Universe ORIGIN will use Gamma-Ray Bursts (GRBs) to study their local environments in their host galaxies. This requires the capability to slew the satellite in less than a minute to the GRB location. By studying the chemical composition and properties of clusters of galaxies we can extend the range of exploration to lower redshifts (z ~ 0.2). For this task we need a high-resolution spectral imaging instrument with a large field of view. Using the same instrument, we can also study the so far only partially detected baryons in the Warm-Hot Intergalactic Medium (WHIM). The less dense part of the WHIM will be studied using absorption lines at low redshift in the spectra for GRBs.

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.

We investigate on the spectral properties of an active black hole, member of a massive (10^7 - 10^9 Msun) sub-parsec black hole binary. We work under the hypothesis that the binary, surrounded by a circum-binary disc, has cleared a gap, and that accretion occurs onto the secondary black hole fed by material closer to the inner edge of the disc. Broad line emission clouds orbit around the active black hole and suffer erosion due to tidal truncation at the Roche Lobe surface, following gap opening and orbital decay. We consider three of the most prominent broad emission lines observed in the spectra of AGNs, i.e. CIV, MgII and H{\beta}, and compute the flux ratios between the lines of MgII and CIV (FMgII/FCIV) and those of MgII and H{\beta} (FMgII/FH{\beta}). We find that close black hole binaries have FMgII/FCIV up to one order of magnitude smaller than single black holes. By contrast FMgII/FH{\beta} may be significantly reduced only at the shortest separations. Peculiarly low values of line flux ratios together with large velocity offsets between the broad and narrow emission lines and/or periodic variability in the continuum (on timescales >= years) would identify genuine sub-pc binary candidates.

With the release of the first year Fermi catalogue, the number of blazars detected above 100 MeV lying at high redshift has been largely increased. There are 28 blazars at z>2 in the "clean" sample. All of them are Flat Spectrum Radio Quasars (FSRQs). We study and model their overall spectral energy distribution in order to find the physical parameters of the jet emitting region, and for all of them we estimate their black hole masses and accretion rates. We then compare the jet with the accretion disk properties, setting these sources in the broader context of all the other bright gamma-ray or hard X-ray blazars. We confirm that the jet power correlates with the accretion luminosity. We find that the high energy emission peak shifts to smaller frequencies as the observed luminosity increases, according to the blazar sequence, making the hard X-ray band the most suitable for searching the most luminous and distant blazars.

We investigate the physical properties of the 10 blazars at redshift greater than 2 detected in the 3-years all sky survey performed by the Burst Alert Telescope (BAT) onboard the Swift satellite. We find that the jets of these blazars are among the most powerful known. Furthermore, the mass of their central black hole, inferred from the optical-UV bump, exceeds a few billions of solar masses, with accretion luminosities being a large fraction of the Eddington one. We compare their properties with those of the brightest blazars of the 3-months survey performed by the Large Area Telescope (LAT) onboard the Fermi satellite. We find that the BAT blazars have more powerful jets, more luminous accretion disks and larger black hole masses than LAT blazars. These findings can be simply understood on the basis of the blazar sequence, that suggests that the most powerful blazars have a spectral energy distribution with a high energy peak at MeV (or even sub-MeV) energies. This implies that the most extreme blazars can be found more efficiently in hard X-rays, rather than in the high energy gamma-ray band. We then discuss the implications of our findings for future missions, such as the New Hard X-ray Mission (NHXM) and especially the Energetic X-ray Imaging Survey Telescope (EXIST) mission which, during its planned 2 years all sky survey, is expected to detect thousands of blazars, with a few of them at z greater than 6.

Galactic nuclei host central massive objects either in the form of supermassive black holes or nuclear stellar clusters. Recent investigations have shown that both components co-exist in at least a few galaxies. In this paper we explore the possibility of a connection between nuclear star clusters and black holes that establishes at the moment of their formation. We here model the evolution of high redshift discs, hosted in dark matter halos with virial temperatures 10^4 K, whose gas has been polluted with metals just above the critical metallicity for fragmentation. A nuclear cluster forms as a result of a central starburst from gas inflowing from the unstable disc. The nuclear stellar cluster provides a suitable environment for the formation of a black hole seed, ensuing from runaway collisions among the most massive stars. Typical masses for the nuclear stellar clusters at the time of black hole formation (z~10) are inthe range 10^4-10^6 solar masses and have half mass radii < 0.5 pc. The black holes forming in these dense, high redshift clusters can have masses in the range ~300-2000 solar masses.

Using high resolution hydrodynamical simulations, we explore the spin evolution of massive dual black holes orbiting inside a circumnuclear disc, relic of a gas-rich galaxy merger. The black holes spiral inwards from initially eccentric co or counter-rotating coplanar orbits relative to the disc's rotation, and accrete gas that is carrying a net angular momentum. As the black hole mass grows, its spin changes in strength and direction due to its gravito-magnetic coupling with the small-scale accretion disc. We find that the black hole spins loose memory of their initial orientation, as accretion torques suffice to align the spins with the angular momentum of their orbit on a short timescale (<1-2 Myr). A residual off-set in the spin direction relative to the orbital angular momentum remains, at the level of <10 degrees for the case of a cold disc, and <30 degrees for a warmer disc. Alignment in a cooler disc is more effective due to the higher coherence of the accretion flow near each black hole that reflects the large-scale coherence of the disc's rotation. If the massive black holes coalesce preserving the spin directions set after formation of a Keplerian binary, the relic black hole resulting from their coalescence receives a relatively small gravitational recoil. The distribution of recoil velocities inferred from a simulated sample of massive black hole binaries has median <70 km/s much smaller than the median resulting from an isotropic distribution of spins.

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.

A strong hard X-ray luminosity from a blazar flags the presence of a very powerful jet. If the jet power is in turn related to the mass accretion rate, the most luminous hard X-ray blazars should pinpoint the largest accretion rates, and therefore the largest black hole masses. These ideas are confirmed by the Swift satellite observations of the blazar S5 0014+813, at the redshift z=3.366. Swift detected this source with all its three instruments, from the optical to the hard X-rays. Through the construction of its spectral energy distribution we are confident that its optical-UV emission is thermal in origin. Associating it to the emission of a standard optically thick geometrically thin accretion disk, we find a black hole mass of 40 billion solar masses, radiating at 40% the Eddington value. The derived mass is among the largest ever found. Super-Eddington slim disks or thick disks with the presence of a collimating funnel can in principle reduce the black hole mass estimate, but tends to produce spectra bluer than observed.

In this Letter we explore the hypothesis that the quasar SDSSJ092712.65+294344.0 is hosting a massive black hole binary embedded in a circumbinary disc. The lightest, secondary black hole is active, and gas orbiting around it is responsible for the blue-shifted broad emission lines with velocity off-set of 2650 km/s, relative to the galaxy rest frame. As the tidal interaction of the binary with the outer disc is expected to excavate a gap, the blue-shifted narrow emission lines are consistent with being emitted from the low-density inhomogeneous gas of the hollow region. From the observations we infer a binary mass ratio q ~ 0.3, a mass for the primary of M1 ~ 2 billion Msun and a semi-major axis of 0.34 pc, corresponding to an orbital period of 370 years. We use the results of cosmological merger trees to estimate the likely-hood of observing SDSSJ092712.65+294344.0 as recoiling black hole or as a binary. We find that the binary hypothesis is preferred being one hundred times more probable than the ejection hypothesis. If SDSSJ092712.65+294344.0 hosts a binary, it would be the one closest massive black hole binary system ever discovered.

We study the orbital evolution and accretion history of massive black hole (MBH) pairs in rotationally supported circumnuclear discs up to the point where MBHs form binary systems. Our simulations have high resolution in mass and space which, for the first time, makes it feasible to follow the orbital decay of a MBH either counter- or co-rotating with respect to the circumnuclear disc. We show that a moving MBH on an initially counter-rotating orbit experiences an "orbital angular momentum flip" due to the gas-dynamical friction, i.e., it starts to corotate with the disc before a MBH binary forms. We stress that this effect can only be captured in very high resolution simulations. Given the extremely large number of gas particles used, the dynamical range is sufficiently large to resolve the Bondi-Hoyle-Lyttleton radii of individual MBHs. As a consequence, we are able to link the accretion processes to the orbital evolution of the MBH pairs. We predict that the accretion rate is significantly suppressed and extremely variable when the MBH is moving on a retrograde orbit. It is only after the orbital angular momentum flip has taken place that the secondary rapidly "lights up" at which point both MBHs can accrete near the Eddington rate for a few Myr. The separation of the double nucleus is expected to be around ~10 pc at this stage. We show that the accretion rate can be highly variable also when the MBH is co-rotating with the disc (albeit to a lesser extent) provided that its orbit is eccentric. Our results have significant consequences for the expected number of observable double AGNs at separations of <100 pc.

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.

All the physical processes involved in the formation, merging, and accretion history of massive black holes along the hierarchical build--up of cosmic structures are likely to leave an imprint on the gravitational waves detectable by future space--borne missions, such as LISA. We report here the results of recent studies, carried out by means of dedicated simulations of black hole build--up, aiming at understanding the impact on LISA observations of two ingredients that are crucial in every massive black hole formation scenario, namely: (i) the nature and abundance of the first black hole seeds and (ii) the large gravitational recoils following the merger of highly spinning black holes. We predict LISA detection rates spanning two order of magnitude, in the range 3-300 events per year, depending on the detail of the assumed massive black hole seed model. On the other hand, large recoil velocities do not dramatically compromise the efficiency of LISA observations. The number of detections may drop substantially (by ~60%), in scenarios characterized byabundant light seeds, but if seeds are already massive and/or relatively rare, the detection rate is basically unaffected.

The discovery of hypervelocity binary stars (HVBs) in the Galactic halo would provide definite evidence of the existence of a massive black hole companion to Sgr A*. Here we use an hybrid approach to compute the rate of ejection and the total number of HVBs produced by a hypothetical intermediate-mass black hole (IMBH, M_2<10^5\msun) orbiting Sgr A*. Depending on the mass of M_2 and on the properties of binary stars in the central parsec of the Milky Way, we show that the number of undisrupted HVBs expected to be expelled from the Galactic Center before binary black hole coalescence ranges from zero to a few dozens at most. Therefore, the non-detection of stellar binaries in a complete survey of hypervelocity stars would not rule out the occurrence of an IMBH-Sgr A* in-spiralling event within the last few 10^8 years.

We study the evolution of a massive black hole pair in a rotationally supported nuclear disc. The distributions of stars and gas mimic the nuclear region of a gas-rich galaxy merger remnant. Using high-resolution SPH simulations, we follow the black hole dynamics and trace the evolution of the underlying background, until the black holes form a binary. We find that the gravitational perturbation of the pair creates a core in the disc density profile, hence decreasing the gas-dynamical drag. This leads the newly formed binary to stall at a separation of ~5 pc. In the early phases of the sinking, black holes lose memory of their initial orbital eccentricity if they co-rotate with the disc, as rotation of the gaseous background promotes circularization of the black hole orbits. Circularization is efficient until the black holes bind in a binary, though in the latest stages of the simulations a residual eccentricity > 0.1 is still present. Black holes are treated as sink particles, allowing for gas accretion. We find that accretion strongly depends on the dynamical properties of the black holes, and occurs preferentially after circularization.

We develop a formalism for studying the dynamics of massive black hole binaries embedded in gravitationally-bound stellar cusps, and study the binary orbital decay by three-body interactions, the impact of stellar slingshots on the density profile of the inner cusp, and the properties of the ejected hypervelocity stars (HVSs). We find that the scattering of bound stars shrinks the binary orbit and increases its eccentricity more effectively than that of unbound ambient stars. Binaries with initial eccentricities e>0.3 and/or unequal-mass companions (M_2/M_1<0.1) can decay by three-body interactions to the gravitational wave emission regime in less than a Hubble time. The stellar cusp is significantly eroded, and cores as shallow as \rho\propto r^-0.7 may develop from a pre-existing singular isothermal density profile. A population of HVSs is ejected in the host galaxy halo, with a total mass ~M_2. We scale our results to the scattering of stars bound to Sgr A*, the massive black hole in the Galactic Center, by an inspiraling companion of intermediate mass. Depending on binary mass ratio, eccentricity, and initial slope of the stellar cusp, a core of radius ~0.1 pc typically forms in 1-10 Myr. On this timescale about 500-2500 HVSs are expelled with speeds sufficiently large to escape the gravitational potential of the Milky Way.

We investigate the distribution of massive black holes (MBHs) in the Virgo cluster. Observations suggest that AGN activity is widespread in massive galaxies (M>1e10 solar masses), while at lower galaxy masses star clusters are more abundant, which might imply a limited presence of central black holes in these galaxy-mass regimes. We explore if this possible threshold in MBH hosting, is linked to nature, nurture, or a mixture of both. The nature scenario arises naturally in hierarchical cosmologies, as MBH formation mechanisms typically are efficient in biased systems, which would later evolve into massive galaxies. Nurture, in the guise of MBH ejections following MBH mergers, provides an additional mechanism that is more effective for low mass, satellite galaxies. The combination of inefficient formation, and lower retention of MBHs, leads to the natural explanation of the distribution of compact massive ob jects in Virgo galaxies. If MBHs arrive to the correlation with the host mass and velocity dispersion during merger-triggered accretion episodes, sustained tidal stripping of the host galaxies creates a population of MBHs which lie above the expected scaling between the holes and their host mass, suggesting a possible environmental dependence.

Narrow Line Seyfert 1 galaxies (NLS1s) are generally considered peculiar objects among the broad class of Type 1 active galactic nuclei, due to the relatively small width of the broad lines, strong X-ray variability, soft X-ray continua, weak [OIII], and strong FeII line intensities. The mass M_BH of the central massive black hole (MBH) is claimed to be lighter than expected from known MBH-host galaxy scaling relations, while the accretion rate onto the MBH larger than the average value appropriate to Seyfert 1 galaxies. In this Letter, we show that NLS1 peculiar M_BH and L/L_Edd turn out to be fairly standard, provided that the broad line region is allowed to have a disc-like, rather than isotropic, geometry. Assuming that NLS1s are rather ``normal'' Seyfert 1 objects seen along the disc axis, we could estimate the typical inclination angles from the fraction of Seyfert 1 classified as NLS1s, and compute the geometrical factor relating the observed FWHM of broad lines to the virial mass of the MBH. We show that the geometrical factor can fully account for the "black hole mass deficit" observed in NLS1s, and that L/L_Edd is (on average) comparable to the value of the more common broad line Seyfert 1 galaxies.

The dynamics of two massive black holes in a rotationally supported nuclear disc of 10^8 solar masses is explored using N-Body/SPH simulations. Gas and star particles are co-present in the disc. Described by a Mestel profile, the disc has a vertical support provided by turbulence of the gas, and by stellar velocity dispersion. A primary black hole of 4 million solar masses is placed at the centre of the disc, while a secondary black hole is set initially on an eccentric co-rotating orbit in the disc plane. Its mass is in a 1 to 1, 1 to 4, and 1 to 10 ratio, relative to the primary. With this choice, we mimic the dynamics of black hole pairs released in the nuclear region at the end of a gas-rich galaxy merger. It is found that, under the action of dynamical friction, the two black holes form a close binary in ~10 Myrs. The inspiral process is insensitive to the mass fraction in stars and gas present in the disc and is accompanied by the circularization of the orbit. We detail the gaseous mass profile bound to each black hole that can lead to the formation of two small Keplerian discs, weighing ~2 % of the black hole mass, and of size ~0.01 pc. The mass of the tightly (loosely) bound particles increases (decreases) with time as the black holes spiral into closer and closer orbits. Double AGN activity is expected to occur on an estimated timescale of ~10 Myrs, comparable to the inspiral timescale. The double nuclear point-like sources that may appear during dynamical evolution have typical separations of ~10 pc.

Hypervelocity stars (HVSs) are a natural consequence of the presence of a massive nuclear black hole (Sgr A*) in the Galactic Center. Here we use the Brown et al. sample of unbound and bound HVSs together with numerical simulations of the propagation of HVSs in the Milky Way halo to constrain three plausible ejection mechanisms: 1) the scattering of stars bound to Sgr A* by an inspiraling intermediate-mass black hole (IMBH); 2) the disruption of stellar binaries in the tidal field of Sgr A*; and 3) the two-body scattering of stars off a cluster of stellar-mass black holes orbiting Sgr A*. We compare the predicted radial and velocity distributions of HVSs with the limited-statistics dataset currently available, and show that the IMBH model appears to produce a spectrum of ejection velocities that is too flat. Future astrometric and deep wide-field surveys of HVSs should shed unambiguous light on the stellar ejection mechanism and probe the Milky Way potential on scales as large as 200 kpc.

The formation, merging, and accretion history of massive black holes along the hierarchical build--up of cosmic structures leaves a unique imprint on the background of gravitational waves at mHz frequencies. We study here, by means of dedicated simulations of black hole build--up, the possibility of constraining different models of black hole cosmic evolution using future gravitational wave space--borne missions, such as LISA. We consider two main scenarios for black hole formation, namely, one where seeds are light (~10^2 \msun, remnant of Population III stars), and one where seeds are heavy (>~10^4 \msun, direct collapse). In all the models we have investigated, massive black hole binary coalescences do not produce a stochastic GW background, but rather, a set of individual resolved events. Detection of several hundreds merging events in 3 year LISA mission will be the sign of a heavy seed scenario with efficient formation of black hole seeds in a large fraction of high redshift halos. On the other extreme, a low event rate, about few tens in 3 years, is peculiar of scenarios where either the seeds are light, and many coalescences do not fall into the LISA band, or seeds are massive, but rare. In this case a decisive diagnostic is provided by the shape of the mass distribution of detected events. Light binaries (m<10^4\msun) are predicted in a fairly large number in Population III remnants models, but are totally absent in direct collapse models. Finally, a further, helpful diagnostic of black hole formation models lies in the distribution of the mass ratios in binary coalescences. While heavy seed models predict that most of the detected events involve equal mass binaries, in the case of light seeds, mass ratios are equally distributed in the range 0.1-1.

We study the long-term evolution of massive black hole binaries (MBHBs) at the centers of galaxies using detailed scattering experiments to solve the full three-body problem. Ambient stars drawn from a isotropic Maxwellian distribution unbound to the binary are ejected by the gravitational slingshot. We construct a minimal, hybrid model for the depletion of the loss cone and the orbital decay of the binary, and show that secondary slingshots - stars returning on small impact parameter orbits to have a second super-elastic scattering with the MBHB - may considerably help the shrinking of the pair in the case of large binary mass ratios. In the absence of loss-cone refilling by two-body relaxation or other processes, the mass ejected before the stalling of a MBHB is half the binary reduced mass. About 50% of the ejected stars are expelled ejected in a "burst" lasting ~1E4 yrs M_6^1/4, where M_6 is the binary mass in units of 1E6 Msun. The loss cone is completely emptied in a few bulge crossing timescales, 1E7 yrs M_6^1/4. Even in the absence of two-body relaxation or gas dynamical processes, unequal mass and/or eccentric binaries with M_6 >0.1 can shrink to the gravitational wave emission regime in less than a Hubble time, and are therefore "safe" targets for the planned Laser Interferometer Space Antenna (LISA).

We study, by means of dedicated simulations of massive black hole build-up, the possibility to constraint the existence and nature of the AGN population at z>6 with available and planned X-ray and near infrared space telescopes. We find that X-ray deep-field observations can set important constraints to the faint-end of the AGN luminosity function at very high redshift. Planned X-ray telescopes should be able to detect AGN hosting black holes with masses down to >10^5 Msun (i.e., X-ray luminosities in excess of 10^42 erg s^-1), and can constrain the evolution of the population of massive black hole at early times (6<z<10). We find that this population of AGN should contribute substantially (~25%) to the unresolved fraction of the cosmic X-ray background in the 0.5-10 keV range, and that a significant fraction (~3-4%) of the total background intensity would remain unaccounted even after future X-ray observations. As byproduct, we compute the expected UV background from AGN at z>6 and we discuss the possible role of AGN in the reionization of the Universe at these early epochs, showing that AGN alone can provide enough ionizing photons only in the (improbable) case of an almost completely homogeneous inter-galactic medium. Finally, we show that super-Eddington accretion, suggested by the observed QSOs at z>6, must be a very rare event, confined to black holes living in the highest density peaks.

We investigate how hierarchical models for the co-evolution of the massive black hole (MBH) and AGN population can reproduce the observed faint X-ray counts. We find that the main variable influencing the theoretical predictions is the Eddington ratio of accreting sources. We compare three different models proposed for the evolution of AGN Eddington ratio, f_Edd: constant f_Edd=1, f_Edd decreasing with redshift, and f_Edd depending on the AGN luminosity, as suggested by simulations of galactic mergers including BHs and AGN feedback. We follow the full assembly of MBHs and host halos from early times to the present in a LambdaCDM cosmology. AGN activity is triggered by halo major mergers and MBHs accrete mass until they satisfy the observed correlation with velocity dispersion. We find that all three models can reproduce fairly well the total faint X-ray counts. The redshift distribution is however poorly matched in the first two models. The Eddington ratios suggested by merger simulations predicts no turn-off of the faint end of the AGN optical luminosity function at redshifts z>=1 down to very low luminosity.

We explore the nature of possible electromagnetic counterparts of coalescences of massive black hole binaries at cosmological distances detectable by the Laser Interferometer Space Antenna (LISA). An electromagnetic precursor, during the last year of gravitational wave (GW)-driven inspiral, or an afterglow within few years after coalescence, may highlight the position in the sky of galaxies hosting LISA sources. We show that observations of precursors and afterglows are mutually exclusive, depending on the mass of the primary black hole. Precursors are expected to occur in binaries where the primary (more massive) black hole is heavier than ~10^7 \Msun. They may correspond to on-off states of accretion, i.e., to a bright X-ray source decaying into quiescence before black hole coalescence, and are likely associated to disturbed galaxies showing signs of ongoing starbursts. Coalescences of lighter binaries, with masses <5x10^6 \Msun, lack of any precursor, as gas is expected to be consumed long before the GW-driven orbital decay. Such events would not be hosted by (massive) galaxies with an associated starburst, given the slow binary inspiral time compared to the typical time scale of starbursts. By contrast, coalescence, for such light binaries, is followed by an electromagnetic afterglow, i.e., an off-on accretion state rising in <20 yrs. Using a cosmological merger tree algorithm, we show that future X-ray missions such as XEUS will be able to identify, in 20 yrs operation, almost all the massive BH binary detectable by LISA, and, in only 5 yrs, all the LISA sources at z>6.

We use full three-body scattering experiments to study the ejection of hypervelocity stars (HVSs) by massive black hole binaries (MBHBs) at the center of galaxies. Ambient stars drawn from a Maxwellian distribution unbound to the binary are expelled by the gravitational slingshot. Accurate measurements of thermally averaged hardening, mass ejection, and eccentricity growth rates (H, J, and K) for MBHBs in a fixed stellar background are obtained by numerical orbit integration from initial conditions determined by Monte Carlo techniques. Three-body interactions create a subpopulation of HVSs on nearly radial, corotating orbits, with a spatial distribution that is initially highly flattened in the inspiral plane of the MBHB, but becomes more isotropic with decreasing binary separation. The degree of anisotropy is smaller for unequal mass binaries and larger for stars with higher kick velocities. Eccentric MBHBs produce a more prominent tail of high-velocity stars and break planar symmetry, ejecting HVSs along a broad jet perpendicular to the semimajor axis. The jet two-sidedness decreases with increasing binary mass ratio, while the jet opening-angle increases with decreasing kick velocity and orbital separation. The detection of a numerous population of HVSs in the halo of the Milky Way by the next generation of large astrometric surveys like GAIA may provide a unique signature of the history, nature, and environment of the MBH at the Galactic center.

We study the long-term evolution of massive black hole binaries (MBHBs) at the centers of galaxies using detailed full three-body scattering experiments. Stars, drawn from a distribution unbound to the binary, are ejected by the gravitational slingshot. We quantify the effect of secondary slingshots -- stars returning on small impact parameter orbits to have a second super-elastic scattering with the MBHB -- on binary separation. Even in the absence of two-body relaxation or gas dynamical processes, very unequal mass binaries of mass M=10^7 solar masses can shrink to the gravitational wave emission regime in less than a Hubble time, and are therefore a target for the planned Laser Interferometer Space Antenna (LISA). Three-body interactions create a subpopulation of hypervelocity stars on nearly radial, corotating orbits, with a spatial distribution that is initially highly flattened in the inspiral plane of the MBHB, but becomes more isotropic with decreasing binary separation. The mass ejected is ~0.7 times the binary reduced mass, and most of the stars are ejected in an initial burst lasting much less than a bulge crossing time.

We study the inspiral of double black holes, with masses in the LISA window of detectability, orbiting inside a massive circum-nuclear disc. Using high-resolution SPH simulations, we follow the black hole dynamics in the early phase when gas-dynamical friction acts on the black holes individually, and continue our simulation until they form a close binary. We find that in the early sinking the black holes lose memory of their initial orbital eccentricity if they co-rotate with the gaseous disc, forming a binary with a low eccentricity, consistent with zero within our numerical resolution limit. The cause of circularization resides in the rotation present in the gaseous background where dynamical friction operates. Circularization may hinder gravitational waves from taking over and leading the binary to coalescence. In the case of counter-rotating orbits the initial eccentricity does not decrease, and the black holes may bind forming an eccentric binary. When dynamical friction has subsided, for equal mass black holes and regardless their initial eccentricity, angular momentum loss, driven by the gravitational torque exerted on the binary by surrounding gas, is nevertheless observable down to the smallest scale probed. In the case of unequal masses, dynamical friction remains efficient down to our resolution limit, and there is no sign of formation of any ellipsoidal gas distribution that may further harden the binary. During inspiral, gravitational capture of gas by the black holes occurs mainly along circular orbits: eccentric orbits imply high relative velocities and weak gravitational focusing. Thus, AGN activity may be excited during the black hole pairing process and double active nuclei may form when circularization is completed, on distance-scales of tens of pcs.

A large population of Intermediate Mass Black Holes (IMBHs) might be produced at early cosmic times as a left over of the evolution of the very massive first stars. Accretion onto IMBHs provides a source of (re)ionizing radiation. We show that the baryon mass fraction locked into IMBHs and their growth is strongly constrained by the observed residual Soft X-ray Background (SXRB) intensity. Thus, unless they are extremely X-ray quiet, miniquasars must be quite rare and/or have a short shining phase. As a byproduct, we find that miniquasars can not be the only source of reionization and that their alleged contribution to the near infrared bands is completely negligible.

We compute the expected low-frequency gravitational wave signal from coalescing massive black hole (MBH) binaries at the center of galaxies. We follow the merging history of halos and associated holes via cosmological Monte Carlo realizations of the merger hierarchy from early times to the present in a LCDM cosmology. MBHs get incorporated through a series of mergers into larger and larger halos, sink to the centre owing to dynamical friction, accrete a fraction of the gas in the merger remnant to become more massive, and form a binary system. Stellar dynamical processes dominates the orbital evolution of the binary at large separations, while gravitational wave emission takes over at small radii, causing the final coalescence of the system. We discuss the observability of inspiraling MBH binaries by a low-frequency gravitational wave experiment such as the planned Laser Interferometer Space Antenna (LISA), discriminating between resolvable sources and unresolved confusion noise. Over a 3-year observing period LISA should resolve this GWB into discrete sources, detecting ~90 individual events above a S/N=5 confidence level, while expected confusion noise is well below planned LISA capabilities.

Massive black hole binaries, with masses in the range 1E3-1E8 Msun, are expected to be the most powerful sources of gravitational radiation at mHz frequencies, and hence are among the primary targets for the planned Laser Interferometer Space Antenna (LISA). We extend and refine our previous analysis (Sesana et al. 2004), detailing the gravitational wave signal expected from a cosmological population of massive black hole binaries. As done in our previous paper, we follow the merger history of dark matter halos, the dynamics of the massive black holes they host, and their growth via gas accretion and binary coalescences in a LCDM cosmology. Stellar dynamical processes dominates the orbital evolution of black hole binaries at large separations, while gravitational wave emission takes over at small radii, causing the final coalescence of the pairs. We show that the GW signal from this population, in a 3 year LISA observation, will be resolved into approx 90 discrete events with S/N>5, among which approx. 35 will be observed above threshold until coalescence. These "merging events" involve relatively massive binaries, M=10E5 Msun, in the redshift range 2<z<6. The remaining approx. 55 events come from higher redshift, less massive binaries (M=5E3 Msun at z>6) and, although their S/N integrated over the duration of the observation can be substantial, the final coalescence phase is at too high frequency to be directly observable by space-based interferometers such as LISA. LISA will be able to detect a fraction approx. 90% of all the coalescences of massive black hole binaries occurring at z=5. The residual confusion noise from unresolved massive black hole binaries is expected to be at least an order of magnitude below the estimated stochastic noise.

We compute the expected gravitational wave signal from coalescing massive black hole (MBH) binaries at the center of galaxies in a hierarchical structure formation scenario in which seed holes of intermediate mass form far up in the dark halo merger tree. The merger history of DM halos and MBHs is followed from z=20 to the present in a LCDM cosmology. MBHs get incorporated through halo mergers into larger and larger structures, sink to the center owing to dynamical friction against the DM background, accrete cold material in the merger remnant, and form MBH binary systems. Stellar dynamical interactions cause the hardening of the binary at large separations, while gravitational wave emission takes over at small radii and leads to the final coalescence of the pair. The integrated emission from inspiraling MBH binaries results in a gravitational wave background (GWB). The characteristic strain spectrum has the standard h_c(f)\propto f^{-2/3} behavior only in the range 1E-9<f<1E-6 Hz. At lower frequencies the orbital decay of MBH binaries is driven by the ejection of background stars, and h_c(f) \propto f. At higher frequencies, f>1E-6 Hz, the strain amplitude is shaped by the convolution of last stable circular orbit emission. We discuss the observability of inspiraling MBH binaries by the planned LISA. Over a 3-year observing period LISA should resolve this GWB into discrete sources, detecting ~60 (~250) individual events above a S/N=5 (S/N=1) confidence level. (Abridged)

Motivated by the recent detection by WMAP of a large optical depth to Thomson scattering -- implying a very early reionization epoch -- we assess a scenario where the universe was reionized by `miniquasars' powered by intermediate-mass black holes (IMBHs), the remnants of the first generation of massive stars. Pregalactic IMBHs form within minihalos above the cosmological Jeans mass collapsing at z=24, get incorporated through mergers into larger and larger systems, sink to the center owing to dynamical friction, and accrete cold material. The merger history of dark halos and associated IMBHs is followed by Monte Carlo realizations of the merger hierarchy in a LCDM cosmology. While seed IMBHs that are as rare as the 3.5-sigma peaks of the primordial density field evolve largely in isolation, a significant number of black hole binary systems will form if IMBHs populate the more numerous 3-sigma peaks instead. In the case of rapid binary coalescence a fraction of IMBHs will be displaced from galaxy centers and ejected into the IGM by the `gravitational rocket' effect, rather than accrete and shine as miniquasars. We show that, under a number of plausible assumptions for the amount of gas accreted onto IMBHs and their emission spectrum, miniquasars powered by IMBHs may be responsible for cosmological reionization at z~15. Reionization by miniquasars with a hard spectrum may be more `economical' than stellar reionization, as soft X-rays escape more easily from the dense sites of star formation and travel further than EUV radiation. Energetic photons will make the diffuse IGM warm and weakly ionized prior to the epoch of reionization breakthrough, set an entropy floor, and reduce gas clumping. Future 21 cm observations may detect a preheated, weakly-ionized IGM in emission against the CMB. (abridged)

We investigate a hierarchical structure formation scenario describing the evolution of a Super Massive Black Holes (SMBHs) population. The seeds of the local SMBHs are assumed to be 'pregalactic' black holes, remnants of the first POPIII stars. As these pregalactic holes become incorporated through a series of mergers into larger and larger halos, they sink to the center owing to dynamical friction, accrete a fraction of the gas in the merger remnant to become supermassive, form a binary system, and eventually coalesce. A simple model in which the damage done to a stellar cusps by decaying BH pairs is cumulative is able to reproduce the observed scaling relation between galaxy luminosity and core size. An accretion model connecting quasar activity with major mergers and the observed BH mass-velocity dispersion correlation reproduces remarkably well the observed luminosity function of optically-selected quasars in the redshift range 1<z<5. We finally asses the potential observability of the gravitational wave background generated by the cosmic evolution of SMBH binaries by the planned space-born interferometer LISA.

We propose that radio-quiet quasars and Seyfert galaxies have central black holes powering outflows and jets which propagate only for a short distance, because the velocity of the ejected material is smaller than the escape velocity. We call them "aborted" jets. If the central engine works intermittently, blobs of material may be produced, which can reach a maximum radial distance and then fall back, colliding with the blobs produced later and still moving outwards. These collisions dissipate the bulk kinetic energy of the blobs by heating the plasma, and can be responsible (entirely or at least in part) for the generation of the high energy emission in radio-quiet objects. This is alternative to the more conventional scenario in which the X-ray spectrum of radio-quiet sources originates in a hot (and possibly patchy) corona above the accretion disk. In the latter case the ultimate source of energy of the emission of both the disk and the corona is accretion. Here we instead propose that the high energy emission is powered also by the extraction of the rotational energy of the black hole (and possibly of the disk). By means of Montecarlo simulations we calculate the time dependent spectra and light curves, and discuss their relevance to the X-ray spectra in radio-quiet AGNs and galactic black hole sources. In particular, we show that time variability and spectra are similar to those observed in Narrow Line Seyfert 1 galaxies.

We present a re-analysis of the simultaneous 30 day IUE/XTE observation of NGC 7469 done in 1996. Our main progress, in comparison to previous spectral analysis, is to adopt and fit directly to the data a detailed Comptonization model. This firstly allows to fit simultaneously the data from the UV to the hard X-ray band in a self-consistent way and secondly it gives direct constraints on the physical parameters of the disc-corona system. Our results are completely consistent with a slab geometry where all the observed UV emission is supposed to cross the corona but more photon-starved geometries also give acceptable fits. Whatever the geometry is, the UV seed photon emission appears to be dominated by the reprocessing of the X-rays. We also found very interesting correlations between the different model parameters, the most important one being the anticorrelation between the corona temperature and the UV flux. Such an anticorrelation is clearly inconsistent with a fixed disc-corona configuration and suggests a variation of the geometry of the system. We also find a correlation between the corona optical depth and the X-ray flux which may reflect processes linked to the corona formation. During the observations, NGC 7469 appears to accrete near its Eddington limit. This source could then be an example of magnetically dominated disc-corona system as recently proposed by Merloni (2003).

Mkn 841 has been observed simultaneously by XMM and BeppoSAX in January 2001. Due to operational contingency, the 30ks XMM observation was split into two parts, separated by about 15 hours. We first report the presence of a narrow iron line which appears to be rapidly variable between the two pointings, requiring a non-standard interpretation. We then focus on the analysis of the broad band (0.3-200 keV) continuum using the XMM/EPIC, RGS and SAX/PDS data. The Mkn 841 spectrum is well fitted by a comptonization model in a geometry more photon-fed than a simple slab geometry above a passive disk. It presents a relatively large reflection (R>2) which does not agree with an apparently weak iron line. It also show the presence of a strong soft excess wellfitted by a comptonized spectrum in a cool plasma, suggesting the presence of a multi-temperature corona.

We report on four observations of the black hole candidate GX 339-4 collected in the spectral range 0.4-140 keV with the BeppoSAX satellite during 1997. All our spectra are typical of a black hole candidate in its low-hard state, i.e. are dominated by the hard tail of the emission from the corona. Although the direct emission from the accretion disk in the soft X-ray range is weak, spectra are well described only with models that take into consideration the Compton reflection of hard X-rays from the corona onto the surface of an accretion disk of non negligible extent.

We investigate a hierarchical structure formation scenario in which galaxy stellar cores are created from the binding energy liberated by shrinking supermassive black hole (SMBH) binaries. The binary orbital decay heats the surrounding stars, eroding a preexisting 1/r^2 stellar cusp. We follow the merger history of dark matter halos and associated SMBHs via cosmological Monte Carlo realizations of the merger hierarchy from early times to the present in a LCDM cosmology. Massive black holes get incorporated through a series of mergers into larger and larger halos, sink to the center owing to dynamical friction, accrete a fraction of the gas in the merger remnant to become supermassive, and form a binary system. Stellar dynamical processes drive the binary to harden and eventually coalesce. A simple scheme is applied in which the loss cone is constantly refilled and a constant density core forms due to the ejection of stellar mass. We find that a model in which the effect of the hierarchy of SMBH interactions is cumulative and cores are preserved during galaxy mergers produces at the present epoch a correlation between the `mass deficit' (the mass needed to bring a flat inner density profile to a 1/r^2 cusp) and the mass of the nuclear SMBH, with a normalization and slope comparable to the observed relation. Models in which the mass displaced by the SMBH binary is replenished after every major galaxy merger appear instead to underestimate the mass deficit observed in `core' galaxies.

We assess models for the assembly of supermassive black holes (SMBHs) at the center of galaxies that trace their hierarchical build-up far up in the dark halo `merger tree'. We assume that the first `seed' black holes (BHs) formed in (mini)halos collapsing at z=20 from high-sigma density fluctuations. As these pregalactic holes become incorporated through a series of mergers into larger and larger halos, they sink to the center owing to dynamical friction, accrete a fraction of the gas in the merger remnant to become supermassive, form a binary system, and eventually coalesce. The merger history of dark matter halos and associated BHs is followed by cosmological Monte Carlo realizations of the merger hierarchy. A simple model, where quasar activity is driven by major mergers and SMBHs accrete at the Eddington rate a mass that scales with the fifth power of the velocity dispersion, is shown to reproduce the optical LF of quasars in the redshift range 1<z<4. Binary and triple BH interactions are followed in our merger tree. The assumptions underlying our scenario lead to the prediction of a population of massive BHs wandering in galaxy halos and the intergalactic medium at the present epoch, and contributing <10% to the total BH mass density. At all epochs the fraction of binary SMBHs in galaxy nuclei is of order 10%, while the fraction of binary quasars is less than 0.3%

We report on the detection of a rapidly variable narrow Fe K$\alpha$ line in Mkn 841. The source has been observed two times by XMM-Newton and simultaneously with BeppoSAX. The two observations, of about 10ks long each, were separated by $\sim$ 15 hours. The line flux reaches a maximum during the first observation and is significantly reduced in the second one. The continuum shape and flux, instead, keep roughly constant between the two pointings. Such rapid variability of a narrow (unresolved by the XMM-pn instrument) line has never been reported in the past. These results are not easily explained in the standard cold reflection model where the narrow line component is supposed to be produced far from the primary X-ray source (e.g. from the torus) and is thus not expected to vary rapidly. Different interpretations are discussed.

The evolution of radiation emitted during the dynamical collapse of metal-free protostellar clouds is investigated within a spherically symmetric hydrodynamical scheme that includes the transfer of radiation and the chemistry of the primordial gas. The cloud centre collapses on a time scale of about 10^5-10^6 years, thanks to line cooling from molecular hydrogen (H2). For most of the collapse time, when the evolution proceeds self-similarly, the luminosity slowly rises up to about 10^36 erg/s and is essentially due to H2 IR line emission. Later, continuum IR radiation provides an additional contribution, which is mostly due to the accretion of an infalling envelope upon a small hydrostatic protostellar core which develops in the centre. We follow the beginning of the accretion phase, when the enormous accretion rate (~ 0.1 Msun/yr) produces a very high continuum luminosity of about 10^36 erg/s. Despite the high luminosities, the radiation field is unable to affect the gas dynamics during the collapse and the first phases of accretion, because the opacity of the infalling gas is too small; this is quite different from present-day star formation. We also find that the protostellar evolution is similar among clouds with different initial configurations, including those resulting from 3D cosmological simulations of primordial objects; in particular, the shape of the molecular spectra is quite universal. Finally, we briefly discuss the detectability of this pristine cosmic star formation activity.

We present results from RXTE data obtained during a systematic monitoring program of four Seyfert galaxies (NGC5548, NGC5506, MCG-6-30-15 and NGC4051). We studied the variability of three hardness ratios derived from the light curves in four energy bands. All the objects show similar spectral variations in all ratios. In order to interpret the results we computed the hardness ratios corresponding to a simple spectral model of a power law plus iron line plus reflection component, considering two possibilities: a)variations with constant line equivalent width and reflection parameter R and b) variations with constant line and reflection flux. The overall, mean observed trends can be explained by spectral slope variations and a constant flux Fe line and reflection component, although the existence of a line component which is variable on short time scales cannot be excluded. Finally, we find that the data are not consistent with an increase of R with flux for individual sources, indicating that, as a single source varies, softer spectra do not correspond to larger R values.

We present photometric and spectroscopic observations of the late afterglow of GRB 000911, starting ~1 day after the burst event and lasting ~8 weeks. We detect a moderately significant re-brightening in the R, I and J lightcurves, associated with a sizable reddening of the spectrum. This can be explained through the presence of an underlying supernova, outshining the afterglow ~30 days after the burst event. Alternative explanations are discussed.

We discuss the effects of a comptonizing corona on the appearance of the reflection components, and in particular of the reflection hump, in the X-rays spectra of accreting black holes. Indeed, in the framework of a thermal corona model, we expect that part (or even all, depending on the coronal covering factor) of the reflection features should cross the hot plasma, and thus suffer Compton scattering, before being observed. We have studied in detail the dependence of these effects on the physical (i.e. temperature and optical depth) and geometrical (i.e. inclination angle) parameters of the corona, concentrating on the slab geometry . Due to the smoothing and shifting towards high energies of the comptonized reflection hump, the main effects on the emerging spectra appear above 100 keV. We have also investigated the importance of such effects on the interpretation of the results obtained with the standard fitting procedures. We found that fitting Comptonization models, taking into account comptonized reflection, by the usual cut-off power law + uncomptonized reflection model, may lead to an underestimation of the reflection normalization and an overestimation of the high energy cut-off. We discuss and illustrate the importance of these effects by analysing recent observational results as those of the galaxy NGC 4258. We also find that the comptonizing corona can produce and/or emphasize correlations between the reflection features characteristics (like the iron line equivalent width or the covering fraction) and the X-ray spectral index similar to those recently reported in the literature. We also underline the importance of these effects when dealing with accurate spectral fitting of the X-ray background.

We report on results of a target of opportunity observation of the X-ray transient XTE J1118+480 performed on 2000 April 14-15 with the Narrow Field Instruments (0.1-200 keV) of the SAX satellite. The measured spectrum is a power law with a photon index of ~1.7 modified by an ultrasoft X-ray excess and a high-energy cutoff above ~100 keV. The soft excess is consistent with a blackbody with temperature of ~40 eV and a low flux, while the cut-off power law is well fitted by thermal Comptonization in a plasma with an electron temperature of 100 keV and an optical depth of order of unity. Consistent with the weakness of the blackbody, Compton reflection is weak. Though the data are consistent with various geometries of the hot and cold phases of the accreting gas, we conclude that a hot accretion disk is the most plausible model. The Eddington ratio implied by recent estimates of the mass and the distance is about 10^{-3}, which may indicate that advection is probably not the dominant cooling mechanism. We finally suggest that the reflecting medium has a low metallicity, consistent with location of the system in the halo.

We present and discuss BeppoSAX observations of Mrk841 and Mrk335, two Seyfert 1 galaxies in which previous observations have established the presence of soft excesses. We confirm the soft excess in both sources, even if for Mrk~841 a warm absorber provides a fit almost as good as the one with a true excess. As far as the hard X-ray continuum is concerned, a Comptonization model provides a fit as good as a power law and a physically sound solution for Mrk841. For Mrk335, the Comptonization model gives a result which is somewhat better on statistical ground, but rather problematic on physical ground. The most interesting results regard the reprocessing components. For Mrk841 we find a very large reflection continuum but an almost normal iron line equivalent width even if, within the errors, a solution in which both components are a factor ~2 larger than expected for an accretion disc is still marginally acceptable. If this is the case, an anisotropy of the primary emission seems the best explanation. On the contrary, in Mrk335 we find a very large iron line EW but a reflection component not accordingly large. In this case, the best solution seems to be in terms of reflection from an ionized disc.

In this paper, I will describe the features of the numerical code CUBA, aimed at the solution of the radiative transfer equation in a cosmological context. CUBA will be soon available for public use at the URL http://pitto.mib.infn.it/~haardt/, allowing for several user-supplied input parameters, such as favourite cosmology, luminosity functions, Type II object evolution, stellar spectra, and many others. I will also present some new results of the UV/X-ray cosmic background as produced by the observed populations of QSOs and star forming galaxies, updating and extending our previous works. The background evolution is complemented with a number of derived quantities such as the ionization and thermal state of the IGM, the HeII opacity, the HI and HeII ionization rates, and the HI, HeII and Compton heating rates.

We have used realistic Comptonization models to fit high quality BeppoSAX data of 6 Seyfert galaxies. Our main effort was to adopt a Comptonization model taking into account the anisotropy of the soft photon field. The most important consequence is a reduction of the first scattering order, which produces a break (the so-called anisotropy break) in the outgoing spectra. The physical parameters of the hot corona (i.e. the temperature and optical depth) obtained fitting this class of models to broad band X-ray spectra are substantially different from those derived fitting the same data with the power law + cut-off model commonly used in the literature. In particular, our best fits with Comptonization models in slab geometry give a temperature generally much larger and an optical depth much smaller than derived from the power law + cut-off fits, using standard Comptonization formulae. The estimate of the reflection normalization is also larger with the slab corona model. For most objects of our sample, both models give Compton parameter values larger than expected in a slab corona geometry, suggesting a more ``photon starved'' X-ray source configuration. Finally, the two models provide different trends and correlation between the physical parameters which has major consequences for the physical interpretation of the data.

The nature of Gamma-Ray Burst progenitors is still a debated issue, but consensus is growing on the association of GRBs with massive stars. Furthermore, current models for the reionization of the universe consider massive Pop--III stars as the sources of the ionizing photons. There could then be a natural link between GRBs and reionization. The reionization epoch can be measured through prompt IR spectral observations of high redshift GRBs. For this, GRBs are better than quasars: they produce a smaller HII region even if they are much brighter than quasars (but for a much shorter time) and then, contrary to quasars, they do not modify the absorption properties of the surrounding IGM.

We used high quality BeppoSAX data of 6 Seyfert galaxies to test realistic thermal Comptonization models. Our main effort was to adopt a Comptonization model taking into account the anisotropy of the soft photon field. The best fit parameter values of the temperature and optical depth of the corona and of the reflection normalization obtained fitting this class of models to the data are substantially different from those derived fitting the same data with the power law + cut--off model commonly used. The two models also provide different trends and correlation between the physical parameters, which has major consequences for the physical interpretation of the data

We report on some results from a three-year program of RXTE observations of 4 Seyfert galaxies: NGC 4051, NGC 5548, MCG-6-30-15 and NGC 5506. We focus here on color-color and color-flux diagrams obtained using the count-rates of the PCA instrument in different energy ranges: 3--5, 5--7 and 7--10 keV. The data of the different sources show interesting trends in the different color-color and color-flux planes. These trends are quite similar from one source to the other, even in the case of the narrow line Seyfert 1 NGC 4051 which exhibits strong flux and spectral variability. We discuss these results in term of a simple spectral model: a power law with high energy cut-off + reflection + Iron line, with the scope of understanding the relations between the observed variations of different components and the physical causes of the variability.

The optical/far-IR extragalactic background light (EBL) from both resolved and unresolved extragalactic sources is an indicator of the total luminosity of cosmic structures, as the cumulative emission from young and evolved galactic systems, as well as from active galactic nuclei (AGNs), is recorded in this radiation. This is a brief review of some of the implications of the observed brightness of the night sky for the stellar mass density and average metallicity of the universe today, and of the possible contribution of MACHO progenitors and QSOs to the EBL. Assuming a Salpeter initial mass function with a cutoff below 0.6 M_sun, a lower limit of Omega_{g+s}h^2>0.0015 I_{60} can be derived to the visible (recycled gas + stars) mass density required to generate an EBL at a level of I_{EBL}=60 I_{60} nW/m^2/sr: our latest, best-guess estimate is Omega_{g+s}h^2 ~ 0.0023 I_{60}, which implies a mean metallicity at the present-epoch of y_Z Omega_{g+s}/Omega_b ~ 0.2 Z_sun. If massive dark halos around spiral galaxies are partially composed of faint, old white dwarfs, i.e. if a non-negligible fraction (a few percent) of the nucleosynthetic baryons is locked in the remnants (MACHOs) of intermediate-mass stars forming at very high redshifts, then the bright early phases of such halos should contribute significantly to the observed EBL. Assuming a standard black hole accretion model for quasar activity and using recent observations of the quasar population and new synthesis models for the cosmic X-ray background, we estimate a present mass density of QSO remnants of rho_{BH} ~ 3e5 M_sun/Mpc^3 for a 10% efficiency of accreted mass-to-radiation conversion. The quasar contribution to the brightness of the night sky is I_QSO ~ 2 nW/m^2/sr.

We discuss implications of a stellar-dominated UV background at high redshifts. The composite spectrum of 29 LBGs evaluated by Steidel etal (2000) at <z>=3.4 can be well fit by a stellar population with ongoing star formation, a Salpeter initial mass function, modest or negligible dust reddening, and no intrinsic HI photoelectric absorption. Fading starbursts in which star formation has ceased for 10^7 yr or more cannot reproduce the observed flux shortward of 1 Ryd.The escape fraction of H-ionizing photons must be close to 100 percent for the observed sample of LBGs. The spectrum of ionizing photons produced by a stellar population with ongoing star formation is similar to that of QSOs between 1 and 3 Ryd, becomes softer between 3 and 4 Ryd and drops sharply shortward of 4 Ryd. A galaxy-dominated UV background appears inconsistent with the observed HeII/HI opacity ratio at z=2.4, but might be able to explain the SiIV/CIV abundances measured at z>3 in absorption spectra. A scenario may be emerging where star-forming galaxies reionized intergalactic hydrogen at $z>6$ and dominate the 1 Ryd metagalactic flux at z>3, with quasi-stellar sources taking over at low redshifts. If the large amplitude of the H-ionizing flux estimated by Steidel etal is correct, hydrodynamical simulations of structure formation in the IGM require a baryon density (to explain the observed Lyalpha opacity in absorption spectra) which is similar to or larger than that favoured by recent CMB experiments, and is inconsistent with standard nucleosynthesis values.

We report on observations of the two persistent black hole candidates LMC X-3 and LMCX-1 performed with \BS in October 1997. The flux of LMC X-1 was possibly measured up to 60 keV, but there is a possible confusion with PSR 0540-69. Fits with an absorbed multicolor disk black body are not satisfactory, while the superposition of this model with a power law is acceptable. The sources showed little variations during the observations. However in LMC X-1 some X-ray color dependence on intensity is apparent, indicating a hardening of the spectrum in the second half of the observation. The inner disk radius and temperature change, featuring the same (anti)correlation found in {\it RXTE} data (Wilms et al. 2000). QPOs were searched for. In LMC X-3 none was detected; in LMCX-1 a 3 sigma upper (9% rms) limit is given at 0.07 Hz, the frequency of the QPO discovered with Ginga.

We report on the 0.5--200 keV spectral properties of Cyg X-1 observed at different epochs with the Narrow Field Instruments of the BeppoSAX satellite. The source was in its soft state during the first observation of 1996 June. In the second observation of 1996 September, the source had parameters characteristic to its hard state. A soft X-ray excess, a broad Fe K$\alpha$ line and Compton reflection are clearly detected in both states. The soft-state broad-band continuum is well modeled by a disk blackbody (accounting for the soft excess) and Compton upscattering of the disk photons by a hybrid, thermal/non-thermal, plasma, probably forming a corona above the disk (also giving rise to the Compton-reflection component). In the hard state, the primary hard X-ray spectrum can be well modeled by Compton upscattering of a weak blackbody emission by a thermal plasma at a temperature of $\sim 60$ keV. The soft excess is then explained by thermal Comptonization of the same blackbody emission by another hot plasma cloud characterized by a low value of its Compton parameter. Finally, we find the characteristic ratio of the bolometric flux in the soft state to that in the hard state to be about 3. This value is much more compatible with theories of state transitions than the previously reported (and likely underestimated) value of 1.5.

We test accurate models of Comptonisation spectra over BeppoSAX observations of Seyfert galaxies focusing on the long look at NGC 5548. The hot plasma temperature derived with these models is significantly higher than that obtained fitting the same data with a power law plus high energy cut off model for the continuum. This is due to the fact that in anisotropic geometries Comptonisation spectra show "intrinsic" curvature which moves the fitted high energy cut-off to higher energies. We also show preliminary results of our analysis of a sample of 5 other Seyfert galaxies.

BeppoSAX observations of the Seyfert galaxies Mrk 509 and MCG+8-11-11 are presented. Earlier evidence of a soft excess in Mrk 509 is confirmed. This excess is found to be better represented by a power law than by a black body: with a photon slope of 2.5, its extrapolation matches the flux recorded in the far UV. An ASCA observation, which appeared to exclude the presence of the excess while showing instead evidence of a warm absorber, turns out to be compatible with the coexistence of the excess seen with BeppoSAX and of the warm absorber. The hard power law of Mrk 509 is seen for the first time to be affected by a cut-off at high energies, with an e-folding energy of about 70 keV. In MCG+8-11-11 the cut-off is found at about 170 keV, consistent within the combined errors with a previous estimate from a ASCA+OSSE/CGRO observation. In both objects the reflection component is clearly detected. In Mrk 509 its strength, together with that of the iron K line, indicates a solid angle Omega, subtended by the reprocessing gas in the accretion disk, much less than 2pi, but lacking a valid constraint on the inclination angle this evidence is not as convincing as that found with BeppoSAX in IC 4329A. In MCG+8-11-11 the same parameters are instead consistent with Omega=2pi, but the comparison with an ASCA observation, when the flux level was about 2.5 times weaker, suggests that a substantial fraction of the angle might be associated with gas farther out than the accretion disk.

We test accurate models of Comptonization spectra over the high quality data of the BeppoSAX long look at NGC 5548. The data are well represented by a plane parallel corona with an inclination angle of 30$^{\circ}$, a soft photon temperature of 5 eV and a hot plasma temperature and optical depth of $kT_{\rm e}\simeq$ 360 keV and $\tau\simeq$ 0.1, respectively. If energy balance applies, such values suggest that a more ``photon-starved'' geometry (e.g. a hemispheric region) is necessary. The spectral softening detected during a flare, appears to be associated to a decrease of the heating-to-cooling ratio, indicating a geometric and/or energetic modification of the disk plus corona system. The hot plasma temperature derived with the models above is significantly higher than that obtained fitting the same data with a power law plus high energy cut off model for the continuum. This is due to the fact that in anisotropic geometries Comptonization spectra show "intrinsic" curvature which moves the fitted high energy cut-off to higher energies.

We test accurate models of Comptonization spectra over the high quality data of the BeppoSAX long look at NGC 5548, allowing for different geometries of the scattering region, different temperatures of the input soft photon field and different viewing angles. We find that the BeppoSAX data are well represented by a plane parallel or hemispherical corona viewed at an inclination angle of 30$^{\circ}$. For both geometries the best fit temperature of the soft photons is close to 15$^{+3}_{-9}$ eV. The corresponding best fit values of the hot plasma temperature and optical depth are $kT_{\rm e}\simeq$ 250--260 keV and $\tau\simeq$ 0.16--0.37 for the slab and hemisphere respectively. These values are substantially different from those derived fitting the data with a power-law + cut off approximation to the Comptonization component ($kT_{\rm e}\lta$ 60 keV, $\tau\simeq$ 2.4). This is due to the fact that accurate Comptonization spectra in anisotropic geometries show "intrinsic" curvature which reduces the necessity of a high energy cut-off. The Comptonization parameter derived for the slab model {is} larger than predicted for a two phase plane parallel corona in energy balance, suggesting that a more ``photon-starved'' geometry is necessary. The spectral softening detected during a flare which occurred in the central part of the observation corresponds to a decrease of the Comptonization parameter, probably associated with an increase of the soft photon luminosity, the {hard} photon luminosity remaining constant.

We present $Beppo$SAX observations of the $\gamma $-ray emitting quasars 0836+710, 1510-089 and 2230+114. All the objects have been detected in the PDS up to 100 keV and have extremely flat power-law spectra above 2 keV ($\alpha _x$=0.3--0.5). 0836+710 shows absorption higher than the galactic value and marginal evidence for the presence of the redshifted 6.4 keV Iron line. 1510-089 shows a spectral break around 1 keV, with the low energy spectrum steeper ($\alpha_l$=1.6) than the high energy power-law ($\alpha_h$=0.3). The data are discussed in the light of current Inverse Compton models for the high energy emission.

The X-ray spectrum from a cool disk embedded in an ion supported torus is computed. The interaction of the hot ions with the disk increases the hard X-ray luminosity of the system}. A surface layer of the disk is heated by the protons from the torus. The Comptonized spectrum produced by this layer has a shape that depends only weakly on the incident energy flux and the distance from the accreting compact object. It consists of a `blue bump' of unComptonized soft photons and a flat high energy tail, reminiscent of the observed spectra. The hard tail becomes flatter as the thermalization depth in the cool disk is increased. Further evidence for ion illumination are the Li abundance in the secondaries of low mass X-ray binaries and the 450 keV lines sometimes seen in black-hole transient spectra.

NGC 5548 was observed by BeppoSAX in a single long (8 day) observation from 0.2 to 200 keV. We find (1) the spectral variation of the source is produced by a change of the intrinsic power law slope; (2) a high energy cut-off at $E_c= 115^{+39}_{-27}$ keV with a hint of change of $E_c$ with flux; (3) OVII and OVIII absorption K edges, and a possible blended OVII-OVIII K$\alpha,\beta$ emission feature at $0.54^{+0.07}_{-0.06}$ keV, inconsistent with a purely photoionized gas in equilibrium. We propose that the temperature of the absorbing and emitting gas is $\sim 10^6$ K so that both collisional ionization and photoionization contribute.

We will present some results on the broad--band observations of BeppoSAX of the bright Seyfert galaxies NGC 4151 and NGC 5548.

If gamma-ray bursts originate in dense stellar forming regions, the interstellar material can imprint detectable absorption features on the observed X-ray spectrum. Such features can be detected by existing and planned X-ray satellites, as long as the X-ray afterglow is observed after a few minutes from the burst. Detection of these X-ray features will make possible the determination of the redshift of gamma-ray bursts even when their optical afterglows are severely dimmed by extinction.

Models for gamma-ray burst afterglows envisage an hyper-relativistic fireball that is decelerated in the ambient medium around the explosion site. This interaction produces a shock wave which amplifies the magnetic field and accelerates electrons to relativistic energies, setting the conditions for an efficient production of synchrotron photons. If produced in a region of large-scale ordered magnetic field, synchrotron radiation can be highly polarized. The optical transient associated with GRB 990510 was observed ~18.5 hr after the event and linear polarization in the R band was measured at a level of 1.7 +/- 0.2 %. This is the first detection of linear polarization in the optical afterglow of a gamma-ray burst. We exclude that this polarization is due to dust in the interstellar material, either in our Galaxy or in the host galaxy of the gamma-ray burst. These results provide important new evidence in favor of the synchrotron origin of the afterglow emission, and constrains the geometry of the fireball and/or magnetic field lines.

We describe BeppoSAX observations of the black hole candidates LMC X--1 and LMC X--3 performed in Oct. 1997. Both sources can be modelled by a multicolor accretion disk spectrum, with temperature $\sim 1$ keV. However, there is some evidence that a thin emitting component coexists with the thick disk at these temperatures. In the direction of LMC X--1, we detected a significant emission above 10 keV, which we suspect originates from the nearby source PSR 0540-69. For LMC X--1, we estimate an absorbing column density of $\simeq 6\times 10^{21}$ cm$^{-2}$, which is almost ten times larger than that found for LMC X--3. In both sources, we find no indication of emission or absorption features whatsoever.

A BeppoSAX observation of the BLRG 3C390.3 is reported. For the first time, both the 6.4 keV iron line and a strong reflection hump, produced by the illumination of the primary X-ray emission on cold matter, are detected in this source. The 0.1-100 keV continuum is modeled by an absorbed power law (\Gamma~ 1.8) reflected at high energies by material with a fairly large covering factor (\Omega/2\pi~1). The iron line is intrinsically narrow ($\sigma=73(+207,-73) eV) and has an equivalent width of ~140$ eV. We discuss the results in the context of current models for AGNs and suggest that the primary X-ray power law continuum is probably produced by a hot inner flow, while the reprocessed radiation comes from an outer cold thin disk, and/or from a thick torus at even larger radii. Further observations with BeppoSAX could distinguish between the latter two cases. Beamed radiation associated to the radio jet is unlikely to contribute significantly to the X-ray emission. Finally, an historical study of the column density NH, also reported here, shows that the absorption along the line of sight changes in time. The NH time variability, which is not correlated with that of the primary continuum, seems to imply variations of the geometry of the absorber rather than variations in the ionization state of the gas.

We report broad band (0.1-200 keV) X-ray observations, made by BeppoSAX, of a sample of bright Seyfert 2 galaxies: NGC7172, NGC2110, NGC4507, Mkn 3 and NGC7674. These spectra provide a better understanding of the effects of X-ray reprocessing by cold material in the source and allow to put tighter constraints on the various spectral parameters involved. In particular, the data are used to determine, with less ambiguities than in the past, the shape of the intrinsic continuum emission by means of the high energy data. Within the small sample both Compton thin and Compton thick sources are found according to the expectations of the unified theory.

The history of the transition from a neutral intergalactic medium (IGM) to one that is almost fully ionized can reveal the character of cosmological ionizing sources. We study the evolution of the volume filling factor of HII and HeIII regions in a clumpy IGM, and discuss the implications for rival reionization scenarios of the rapid decline in the space density of radio-loud quasars and of the large population of star-forming galaxies recently observed at z>3. The hydrogen component in a highly inhomogeneous universe is completely reionized when the number of photons emitted above 1 ryd in one recombination time equals the mean number of hydrogen atoms. If stellar sources are responsible for keeping the IGM ionized at z=5, the rate of star formation at this epoch must be comparable or greater than the one inferred from optical observations of galaxies at z=3, and the mean metallicity per baryon in the universe of order 1/500 solar. An early generation of stars in dark matter halos with circular velocities v_circ=50 km/s, possibly one of the main source of UV photons at high-z, could be detectable with the Next Generation Space Telescope. Models in which the quasar emissivity declines rapidly at z>3 predict a late HeII reionization epoch, a feature that could explain the recent detection of patchy HeII Lyman-alpha at z=2.9 by Reimers et al. (1997) and the abrupt change observed by Songaila (1998) at about the same epoch of the SiIV/CIV ratio, but appear unable to provide the required number of hydrogen-ionizing photons at z=5.

If gamma-ray bursts originate in dense stellar forming regions, the interstellar material can imprint detectable absorption features on the observed X-ray spectrum. Such features can be detected by existing and planned X-ray satellites, as long as the X-ray afterglow is observed after a few minutes from the burst. If the column density of the interstellar material exceeds ~10^{23} cm^{-2} there exists the possibility to detect the K_alpha fluorescent iron line, which should be visible for more than one year, long after the X-ray afterglow continuum has faded away. Detection of these X-ray features will make possible the determination of the redshift of gamma-ray bursts even when their optical afterglow is severely dimmed by extinction.

The apparency of low energy absorption features in Blazar LECS spectra has been reported during the early phase of the BeppoSAX mission (e.g. by Grandi et al. 1997 for 3C~273 SVP data; Giommi et al. 1997 for PKS 2155--304 SVP data), with data reduction and analysis tools released prior to July 1997. We set up a robust method to investigate the reliability of absorption features detections in LECS data of bright blazars. The analysis shows that similar features are common to all datasets (with possibly the exception of 3C~273 SVP data). This suggests that most of the features observed from 0.2 to 1 keV are possibly due to miscalibration of the instrument response. The reported problems have been overcome with the latest release of calibrations (September 1997).

GX 339-4 has been observed by BeppoSAX twice in spring 1997 as part of a longer monitoring program. The source was close to the highest levels (50 mCrab) of the extended low state (as measured by the XTE ASM during the last 2 years). Its spectrum was quite hard, similar to the Exosat 1984 off state, but 40 times stronger. The source is detected up to more than 120 keV, enabling the possibility to study its high energy spectrum

We present preliminary results of BeppoSAX AO1 observations of 3C 273 performed in January 1997. We also present a close comparison with data obtained during the satellite SVP, in July 1996. On average, the AO1 flux is about a factor 2 higher than the flux detected during the SVP, and roughly on the middle of the historical X-ray flux range. Power law fits with galactic absorption to all observations yield spectral indices in the range 1.53-1.6, with the spectrum extending from 0.2 to at least up to 200 keV without any significant slope change. The broad band spectrum appears basically featureless, marking a clear difference from the SVP data, where an absorption feature at low energy and a fluorescence iron emission line are present. The lack of cold/warm matter signatures in our data may indicate that, at this "high" level of luminosity, the featureless continuum produced in a relativistic jet overwhelms any thermal and/or reprocessing component, while the two components were at least comparable during the "low" state of July 1996.

The high energy continuum in Seyfert galaxies and galactic black hole candidates is likely to be produced by a thermal plasma. There are difficulties in understanding what can keep the plasma thermal, especially during fast variations of the emitted flux. Particle-particle collisions are too inefficient in hot and rarefied plasmas, and a faster process is called for. We show that cyclo-synchrotron absorption can be such a process: mildly relativistic electrons thermalize in a few synchrotron cooling times by emitting and absorbing cyclo-synchrotron photons. The resulting equilibrium function is a Maxwellian at low energies, with a high energy tail when Compton cooling is important. Assuming that electrons emit completely self absorbed synchrotron radiation and at the same time Compton scatter their own cyclo-synchrotron radiation and ambient UV photons, we calculate the time dependent behavior of the electron distribution function, and the final radiation spectra. In some cases, the 2-10 keV spectra are found to be dominated by thermal synchrotron self-Compton process rather than by thermal Comptonization of UV disk radiation.

We present the results from ASCA and OSSE simultaneous observations of the Seyfert 1.5 galaxy MCG+8-11-11 performed in August-September 1995. The ASCA observations indicate a modest flux increase (20%) in 3 days, possibly correlated to a softening of the 0.6-9 keV spectrum. The spectrum is well described by a hard power law (Gamma=1.64) absorbed by a column density slightly larger than the Galactic value, with an iron line at 6.4 keV of EW=400 eV. The simultaneous OSSE data are characterized by a much softer power law with photon index Gamma=3.0, strongly suggesting the presence of a spectral break in the hard X/soft gamma-ray band. A joint fit to OSSE and ASCA data clearly shows an exponential cut-off at about 300 keV, and strong reflection component. MCG+8-11-11 features a spectral break in the underlying continuum unambiguously. This, together with the inferred low compactness of this source, favours thermal or quasi-thermal electron Comptonization in a structured Corona as the leading process of high energy radiation production.

We report on the first detection above ~10 keV of the archetypal Compton-thick Seyfert 2 galaxy NGC 1068. This detection, obtained with the PDS instrument onboard BeppoSAX in the 20-100 keV range, confirms the hardness of the X--ray spectrum above a few keV (as indicated by ASCA observations and supports models envisaging a mixture of both neutral and ionized reflections of an otherwise invisible nuclear continuum, the neutral reflection component being the dominant one in hard X-rays.

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.

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.

The current status of understanding of the X-ray emission from Seyfert galaxies involves thermal Comptonization of soft photons by mildrelativistic electrons and positrons. I review observational and theoretical arguments supporting such a view, discussing current models proposed for the structure of the innermost part of the accretion flow: extended coronae, small scale flaring blobs, advection dominated accretion disks and small scale cold cloudlets.

We present the analysis of the X-ray data of the young active star HD35850 obtained with ASCA and ROSAT. Our main goal was to see if there is a difference in the elemental abundances of active stars between young and more evolved objects. A two temperature plasma with subsolar abundances, of the order of Z = 0.15 - 0.3, is required to fit the SIS spectra. Similar results are obtained from a ROSAT PSPC observation. Metal abundances of 0.2 - 0.4 the solar value are required to fit both the ASCA and ROSAT data together. From a simultaneous SIS0+SIS1 spectral fit, with 2T plasma models and abundances free to vary in non-solar proportions, we find that, besides N, O and Ne for which we find very low values, all other elements have values relative to solar abundances around 0.2-0.3. These subsolar abundances are in line with those typically observed in more evolved, active stars like RS CVn and Algol-type binaries. The two temperature values required to fit the ASCA SIS spectra are about 0.5 and 1.0 keV. These temperatures, especially the higher one, are lower with respect to the values found for the RS CVn and Algol binaries or for the young star AB Dor, but higher than other single G/K stars. All our data show that this single, late F-type star is actually a very active source, indirectly confirming that this fast rotating star is probably a young object. In the simultaneous fit of the ASCA+ROSAT data, a third temperature is required. However this is not just an addition of a softer component, but is more a redistribution of the dominant temperatures. Indeed, the range spanned by the three temperatures, from 5 to 15 million degrees, is not very large.

The current status of understanding of the X-ray emission from Seyfert galaxies involves Comptonization of soft photons by hot subrelativistic electrons. After briefly reviewing the early theoretical basis for the presence of hot optically thin plasma in or around accretion disks and the key observations that led to develop the presently popular model of an accretion disk with a hot corona, we summarize recent progress in accretion models that take into account energy dissipation and/or angular momentum transport in the corona. Finally, adopting the simple scheme of a homogeneous plane parallel corona, we discuss in detail the dependence of the X-ray spectrum on the coronal parameters. Despite the strong coupling between optical depth and temperature which determines a spectral shape insensitive to their precise values, moderate spectral changes are possible. The spectral variability patterns can be used as diagnostics for coronal physics and should allow to determine whether the optical depth of the corona is dominated by e+e- pairs.

GRS 1758-258 is one of the few persistent hard X-ray emitters (E>100 keV) in the Galaxy. Using the ASCA satellite, we have obtained the first detailed data on GRS 1758-258 in the 1-10 keV range, where previous observations were affected by confusion problems caused by the nearby strong source GX5-1. The spectrum is well described by a power law with photon index 1.7 without strong Fe emission lines. A prominent soft excess, as observed with ROSAT when the hard X-ray flux was in a lower intensity state, was not detected. However, the presence of a soft spectral component, accounting for at most 5% of the 0.1-300 keV flux, cannot be excluded. The accurate measurement of interstellar absorption (N_H=(1.5+-0.1) x 10^22 cm -2) corresponds to an optical extinction which definitely excludes the presence of a massive companion.

We discuss in detail the broad band X-ray variability expected in Seyfert galaxies in the framework of a Compton cooled corona model. Variations in the optical depth tau of the corona cause spectral changes in the Comptonized emission in the sense that the spectrum steepens and the temperature decreases with increasing tau. If the corona is pair dominated, tau is determined by the compactness (proportional to luminosity for fixed size). This yields a definite relation between spectral shape and intensity implying a variation in intensity of a factor 10 for a change of 0.2 in the 2-10 keV spectral index. For low pair density coronae, no definite prediction is possible. In the absence of substantial changes in the structure of the corona, the soft thermal component correlates with the medium-hard Comptonized component. Examples are computed using the ROSAT response matrix, showing that current observations of NGC 5548 and Mrk 766 showing different behaviours can be reconciled with the model. The spectral behaviour of Mrk 766 strongly suggests that the coronal plasma in this source is not pair dominated.

The detection of a HeII absorption trough in the spectra of three high-redshift quasars provides new constraints on the spectral shape of the UV extragalactic background. The relative strengths of the observed flux decrements at the rest-frame wavelengths of 1216 and 304 A require a relatively soft radiation field at 4 ryd compared to 1 ryd. If the ionizing metagalactic flux is dominated by the integrated light from QSOs, and the HeII reionization of the universe was completed well before z~3, then the UV/soft X-ray spectrum of individual quasars at these epochs must, once the cosmological ``filtering'' through material along the line of sight is taken into account, satisfy similar constraints on the average. We model the propagation of AGN-like ionizing radiation through the intergalactic medium using CUBA, a numerical code developed in our prior work. We show that, in order to explain the reported HeII absorption, any thermal component responsible for the ``soft X-ray excess'' observed in the spectra of AGNs at z<0.3 must, at z~3, have a typical temperature >80 eV, and a luminosity not exceeding 20% of that of the ``UV bump''.

We present a detailed calculation of the propagation of AGN-like ionizing radiation through the intergalactic space. We model the ionization state of absorbing clouds, and show that the universe will be more opaque above $4\,$Ryd than previously estimated. Singly ionized helium in \Lya forest clouds and Lyman-limit systems is found to be very efficient in reprocessing soft X-ray, helium-ionizing photons into ultraviolet, hydrogen-ionizing ones. We demonstrate that a significant fraction of the absorbed primary photons (emitted, e.g., by quasar sources) will be reradiated by the photoionized gas through \Lya line emission, two-photon continuum, and recombination continuum radiation. In the light of new data and recent studies, we also reassess the contribution of the QSOs observed in optical surveys to the UV extragalactic background, and find that the stochastic reprocessing of quasar Lyman continuum radiation by hydrogen and helium along the line of sight will significantly affect the amplitude, spectral shape, and fluctuations properties of the metagalactic flux. In a scenario in which QSOs are the primary source of ionizing photons in the universe, the integrated \HI \Lya emission at $z=0$ from photoionized \Lya clouds and Lyman-limit systems is found to be at a level of less than 5\% of current observational limits on the far-UV extragalactic radiation flux. We show that $J_{912}$ increases from $\approx 10^{-23} \uvunits$ at the present epoch to $\approx 5\times 10^{-22}\uvunits$ at $z=2.5$. The attenuated direct quasar emission plus recombination radiation from intergalactic gas appears to provide enough hydrogen-ionizing photons to satisfy the proximity effect at large redshift. The \HeII/\HI ratio in the diffuse intergalactic medium and the \Lya clouds increases from $\approx25$ at

We propose that the X--ray emission from radio quiet AGN and galactic black holes is due to Comptonization of soft thermal photons emitted by the underlying accretion disk in localized structures (blobs). The power per unit area produced by the blobs, impinging on the disk, can easily dominate the radiation internally produced by the disk. In this case the electron temperature and the high energy spectrum can be determined in a similar way as in the previously studied homogeneous model (Haardt \& Maraschi 1991). However in the present model: a) the emitted spectrum is largely independent of the {\it fraction} of gravitational power dissipated in the blobs; b) the X--ray spectrum can be harder depending on a form factor of the blobs; c) the UV (or soft X--ray for galactic objects) luminosity that is not intercepted by the blobs can be larger than the X--ray luminosity. In the framework of a simplified accretion disk $\alpha-\Omega$ dynamo model, we make order of magnitude estimates of the number of active blobs, their size, luminosity and hence their compactness, finding values in agreement with what is observed. The expected UV to X--ray spectra and correlations of X--ray and UV light curves are discussed.

We study $e^{\pm}$ pair plasmas in pair equilibrium, which emit high energy radiation by thermal Comptonization of soft photons. We find that the maximum luminosity to size ratio of the source (i.e. the compactness) depends not only on the hot plasma temperature, but also on the spectral index of the resulting Comptonized spectrum. In the observationally interesting range, sources of same compactness can be hotter if their spectrum is steeper. Instruments observing in the 50--500 keV energy range, such as OSSE on board CGRO, and especially the future SAX satellite, can be more successful in detecting sources moderately steep, the flattest sources being characterized by an high energy cut--off at too low frequencies. For any given pair of values of spectral index and temperature, Comptonization theory alone fixes the $ratio$ of the compactnesses in hard and soft photons. However, if the source is pair dominated, the absolute values of the two compactnesses are fixed. Therefore there is a a one--to--one correspondence between the physical parameters of the source (the compactnesses in soft and hard photons) and the observable quantities (spectral index and temperature). This correspondence can be extremely useful in interpreting the physical behaviour of the sources, especially during variations, and can help discriminating between different models for the high energy emission of compact sources.

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.