We present $NuSTAR$ observations of four active galactic nuclei (AGN) located within 15 Mpc. These AGN, namely ESO 121-G6, NGC 660, NGC 3486 and NGC 5195, have observed X-ray luminosities of $L_{\rm 2-10\ keV, obs} \lesssim$ 10$^{39}$ erg s$^{-1}$, classifying them as low luminosity AGN (LLAGN). We perform broadband X-ray spectral analysis for the AGN by combining our $NuSTAR$ data with $Chandra$ or $XMM-Newton$ observations to directly measure their column densities ($N_{\rm H}$) and infer their intrinsic power. We complement our X-ray data with archival and new high angular resolution mid-infrared (mid-IR) data for all objects, except NGC 5195. Based on our X-ray spectral analysis, we found that both ESO 121-G6 and NGC 660 are heavily obscured ($N_{\rm H}$ > 10$^{23}$ cm$^{-2}$; $L_{\rm 2-10\ keV,\ int} \sim$ 10$^{41}$ erg s$^{-1}$), and NGC 660 may be Compton-thick. We also note that the X-ray flux and spectral slope for ESO 121-G6 have significantly changed over the last decade, indicating significant changes in the obscuration and potentially accretion rate. On the other hand, NGC 3486 and NGC 5195 appear to be unobscured and just mildly obscured, respectively, with $L_{\rm 2-10\ keV,\ int} <$ 10$^{39}$ erg s$^{-1}$; i.e., genuine LLAGN. Both of the heavily obscured AGN have $L_{\rm bol} >$ 10$^{41}$ erg s$^{-1}$ and $\lambda_{\rm Edd} \gtrsim$ 10$^{-3}$, and are detected in high angular resolution mid-IR imaging, indicating the presence of obscuring dust on nuclear scale. NGC 3486 however, is undetected in high-resolution mid-IR imaging, and the current data do not provide stringent constraints on the presence or absence of obscuring nuclear dust in the AGN.

Hyperluminous quasars ($L_{\rm bol}\gtrsim 10^{47}$ erg s$^{-1}$) are ideal laboratories to study the interaction and impact of extreme radiative field and the most powerful winds in the AGN nuclear regions. They typically exhibit low coronal X-ray luminosity ($L_{\rm X}$) compared to the UV and MIR radiative outputs ($L_{\rm UV}$ and $L_{\rm MIR}$) with a non-negligible fraction of them reporting even $\sim$1 dex weaker $L_{\rm X}$ compared to the prediction of the well established $L_{\rm X}$-$L_{\rm UV}$ and $L_{\rm X}$-$L_{\rm MIR}$ relations followed by the bulk of the AGN population. We report in our WISE/SDSS-selected Hyperluminous (WISSH) $z=2-4$ broad-line quasar sample, the discovery of a dependence between the intrinsic 2-10 keV luminosity ($L_{\rm 2-10}$) and the blueshifted velocity of the CIV emission line ($v_{\rm CIV}$) indicative of accretion disc winds. In particular, sources with fastest winds ($v_{\rm CIV}\gtrsim 3000~\rm km s^{-1}$) possess $\sim$0.5-1 dex lower $L_{\rm 2-10}$ than sources with negligible $v_{\rm CIV}$. No similar dependence is found on $L_{\rm UV}$, $L_{\rm MIR}$, $L_{\rm bol}$, photon index and absorption column density. We interpret these findings in the context of accretion disc wind models. Both magnetohydrodynamic and line-driven models can qualitatively explain the reported relations as a consequence of X-ray shielding from the inner wind regions. In case of line-driven winds, the launch of fast winds is favoured by a reduced X-ray emission, and we speculate that these winds may play a role in directly limiting the coronal hard X-ray production.

During the last years, Ly$\alpha$ nebulae have been routinely detected around high-z, radio-quiet quasars (RQQs) thanks to the advent of sensitive integral field spectrographs. Constraining the physical properties of the Ly$\alpha$ nebulae is crucial for a full understanding of the circum-galactic medium (CGM), which is a venue of feeding and feedback processes. The most luminous quasars are privileged test-beds to study these processes, given their large ionizing fluxes and dense CGM environments in which they are expected to be embedded. We aim at characterizing the rest-frame UV emission lines in the CGM around a hyper-luminous, broad emission line, RQQ at z~3.6, that exhibits powerful outflows at both nuclear and host galaxy scales. We analyze VLT/MUSE observations of the quasar J1538+08 and perform a search for extended UV emission lines to characterize its morphology, emissivity, kinematics and metal content. We report the discovery of a very luminous ($\sim2 \times10^{44}~erg~s^{-1}$), giant Ly$\alpha$ nebula and a likely associated extended CIV nebula. The Ly$\alpha$ nebula emission exhibits moderate blueshift compared with the quasar systemic redshift and large average velocity dispersion ($\sigma_{v}$ ~700 $km~s^{-1}$) across the nebula, while the CIV nebula shows $\sigma_{v}$~$350~km~s^{-1}$. The Ly$\alpha$ line profile exhibits a significant asymmetry towards negative velocity values at 20-30 kpc south of the quasar and is well parameterized by two Gaussian components: a narrow ($\sigma$~$470~km~s^{-1}$) systemic one plus a broad ($\sigma$~1200 $km~s^{-1}$), blueshifted (~1500 $km~s^{-1}$) one. Our analysis of the MUSE observation of J1538+08 reveals metal-enriched CGM around this hyper-luminous quasar and our detection of blueshifted emission in the line profile of the Ly$\alpha$ nebula suggests that powerful nuclear outflows can propagate through the CGM over tens of kpc.

Outflows are observed in a variety of astrophysical sources. Remarkably, ultra-fast ($v\geq 0.1c$), outflows in the UV and X-ray bands are often seen in AGNs. Depending on their energy and mass outflow rate, respectively $\dot{E}_{out}, \dot{M}_{out}$, such outflows may play a key role in regulating the AGN-host galaxy co-evolution process through cosmic time. It is therefore crucial to provide accurate estimates of the wind properties. Here, we concentrate on special relativistic effects concerning the interaction of light with matter moving at relativistic speed relatively to the source of radiation. Our aim is to assess the impact of these effects on the observed properties of the outflows and implement a relativistic correction in the existing spectral modelling routines. We define a simple procedure to incorporate relativistic effects in radiative transfer codes. Following this procedure, we run a series of simulations to explore the impact of these effects on the simulated spectra, for different $v$ and column densities of the outflow. The observed optical depth is usually considered a proxy for the wind $N_H$, independently on its velocity. However, our simulations show that the observed optical depth of an outflow with a given column density $N_H$ decreases rapidly as the velocity of the wind approaches relativistic values. This, in turn, implies that when estimating $N_H$ from the optical depth, it is necessary to include a velocity-dependent correction, already for moderate velocities (e.g. $v \geq 0.05c$). This correction linearly propagates to the derived $\dot{M}_{out}, \dot{E}_{out}$. As an example of these effects, we calculate the relativistically corrected values of $\dot{M}_{out}$ and $\dot{E}_{out}$ for a sample of $\sim 30$ Ultra-Fast Outflows taken from the literature, and find correction factors of $20-120 \%$ within the observed range of outflowing velocities.

Metals form an essential part of the Universe at all scales. Without metals we would not exist, and the Cosmos would look completely different. Metals are primarily born through nuclear processes in stars. They leave their cradles through winds or explosions, and then start their journey through space. This can lead them in and out of astronomical objects on all scales, ranging from comets, planets, stars, entire galaxies, groups and clusters of galaxies to the largest structures of the Universe. Their wanderings are fundamental in determining how these objects, and the entire universe, evolve. In addition, their bare presence can be used to trace what these structures look like. The scope of this paper is to highlight the most important open astrophysical problems that will be central in the next decades and for which a deep understanding of the Universe-wandering metals, their physical and kinematical states and their chemical composition represents the only viable solution. The majority of these studies can only be efficiently performed through High Resolution Spectroscopy in the soft X-ray band.

The WISSH quasars project aims at studying the nuclear and host galaxy properties of the most luminous quasars ($L_{bol}>10^{47}$ erg/s, $1.8<z<4.6$). Nuclear winds are manifested as UV broad ($\geq$ 2,000 km/s) absorption lines (BAL) in $\sim$ 15\% of quasars. We aim at studying the incidence and properties of such winds in the WISSH sample, to investigate possible differences with respect to lower luminosity AGN regimes. We collected optical spectra from the SDSS data release 12, and identified those showing absorption troughs in the region between \siiv and \civ emission lines. We find a higher observed fraction of \civ BAL quasars in the WISSH sample (24\%) with respect to previous catalogues (10-15\%). These WISSH BAL quasars are also characterized by a larger average BI ($\sim$4,000 km/s) and maximum velocity ($\sim$17,000 km/s). Moreover, for two objects we discovered BAL features bluewards of the \siiv peak, which can be associated to \civ absorption with velocity of 0.15c. Finally, we updated previous studies on the dependence of maximum outflow velocity upon bolometric luminosity, showing that BAL winds have intermediate properties compared to molecular/ionized winds and ultra fast outflows (UFOs). Finally, the radio properties of the WISSH BAL quasars as a whole are in line with those of samples at lower luminosities from previous studies. Our results suggest that the higher $L_{\rm bol}$ of the WISSH quasars likely favours the acceleration of BAL outflows and that their most likely driving mechanism is radiation pressure. Furthermore we estimate that the kinetic power associated to these winds in hyperluminous quasars is sufficient, for the highest column density and fastest winds, to provide efficient feedback onto the host galaxy.

We report on the first ALMA observation of the CO(3$-$2) and rest-frame ~340 GHz continuum emission in PDS 456, which is the most luminous, radio-quiet QSO in the local Universe ($z$~0.18), with a bolometric luminosity $L_{\rm Bol}\sim10^{47}$ erg s$^{-1}$. ALMA angular resolution allowed us to map scales as small as ~700 pc. The molecular gas reservoir, traced by the core of the very bright CO(3$-$2) emission line, is distributed in a compact rotating disk, with size of ~1.3 kpc, seen close to face-on ($i$~25 deg). Fast CO(3$-$2) emission in the velocity range $v\in[-1000,+500]$ km s$^{-1}$ is also present. Specifically, we detect several blue-shifted clumps out to ~5 kpc from the nucleus, in addition to a compact ($R\lesssim1.2$ kpc), broad emission component. These components reveal a galaxy-wide molecular outflow, with a total mass $M_{\rm mol}^{\rm out}\sim2.5\times10^8$ $M_{\odot}$ and a mass outflow rate $\dot{M}_{\rm mol}\sim290$ $M_{\odot}$ yr$^{-1}$. The corresponding depletion time is ~8 Myr, shorter than the rate at which the molecular gas is converted into stars, indicating that the detected outflow is potentially able to quench star-formation in the host. The momentum flux of the molecular outflow normalised to the radiative momentum output (i.e. $L_{\rm Bol}/c$) is $\lesssim1$, comparable to that of the X-ray ultra-fast outflow (UFO) detected in PDS 456. This is at odds with the expectations for an energy-conserving expansion suggested for most of the large-scale outflows detected in low-luminosity AGN so far. We suggest three possible scenarios that may explain this observation: (i) in very luminous AGN such as our target the molecular gas phase is tracing only a fraction of the total outflowing mass; (ii) a small coupling between the shocked gas by the UFO and the host-galaxy ISM (iii) AGN radiation pressure may play an important role in driving the outflow.

We present an analysis of a sample of 69 local obscured Swift/Burst Alert Telescope active galactic nuclei (AGNs) with X-ray spectra from NuSTAR and infrared (IR) spectral energy distributions from Herschel and WISE. We combine this X-ray and IR phenomenological modeling and find a significant correlation between reflected hard X-ray emission and IR AGN emission, with suggestive indications that this correlation may be stronger than the one between intrinsic hard X-ray and IR emissions. This relation between the IR and reflected X-ray emission suggests that both are the result of the processing of intrinsic emission from the corona and accretion disk by the same structure. We explore the resulting implications on the underlying distribution of covering fraction for all AGNs, by generating mock observables for the reflection parameter and IR luminosity ratio using empirical relations found for the covering fraction with each quantity. We find that the observed distributions of the reflection parameter and IR-to-X-ray ratio are reproduced with broad distributions centered around covering fractions of at least ~40%-50%, whereas narrower distributions match our observations only when centered around covering fractions of ~70%-80%. Our results are consistent with both independent estimates of the covering fractions of individual objects and the typical covering fraction obtained on the basis of obscured fractions for samples of AGNs. These results suggest that the level of reprocessing in AGNs, including X-ray reflection, is related in a relatively straightforward way to the geometry of the obscuring material.

Outflows from active galactic nuclei (AGN) are often invoked to explain the co-evolution of AGN and their host galaxies, and the scaling relations between the central black hole mass and the bulge velocity dispersion. Nuclear winds are often seen in the X-ray spectra through Fe K shell transitions and some of them are called ultra fast outflows (UFOs) due to their high velocities, up to some fractions of the speed of light. If they were able to transfer some percentage of the AGN luminosity to the host galaxy, this might be enough to trigger an efficient feedback mechanism. We aim to establish new constraints on the covering fraction and on the kinematic properties of the UFO in the powerful (L(bol) ~ 10^(47) erg/s) quasar PDS 456, an established Rosetta stone for studying AGN feedback from disk winds. This will allow us to estimate the mass outflow rate and the energy transfer rate of the wind, which are key quantities to understand the potential impact on the host galaxy. We analyze two sets of simultaneous XMM-Newton and NuSTAR observations taken in September 2013 and reported in Nardini et al. (2015) as having similar broadband spectral properties. We fit the Fe K features with a P-Cygni profile between 5 and 14 keV, using a novel Monte Carlo model for the WINd Emission (WINE). We find an outflow velocity ranging from 0.17 to 0.28 c, with a mean value of 0.23 c. We obtain an opening angle of the wind of 71(+13,-8) deg and a covering fraction of 0.7(+0.2,-0.3), suggesting a wide-angle outflow. We check the reliability of the WINE model by performing extensive simulations of joint XMM-Newton and NuSTAR observations. Furthermore, we test the accuracy of the WINE model in recovering the geometrical properties of UFOs by simulating observations with the forthcoming X-ray observatory ATHENA.

Theoretical models of galaxy formation suggest that the presence of an active galactic nucleus (AGN) is required to regulate the growth of its host galaxy through feedback mechanisms, produced by e.g. AGN-driven outflows. Although such outflows are common both at low and high redshift, a comprehensive picture is still missing. The peak epoch of galaxy assembly (1<z<3) has been poorly explored so far, and current observations in this redshift range are mostly limited to targets with high chances to be in an outflowing phase. This paper introduces SUPER (a SINFONI Survey for Unveiling the Physics and Effect of Radiative feedback), an ongoing ESO's VLT/SINFONI Large Programme. SUPER will perform the first systematic investigation of ionized outflows in a sizeable and blindly-selected sample of 39 X-ray AGN at z~2, which reaches high spatial resolutions (~2 kpc) thanks to the adaptive optics-assisted IFU observations. The outflow morphology and star formation in the host galaxy will be mapped through the broad component of [OIII] and the narrow component of Ha emission lines. The main aim of our survey is to infer the impact of outflows on the on-going star formation and to link the outflow properties to a number of AGN and host galaxy properties. We describe here the survey characteristics and goals, as well as the selection of the target sample. Moreover, we present a full characterization of its multi-wavelength properties: we measure, via spectral energy distribution fitting of UV-to-FIR photometry, stellar masses (4x10^9-2x10^11 Msun), star formation rates (25-680 Msun yr^-1) and AGN bolometric luminosities (2x10^44-8x10^47 erg s^-1), along with obscuring column densities (<2x10^24 cm^-2) and 2-10 keV luminosities (2x10^43-6x10^45 erg s^-1) derived through X-ray spectral analysis. Finally, we classify our AGN as jetted or non-jetted according to their radio and FIR emission.

We present ALMA observations of the CO(6-5) and [CII] emission lines and the sub-millimeter continuum of the $z\sim6$ quasi-stellar object (QSO) SDSS J231038.88+185519.7. Compared to previous studies, we have analyzed a synthetic beam that is ten times smaller in angular size, we have achieved ten times better sensitivity in the CO(6-5) line, and two and half times better sensitivity in the [CII] line, enabling us to resolve the molecular gas emission. We obtain a size of the dense molecular gas of $2.9\pm0.5$ kpc, and of $1.4\pm0.2$ kpc for the 91.5 GHz dust continuum. By assuming that CO(6-5) is thermalized, and by adopting a CO--to--$H_2$ conversion factor $\rm \alpha_{CO} = 0.8~ M_{\odot}~K^{-1}~ (km/s)^{-1} ~pc^{2}$, we infer a molecular gas mass of $\rm M(H_2)=(3.2 \pm0.2) \times 10^{10}\rm M_{\odot}$. Assuming that the observed CO velocity gradient is due to an inclined rotating disk, we derive a dynamical mass of $\rm M_{dyn}~sin^2(i) = (2.4\pm0.5) \times 10^{10}~ M_{\odot}$, which is a factor of approximately two smaller than the previously reported estimate based on [CII]. Regarding the central black hole, we provide a new estimate of the black hole mass based on the C~IV emission line detected in the X-SHOOTER/VLT spectrum: $\rm M_{BH}=(1.8\pm 0.5) \times 10^{9}~ M_{\odot}$. We find a molecular gas fraction of $\rm \mu=M(H_2)/M^*\sim4.4$, where $\rm M^*\approx M_{dyn} - M(H_2)-M(BH)$. We derive a ratio $v_{rot}/\sigma \approx 1-2$ suggesting high gas turbulence, outflows/inflows and/or complex kinematics due to a merger event. We estimate a global Toomre parameter $Q\sim 0.2-0.5$, indicating likely cloud fragmentation. We compare, at the same angular resolution, the CO(6-5) and [CII] distributions, finding that dense molecular gas is more centrally concentrated with respect to [CII]. We find that the current BH growth rate is similar to that of its host galaxy.

We present a 155ks NuSTAR observation of the $z\sim2$ hot dust-obscured galaxy (hot DOG) W1835+4355. We extracted spectra from the two NuSTAR detectors and analyzed them jointly with the archival XMM PN and MOS spectra. We performed a spectroscopic analysis based on both phenomenological and physically motivated models employing toroidal and spherical geometry for the obscurer. In all the modelings, the source exhibits a Compton-thick column density $N_{\rm H} \gtrsim 10^{24}$ cm$^{-2}$, a 2-10 keV luminosity $L_{2-10}\approx2\times10^{45}$ erg s$^{-1}$ , and a prominent soft excess ($\sim5-10$ % of the primary radiative output), which translates into a luminosity $\sim10^{44}$ erg s$^{-1}$. We modeled the spectral energy distribution from 1.6 to 850 $\mu m$ using a clumpy two-phase dusty torus model plus a modified blackbody to account for emission powered by star formation in the far-infrared. We employed several geometrical configurations consistent with those applied in the X-ray analysis. In all cases we obtained a bolometric luminosity $L_{\rm bol}\approx3-5\times10^{47}$ erg s$^{-1}$, which confirms the hyperluminous nature of this active galactic nucleus. Finally, we estimate a prodigious star formation rate of $\sim$3000 $M_{\odot}\,yr^{-1}$, which is consistent with the rates inferred for $z\approx2-4$ hyperluminous type I quasars. The heavily obscured nature, together with $L_{\rm bol}$, the ratio of X-ray to mid-infrared luminosity, the rest-frame optical morphology, and the host star formation rate are indicative of its evolutionary stage. We can interpret this as a late-stage merger event in the transitional, dust-enshrouded, evolutionary phase eventually leading to an optically bright AGN.

It has been known for decades that the observed number of baryons in the local universe falls about 30-40% short of the total number of baryons predicted by Big-Bang Nucleosynthesis, as inferred from density fluctuations of the Cosmic Microwave Background and seen during the first 2-3 billion years of the universe in the so called Lyman-alpha Forest. A theoretical solution to this paradox locates the missing baryons in the hot and tenuous filamentary gas between galaxies, known as the warm-hot intergalactic medium. However, it is difficult to detect them there because the largest by far constituent of this gas - hydrogen - is mostly ionized and therefore almost invisible in far-ultraviolet spectra with typical signal-to-noise ratios. Indeed, despite the large observational efforts, only a few marginal claims of detection have been made so far. Here we report observations of two absorbers of highly ionized oxygen (OVII) in the high signal-to-noise-ratio X-ray spectrum of a quasar at redshift >0.4. These absorbers show no variability over a 2-year timescale and have no associated cold absorption, making the assumption that they originate from the quasar's intrinsic outflow or the host galaxy's interstellar medium implausible. The OVII systems lie in regions characterized by large (x4 compared to average) galaxy over-densities and their number (down to the sensitivity threshold of our data), agrees well with numerical simulation predictions for the long-sought warm-hot intergalactic medium (WHIM). We conclude that the missing baryons have been found.

We present an ALMA high-resolution observation of the 840 um continuum and [CII] line emission in the WISE-SDSS selected hyper-luminous (WISSH) QSO J1015+0020 at z~4.4. Our analysis reveals an exceptional overdensity of [CII]-emitting companions with a very small (<150 km/s) velocity shift with respect to the QSO redshift. We report the discovery of the closest companion observed so far in submillimetre observations of high-z QSOs. It is only 2.2 kpc distant and merging with J1015+0020, while two other [CII] emitters are found at 8 and 17 kpc. Two strong continuum emitters are also detected at <3.5 arcsec. They are likely associated to the same overdense structure of J1015+0020, as they exceed by a factor of 100 the number of expected sources, considering the Log(N)-Log(S) at 850 um. The host galaxy of J1015+0020 shows a SFR of about 100 Msun/yr while the total SFR of the QSO and its companion galaxies is a factor of 10 higher, indicating that substantial stellar mass assembly at early epochs may have taken place in the QSO satellites. For J1015+0020 we compute a SMBH mass MBH~6E9 Msun and a dynamical mass Mdyn~4E10 Msun . This translates into an extreme ratio Mdyn/MBH~7. The total stellar mass of the QSO host galaxy plus the [CII] emitters already exceeds 1E11 Msun at z~4.4. These sources will likely merge and develop into a giant galaxy of 1.3E12 Msun. Under the assumption of constant mass or Eddington accretion rate equal to the observed values, we find that the growth timescale of the host galaxy is comparable or even shorter than that inferred for the SMBH.

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.

Observing the light emitted by the first accreting black holes (BHs) would dramatically improve our understanding of the formation of quasars at z > 6, possibly unveiling the nature of their supermassive black hole (SMBH) seeds. In previous works we explored the relative role of the two main competing BH seed formation channels, Population III remnants (low-mass seeds) and direct collapse BHs (high-mass seeds), investigating the properties of their host galaxies in a cosmological context. Building on this analysis, we predict here the spectral energy distribution and observational features of low- and high-mass BH seeds selected among the progenitors of a z~6 SMBH. We derive the processed emission from both accreting BHs and stars by using the photo-ionization code Cloudy, accounting for the evolution of metallicity and dust-to-gas mass ratio in the interstellar medium of the host galaxies, as predicted by the cosmological data- constrained model GAMETE/QSOdust. We show how future missions like JWST and ATHENA will be able to detect the light coming from SMBH progenitors already at z~16. We build upon previous complementary studies and propose a method based on the combined analysis of near infrared (NIR) colors, IR excess (IRX) and UV continuum slopes (i.e. color-color and IRX-Beta diagrams) to distinguish growing seed BH host galaxies from starburst-dominated systems in JWST surveys. Sources selected through this criterion would be the best target for follow-up X-ray observations.

We discuss the spectral analysis of a sample of 63 Active Galactic Nuclei (AGN) detected above a limiting flux of $S(8-24$ keV$)=7\times10^{-14}$ erg/s/cm$^2$ in the multi-tiered NuSTAR Extragalactic Survey program. The sources span a redshift range z=0-2.1 (median $\langle$z$\rangle=$0.58). The spectral analysis is performed over the broad 0.5-24 keV energy range, combining NuSTAR with Chandra and/or XMM-Newton data and employing empirical and physically motivated models. This constitutes the largest sample of AGN selected at $ >10$ keV to be homogeneously spectrally analyzed at these flux levels. We study the distribution of spectral parameters such as photon index, column density ($N_{\rm H}$), reflection parameter ($R$) and 10-40 keV luminosity ($L_{X}$). Heavily obscured ($\log [N_{\rm H}/cm^{-2}]\ge23$) and Compton Thick (CT; $\log [N_{\rm H}/cm^{-2}]\ge24$) AGN constitute $\sim$25% (15-17 sources) and $\sim$2-3% ( 1-2 sources) of the sample, respectively. The observed $N_{\rm H}$ distribution fairly agrees with predictions of Cosmic X-ray Background population synthesis models (CXBPSM). We estimate the intrinsic fraction of AGN as a function of $N_{H}$, accounting for the bias against obscured AGN in a flux-selected sample. The fraction of CT AGN relative to $\log [N_{\rm H}/cm^{-2}]=20-24$ AGN is poorly constrainted, formally in the range 2-56% (90% upper limit of 66%). We derived a fraction ($f_{abs}$) of obscured AGN ($\log [N_{\rm H}/cm^{-2}]=22-24$) as a function of $L_{X}$ in agreement with CXBPSM and previous z < 1 X-ray determinations. Furthermore $f_{abs}$ at z=0.1-0.5 and $\log (L_{X}/erg$ $s^{-1})\approx43.6-44.3$ agrees with observational measurements/trends obtained over larger redshift intervals. We report a significant anti-correlation of $R$ with $L_{X}$ (confirmed by our companion paper on stacked spectra) with considerable scatter around the median $R$ values.

We present the results and the source catalog of the NuSTAR survey in the UKIDSS Ultra Deep Survey (UDS) field, bridging the gap in depth and area between NuSTAR's ECDFS and COSMOS surveys. The survey covers a $\sim 0.6$ deg$^2$ area of the field for a total observing time of $\sim$ 1.75 Ms, to a half-area depth of $\sim$ 155 ks corrected for vignetting at $3-24$ keV, and reaching sensitivity limits at half-area in the full ($3-24$ keV), soft ($3-8$ keV) and hard ($8-24$ keV) bands of $2.2 \times 10^{-14}$ erg cm$^{-2}$ s$^{-1}$, $1.0 \times 10^{-14}$ erg cm$^{-2}$ s$^{-1}$, and $2.7 \times 10^{-14}$ erg cm$^{-2}$ s$^{-1}$, respectively. A total of 67 sources are detected in at least one of the three bands, 56 of which have a robust optical redshift with a median of $\langle z\rangle \sim 1.1$. Through a broadband ($0.5-24$ keV) spectral analysis of the whole sample combined with the NuSTAR hardness ratios, we compute the observed Compton-thick (CT; $N_{\rm H} > 10^{24}$ cm$^{-2}$) fraction. Taking into account the uncertainties on each $N_{\rm H}$ measurement, the final number of CT sources is $6.8\pm1.2$. This corresponds to an observed CT fraction of $(11.5\pm2.0)\%$, providing a robust lower limit to the intrinsic fraction of CT AGN and placing constraints on cosmic X-ray background synthesis models.

Supermassive black holes (SMBHs) of 10^9-10^10 Msun were already in place ~13 Gyr ago, at z>6. Super-Eddington growth of low-mass BH seeds (~100 Msun) or less extreme accretion onto ~10^5 Msun seeds have been recently considered as the main viable routes to these SMBHs. Here we study the statistics of these SMBH progenitors at z~6. The growth of low- and high-mass seeds and their host galaxies are consistently followed using the cosmological data constrained model GAMETE/QSOdust, which reproduces the observed properties of high-z quasars, like SDSS J1148+5251. We show that both seed formation channels can be in action over a similar redshift range, 15 < z < 18 and are found in dark matter halos with comparable mass, ~5x10^7 Msun. However, as long as the systems evolve in isolation (i.e. no mergers occur), noticeable differences in their properties emerge: at z>= 10 galaxies hosting high-mass seeds have smaller stellar mass and metallicity, the BHs accrete gas at higher rates and star formation proceeds less efficiently than in low-mass seeds hosts. At z<10 these differences are progressively erased, as the systems experience minor or major mergers and every trace of the BH origin gets lost.

We present a study of the average X-ray spectral properties of the sources detected by the NuSTAR extragalactic survey, comprising observations of the E-CDFS, EGS and COSMOS fields. The sample includes 182 NuSTAR sources (64 detected at 8-24 keV), with 3-24 keV fluxes ranging between $f_{\rm 3-24 keV}\approx10^{-14}$ and $6\times10^{-13}$ erg/cm$^2$/s ($f_{\rm 8-24 keV}\approx3\times10^{-14}-3\times10^{-13}$ erg/cm$^2$/s) and redshifts of $z=0.04-3.21$. We produce composite spectra from the Chandra+NuSTAR data ($E\approx2-40$ keV, rest frame) for all the sources with redshift identifications (95%) and investigate the intrinsic, average spectra of the sources, divided into broad-line (BL) and narrow-line (NL) AGN, and also in different bins of X-ray column density and luminosity. The average power-law photon index for the whole sample is $\Gamma=1.65_{-0.03}^{+0.03}$, flatter than $\Gamma\approx1.8$ typically found for AGN. While the spectral slope of BL and X-ray unabsorbed AGN is consistent with typical values ($\Gamma=1.79_{-0.01}^{+0.01}$), a significant flattening is seen in NL AGN and heavily-absorbed sources ($\Gamma=1.60_{-0.05}^{+0.08}$ and $\Gamma=1.38_{-0.12}^{+0.12}$, respectively), likely due to the effect of absorption and to the contribution from Compton reflection to the high-energy flux (E>10 keV). We find that the typical reflection fraction in our spectra is $R\approx0.5$ (for $\Gamma=1.8$), with a tentative indication of an increase of the reflection strength with column density. While there is no significant evidence for a dependence of the photon index with X-ray luminosity in our sample, we find that $R$ decreases with luminosity, with relatively high levels of reflection ($R\approx1.2$) for $L_{\rm 10-40 keV}<10^{44}$ erg/s and $R\approx0.3$ for $L_{\rm 10-40 keV}>10^{44}$ erg/s AGN, assuming $\Gamma=1.8$.

We have performed a high sensitivity observation of the UFO/BAL quasar APM 08279+5255 at z=3.912 with NOEMA at 3.2 mm, aimed at detecting fast moving molecular gas. We report the detection of blueshifted CO(4-3) with maximum velocity (v95\%) of $-1340$ km s$^{-1}$, with respect to the systemic peak emission, and a luminosity of $L' = 9.9\times 10^9 ~\mu^{-1}$ K km s$^{-1}$ pc$^{-2}$ (where $\mu$ is the lensing magnification factor). We discuss various scenarios for the nature of this emission, and conclude that this is the first detection of fast molecular gas at redshift $>3$. We derive a mass flow rate of molecular gas in the range $\rm \dot M=3-7.4\times 10^3$ M$_\odot$/yr, and momentum boost $\dot P_{OF} / \dot P_{AGN} \sim 2-6$, therefore consistent with a momentum conserving flow. For the largest $\dot P_{OF}$ the scaling is also consistent with a energy conserving flow with an efficiency of $\sim$10-20\%. The present data can hardly discriminate between the two expansion modes. The mass loading factor of the molecular outflow $\eta=\dot M_{OF}/SFR$ is $>>1$. We also detect a molecular emission line at a frequency of 94.83 GHz, corresponding to a rest frame frequency of 465.8 GHz, which we tentatively identified with the cation molecule $\rm N_2H^+$(5-4), which would be the first detection of this species at high redshift. We discuss the alternative possibility that this emission is due to a CO emission line from the, so far undetected, lens galaxy. Further observations of additional transitions of the same species with NOEMA can discriminate between the two scenarios.

We identify sources with extremely hard X-ray spectra (i.e., with photon indices of Gamma<0.6 in the 13 sq. deg. NuSTAR serendipitous survey, to search for the most highly obscured AGNs detected at >10 keV. Eight extreme NuSTAR sources are identified, and we use the NuSTAR data in combination with lower energy X-ray observations (from Chandra, Swift XRT, and XMM-Newton) to characterize the broad-band (0.5-24 keV) X-ray spectra. We find that all of the extreme sources are highly obscured AGNs, including three robust Compton-thick (CT; N_H > 1.5e24 cm^-2) AGNs at low redshift (z<0.1), and a likely-CT AGN at higher redshift (z=0.16). Most of the extreme sources would not have been identified as highly obscured based on the low energy (<10 keV) X-ray coverage alone. The multiwavelength properties (e.g., optical spectra and X-ray/MIR luminosity ratios) provide further support for the eight sources being significantly obscured. Correcting for absorption, the intrinsic rest-frame 10-40 keV luminosities of the extreme sources cover a broad range, from ~ 5 x 10^42 to 10^45 erg s^-1. The estimated number counts of CT AGNs in the NuSTAR serendipitous survey are in broad agreement with model expectations based on previous X-ray surveys, except for the lowest redshifts (z<0.07) where we measure a high CT fraction of f_CT^obs = 30 (+16 -12) %. For the small sample of CT AGNs, we find a high fraction of galaxy major mergers (50 +/- 33%) compared to control samples of "normal" AGNs.

We perform a survey of the X-ray properties of 41 objects from the WISE/SDSS selected Hyper-luminous (WISSH) quasars sample, composed by 86 broad-line quasars (QSOs) with bolometric luminosity $L_{Bol}\geq 2\times 10^{47}\,erg\, s^{-1}$, at z~2-4. All but 3 QSOs show unabsorbed 2-10 keV luminosities $L_{2-10}\geq10^{45} \,erg \,s^{-1}$. Thanks to their extreme radiative output across the Mid-IR-to-X-ray range, WISSH QSOs offer the opportunity to significantly extend and validate the existing relations involving $L_{2-10}$. We study $L_{2-10}$ as a function of (i) X-ray-to-Optical (X/O) flux ratio, (ii) mid-IR luminosity ($L_{MIR}$), (iii) $L_{Bol}$ as well as (iv) $\alpha_{OX}$ vs. the 2500$\mathring{A}$ luminosity. We find that WISSH QSOs show very low X/O(<0.1) compared to typical AGN values; $L_{2-10}/L_{MIR}$ ratios significantly smaller than those derived for AGN with lower luminosity; large X-ray bolometric corrections $k_{\rm Bol,X}\sim$ 100-1000; and steep $-2<\alpha_{OX}<-1.7$. These results lead to a scenario where the X-ray emission of hyper-luminous quasars is relatively weaker compared to lower-luminosity AGN. Models predict that such an X-ray weakness can be relevant for the acceleration of powerful high-ionization emission line-driven winds, commonly detected in the UV spectra of WISSH QSOs, which can in turn perturb the X-ray corona and weaken its emission. Accordingly, hyper-luminous QSOs represent the ideal laboratory to study the link between the AGN energy output and wind acceleration. Additionally, WISSH QSOs show very large BH masses ($\log[M_{\rm BH}/M_{\odot}]$>9.5). This enables a more robust modeling of the $\Gamma-M_{BH}$ relation by increasing the statistics at high masses. We derive a flatter $\Gamma$ dependence than previously found over the broad range 5 <$\log(M_{\rm BH}/M_{\odot})$ < 11.

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 present the 78-ks Chandra observations of the $z=6.4$ quasar SDSS J1148+5251. The source is clearly detected in the energy range 0.3-7 keV with 42 counts (with a significance $\gtrsim9\sigma$). The X-ray spectrum is best-fitted by a power-law with photon index $\Gamma=1.9$ absorbed by a gas column density of $\rm N_{\rm H}=2.0^{+2.0}_{-1.5}\times10^{23}\,\rm cm^{-2}$. We measure an intrinsic luminosity at 2-10 keV and 10-40 keV equal to $\sim 1.5\times 10^{45}~\rm erg~s^{-1}$, comparable with luminous local and intermediate-redshift quasar properties. Moreover, the X-ray to optical power-law slope value ($\alpha_{\rm OX}=-1.76\pm 0.14$) of J1148 is consistent with the one found in quasars with similar rest-frame 2500 \AA ~luminosity ($L_{\rm 2500}\sim 10^{32}~\rm erg~s^{-1}$\AA$^{-1}$). Then we use Chandra data to test a physically motivated model that computes the intrinsic X-ray flux emitted by a quasar starting from the properties of the powering black hole and assuming that X-ray emission is attenuated by intervening, metal-rich ($Z\geq \rm Z_{\odot}$) molecular clouds distributed on $\sim$kpc scales in the host galaxy. Our analysis favors a black hole mass $M_{\rm BH} \sim 3\times 10^9 \rm M_\odot$ and a molecular hydrogen mass $M_{\rm H_2}\sim 2\times 10^{10} \rm M_\odot$, in good agreement with estimates obtained from previous studies. We finally discuss strengths and limits of our analysis.

Feedback from accreting SMBHs is often identified as the main mechanism responsible for regulating star-formation in AGN host galaxies. However, the relationships between AGN activity, radiation, winds, and star-formation are complex and still far from being understood. We study scaling relations between AGN properties, host galaxy properties and AGN winds. We then evaluate the wind mean impact on the global star-formation history, taking into account the short AGN duty cycle with respect to that of star-formation. We first collect AGN wind observations for 94 AGN with detected massive winds at sub-pc to kpc spatial scales. We then fold AGN wind scaling relations with AGN luminosity functions, to evaluate the average AGN wind mass-loading factor as a function of cosmic time. We find strong correlations between the AGN molecular and ionised wind mass outflow rates and the AGN bolometric luminosity. The power law scaling is steeper for ionised winds (slope 1.29+/-0.38) than for molecular winds (0.76+/-0.06), meaning that the two rates converge at high bolometric luminosities. The molecular gas depletion timescale and the molecular gas fraction of galaxies hosting powerful AGN winds are 3-10 times shorter and smaller than those of main-sequence galaxies with similar SFR, stellar mass and redshift. These findings suggest that, at high AGN bolometric luminosity, the reduced molecular gas fraction may be due to the destruction of molecules by the wind, leading to a larger fraction of gas in the atomic ionised phase. The AGN wind mass-loading factor $\eta=\dot M_{OF}/SFR$ is systematically higher than that of starburst driven winds. Our analysis shows that AGN winds are, on average, powerful enough to clean galaxies from their molecular gas only in massive systems at z<=2, i.e. a strong form of co-evolution between SMBHs and galaxies appears to break down for the least massive galaxies.

We present the first full catalog and science results for the NuSTAR serendipitous survey. The catalog incorporates data taken during the first 40 months of NuSTAR operation, which provide ~20Ms of effective exposure time over 331 fields, with an areal coverage of 13 sq deg, and 497 sources detected in total over the 3-24 keV energy range. There are 276 sources with spectroscopic redshifts and classifications, largely resulting from our extensive campaign of ground-based spectroscopic followup. We characterize the overall sample in terms of the X-ray, optical, and infrared source properties. The sample is primarily comprised of active galactic nuclei (AGNs), detected over a large range in redshift from z = 0.002 - 3.4 (median of <z> = 0.56), but also includes 16 spectroscopically confirmed Galactic sources. There is a large range in X-ray flux, from log( f_3-24keV / erg s^-1 cm^-2 ) ~ -14 to -11, and in rest-frame 10-40 keV luminosity, from log( L_10-40keV / erg s^-1 ) ~ 39 to 46, with a median of 44.1. Approximately 79% of the NuSTAR sources have lower energy (<10 keV) X-ray counterparts from XMM-Newton, Chandra, and Swift/XRT. The mid-infrared (MIR) analysis, using WISE all-sky survey data, shows that MIR AGN color selections miss a large fraction of the NuSTAR-selected AGN population, from ~15% at the highest luminosities (Lx > 10^44 erg s^-1) to ~80% at the lowest luminosities (Lx < 10^43 erg s^-1). Our optical spectroscopic analysis finds that the observed fraction of optically obscured AGNs (i.e., the Type 2 fraction) is F_Type2 = 53(+14-15)%, for a well-defined subset of the 8-24 keV selected sample. This is higher, albeit at a low significance level, than the Type 2 fraction measured for redshift- and luminosity-matched AGNs selected by <10 keV X-ray missions.

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.

Hot, Dust-Obscured Galaxies (Hot DOGs), selected from the WISE all sky infrared survey, host some of the most powerful Active Galactic Nuclei (AGN) known, and might represent an important stage in the evolution of galaxies. Most known Hot DOGs are at $z> 1.5$, due in part to a strong bias against identifying them at lower redshift related to the selection criteria. We present a new selection method that identifies 153 Hot DOG candidates at $z\sim 1$, where they are significantly brighter and easier to study. We validate this approach by measuring a redshift $z=1.009$, and an SED similar to higher redshift Hot DOGs for one of these objects, WISE J1036+0449 ($L_{\rm\,Bol}\simeq 8\times 10^{46}\rm\,erg\,s^{-1}$), using data from Keck/LRIS and NIRSPEC, SDSS, and CSO. We find evidence of a broadened component in MgII, which, if due to the gravitational potential of the supermassive black hole, would imply a black hole mass of $M_{\rm\,BH}\simeq 2 \times 10^8 M_{\odot}$, and an Eddington ratio of $\lambda_{\rm\,Edd}\simeq 2.7$. WISE J1036+0449 is the first Hot DOG detected by NuSTAR, and the observations show that the source is heavily obscured, with a column density of $N_{\rm\,H}\simeq(2-15)\times10^{23}\rm\,cm^{-2}$. The source has an intrinsic 2-10 keV luminosity of $\sim 6\times 10^{44}\rm\,erg\,s^{-1}$, a value significantly lower than that expected from the mid-infrared/X-ray correlation. We also find that the other Hot DOGs observed by X-ray facilities show a similar deficiency of X-ray flux. We discuss the origin of the X-ray weakness and the absorption properties of Hot DOGs. Hot DOGs at $z\lesssim1$ could be excellent laboratories to probe the characteristics of the accretion flow and of the X-ray emitting plasma at extreme values of the Eddington ratio.

We compute the non-thermal emissions produced by relativistic particles accelerated by the AGN-driven shocks in NGC 1068, and we compare the model predictions with the observed gamma-ray and radio spectra . The former is contributed by pion decay, inverse Compton scattering, and bremsstrahlung, while the latter is produced by synchrotron radiation. We derive the gamma-ray and radio emissions by assuming the standard acceleration theory, and we discuss how our results compare with those corresponding to other commonly assumed sources of gamma-ray and radio emissions, like Supernova remnants (SNR) or AGN jets. We find that the AGN-driven shocks observed in the circumnuclear molecular disk of such a galaxy provide a contribution to the gamma-ray emission comparable to that provided by the starburst activity when standard particle acceleration efficiencies are assumed, while they can yield the whole gamma-ray emission only when the parameters describing the acceleration efficiency and the proton coupling with the molecular gas are tuned to values larger than those assumed in standard, SNR-driven shocks. We discuss the range of acceleration efficiencies (for protons and electrons) and of proton calorimetric fractions required to account for the observed gamma-ray emission in the AGN outflow model. We further compare the neutrino flux expected in our model with constraints from current experiments, and we provide predictions for the detections by the upcoming KM3NeT neutrino telescope. This analysis strongly motivates observations of NGC 1068 at >TeV energies with current and future Cherenkov telescopes in order to gain insight into the nature of the gamma-rays source.

We present the 3-8 keV and 8-24 keV number counts of active galactic nuclei (AGN) identified in the NuSTAR extragalactic surveys. NuSTAR has now resolved 33-39% of the X-ray background in the 8-24 keV band, directly identifying AGN with obscuring columns up to approximately 1e25 /cm2. In the softer 3-8 keV band the number counts are in general agreement with those measured by XMM-Newton and Chandra over the flux range 5e-15 < S(3 - 8 keV)/(erg/cm2/s) < 1e-12 probed by NuSTAR. In the hard 8-24 keV band NuSTAR probes fluxes over the range 2e-14 < S(8-24 keV)/(erg/cm2/s) < 1e-12, a factor of approximately 100 fainter than previous measurements. The 8-24 keV number counts match predictions from AGN population synthesis models, directly confirming the existence of a population of obscured and/or hard X-ray sources inferred from the shape of the integrated cosmic X-ray background. The measured NuSTAR counts lie significantly above simple extrapolation with a Euclidian slope to low flux of the Swift/BAT 15-55 keV number counts measured at higher fluxes S(15-55 keV) > 1e-11 erg/cm2/s, reflecting the evolution of the AGN population between the Swift/BAT local (z<0.1) sample and NuSTAR's z~1 sample. CXB synthesis models, which account for AGN evolution, lie above the Swift/BAT measurements, suggesting that they do not fully capture the evolution of obscured AGN at low redshifts.

We present the first direct measurements of the rest-frame 10-40 keV X-ray luminosity function (XLF) of Active Galactic Nuclei (AGNs) based on a sample of 94 sources at 0.1 < z <3, selected at 8-24 keV energies from sources in the NuSTAR extragalactic survey program. Our results are consistent with the strong evolution of the AGN population seen in prior, lower-energy studies of the XLF. However, different models of the intrinsic distribution of absorption, which are used to correct for selection biases, give significantly different predictions for the total number of sources in our sample, leading to small, systematic differences in our binned estimates of the XLF. Adopting a model with a lower intrinsic fraction of Compton-thick sources and a larger population of sources with column densities N_H ~ 10^{23-24} /cm2 or a model with stronger Compton reflection component (with a relative normalization of R ~ 2 at all luminosities) can bring extrapolations of the XLF from 2-10 keV into agreement with our NuSTAR sample. Ultimately, X-ray spectral analysis of the NuSTAR sources is required to break this degeneracy between the distribution of absorbing column densities and the strength of the Compton reflection component and thus refine our measurements of the XLF. Furthermore, the models that successfully describe the high-redshift population seen by NuSTAR tend to over-predict previous, high-energy measurements of the local XLF, indicating that there is evolution of the AGN population that is not fully captured by the current models.

To provide the census of the sources contributing to the X-ray background peak above 10 keV, NuSTAR is performing extragalactic surveys using a three-tier "wedding cake" approach. We present the NuSTAR survey of the COSMOS field, the medium sensitivity and medium area tier, covering 1.7 deg2 and overlapping with both Chandra and XMM-Newton data. This survey consists of 121 observations for a total exposure of ~3 Ms. To fully exploit these data, we developed a new detection strategy, carefully tested through extensive simulations. The survey sensitivity at 20% completeness is 5.9, 2.9 and 6.4 x 10^-14 erg/cm2/s in the 3-24 keV, 3-8 keV and 8-24 keV bands, respectively. By combining detections in 3 bands, we have a sample of 91 NuSTAR sources with 10^42 -10^45.5 erg/s luminosities and redshift z=0.04-2.5. Thirty two sources are detected in the 8-24 keV band with fluxes ~100 times fainter than sources detected by Swift-BAT. Of the 91 detections, all but four are associated with a Chandra and/or XMM-Newton point-like counterpart. One source is associated with an extended lower energy X-ray source. We present the X-ray (hardness ratio and luminosity) and optical-to-X-ray properties. The observed fraction of candidate Compton-thick AGN measured from the hardness ratio is between 13-20%. We discuss the spectral properties of NuSTAR J100259+0220.6 (ID 330) at z=0.044, with the highest hardness ratio in the entire sample. The measured column density exceeds 10^24 /cm2, implying the source is Compton-thick. This source was not previously recognized as such without the >10 keV data.

We present initial results and the source catalog from the NuSTAR survey of the Extended Chandra Deep Field South (hereafter, ECDFS) - currently the deepest contiguous component of the NuSTAR extragalactic survey program. The survey covers the full ~30 arcmin x 30 arcmin area of this field to a maximum depth of ~360 ks (~220 ks when corrected for vignetting at 3-24 keV), reaching sensitivity limits of ~1.3 x 10^-14 erg/cm2/s (3-8 keV), ~3.4 x 10^-14 erg/cm2/s (8-24 keV) and ~3.0 x 10^-14 erg/cm2/s (3-24 keV). Fifty four (54) sources are detected over the full field, although five of these are found to lie below our significance threshold once contaminating flux from neighboring (i.e., blended) sources is taken into account. Of the remaining 49 that are significant, 19 are detected in the 8-24 keV band. The 8-24 keV to 3-8 keV band ratios of the twelve sources that are detected in both bands span the range 0.39-1.7, corresponding to a photon index range of Gamma ~ 0.5-2.3, with a median photon index of 1.70 +/- 0.52. The redshifts of the 49 sources in our main sample span the range z = 0.21-2.7, and their rest-frame 10-40 keV luminosities (derived from the observed 8-24 keV fluxes) span the range L(10-40 keV) ~ (0.7-300) x 10^43 erg/s, sampling below the "knee" of the X-ray luminosity function out to z ~ 0.8-1. Finally, we identify one NuSTAR source that has neither a Chandra nor an XMM-Newton counterpart, but that shows evidence of nuclear activity at infrared wavelengths, and thus may represent a genuine, new X-ray source detected by NuSTAR in the ECDFS.

We present a broad-band (~0.3-70 keV) spectral and temporal analysis of NuSTAR observations of the luminous infrared galaxy NGC 6240, combined with archival Chandra, XMM-Newton and BeppoSAX data. NGC 6240 is a galaxy in a relatively early merger state with two distinct nuclei separated by ~1."5. Previous Chandra observations have resolved the two nuclei, showing that they are both active and obscured by Compton-thick material. Although they cannot be resolved by NuSTAR, thanks to the unprecedented quality of the NuSTAR data at energies >10 keV, we clearly detect, for the first time, both the primary and the reflection continuum components. The NuSTAR hard X-ray spectrum is dominated by the primary continuum piercing through an absorbing column density which is mildly optically thick to Compton scattering (tau ~ 1.2, N_H ~ 1.5 x 10^(24) cm^-2). We detect moderate hard X-ray (> 10 keV) flux variability up to 20% on short (15-20 ksec) timescales. The amplitude of the variability is maximum at ~30 keV and is likely to originate from the primary continuum of the southern nucleus. Nevertheless, the mean hard X-ray flux on longer timescales (years) is relatively constant. Moreover, the two nuclei remain Compton-thick, although we find evidence of variability of the material along the line of sight with column densities N_H <~ 2 x 10^(23) cm-2 over long (~3-15 years) timescales. The observed X-ray emission in the NuSTAR energy range is fully consistent with the sum of the best-fit models of the spatially resolved Chandra spectra of the two nuclei.

We present two NuSTAR observations of the local Seyfert 2 active galactic nucleus (AGN) and an ultraluminous X-ray source (ULX) candidate in NGC 5643. Together with archival data from Chandra, XMM-Newton and Swift-BAT, we perform a high-quality broadband spectral analysis of the AGN over two decades in energy ($\sim$0.5-100 keV). Previous X-ray observations suggested that the AGN is obscured by a Compton-thick (CT) column of obscuring gas along our line-of-sight. However, the lack of high-quality $\gtrsim$ 10 keV observations, together with the presence of a nearby X-ray luminous source, NGC 5643 X-1, had left significant uncertainties in the characterization of the nuclear spectrum. NuSTAR now enables the AGN and NGC 5643 X-1 to be separately resolved above 10 keV for the first time and allows a direct measurement of the absorbing column density toward the nucleus. The new data show that the nucleus is indeed obscured by a CT column of $N_{\rm{H}}$ $\gtrsim$ 5 $\times$ 10$^{24}$ cm$^{-2}$. The range of 2-10 keV absorption-corrected luminosity inferred from the best fitting models is $L_{2-10,\rm{int}} =$ (0.8-1.7) $\times$ 10$^{42}$ erg s$^{-1}$, consistent with that predicted from multiwavelength intrinsic luminosity indicators. We also study the NuSTAR data for NGC 5643 X-1, and show that it exhibits evidence for a spectral cut-off at energy, $E$ $\sim$ 10 keV, similar to that seen in other ULXs observed by NuSTAR. Along with the evidence for significant X-ray luminosity variations in the 3-8 keV band from 2003-2014, our results further strengthen the ULX classification of NGC 5643 X-1.

We present the best sensitivity and angular resolution maps of the molecular disk and outflow of Mrk 231, as traced by CO observations obtained with IRAM/PdBI, and we analyze archival Chandra and NuSTAR observations. We constrain the physical properties of both the molecular disk and outflow, the presence of a highly-ionized ultra-fast nuclear wind, and their connection. The molecular outflow has a size of ~1 kpc, and extends in all directions around the nucleus, being more prominent along the south-west to north-east direction, suggesting a wide-angle biconical geometry. The maximum projected velocity of the outflow is nearly constant out to ~1 kpc, thus implying that the density of the outflowing material decreases from the nucleus outwards as $r^{-2}$. This suggests that either a large part of the gas leaves the flow during its expansion or that the bulk of the outflow has not yet reached out to ~1 kpc, thus implying a limit on its age of ~1 Myr. We find $\dot M_{OF}=[ 500-1000]~ M_{\odot}~yr^{-1}$ and $\dot E_{kin,OF}=[7-10]\times 10^{43}$ erg s$^{-1}$. Remarkably, our analysis of the X-ray data reveals a nuclear ultra-fast outflow (UFO) with velocity -20000 km s$^{-1}$, $\dot M_{UFO}=[0.3- 2.1] ~M_\odot yr^{-1}$, and momentum load $\dot P_{UFO}/\dot P_{rad}=[0.2-1.6]$.We find $\dot E_{kin,UFO}\sim \dot E_{kin,OF}$ as predicted for outflows undergoing an energy conserving expansion. This suggests that most of the UFO kinetic energy is transferred to mechanical energy of the kpc-scale outflow, strongly supporting that the energy released during accretion of matter onto super-massive black holes is the ultimate driver of giant massive outflows. We estimate a momentum boost $\dot P_{OF}/\dot P_{UFO}\approx [30-60]$. The ratios $\dot E_{kin, UFO}/L_{bol,AGN} =[ 1-5]\%$ and $\dot E_{kin,OF}/L_{bol,AGN} = [1-3]\%$ agree with the requirements of the most popular models of AGN feedback.

The intrinsic column density (NH) distribution of quasars is poorly known. At the high obscuration end of the quasar population and for redshifts z<1, the X-ray spectra can only be reliably characterized using broad-band measurements which extend to energies above 10 keV. Using the hard X-ray observatory NuSTAR, along with archival Chandra and XMM-Newton data, we study the broad-band X-ray spectra of nine optically selected (from the SDSS), candidate Compton-thick (NH > 1.5e24 cm^-2) type 2 quasars (CTQSO2s); five new NuSTAR observations are reported herein, and four have been previously published. The candidate CTQSO2s lie at z<0.5, have observed [OIII] luminosities in the range 8.4 < log (L_[OIII]/L_solar) < 9.6, and show evidence for extreme, Compton-thick absorption when indirect absorption diagnostics are considered. Amongst the nine candidate CTQSO2s, five are detected by NuSTAR in the high energy (8-24 keV) band: two are weakly detected at the ~ 3 sigma confidence level and three are strongly detected with sufficient counts for spectral modeling (>~ 90 net source counts at 8-24 keV). For these NuSTAR-detected sources direct (i.e., X-ray spectral) constraints on the intrinsic AGN properties are feasible, and we measure column densities ~2.5-1600 times higher and intrinsic (unabsorbed) X-ray luminosities ~10-70 times higher than pre-NuSTAR constraints from Chandra and XMM-Newton. Assuming the NuSTAR-detected type 2 quasars are representative of other Compton-thick candidates, we make a correction to the NH distribution for optically selected type 2 quasars as measured by Chandra and XMM-Newton for 39 objects. With this approach, we predict a Compton-thick fraction of f_CT = 36^{+14}_{-12} %, although higher fractions (up to 76%) are possible if indirect absorption diagnostics are assumed to be reliable.

We present the first X-ray spectrum of a Hot dust-obscured galaxy (DOG), namely W1835+4355 at z ~ 2.3. Hot DOGs represent a very rare population of hyperluminous (>= 10^47 erg/s), dust-enshrouded objects at z > 2 recently discovered in the WISE All Sky Survey. The 40 ks XMM-Newton spectrum reveals a continuum as flat (Gamma ~ 0.8) as typically seen in heavily obscured AGN. This, along with the presence of strong Fe Kalpha emission, clearly suggests a reflection-dominated spectrum due to Compton-thick absorption. In this scenario, the observed luminosity of L(2-10 keV) ~ 2 x 10^44 erg/s is a fraction (<10%) of the intrinsic one, which is estimated to be >~ 5 x 10^45 erg/s by using several proxies. The Herschel data allow us to constrain the SED up to the sub-mm band, providing a reliable estimate of the quasar contribution (~ 75%) to the IR luminosity as well as the amount of star formation (~ 2100 Msun/yr). Our results thus provide additional pieces of evidence that associate Hot DOGs with an exceptionally dusty phase during which luminous quasars and massive galaxies co-evolve and a very efficient and powerful AGN-driven feedback mechanism is predicted by models.

Understanding the relationship between the formation and evolution of galaxies and their central super massive black holes (SMBH) is one of the main topics in extragalactic astrophysics. Links and feedback may reciprocally affect both black hole and galaxy growth. Observations of the CO line at redshifts of 2-4 are crucial to investigate the gas mass, star formation activity and accretion onto SMBHs, as well as the effect of AGN feedback. Potential correlations between AGN and host galaxy properties can be highlighted by observing extreme objects. Despite their luminosity, hyper-luminous QSOs at z=2-4 are still little studied at mm wavelengths. We targeted CO(3-2) in ULAS J1539+0557, an hyper-luminos QSO (Lbol> 10^48 erg/s) at z=2.658, selected through its unusual red colors in the UKIDSS Large Area Survey (ULAS). We find a molecular gas mass of 4.1+-0.8 10^10 Msun, and a gas fraction of 0.4-0.1, depending mostly on the assumed source inclination. We also find a robust lower limit to the star-formation rate (SFR=250-1600 Msun/yr) and star-formation efficiency (SFE=25-350 Lsun/(K km s-1 pc2) by comparing the observed optical-near-infrared spectral energy distribution with AGN and galaxy templates. The black hole gas consumption timescale, M(H_2)/dM(accretion)/dt, is ~160 Myr, similar or higher than the gas consumption timescale. The gas content and the star formation efficiency are similar to those of other high-luminosity, highly obscured QSOs, and at the lower end of the star-formation efficiency of unobscured QSOs, in line with predictions from AGN-galaxy co-evolutionary scenarios. Further measurements of the (sub)-mm continuum in this and similar sources are mandatory to obtain a robust observational picture of the AGN evolutionary sequence.

We present the first results from our pilot 500 ks Chandra-LETG Large Program observation of the soft X-ray brightest source in the z>=0.4 sky, the blazar 1ES 1553+113, aimed to secure the first uncontroversial detections of the missing baryons in the X-rays. We identify a total of 11 possible absorption lines, with single-line statistical significances between 2.2-4.1\sigma. Six of these lines are detected at high single-line statistical significance (3.6 <= sigma <= 4.1), while the remaining five are regarded as marginal detections in association with either other X-ray lines detected at higher significance and/or Far-Ultraviolet (FUV) signposts. In particular, five of these possible intervening absorption lines, are identified as CV and CVI K\alpha absorbers belonging to three WHIM systems at z_X = 0.312, z_X = 0.237 and <z_X> = 0.133, which also produce broad HI (and OVI for the z_X = 0.312 system) absorption in the FUV. For two of these systems (z_X = 0.312 and 0.237), the Chandra X-ray data led the a-posteriori discovery of physically consistent broad HI associations in the FUV, so confirming the power of the X-ray-FUV synergy for WHIM studies. The true statistical significances of these three X-ray absorption systems, after properly accounting for the number of redshift trials, are 5.8 sigma (z_X = 0.312; 6.3 sigma if the low-significance OV and CV K-beta associations are considered), 3.9 sigma (z_X = 0.237), and 3.8 sigma (\langle z_X \rangle = 0.133), respectively.

We report a detection of an absorption line at ~44.8 {\AA} in a > 500 ks Chandra HRC-S/LETG X-ray grating spectrum of the blazar H 2356-309. This line can be identified as intervening CV-K{\alpha} absorption, at z\approx0.112, produced by a warm (log T = 5.1 K) intergalactic absorber. The feature is significant at a 2.9{\sigma} level (accounting for the number of independent redshift trials). We estimate an equivalent hydrogen column density of log N_H=19.05 (Z/Zsun)^-1 cm^-2. Unlike other previously reported FUV/X-ray metal detections of warm-hot intergalactic medium (WHIM), this CV absorber lies in a region with locally low galaxy density, at ~2.2 Mpc from the closest galaxy at that redshift, and therefore is unlikely to be associated with an extended galactic halo. We instead tentatively identify this absorber with an intervening Warm-Hot Intergalactic Medium filament possibly permeating a large-scale, 30 Mpc extended, structure of galaxies whose redshift centroid, within a cylinder of 7.5 Mpc radius centered on the line of sight to H 2356-309, is marginally consistent (at a 1.8{\sigma} level) with the redshift of the absorber.

We make use of a 500ks Chandra HRC-S/LETG spectrum of the blazar H2356-309, combined with a lower S/N spectrum of the same target, to search for the presence of warm-hot absorbing gas associated with two Large-Scale Structures (LSSs) crossed by this sightline at z=0.062 (the Pisces-Cetus Supercluster, PCS) and at z=0.128 ("Farther Sculptor Wall", FSW). No statistically significant (>=3sigma) individual absorption is detected from any of the strong He- or H-like transitions of C, O and Ne at the redshifts of the structures. However we are still able to constrain the physical and geometrical parameters of the associated putative absorbing gas, by performing joint spectral fit of marginal detections and upper limits of the strongest expected lines with our self-consistent hybrid ionization WHIM spectral model. At the redshift of the PCS we identify a warm phase with logT=5.35_-0.13^+0.07 K and log N_H =19.1+/-0.2 cm^-2 possibly coexisting with a hotter and less significant phase with logT=6.9^+0.1_-0.8 K and log N_H=20.1^+0.3_-1.7 cm^-2 (1sigma errors). For the FSW we estimate logT=6.6_-0.2^+0.1 K and log N_H=19.8_-0.8^+0.4 cm^-2. Our constraints allow us to estimate the cumulative number density per unit redshifts of OVII WHIM absorbers. We also estimate the cosmological mass density obtaining Omega_b(WHIM)=(0.021^+0.031_-0.018) (Z/Z_sun)^-1, consistent with the mass density of the intergalactic 'missing baryons' for high metallicities.

We present a possible OVIII X-ray absorption line at $z=0.117 \pm 0.001$ which, if confirmed, will be the first one associated with a broad HI Ly$\beta$ (BLB: FWHM=$160^{+50}_{-30}$ km s$^{-1}$) absorber. The absorber lies along the line of sight to the nearby ($z=0.1372$) Seyfert 1 galaxy PKS~0558-504, consistent with being a WHIM filament. The X-ray absorber is marginally detected in two independent XMM-Newton spectra of PKS~0558-504, a long $\sim 600$ ks Guest-Observer observation and a shorter, $\sim 300$ ks total, calibration observation, with a combined single line statistical significance of 2.8$\sigma$ (2.7$\sigma$ and 1.2$\sigma$ in the two spectra, respectively). When fitted with our self-consistent hybrid-photoionization WHIM models, the combined XMM-{\em Newton} spectrum is consistent with the presence of OVIII K$\alpha$ at $z=(0.117 \pm 0.001)$. This model gives best fitting temperature and equivalent H column density of the absorber of log$T=6.56_{-0.17}^{+0.19}$ K, and logN$_H=(21.5 \pm 0.3) (Z/Z_{0.01\odot})^{-1}$ cm$^{-2}$. The statistical sigificance of this single X-ray detection is increased by the detection of broad and complex HI Ly$\beta$ absorption in archival FUSE spectra of PKS~0558-504, at redshifts $z=0.1183 \pm 0.0001$ consistent with the best-fitting redshift of the X-ray absorber. The single line statistical significance of this line is 4.1$\sigma$ (3.7$\sigma$ if systematics are considered), and thus the combined (HI+OVIII) statistical significance of the detection is of 5.0$\sigma$. The detection of both metal and H lines at a consistent redshift, in this hot absorbing system, allows us to speculate on its metallicity. By associating the bulk of the X-ray absorber with the BLB line detected in the FUSE spectrum at $z_{BLB}=0.1183 \pm 0.0001$, we obtain a metallicity of 1-4\% Solar.

In a previous paper we reported a 3-sigma detection of an absorption line from the Warm-Hot Intergalactic Medium (WHIM) using the Chandra and XMM X-ray grating spectra of the blazar H2356-309, the sight-line of which intercepts the Sculptor Wall, a large-scale superstructure of galaxies at z ~ 0.03. To verify our initial detection, we obtained a deep (500 ks), follow-up exposure of H2356-309 as part of the Cycle-10 Chandra Large Project Program. From a joint analysis of the Cycle-10 and previous (Cycle-8) Chandra grating data we detect the redshifted OVII WHIM line at a significance level of 3.4-sigma, a substantial improvement over the 1.7-sigma level reported previously when using only the Cycle-8 data. The significance increases to 4.0-sigma when the existing XMM grating data are included in the analysis, thus confirming at higher significance the existence of the line at the redshift of the Sculptor Wall with an equivalent width of 28.5+/-10.5 mA (90% confidence). We obtain a 90% lower limit on the OVII column density of 0.8 10^16 cm^-2 and a 90% upper limit on the Doppler-b parameter of 460 km/s. Assuming the absorber is uniformly distributed throughout the ~ 15 Mpc portion of the blazar's sight-line that intercepts the Sculptor Wall, that the OVII column density is ~ 2 10^16 cm^-2 (corresponding to b > 150 km/s where the inferred column density is only weakly dependent on b), and that the oxygen abundance is 0.1 solar, we estimate a baryon over-density of ~ 30 for the WHIM, which is consistent with the peak of the WHIM mass fraction predicted by cosmological simulations. The clear detection of OVII absorption in the Sculptor Wall demonstrates the viability of using current observatories to study WHIM in the X-ray absorption spectra of blazars behind known large-scale structures.

We present XMM RGS and Chandra LETG observations of the blazar, H 2356-309, located behind the Sculptor Wall, a large-scale galaxy structure expected to harbor high-density Warm-Hot Intergalactic Medium (WHIM). Our simultaneous analysis of the RGS and LETG spectra yields a 3-sigma detection of the crucial redshifted O vii K-alpha line with a column density (>~ 10^{16} cm^{-2}) consistent with similar large-scale structures produced in cosmological simulations. This represents the first detection of non-local WHIM from X-ray absorption studies where XMM and Chandra data are analyzed simultaneously and the absorber redshift is already known, thus providing robust evidence for the expected repository of the "missing baryons".

We present radial mass profiles within 0.3 r_vir for 16 relaxed galaxy groups-poor clusters (kT range 1-3 keV) selected for optimal mass constraints from the Chandra and XMM data archives. After accounting for the mass of hot gas, the resulting mass profiles are described well by a two-component model consisting of dark matter (DM), represented by an NFW model, and stars from the central galaxy. The stellar component is required only for 8 systems, for which reasonable stellar mass-to-light ratios (M/L_K) are obtained, assuming a Kroupa IMF. Modifying the NFW dark matter halo by adiabatic contraction does not improve the fit and yields systematically lower M/L_K. In contrast to previous results for massive clusters, we find that the NFW concentration parameter (c_vir) for groups decreases with increasing M_vir and is inconsistent with no variation at the 3 sigma level. The normalization and slope of the c_vir-M_vir relation are consistent with the standard LambdaCDM cosmological model with sigma_8 = 0.9. The small intrinsic scatter measured about the c_vir-M_vir relation implies the groups represent preferentially relaxed, early forming systems. The mean gas fraction (f =0.05 +/- 0.01) of the groups measured within an overdensity Delta=2500 is lower than for hot, massive clusters, but the fractional scatter (sigma_f/f=0.2) for groups is larger, implying a greater impact of feedback processes on groups, as expected.

We present the concentration (c)-virial mass (M) relation of 39 galaxy systems ranging in mass from individual early-type galaxies up to the most massive galaxy clusters, (0.06-20) x 10^{14} M_sun. We selected for analysis the most relaxed systems possessing the highest quality data currently available in the Chandra and XMM public data archives. A power-law model fitted to the X-ray c-M relation requires at high significance (6.6 sigma) that c decreases with increasing M, which is a general feature of CDM models. The median and scatter of the c-M relation produced by the flat, concordance LCDM model (Omega_m=0.3, sigma_8=0.9) agrees with the X-ray data provided the sample is comprised of the most relaxed, early forming systems, which is consistent with our selection criteria. Holding the rest of the cosmological parameters fixed to those in the concordance model the c-M relation requires 0.76< sigma_8 <1.07 (99% conf.), assuming a 10% upward bias in the concentrations for early forming systems. The tilted, low-sigma_8 model suggested by a new WMAP analysis is rejected at 99.99% confidence, but a model with the same tilt and normalization can be reconciled with the X-ray data by increasing the dark energy equation of state parameter to w ~ -0.8. When imposing the additional constraint of the tight relation between sigma_8 and Omega_m from studies of cluster abundances, the X-ray c-M relation excludes (>99% conf.) both open CDM models and flat CDM models with Omega_m ~1. This result provides novel evidence for a flat, low-Omega_m universe with dark energy using observations only in the local (z << 1) universe. Possible systematic errors in the X-ray mass measurements of a magnitude ~10% suggested by CDM simulations do not change our conclusions.

We present the results of spectroscopic observations of galaxies associated with the diffuse X-ray emitting structure discovered by Zappacosta et al. (2002). After measuring the redshifts of 161 galaxies, we confirm an overdensity of galaxies with projected dimensions of at least 2 Mpc, determine its spectroscopic redshift in z=0.401+/-0.002, and show that it is spatially coincident with the diffuse X-ray emission. This confirms the original claim that this X-ray emission has an extragalactic nature and is due to the Warm-Hot Intergalactic Medium (WHIM). We used this value of the redshift to compute the temperature of the emitting gas. The resulting value depends on the metallicity that is assumed for the IGM, and is constrained to be between 0.3 and 0.6 keV for metallicities between 0.05 and 0.3 solar, in good agreement with the expectations from the WHIM.

We report the discovery of a galaxy cluster serendipitously detected as an extended X-ray source in an offset observation of the group NGC 5044. The cluster redshift, z=0.281, determined from the optical spectrum of the brightest cluster galaxy, agrees with that inferred from the X-ray spectrum using the Fe K alpha complex of the hot ICM (z=0.27 +/- 0.01). Based on the 50 ks XMM observation, we find that within a radius of 383 kpc the cluster has an unabsorbed X-ray flux, f_X (0.5-2 keV) = 3.34 (+0.08, -0.13) x 10^{-13} erg/cm^2/s, a bolometric X-ray luminosity, L_X = 2.21 (+0.34, -0.19) x 10^{44} erg/s, kT = 3.57 +/- 0.12 keV, and metallicity, 0.60 +/- 0.09 solar. The cluster obeys the scaling relations for L_X and T observed at intermediate redshift. The mass derived from an isothermal NFW model fit is, M_vir = 3.89 +/- 0.35 x 10^{14} solar masses, with a concentration parameter, c = 6.7 +/- 0.4, consistent with the range of values expected in the concordance cosmological model for relaxed clusters. The optical properties suggest this could be a ``fossil cluster''.

The enormous range of X-ray luminosities among elliptical galaxies with similar optical luminosities can be attributed to large scatter and systematic variations in the dark halo mass determined from X-ray observations. The mean halo mass decreases sharply with both the X-ray and K-band luminosities. Smaller halos contain less diffuse hot gas and have lower X-ray luminosities. In addition the hot gas in low mass halos is more easily depleted by supernova-driven outflows, further lowering the X-ray luminosity.

We present an X-ray analysis of the radial mass profile of the radio-quiet galaxy cluster A2589 between 0.015-0.25 r_vir using an XMM-Newton observation. Except for a ~16 kpc shift of the X-ray center of the R=45-60 kpc annulus, A2589 possesses a remarkably symmetrical X-ray image and is therefore an exceptional candidate for precision studies of its mass profile by applying hydrostatic equilibrium. The total gravitating matter profile is well described by the NFW model (fractional residuals <~10%) with c_vir=6.1 +/- 0.3 and M_vir = 3.3 +/- 0.3 x 10^{14} M_sun (r_vir = 1.74 +/- 0.05 Mpc) in excellent agreement with LCDM. When the mass of the hot ICM is subtracted from the gravitating matter profile, the NFW model fitted to the resulting dark matter (DM) profile produces essentially the same result. However, if a component accounting for the stellar mass (M_*) of the cD galaxy is included, then the NFW fit to the DM profile is substantially degraded in the central r ~50 kpc for reasonable M_*/L_V. Modifying the NFW DM halo by adiabatic contraction arising from the early condensation of stellar baryons in the cD galaxy further degrades the fit. The fit is improved substantially with a Sersic-like model recently suggested by high resolution N-body simulations but with an inverse Sersic index, alpha ~0.5, a factor of ~3 higher than predicted. We argue that neither random turbulent motions nor magnetic fields can provide sufficient non-thermal pressure support to reconcile the XMM mass profile with adiabatic contraction of a CDM halo assuming reasonable M_*/L_V. Our results support the scenario where, at least for galaxy clusters, processes during halo formation counteract adiabatic contraction so that the total gravitating mass in the core approximately follows the NFW profile.

We present a Chandra study of mass profiles in 7 elliptical galaxies, of which 3 have galaxy-scale and 4 group-scale halos demarcated at 1E13Msun. These represent the best available data for nearby objects with comparable X-ray luminosities. We measure ~flat mass-to-light (M/L) profiles within an optical half-light radius (Reff), rising by an order of magnitude at ~10Reff, which confirms the presence of dark matter (DM). The data indicate hydrostatic equilibrium, which is also supported by agreement with studies of stellar kinematics in elliptical galaxies. The data are well-fitted by a model comprising an NFW DM profile and a baryonic component following the optical light. The distribution of DM halo concentration parameters (c) versus Mvir agrees with LCDM predictions and our observations of bright groups. Concentrations are slightly higher than expected, which is most likely a selection effect. Omitting the stellar mass drastically increases c, possibly explaining large concentrations found by some past observers. The stellar M/LK agree with population synthesis models, assuming a Kroupa IMF. Allowing adiabatic compression (AC) of the DM halo by baryons made M/L more discrepant, casting some doubt on AC. Our best-fitting models imply total baryon fractions \~0.04--0.09, consistent with models of galaxy formation incorporating strong feedback. The groups exhibit positive temperature gradients, consistent with the "Universal" profiles found in other groups and clusters, whereas the galaxies have negative gradients, suggesting a change in the evolutionary history of the systems around Mvir=1E13 Msun.

The Chandra and XMM data for a sample of 19 relaxed groups/poor clusters, covering the temperature range 1-3 keV and selected to have the best available data, reveal a remarkable similarity in their temperature profile: cool cores outside of which the temperatures reach a peak for radii less than 0.1 of the virial radius and then decline. We fitted the derived mass profiles using an NFW model, which provides a good fit to the data when accounting for the central galaxy in the inner region. The concentration parameters and virial masses are in the range c = 5-22 and M = 2 x 10^13 - 4 x 10^14 solar masses, in general agreement with the concentrations found in numerical simulations.

Cosmological simulations of galaxy formation predict a universal form for the mass profile of dark matter (DM) halos from cluster to galaxy scales. Remarkably few interesting constraints exist, however, on DM halos in early-type galaxies. Using Chandra we present the temperature, density and mass profiles of a small sample of early-type galaxies, revealing significant DM in each case. When a component is included to account for stellar mass and the DM halo is allowed to respond adiabatically to the baryonic condensation into stars, the mass profiles are well-fitted by the universal profile, with Virial masses and concentrations in agreement with simulations. However, only ~half, or less, of the mass within R_e seems attributable to the stars, implying stellar M_*/L_B ~1-5.

We have analyzed an XMM-Newton observation of the cluster Abell 2589. Apart from a low-level asymmetry in the central region, the cluster appears very relaxed and does not show presence of a central AGN. We derived constraints for the radial temperature, density and, assuming hydrostatic equilibrium, mass profiles. We find that the best fit to the dark matter profile is given by the Sersic-like profile proposed by Navarro et al. (2004). The NFW model does not provide a good fit. We also tested whether the central stellar component could affect the profile through the adiabatic contraction model but were unable to distinguish it from a simple "dark matter + stars" modeling.

Q0056-363 is the most powerful X-ray quasar known to exhibit a broad, likely relativistic iron line (Porquet & Reeves 2003). It has been observed twice by XMM-$Newton$, three and half years apart (July 2000 and December 2003). In the second observation, the UV and soft X-ray fluxes were fainter, the hard X-ray power law flatter, and the iron line equivalent width (EW) smaller than in the 2000 observation. These variations can all be explained, at least qualitatively, if the disc is truncated in the second observation. We report also on the possible detection of a transient, redshifted iron absorption line during the 2003 observation.

We have identified a large-scale structure traced by galaxies at z=0.8, within the Lockman Hole, by means of multi-object spectroscopic observations. By using deep XMM images we have investigated the soft X-ray emission from the Warm-Hot Intergalactic Medium (WHIM) expected to be associated with this large-scale structure and we set a tight upper limit to its flux in the very soft 0.2-0.4 keV band. The non-detection requires the WHIM at these redshifts to be cooler than 0.1 keV. Combined with the WHIM emission detections at lower redshift, our result indicates that the WHIM temperature is rapidly decreasing with redshift, as expected in popular cosmological models.

We have analyzed the soft X-ray emission in a wide area of the Sculptor supercluster by using overlapping ROSAT PSPC pointings. After subtraction of the point sources we have found evidence for extended, diffuse soft X-ray emission. We have investigated the nature of such extended emission through the cross-correlation with the density of galaxies as inferred from the Muenster Redshift Survey. In particular we have analyzed the correlation as a function of the temperature of the X-ray emitting gas. We have found a significant correlation of the galaxy distribution only with the softest X-ray emission (0.1--0.3 keV) and only for gas temperatures kT < 0.5 keV. We have excluded that this soft X-ray diffuse emission, and its correlation with the galaxy distribution, is significantly contributed by unresolved AGN, group of galaxies or individual galaxies. The most likely explanation is that the soft, diffuse X-ray emission is tracing Warm-Hot Intergalactic Medium, with temperatures below 0.5 keV, associated with the large-scale structures in the Sculptor supercluster.

Several popular cosmological models predict that most of the baryonic mass in the local universe is located in filamentary and sheet-like structures associated with groups and clusters of galaxies. This gas is expected to be gravitationally heated to ~10^6 K and therefore emitting in the soft X-rays. We have investigated three fields with large scale structures of galaxies at redshifts 0.1, 0.45, 0.79 and found signatures of warm-hot thermal emission (kT< 1 keV) correlated with the distribution of galaxies for the first two. The correlation and the properties of both X-ray and galaxy distribution strongly suggest that the diffuse X-ray flux is due to extragalactic emission by the Warm-Hot Intergalactic Medium (WHIM) predicted by cosmological models.

Several popular cosmological models predict that most of the baryonic mass in the local universe is located in filamentary and sheet-like structures associated with galaxy overdensities. This gas is expected to be gravitationally heated to ~10^6 K and therefore emitting in the soft X-rays. We have detected diffuse soft X-ray structures in a high Galactic latitude ROSAT field after point source subtraction and correction for Galactic absorption. These diffuse structures have an X-ray energy distribution that is much softer than expected from clusters, groups or unresolved emission from AGNs, but are consistent with that expected from a diffuse warm intergalactic medium. To discriminate between a Galactic or extragalactic nature of the diffuse gas we have correlated the soft X--map with multiband optical images in this field. We have found a significant overdensity of galaxies in correspondence with the strongest diffuse X-ray structure. The photometric redshift distribution of the galaxies over the X-ray peak has an excess over field galaxies at z~0.45. This result strongly suggests that the diffuse X-ray flux is due to extragalactic emission by warm gas associated with an overdense galaxy region at z~0.45.