The fourth Fermi Large Area Telescope Source Catalog (4FGL) lists over 5000 gamma-ray sources with statistical significance above 4$\sigma$. About 23% of the sources listed in this catalog are unidentified/unassociated gamma-ray sources while ~26% of the sources are classified as blazar candidates of uncertain type (BCUs), lacking optical spectroscopic information. To probe the blazar nature of candidate counterparts of UGSs and BCUs, we started our optical spectroscopic follow up campaign in 2012, which up to date account for more than 350 observed sources. In this paper, the tenth of our campaign, we report on the spectroscopic observations of 37 sources, mostly BCUs, whose observations were carried out predominantly at the Observatorio Astron\'omico Nacional San Pedro M\'artir and the Southern Astrophysical Research Observatory between August 2018 and September 2019. We confirm the BL Lac nature of 27 sources and the flat spectrum radio quasar nature of three sources. The remaining ones are classified as six BL Lacs galaxy-dominated and one normal galaxy. We were also able to measure the redshifts for 20 sources, including 10 BL Lacs. As in previous analyses, the largest fraction of BCUs revealed to be BL Lac objects.

Nearly one-third of the sources in the $Fermi$-LAT catalogs lack a lower energy counterpart, hence being referred as unidentified/unassociated gamma-ray sources (UGSs). In order to firmly classify them, dedicated multifrequency follow-up campaigns are necessary. These will permit to unveil their nature and identify the fraction that could belong to the class of active galaxies known as blazars that is the largest population of extragalactic $\gamma$-ray sources. In $Fermi$-LAT catalogs there are also gamma-ray sources associated with multifrequency blazar-like objects known as Blazars Candidates of Uncertain type (i.e., BCUs) for which follow up spectroscopic campaigns are mandatory to confirm their blazar nature. Thus, in 2013 we started an optical spectroscopic campaign to identify blazar-like objects potential counterparts of UGSs and BCUs. Here we report the spectra of 31 additional targets observed as part of our follow up campaign. Thirteen of them are BCUs for which we acquired spectroscopic observations at Observatorio Astrof\'isico Guillermo Haro (OAGH) and at Southern Astrophysical Research Observatory (SOAR) telescopes, while the rest has been identified thanks to the archival observations available from the Sloan Digital Sky Survey (SDSS). We confirm the blazar nature of all BCUs: three of them are in blazar of quasar type (BZQs) while the remaining ones can be spectroscopically classified as BL Lac objects (BZBs). Then we also discovered 18 BL Lac objects lying within the positional uncertainty regions of UGSs that could be their potential counterparts.

The AGN bolometric correction is a key element to understand BH demographics and compute accurate BH accretion histories from AGN luminosities. However, current estimates still differ from each other by up to a factor of two to three, and rely on extrapolations at the lowest and highest luminosities. Here we revisit this fundamental issue presenting general hard X-ray ($K_{X}$) and optical ($K_{O}$) bolometric corrections, computed combining several AGN samples spanning the widest (about 7 dex) luminosity range ever used for this kind of studies. We analysed a total of $\sim 1000$ type 1 and type 2 AGN for which a dedicated SED-fitting has been carried out. We provide a bolometric correction separately for type 1 and type 2 AGN; the two bolometric corrections results to be in agreement in the overlapping luminosity range and therefore, for the first time, a universal bolometric correction for the whole AGN sample (both type 1 and type 2) has been computed. We found that $K_{X}$ is fairly constant at $log(L_{BOL}/L_{\odot}) < 11$, while it increases up to about one order of magnitude at $log(L_{BOL}/L_{\odot}) \sim 14.5$. A similar increasing trend has been observed when its dependence on either the Eddington ratio or the BH mass is considered, while no dependence on redshift up to $z\sim3.5$ has been found. On the contrary, the optical bolometric correction appears to be fairly constant (i.e. $K_{O} \sim 5$) whatever is the independent variable. We also verified that our bolometric corrections correctly predict the AGN bolometric luminosity functions. According to this analysis, our bolometric corrections can be applied to the whole AGN population in a wide range of luminosity and redshift.

The masses of supermassive black holes at the centres of local galaxies appear to be tightly correlated with the mass and velocity dispersions of their galactic hosts. However, the local Mbh-Mstar relation inferred from dynamically measured inactive black holes is up to an order-of-magnitude higher than some estimates from active black holes, and recent work suggests that this discrepancy arises from selection bias on the sample of dynamical black hole mass measurements. In this work we combine X-ray measurements of the mean black hole accretion luminosity as a function of stellar mass and redshift with empirical models of galaxy stellar mass growth, integrating over time to predict the evolving Mbh-Mstar relation. The implied relation is nearly independent of redshift, indicating that stellar and black hole masses grow, on average, at similar rates. Matching the de-biased local Mbh-Mstar relation requires a mean radiative efficiency ~0.15, in line with theoretical expectations for accretion onto spinning black holes. However, matching the "raw" observed relation for inactive black holes requires a mean radiative efficiency around 0.02, far below theoretical expectations. This result provides independent evidence for selection bias in dynamically estimated black hole masses, a conclusion that is robust to uncertainties in bolometric corrections, obscured active black hole fractions, and kinetic accretion efficiency. For our fiducial assumptions, they favour moderate-to-rapid spins of typical supermassive black holes, to achieve a mean radiative efficiency ~0.12-0.20. Our approach has similarities to the classic Soltan analysis, but by using galaxy-based data instead of integrated quantities we are able to focus on regimes where observational uncertainties are minimized.

A supermassive black hole has been found at the centre of nearly every galaxy observed with sufficient sensitivity. The masses of these black holes are observed to increase with either the total mass or the mean (random) velocity of the stars in their host galaxies. The origin of these correlations remains elusive. Observational systematics and biases severely limit our knowledge of the local demography of supermassive black holes thus preventing accurate model comparisons and progress in this field. Here we show that the large-scale spatial distribution of local active galactic nuclei (AGN), believed to be accreting supermassive black holes, can constrain the shape and normalization of the black hole-stellar mass relation thus bypassing resolution-related observational biases. In turn, our results can set more stringent constraints on the so-called "radiative efficiency", a fundamental parameter describing the inner physics of supermassive black holes that is closely linked to their spin, geometry, and ability to release energy. The mean value of the radiative efficiency can be estimated by comparing the average total luminous output of AGN with the relic mass density locked up in quiescent supermassive black holes at galaxy centres today. For currently accepted values of the AGN obscured fractions and bolometric corrections, our newest estimates of the local supermassive black hole mass density favour mean radiative efficiencies of ~10-20\%, suggesting that the vast majority of supermassive black holes are spinning moderately to rapidly. With large-scale AGN surveys coming online, our novel methodology will enable even tighter constraints on the fundamental parameters that regulate the growth of supermassive black holes.

NGC 3147 has been considered the best case of a true type 2 AGN: an unobscured AGN, based on the unabsorbed compact X-ray continuum, which lacks a broad line region (BLR). However, the very low luminosity of NGC 3147 implies a compact BLR, which produces very broad lines, hard to detect against the dominant background host galaxy. Narrow (0.1"x0.1") slit HST spectroscopy allowed us to exclude most of the host galaxy light, and revealed an H$\alpha$ line with an extremely broad base (FWZI$\sim27\,000$ km s$^{-1}$). The line profile shows a steep cutoff blue wing and an extended red wing, which match the signature of a mildly relativistic thin accretion disk line profile. It is indeed well fit with a nearly face on thin disk, at $i\sim23^\circ$, with an inner radius at $77\pm15$ r$_g$, which matches the prediction of $62^{+18}_{-14}$ r$_g$ from the $R_{\rm BLR} \sim L^{1/2}$ relation. This result questions the very existence of true type 2 AGN. Moreover, the detection of a thin disk, which extends below 100 r$_g$ in an $L/L_{\rm Edd}\sim10^{-4}$ system, contradicts the current view of the accretion flow configuration at extremely low accretion rates.

We present two catalogs of radio-loud candidate blazars whose WISE mid-infrared colors are selected to be consistent with the colors of confirmed gamma-ray emitting blazars. The first catalog is the improved and expanded release of the WIBRaLS catalog presented by D'Abrusco et al. (2014): it includes sources detected in all four WISE filters, spatially cross-matched with radio source in one of three radio surveys and radio-loud based on their q22 spectral parameter. WIBRaLS2 includes 9541 sources classified as BL Lacs, FSRQs or mixed candidates based on their WISE colors. The second catalog, called KDEBLLACS, based on a new selection technique, contains 5579 candidate BL Lacs extracted from the population of WISE sources detected in the first three WISE passbands ([3.4], [4.6] and [12]) only, whose mid-infrared colors are similar to those of confirmed, gamma-ray BL Lacs. KDBLLACS members area also required to have a radio counterpart and be radio-loud based on the parameter q12, defined similarly to q22 used for the WIBRaLS2. We describe the properties of these catalogs and compare them with the largest samples of confirmed and candidate blazars in the literature. We crossmatch the two new catalogs with the most recent catalogs of gamma-ray sources detected by Fermi LAT instrument. Since spectroscopic observations of candidate blazars from the first WIBRaLS catalog within the uncertainty regions of gamma-ray unassociated sources confirmed that ~90% of these candidates are blazars, we anticipate that these new catalogs will play again an important role in the identification of the gamma-ray sky.

Local samples of quiescent galaxies with dynamically measured black hole masses (Mbh) may suffer from an angular resolution-related selection effect, which could bias the observed scaling relations between Mbh and host galaxy properties away from the intrinsic relations. In particular, previous work has shown that the observed Mbh-Mstar (stellar mass) relation is more strongly biased than the Mbh-sigma (velocity dispersion) relation. Local samples of active galactic nuclei (AGN) do not suffer from this selection effect, as in these samples Mbh is estimated from megamasers and/or reverberation mapping-based techniques. With the exception of megamasers, Mbh-estimates in these AGN samples are proportional to a virial coefficient fvir. Direct modelling of the broad line region suggests that fvir~3.5. However, this results in a Mbh-Mstar relation for AGN which lies below and is steeper than the one observed for quiescent black hole samples. A similar though milder trend is seen for the Mbh-sigma relation. Matching the high-mass end of the Mbh-Mstar and Mbh-sigma relations observed in quiescent samples requires fvir~15 and fvir~7, respectively. On the other hand, fvir~3.5 yields Mbh-sigma and Mbh-Mstar relations for AGN which are remarkably consistent with the expected `intrinsic' correlations for quiescent samples (i.e., once account has been made of the angular resolution-related selection effect), providing additional evidence that the sample of local quiescent black holes is biased. We also show that, as is the case for quiescent black holes, the Mbh-Mstar scaling relation of AGN is driven by velocity dispersion, thus providing additional key constraints to black hole-galaxy co-evolution models.

Ionized outflows, revealed by broad asymmetric wings of the [OIII] line, are commonly observed in AGN but the low intrinsic spatial resolution of observations has generally prevented a detailed characterization of their properties. The MAGNUM survey aims at overcoming these limitations by focusing on the nearest AGN, including NGC 1365, a nearby Seyfert galaxy (D~17 Mpc), hosting a low-luminosity AGN (Lbol ~ 2x10^43 erg/s). We want to obtain a detailed picture of the ionized gas in the central ~5 kpc of NGC 1365 in terms of physical properties, kinematics, and ionization mechanisms. We also aim to characterize the warm ionized outflow as a function of distance from the nucleus and its relation with the nuclear X-ray wind. We employed VLT/MUSE optical integral field spectroscopic observations to investigate the warm ionized gas and Chandra ACIS-S X-ray data for the hot highly-ionized phase. We obtained flux, kinematic, and diagnostic maps of the optical emission lines, which we used to disentangle outflows from disk motions and measure the gas properties down to a spatial resolution of ~70 pc. [OIII] emission mostly traces an AGN-ionized kpc-scale biconical outflow with velocities up to ~200 km/s. H{\alpha} emission traces instead star formation in a circumnuclear ring and along the bar, where we detect non-circular motions. Soft X-rays are mostly due to thermal emission from the star-forming regions, but we could isolate the AGN photoionized component which matches the [OIII] emission. The mass outflow rate of the extended ionized outflow matches that of the nuclear X-ray wind and then decreases with radius. However, the hard X-ray emission from the circumnuclear ring suggests that star formation might contribute to the outflow. The integrated mass outflow rate, kinetic energy rate, and outflow velocity are broadly consistent with the typical relations observed in more luminous AGN.

The active galaxy NGC 1275 lies at the center of the Perseus cluster of galaxies, being an archetypal BH-galaxy system that is supposed to fit well with the M_BH-host scaling relations obtained for quiescent galaxies. Since it harbours an obscured AGN, only recently our group has been able to estimate its black hole mass. Here our aim is to pinpoint NGC 1275 on the less dispersed scaling relations, namely the M_BH-sigma_star and M_BH-L_bul planes. Starting from our previous work Ricci et al. 2017b, we estimate that NGC 1275 falls well outside the intrinsic dispersion of the M_BH-sigma_star plane being ~1.2 dex (in black hole mass) displaced with respect to the scaling relations. We then perform a 2D morphological decomposition analysis on Spitzer/IRAC images at 3.6 mic and find that, beyond the bright compact nucleus that dominates the central emission, NGC 1275 follows a de Vaucouleurs profile with no sign of significant star formation nor clear merger remnants. Nonetheless, its displacement on the M_BH-L_(3.6,bul) plane with respect to the scaling relation is as high as observed in the M_BH-sigma_star.

We have undertaken a multi-band observing program aimed at obtaining a complete census of winds in a sample of WISE/SDSS selected hyper-luminous (WISSH) QSOs at z~2-4. We have analyzed the rest-frame optical (LBT/LUCI and VLT/SINFONI) and UV (SDSS) spectra of 18 randomly selected WISSH QSOs to measure the SMBH mass and study the properties of winds both in the NLR and BLR traced by blueshifted/skewed [OIII] and CIV emission lines, respectively. These WISSH QSOs are powered by SMBH with masses $\ge$10$^9$ Msun accreting at 0.4<$\lambda_{Edd}$<3.1. We have found the existence of two sub-populations characterized by the presence of outflows at different distances from the SMBH. One population ([OIII] sources) exhibits powerful [OIII] outflows, rest-frame EW (REW) of the CIV emission REW$_{CIV}\approx$20-40 A and modest CIV velocity shift (v$_{CIV}^{peak}$) with respect to the systemic redshift (<=2000 km/s). The second population (Weak [OIII] sources), representing ~70% of the analyzed WISSH QSOs, shows weak/absent [OIII] emission and an extremely large v$_{CIV}^{peak}$ (up to ~8000 km/s and REW$_{CIV}$<=20 A). We propose two explanations for the observed behavior of the strength of the [OIII] emission in terms of orientation effects of the line of sight and ionization cone. The dichotomy in the presence of BLR and NLR winds could be likely due to inclination effects considering a polar geometry scenario for the BLR winds. We find a strong correlation with L$_{Bol}$ and an anti-correlation with $\alpha_{ox}$, whereby the higher L$_{Bol}$, the steeper $\alpha_{ox}$ and the larger is the v$_{CIV}^{peak}$. Finally, the observed dependence v$_{CIV}^{peak}\propto L_{Bol}^{0.28\pm0.04}$ is consistent with radiatively driven winds scenario, where strong UV continuum is necessary to launch the wind and a weakness of the X-ray emission is fundamental to prevent overionization of the wind itself.

We imaged with ALMA and ARGOS/LUCI the molecular gas and the dust and stellar continuum in XID2028, an obscured QSO at z=1.593, where the presence of a massive outflow in the ionized gas component traced by the [O III]5007 emission has been resolved up to 10 kpc. This target represents a unique test case to study QSO 'feedback in action' at the peak epoch of AGN-galaxy coevolution. The QSO has been detected in the CO(5-4) transition and in the 1.3mm continuum, at ~30 and ~20 {\sigma} significance respectively, with both emissions confined in the central (<4 kpc) radius area. Our analysis suggests the presence of a fast rotating molecular disc (v~400 km/s) on very compact scales, and well inside the galaxy extent seen in the rest-frame optical light (~10 kpc, as inferred from the LUCI data). Adding available measurements in additional two CO transitions, CO(2-1) and CO(3-2), we could derive a total gas mass of ~10$^{10}$ M$_\odot$, thanks to a critical assessment of CO excitation and the comparison with Rayleigh-Jeans continuum estimate. This translates into a very low gas fraction (<5%) and depletion time scales of 40-75 Myr, reinforcing the result of atypical gas consumption conditions in XID2028, possibly due to feedback effects on the host galaxy. Finally, we also detect at ~5{\sigma} the presence of high velocity CO gas, which we interpret as a signature of galaxy-scale molecular outflow, spatially coincident with the ionised gas outflow. XID2028 represents therefore a unique case where the measurement of total outflowing mass (~500-800 M$_\odot$/yr) including the molecular and atomic components, in both the ionised and neutral phases, has been attempted for a high-z QSO.

The XXL survey currently covers two 25 sq. deg. patches with XMM observations of ~10ks. We summarise the scientific results associated with the first release of the XXL data set, that occurred mid 2016. We review several arguments for increasing the survey depth to 40 ks during the next decade of XMM operations. X-ray (z<2) cluster, (z<4) AGN and cosmic background survey science will then benefit from an extraordinary data reservoir. This, combined with deep multi-$\lambda$ observations, will lead to solid standalone cosmological constraints and provide a wealth of information on the formation and evolution of AGN, clusters and the X-ray background. In particular, it will offer a unique opportunity to pinpoint the z>1 cluster density. It will eventually constitute a reference study and an ideal calibration field for the upcoming eROSITA and Euclid missions.

We present new ALMA observations aimed at mapping molecular gas reservoirs through the CO(3-2) transition in three quasars at $z\simeq2.4$, LBQS 0109+0213, 2QZ J002830.4-281706, and [HB89] 0329-385. Previous [OIII]5007 observations of these quasars showed evidence for ionised outflows quenching star formation in their host galaxies. Systemic CO(3-2) emission has been detected only in one quasar, LBQS 0109+0213, where the CO(3-2) emission is spatially anti-correlated with the ionised outflow, suggesting that most of the molecular gas may have been dispersed or heated in the region swept by the outflow. In all three sources, including the one detected in CO, our constraints on the molecular gas mass indicate a significantly reduced reservoir compared to main-sequence galaxies at the same redshift, supporting a negative feedback scenario. In the quasar 2QZ J002830.4-281706, we tentatively detect an emission line blob blue-shifted by $v\sim-2000$ km/s with respect to the galaxy systemic velocity and spatially offset by 0.2 arcsec (1.7 kpc) with respect to the ALMA continuum peak. Interestingly, such emission feature is coincident in both velocity and space with the ionised outflow as seen in [OIII]5007. This tentative detection must be confirmed with deeper observations but, if real, it could represent the molecular counterpart of the ionised gas outflow driven by the AGN. Finally, in all ALMA maps we detect the presence of serendipitous line emitters within a projected distance $\sim 160$ kpc from the quasars. By identifying these features with the CO(3-2) transition, the serendipitous line emitters would be located within |$\Delta v$|$<$500 km/s from the quasars, hence suggesting an overdensity of galaxies in two out of three quasars.

Type 2 active galactic nuclei (AGN) represent the majority of the AGN population. However, due to the difficulties in measuring their black hole (BH) masses, it is still unknown whether they follow the same BH mass-host galaxy scaling relations valid for quiescent galaxies and type 1 AGN. Here we present the locus of type 2 AGN having virial BH mass estimates in the $M_{BH}-\sigma_\star$ plane. Our analysis shows that the BH masses of type 2 AGN are $\sim0.9$ dex smaller than type 1 AGN at $\sigma_\star\sim 185$ km s$^{-1}$, regardless of the (early/late) AGN host galaxy morphology. Equivalently, type 2 AGN host galaxies have stellar velocity dispersions $\sim 0.2$ dex higher than type 1 AGN hosts at $M_{BH}\sim10^7$ M$_\odot$.

Studying the coupling between the energy output produced by the central quasar and the host galaxy is fundamental to fully understand galaxy evolution. Quasar feedback is indeed supposed to dramatically affect the galaxy properties by depositing large amounts of energy and momentum into the ISM. In order to gain further insights on this process, we study the SEDs of sources at the brightest end of the quasar luminosity function, for which the feedback mechanism is supposed to be at its maximum. We model the rest-frame UV-to-FIR SEDs of 16 WISE-SDSS Selected Hyper-luminous (WISSH) quasars at 1.8 < z < 4.6 disentangling the different emission components and deriving physical parameters of both the nuclear component and the host galaxy. We also use a radiative transfer code to account for the contribution of the quasar-related emission to the FIR fluxes. Most SEDs are well described by a standard combination of accretion disk+torus and cold dust emission. However, about 30% of them require an additional emission component in the NIR, with temperatures peaking at 750K, which indicates the presence of a hotter dust component in these powerful quasars. We measure extreme values of both AGN bolometric luminosity (LBOL > 10^47 erg/s) and SFR (up to 2000 Msun/yr). A new relation between quasar and star-formation luminosity is derived (LSF propto LQSO^(0.73)) by combining several Herschel-detected quasar samples from z=0 to 4. Future observations will be crucial to measure the molecular gas content in these systems, probe the impact between quasar-driven outflows and on-going star-formation, and reveal the presence of merger signatures in their host galaxies.

We report the virial measurements of the BH mass of a sample of 17 type 2 AGN, drawn from the Swift/BAT 70-month 14-195 keV hard X-ray catalogue, where a faint BLR component has been measured via deep NIR (0.8-2.5 $\mu$m) spectroscopy. We compared the type 2 AGN with a control sample of 33 type 1 AGN. We find that the type 2 AGN BH masses span the 5$<$ log(M$_{BH}$ /M$_{\odot}$) $< $7.5 range, with an average log(M$_{BH}$/M$_{\odot}$) = 6.7, which is $\sim$ 0.8 dex smaller than found for type 1 AGN. If type 1 and type 2 AGN of the same X-ray luminosity log($L_{14-195}$/erg s$^{-1}$) $\sim$ 43.5 are compared, type 2 AGN have 0.5 dex smaller BH masses than type 1 AGN. Although based on few tens of objects, this result disagrees with the standard AGN unification scenarios in which type 1 and type 2 AGN are the same objects observed along different viewing angles with respect to a toroidal absorbing material.

We present the first NuSTAR observation of a 'true' Type 2 Seyfert galaxy. The 3-40 keV X-ray spectrum of NGC3147 is characterised by a simple power-law, with a standard {\Gamma}~1.7 and an iron emission line, with no need for any further component up to ~40 keV. These spectral properties, together with significant variability on time-scales as short as weeks (as shown in a 2014 Swift monitoring campaign), strongly support an unobscured line-of-sight for this source. An alternative scenario in terms of a Compton-thick source is strongly disfavoured, requiring an exceptional geometrical configuration, whereas a large fraction of the solid angle to the source is filled by a highly ionised gas, whose reprocessed emission would dominate the observed luminosity. Moreover, in this scenario the implied intrinsic X-ray luminosity of the source would be much larger than the value predicted by other luminosity proxies, like the [OIII]{\lambda}5007 emission line extinction-corrected luminosity. Therefore, we confirm with high confidence that NGC3147 is a true Type 2 Seyfert galaxy, intrinsically characterised by the absence of a BLR.

We present dust yields for asymptotic giant branch (AGB) and super--asymptotic giant branch (SAGB) stars of solar metallicity. Stars with initial mass $1.5~M_{\odot} \leq M_{\rm ini} \leq 3~M_{\odot}$ reach the carbon star stage during the AGB phase and produce mainly solid carbon and SiC. The size and the amount of the carbon particles formed follows a positive trend with themass of the star; the carbon grains with the largest size ($a_{\rm C} \sim 0.2\mu$m) are produced by AGB stars with $M_{\rm ini} = 2.5-3~M_{\odot}$, as these stars are those achieving the largest enrichment of carbon in the surface regions. The size of SiC grains, being sensitive to the surface silicon abundance, keeps around $a_{\rm SiC} \sim 0.1\mu$m. The mass of carbonaceous dust formed is in the range $10^{-4} - 5\times 10^{-3}~M_{\odot}$, whereas the amount of SiC produced is $2\times 10^{-4} - 10^{-3}~M_{\odot}$. Massive AGB/SAGB stars with $M_{\rm ini} > 3~M_{\odot}$ experience HBB, that inhibits formation of carbon stars. The most relevant dust species formed in these stars are silicates and alumina dust, with grain sizes in the range $0.1\mu m < a_{\rm ol} < 0.15\mu$m and $a_{\rm Al_2O_3} \sim 0.07\mu$m, respectively. The mass of silicates produced spans the interval $3.4\times 10^{-3}~M_{\odot} \leq M_{\rm dust} \leq 1.1\times 10^{-2}~M_{\odot}$ and increases with the initial mass of the star.

Variability in the X-rays is a key ingredient in understanding and unveiling active galactic nuclei (AGN) properties. In this band flux variations occur on short time scales (hours) as well as on larger times scales. While short time scale variability is often investigated in single source studies, only few works are able to explore flux variation on very long time scales.This work provides a statistical analysis of the AGN long term X-ray variability. We study variability on the largest time interval ever investigated for the 0.2-2 keV band, up to $\sim 20$ years rest-frame for a sample of 220 sources. Moreover, we study variability for 2,700 quasars up to $\sim 8$ years rest-frame in the same (soft) band.We build our source sample using the 3XMM serendipitous source catalogue data release 5, and data from ROSAT All Sky Survey Bright and Faint source catalogues. In order to select only AGN we use the Sloan Digital Sky Survey quasar catalogues data releases 7 and 12. Combining ROSAT and XMM-Newton observations, we investigate variability using the structure function analysis which describes the amount of variability as a function of the lag between the observations.Our work shows an increase of the structure function up to 20 years. We do not find evidence of a plateau in the structure function on these long time scales.The increase of the structure function at long time lags suggests that variability in the soft X-rays can be influenced by flux variations originated in the accretion disk or that they take place in a region large enough to justify variation on such long time scales.

We discuss the broadband X-ray spectrum of GRS 1734-292 obtained from non-simultaneous XMM-Newton and NuSTAR observations, performed in 2009 and 2014, respectively. GRS1734-292 is a Seyfert 1 galaxy, located near the Galactic plane at $z=0.0214$. The NuSTAR spectrum ($3-80$ keV) is dominated by a primary power-law continuum with $\Gamma=1.65 \pm 0.05$ and a high-energy cutoff $E_c=53^{+11}_{-8}$ keV, one of the lowest measured by NuSTAR in a Seyfert galaxy. Comptonization models show a temperature of the coronal plasma of $kT_e=11.9^{+1.2}_{-0.9}$ keV and an optical depth, assuming a slab geometry, $\tau=2.98^{+0.16}_{-0.19}$ or a similar temperature and $\tau=6.7^{+0.3}_{-0.4}$ assuming a spherical geometry. The 2009 XMM-Newton spectrum is well described by a flatter intrinsic continuum ($\Gamma=1.47^{+0.07}_{-0.03}$) and one absorption line due to Fe\textsc{XXV} K$\alpha$ produced by a warm absorber. Both data sets show a modest iron K$\alpha$ emission line at $6.4$ keV and the associated Compton reflection, due to reprocessing from neutral circumnuclear material.

Models and observations suggest that both power and effects of AGN feedback should be maximised in hyper-luminous (L_Bol>10^47 erg/s) quasars, i.e. objects at the brightest end of the AGN luminosity function. We present the first results of a multi-wavelength observing program, focusing on a sample of WISE/SDSS selected hyper-luminous (WISSH) broad-line quasars at z~1.5-5. The WISSH quasars project has been designed to reveal the most energetic AGN-driven outflows, estimate their occurrence at the peak of quasar activity and extend the study of correlations between outflows and nuclear properties up to poorly-investigated extreme AGN luminosities (L_Bol~10^47 -10^48 erg/s). We present NIR, long-slit LBT/LUCI1 spectroscopy of five WISSH quasars at z~2.3-3.5 showing prominent [OIII] emission lines with broad (FWHM~1200-2200 km/s) and skewed profiles. The luminosities of the broad [OIII] wings are the highest measured so far (L_[OIII]^broad >~5x10^44 erg/s) and reveal the presence of powerful ionised outflows with mass outflow rates Mdot >~1700 M_Sun/yr and kinetic powers Edot >~10^45 erg/s. Although these estimates are affected by large uncertainties, due to the use of [OIII] as tracer of ionized outflows and the very basic outflow model we assume, these results suggest that the AGN is highly efficient in pushing outwards large amounts of ionised gas in our targets. The mechanical outflow luminosities for WISSH quasars correspond to higher fractions (~1-3%) of L_Bol than those derived for lower L_Bol AGN. Our targets host very massive (M_BH>~2x10^9 M_Sun) black holes which are still accreting at a high rate (i.e. a factor of ~0.4-3 of the Eddington limit). These findings demonstrate that WISSH quasars offer the opportunity of probing the extreme end of both luminosity and SMBH mass functions and revealing powerful ionised outflows able to affect the evolution of their host galaxies.

We present medium resolution near infrared spectroscopic observations of 41 obscured and intermediate class AGN (type 2, 1.9 and 1.8; AGN2) with redshift $z \lesssim$0.1, selected from the Swift/BAT 70-month catalogue. The observations have been carried out in the framework of a systematic study of the AGN2 near infrared spectral properties and have been executed using ISAAC/VLT, X-shooter/VLT and LUCI/LBT, reaching an average S/N ratio of $\sim$30 per resolution element. For those objects observed with X-shooter we also obtained simultaneous optical and UV spectroscopy. We have identified a component from the broad line region in 13 out of 41 AGN2, with FWHM ${\rm > 800 }$ km/s. We have verified that the detection of the broad line region components does not significantly depend on selection effects due to the quality of the spectra, the X-ray or near infrared fluxes, the orientation angle of the host galaxy or the hydrogen column density measured in the X-ray band. The average broad line region components found in AGN2 has a significantly (a factor 2) smaller FWHM if compared with a control sample of type 1 AGN.

Accurately weigh the masses of SMBH in AGN is currently possible for only a small group of local and bright broad-line AGN through reverberation mapping (RM). Statistical demographic studies can be carried out considering the empirical scaling relation between the size of the BLR and the AGN optical continuum luminosity. However, there are still biases against low-luminosity or reddened AGN, in which the rest-frame optical radiation can be severely absorbed/diluted by the host and the BLR emission lines could be hard to detect. Our purpose is to widen the applicability of virial-based SE relations to reliably measure the BH masses also for low-luminosity or intermediate/type 2 AGN that are missed by current methodology. We achieve this goal by calibrating virial relations based on unbiased quantities: the hard X-ray luminosities, in the 2-10 keV and 14-195 keV bands, that are less sensitive to galaxy contamination, and the FWHM of the most important rest-frame NIR and optical BLR emission lines. We built a sample of RM AGN having both X-ray luminosity and broad optical/NIR FWHM measurements available in order to calibrate new virial BH mass estimators. We found that the FWHM of the H$\alpha$, H$\beta$ and NIR lines (i.e. Pa$\alpha$, Pa$\beta$ and HeI$\lambda$10830) all correlate each other having negligible or small offsets. This result allowed us to derive virial BH mass estimators based on either the 2-10 keV or 14-195 keV luminosity. We took also into account the recent determination of the different virial coefficients $f$ for pseudo and classical bulges. By splitting the sample according to the bulge type and adopting separate $f$ factors we found that our virial relations predict BH masses of AGN hosted in pseudobulges $\sim$0.5 dex smaller than in classical bulges. Assuming the same average $f$ factor for both populations, a difference of $\sim$0.2 dex is still found.

The cosmological process of hydrogen (HI) reionization in the intergalactic medium is thought to be driven by UV photons emitted by star-forming galaxies and ionizing active galactic nuclei (AGN). The contribution of QSOs to HI reionization at $z>4$ has been traditionally believed to be quite modest. However, this view has been recently challenged by new estimates of a higher faint-end UV luminosity function (LF). To set firmer constraints on the emissivity of AGN at $z<6$, we here make use of complete X-ray selected samples including deep Chandra and new COSMOS data, capable to efficiently measure the 1 ryd comoving AGN emissivity up to $z\sim5-6$ and down to five magnitudes fainter than probed by current optical surveys, without any luminosity extrapolation. We find good agreement between the logN$\rm{_H}\lesssim21-22$ cm$^{-2}$ X-ray LF and the optically-selected QSO LF at all redshifts for $M_{1450}\leq -23$. The full range of the logN$\rm{_H}\lesssim21-22$ cm$^{-2}$ LF ($M_{1450} \leq -17$) was then used to quantify the contribution of AGN to the photon budget critical value needed to keep the Universe ionized. We find that the contribution of ionizing AGN at $z = 6$ is as small as 1\% - 7\%, and very unlikely to be greater than 30\%, thus excluding an AGN-dominated reionization scenario.

We recently developed a procedure to recognize gamma-ray blazar candidates within the positional uncertainty regions of the unidentified/unassociated gamma-ray sources (UGSs). Such procedure was based on the discovery that Fermi blazars show peculiar infrared colors. However, to confirm the real nature of the selected candidates, optical spectroscopic data are necessary. Thus, we performed an extensive archival search for spectra available in the literature in parallel with an optical spectroscopic campaign aimed to reveal and confirm the nature of the selected gamma-ray blazar candidates. Here, we first search for optical spectra of a selected sample of gamma-ray blazar candidates that can be potential counterparts of UGSs using the Sloan Digital Sky Survey (SDSS DR12). This search enables us to update the archival search carried out to date. We also describe the state-of-art and the future perspectives of our campaign to discover previously unknown gamma-ray blazars.

Despite the fact that blazars constitute the rarest class among active galactic nuclei (AGNs) they are the largest known population of associated $\gamma$-ray sources. Many of the $\gamma$-ray objects listed in the Fermi-Large Area Telescope Third Source catalog (3FGL) are classified as blazar candidates of uncertain type (BCUs), either because they show multifrequency behaviour similar to blazars but lacking optical spectra in the literature, or because the quality of such spectra is too low to confirm their nature. Here we select, out of 585 BCUs in the 3FGL, 42 BCUs which we identify as probable blazars by their WISE infrared colors and which also have optical spectra that are available in the Sloan Digital Sky Survey (SDSS) and/or Six-Degree Field Galaxy Survey Database (6dFGS). We confirm the blazar nature of all of the sources. We furthermore conclude that 28 of them are BL Lacs, 8 are radio-loud quasars with flat radio spectrum and 6 are BL Lac whose emission is dominated by their host galaxy.

Blazars, one of the most extreme class of active galaxies, constitute so far the largest known population of $\gamma$-ray sources and their number is continuously growing in the Fermi catalogs. However in the latest release of the Fermi catalog there is still a large fraction of sources that are classified as blazar candidate of uncertain type (BCUs) for which optical spectroscopic observations are necessary to confirm their nature and their associations. In addition about 1/3 of the $\gamma$-ray point sources listed in the Third Fermi-LAT Source Catalog (3FGL) are still unassociated and lacking an assigned lower energy counterpart. Since 2012 we have been carrying out an optical spectroscopic campaign to observe blazar candidates to confirm their nature. In this paper, the sixth of the series, we present optical spectroscopic observations for 30 $\gamma$-ray blazar candidates from different observing programs we carried out with the TNG, WHT, OAN, SOAR and Magellan telescopes. We found that 21 out of 30 sources investigated are BL Lac objects while the remaining targets are classified as flat spectrum radio quasars showing the typical broad emission lines of normal quasi stellar objects. We conclude that our selection of $\gamma$-ray blazars candidates based on their multifrequency properties continues to be a successful way to discover potential low-energy counterparts of the FermiUnidentified Gamma-ray Sources (UGSs) and to confirm the nature of BCUs.

The extragalactic $\gamma$-ray sky is dominated by emission from blazars, a peculiar class of active galactic nuclei (AGNs). Many of the $\gamma$-ray sources included in Fermi -Large Area Telescope Third Source catalog (3FGL) are classified as a blazar candidate of uncertain type (BCU) because there is no optical spectra available in the literature to confirm their nature. In 2013 we started a spectroscopic campaign to look for the optical counterparts of the BCUs and of the Unidentified $\gamma$-ray Sources. The main goal of our investigation is to confirm the blazar nature of these sources having peculiar properties as compact radio emission and/or selected on the basis of their infrared (IR) colors. Whenever possible we also determine their redshifts. Here we present the results of the observations carried out in the Northern hemisphere in 2013 and 2014 at Telescopio Nazionale Galilleo (TNG), Kitt Peak National Observatory (KPNO) and Observatorio Astron\'omico Nacional (OAN) in San Pedro M\'artir. In this paper we describe the optical spectra of 25 sources. We confirmed that all the 15 BCUs observed in our campaign and included in our sample are blazars and we estimated the redshift for 3 of them. In addition, we present the spectra for 3 sources classified as BL Lacs in the literature but with no optical spectra available to date. We found that one of them is a quasar (QSO) at a redshift $z$ = 0.208 and the other 2 are BL Lacs. Moreover, we also present 7 new spectra for known blazars listed in the Roma-BZCAT having an uncertain redshift or being classified as BL Lac candidates. We found that one of them, 5BZB J0724+2621 is a `changing look' blazar. According to the spectrum available in the literature it was classified as a BL Lac but in our observation we clearly detected a broad emission line that lead to classify this source as a QSO at $z$=1.17.

In order to investigate the role of absorption in AGN with jets, we have studied the column density distribution of a hard X-ray selected sample of radio galaxies, derived from the INTEGRAL/IBIS and Swift/BAT AGN catalogues (~7-10% of the total AGN population). The 64 radio galaxies have a typical FRII radio morphology and are characterized by high 20-100 keV luminosities (from 10^42 to 10^46 erg/s) and high Eddington ratios (Log L_{Bol}/L_{Edd} typically larger than ~0.01). The observed fraction of absorbed AGN (N_{H} > 10^{22} cm^{-2}) is around 40% among the total sample, and ~75% among type 2 AGN. The majority of obscured AGN are narrow line objects, while unobscured AGN are broad line objects, obeying to the zeroth-order predictions of unified models. A significant anti-correlation between the radio core dominance parameter and the X-ray column density is found. The observed fraction of Compton thick AGN is ~2-3%, in comparison with the 5-7% found in radio-quiet hard X-ray selected AGN. We have estimated the absorption and Compton thick fractions in a hard X-ray sample containing both radio galaxies and non-radio galaxies and therefore affected by the same selection biases. No statistical significant difference was found in the absorption properties of radio galaxies and non radio galaxies sample. In particular, the Compton thick objects are likely missing in both samples and the fraction of obscured radio galaxies appears to decrease with luminosity as observed in hard X-ray non-radio galaxies.

In the last decade, using single epoch virial based techniques in the optical band, it has been possible to measure the central black hole mass on large AGN1 samples. However these measurements use the width of the broad line region as a proxy of the virial velocities and are therefore difficult to be carried out on those obscured (type 2) or low luminosity AGN where the nuclear component does not dominate in the optical. Here we present the optical and near infrared spectrum of the starburst/Seyfert galaxy NGC 6221, observed with X-shooter/VLT. Previous observations of NGC 6221 in the X-ray band show an absorbed (N_H=8.5 +/- 0.4 x 10^21 cm^-2) spectrum typical of a type 2 AGN with luminosity log(L_14-195 keV) = 42.05 erg/s, while in the optical band its spectrum is typical of a reddened (A_V=3) starburst. Our deep X-shooter/VLT observations have allowed us to detect faint broad emission in the H_alpha, HeI and Pa_beta lines (FWHM ~1400-2300 km/s) confirming previous studies indicating that NGC 6221 is a reddened starburst galaxy which hosts an AGN. We use the measure of the broad components to provide a first estimate of its central black hole mass (M_BH = 10^(6.6 +/- 0.3) Msol, lambda_Edd=0.01-0.03), obtained using recently calibrated virial relations suitable for moderately obscured (N_H<10^24 cm^-2) AGN.

Negative feedback from active galactic nuclei (AGN) is considered a key mechanism in shaping galaxy evolution. Fast, extended outflows are frequently detected in the AGN host galaxies at all redshifts and luminosities, both in ionised and molecular gas. However, these outflows are only "potentially" able to quench star formation and we are still missing a decisive evidence of negative feedback in action. Here we present Spectrograph for INtegral Field Observations in the Near Infrared (SINFONI) H- and K-band integral-field spectroscopic observations of two quasars at $z\sim$2.4 characterised by fast, extended outflows detected through the [OIII]$\lambda$5007 line (Carniani et al. 2015). The high signal-to-noise ratio of our observations allows us to identify faint narrow (FWHM $< 500$ km/s), and spatially extended components in [OIII]$\lambda$5007 and H$\alpha$ emission associated with star formation in the host galaxy. Such star-formation powered emission is spatially anti-correlated with the fast outflows. The ionised outflows therefore appear to be able to suppress star formation in the region where the outflow is expanding. However the detection of narrow, spatially-extended H$\alpha$ emission indicates star formation rates of at least $\sim 50-90$ M$_{\odot}$/yr, suggesting either that AGN feedback does not affect the whole galaxy or that many feedback episodes are required before star formation is completely quenched. On the other hand, the narrow H$\alpha$ emission extending along the edges of the outflow cone may lead also to a positive feedback interpretation. Our results highlight the possible twofold role of galaxy-wide outflows in host galaxy evolution.

We investigate the role of supermassive black holes in the global context of galaxy evolution by measuring the host galaxy stellar mass function (HGMF) and the specific accretion rate i.e., lambda_SAR, distribution function (SARDF) up to z~2.5 with ~1000 X-ray selected AGN from XMM-COSMOS. Using a maximum likelihood approach, we jointly fit the stellar mass function and specific accretion rate distribution function, with the X-ray luminosity function as an additional constraint. Our best fit model characterizes the SARDF as a double power-law with mass dependent but redshift independent break whose low lambda_SAR slope flattens with increasing redshift while the normalization increases. This implies that, for a given stellar mass, higher lambda_SAR objects have a peak in their space density at earlier epoch compared to the lower lambda_SAR ones, following and mimicking the well known AGN cosmic downsizing as observed in the AGN luminosity function. The mass function of active galaxies is described by a Schechter function with a almost constant Mstar* and a low mass slope alpha that flattens with redshift. Compared to the stellar mass function, we find that the HGMF has a similar shape and that, up to log((Mstar/Msun)~11.5 the ratio of AGN host galaxies to star forming galaxies is basically constant (~10%). Finally, the comparison of the AGN HGMF for different luminosity and specific accretion rate sub-classes with the phenomenological model prediction by Peng et al. (2010) for the "transient" population, i.e. galaxies in the process of being mass-quenched, reveals that low-luminosity AGN do not appear to be able to contribute significantly to the quenching and that at least at high masses, i.e. Mstar>10^(10.7) Msun , feedback from luminous AGN (log(Lbol)>~46 [erg/s]) may be responsible for the quenching of star formation in the host galaxy.

We analyzed a large sample of radio-loud and radio-quiet quasar spectra at redshift 1.0 < z < 1.2 to compare the inferred underlying quasar continuum slopes (after removal of the host galaxy contribution) with accretion disk models. The latter predict redder (decreasing) alpha_3000 continuum slopes (L_\nu~\nu^alpha at 3000Ang) with increasing black hole mass, bluer alpha_3000 with increasing luminosity at 3000Ang, and bluer alpha_3000 with increasing spin of the black hole, when all other parameters are held fixed. We find no clear evidence for any of these predictions in the data. In particular we find that: (i) alpha_3000 shows no significant dependence on black hole mass or luminosity. Dedicated Monte Carlo tests suggest that the substantial observational uncertainties in the black hole virial masses can effectively erase any intrinsic dependence of alpha_3000 on black hole mass, in line with some previous studies. (ii) The mean slope alpha_3000 of radio-loud sources, thought to be produced by rapidly spinning black holes, is comparable to, or even redder than, that of radio-quiet quasars. Indeed, although quasars appear to become more radio loud with decreasing luminosity, we still do not detect any significant dependence of alpha_3000 on radio loudness. The predicted mean alpha_3000 slopes tend to be bluer than in the data. Disk models with high inclinations and dust extinction tend to produce redder slopes closer to empirical estimates. Our mean alpha_3000 values are close to the ones independently inferred at z<0.5 suggesting weak evolution with redshift, at least for moderately luminous quasars.

AGN-driven outflows are invoked by galaxy evolutionary models to quench star formation and to explain the origin of the relations observed locally between super massive black holes and their host galaxies. This work aims to detect the presence of extended ionised outflows in luminous quasars where we expect the maximum activity both in star formation and in black hole accretion. Currently, there are only a few studies based on spatially resolved observations of outflows at high redshift, $z>2$. We analyse a sample of six luminous (${\rm L>10^{47} \ erg/s}$) quasars at $z\sim2.4$, observed in H-band using the near-IR integral field spectrometer SINFONI at VLT. We perform a kinematic analysis of the [OIII] emission line at $\lambda = 5007\AA$. [OIII] has a complex gas kinematic, with blue-shifted velocities of a few hundreds of km/s and line widths up to 1500 km/s. Using the spectroastrometric method we infer size of the ionised outflows of up to $\sim$2 kpc. The properties of the ionised outflows, mass outflow rate, momentum rate and kinetic power, are correlated with the AGN luminosity. The increase in outflow rate with increasing AGN luminosity is consistent with the idea that a luminous AGN pushes away the surrounding gas through fast outflows driven by radiation pressure, which depends on the emitted luminosity. We derive mass outflow rates of about 6-700 M$_{\odot}$/yr for our sample, which are lower than those observed in molecular outflows. Indeed physical properties of ionised outflows show dependences on AGN luminosity which are similar to those of molecular outflows but indicating that the mass of ionised gas is smaller than that of the molecular one. Alternatively, this discrepancy between ionised and molecular outflows could be explained with different acceleration mechanisms.

In the last decade, using single epoch (SE) virial based spectroscopic optical observations, it has been possible to measure the black hole (BH) mass on large type 1 Active Galactic Nuclei (AGN) samples. However this kind of measurements can not be applied on those obscured type 2 and/or low luminosity AGN where the nuclear component does not dominate in the optical. We have derived new SE relationships, based on the FWHM and luminosity of the broad line region component of the Pabeta emission line and/or the hard X-ray luminosity in the 14-195 keV band, which have the prospect of better working with low luminosity or obscured AGN. The SE relationships have been calibrated in the 10^5-10^9 M_sol mass range, using a sample of AGN whose BH masses have been previously measured using reverberation mapping techniques. Our tightest relationship between the reverberation-based BH mass and the SE virial product has an intrinsic spread of 0.20 dex. Thanks to these SE relations, in agreement with previous estimates, we have measured a BH mass of M_BH =1.7^+1.3_-0.7 X 10^5 M_sol for the low luminosity, type 1, AGN NGC 4395 (one of the smallest active galactic BH known). We also measured, for the first time, a BH mass of M_BH = 1.5^+1.1_-0.6 X 10^7 M_sol for the Seyfert 2 galaxy MCG -01-24-012.

This is a White Paper in support of the mission concept of the Large Observatory for X-ray Timing (LOFT), proposed as a medium-sized ESA mission. We discuss the potential of LOFT for the study of active galactic nuclei. For a summary, we refer to the paper.

We present X-shooter at the Vewry Large Telescope observations of a sample of 10 luminous, X-ray obscured QSOs at z$\sim1.5$ from the XMM-COSMOS survey, expected to be caught in the transitioning phase from starburst to active galactic nucleus (AGN)-dominated systems. The main selection criterion is X-ray detection at bright fluxes (L$>=10^{44}$ erg s$^{-1}$) coupled to red optical-to-NIR-to-MIR colors. Thanks to its large wavelength coverage, X-shooter allowed us to determine accurate redshifts from the presence of multiple emission lines for five out of six targets for which we had only a photometric redshift estimate, with an 80 percent success rate, significantly larger than what is observed in similar programs of spectroscopic follow-up of red QSOs. We report the detection of broad and shifted components in the [OIII]$\lambda\lambda$5007,4959 complexes for six out of eight sources with these lines observable in regions free from strong atmospheric absorptions. The full width half-maximum (FWHM) associated with the broad components are in the range FWHM$\sim900-1600$ km s$^{-1}$, larger than the average value observed in SDSS Type 2 AGN samples at similar observed [OIII] luminosity, but comparable to those observed for QSO/ultraluminous infrared galaxies systems for which the presence of kpc scale outflows have been revealed through integral field unit spectroscopy. Although the total outflow energetics (inferred under reasonable assumptions) may be consistent with winds accelerated by stellar processes, we favour an AGN origin for the outflows given the high outflow velocities oberved (v$>1000$ km s$^{-1}$) and the presence of strong winds also in objects undetected in the far infrared.

We present near-infrared spectra, obtained with SINFONI and XShooter observations at ESO VLT, of nine dusty, red QSOs at 1.2<z<2.6. The sources are hard X-ray detected, characterized by cold absorption (N_H>10^{21} - 10^{22} cm^{-2}) and show a broad Ha component in the NIR spectra. We complement this sample with twelve additional sources taken from the literature with similar properties resulting in a total sample of 21 X-ray obscured, intermediate type (1.8-1.9), dusty reddened QSOs. From the broad Ha line we have computed the BH masses through the virial formula and derived Eddington ratios. Moreover, from optical/IR multi-component SED fitting we have derived the stellar mass of their host galaxies and their SFRs. We find that most of the sources in our sample are hosted in starburst and main sequence star-forming galaxies with Eddington ratios lambda>0.1. We find a strong trend with the BH mass i.e. less massive objects are scattered below and above the local relation while the most massive ones are mainly located above it. We also studied the evolution of these sources on the M_BH-M_star plane compared to a sample of optically blue type--1 QSOs and we find that obscured red QSOs show a ratio of M_BH to M_star that increases with redshift which is consistent with or slightly lower than what has been found for blue QSOs. These sources may represent the blow-out phase at the end of the rapid BH growth and immediately preceding the classical blue QSOs typically sampled in optical surveys. They in fact show evidence of outflows in the ionized gas component, but their BH has already fully formed.

We report the discovery of a luminosity distance estimator using Active Galactic Nuclei (AGN). We combine the correlation between the X-ray variability amplitude and the Black Hole (BH) mass with the single epoch spectra BH mass estimates which depend on the AGN luminosity and the line width emitted by the broad line region. We demonstrate that significant correlations do exist which allows one to predict the AGN (optical or X-ray) luminosity as a function of the AGN X-ray variability and either the HBeta or the PaBeta line widths. In the best case, when the PaBeta is used, the relationship has an intrinsic dispersion of ~0.6 dex. Although intrinsically more disperse than Supernovae Ia, this relation constitutes an alternative distance indicator potentially able to probe, in an independent way, the expansion history of the Universe. With this respect, we show that the new mission concept Athena should be able to measure the X-ray variability of hundreds of AGN and then constrain the distance modulus with uncertainties of 0.1 mag up to z~0.6. We also discuss how, using a new dedicated wide field X-ray telescope able to measure the variability of thousands of AGNs, our estimator has the prospect to become a cosmological probe even more sensitive than current Supernovae Ia samples.

The dust formation process in the winds of Asymptotic Giant Branch stars is discussed, based on full evolutionary models of stars with mass in the range $1$M$_{\odot} \leq$M$\leq 8$M$_{\odot}$, and metallicities $0.001 < Z <0.008$. Dust grains are assumed to form in an isotropically expanding wind, by growth of pre--existing seed nuclei. Convection, for what concerns the treatment of convective borders and the efficiency of the schematization adopted, turns out to be the physical ingredient used to calculate the evolutionary sequences with the highest impact on the results obtained. Low--mass stars with M$\leq 3$M$_{\odot}$ produce carbon type dust with also traces of silicon carbide. The mass of solid carbon formed, fairly independently of metallicity, ranges from a few $10^{-4}$M$_{\odot}$, for stars of initial mass $1-1.5$M$_{\odot}$, to $\sim 10^{-2}$M$_{\odot}$ for M$\sim 2-2.5$M$_{\odot}$; the size of dust particles is in the range $0.1 \mu$m$\leq a_C \leq 0.2\mu$m. On the contrary, the production of silicon carbide (SiC) depends on metallicity. For $10^{-3} \leq Z \leq 8\times 10^{-3}$ the size of SiC grains varies in the range $0.05 \mu {\rm m} < {\rm a_{SiC}} < 0.1 \mu$m, while the mass of SiC formed is $10^{-5}{\rm M}_{\odot} < {\rm M_{SiC}} < 10^{-3}{\rm M}_{\odot}$. Models of higher mass experience Hot Bottom Burning, which prevents the formation of carbon stars, and favours the formation of silicates and corundum. In this case the results scale with metallicity, owing to the larger silicon and aluminium contained in higher--Z models. At Z=$8\times 10^{-3}$ we find that the most massive stars produce dust masses $m_d \sim 0.01$M$_{\odot}$, whereas models of smaller mass produce a dust mass ten times smaller. The main component of dust are silicates, although corundum is also formed, in not negligible quantities ($\sim 10-20\%$).

We present X-ray (Chandra), radio (EVLA and EVN), and archival optical data of the triple-merging system in NGC3341. Our panchromatic analysis confirms the presence of a Seyfert 2 AGN in NGC3341B, one of the secondary dwarf companions. On the other hand, the nucleus of the primary galaxy, consistent with a star-forming region of a few solar masses per year, and NGC3341C are very unlikely to host an AGN. We therefore suggest that NGC3341 is an exceptional example of an AGN triggered in the satellite galaxy of a minor-merging system. The existence of such a system can have important implications in the models of hierarchical growth of structures. Further observational and theoretical efforts on NGC3341 and potentially similar sources are needed in order to understand the role of minor mergers on the onset of AGN activity, and in the evolution of massive galaxies.

We present new models for the evolution of stars with mass in the range 1Msun < M < 7.5Msun, followed from the pre-main-sequence through the asymptotic giant branch phase. The metallicity adopted is $Z=3*10^{-4} (which, with an alpha-enhancement of +0.4, corresponds to [Fe/H]=-2). Dust formation is described by following the growth of dust grains of various types as the wind expands from the stellar surface. Models with mass M>3Msun experience Hot Bottom Burning, thus maintaining the surface C/O below unity. Unlike higher Z models, the scarcity of silicon available in the envelope prevents the formation of silicates in meaningful quantities, sufficient to trigger the acceleration of the wind via radiation pressure on the dust grains formed. No silicate formation occurs below a threshold metallicity of Z=10^{-3}. Low--mass stars, with M< 2.5Msun become carbon stars, forming solid carbon dust in their surroundings. The total dust mass formed depends on the uncertain extent of the inwards penetration of the convective envelope during the Third Dredge--Up episodes following the Thermal Pulses. Carbon grains have sizes 0.08 micron < a_C < 0.12 micron and the total amount of dust formed (increasing with the mass of the star) is M_C=(2-6)*10^{-4}Msun. Our results imply that AGB stars with Z=3*10^{-4} can only contribute to carbon dust enrichment of the interstellar medium on relatively long timescales, > 300 Myr, comparable to the evolutionary time of a 3Msun star. At lower metallicities the scarcity of silicon available and the presence of Hot Bottom Burning even in M< 2Msun, prevents the formation of silicate and carbon grains. We extrapolate our conclusion to more metal--poor environments, and deduce that at Z < 10^{-4} dust enrichment is mostly due to metal condensation in supernova ejecta.

We have studied the dependence of the AGN nuclear radio (1.4 GHz) luminosity on both the AGN 2-10 keV X-ray and the host-galaxy K-band luminosity. A complete sample of 1268 X-ray selected AGN (both type 1 and type 2) has been used, which is the largest catalogue of AGN belonging to statistically well defined samples where radio, X and K band information exists. At variance with previous studies, radio upper limits have been statistically taken into account using a Bayesian Maximum Likelihood fitting method. It resulted that a good fit is obtained assuming a plane in the 3D L_R-L_X-L_K space, namely logL_R= xi_X logL_X + xi_K logL_K + xi_0, having a ~1 dex wide (1 sigma) spread in radio luminosity. As already shown, no evidence of bimodality in the radio luminosity distribution was found and therefore any definition of radio loudness in AGN is arbitrary. Using scaling relations between the BH mass and the host galaxy K-band luminosity, we have also derived a new estimate of the BH fundamental plane (in the L_5GHz -L_X-M_BH space). Our analysis shows that previous measures of the BH fundamental plane are biased by ~0.8 dex in favor of the most luminous radio sources. Therefore, many AGN studies, where the BH fundamental plane is used to investigate how AGN regulate their radiative and mechanical luminosity as a function of the accretion rate, or many AGN/galaxy co-evolution models, where radio-feedback is computed using the AGN fundamental plane, should revise their conclusions.

We present the results of a campaign of simultaneous X-ray and optical observations of `true' Type 2 Seyfert galaxies candidates, i.e. AGN without a Broad Line Region (BLR). Out of the initial sample composed by 8 sources, one object, IC1631, was found to be a misclassified starburst galaxy, another, Q2130-431, does show broad optical lines, while other two, IRAS01428-0404 and NGC4698, are very likely absorbed by Compton-thick gas along the line of sight. Therefore, these four sources are not unabsorbed Seyfert 2s as previously suggested in the literature. On the other hand, we confirm that NGC3147, NGC3660, and Q2131-427 belong to the class of true Type 2 Seyfert galaxies, since they do not show any evidence for a broad component of the optical lines nor for obscuration in their X-ray spectra. These three sources have low accretion rates ($\dot{m}= L_{bol}/L_{Edd} \la0.01$), in agreement with theoretical models which predict that the BLR disappears below a critical value of $L_{bol}/L_{Edd}$. The last source, Mrk273x, would represent an exception even of this accretion-dependent versions of the Unification Models, due to its high X-ray luminosity and accretion rate, and no evidence for obscuration. However, its optical classification as a Seyfert 2 is only based on the absence of a broad component of the H\beta, due to the lack of optical spectra encompassing the H\alpha\ band.

Based on previous spectral and temporal optical studies, Q2122-444 has been classified as a naked AGN or true type 2 AGN, that is, an AGN that genuinely lacks a broad line region (BLR). Its optical spectrum seemed to possess only narrow forbidden emission lines that are typical of type 2 (obscured) AGNs, but the long-term optical light curve, obtained from a monitoring campaign over more than two decades, showed strong variability, apparently ruling out the presence of heavy obscuration. Here, we present the results from a ~40 ks XMM-Newton observation of Q2122-444 carried out to shed light on the energetics of this enigmatic AGN. The X-ray analysis was complemented with ATCA radio data to assess the possible presence of a jet, and with new NTT/EFOSC2 optical spectroscopic data to verify the actual absence of a BLR. The higher-quality optical data revealed the presence of strong and broad Balmer lines that are at odds with the previous spectral classification of this AGN. The lack of detection of radio emission rules out the presence of a jet. The X-ray data combined with simultaneous UV observations carried out by the OM aboard XMM-Newton confirm that Q2122-444 is a typical type-1 AGN without any significant intrinsic absorption. New estimates of the black hole mass independently obtained from the broad Balmer lines and from a new scaling technique based on X-ray spectral data suggest that Q2122-444 is accreting at a relatively high rate in Eddington units.

We report the detection of broad Halpha emission in three X-ray selected obscured AGNs at z=1-2. By exploiting the Halpha width and the intrinsic X-ray luminosity, we estimate their black hole masses, which are in the range 0.1-3x10^9 Msun. By means of multi-band photometric data, we measure the stellar mass of their host galaxy and, therefore, infer their M_BH/M_star ratio. These are the first obscured AGNs at high-z, selected based on their black hole accretion (i.e. on the basis of their X-ray luminosity), that can be located on the M_BH-M_star relation at high-z. All of these obscured high-z AGNs are fully consistent with the local M_BH-M_star relation. This result conflicts with those for other samples of AGNs in the same redshift range, whose M_BH/M_star ratio departs significantly from the value observed in local galaxies. We suggest that the obscured AGNs in our sample are in an advanced evolutionary stage, have already settled onto the local M_BH-M_star relation, and whose nuclear activity has been temporarily revived by recent galaxy interactions.

We studied the Active Galactic Nuclei (AGN) radio emission from a compilation of hard X-ray selected samples, all observed in the 1.4 GHz band. A total of more than 1600 AGN with 2-10 keV de-absorbed luminosities higher than 10^42 erg/s were used. For a sub-sample of about 50 z\lsim 0.1 AGN it was possible to reach a ~80% fraction of radio detections and therefore, for the first time, it was possible to almost completely measure the probability distribution function of the ratio between the radio and the X-ray luminosity Rx=log[L(1.4)/Lx]. The probability distribution function of Rx was functionally fitted as dependent on the X-ray luminosity and redshift, P(Rx|Lx,z). It roughly spans over 6 decades (-7<Rx<-1), and does not show any sign of bi-modality. It resulted that the probability of finding large values of the Rx ratio increases with decreasing X-ray luminosities and (possibly) with increasing redshift. No statistical significant difference was found between the radio properties of the X-ray absorbed and unabsorbed AGN. The measure of the probability distribution function of Rx allowed us to compute the kinetic luminosity function and the kinetic energy density which, at variance with what assumed in many galaxy evolution models, is observed to decrease of about a factor of five at redshift below 0.5. About half of the kinetic energy density results to be produced by the more radio quiet (Rx<-4) AGN. In agreement with previous estimates, the AGN efficiency in converting the accreted mass energy into kinetic power is, on average, ~5x10-3.

We present a catalog of optical spectroscopic identifications of sources detected by Spitzer at 3.6 or 24 micron down to 10 and 280 microJy, respectively, in the SWIRE/XMM-Newton/ELAIS-S1 field and classified via line width analysis and diagnostic diagrams. A total of 1376 sources down to R~24.2 mag have been identified (1362 detected at 3.6 micron, 419 at 24 micron, and 405 at both) by low-resolution optical spectroscopy carried out with FORS2, VIMOS, and EFOSC2 at the Very Large Telescope and 3.6m ESO telescopes. The spectroscopic campaigns have been carried out over the central 0.6 square degrees area of ELAIS-S1 which, in particular, has also been observed by XMM-Newton and Chandra. We find the first direct optical spectroscopic evidence that the fraction of active galactic nuclei (AGN; mostly AGN2) increases with increasing F(24 micron)/F(R) ratio, reaching values of 70(+/-20)% in the range 316<F(24 micron)F(R)<1000. We present an IRAC-MIPS color-color diagram able to separate AGN1 from obscured AGN2 candidates. After having corrected for the spectroscopic incompleteness of our sample, it results that the AGN fraction at F(24 micron)=0.8 mJy is ~22(+/-7)% and decreases slowly to ~19(+/-5)% down to F(24 micron)=0.3 mJy.

We present multi-wavelength observations (from optical to sub-millimeter, including Spitzer and SCUBA) of H2XMMJ 003357.2-120038 (also GD158_19), an X-ray selected, luminous narrow-line (Type 2) quasar at z=1.957 selected from the HELLAS2XMM survey. Its broad-band properties can be reasonably well modeled assuming three components: a stellar component to account for the optical and near-IR emission, an AGN component (i.e., dust heated by an accreting active nucleus), dominant in the mid-IR, with an optical depth at 9.7 micron along the line of sight (close to the equatorial plane of the obscuring matter) of tau(9.7)=1 and a full covering angle of the reprocessing matter (torus) of 140 degrees, and a far-IR starburst component (i.e., dust heated by star formation) to reproduce the wide bump observed longward of 70 micron. The derived star-formation rate is about 1500 solar masses per year. The overall modeling indicates that GD158_19 is a high-redshift X-ray luminous, obscured quasar with coeval powerful AGN activity and intense star formation. It is probably caught before the process of expelling the obscuring gas has started, thus quenching the star formation.

We present results on low-resolution mid-infrared (MIR) spectra of 70 infrared-luminous galaxies obtained with the Infrared Spectrograph (IRS) onboard Spitzer. We selected sources from the European Large Area Infrared Survey (ELAIS) with S15 > 0.8 mJy and photometric or spectroscopic z > 1. About half of the sample are QSOs in the optical, while the remaining sources are galaxies, comprising both obscured AGN and starbursts. We classify the spectra using well-known infrared diagnostics, as well as a new one that we propose, into three types of source: those dominated by an unobscured AGN (QSOs), obscured AGN, and starburst-dominated sources. Starbursts concentrate at z ~ 0.6-1.0 favored by the shift of the 7.7-micron PAH band into the selection 15 micron band, while AGN spread over the 0.5 < z < 3.1 range. Star formation rates (SFR) are estimated for individual sources from the luminosity of the PAH features. An estimate of the average PAH luminosity in QSOs and obscured AGN is obtained from the composite spectrum of all sources with reliable redshifts. The estimated mean SFR in the QSOs is 50-100 Mo yr^-1, but the implied FIR luminosity is 3-10 times lower than that obtained from stacking analysis of the FIR photometry, suggesting destruction of the PAH carriers by energetic photons from the AGN. The SFR estimated in obscured AGN is 2-3 times higher than in QSOs of similar MIR luminosity. This discrepancy might not be due to luminosity effects or selection bias alone, but could instead indicate a connection between obscuration and star formation. However, the observed correlation between silicate absorption and the slope of the near- to mid-infrared spectrum is compatible with the obscuration of the AGN emission in these sources being produced in a dust torus.

We present optical identifications and a multi-band catalogue of a sample of 478 X-ray sources in the XMM and Chandra surveys of the central 0.6 deg^2 of the ELAIS-S1 field. The optical/infrared counterpart of each X-ray source was identified using R and IRAC 3.6 um bands. This method was complemented by the precise positions obtained through Chandra observations. Approximately 94% of the counterparts are detected in the R band, while the remaining are blank fields in the optical down to R~24.5, but have a near-infrared counterpart detected by IRAC within 6 arcsec from the XMM centroid. The multi-band catalogue contains photometry in ten photometric bands (B to the MIPS 24 um). We determined redshift and classification for 237 sources (~50% of the sample) brighter than R=24. We classified 47% of the sources with spectroscopic redshift as broad-line active galactic nuclei (BL AGNs) with z=0.1-3.5, while sources without broad-lines are about 46% of the spectroscopic sample and are found up to z=2.6. We identified 11 type 2 QSOs among the sources with X/O>8, with z=0.9-2.6, high 2-10 keV luminosity (log(L2-10keV)>=43.8 erg/s) and hard X-ray colors suggesting large absorbing columns at the rest frame (logN_H up to 23.6 cm^-2). BL AGNs show on average blue optical-to-near-infrared colors, softer X-ray colors and X-ray-to-optical colors typical of optically selected AGNs. Conversely, narrow-line sources show redder optical colors, harder X-ray flux ratio and span a wider range of X-ray-to-optical colors. On average the SEDs of high-luminosity BL AGNs resemble the power-law typical of unobscured AGNs. The SEDs of NOT BL AGNs are dominated by the galaxy emission in the optical/near-infrared, and show a rise in the mid-infrared which suggests the presence of an obscured active nucleus.

(abridged) The ESIS survey is the optical follow up of the SWIRE/Spitzer in the ELAIS-S1 region of the sky. In the era of observational cosmology, the main efforts are focused on the study of galaxy evolution and its environmental dependence. Wide area, multiwavelength, extragalactic surveys are needed in order to probe sufficiently large volumes, minimize cosmic variance and find significant numbers of rare objects. We present VIMOS I and z band imaging belonging to the ESIS survey. A total of ~4 deg2 were targeted in I and ~1 deg2 in z. More than 300000 galaxies have been detected in the I band and ~50000 in the z band. Object coordinates are defined within an uncertainty of ~0.2 arcsec r.m.s., with respect to GSC 2.2. We reach a 90% average completeness at 23.1 and 22.5 mag (Vega) in the I and z bands, respectively. On the basis of IRAC colors, we identified galaxies having the 1.6 um stellar peak shifted to z=1-3. The new I, z band data provide reliable constraints to avoid low-redshift interlopers and reinforce this selection. Roughly 1000 galaxies between z=2-3 were identified over the ESIS ~4 deg2, at the SWIRE 5.8 um depth (25.8 uJy at 3sigma). These are the best galaxy candidates to dominate the massive tail (M>1e11 Msun) of the z>2 mass function.

The X-ray spectrum of NGC7213 is known to present no evidence for Compton reflection, a unique result among bright Seyfert 1s. The observed neutral iron K$\alpha$ line, therefore, cannot be associated with a Compton-thick material, like the disc or the torus, but is due to Compton-thin gas, with the Broad Line Region (BLR) as the most likely candidate. To check this hypothesis, a long Chandra HETG observation, together with a quasi-simultaneous optical spectroscopic observation at the ESO NTT EMMI were performed. We found that the iron line is resolved with a FWHM=$2 400^{+1 100}_{-600}$ km/s, in perfect agreement with the value measured for the broad component of the H$\alpha$, $2640^{+110}_{-90}$ km/s. Therefore, NGC7213 is the only Seyfert 1 galaxy whose iron K$\alpha$ line is unambiguously produced in the BLR. We also confirmed the presence of two ionised iron lines and studied them in greater detail than before. The resonant line is the dominant component in the Fe XXV triplet, therefore suggesting an origin in collisionally ionised gas. If this is the case, the blueshift of around 1000 km/s of the two ionised iron lines could be the first measure of the velocity of a starburst wind from its X-ray emission.

We present the broad-band Spectral Energy Distributions (SEDs) of the largest available highly (72%) complete spectroscopic sample of mid-infrared (MIR) selected galaxies and AGN at intermediate redshift. The sample contains 203 extragalactic sources from the 15-micron survey in the ELAIS-SWIRE field S1, all with measured spectroscopic redshift. Most of these sources have full multi-wavelength coverage from the far-UV to the far-infrared and lie in the redshift range 0.1<z<1.3. Due to its size, this sample allows us for the first time to characterise the spectral properties of the sources responsible for the strong evolution observed in the MIR. Based on SED-fitting technique we have classified the MIR sources, identifying AGN signatures in about 50% of them. This fraction is significantly higher than that derived from optical spectroscopy (~29%) and is due in particular to the identification of AGN activity in objects spectroscopically classified as galaxies. It is likely that in most of our objects, the AGN is either obscured or of low-luminosity, and thus it does not dominate the energetic output at any wavelength, except in the MIR, showing up just in the range where the host galaxy SED has a minimum. The fraction of AGN strongly depends on the flux density, with that derived through the SED-fitting being about 20% at S(15)~0.5-1 mJy and gradually increasing up to 100% at S(15)>10 mJy, while that obtained from optical spectroscopy never being >30%, even at the higher flux densities. The results of this work will be very useful for updating all the models aimed at interpreting the deep infrared survey data and, in particular, for constraining the nature and the role of dust-obscured systems in the intermediate/high-redshift Universe.

(Abridged) We selected high-z massive galaxies at 5.8 microns, in the SWIRE ELAIS-S1 field (1 sq. deg.). Galaxies with the 1.6 microns stellar peak redshifted into the IRAC bands (z~1-3, called ``IR-peakers'') were identified. Stellar masses were derived by means of spectro-photometric fitting and used to compute the stellar mass function (MF) at z=1-2 and 2-3. A parametric fit to the MF was performed, based on a Bayesian formalism, and the stellar mass density of massive galaxies above z=2 determined. We present the first systematic study of the very-massive tail of the galaxy stellar mass function at high redshift. A total of 326 sources were selected. The majority of these galaxies have stellar masses in excess of 1e11 Msun and lie at z>1.5. The availability of mid-IR data turned out to be a valuable tool to constrain the contribution of young stars to galaxy SEDs, and thus their M(stars)/L ratio. The influence of near-IR data and of the chosen stellar library on the SED fitting are also discussed. A significant evolution is found not only for galaxies with M~1e11 Msun, but also in the highest mass bins considered. The comoving number density of these galaxies was lower by more than a factor of 10 at z=2-3, with respect to the local estimate. SWIRE 5.8 micron peakers more massive than 1.6x1e11 Msun provide 30-50% of the total stellar mass density in galaxies at z=2-3.

X-ray Bright Optically Normal Galaxies (XBONGs) constitute a small but not negligible fraction of hard X-ray selected sources in recent Chandra and XMM-Newton surveys. Even though several possibilities were proposed to explain why a relatively luminous hard X-ray source does not leave any significant signature of its presence in terms of optical emission lines, the nature of XBONGs is still subject of debate. We aim to a better understanding of their nature by means of a multiwavelength and morphological analysis of a small sample of these sources. Good-quality photometric near-infrared data (ISAAC/VLT) of four low-redshift (z=0.1-0.3) XBONGs, selected from the HELLAS2XMM survey, have been used to search for the presence of the putative nucleus, applying the surface-brightness decomposition technique through the least-squares fitting program GALFIT. The surface brightness decomposition allows us to reveal a nuclear point-like source, likely to be responsible of the X-ray emission, in two out of the four sources. The results indicate that moderate amounts of gas and dust, covering a large solid angle (possibly 4pi) at the nuclear source, combined with the low nuclear activity, may explain the lack of optical emission lines. The third XBONG is associated with an X-ray extended source and no nuclear excess is detected in the near infrared at the limits of our observations. The last source is associated to a close (d< 1 arcsec) double system and the fitting procedure cannot achieve a firm conclusion.

Context: The counts of galaxies and AGN in the mid infra-red (MIR) bands are important instruments for studying their cosmological evolution. However, the classic spectral line ratios techniques can become misleading when trying to properly separate AGN from starbursts or even from apparently normal galaxies. Aims: We use X-ray band observations to discriminate AGN activity in previously classified MIR-selected starburst galaxies and to derive updated AGN1 and (Compton thin) AGN2 counts at 15 um. Methods: XMM observations of the ELAIS-S1 15um sample down to flux limits ~2x10^-15 erg cm^-2 s^-1 (2-10 keV band) were used. We classified as AGN all those MIR sources with a unabsorbed 2-10 keV X-ray luminosity higher that ~10^42 erg/s. Results: We find that at least about 13(+/-6) per cent of the previously classified starburst galaxies harbor an AGN. According to these figures, we provide an updated estimate of the counts of AGN1 and (Compton thin) AGN2 at 15 um. It turns out that at least 24% of the extragalactic sources brighter than 0.6 my at 15 um are AGN (~13% contribution to the extragalactic background produced at fluxes brighter than 0.6 mJy).

Aims: We aim at estimating the spectral energy distributions (SEDs) and the physical parameters related to the black holes harbored in eight high X-ray-to-optical (F_X/F_R>10) obscured quasars at z>0.9 selected in the 2--10 keV band from the HELLAS2XMM survey. Methods: We use IRAC and MIPS 24 micron observations, along with optical and Ks-band photometry, to obtain the SEDs of the sources. The observed SEDs are modeled using a combination of an elliptical template and torus emission (using the phenomenological templates of Silva et al. 2004) for six sources associated with passive galaxies; for two point-like sources, the empirical SEDs of red quasars are adopted. The bolometric luminosities and the M_BH-L_K relation are used to provide an estimate of the masses and Eddington ratios of the black holes residing in these AGN. Results: All of our sources are detected in the IRAC and MIPS (at 24 micron) bands. The SED modeling described above is in good agreement with the observed near- and mid-infrared data. The derived bolometric luminosities are in the range ~10^45-10^47 erg s^-1, and the median 2--10 keV bolometric correction is ~25, consistent with the widely adopted value derived by Elvis et al. (1994). For the objects with elliptical-like profiles in the K_s band, we derive high stellar masses (0.8-6.2)X10^11 Mo, black hole masses in the range (0.2-2.5)X10^9 Mo, and Eddington ratios L/L_Edd<0.1, suggesting a low-accretion phase.

(Abridged) We present the results of the photometric and spectroscopic identification of 110 hard X-ray selected sources (6e-15<F(2-10 keV)<4e-14 cgs) from 5 additional XMM-Newton fields, nearly doubling the original HELLAS2XMM sample. We spectroscopically identified 59 new sources, bringing the completeness of the full HELLAS2XMM sample to ~70%. We found optical counterparts for 214 out of the 232 X-ray sources of the full sample down to R~25. We measure the flux and luminosity of the [OIII](5007) emission line for 59 such sources. Results. We use the full HELLAS2XMM and the CDF samples to estimate the LogN-LogS and the obscured QSO surface densities down to F(2-10 keV)~1e-14 and 1e-15 cgs. At these flux limits the fraction of X-ray selected obscured and unobscured QSO turns out to be similar. Since X-ray selection misses most Compton thick AGN, the number of obscured QSO may well outnumber that of unobscured QSOs. We find that hard X-ray selected AGNs with detected [OIII] emission span a wide range of L(2-10 keV)/L[OIII] with a median higher than that of optically selected AGNs, suggesting that optically selected samples are at least partly incomplete, and/or [OIII] emission is not a perfect isotropic indicator of the nuclear power. Differently from narrow-line AGNs in the same z interval, the 7 XBONGs candidates in the sample have L(2-10 keV)/L[OIII]> 1000, while their Lx, L(R) and Nh are similar. This suggests that while the central engine of narrow-line AGNs and XBONGs looks similar, the narrow-line region in XBONGs could be strongly inhibited or obscured.

We report results from a Spitzer GO-1 program of IRS spectroscopy of a large sample of Luminous Infrared Galaxies and quasars selected from the European Large Area ISO Survey (ELAIS) which have a wide multiwavelength coverage, including ISOCAM, ISOPHOT, IRAC and MIPS (from SWIRE), and optical photometry. We present the sample selection and results from the IRS spectroscopy.

We have surveyed with XMM-Newton the central ~0.6 deg2 region of the ELAIS-S1 field down to flux limits of ~5.5X10-16 cgs (0.5-2 keV, S band), ~2X10-15 cgs (2-10 keV, H band), and ~4X10-15 cgs (5-10 keV, HH band). We detect a total of 478 sources, 395 and 205 of which detected in the S and H bands respectively. We identified 7 clearly extended sources and estimated their redshift through X-ray spectral fits with thermal models. In four cases the redshift is consistent with z=0.4. We have computed the angular correlation function of the sources in the S and H bands, finding best fit correlation angles theta_0=5.2+/-3.8 arcsec and theta_0=12.8+/-7.8 arcsec respectively. A rough estimate of the present-day correlation length r_0 can be obtained inverting the Limber equation and assuming an appropriate redshift distribution dN/dz. The results range between 12.8 and 9.8 h-1 Mpc in the S band and between 17.9 and 13.4 h-1 Mpc in the H band, with 30-40% statistical errors, assuming either smooth redshift distributions or redshift distributions with spikes accounting for the presence of a structure at z=0.4. The relative density of the S band sources is higher near the clusters and groups at z~0.4 and extends toward East and toward South/West. This suggests that the structure is complex, with a size comparable to the full XMM-Newton field. Conversely, the highest relative source densities of the H band sources are located in the central-west region of the field.

The ESO-Spitzer extragalactic Imaging Survey (ESIS) is the optical follow up of the Spitzer Wide-Area InfraRed Extragalactic (SWIRE) survey in the ELAIS-S1 area. This paper presents B, V, R Wide Field Imager observations of the first 1.5 square degree of the ESIS survey. Data reduction is described including astrometric calibration, illumination and color corrections, completeness and photometric accuracy estimates. Number counts and color distributions are compared to literature observational and theoretical data, including non-evolutionary, PLE, evolutionary and semi-analytic Lambda-CDM galaxy models, as well as Milky Way stellar predictions. ESIS data are in good agreement with previous works and are best reproduced by evolutionary and hierarchical Lambda-CDM scenarios. The ELAIS-S1 area benefits from extensive follow-up from X-ray to radio frequencies: some potential uses of the multi-wavelength observations are illustrated. Optical-Spitzer color-color plots promise to be very powerful tools to disentangle different classes of sources (e.g. AGNs, starbursts, quiescent galaxies). Ultraviolet GALEX data are matched to optical and Spitzer samples, leading to a discussion of galaxy properties in the UV-to-24 microns color space. The spectral energy distribution of a few objects, from the X-rays to the far-IR are presented as examples of the multi-wavelength study of galaxy emission components in different spectral domains.

We study the evolution of the luminosity function (LF) of type-1 and type-2 AGN in the mid-infrared, and derive their contribution to the Cosmic InfraRed Background (CIRB) and the expected deep source counts to be observed by Spitzer at 24 micron. The sample of type-1 and type-2 AGN was selected at 15 micron (ISO) and 12 micron (IRAS), and classified on the basis of their optical spectra. Local templates of type-1 and type-2 AGN have been used to derive the intrinsic 15 micron luminosities. We adopted an evolving smooth two-power law shape of the LF, whose parameters have been derived using an un-binned maximum likelihood method. We find that the LF of type-1 AGN is compatible with a pure luminosity evolution (L(z)=L(0)(1+z)^k_L) model where k_L~2.9. A small flattening of the faint slope of the LF with increasing redshift is favoured by the data. A similar evolutionary scenario is found for the type-2 population with a rate k_L ranging from ~1.8 to 2.6, depending significantly on the adopted mid-infrared spectral energy distribution. Also for type-2 AGN a flattening of the LF with increasing redshift is suggested by the data, possibly caused by the loss of a fraction of type-2 AGN hidden within the optically classified starburst and normal galaxies. The type-1 AGN contribution to the CIRB at 15 micron is (4.2-12.1) x 10e-11 W m^-2 sr^-1, while the type-2 AGN contribution is (5.5-11.0) x 10e-11 W m^-2 sr^-1. We expect that Spitzer will observe, down to a flux limit of S_24 = 0.01 mJy, a density of ~1200 deg^-2 type-1 and ~1000 deg^-2 type-2 optically classified AGN. The derived total contribution of the AGN galaxies to the CIRB (4-10%) and Spitzer counts should be considered as lower limits, because of a possible loss of type-2 sources caused by the optical classification.

We present VLT near-IR spectroscopic observations of three X-ray sources characterized by extremely high X-ray-to-optical ratios (X/O>40), extremely red colors (6.3<R-K<7.4, i.e. EROs) and bright infrared magnitudes (17.6<K<18.3). These objects are very faint in the optical, making their spectroscopic identification extremely challenging. Instead, our near-IR spectroscopic observations have been successful in identifying the redshift of two of them (z=2.08 and z=1.35), and tentatively even of the third one (z=2.13). When combined with the X-ray properties, our results clearly indicate that all these objects host obscured QSOs (4e44 < L(2-10keV) < 1.5e45 erg/s, 2e22 < N_H < 4e23 cm-2) at high redshift. The only object with unresolved morphology in the K band shows broad Halpha emission, but not broad Hbeta, implying a type 1.9 AGN classification. The other two objects are resolved and dominated by the host galaxy light in the K band, and appear relatively quiescent: one of them has a LINER-like emission line spectrum and the other presents only a single, weak emission line which we tentatively identify with Halpha. The galaxy luminosities for the latter two objects are an order of magnitude brighter than typical local L* galaxies and the derived stellar masses are well in excess of 10^11 Msun. For these objects we estimate black hole masses higher than 10^9 Msun and we infer that they are radiating at Eddington ratios L/L_Edd < 0.1. We discuss the implications of these findings for the coevolution of galaxies and black hole growth. Our results provide further support that X-ray sources with high X/O ratios and very red colors tend to host obscured QSO in very massive galaxies at high redshift.

We have determined the cosmological evolution of the density of AGN and of their Nh distribution as a function of the un-absorbed 2-10 keV luminosity up to redshift 4. We used the HELLAS2XMM sample combined with other published catalogs, yielding a total of 508 AGN. Our best fit is obtained with a luminosity-dependent density evolution (LDDE) model where low luminosity (Lx~10^43 erg s^-1) AGN peak at z~0.7, while high luminosity AGN (Lx>10^45 erg s^-1) peak at z~2.0. A pure luminosity evolution model (PLE) can instead be rejected. There is evidence that the fraction of absorbed (Nh>10^22 cm^-2) AGN decreases with the intrinsic X-ray luminosity, and increases with the redshift. Our best fit solution provides a good fit to the observed counts, the cosmic X-ray background, and to the observed fraction of absorbed AGN as a function of the flux in the 10^-15<S(2-10)<10^-10 erg s^-1 cm^-2 range. We find that the absorbed, high luminosity (Lx>10^44 erg s^-1) AGN have a density of 267 deg^-2 at fluxes S(2-10)>10^-15 erg s^-1 cm^-2. Using these results, we estimate a density of supermassive black holes in the local Universe of 3.2 (h^2_70) x 10^5 M_sol Mpc^-3, which is consistent with the recent measurements of the black hole mass function in the local galaxies.

We investigate the clustering properties of galaxies in the recently completed ELAIS-S1 redshift survey through their spatial two point autocorrelation function. We used a sub-sample of the ELAIS-S1 catalog covering approximately 4 deg^2 and consisting of 148 objects selected at 15 micron with a flux >0.5 mJy and redshift z<0.5. We detected a positive signal in the correlation function that, in the range of separations 1-10 h mpc is well approximated by a power law with a slope gamma = 1.4 and a correlation length s =5.4 h mpc. This result is in good agreement with the redshift-space correlation function measured in more local samples of mid infrared selected galaxies like the IRAS PSC$z$ redshift survey. This suggests a lack of significant clustering evolution of infrared selected objects out to $z=0.5$ that is further confirmed by the consistency found between the correlation functions measured in a local (z<0.2) and a distant (0.2<z<0.5) subsample of ELAIS-S1 galaxies. We also confirm that optically selected galaxies in the local redshift surveys, especially those of the SDSS sample, are significantly more clustered than infrared objects.

We present the Final Analysis of the European Large Area ISO Survey (ELAIS) 15 $\mu$m observations, carried out with the ISOCAM instrument on board the Infrared Space Observatory (ISO). The data reduction method, known as LARI method, is based on a mathematical model of the detector's behaviour and was specifically designed for the detection of faint sources in ISO-CAM/PHOT data. The method is fully interactive and leads to very reliable and complete source lists. The resulting catalogue includes 1923 sources detected with $S/N > 5$ in the {0.5 -- 100 mJy} flux range and over an area of 10.85 {deg$^2$} split into four fields, making it the largest non-serendipitous extragalactic source catalogue obtained to date from ISO data. This paper presents the concepts underlying the data reduction method together with its latest enhancements. The data reduction process, the production and basic properties of the resulting catalogue are discussed. The catalogue quality is assessed by means of detailed simulations, optical identifications and comparison with previous analyses.

We present the 90 micron luminosity function of the Final Analysis of the European Large Area ISO Survey (ELAIS), extending the sample size of our previous analysis (paper IV) by about a factor of 4. Our sample extends to z=1.1, around 50 times the comoving volume of paper IV, and 10^{7.7} < h^{-2}L/Lsun < 10^{12.5}. From our optical spectroscopy campaigns of the northern ELAIS 90 mircon survey (7.4 deg^2 in total, to S(90um)>70mJy), we obtained redshifts for 61% of the sample (151 redshifts) to B<21 identified at 7 microns, 15 microns, 20cm or with bright (B<18.5) optical identifications. The selection function is well-defined, permitting the construction of the 90 micron luminosity function of the Final Analysis catalogue in the ELAIS northern fields, which is in excellent agreement with our Preliminary Analysis luminosity function in the ELAIS S1 field from paper IV. The luminosity function is also in good agreement with the IRAS-based prediction of Serjeant & Harrison (2004), which if correct requires luminosity evolution of (1+z)^{3.4 +/- 1.0} for consistency with the source counts. This implies an evolution in comoving volume averaged star formation rate at z<~1 consistent with that derived from rest-frame optical and ultraviolet surveys.

The spectroscopic analysis of 117 serendipitous sources in the HELLAS2XMM 1df (1 degree field) survey is described. Of these, 106 sources, of which 86% have a spectroscopic redshift, are used to evaluate the fraction of X-ray absorbed (log NH>22) Active Galactic Nuclei (AGN) in the 2--10 keV flux range 0.8-20E-14 erg/cm2/s. This fraction turns out lower than what is predicted by two well known Cosmic X-Ray Background synthesis models, and the discrepancy is significant at the 99.999% level. This result consolidates the findings recently obtained by other authors. In the flux interval explored, the data are consistent with an intrinsic distribution of the absorbing columns (flat per decade above logNH>21) independent of luminosity and redshift, together with an AGN luminosity function evolving purely in luminosity. It is shown that, on the other hand, extrapolation to lower fluxes fails to reproduce the results inferred from the Chandra Deep Field North survey. It is found that about 40% of the high luminosity sources in our sample have best fit logNH>22, and the surface density of these X-ray obscured QSOs can then be estimated at about 48 per square degree, at the flux limit of ~1E-14} erg/cm2/s of the HELLAS2XMM 1df survey. As a side issue, 5 or 6 out of 60 sources, that is about 10%, identified with broad line AGN, turn out to be affected by logNH>22 absorption.

We present the final band-merged ELAIS catalogue at 6.7, 15, 90, and 175 $\mu$m, and the associated data at u,g,r,i,z,J,H,K, and 20cm. The origin of the survey, infrared and radio observations, data-reduction and optical identifications are briefly reviewed, and a summary of the area covered, and completeness limit for each infrared band is given. A detailed discussion of the band-merging and optical association strategy is given. The total catalogues consists of 2860 sources. For extragalactic sources observed in 3 or more infrared bands, colour-colour diagrams are presented and discussed in terms of the contributing infrared populations. Spectral energy distributions are shown for selected sources and compared with cirrus, M82 and Arp220 starburst, and AGN dust torus models. Spectroscopic redshifts are tabulated, where available. For the N1 and N2 areas, the INT ugriz Wide Field Survey, permits photometric redshifts to be estimated for galaxies and quasars. These agree well with the spectroscopic redshifts, within the uncertainty of the photometric method. The redshift distribution is given for selected ELAIS bands and colour-redshift diagrams are discussed. There is a high proportion of ultraluminous infrared galaxies in the ELAIS Catalogue ($> 10 %$ of 15 $\mu$m sources), many with Arp220-like colours. 10 hyperluminous infrared galaxies and 10 EROs are found in the survey. The large numbers of ultraluminous galaxies imply very strong evolution in the star-formation rate between z = 0 and 1.

We present a multi-wavelength catalog (15 um, R, K-band, 1.4 GHz flux) plus spectroscopic identifications for 406 15 um sources detected in the ELAIS region S1, over the flux density range 0.5<S<150 mJy. 332 (~82%) sources are optically identified down to R~23.0. Spectra or bona fide stellar identifications are obtained for 290 objects. The areal coverage, MIR and optical completeness of the sample are discussed in order to allow statistical and evolutionary analyses. Two main spectroscopic classes have been found to dominate the MIR extragalactic population: z<0.5 star-forming galaxies (from absorbed to extreme starbursts: nuL_nu~10^8-10^11 L_odot), which account for \~75% of the sources, and AGN (both type 1 and 2), which account for ~25% of the sources. About 20% of the extragalactic sources are dust-enshrouded starburst galaxies [e(a) spectra], and all the starburst galaxies appear more dust extincted in the optical than nearby normal galaxies. We also identified 91 stellar objects (~22% of the MIR sources). The counts for starburst galaxies and AGN down to 0.6 mJy have been derived. A general trend is found in the optical-MIR SED of the galaxies, where the MIR-luminous objects have larger MIR to optical luminosity ratios.

We present the first determination of the 15micron luminosity functio galaxies from the European Large Area ISO survey (ELAIS) southern fields. We have adopted a new criterion to separate the quiescent, non-evolving and the starburst, evolving populations based on the ratio of mid-infrared to optical luminosities. Strong evolution is suggested by our data for the starburst galaxy population, while normal spiral galaxies are consistent with no evolution. The starburst population must evolve both in luminosity and in density with rates of the order L(z) (1+z)^(3.5) and rho(z) (1+z)^3.8 up to z~1. The evolutionary parameters of our model have been tested by comparing the model predictions with other observables, like source counts at all flux density levels (from 0.1 to 300 mJy) and redshift distributions and luminosity functions at high-z (0.7 < z < 1.0 from HDF-N data). The agreement between our model predictions and the observed data is remarkably good. We use our data to estimate the star-formation density of the Universe up to z=0.4 and we use the luminosity function model to predict the trend of the star-formation history up to z=1.

We present an analysis of the two-point angular correlation function of the ELAIS S1 survey. The survey covers 4 deg$^2$ and contains 462 sources detected at 15$\mu$m to a 5$\sigma$ flux limit of 0.45 mJy. Using the 329 extragalactic sources not repeated in different observations, we detect a significant clustering signal; the resulting angular correlation function can be fitted by a exponential law $w(\theta) = A \theta^{1-\gamma}$ with $A = 0.014\pm0.005$ and $\gamma = 2.04\pm0.18$. Assuming a redshift distribution of the objects, we invert Limber's equation and deduce a spatial correlation length $r_0=4.3^{+0.4}_{-0.7} h^{-1}$ Mpc. This is smaller than that obtained from optical surveys but it is in agreement with results from IRAS. This extends to higher redshift the observational evidence that infrared selected surveys show smaller correlation lengths (i.e. reduced clustering amplitudes) than optical surveys.

We present the results of deep near-infrared observations of 11 hard X-ray selected sources in the Hellas2XMM survey, with faint optical magnitude (R>24) and high X-ray-to-optical flux ratio. All but one of the sources have been detected in the Ks band, with bright counterparts (Ks<19) and red colors (R-K>5), and therefore belong to the ERO population. A detailed analysis of the surface brightness profiles allows us to classify all of the NIR counterparts. There are 2 point-like objects, 7 elliptical galaxies and one source with a disky profile. None of the extended sources shows any evidence for the presence of a central unresolved object tracing the putative X-ray emitting AGN. Using both the R-K colors and the morphological information, we have estimated for all the sources a ``minimum photometric redshift'', ranging between 0.8 and 2.4; the elliptical hosts have zmin=0.9-1.4. We computed the X-ray properties using these redshifts: most of the sources have NH>10^{22}, with unabsorbed X-ray luminosities up to 10^{45}erg s^{-1}. These objects therefore belong to the population of obscured (Type II) quasars and, from a statistical point of view, they turn out to be a non-negligible fraction (~10%) of the most luminous AGN. Selecting the high X/O sources for a follow-up study in the NIR is therefore a powerful technique aimed at studying at high redshift the hosts of Type II AGN. Overall, our results seem to indicate that the hosts are mostly elliptical galaxies at z~1, and that these near-IR bright objects would be among the most massive spheroids at these epochs.

We present results from the photometric and spectroscopic identification of 122 X-ray sources discovered by XMM-Newton in the 2-10 keV band (the HELLAS2XMM 1dF sample). One of the most interesting results (which is found also in deeper sourveys) is that ~20% of the sources have an X-ray to optical flux ratio (X/O) ten times or more higher than that of optically selected AGN. Unlike the faint sources found in the ultra-deep Chandra and XMM-Newton surveys, which reach X-ray (and optical) fluxes >~10 times lower than in the HELLAS2XMM sample, many of the extreme X/O sources in our sample have R~25 and are therefore accessible to optical spectroscopy. We report the identification of 13 sources with extreme X/O values. While four of these sources are broad line QSO, eight of them are narrow line QSO, seemingly the extension to very high luminosity of the type 2 Seyfert galaxies.

We present results from the photometric and spectroscopic identification of 122 X-ray sources recently discovered by XMM-Newton in the 2-10 keV band (the HELLAS2XMM 1dF sample). Their flux cover the range 8E-15-4E-13 cgs and the total area surveyed is 0.9 deg2. About 20% of the hard X-ray selected sources have an X-ray to optical flux ratio (X/O) ten times or more higher than that of optically selected AGN. Unlike the faint sources found in the ultra-deep Chandra and XMM-Newton surveys, which reach X-ray (and optical) fluxes more than one order of magnitude lower than the HELLAS2XMM survey sources, many of the extreme X/O sources in our sample have R<=25 and are therefore accessible to optical spectroscopy. We report the identification of 13 sources with X/O>10: 8 are narrow line QSO (i.e. QSO2), four are broad line QSO. We use a combined sample of 317 hard X-ray selected sources (HELLAS2XMM 1dF, CDFN 1Msec, SSA13 and Lockman Hole flux limited samples), 221 with measured z, to evaluate the cosmological evolution of the hard X-ray source's number and luminosity densities. Looking backward in time, the low luminosity sources (logL(2-10keV) = 43-44 erg/s) increase in number at a rate different than the high luminosity sources (logL(2-10keV)>44.5 erg/s), reaching a maximum around z=1 and then levelling off beyond z=2. This translates into an accretion driven luminosity density which is dominated by sources with logL(2-10keV) < 44.5 erg/s up to at least z=1, while the contribution of the same sources and of those with logL(2-10keV)>44.5 erg/s appear to be comparable between z=2 and 4.

We present extensive optical, radio and infrared follow-up observations of a sample of 35 hard X-ray (2-10 keV) selected sources discovered serendipitously in the PV XMM-Newton observation of the radio-loud quasar PKS 0312-77 field, for which also an archival Chandra observation is available. The observations have been carried out as part of the HELLAS2XMM survey, a program aimed to understand the nature of the sources responsible for the bulk of the hard X-ray Background (XRB). The identification of the optical counterparts greatly benefits from the positional accuracy obtained from Chandra and radio observations. As a consequence, the spectroscopic completeness of the present sample (80%) is limited only by the faintness of the optical counterparts. The multiwavelength coverage of our survey allows us to unveil a large spread in the overall properties of hard X-ray selected sources. At low redshift (z<1), the source breakdown includes Broad Line AGN, Narrow Emission-Line Galaxies, and optically ``normal'' galaxies. All the ten sources at z>1 are spectroscopically classified as Broad Line AGNs. A few of them show significant intrinsic X-ray absorption (N(H)>10^22 cm^-2), further supporting previous evidence of a decoupling between optical and X-ray properties at high luminosities and redshifts. Finally, a non negligible fraction (15%) of the hard X-ray sources are not detected down to the limiting magnitude of the optical images. The corresponding high X-ray to optical flux ratio, X-ray and optical-infrared colors strongly suggest that they are high redshift, obscured AGN.

The radio counterparts to the 15-micron sources in the European Large Area ISO Survey southern fields are identified in 1.4-GHz maps down to about 80 microJy. The radio-MIR correlation is investigated and derived for the first time at these flux densities for a sample of this size. Our results show that radio and MIR luminosities correlate almost as well as radio and FIR, at least up to z ~ 0.6. Using the derived relation and its spread together with the observed 15-micron counts, we have estimated the expected contribution of the 15-micron extragalactic populations to the radio source counts and the role of MIR starburst galaxies in the well known 1.4-GHz source excess observed at sub-mJy levels. Our analysis demonstrates that IR emitting starburst galaxies do not contribute significantly to the 1.4-GHz counts for strong sources, but start to become a significant fraction of the radio source population at flux densities < 0.5 - 0.8 mJy. They are expected to be responsible for more than 60% of the observed radio counts at < 0.05 mJy. These results are in agreement with the existing results on optical identifications of faint radio sources.

We present results of a complete radio follow-up obtained with the VLA and ATCA radio telescopes down to a 6 cm flux limit of about 0.3 mJy of all the 147 X-ray sources detected in the BeppoSAX HELLAS survey. We found 53 X-ray/radio likely associations, corresponding to about one third of the X-ray sample. Using the two point spectral index alpha_ro=0.35 we divided all the HELLAS X-ray sources in radio quiet and radio loud. We have 26 sources classified as radio-loud objects, corresponding to about 18% of the HELLAS sample. In agreement with previous results, the identified radio-loud sources are associated mainly with Type 1 AGNs with L(5-10 keV) > 10^44 erg/s, while all the identified Type 2 AGNs and Emission Line Galaxies are radio quiet objects with L(5-10 keV) < 10^44 erg/s. The analysis of the radio spectral index suggests that Type 1 AGNs have a mean radio spectral index flatter than Type 2 AGNs and Emission Line Galaxies.

We present the ISO far-IR photometry of a complete sub-sample of optically selected bright quasars belonging to two complete surveys selected through multicolour (U,B,V,R,I) techniques. The ISOPHOT camera on board of the ISO Satellite was used to target these quasars at wavelengths of 7.3, 11.5, 60, 100 and 160 micron. Almost two thirds of the objects were detected at least in one ISOPHOT band. The detection rate is independent of the source redshift, very likely due to the negative K-correction of the far-IR thermal emission. More than a half of the optically selected QSOs show significant emission between 4 and 100 micron in the quasar rest-frame. These fluxes have a very likely thermal origin, although in a few objects an additional contribution from a non-thermal component is plausible in the long wavelength bands. In a colour-colour diagram these objects span a wide range of properties from AGN-dominated to ULIRG-like. The far-IR composite spectrum of the quasar population presents a broad far-IR bump between 10 and 30 micron and a sharp drop at wavelengths greater than 100 micron in the quasar restframe. The amount of energy emitted in the far-IR, is on average a few times larger than that emitted in the blue and the ratio L(FIR)/L(B) increases with the bolometric luminosity. Objects with fainter blue magnitudes have larger ratios between the far-IR (wavelengths > 60 micron) fluxes and the blue band flux, which is attributed to extinction by dust around the central source. No relation between the blue absolute magnitude and the dust colour temperature is seen, suggesting that the dominant source of FIR energy could be linked to a concurrent starburst rather than to gravitational energy produced by the central engine.

We present the 15 micron luminosity function of type 1 AGN (QSO + Seyfert 1). Our sample of 21 high-redshift sources is selected from the Preliminary Analysis catalogue in the S1 field of the European Large Area ISO Survey (ELAIS). To study the cosmological evolution of the AGN1 luminosity function, our sample has been combined with a local sample of 41 sources observed by IRAS. We find that the luminosity function of AGN1 at 15 micron is fairly well represented by a double-power-law-function. There is evidence for significant cosmological evolution consistent with a pure luminosity evolution (PLE) model L(z) (1+z)^k, with k=3.0-3.3. The value of k depends on the existence or not of an evolution cut-off at redshift ~2, and on the adopted cosmology. From the luminosity function and its evolution we estimate a contribution of AGN1 to the Cosmic Infrared Background (CIRB) of nuI_nu ~ 6 x 10^{-11}W m^{-2} sr^{-1} at 15 micron. This corresponds to ~2-3% of the total observed CIRB at this wavelength. Under the usual assumptions of unified models for AGN, the expected contribution of the whole AGN population to the CIRB at 15 micron is 10-15%.

Multiwavelength observations of the hard X-ray selected sources by Chandra and XMM-Newton surveys have significantly improved our knowledge of the objects responsible of the hard X-ray background. A surprising finding is the discovery of a population of optically dull, X-ray bright galaxies emerging at 2-10 keV fluxes of the order of 10^{-14} erg cm^{-2} s^{-1}. We present preliminary results of multiwavelength observations of a few objects serendipitously discovered in the field of XMM-Newton and Chandra observations and intensively studied at longer wavelengths.

We present optical spectroscopic identifications of hard X-ray (5-10 keV) selected sources belonging to the HELLAS sample obtained with BeppoSAX down to a 5-10 keV flux limit of f_{5-10 keV}~3x10^-14 erg cm^-2 s^-1. The sample consists of 118 sources. 25 sources have been identified trough correlations with catalogues of known sources. 49 have been searched for spectroscopic identification at the telescope. 13 fields resulted empty down to R=21. 37 sources have been identified with type 1 AGN and 9 with type 2 AGN. The remaining are: 5 narrow emission line galaxies, 6 Clusters, 2 BL Lac, 1 Radio Galaxy and 1 Star. Combining these objects with other hard X-ray selected AGNs from ASCA and HEAO1, we find that the local luminosity function of type 1 AGN in the 2-10 keV band is fairly well represented by a double-power-law-function. There is evidence for significant cosmological evolution according to a pure luminosity evolution (PLE) model L_X(z)~(1+z)^k, with k=2.12 and k=2.22 in a (Omega_m,Omega_lambda)=(1.0,0.0) and in a (Omega_m,Omega_lambda)=(0.3,0.7) cosmology, respectively. The data show an excess of faint high redshift type 1 AGN which is well modeled by a luminosity dependent density evolution (LDDE), similarly to what observed in the soft X-rays. However, in both cosmologies, the statistic is not significant enough to distinguish between the PLE and LDDE models. The fitted models imply a contribution of AGN1 to the 2-10 keV X-ray background from 35% up to 60%.

Recent X-ray surveys have clearly demonstrated that a population of optically dull, X-ray bright galaxies is emerging at 2-10 keV fluxes of the order of 10^{-14} erg cm^{-2} s^{-1}. Although they might constitute an important fraction of the sources responsible for the hard X-ray background, their nature is still unknown. With the aim to better understand the physical mechanisms responsible for the observed properties, we have started an extensive program of multiwavelength follow-up observations of hard X-ray, optically quiet galaxies discovered with XMM-Newton. Here we report the results of what can be considered the first example of this class of objects: CXOUJ031238.9-765134, originally discovered by Chandra, and optically identified by Fiore et al. (2000) with an apparently normal early-type galaxy at z=0.159, usually known as "FIORE P3". The analysis of the broad-band energy distribution suggests the presence of a heavily obscured active nucleus.

We present the estimate of the evolution of type 1 AGNs in the hard (2-10 keV) X-rays drawn from the HELLAS survey, and in the IR (15um) obtained from the ELAIS survey. We find that the local luminosity function (LF) of AGN1 in the 2-10 keV band is fairly well represented by a double-power-law-function. There is evidence for significant cosmological evolution according to a pure luminosity evolution model Lx(z)~(1+z)^k, with k=2.12 and k=2.19 in a (Omega_m,Omega_lambda)=(1.0,0.0) and in a (Omega_m,Omega_lambda)=(0.3,0.7) cosmology respectively. In a (Omega_m,Omega_lambda)=(1.0,0.0) Universe the data show an excess of faint high redshift type 1 AGN which is well modeled by a luminosity dependent density evolution, similarly to what observed in the soft X-rays. In the IR band, with a (Omega_m,Omega_lambda)=(1.0,0.0) cosmology, the evolution found is similar to what observed in other wavebands, the LF is a double-power-law-function with a bright slope 2.9 and a faint slope 1.1, following a luminosity evolution model L(z)~(1+z)^3.0.

We present preliminary results on the optical identifications of sources from the High Energy Large Area Survey performed with XMM-Newton (HELLAS2XMM). This survey covers about 3 square degrees of sky down to a 2-10 keV flux of 7 x 10^-15 erg cm-2 s-1 (Baldi & Molendi these proceedings, Baldi et al. 2001). The survey good sensitivity over a large area allows us to confirm and extend previous Chandra results about a population of X-ray luminous but optically dull galaxies emerging at 2-10 keV fluxes of about 10-14 erg cm-2 s-1. Although the statistics are still rather poor, it appears that at these X-ray fluxes the fraction of these galaxies is similar or even higher than that of narrow line AGN.

Recent X-ray surveys have clearly demonstrated that a population of optically dull, X-ray bright galaxies is emerging at X-ray fluxes of the order of 10-14 erg cm-2 s-1. The nature of these objects is still unknown. We present the results of multiwavelength observations of what can be considered the best studied example: the Chandra source CXOUJ031238.9-765134 optically identified by Fiore et al. (2000) with an apparently normal galaxy at z=0.159 and called ``Fiore P3''.

The BeppoSAX High Energy Large Area Survey (HELLAS) has surveyed several tens of square degrees of the sky in the 5--10 keV band down to a flux of about 5 10^-14 erg cm-2 s-1. The extrapolation of the HELLAS logN--logS towards fainter fluxes with an euclidean slope is consistent with the first XMM measurements, in the same energy band, which are a factor 20 more sensitive. The source counts in the hardest band so far surveyed by X-ray satellites are used to constrain XRB models. It is shown that in order to reproduce the 5--10 keV counts over the range of fluxes covered by BeppoSAX and XMM a large fraction of highly absorbed (logN_H = 23--24 cm-2), luminous (L_X > 10^44 erg s-1) AGN is needed. A sizeable number of more heavily obscured, Compton thick, objects cannot be ruled out but it is not required by the present data. The model predicts an absorption distribution consistent with that found from the hardness ratios analysis of the so far identified HELLAS sources. Interestingly enough, there is evidence of a decoupling between X-ray absorption and optical reddening indicators especially at high redshifts/luminosities where several broad line quasars show hardness ratios typical of absorbed power law models with logN_H=22--24 cm-2.

The BeppoSAX High Energy Large Area Survey (HELLAS) has surveyed about 85 deg^2 of sky in the 5-10 keV band down to a flux of 4-5 10^-14 erg cm-2 s-1. The source surface density of 16.9+-6.4 deg$^{-2}$ at the survey limit corresponds to a resolved fraction of the 5-10 keV X-ray background (XRB) of the order of 20-30 %. Hardness ratios analysis indicates that the spectra of a substantial fraction of the HELLAS sources (at least one third) are harder than a alpha_E = 0.6 power law. This hardness may be due to large absorbing columns. The hardness ratio analysis also indicates that many HELLAS sources may have a spectrum more complex than a single absorbed power law. A soft component, superimposed to a strongly cut-off power law, is likely to be present in several sources.

We present a comprehensive study of the soft X-ray properties of the BeppoSAX High-Energy Large Area Survey (HELLAS) sources. A large fraction (about 2/3) of the hard X-ray selected sources is detected by ROSAT. The soft X-ray colors for many of these objects, along with the 0.5-2 keV flux upper limits for those undetected in the ROSAT band, do imply the presence of absorption. The comparison with the ROSAT Deep Survey sources indicates that a larger fraction of absorbed objects among the HELLAS sources is present, in agreement with their hard X-ray selection and the predictions of the X-ray background synthesis models. Another striking result is the presence of a soft (additional) X-ray component in a significant fraction of absorbed objects.

Optical identification of hard X-ray selected BeppoSAX and Chandra sources indicates that a large fraction of the sources are ``intermediate'' AGN, i.e. type 1.8-1.9 AGN, broad-line quasars and even X-ray loud but optically dull (apparently normal) galaxies, all obscured in X-rays by columns of the order of logNH=22-23.5 cm-2. Because of this obscuring matter, these sources are more difficult to detect or select at other wavelengths, implying that a fraction of the accretion power in the Universe may have been missed so far.

We present BeppoSAX observations of the Southern S1 region in the European Large Area ISO Survey (ELAIS). These observations cover an area of ~1.7 sq. deg. and reach an on-axis (~0.7 sq. deg) 2-10 keV (HX) sensitivity of ~10E-13 cgs. This is the first HX analysis of an ISOCAM survey. We detect 9 sources with a signal to noise ratio SNR_HX>3, 4 additional sources with a 1.3 to 10 keV (T) SNR_T>3 and 2 additional sources which appear to be associated with QSOs with SNR_T>2.9. The number densities of the SNR_HX>3 sources are consistent with the ASCA and BeppoSAX logN-logS functions. Six BeppoSAX sources have reliable ISOCAM 15 micron counterparts within ~60 arcsec. All these ISOCAM sources have optical counterparts of R<20 mags. Five of these sources have been previously optically classified giving 4 QSOs and 1 BALQSO at z=2.2. The remaining unclassified source has X-ray and photometric properties consistent with that of a nearby Seyfert galaxy. One further HX source has a 2.6sigma ISOCAM counterpart associated with a galaxy at z=0.325. If this ISOCAM source is real, the HX/MIR properties suggest either an unusual QSO or a cD cluster galaxy. We have constructed MIR and HX spectral energy distributions to compute the expected HX/MIR ratios for these classes of objects up to z=3.2 and assess the HX/MIR survey depth. The BALQSO has an observed X-ray softness ratio and HX/MIR flux ratio similar to QSOs but different to those found for low redshift BALQSOs. This difference can be explained in terms of absorption, and suggests that high redshift BALQSOs should be comparatively easy to detect in the HX band, allowing their true fraction in the high redshift QSO population to be determined (ABRIDGED).

We present optical spectroscopic identification of sources identified in the BeppoSAX MECS fields of the High Energy LLarge Area Survey (HELLAS). In total 62 sources from a sample of 115 brighter than F_{5-10keV}>5x10^{-14} cgs have been identified. We find a density of 13-21 sources/deg^2 at F_{5-10keV}>5x10^{-14} cgs, which contribute to 20-30 % of the hard X-ray background. Evidences are found for type 1 AGN being more absorbed with increasing redshift or luminosity. The low redshift (z<0.2) ratio of type 2 to type 1 AGN is 3+/-1.5, in agreement with the unified models for AGN. The luminosity function of type 1 AGN in the 2-10 keV band is preliminary fitted by a two-power law function evolving according to a pure luminosity evolution model: L~(1+z)^{2.2}.

We present the luminosity function of 90um selected galaxies from the European Large Area ISO Survey (ELAIS), extending to z=0.3. Their luminosities are in the range 10^9 < h_65^-2 L/Lsun < 10^12, i.e. non-ultraluminous. From our sample of 37 reliably detected galaxies in the ELAIS S1 region from the Efstathiou et al. (2000) S_90 >= 100mJy database, we found optical, 15um or 1.4GHz identifications for 24 (65%). We have obtained 2dF and UK Schmidt FLAIR spectroscopy of 89% of IDs to rigid multivariate flux limits. We construct a luminosity function assuming (a) our spectroscopic subset is an unbiased sparse sample, and (b) there are no galaxies which would not be represented in our spectroscopic sample at {\it any} redshift. We argue that we can be confident of both assumptions. We find the luminosity function is well-described by the local 100um luminosity function of Rowan-Robinson, Helou & Walker (1987). {\it Assuming} this local normalisation, we derive luminosity evolution of (1+z)^{2.45\pm0.85} (95% confidence). We argue that star formation dominates the bolometric luminosities of these galaxies and we derive comoving star formation rates in broad agreement with the Flores et al. (1999) and Rowan-Robinson et al. (1997) mid-IR-based estimates.

We present observational evidences that dust in the circumnuclear region of AGNs has different properties than in the Galactic diffuse interstellar medium. By comparing the reddening of optical and infrared broad lines and the X-ray absorbing column density we find that the E(B-V)/N_H ratio is nearly always lower than Galactic by a factor ranging from ~3 up to ~100. Other observational results indicate that the Av/N_H ratio is significantly lower than Galactic in various classes of AGNs including intermediate type 1.8-1.9 Seyferts, hard X-ray selected and radio selected quasars, broad absorption line QSOs and grism selected QSOs. The lack of prominent absorption features at 9.7um (silicates) and at 2175A (carbon dip) in the spectra of Seyfert 2s and of reddened Seyfert 1s, respectively, add further evidence for dust in the circumnuclear region of AGNs being different from Galactic. These observational results indicate that the dust composition in the circumnuclear region of AGNs could be dominated by large grains, which make the extinction curve flatter, featureless and are responsible for the reduction of the E(B-V)/N_H and Av/N_H ratios. Regardless of the physical origin of these phenomena, the reduced dust absorption with respect to what expected from the gaseous column density should warn about a mismatch between the optical and the X-ray classification of the active galactic nuclei in terms of their obscuration.

The BeppoSAX 4.5-10 keV High Energy Large Area Survey has covered about 80 square degrees of sky down to a flux of F(5-10keV)~5E-14 cgs. Optical spectroscopic identification of about half of the sources in the sample (62) shows that many (~50%) are highly obscured AGN, in line with the predictions of AGN synthesis models for the hard X-ray background (XRB, see e.g. Comastri et al. 1995). The X-ray data, complemented by optical, near-IR and radio follow-up, indicate that the majority of these AGN are ``intermediate'' objects, i.e. type 1.8-1.9 AGN,`red' quasars, and even a few broad line, blue continuum quasars, obscured in X-rays by columns of the order of logNH=22.5-23.5 cm-2, but showing a wide dispersion in optical extinction. The optical and near-IR photometry of the obscured objects are dominated by galaxy starlight, indicating that a sizeable fraction of the accretion power in the Universe may actually have been missed in optical color surveys. This also implies that multicolor photometry techniques may be efficiently used to assess the redshift of the hard X-ray selected sources.

The energy density of the Cosmic X-ray background (XRB) peaks around 30 keV (see Figure 1), an energy not yet probed by focussing imaging instruments. The first hard X-ray telescope due to fly on a space mission will be that on board Constellation-X. The imaging capability, besides providing an improvement of several orders of magnitude in sensitivity over current passively collimated detectors, will permit for the first time to resolve a fraction of the XRB at this most crucial energy. Synthesis models of the XRB based on obscured AGN predict that at least 40% of the 10-40 keV XRB will be resolved by Constellation-X.

Obscured AGNs are thought to contribute a large fraction of the hard X-ray background (2-10 keV), and have also been proposed as the powerhouse of a fraction of the SCUBA sources which make most of the background at 850um, thus providing a link between the two spectral windows. We have tackled this issue by comparing data at 2-10 keV and at 850um for a sample of 34 sources at fluxes (or limiting fluxes) which resolve most of the background in the two bands. We present here new SCUBA observations, and new correlations between separate data sets retrieved from the literature. Similar correlations presented by others are added for completeness. None of the 11 hard X-ray (2-10 keV) sources has a counterpart at 850um, with the exception of a Chandra source in the SSA13 field, which is a candidate type 2, heavily absorbed QSO at high redshift. The ratios F(850um)/F(5keV) (mostly upper limits) of the X-ray sources are significantly lower than the value observed for the cosmic background. In particular, we obtain that 2-10 keV sources brighter than 10^-15 erg s^-1 cm^-2, which make at least 75% of the background in this band, contribute for less than 7% to the submillimeter background. Out of the 24 SCUBA sources, 23 are undetected by Chandra. The ratios F(850um)/F(5keV) (mostly lower limits) of these SCUBA sources indicate that most of them must be powered either by starburst activity, or by an AGN which is obscured by a column Nh > 10^25 cm^-2, with a reflection efficiency in the hard X rays significantly lower than 1% in most cases. However, AGNs of this type could not contribute significantly to the 2-10 keV background.

We present the first preliminary estimate of the evolution of type 1 AGNs at 15 um. A new sample of sources selected from the ELAIS-S1 survey has been used together with the RMS sample. The ELAIS-S1 survey covers an effective area of 2.2 deg^2 (total: 4 deg^2) down to 1.2 mJy at 15 um. Spectroscopic identifications have been carried out in the range 17.2<R<20.0 at the 3.6m/ESO, NTT/ESO and 2dF/AAT. The luminosity function is steeper at bright luminosities in comparison to the previous estimates of Rush, Malkan & Spinoglio (1993). The evolution is compatible with a Pure Luminosity Evolution model with L(z)=L(0)(1+z)^3.4.

We present some recent results from the European Large Area ISO Survey (ELAIS). This survey was the largest non-serendipitous ISO field survey. A preliminary reduction has recently been completed and catalogues of sources released to the community. Early results show strongly evolving source counts. A comprehensive identification programme is underway and a number of extremely luminous objects have already been discovered. This survey provides an exciting legacy from the ISO mission and (amongst many goals) will allow us to provide important constraints on the obscured star-formation history of the Universe.

We report optical spectroscopic identifications of 10 hard (2-10 keV) X-ray selected sources discovered by Chandra. The X-ray flux of the sources ranges between 1.5 and 25 x 10-14 cgs, the lower value being 3 times fainter than in previous BeppoSAX and ASCA surveys. Their R band magnitudes are in the range 12.8-22. Six of the Chandra sources are broad line quasars with redshifts between 0.42 and 1.19, while the optical identification of the remaining four is quite varied: two are X-ray obscured, emission line AGN at z=0.272 and z=0.683, one is a starburst galaxy at z=0.016 and one, most unusually, is an apparently normal galaxy at z=0.158. These findings confirm and extend down to fainter X-ray fluxes the BeppoSAX results, in providing samples with a wide range of X-ray and optical properties. The ratio between the soft X-ray and the optical luminosity of the z=0.158 galaxy is a factor at least 30 higher than that of normal galaxies, and similar to those of AGN. The high X-ray luminosity and the lack of optical emission lines suggest an AGN in which either continuum beaming dominates, or emission lines are obscured or not efficiently produced.

We present near-IR photometry and imaging observations of a small sample of sources identified in the BeppoSAX 5-10 keV survey (HELLAS) which resolves ~ 20-30% of the X-ray background at these energies. The near-IR data are combined with optical spectra and photometry. Only 40% of the sources in our sample have the blue, power law continuum typical of color-selected QSOs. The remaining 60% are dominated by a galactic component which, on the basis of the continuum colors and shape, have ages ranging from 10^9 to 10^10 years. The images show that the blue QSOs are pointlike at our angular resolution, while all the other sources are extended, consistent with their spectral appearance and low redshift. Since down to R=20 only about two thirds of the HELLAS sources have a counterpart, the preliminary HELLAS census comprises in roughly equal parts: i) blue QSOs (mostly at high redshifts); ii) optically dim, galaxy-dominated active nuclei (mostly at modest redshifts); and iii) empty fields (possibly highly absorbed QSOs at high redshifts).

The integrated emission of highly obscured AGN is expected to provide a major contribution to the X-ray energy density in the Universe: the X-ray background (XRB). The study of these objects is possible only at energies where the effects of absorption are less severe. For this reason we have carried out the BeppoSAX High Energy LLarge Area Survey in the hardest band (5-10 keV) accessible so far with imaging X-ray instruments. The source surface density at the survey limiting flux accounts for a significant fraction (20-30%) of the hard XRB. The X-ray data complemented by multiwavelength follow--up observations suggest that a large fraction of the hard sources are AGN and that X-ray absorption with column densities in the range 10^22-23.5 cm-2 is common among them. The great diversity in their optical-near-IR properties suggests that the optical appearance of obscured sources is a function of the X-ray luminosity. We briefly discuss the implications of these findings for the XRB models.

A synthesis model for the cosmic X-ray Background (XRB) is presented, which includes a proper treatment of Compton scattering in the absorbing matter for type 2 AGN. Evidence for a decrease of the relative importance of type 2 AGN at high redshift is found, which may be due either to a decrease of the relative number of obscured sources, or (more plausibly) to an increase of the fraction of Compton-thick absorbed sources. The XRB spectrum, soft X-rays and hard X-rays source counts can be simultaneously fitted only if the XRB normalization as derived from BeppoSAX/MECS measurements is adopted.

Spectra, position, magnitudes and colors are presented for 485 faint (B<20.5) emission line objects selected with the ultraviolet-excess (UVX) criterion on a area of 24.6 sq. deg in the South Galactic Pole. The objects were selected from the analysis of pixel-to-pixel stacking of COSMOS scans of UKST U, J and R plates. The candidates were observed with the Meudon-ESO Fiber Optics System (MEFOS) at the ESO 3.6m telescope. 429 type 1 AGNs have been identified (373 in the redshift range 0.3<z<2.2). This sample has allowed the measure of a difference on the QSO clustering evolution in comparison with that found for galaxies (La Franca et al 1998). The region is part of the ESO Imaging Survey (EIS) and of the 2dF QSO redshift survey.

We discuss the constraints on the AGN evolution from the cosmic X-ray background and source counts. A synthesis model to fit the X-ray background is presented. In the model, the spectrum of type 2 AGN has been modeled including Compton down--scattering within the absorbing material. Besides, we introduced a dependence on redshift of the relative number of obscured sources and found a decrease of the fraction of type 2 AGN at redshifts larger than 2.

X-ray spectra transmitted through matter which is optically thick to Compton scattering are computed by means of Monte Carlo simulations. Applications to the BeppoSAX data of the Seyfert 2 galaxy in Circinus, and to the spectral modeling of the Cosmic X-ray Background, are discussed.

Hard X-ray selection is the most efficient way to discriminate between accretion-powered sources, such as AGN, from sources dominated by starlight. Hard X-rays are also less affected than other bands by obscuration. We have then carried out the BeppoSAX High Energy Large Area Survey (HELLAS) in the largely unexplored 5-10 keV band, finding 180 sources in ~50 deg^2 of sky with flux >5E-14 erg cm-2 s-1. After correction for the non uniform sky coverage this corresponds to resolving about 30 % of the hard Cosmic X-ray Background (XRB). Here we report on a first optical spectroscopic identification campaign, finding 12 AGN out of 14 X-ray error-boxes studied. Seven AGN show evidence for obscuration in X-ray and optical bands, a fraction higher than in previous ROSAT or ASCA-ROSAT surveys (at a 95-99 % and 90 % confidence level respectively), thus supporting the scenario in which a significant fraction of the XRB is made by obscured AGN.

First results on a medium-deep X-ray survey in the "new" 5-10 keV band carried out with the MECS detectors onboard BeppoSAX are presented. The High Energy Llarge Area Survey (HELLAS) is aimed to directly explore a band where the energy density of the X-ray background is more than twice than that in the soft (0.5-2.0 keV) band. The optical identification follow-up of the first ten HELLAS hard X-ray sources indicate that Active Galactic Nuclei are the dominant population at 5-10 keV fluxes of the order of 10e-13 cgs. We discuss the implications of these findings for the AGN synthesis models for the XRB.

We have surveyed about 50 deg-2 of sky in the largely unexplored 5-10 keV band using the BeppoSAX MECS instrument, finding 180 sources. After correction for the non uniform sky coverage we find about 18 sources deg-2 with F(5-10keV) >= 5E-14 cgs, and resolve 30-40 % of the hard 5-10 keV Cosmic X-ray Background (XRB). Optical identification of a first small sample of sources show that most (11 out of 14) are AGN. Six of these show evidence of absorption/extinction in X-ray/optical, thus providing support to the scenario in which the hard XRB is largely made by obscured AGN (Setti & Woltjer 1989, Comastri et al. 1995).

Recent results have questioned the description of the QSO luminosity function in terms of a pure luminosity evolution and call for a luminosity dependent luminosity evolution. Measurements of the QSO clustering amplitude and evolution allow further distinguishing among the various physical scenarios proposed to interpret the QSO phenomenon. The general properties of the QSO population would arise naturally if quasars are short-lived events connected to a characteristic halo mass $\sim 5 \cdot 10^{12}$ M$_\odot$. This is the typical mass of groups of galaxies in which the interactions triggering the QSO activity preferentially take place.

The evolution of QSO clustering is investigated with a new sample of 388 QSOs with 0.3<z<2.2, B<20.5 and M_B<-23. Evidence is found for an increase of the clustering amplitude with increasing redshift. These measurements allow to further distinguish among the various physical scenarios proposed to interpret the QSO phenomenon. A single population model is inconsistent with the observations. The general properties of the QSO population would arise naturally if quasars are short-lived events connected to a characteristic halo mass ~5x10^{12} M_sun.

An ESO Key programme dedicated to an Homogeneous Bright QSO Survey (HBQS) has been completed. 327 QSOs (Mb<-23, 0.3<z<2.2) have been selected over 555 deg^2 with 15<B<18.75. For B<16.4 the QSO surface density turns out to be a factor 2.2 higher than what measured by the PG survey, corresponding to a surface density of 0.013+/-.006 deg^{-2}. If the Edinburgh QSO Survey is included, an overdensity of a factor 2.7 is observed, corresponding to a density of 0.016+/-0.005 deg^{-2}. In order to derive the QSO optical luminosity function (LF) we used Monte Carlo simulations that take into account of the selection criteria, photometric errors and QSO spectral slope distribution. The LF can be represented with a Pure Luminosity Evolution (L(z)\propto(1+z)^k) of a two power law both for q_0=0.5 and q_0=0.1. For q_0=0.5 k=3.26, slower than the previous Boyle's (1992) estimations of k=3.45. A flatter slope beta=-3.72 of the bright part of the LF is also required. The observed overdensity of bright QSOs is concentrated at z<0.6. It results that in the range 0.3<z<0.6 the luminosity function is flatter than observed at higher redshifts. In this redshift range, for Mb<-25, 32 QSOs are observed instead of 19 expected from our best-fit PLE model. This feature requires a luminosity dependent luminosity evolution in order to satisfactorily represent the data in the whole 0.3<z<2.2 interval.

We present BeppoSAX observations of the peculiar long period polar system V1309 Ori (RXJ0515.6+0105). The source was detected simultaneously at soft and, for the first time, at hard X-rays with the LECS and the MECS detectors. Both, the LECS and the MECS light curves are irregular with a bursting/flaring type behaviour indicating inhomogeneous accretion onto the white dwarf. This peculiar variability, together with an extreme high soft-to-hard X-ray luminosity ratio, indicates that in V1309 Ori accretion occurs predominantly in highly compressed chunks or ``blobs'' of matter. From coordinated ESO optical spectroscopy, we find indications that the magnetic field strength of the white dwarf is < 70 MG, not expected either from the 8 hr orbital period synchronism or from the strong soft-to-hard X-ray ratio suggesting alternative solutions for sustaining synchronism in this system.

The average optical-UV continuum shape of QSOs has been investigated using spectra of 62 QSOs having good relative photometric calibrations. The QSO spectra were extracted from two complete color selected samples in the magnitude intervals B ~ 16-20. The analysis was performed fitting power-law continua (f proportional to nu^(alpha)) in well defined rest-frame wavelength intervals after removing regions of the spectrum affected by strong emission lines or weak emission bumps. The average slope in the rest-frame optical-UV region 1200 - 5500 A shows a rapid change around the 3000 A emission bump with alpha=0.15 longward of it and alpha=-0.65 at shorter wavelengths. Although these average slopes have been obtained using spectra of QSOs with different luminosities and redshifts, there are no significant correlations of the average spectral index with these quantities. For a few QSOs in the sample we were able to measure the same softening of the spectral shape within the individual spectrum. These results have significant consequences on the estimate of the cosmological evolution of the optically selected QSOs as they affect, for instance, the k-corrections. New k-corrections in the B, V, R and Gr bands were computed. The derived average spectral shape in the optical-UV band puts interesting constraints on the expected emission mechanisms.

The evolution of quasar clustering is investigated with a new sample of 388 quasars with 0.3<z<=2.2, B<=20.5 and Mb<-23, selected over an area of 24.6 sq. deg. in the South Galactic Pole. Assuming a two-point correlation function of the form xi(r) = (r/r_o)^-1.8, we detect clustering with r_0=6.2 +/- 1.6 h^-1 comoving Mpc at an average redshift of z=1.3. We find a 2 sigma significant increase of the quasar clustering between z=0.95 and z=1.8, independent of the quasar absolute magnitude and inconsistent with recent evidence on the evolution of galaxy clustering. If other quasar samples are added (resulting in a total data-set of 737 quasars) the increase of the quasar clustering is still favoured although it becomes less significant. We find epsilon=-2.5. Evolutionary parameters epsilon>0.0 are excluded at a 0.3% probability level, to be compared with epsilon=0.8 found for galaxies. The observed clustering properties appear qualitatively consistent with a scenario of Omega=1 CDM in which a) the difference between the quasar and the galaxy clustering can be explained as a difference in the effective bias and redshift distributions, and b) the quasars, with a lifetime of t~10^8 yr, sparsely sample halos of mass greater than M_min~10^12-10^13 h^-1 M_sun. We discuss also the possibility that the observed change in the quasar clustering is due to an increase in the fraction of early-type galaxies as quasar hosts at high z.

An ESO Key programme dedicated to an Homogeneous Bright QSO Survey (HBQS) has been completed. 327 QSOs (Mb<-23, 0.3<z<2.2) have been selected over 555 deg^2 with 15<B<18.75. For B<16.4 the QSO surface density turns out to be a factor 2.2 higher than what measured by the PG survey, corresponding to a surface density of 0.013+/-.006 deg^{-2}. If the Edinburgh QSO Survey is included, an overdensity of a factor 2.7 is observed, corresponding to a density of 0.016+/-0.005 deg^{-2}. In order to derive the QSO optical luminosity function (LF) we used Monte Carlo simulations that take into account of the selection criteria, photometric errors and QSO spectral slope distribution. The LF can be represented with a Pure Luminosity Evolution (L(z)\propto(1+z)^k) of a two power law both for q_0=0.5 and q_0=0.1. For q_0=0.5 k=3.26, slower than the previous Boyle's (1992) estimations of k=3.45. A flatter slope beta=-3.72 of the bright part of the LF is also required. The observed overdensity of bright QSOs is concentrated at z<0.6. It results that in the range 0.3<z<0.6 the luminosity function is flatter than observed at higher redshifts. In this redshift range, for Mb<-25, 32 QSOs are observed instead of 19 expected from our best-fit PLE model. This feature requires a luminosity dependent luminosity evolution in order to satisfactorily represent the data in the whole 0.3<z<2.2 interval.

The long-term variability of a sample of 149 optically selected QSOs in the field of the Selected Area 94 has been studied in the R-band. The relations between variability and luminosity and between variability and redshift have been investigated by means of ``robust'' statistical estimators, allowing to disentangle the effects of the measurement errors. The results are compared with the corresponding properties of the variability in the B-band for the same sample. An anti-correlation between the R-band variability and the intrinsic luminosity is found, analogously to what is observed in the B-band. The amplitude of the R-band variability turns out to be smaller (of a factor $1.13 \pm 0.05$) than the B-band variability. The implications in terms of the black-hole, starburst and microlensing models are discussed.

An ESO Key programme dedicated to an Homogeneous Bright QSO Survey (HBQS) has been completed. 327 QSOs (Mb<-23, 0.3<z<2.2) have been selected over 555 deg^2 with 15<B<18.75. For B<16.4 the QSO surface density turns out to be a factor 2.2 higher than what measured by the PG survey, corresponding to a surface density of 0.013+/-.006 deg^{-2}. If the Edinburgh QSO Survey is included, an overdensity of a factor 2.5 is observed, corresponding to a density of 0.016+/-0.005 deg^{-2}. In order to derive the QSO optical luminosity function (LF) we used Monte Carlo simulations that take into account of the selection criteria, photometric errors and QSO spectral slope distribution. The LF can be represented with a Pure Luminosity Evolution (L(z)\propto(1+z)^k) of a two power law both for q_0=0.5 and q_0=0.1. For q_0=0.5 k=3.26, slower than the previous Boyle's (1992) estimations of k=3.45. A flatter slope beta=-3.72 of the bright part of the LF is also required. The observed overdensity of bright QSOs is concentrated at z<0.6. It results that in the range 0.3<z<0.6 the luminosity function is flatter than observed at higher redshifts. In this redshift range, for Mb<-25, 32 QSOs are observed instead of 19 expected from our best-fit PLE model. This feature requires a luminosity dependent luminosity evolution in order to satisfactorily represent the data in the whole 0.3<z<2.2 interval.

The long-term variability of a sample of 486 optically selected QSOs in the fields of the SA94, SA57 and SGP has been studied. The relations of variability with luminosity and redshift have been investigated by means of statistical estimators that are ``robust'' and allow at the same time to eliminate the influence of the measurement errors. The analysis in the QSOs rest frame of both the ensemble structure function ($SF$) and individual variability indices show: 1) a negative correlation between variability and luminosity; 2) a significant positive correlation between variability and $z$; 3) such correlations may be equally well parameterized either with a model in which the timescale of the variability is fixed for all the QSOs ($\tau \sim 2.4$ yr), while the amplitude linearly increases with the absolute magnitude and $z$, or with a model in which the timescale of the variability linearly depends on the absolute magnitude and the amplitude is only a function of $z$. The same analysis carried out in the observer's frame shows: 1) a negative correlation between variability and luminosity; 2) the timescale of variability does not depend significantly either on the absolute magnitude or on $z$; 3) the ensemble structure function is well represented by a parameterization in which, with a fixed timescale of about 5.5 yr, the amplitude linearly increases with the absolute magnitude; 4) although the general behaviour of the $SF$ does not show a systematic variation of the timescale and/or amplitude with $z$, if we examine the average variability index for objects with $-25>M_B>-27$, we find that below $z=1$ QSOs are significantly less variable than at higher redshift. Implications in terms of the black-hole, starburst and microlensing models are briefly discussed.

We have investigated the spatial distribution of quasars and its relationship with redshift by using the two-point correlation function, the variance of cell counts and the conditional density as a function of redshift. By comparing our results with those found by other authors we conclude that quasars have a correlation function with amplitude similar to or greater than that of galaxies, but the issue of its evolution still needs to be settled. We discuss how this can be achieved by adding new observations and comparing real data with survey simulations.

This is the first paper in a series aimed at defining a statistically significant sample of QSOs in the range $ 15 < B < 18.75$ and $ 0.3 < z < 2.2$. The selection is carried out using direct plates obtained at the ESO and UK Schmidt Telescopes, scanned with the COSMOS facility and searched for objects with an ultraviolet excess. Follow-up spectroscopy, carried out at ESO La Silla, is used to classify each candidate. In this initial paper, we describe the scientific objectives of the survey; the selection and observing techniques used. We present the first sample of 285 QSOs ($M_B < -23$) in a 153 deg$^2$ area, covered by the six ``deep'' fields, intended to obtain significant statistics down $B \simeq 18.75$ with unprecedented photometric accuracy. From this database, QSO counts are determined in the magnitude range $ 17 < B < 18.75$.

The issue of the X-ray to optical luminosity relationship ($L_o-L_x$) is addressed for both optically and X-ray selected quasar samples. We have applied a generalized regression algorithm for the case of samples involving censored data, with errors on both the dependent and the independent variable. Contrary to some previous results, we find that such relationship is consistent with being a linear one ($L_x\propto L_o$). We argue that previous reports of non-linear relationships (i.e. $L_x\propto L_o^e$ with $e < 1$) are due to the neglect of the influence of the photometric errors, the precise knowledge of which strongly influences the reliability of the results. Further progresses in the determination of the $L_o-L_x$ relationship can be achieved with ROSAT observations of the new generation of large bright quasar surveys.

The present status of the ESO key-programme ``A Homogeneous Bright QSO Survey'' is described together with the first results concerning the QSO counts and the optical luminosity function. The related analysis of the evolution of the radio and X-ray selected QSOs is presented.

The sky subtraction performances of multi-fiber spectrographs are discussed, analysing in detail the case of the OPTOPUS system at the 3.6 meter ESO telescope at La Silla. A standard technique, based on flat-fields obtained with a uniformly illuminated screen on the dome, provides poor results. A new method has been developed, using the [OI] emission line at 5577 \AA\ as a calibrator of the fiber transmittance, taking into account the diffuse light and the influence of each fiber on the adjacent ones, and correcting for the effects of the image distortions on the sky sampling. In this way the accuracy of the sky subtraction improves from $2-8 \%$ to $1.3-1.6 \%$.

The average flux decrement shortward the Ly$_{\alpha}$ emission, due to the well-known ``forest'' of absorptions, has been measured in the spectra of 8 quasars. Quasi-simultaneous optical and IUE observations of the two low redshift quasars PKS 0637--75 (z=0.654) and MC 1104+16 (z=0.632) have been carried out, obtaining relatively high S/N, spectrophotometrically calibrated data on their energy distribution from the rest frame H$_{\beta}$ to the Lyman continuum. Six more quasars in the redshift range 2.5-3.4 have been observed in the optical domain. For all the quasars the ``intrinsic'' continuum slope and normalization have been estimated longward the Ly$_{\alpha}$ emission and extrapolated towards the Lyman continuum to measure the average depressions, which have been compared with the model statistics of the Ly$_{\alpha}$ clouds. When all the known classes of absorbers are taken into account with plausible values for their equivalent width distribution and evolution, a good agreement is obtained with the observations. The results for the observed continuum decrement at $z \sim 0.65$ are identical to those predicted by the evolution with redshift of the number of Ly$_{\alpha}$ forest systems including the HST data and within $2\sigma$ of the predicted value using the ``standard'' Ly$_{\alpha}$ evolution (as determined only at high z).