Quasars are the most luminous of active galactic nuclei (AGN), and are perhaps responsible for quenching star formation in their hosts. The Stripe 82X catalog covers 31.3 deg$^2$ of the Stripe 82 field, of which the 15.6 deg$^2$ covered with XMM-Newton is also covered by Herschel/SPIRE. We have 2500 X-ray detected sources with multi-wavelength counterparts, and 30% of these are unobscured quasars, with $L_X > 10^{44}\,$erg/s and $M_B < -23$. We define a new population of quasars which are unobscured, have X-ray luminosities in excess of $10^{44}\,$erg/s, have broad emission lines, and yet are also bright in the far-infrared, with a 250$\mu$m flux density of $S_{\rm 250}>30$mJy. We refer to these Herschel-detected, unobscured quasars as "Cold Quasars". A mere 4% (21) of the X-ray- and optically-selected unobscured quasars in Stripe 82X are detected at 250$\mu$m. These Cold Quasars lie at $z\sim1-3$, have $L_{\rm IR}>10^{12}\,L_\odot$, and have star formation rates of $\sim200-1400\,M_\odot$/yr. Cold Quasars are bluer in the mid-IR than the full quasar population, and 72% of our Cold Quasars have WISE W3 $<$ 11.5 [Vega], while only 19% of the full quasar sample meets this criteria. Crucially, Cold Quasars have on average $\sim9\times$ as much star formation as the main sequence of star forming galaxies at similar redshifts. Although dust-rich, unobscured quasars have occasionally been noted in the literature before, we argue that they should be considered as a separate class of quasars due to their high star formation rates. This phase is likely short-lived, as the central engine and immense star formation consume the gas reservoir. Cold Quasars are type-1 blue quasars that reside in starburst galaxies.

We report on the investigation of a very high energy (VHE), Galactic gamma-ray source recently discovered at >50GeV using the Large Area Telescope (LAT) on board the Fermi Gamma-Ray Space Telescope. This object, 2FHL J1703.4-4145, displays a very hard >50GeV spectrum with a photon index ~1.2 in the 2FHL catalog and, as such, is one of the most extreme sources in the 2FHL sub-sample of Galactic objects. A detailed analysis of the available multi-wavelength data shows that this source is located on the western edge of the supernova remnant (SNR) G344.7--0.1, along with extended TeV source, HESS J1702-420. The observations and the spectral energy distribution modeling support a scenario where this gamma-ray source is the byproduct of the interaction between the SNR shock and the dense surrounding medium, with escaping cosmic rays (CRs) diffusing into the dense environment and interacting with a large local cloud, generating the observed TeV emission. If confirmed, an interaction between the SNR CRs and a nearby cloud would make 2FHL J1703.4-4145 another promising candidate for efficient particle acceleration of the 2FHL Galactic sample, following the first candidate from our previous investigation of a likely shock-cloud interaction occurring on the West edge of the Vela SNR.

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.

We report a detailed study of an unidentified gamma-ray source located in the region of the compact stellar cluster NGC 3603. This is a star-forming region (SFR) powered by a massive cluster of OB stars. A dedicated analysis of about 10 years of data from 10 GeV - 1 TeV, provided by the Large Area Telescope (LAT) onboard the Fermi Gamma-ray Space Telescope, yields the detection of a pointlike source at a significance of 9$\sigma$. The source photon spectrum can be described by a power-law model with best fit spectral index of $2.35 \pm 0.03$. In addition, the analysis of a deep Chandra image in the 0.5 - 7 keV band reliably rules out an extragalactic origin for the gamma rays. We also conclude that the broadband spectral energy distribution of the point source can be explained well with both leptonic and hadronic models. No firm evidence of association with any other classes of known gamma-ray emitters is found, therefore we speculate that 4FGL J1115.1-6118 is a gamma-ray emitting SFR.

In modeling the X-ray spectra of active galactic nuclei (AGNs), the inclination angle is a parameter that can play an important role in analyzing the X-ray spectra of AGN, but it has never been studied in detail. We present a broadband X-ray spectral analysis of the joint NuSTAR-XMM-Newton observations of 13 sources with [OIII] measured inclinations determined by Fischer et al. 2013. By freezing the inclination angles at the [OIII] measured values when modeling the observations, the spectra are well fitted and the geometrical properties of the obscuring structure of the AGNs are slightly better constrained than those fitted when the inclination angles are left free to vary. We also test if one could freeze the inclinations at other specific angles in fitting the AGN X-ray spectra as commonly did in the literatures. We find that one should always let the inclination angle free to vary in modeling the X-ray spectra of AGNs, while fixing the inclination angle at [OIII] measured values and fixing the inclination angle at 60$^\circ$ also present correct fits of the sources in our sample.Correlations between the covering factor and the average column density of the obscuring torus with respect to the Eddington ratio are also measured, suggesting that the distribution of the material in the obscuring torus is regulated by the Eddington ratio, which is in agreement with previous studies. In addition, no geometrical correlation is found between the narrow line region of the AGN and the obscuring torus, suggesting that the geometry might be more complex than what is assumed in the simplistic unified model.

What triggers collimated relativistic outflows or jets, from the centers of galaxies remains a fundamental question in astrophysics. The merging of two galaxies has been proposed to realize the conditions to successfully launch and drive such jets into the intergalactic medium. However, evidences for the operation of this mechanism are scarce. Here we report the first unambiguous detection of an ongoing merger of a narrow-line Seyfert 1 galaxy, TXS 2116$-$077, hosting a closely aligned, $\gamma$-ray emitting relativistic jet with a Seyfert 2 galaxy at a separation of $\sim$12 kpc, using the observations taken with 8.2 m Subaru telescope. Our subsequent followup observations with 10.4 m Gran Telescopio Canarias, 4.2 m William Herschel Telescope, and Chandra X-ray observatory have provided what is likely to be the first glimpse of the merging environment hosting a closely aligned relativistic jet. Our finding that the jet is considerably younger than the merger demonstrates that jet activity can be triggered by galaxy mergers and that $\gamma$-ray detected narrow-line Seyfert 1 galaxies represent the beginning phase of that activity. These results also highlight the crucial role of mergers in shaping the fate of galaxies in their cosmological evolution and are consistent with recent studies focused on the host galaxy imaging of this enigmatic class of active galactic nuclei.

We present the X-ray source catalog for the 479 ks Chandra exposure of the SDSS J1030+0524 field, that is centered on a region that shows the best evidence to date of an overdensity around a z > 6 quasar, and also includes a galaxy overdensity around a Compton-thick Fanaroff-Riley type II radio galaxy at z = 1.7. Using wavdetect for initial source detection and ACIS Extract for source photometry and significance assessment, we create preliminary catalogs of sources that are detected in the full, soft, and hard bands, respectively. We produce X-ray simulations that mirror our Chandra observation to filter our preliminary catalogs and get a completeness level of > 91% and a reliability level of 95% in each band. The catalogs in the three bands are then matched into a final main catalog of 256 unique sources. Among them, 244, 193, and 208 are detected in the full, soft, and hard bands, respectively. This makes J1030 field the fifth deepest extragalactic X-ray survey to date. The field is part of the Multiwavelength Survey by Yale-Chile (MUSYC), and is also covered by optical imaging data from the Large Binocular Camera (LBC) at the Large Binocular Telescope, near-IR imaging data from the Canada France Hawaii Telescope WIRCam, and Spitzer IRAC. Thanks to its dense multi-wavelength coverage, J1030 represents a legacy field for the study of large-scale structures around distant accreting supermassive black holes. Using a likelihood ratio analysis, we associate multi-band counterparts for 252 (98.4%) of the 256 Chandra sources, with an estimated reliability of 95%. Finally, we compute the cumulative number of sources in each X-ray band, finding that they are in general agreement with the results from the Chandra Deep Fields.

With bolometric luminosities exceeding $10^{48}$ erg s$^{-1}$, powerful jets and supermassive black holes at their center, MeV blazars are some of the most extreme sources in the Universe. Recently, the Fermi-Large Area Telescope detected five new $\gamma$-ray emitting MeV blazars beyond redshift $z=3.1$. With the goal of precisely characterizing the jet properties of these extreme sources, we started a multiwavelength campaign to follow them up with joint NuSTAR, Swift and SARA observations. We observe six high-redshift quasars, four of them belonging to the new $\gamma$-ray emitting MeV blazars. Thorough X-ray analysis reveals spectral flattening at soft X-ray for three of these objects. The source NVSS J151002$+$570243 also shows a peculiar re-hardening of the X-ray spectrum at energies $E>6\,\rm keV$. Adopting a one-zone leptonic emission model, this combination of hard X-rays and $\gamma$-rays enables us to determine the location of the Inverse Compton peak and to accurately constrain the jet characteristics. In the context of the jet-accretion disk connection, we find that all six sources have jet powers exceeding accretion disk luminosity, seemingly validating this positive correlation even beyond $z>3$. Our six sources are found to have $10^9 \rm M_{\odot}$ black holes, further raising the space density of supermassive black holes in the redshift bin $z=[3,4]$.

We constrain X-ray spectral shapes for the ensemble of AGN based on the shape of the Cosmic X-ray Background (CXB). Specifically, we rule out regions of X-ray spectral parameter space that do not reproduce the CXB in the energy range 1-100 keV. The key X-ray spectral parameters are the photon index, {\Gamma}; the cutoff energy, Ecutoff; and the reflection scaling factor, R. Assuming each parameter follows a Gaussian distribution, we first explore the parameter space using a Bayesian approach and a fixed X-ray luminosity function (XLF). For {\sigma}_E = 36 keV and {\sigma}_R = 0.14, fixed at the observed values from the Swift-BAT 70-month sample, we allow <R>, <Ecutoff > and <{\Gamma}> to vary subject to reproducing the CXB. We report results for {\sigma}_{\Gamma} = 0.1-0.5. In an alternative approach, we define the parameter distributions, then forward model to fit the CXB by perturbing the XLF using a neural network. This approach allows us to rule out parameter combinations that cannot reproduce the CXB for any XLF. The marginalized conditional probabilities for the four free parameters are: <R> = 0.99^{+0.11}_{-0.26}, <Ecutoff> = 118^{+24}_{-23}, {\sigma}_{\Gamma} = 0.101^{+0.097}_{-0.001} and <{\Gamma}> = 1.9^{+0.08}_{-0.09}. We provide an interactive online tool for users to explore any combination of <Ecutoff>, {\sigma}_E, <{\Gamma}>, {\sigma}_{\Gamma}, <R> and {\sigma}_R including different distributions for each absorption bin, subject to the integral CXB constraint. The distributions observed in many AGN samples can be ruled out by our analysis, meaning these samples can not be representative of the full AGN population. The few samples that fall within the acceptable parameter space are hard X-ray-selected, commensurate with their having fewer selection biases.

We present rest-frame far-infrared (FIR) and optical size measurements of AGN hosts and star-forming galaxies in the COSMOS field, enabled by high-resolution ALMA/1 mm (0.1 arcsec - 0.4 arcsec) and HST/F814W imaging (~ 0.1 arcsec). Our sample includes 27 galaxies at z<2.5, classified as infrared-selected AGN (3 sources), X-ray selected AGN (4 sources), and non-AGN star-forming galaxies (20 sources), for which high-resolution Band 6/7 ALMA images are available at 1 mm from our own observing program as well as archival observations. The sizes and SFR surface densities measured from both ALMA/1 mm and HST/F814W images show that obscured AGN host galaxies are more compact than non-AGN star-forming galaxies at similar redshift and stellar mass. This result suggests that the obscured accretion phase may be related to galaxies experiencing a compaction of their gaseous component, which could be associated with enhanced central star formation before a subsequent quenching driving the formation of compact passive galaxies. Moreover, most of the detected and stacked rest-frame FIR sizes of AGNs in our sample are similar or more compact than their rest-frame optical sizes, which is consistent with recent results of ALMA detected sources. This might be explained by the fact that the dusty starbursts take place in the compact regions, and suggests that the star formation mechanisms in the compact regions of AGN hosts are similar to those observed in star-forming galaxies observed with ALMA.

Aims. We perform clustering measurements of 800 X-ray selected Chandra COSMOS Legacy (CCL) Type 2 AGN with known spectroscopic redshift to probe the halo mass dependence on AGN host galaxy properties, such as galaxy stellar mass Mstar, star formation rate (SFR) and specific black hole accretion rate (BHAR), in the redshift range z = [0 - 3]. Methods. We split the sample of AGN with known spectroscopic redshits according to Mstar, SFR and specific BHAR, while matching the distributions in terms of the other parameters, including redshift. We measure the projected two-point correlation function wp(rp) and model it with the 2-halo term to derive the large-scale bias b and the corresponding typical mass of the hosting halo, for the different subsamples. Results. We found no significant dependence of the large-scale bias and typical halo mass on galaxy stellar mass and specific BHAR for CCL Type 2 AGN at mean z~1, while a negative dependence on SFR is observed, with lower SFR AGN residing in richer environment. Mock catalogs of AGN matched to have the same X-ray luminosity, stellar mass, specific BHAR and SFR of CCL Type 2 AGN, almost reproduce the observed Mstar, specific BHAR and SFR-Mh relations, when assuming a fraction of satellite AGN of 15%, which corresponds to a ratio between the probabilities of satellite and central AGN of being active Q = 2. Mock matched normal galaxies follow a slightly steeper Mstar -Mh relation with low mass mock galaxies residing in less massive halos than mock AGN of similar mass, and are less biased than mock AGN with similar specific BHAR and SFR, at least for Q > 1.

We present the joint NuSTAR and XMM-Newton spectral analysis in the 0.6-70 keV band of three candidate Compton thick (CT-) AGN selected in the 100-month Swift-BAT catalog. These objects were previously classified as CT-AGNs based on low quality Swift-XRT and Swift-BAT data, and had soft photon indices (Gamma>2.2) that suggested a potential overestimation of the line of sight column density. Thanks to the high-quality NuSTAR and XMM-Newton data we were able to determine that in all three objects the photon index was significantly overestimated, and two out of three sources are reclassified from CT to Compton thin, confirming a previously observed trend, i.e., that a significant fraction of BAT-selected, candidate CT-AGNs with poor soft X-ray data are reclassified as Compton thin when the NuSTAR data are added to the fit. Finally, thanks to both the good XMM-Newton spatial resolution and the high NuSTAR and XMM-Newton spectral quality, we found that the third object in our sample was associated to the wrong counterpart: the correct one, 2MASX J10331570+5252182, has redshift z=0.14036, which makes it one of the very few candidate CT-AGNs in the 100-month BAT catalog detected at z>0.1, and a rare CT quasar.

The cosmic X-ray background (CXB), which peaks at an energy of ~30 keV, is produced primarily by emission from accreting supermassive black holes (SMBHs). The CXB therefore serves as a constraint on the integrated SMBH growth in the Universe and the accretion physics and obscuration in active galactic nuclei (AGNs). This paper gives an overview of recent progress in understanding the high-energy (>~10 keV) X-ray emission from AGNs and the synthesis of the CXB, with an emphasis on results from NASA's NuSTAR hard X-ray mission. We then discuss remaining challenges and open questions regarding the nature of AGN obscuration and AGN physics. Finally, we highlight the exciting opportunities for a next-generation, high-resolution hard X-ray mission to achieve the long-standing goal of resolving and characterizing the vast majority of the accreting SMBHs that produce the CXB.

We present the X-ray properties of 108 Dust-Obscured Galaxies (DOGs; F$_{24 \mu m}$/F$_{R} >$ 1000) in the COSMOS field, all of which detected in at least three far-infrared bands with the Herschel Observatory. Out of the entire sample, 22 are individually detected in the hard 2-8 keV X-ray band by the Chandra COSMOS Legacy survey, allowing us to classify them as AGN. Of them, 6 (27%) are Compton Thick AGN candidates with column densities N$_{H}$$>$10$^{24}$ cm$^{-2}$ while 15 are moderately obscured AGNs with 10$^{22}$ $<$ N$_{H}$ $<$ 10$^{24}$ cm$^{-2}$. Additionally, we estimate AGN contributions to the IR luminosity (8-1000$\mu$m rest-frame) greater than 20% for 19 DOGs based on SED decomposition using Spitzer/MIPS 24$\mu$m and the five Herschel bands (100-500 $\mu$m). Only 7 of these are detected in X-rays individually. We performed a X-ray stacking analysis for the 86 undetected DOGs. We find that the AGN fraction in DOGs increases with 24$\mu$m flux and that it is higher than that of the general 24$\mu$m population. However, no significant difference is found when considering only X-ray detections. This strongly motivates the combined use of X-ray and far-IR surveys to successfully probe a wider population of AGNs, particularly for the most obscured ones.

The discoveries made over the past 20 years by Chandra and XMM-Newton surveys in conjunction with multiwavelength imaging and spectroscopic data available in the same fields have significantly changed the view of the supermassive black hole (SMBH) and galaxy connection. These discoveries have opened up several exciting questions that are beyond the capabilities of current X-ray telescopes and will need to be addressed by observatories in the next two decades. As new observatories peer into the early Universe, we will begin to understand the physics and demographics of SMBH infancy (at $z>6$) and investigate the influence of their accretion on the formation of the first galaxies ($\S$ 2.1). We will also be able to understand the accretion and evolution over the cosmic history (at $z\sim$1-6) of the full population of black holes in galaxies, including low accretion rate, heavily obscured AGNs at luminosities beyond the reach of current X-ray surveys ($\S$2.2 and $\S$2.3), enabling us to resolve the connection between SMBH growth and their environment.

We investigate the multi-wavelength properties of host galaxies of 3701 X-ray-selected active galactic nuclei (AGNs) out to z~5 in the Chandra-COSMOS Legacy Survey. Thanks to the extensive multi-wavelength photometry available in the COSMOS field, we derive AGN luminosities, host stellar masses, and star formation rates (SFRs) via a multi-component SED fitting technique. Type 1 and Type 2 AGNs follow the same intrinsic Lx-L6um relation, suggesting that mid-infrared emission is a reasonably good measure of the AGN accretion power regardless of obscuration. We find that there is a strong increase in Type 1 AGN fraction toward higher AGN luminosity, possibly due to the fact that Type 1 AGNs tend to be hosted by more massive galaxies. The AGN luminosity and SFR are consistent with an increase toward high stellar mass, while both the Mstellar-dependence is weaker towards the high-mass end, which could be interpreted as a consequence of quenching both star formation and AGN activity in massive galaxies. AGN host galaxies tend to have SFRs that are consistent with normal star-forming galaxies, independent of AGN luminosities. We confirm that black hole accretion rate and SFR are correlated up to z~5, when forming stars. The majority (~73%) of our AGN sample are faint in the far-infrared, implying that the moderate-luminosity AGNs seem to be still active after the star formation is suppressed. It is not certain whether AGN activity plays a role in quenching the star formation. We conclude that both AGN activity and star formation might be more fundamentally related to host stellar mass.

Active galactic Nuclei (AGNs) with their relativistic jets pointed towards the observer, form a subclass of luminous gamma-ray sources commonly known as blazars. The study of blazars is essential to improve our understanding on the AGNs emission mechanisms and evolution, as well as to map the extragalactic background light. To do so, however, one needs to correctly classify and measure a redshift for a large sample of these sources. The Third {\it Fermi}--LAT Catalog of High-Energy Sources (3FHL) contains $\approx1160$ blazars reported at energies greater than $10$\,GeV. However $\sim$25\% of these sources are unclassified and $\sim$50\% lack of redshift information. To increase the spectral completeness of the 3FHL catalog, we are working on an optical spectroscopic follow up campaign using 4--m and 8--m telescopes. In this paper, we present the results of the second part of this campaign, where we observed 23 blazars using the 4\,$m$ telescope at CTIO in Chile. We report all the 23 sources to be classified as BL Lacs, a confirmed redshift measurement for 3 sources, a redshift lower limit for 2 sources and a tentative redshift measurement for 3 sources.

In this data paper we present and characterise the multi-component radio sources identified in the VLA-COSMOS Large Project at 3 GHz (0.75 arcsec resolution, 2.3 {\mu}Jy/beam rms), i.e. the radio sources which are composed of two or more radio blobs.The classification of objects into multi-components was done by visual inspection of 351 of the brightest and most extended blobs from a sample of 10,899 blobs identified by the automatic code blobcat. For that purpose we used multi-wavelength information of the field, such as the 1.4-GHz VLA-COSMOS data and the UltraVISTA stacked mosaic available for COSMOS. We have identified 67 multi-component radio sources at 3 GHz: 58 sources with AGN powered radio emission and 9 star-forming galaxies. We report 8 new detections that were not observed by the VLA-COSMOS Large Project at 1.4 GHz, due to the slightly larger area coverage at 3 GHz. The increased spatial resolution of 0.75 arcsec has allowed us to resolve (and isolate) multiple emission peaks of 28 extended radio sources not identified in the 1.4-GHz VLA-COSMOS map. We report the multi-frequency flux densities (324 MHz, 325 MHz, 1.4 GHz & 3 GHz), star-formation-rates, and stellar masses of these objects. Multi-component objects at 3-GHz VLA-COSMOS inhabit mainly massive galaxies (>10^10.5 Msun). The majority of the multi-component AGN lie below the main-sequence of star-forming galaxies (SFGs), in the green valley and the quiescent region. We provide detailed description of the objects: amongst the AGN there are 2 head-tail, 10 core-lobe, 9 wide-angle-tail (WAT), 8 double-double or Z-/X-shaped, 3 bent-tail radio sources, and 26 symmetric sources, while amongst the SFGs we find the only star-forming ring seen in radio emission in COSMOS. We report a large number (32/58) of disturbed/bent multi-component AGN, 18 of which do not lie within X-ray groups in COSMOS (0.08 < z < 1.53). [abridged]

We present the results of an identification campaign of unassociated sources from the Fermi Large Area Telescope 3FHL catalog. Out of 200 unidentified sources, we selected 110 sources for which archival Swift-XRT observations were available, 52 of which were found to have exactly one X-ray counterpart within the 3FHL 95\% positional uncertainty. In this work, we report the X-ray, optical, IR, and radio properties of these 52 sources using positional associations with objects in various catalogs. The Wide-field Infrared Survey Explorer color-color plot for sources suggests that most of these belong to the blazar class family. The redshift measurements for these objects range from z=0.277 to z=2.1. Additionally, under the assumption that the majority of these sources are blazars, three machine-learning algorithms are employed to classify the sample into flat spectrum radio quasars or BL Lacertae objects. These suggest that the majority of the previously unassociated sources are BL Lac objects, in agreement with the fact the BL Lac objects represent by far the most numerous population detected above 10 GeV in 3FHL.

We present the analysis of a sample of 35 candidate Compton thick (CT-) active galactic nuclei (AGNs) selected in the nearby Universe (average redshift <z>~0.03) with the Swift-BAT 100-month survey. All sources have available NuSTAR data, thus allowing us to constrain with unprecedented quality important spectral parameters such as the obscuring torus line-of-sight column density (N_{H, z}), the average torus column density (N_{H, tor}) and the torus covering factor (f_c). We compare the best-fit results obtained with the widely used MyTorus (Murphy et al. 2009) model with those of the recently published borus02 model (Balokovic et al. 2018) used in the same geometrical configuration of MyTorus (i.e., with f_c=0.5). We find a remarkable agreement between the two, although with increasing dispersion in N_{H, z} moving towards higher column densities. We then use borus02 to measure f_c. High-f_c sources have, on average, smaller offset between N_{H, z} and N_{H, tor} than low-f_c ones. Therefore, low f_c values can be linked to a "patchy torus" scenario, where the AGN is seen through an over-dense region in the torus, while high-f_c objects are more likely to be obscured by a more uniform gas distribution. Finally, we find potential evidence of an inverse trend between f_c and the AGN 2-10 keV luminosity, i.e., sources with higher f_c values have on average lower luminosities.

We present the 2--78 keV spectral analysis of the deep NuSTAR and XMM-Newton observation of a nearby Seyfert 2 galaxy, ESO 116-G018, which is selected as a candidate Compton-thick (CT-) active galactic nucleus (AGN) based on a previous Chandra-Swift/BAT study. Through our analysis, the source is for the first time confirmed to be a CT-AGN at a $>$3 $\sigma$ confidence level, with the 'line-of-sight' column density N$_{\rm H,Z}$ = [2.46-2.76] $\times$ $10^{24}$ cm$^{-2}$. The 'global average' column density of the obscuring torus is N$_{\rm H,S}$ = [0.46-0.62] $\times$ $10^{24}$ cm$^{-2}$, which suggests a clumpy, rather than uniform, distribution of the obscuring material surrounding the accreting supermassive black hole. The excellent-quality data given by the combined NuSTAR and XMM-Newton observations enable us to produce a strong constraint on the covering factor of the torus of ESO 116-G018, which is found to be $f_c$ = [0.13-0.15]. We also estimate the bolometric luminosity from the broad-band X-ray spectrum, being L$\rm _{bol}$ = [2.57-3.41] $\times$ 10$^{44}$ erg s$^{-1}$.

As matter accretes onto the central supermassive black holes in active galactic nuclei (AGN), X-rays are emitted. We present a population synthesis model that accounts for the summed X-ray emission from growing black holes; modulo the efficiency of converting mass to X-rays, this is effectively a record of the accreted mass. We need this population synthesis model to reproduce observed constraints from X-ray surveys: the X-ray number counts, the observed fraction of Compton-thick AGN [$\log$ (N$_{\rm H}$/cm$^{\rm -2}$) $>$ 24] and the spectrum of the cosmic X-ray background (CXB), after accounting for selection biases. Over the past decade, X-ray surveys by {\it XMM-Newton}, {\it Chandra}, \textit{NuSTAR} and \textit{Swift}-BAT have provided greatly improved observational constraints. We find that no existing X-ray luminosity function (XLF) consistently reproduces all these observations. We take the uncertainty in AGN spectra into account, and use a neural network to compute an XLF that fits all observed constraints, including observed Compton-thick number counts and fractions. This new population synthesis model suggests that, intrinsically, 50$\pm$9\% (56$\pm$7\%) of all AGN within z $\simeq$ 0.1 (1.0) are Compton-thick.

We report on the investigation of a very high energy (VHE), Galactic gamma-ray source recently discovered at >50GeV using the Large Area Telescope (LAT) on board the Fermi-LAT. This object, 2FHL J0826.1-4500, displays one of the hardest >50GeV spectra (gamma-ray photon index ~1.6) in the 2FHL catalog, and a follow-up observation with XMM-Newton has uncovered diffuse, soft thermal emission at the position of the gamma-ray source. A detailed analysis of the available multi-wavelength data shows that this source is located on the Western edge of the Vela supernova remnant (SNR): the observations and the spectral energy distribution modelling support a scenario where this gamma-ray source is the byproduct of the interaction between the SNR shock and a neutral Hydrogen cloud. If confirmed, this shock-cloud interaction would make 2FHL J0826.1-4500 a promising candidate for efficient particle acceleration.

We present the combined NuSTATR and XMM-Newton 0.6-79 keV spectral analysis of a Seyfert 2 galaxy, NGC 1358, which we selected as a candidate Compton thick (CT-) active galactic nucleus (AGN) on the basis of previous Swift/BAT and Chandra studies. According to our analysis, NGC 1358 is confirmed to be a CT-AGN using physical motivated models, at >3 $\sigma$ confidence level. Our best-fit shows that the column density along the 'line-of-sight' of the obscuring material surrounding the accreting super-massive black hole is N$\rm _H$ = [1.96--2.80] $\times$ 10$^{24}$ cm$^{-2}$. The high-quality data from NuSTAR gives the best constraints on the spectral shape above $\sim$10 keV to date on NGC 1358. Moreover, by combining NuSTAR and XMM-Newton data, we find that the obscuring torus has a low covering factor ($f_c$ <0.17), and the obscuring material is distributed in clumps, rather than uniformly. We also derive an estimate of NGC 1358's Eddington ratio, finding it to be $\lambda_{\rm Edd}$ $\sim$$4.7_{-0.3}^{+0.3}$ $\times$ 10$^{-2}$, which is in acceptable agreement with previous measurements. Finally, we find no evidence of short-term variability, over a $\sim$100 ks time-span, in terms of both 'line-of-sight' column density and flux.

The existence of a large population of Compton thick (CT, $N_{H}>10^{24} cm^{-2}$) AGN is a key ingredient of most Cosmic X-ray background synthesis models. However, direct identification of these sources, especially at high redshift, is difficult due to the flux suppression and complex spectral shape produced by CT obscuration. We explored the Chandra COSMOS Legacy point source catalog, comprising 1855 sources, to select via X-ray spectroscopy, a large sample of CT candidates at high redshift. Adopting a physical model to reproduce the toroidal absorber, and a Monte-Carlo sampling method, we selected 67 individual sources with >5% probability of being CT, in the redshift range $0.04<z<3.5$. The sum of the probabilities above $N_{H}>10^{24} cm^{-2}$, gives a total of 41.9 effective CT, corrected for classification bias. We derive number counts in the 2-10 keV band in three redshift bins. The observed logN-logS is consistent with an increase of the intrinsic CT fraction ($f_{CT}$) from $\sim0.30$ to $\sim0.55$ from low to high redshift. When rescaled to a common luminosity (log(L$_{\rm X}$/erg/s)$=44.5$) we find an increase from $f_{CT}=0.19_{-0.06}^{+0.07}$ to $f_{CT}=0.30_{-0.08}^{+0.10}$ and $f_{CT}=0.49_{-0.11}^{+0.12}$ from low to high z. This evolution can be parametrized as $f_{CT}=0.11_{-0.04}^{+0.05}(1+z)^{1.11\pm0.13}$. Thanks to HST-ACS deep imaging, we find that the fraction of CT AGN in mergers/interacting systems increases with luminosity and redshift and is significantly higher than for non-CT AGN hosts.

Exploring how radio-emitting active galactic nuclei (AGN) behave and evolve with time is critical for understanding how AGN feedback impacts galaxy evolution. In this work, we investigate the relationship between 1.4 GHz radio continuum AGN luminosity ($L^{\rm AGN}_{\rm 1.4}$), specific black hole accretion rate (s-BHAR, defined as the accretion luminosity relative to the galaxy stellar mass) and redshift, for a luminosity-complete sample of radio-selected AGN in the VLA COSMOS 3 GHz Large Project. The sample was originally selected from radio-continuum observations at 3 GHz, and includes about 1800 radio AGN identified via ($>2\sigma$) radio-excess relative to the infrared-radio correlation of star-forming galaxies. We further select a subsample of over 1200 radio AGN that is complete in $L^{\rm AGN}_{\rm 1.4}$ over different redshift ranges, out to z~4, and use X-ray stacking to calculate the average s-BHAR in each $L^{\rm AGN}_{\rm 1.4}$-$z$ bin. We find that the average s-BHAR is independent of $L^{\rm AGN}_{\rm 1.4}$, at all redshifts. However, we see a strong increase of s-BHAR with redshift, at fixed $L^{\rm AGN}_{\rm 1.4}$. This trend resembles the strong increase in the fraction of star-forming host galaxies (based on the $(NUV-r)$ / $(r-J)$ colours) with redshift, at fixed $L^{\rm AGN}_{\rm 1.4}$. A possible explanation for this similarity might imply a link between average AGN radiative power and availability of cold gas supply within the host galaxy. This study corroborates the idea that radio-selected AGN become more radiatively efficient towards earlier epochs, independently of their radio power.

We present the results of the optical spectroscopy follow-up of a sample of 28 unclassified blazars from the Third Fermi-LAT Catalog of High-Energy Sources (3FHL). All the spectra were taken with the 4m Mayall telescope at Kitt Peak. With this follow-up program we are able to classify 27 out of 28 objects as BL Lacs, while the remaining one is a flat spectrum radio quasar. We determine a redshift (z) for three of these objects and a lower limit on z for other four sources: the farthest object for which we obtain a redshift has z>0.836. These results are part of a more extended campaign of optical spectroscopy follow-up of 3FHL blazars, aimed to obtain a complete sample of blazars at >10 GeV which will then be used to extend our knowledge on blazar emission mechanisms and on the extragalactic background light.

We present an X-ray stacking analysis of $\sim$75,000 star-forming galaxies between $0.1<z<5.0$ using the Chandra COSMOS Legacy survey to study the X-ray emission of low-luminosity active galactic nuclei (AGN) and its connection to host galaxy properties. The stacks at $z<0.9$ have luminosity limits as low as $10^{40}-10^{41}$ erg s$^{-1}$, a regime in which X-ray binaries (XRBs) can dominate the X-ray emission. Comparing the measured luminosities to established XRB scaling relations, we find that the redshift evolution of the luminosity per star formation rate (SFR) of XRBs depends sensitively on the assumed obscuration and may be weaker than previously found. The XRB scaling relation based on stacks from the Chandra Deep Field South overestimates the XRB contribution to the COSMOS high specific SFR (sSFR) stacks, possibly due to a bias affecting the CDF-S stacks because of their small galaxy samples. After subtracting the estimated XRB contribution from the stacks, we find that most stacks at $z>1.3$ exhibit a significant X-ray excess indicating nuclear emission. The AGN emission is strongly correlated with stellar mass but does not exhibit an additional correlation with SFR. The hardness ratios of the high-redshift stacks indicate that the AGN are substantially obscured ($N_{\mathrm{H}}\sim10^{23}$ cm$^{-2}$). These obscured AGN are not identified by IRAC color selection and have $L_X\sim10^{41}-10^{43}$ erg s$^{-1}$, consistent with accretion at an Eddington rate of $\sim10^{-3}$ onto $10^7-10^8 M_{\odot}$ black holes. Combining our results with other X-ray studies suggests that AGN obscuration depends on stellar mass and an additional variable, possibly the Eddington rate.

We present a sample of 40 AGN in dwarf galaxies at redshifts $z \lesssim$ 2.4. The galaxies are drawn from the \textit{Chandra} COSMOS-Legacy survey as having stellar masses $10^{7}\leq M_{*}\leq3 \times 10^{9}$ M$_{\odot}$. Most of the dwarf galaxies are star-forming. After removing the contribution from star formation to the X-ray emission, the AGN luminosities of the 40 dwarf galaxies are in the range $L_\mathrm{0.5-10 keV} \sim10^{39} - 10^{44}$ erg s$^{-1}$. With 12 sources at $z > 0.5$, our sample constitutes the highest-redshift discovery of AGN in dwarf galaxies. The record-holder is cid\_1192, at $z = 2.39$ and with $L_\mathrm{0.5-10 keV} \sim 10^{44}$ erg s$^{-1}$. One of the dwarf galaxies has $M_\mathrm{*} = 6.6 \times 10^{7}$ M$_{\odot}$ and is the least massive galaxy found so far to host an AGN. All the AGN are of type 2 and consistent with hosting intermediate-mass black holes (BHs) with masses $\sim 10^{4} - 10^{5}$ M$_{\odot}$ and typical Eddington ratios $> 1\%$. We also study the evolution, corrected for completeness, of AGN fraction with stellar mass, X-ray luminosity, and redshift in dwarf galaxies out to $z$ = 0.7. We find that the AGN fraction for $10^{9}< M_{*}\leq3 \times 10^{9}$ M$_{\odot}$ and $L_\mathrm{X} \sim 10^{41}-10^{42}$ erg s$^{-1}$ is $\sim$0.4\% for $z \leq$ 0.3 and that it decreases with X-ray luminosity and decreasing stellar mass. Unlike massive galaxies, the AGN fraction seems to decrease with redshift, suggesting that AGN in dwarf galaxies evolve differently than those in high-mass galaxies. Mindful of potential caveats, the results seem to favor a direct collapse formation mechanism for the seed BHs in the early Universe.

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 2-100 keV spectral analysis of 30 candidate Compton thick (CT-) active galactic nuclei (AGN) selected in the Swift-BAT 100-month survey. The average redshift of these objects is $\langle z\rangle\sim$0.03 and they all lie within $\sim$500 Mpc. We used the MyTorus (Murphy et al. 2009) model to perform X-ray spectral fitting both without and with the contribution of the NuSTAR data in the 3-50 keV energy range. When the NuSTAR data are added to the fit, 14 out of 30 of these objects (47% of the whole sample) have intrinsic absorption N$_{\rm H}<10^{24}$ cm$^{-2}$ at the $>$3$\sigma$ confidence level, i.e., they are re-classified from Compton thick to Compton thin. Consequently, we infer an overall observed fraction of CT-AGN with respect to the whole AGN population lower than the one reported in previous works, and as low as $\sim$4%. We find evidence that this over-estimation of N$_{\rm H}$ is likely due to the low quality of a subsample of spectra, either in the 2-10 keV band or in the Swift-BAT one.

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.

We present a study of the connection between black hole accretion, star formation, and galaxy morphology at z<=2.5. We focus on active galactic nuclei (AGNs) selected by their mid-IR power-law emission. By fitting optical to far-IR photometry with state-of-the-art spectral energy distribution (SED) techniques, we derive stellar masses, star formation rates, dust properties, and AGN contributions in galaxies over the whole COSMOS field. We find that obscured AGNs lie within or slightly above the star-forming sequence. We confirm our previous finding about compact host galaxies of obscured AGNs at z~1, and find that galaxies with 20-50% AGN contributions tend to have smaller sizes, by ~25-50%, compared to galaxies without AGNs. Furthermore, we find that a high merger fraction of up to 0.5 is appropriate for the most luminous (log (LIR/Lsun) ~ 12.5) AGN hosts and non-AGN galaxies, but not for the whole obscured AGN sample. Moreover, merger fraction depends on the total and star-forming infrared luminosity, rather than the decomposed AGN infrared luminosity. Our results suggest that major mergers are not the main driver of AGN activity, and therefore obscured AGNs might be triggered by internal mechanisms, such as secular processes, disk instabilities, and compaction in a particular evolutionary stage. We make the SED modeling results publicly available.

Multi-wavelength surveys covering large sky volumes are necessary to obtain an accurate census of rare objects such as high luminosity and/or high redshift active galactic nuclei (AGN). Stripe 82X is a 31.3 deg$^2$ X-ray survey with $Chandra$ and $XMM$-Newton observations overlapping the legacy Sloan Digital Sky Survey (SDSS) Stripe 82 field, which has a rich investment of multi-wavelength coverage from the ultraviolet to the radio. The wide-area nature of this survey presents new challenges for photometric redshifts for AGN compared to previous work on narrow-deep fields because it probes different populations of objects that need to be identified and represented in the library of templates. Here we present an updated X-ray plus multi-wavelength matched catalog, including $Spitzer$ counterparts, and estimated photometric redshifts for 5961 (96% of a total of 6181) X-ray sources, which have a normalized median absolute deviation, $\sigma_{\rm nmad}$ = 0.06 and an outlier fraction, $\eta$ = 13.7%. The populations found in this survey, and the template libraries used for photometric redshifts, provide important guiding principles for upcoming large-area surveys such as $eROSITA$ and 3$XMM$ (in X-ray) and the Large Synoptic Survey Telescope (LSST; optical).

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 present our statistical study of near infrared (NIR) variability of X-ray selected Active Galactic Nuclei (AGN) in the COSMOS field, using UltraVISTA data. This is the largest sample of AGN light curves in YJHKs bands, making possible to have a global description of the nature of AGN for a large range of redshifts, and for different levels of obscuration. To characterize the variability properties of the sources we computed the Structure Function. Our results show that there is an anti-correlation between the Structure Function $A$ parameter (variability amplitude) and the wavelength of emission, and a weak anti-correlation between $A$ and the bolometric luminosity. We find that Broad Line (BL) AGN have a considerably larger fraction of variable sources than Narrow Line (NL) AGN, and that they have different distributions of the $A$ parameter. We find evidence that suggests that most of the low luminosity variable NL sources correspond to BL AGN, where the host galaxy could be damping the variability signal. For high luminosity variable NL, we propose that they can be examples of "True type II" AGN or BL AGN with limited spectral coverage which results in missing the Broad Line emission. We also find that the fraction of variable sources classified as unobscured in the X-ray is smaller than the fraction of variable sources unobscured in the optical range. We present evidence that this is related to the differences in the origin of the obscuration in the optical and X-ray regimes.

Based on current models of the Cosmic X-ray Background (CXB), heavily obscured Active Galactic Nuclei (AGN) are expected to make up ~10% of the peak emission of the CXB and ~20% of the total population of AGN, yet few of these sources have been recorded and characterized in current surveys. Here we present the Chandra follow-up observation of 14 AGN detected by Swift-BAT. For five sources in the sample, NuSTAR observations in the 3-80 keV band are also available. The X-ray spectral fitting over the 0.3-150 keV energy range allows us to determine the main X-ray spectral parameters, such as the photon index and the intrinsic absorption, of these objects, and to make hypotheses on the physical structures responsible for the observed spectra. We find that 13 of the 14 objects are absorbed AGN, and one is a candidate Compton thick AGN, having intrinsic absorption NH>1E24 cm^{-2}. Finally, we verified that the use of NuSTAR observations is strategic to strongly constrain the properties of obscured AGN, since the best-fit values we obtained for parameters such as the power-law photon index and the intrinsic absorption NH changed sometimes significantly fitting the spectra with and without the use of NuSTAR data.

We carried out targeted ALMA observations of 129 fields in the COSMOS region at 1.25 mm, detecting 152 galaxies at S/N$\geq$5 with an average continuum RMS of 150 $\mu$Jy. These fields represent a S/N-limited sample of AzTEC / ASTE sources with 1.1 mm S/N$\geq$4 over an area of 0.72 square degrees. Given ALMA's fine resolution and the exceptional spectroscopic and multiwavelength photometric data available in COSMOS, this survey allows us unprecedented power in identifying submillimeter galaxy counterparts and determining their redshifts through spectroscopic or photometric means. In addition to 30 sources with prior spectroscopic redshifts, we identified redshifts for 113 galaxies through photometric methods and an additional nine sources with lower limits, which allowed a statistically robust determination of the redshift distribution. We have resolved 33 AzTEC sources into multi-component systems and our redshifts suggest that nine are likely to be physically associated. Our overall redshift distribution peaks at $z\sim$2.0 with a high redshift tail skewing the median redshift to $\tilde{z}$=2.48$\pm$0.05. We find that brighter millimeter sources are preferentially found at higher redshifts. Our faintest sources, with S$_{1.25 \rm mm}$<1.25 mJy, have a median redshift of $\tilde{z}$=2.18$\pm$0.09, while the brightest sources, S$_{1.25 \rm mm}$>1.8 mJy, have a median redshift of $\tilde{z}$=3.08$\pm$0.17. After accounting for spectral energy distribution shape and selection effects these results are consistent with several previous submillimeter galaxy surveys, and moreover, support the conclusion that the submillimeter galaxy redshift distribution is sensitive to survey depth.

We present results of a ground-based near-infrared campaign with Palomar TripleSpec, Keck NIRSPEC, and Gemini GNIRS to target two samples of reddened active galactic nucleus (AGN) candidates from the 31 deg$^2$ Stripe 82 X-ray survey. One sample, which is $\sim$89\% complete to $K<16$ (Vega), consists of eight confirmed AGNs, four of which were identified with our follow-up program, and is selected to have red $R-K$ colors ($>4$, Vega). The fainter sample ($K>17$, Vega) represents a pilot program to follow-up four sources from a parent sample of 34 that are not detected in the single-epoch SDSS catalog and have {\it WISE} quasar colors. All twelve sources are broad-line AGNs (at least one permitted emission line has a FWHM exceeding 1300 km s$^{-1}$) and span a redshift range $0.59 < z < 2.5$. Half the ($R-K$)-selected AGNs have features in their spectra suggestive of outflows. When comparing these sources to a matched sample of blue Type 1 AGNs, we find the reddened AGNs are more distant ($z > 0.5$) and a greater percentage have high X-ray luminosities ($L_{\rm X,full} > 10^{44}$ erg s$^{-1}$). Such outflows and high luminosities may be consistent with the paradigm that reddened broad-line AGNs represent a transitory phase in AGN evolution as described by the major merger model for black hole growth. Results from our pilot program demonstrate proof-of-concept that our selection technique is successful in discovering reddened quasars at $z > 1$ missed by optical surveys.

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.

To fully understand cosmic black hole growth we need to constrain the population of heavily obscured active galactic nuclei (AGN) at the peak of cosmic black hole growth ($z\sim$1-3). Sources with obscuring column densities higher than $\mathrm{10^{24}}$ atoms $\mathrm{cm^{-2}}$, called Compton-thick (CT) AGN, can be identified by excess X-ray emission at $\sim$20-30 keV, called the "Compton hump". We apply the recently developed Spectral Curvature (SC) method to high-redshift AGN (2<z<5) detected with Chandra. This method parametrizes the characteristic "Compton hump" feature cosmologically redshifted into the X-ray band at observed energies <10 keV. We find good agreement in CT AGN found using the SC method and bright sources fit using their full spectrum with X-ray spectroscopy. In the Chandra deep field south, we measure a CT fraction of $\mathrm{17^{+19}_{-11}\%}$ (3/17) for sources with observed luminosity $\mathrm{>5\times 10^{43}}$ erg $\mathrm{s^{-1}}$. In the Cosmological evolution survey (COSMOS), we find an observed CT fraction of $\mathrm{15^{+4}_{-3}\%}$ (40/272) or $\mathrm{32\pm11 \%}$ when corrected for the survey sensitivity. When comparing to low redshift AGN with similar X-ray luminosities, our results imply the CT AGN fraction is consistent with having no redshift evolution. Finally, we provide SC equations that can be used to find high-redshift CT AGN (z>1) for current (XMM-Newton) and future (eROSITA and ATHENA) X-ray missions.

Based on a sample of over 1,800 radio AGN at redshifts out to z~5, which have typical stellar masses within ~3x(10^{10}-10^{11}) Msol, and 3 GHz radio data in the COSMOS field, we derived the 1.4 GHz radio luminosity functions for radio AGN (L_1.4GHz ~ 10^{22}-10^{27} W/Hz) out to z~5. We constrained the evolution of this population via continuous models of pure density and pure luminosity evolutions, and we found best-fit parametrizations of Phi*~(1+z)^{(2.00+/-0.18)-(0.60+/-0.14)z}, and L*~(1+z)^{(2.88+/-0.82)-(0.84+/-0.34)z}, respectively, with a turnover in number and luminosity densities of the population at z~1.5. We converted 1.4 GHz luminosity to kinetic luminosity taking uncertainties of the scaling relation used into account. We thereby derived the cosmic evolution of the kinetic luminosity density provided by the AGN and compared this luminosity density to the radio-mode AGN feedback assumed in the Semi-Analytic Galaxy Evolution (SAGE) model, i.e., to the redshift evolution of the central supermassive black hole accretion luminosity taken in the model as the source of heating that offsets the energy losses of the cooling, hot halo gas, and thereby limits further stellar mass growth of massive galaxies. We find that the kinetic luminosity exerted by our radio AGN may be high enough to balance the radiative cooling of the hot gas at each cosmic epoch since z~5. However, although our findings support the idea of radio-mode AGN feedback as a cosmologically relevant process in massive galaxy formation, many simplifications in both the observational and semi-analytic approaches still remain and need to be resolved before robust conclusions can be reached.

We investigate the star formation properties of a large sample of ~2300 X-ray-selected Type 2 Active Galactic Nuclei (AGNs) host galaxies out to z~3 in the Chandra COSMOS Legacy Survey in order to understand the connection between the star formation and nuclear activity. Making use of the existing multi-wavelength photometric data available in the COSMOS field, we perform a multi-component modeling from far-infrared to near-ultraviolet using a nuclear dust torus model, a stellar population model and a starburst model of the spectral energy distributions (SEDs). Through detailed analysis of SEDs, we derive the stellar masses and the star formation rates (SFRs) of Type 2 AGN host galaxies. The stellar mass of our sample is in the range 9 < log M_{stellar}/M_{\odot} < 12 with uncertainties of ~0.19 dex. We find that Type 2 AGN host galaxies have, on average, similar SFRs compared to the normal star-forming galaxies with similar M_{stellar} and redshift ranges, suggesting no significant evidence for enhancement or quenching of star formation. This could be interpreted in a scenario, where the relative massive galaxies have already experienced substantial growth at higher redshift (z>3), and grow slowly through secular fueling processes hosting moderate-luminosity AGNs.

We examine the behaviour of the infrared-radio correlation (IRRC) over the range $0<z<6$ using new, highly sensitive 3GHz observations with the Karl G. Jansky Very Large Array (VLA) and infrared data from the Herschel Space Observatory in the 2deg$^{2}$ COSMOS field. We distinguish between objects where emission is believed to arise solely from star-formation, and those where an active galactic nucleus (AGN) is thought to be present. We account for non-detections in the radio or in the infrared using a doubly-censored survival analysis. We find that the IRRC of star-forming galaxies, quantified by the infrared-to-1.4GHz radio luminosity ratio ($q_{\rm TIR}$), decreases with increasing redshift: $q_{\rm TIR}(z)=(2.88\pm0.03)(1+z)^{-0.19\pm0.01}$. Moderate-to-high radiative luminosity AGN do not follow the same $q_{\rm TIR}$$(z)$ trend, having a lower normalisation and steeper decrease with redshift. We cannot rule out the possibility that unidentified AGN contributions only to the radio regime may be steepening the observed $q_{\rm TIR}(z)$ trend of the star-forming population. An increasing fractional contribution to the observed 3GHz flux by free-free emission of star-forming galaxies may also affect the derived evolution. However, we find that the standard (M82-based) assumption of the typical radio spectral energy distribution (SED) for star-forming galaxies is inconsistent with our results. This suggests a more complex shape of the typical radio SED for star-forming galaxies, and that imperfect $K$ corrections in the radio may govern the derived redshift trend of $q_{\rm TIR}$. Lastly, we present a redshift-dependent relation between rest-frame 1.4GHz radio luminosity and star formation rate taking the derived redshift trend into account.

(abridged) We study the composition of the faint radio population selected from the VLA-COSMOS 3GHz Large Project. The survey covers a 2.6sq.deg. area with a mean rms of ~2.3uJy/b, cataloging 10830 sources (>5sigma). Combining these radio data with optical, near-infrared (UltraVISTA), mid-infrared (Spitzer/IRAC) data, and X-ray data (Chandra), we find counterparts to radio sources for ~93% of the radio sample (in the areas of the COSMOS field not affected by saturated or bright sources in the optical to NIR bands), reaching out to z<6. We further classify the sources as star forming galaxies or AGN based on various criteria, such as X-ray luminosity, observed MIR color, UV-FIR spectral-energy distribution, rest-frame NUV-optical color corrected for dust extinction, and radio-excess relative to that expected from the hosts' star-formation rate. We separate the AGN into sub-samples dominated by low-to-moderate and moderate-to-high radiative luminosity AGN, candidates for high-redshift analogues to local low- and high-excitation emission line AGN, respectively. We study the fractional contributions of these sub-populations down to radio flux levels of ~11uJy at 3GHz (or ~20uJy at 1.4GHz assuming a spectral index of -0.7). We find that the dominant fraction at 1.4GHz flux densities above ~200uJy is constituted of low-to-moderate radiative luminosity AGN. Below densities of ~100uJy the fraction of star-forming galaxies increases to ~60%, followed by the moderate-to-high radiative luminosity AGN (~20%), and low-to-moderate radiative luminosity AGN (~20%). Based on this observational evidence, we extrapolate the fractions down to sensitivities of the SKA. Our estimates suggest that at the faint flux limits to be reached by the SKA1 surveys, a selection based only on radio flux limits can provide a simple tool to efficiently identify samples highly (>75%) dominated by star-forming galaxies.

We explore the multiwavelength properties of AGN host galaxies for different classes of radio-selected AGN out to z$\lesssim$6 via a multiwavelength analysis of about 7700 radio sources in the COSMOS field. The sources were selected with the Very Large Array (VLA) at 3 GHz (10 cm) within the VLA-COSMOS 3 GHz Large Project, and cross-matched with multiwavelength ancillary data. This is the largest sample of high-redshift (z$\lesssim$6) radio sources with exquisite photometric coverage and redshift measurements available. We constructed a sample of moderate-to-high radiative luminosity AGN (HLAGN) via spectral energy distribution (SED) decomposition combined with standard X-ray and mid-infrared diagnostics. Within the remainder of the sample we further identified low-to-moderate radiative luminosity AGN (MLAGN) via excess in radio emission relative to the star formation rates in their host galaxies. We show that AGN power in HLAGN occurs predominantly in radiative form, while MLAGN display a substantial mechanical AGN luminosity component. We found significant differences in the host properties of the two AGN classes, as a function of redshift. At z$<$1.5, MLAGN appear to reside in significantly more massive and less star-forming galaxies compared to HLAGN. At z$>$1.5, we observed a reversal in the behaviour of the stellar mass distributions with the HLAGN populating the higher stellar mass tail. We interpret this finding as a possible hint of the downsizing of galaxies hosting HLAGN, with the most massive galaxies triggering AGN activity earlier than less massive galaxies, and then fading to MLAGN at lower redshifts. Our conclusion is that HLAGN and MLAGN samples trace two distinct galaxy and AGN populations in a wide range of redshifts, possibly resembling the radio AGN types often referred to as radiative- and jet-mode (or high- and low-excitation), respectively.

We make use of the deep Karl G. Jansky Very Large Array (VLA) COSMOS radio observations at 3 GHz to infer radio luminosity functions of star-forming galaxies up to redshifts of z~5 based on approximately 6000 detections with reliable optical counterparts. This is currently the largest radio-selected sample available out to z~5 across an area of 2 square degrees with a sensitivity of rms=2.3 ujy/beam. By fixing the faint and bright end shape of the radio luminosity function to the local values, we find a strong redshift trend that can be fitted with a pure luminosity evolution L~(1+z)^{(3.16 +- 0.2)-(0.32 +- 0.07) z}. We estimate star formation rates (SFRs) from our radio luminosities using an infrared (IR)-radio correlation that is redshift dependent. By integrating the parametric fits of the evolved luminosity function we calculate the cosmic SFR density (SFRD) history since z~5. Our data suggest that the SFRD history peaks between 2<z<3 and that the ultraluminous infrared galaxies (ULIRGs; 100 Msol/yr<SFR<1000 Msol/yr) contribute up to ~25% to the total SFRD in the same redshift range. Hyperluminous infrared galaxies (HyLIRGs; SFR>1000 Msol/yr) contribute an additional <2% in the entire observed redshift range. We find evidence of a potential underestimation of SFRD based on ultraviolet (UV) rest-frame observations of Lyman break galaxies (LBGs) at high redshifts (z>4) on the order of 15-20%, owing to appreciable star formation in highly dust-obscured galaxies, which might remain undetected in such UV observations.

We determine the radio size distribution of a large sample of 152 SMGs in COSMOS that were detected with ALMA at 1.3 mm. For this purpose, we used the observations taken by the VLA-COSMOS 3 GHz Large Project. One hundred and fifteen of the 152 target SMGs were found to have a 3 GHz counterpart. The median value of the major axis FWHM at 3 GHz is derived to be $4.6\pm0.4$ kpc. The radio sizes show no evolutionary trend with redshift, or difference between different galaxy morphologies. We also derived the spectral indices between 1.4 and 3 GHz, and 3 GHz brightness temperatures for the sources, and the median values were found to be $\alpha=-0.67$ and $T_{\rm B}=12.6\pm2$ K. Three of the target SMGs, which are also detected with the VLBA, show clearly higher brightness temperatures than the typical values. Although the observed radio emission appears to be predominantly powered by star formation and supernova activity, our results provide a strong indication of the presence of an AGN in the VLBA and X-ray-detected SMG AzTEC/C61. The median radio-emitting size we have derived is 1.5-3 times larger than the typical FIR dust-emitting sizes of SMGs, but similar to that of the SMGs' molecular gas component traced through mid-$J$ line emission of CO. The physical conditions of SMGs probably render the diffusion of cosmic-ray electrons inefficient, and hence an unlikely process to lead to the observed extended radio sizes. Instead, our results point towards a scenario where SMGs are driven by galaxy interactions and mergers. Besides triggering vigorous starbursts, galaxy collisions can also pull out the magnetised fluids from the interacting disks, and give rise to a taffy-like synchrotron-emitting bridge. This provides an explanation for the spatially extended radio emission of SMGs, and can also cause a deviation from the well-known IR-radio correlation.

The coeval AGN and galaxy evolution and the observed local relations between SMBHs and galaxy properties suggest some connection or feedback between SMBH growth and galaxy build-up. We looked for correlations between properties of X-ray detected AGN and their FIR detected host galaxies, to find quantitative evidences for this connection, highly debated in the latest years. We exploit the rich multi-wavelength data set available in the COSMOS field for a large sample (692 sources) of AGN and their hosts, in the redshift range $0.1<z<4$. We use X-ray data to select AGN and determine their properties (intrinsic luminosity and nuclear obscuration), and broad-band SED fitting to derive host galaxy properties (stellar mass $M_*$ and star formation rate SFR). We find that the AGN 2-10 keV luminosity ($L_{\rm X}$) and the host $8-1000~\mu m$ star formation luminosity ($L_{\rm IR}^{\rm SF}$) are significantly correlated. However, the average host $L_{\rm IR}^{\rm SF}$ has a flat distribution in bins of AGN $L_{\rm X}$, while the average AGN $L_{\rm X}$ increases in bins of host $L_{\rm IR}^{\rm SF}$, with logarithmic slope of $\sim0.7$, in the redshifts range $0.4<z<1.2$. We also discuss the comparison between the distribution of these two quantities and the predictions from hydro-dynamical simulations. Finally we find that the average column density ($N_H$) shows a positive correlation with the host $M_*$, at all redshifts, but not with the SFR (or $L_{\rm IR}^{\rm SF}$). This translates into a negative correlation with specific SFR. Our results are in agreement with the idea that BH accretion and SF rates are correlated, but occur with different variability time scales. The presence of a positive correlation between $N_H$ and host $M_*$ suggests that the X-ray $N_H$ is not entirely due to the circum-nuclear obscuring torus, but may also include a contribution from the host galaxy.

Using {\em Chandra} observations in the 2.15 deg$^{2}$ COSMOS legacy field, we present one of the most accurate measurements of the Cosmic X-ray Background (CXB) spectrum to date in the [0.3-7] keV energy band. The CXB has three distinct components: contributions from two Galactic collisional thermal plasmas at kT$\sim$0.27 and 0.07 keV and an extragalactic power-law with photon spectral index $\Gamma$=1.45$\pm{0.02}$. The 1 keV normalization of the extragalactic component is 10.91$\pm{0.16}$ keV cm$^{-2}$ s$^{-1}$ sr$^{-1}$ keV$^{-1}$. Removing all X-ray detected sources, the remaining unresolved CXB is best-fit by a power-law with normalization 4.18$\pm{0.26}$ keV cm$^{-2}$ s$^{-1}$ sr$^{-1}$ keV$^{-1}$ and photon spectral index $\Gamma$=1.57$\pm{0.10}$. Removing faint galaxies down to i$_{AB}\sim$27-28 leaves a hard spectrum with $\Gamma\sim$1.25 and a 1 keV normalization of $\sim$1.37 keV cm$^{-2}$ s$^{-1}$ sr$^{-1}$ keV$^{-1}$. This means that $\sim$91\% of the observed CXB is resolved into detected X-ray sources and undetected galaxies. Unresolved sources that contribute $\sim 8-9\%$ of the total CXB show a marginal evidence of being harder and possibly more obscured than resolved sources. Another $\sim$1\% of the CXB can be attributed to still undetected star forming galaxies and absorbed AGN. According to these limits, we investigate a scenario where early black holes totally account for non source CXB fraction and constrain some of their properties. In order to not exceed the remaining CXB and the $z\sim$6 accreted mass density, such a population of black holes must grow in Compton-thick envelopes with N$_{H}>$1.6$\times$10$^{25}$ cm$^{-2}$ and form in extremely low metallicity environments $(Z_\odot)\sim10^{-3}$.

We present the combined Chandra and Swift-BAT spectral analysis of seven Seyfert 2 galaxies selected from the Swift-BAT 100-month catalog. We selected nearby (z<=0.03) sources lacking of a ROSAT counterpart and never previously observed with Chandra in the 0.3-10 keV energy range, and targeted these objects with 10 ks Chandra ACIS-S observations. The X-ray spectral fitting over the 0.3-150 keV energy range allows us to determine that all the objects are significantly obscured, having NH>=1E23 cm^(-2) at a >99% confidence level. Moreover, one to three sources are candidate Compton thick Active Galactic Nuclei (CT-AGN), i.e., have NH>=1E24 cm^(-2). We also test the recent "spectral curvature" method developed by Koss et al. (2016) to find candidate CT-AGN, finding a good agreement between our results and their predictions. Since the selection criteria we adopted have been effective in detecting highly obscured AGN, further observations of these and other Seyfert 2 galaxies selected from the Swift-BAT 100-month catalog will allow us to create a statistically significant sample of highly obscured AGN, therefore better understanding the physics of the obscuration processes.

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.

We present a structural study of 182 obscured Active Galactic Nuclei (AGNs) at z<=1.5, selected in the COSMOS field from their extreme infrared to X-ray luminosity ratio and their negligible emission at optical wavelengths. We fit optical to far-infrared spectral energy distributions and analyze deep HST imaging to derive the physical and morphological properties of their host galaxies. We find that such galaxies are more compact than normal star-forming sources at similar redshift and stellar mass, and we show that it is not an observational bias related to the emission of the AGN. Based on the distribution of their UVJ colors, we also argue that this increased compactness is not due to the additional contribution of a passive bulge. We thus postulate that a vast majority of obscured AGNs reside in galaxies undergoing dynamical compaction, similar to processes recently invoked to explain the formation of compact star-forming sources at high redshift.

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 present the measurement of the projected and redshift space 2-point correlation function (2pcf) of the new catalog of Chandra COSMOS-Legacy AGN at 2.9$\leq$z$\leq$5.5 ($\langle L_{bol} \rangle \sim$10$^{46}$ erg/s) using the generalized clustering estimator based on phot-z probability distribution functions (Pdfs) in addition to any available spec-z. We model the projected 2pcf estimated using $\pi_{max}$ = 200 h$^{-1}$ Mpc with the 2-halo term and we derive a bias at z$\sim$3.4 equal to b = 6.6$^{+0.60}_{-0.55}$, which corresponds to a typical mass of the hosting halos of log M$_h$ = 12.83$^{+0.12}_{-0.11}$ h$^{-1}$ M$_{\odot}$. A similar bias is derived using the redshift-space 2pcf, modelled including the typical phot-z error $\sigma_z$ = 0.052 of our sample at z$\geq$2.9. Once we integrate the projected 2pcf up to $\pi_{max}$ = 200 h$^{-1}$ Mpc, the bias of XMM and \textit{Chandra} COSMOS at z=2.8 used in Allevato et al. (2014) is consistent with our results at higher redshift. The results suggest only a slight increase of the bias factor of COSMOS AGN at z$\gtrsim$3 with the typical hosting halo mass of moderate luminosity AGN almost constant with redshift and equal to logM$_h$ = 12.92$^{+0.13}_{-0.18}$ at z=2.8 and log M$_h$ = 12.83$^{+0.12}_{-0.11}$ at z$\sim$3.4, respectively. The observed redshift evolution of the bias of COSMOS AGN implies that moderate luminosity AGN still inhabit group-sized halos at z$\gtrsim$3, but slightly less massive than observed in different independent studies using X-ray AGN at z$\leq2$.

We present the X-ray spectral analysis of the 1855 extragalactic sources in the Chandra COSMOS-Legacy survey catalog having more than 30 net counts in the 0.5-7 keV band. 38% of the sources are optically classified Type 1 active galactic nuclei (AGN), 60% are Type 2 AGN and 2% are passive, low-redshift galaxies. We study the distribution of AGN photon index and of the intrinsic absorption N(H,z) based on the sources optical classification: Type 1 have a slightly steeper mean photon index than Type 2 AGN, which on the other hand have average intrinsic absorption ~3 times higher than Type 1 AGN. We find that ~15% of Type 1 AGN have N(H,z)>1E22 cm^(-2), i.e., are obscured according to the X-ray spectral fitting; the vast majority of these sources have L(2-10keV)>$1E44 erg/s. The existence of these objects suggests that optical and X-ray obscuration can be caused by different phenomena, the X-ray obscuration being for example caused by dust-free material surrounding the inner part of the nuclei. ~18% of Type 2 AGN have N(H,z)<1E22 cm^(-2), and most of these sources have low X-ray luminosities (L(2-10keV)<$1E43 erg/s). We expect a part of these sources to be low-accretion, unobscured AGN lacking of broad emission lines. Finally, we also find a direct proportional trend between N(H,z) and host galaxy mass and star formation rate, although part of this trend is due to a redshift selection effect.

We present the largest high-redshift (3<z<6.85) sample of X-ray-selected active galactic nuclei (AGN) on a contiguous field, using sources detected in the Chandra COSMOS Legacy survey. The sample contains 174 sources, 87 with spectroscopic redshift, the other 87 with photometric redshift (z_phot). In this work we treat z_phot as a probability weighted sum of contributions, adding to our sample the contribution of sources with z_phot<3 but z_phot probability distribution >0 at z>3. We compute the number counts in the observed 0.5-2 keV band, finding a decline in the number of sources at z>3 and constraining phenomenological models of X-ray background. We compute the AGN space density at z>3 in two different luminosity bins. At higher luminosities (logL(2-10 keV) > 44.1 erg/s) the space density declines exponentially, dropping by a factor ~20 from z~3 to z~6. The observed decline is ~80% steeper at lower luminosities (43.55 erg/s < logL(2-10 keV) < 44.1 erg/s), from z~3 to z~4.5. We study the space density evolution dividing our sample in optically classified Type 1 and Type 2 AGN. At logL(2-10 keV) > 44.1 erg/s, unobscured and obscured objects may have different evolution with redshift, the obscured component being three times higher at z~5. Finally, we compare our space density with predictions of quasar activation merger models, whose calibration is based on optically luminous AGN. These models significantly overpredict the number of expected AGN at logL(2-10 keV) > 44.1 erg/s with respect to our data.

We present new Keck/MOSFIRE K-band spectroscopy for a sample of 14 faint, X-ray-selected active galactic nuclei (AGNs) in the COSMOS field. The data cover the spectral region surrounding the broad Balmer emission lines, which enables the estimation of black hole masses (M_BH) and accretion rates (in terms of L/L_Edd). We focus on 10 AGN at z~3.3, where we observe the Hbeta spectral region, while for the other four z~2.4 sources we use the Halpha broad emission line. Compared with previous detailed studies of unobscured AGNs at these high redshifts, our sources are fainter by an order of magnitude, corresponding to number densities of order ~10^-6--10^-5 Mpc^-3. The lower AGN luminosities also allow for a robust identification of the host galaxy emission, necessary to obtain reliable intrinsic AGN luminosities, BH masses and accretion rates. We find the AGNs in our sample to be powered by supermassive black holes (SMBHs) with a typical mass of M_BH~5*10^8 M_Sol -- significantly lower than the higher-luminosity, rarer quasars reported in previous studies. The accretion rates are in the range L/L_Edd~0.1-0.5, with an evident lack of sources with lower L/L_Edd (and higher M_BH), as found in several studies of faint AGNs at intermediate redshifts. Based on the early growth expected for the SMBHs in our sample, we argue that a significant population of faint z~6 AGNs, with M_BH~10^6 M_Sol, should be detectable in the deepest X-ray surveys available, but this is not observed. We discuss several possible explanations for the apparent absence of such a population, concluding that the most probable scenario involves an evolution in source obscuration and/or radiative efficiencies.

We present the COSMOS2015 catalog which contains precise photometric redshifts and stellar masses for more than half a million objects over the 2deg$^{2}$ COSMOS field. Including new $YJHK_{\rm s}$ images from the UltraVISTA-DR2 survey, $Y$-band from Subaru/Hyper-Suprime-Cam and infrared data from the Spitzer Large Area Survey with the Hyper-Suprime-Cam Spitzer legacy program, this near-infrared-selected catalog is highly optimized for the study of galaxy evolution and environments in the early Universe. To maximise catalog completeness for bluer objects and at higher redshifts, objects have been detected on a $\chi^{2}$ sum of the $YJHK_{\rm s}$ and $z^{++}$ images. The catalog contains $\sim 6\times 10^5$ objects in the 1.5 deg$^{2}$ UltraVISTA-DR2 region, and $\sim 1.5\times 10^5$ objects are detected in the "ultra-deep stripes" (0.62 deg$^{2}$) at $K_{\rm s}\leq 24.7$ (3$\sigma$, 3", AB magnitude). Through a comparison with the zCOSMOS-bright spectroscopic redshifts, we measure a photometric redshift precision of $\sigma_{\Delta z/(1+z_s)}$ = 0.007 and a catastrophic failure fraction of $\eta=0.5$%. At $3<z<6$, using the unique database of spectroscopic redshifts in COSMOS, we find $\sigma_{\Delta z/(1+z_s)}$ = 0.021 and $\eta=13.2\% $. The deepest regions reach a 90\% completeness limit of 10$^{10}M_\odot$ to $z=4$. Detailed comparisons of the color distributions, number counts, and clustering show excellent agreement with the literature in the same mass ranges. COSMOS2015 represents a unique, publicly available, valuable resource with which to investigate the evolution of galaxies within their environment back to the earliest stages of the history of the Universe. The COSMOS2015 catalog is distributed via anonymous ftp (ftp://ftp.iap.fr/pub/from_users/hjmcc/COSMOS2015/) and through the usual astronomical archive systems (CDS, ESO Phase 3, IRSA).

We study the environment of 23 submillimetre galaxies (SMGs) drawn from the JCMT/AzTEC 1.1mm S/N-limited sample in the COSMOS field, as well as 4 COSMOS SMGs at z_spec>4.5, and 1 at z_spec=2.49, yielding a sample of 28 SMGs. We search for overdensities using the COSMOS photometric redshifts based on over 30 UV-NIR photometric bands, reaching an accuracy of sigma(Delta z/(1+z))=0.0067 (0.0155) at z<3.5 (>3.5). To identify overdensities we apply the Voronoi tessellation analysis, and estimate the overdensity estimator delta_g as a function of distance from the SMG and/or overdensity center. We test and validate our approach via simulations, X-ray detected groups, and spectroscopic verifications using VUDS and zCOSMOS catalogues showing that even with photometric redshifts in COSMOS we can efficiently retrieve overdensities out to z~5. Our results yield that 11/23 (48%) JCMT/AzTEC 1.1mm SMGs occupy overdense environments. Considering the entire JCMT/AzTEC 1.1mm S/N>4 sample, and accounting for the expected fraction of spurious detections, yields that 35-61% of the SMGs in the S/N-limited sample occupy overdense environments. We perform an X-ray stacking analysis in the 0.5-2keV band using a 32" aperture and our SMG positions, and find statistically significant detections. For our z<2 [z>2] subsample we find an average flux of (4.0+/-0.8)x10^{-16} [(1.3+/-0.5)x10^{-16}] erg/s/cm^2 and a corresponding total mass of M_200 = 2.8x10^{13} [2x10^{13}] MSol. Our results suggest a higher occurrence of SMGs occupying overdense environments at z>3, than at z<3. This may be understood if highly star forming galaxies can only be formed in the highest peaks of the density field tracing the most massive dark matter haloes at early cosmic epochs, while at later times cosmic structure may have matured sufficiently that more modest overdensities correspond to sufficiently massive haloes to form SMGs.

We study the composition of the faint radio population selected from the VLA-COSMOS 3GHz Large Project, a radio continuum survey performed at 10 cm wavelength. The survey covers the full 2 square degree COSMOS field with mean $rms\sim2.3$ $\mu$Jy/beam, cataloging 10,899 source components above $5\times rms$. By combining these radio data with UltraVISTA, optical, near-infrared, and Spitzer/IRAC mid-infrared data, as well as X-ray data from the Chandra Legacy, and Chandra COSMOS surveys, we gain insight into the emission mechanisms within our radio sources out to redshifts of $z\sim5$. From these emission characteristics we classify our souces as star forming galaxies or AGN. Using their multi-wavelength properties we further separate the AGN into sub-samples dominated by radiatively efficient and inefficient AGN, often referred to as high- and low-excitation emission line AGN. We compare our method with other results based on fitting of the sources' spectral energy distributions using both galaxy and AGN spectral models, and those based on the infrared-radio correlation. We study the fractional contributions of these sub-populations down to radio flux levels of $\sim$10 $\mu$Jy. We find that at 3 GHz flux densities above $\sim$400 $\mu$Jy quiescent, red galaxies, consistent with the low-excitation radio AGN class constitute the dominant fraction. Below densities of $\sim$200 $\mu$Jy star-forming galaxies begin to constitute the largest fraction, followed by the low-excitation, and X-ray- and IR-identified high-excitation radio AGN.

We study a sample of $\sim$50,000 dwarf starburst and late-type galaxies drawn from the COSMOS survey with the aim of investigating the presence of nuclear accreting black holes (BHs) as those seed BHs from which supermassive BHs could grow in the early Universe. We divide the sample into five complete redshift bins up to $z=1.5$ and perform an X-ray stacking analysis using the \textit{Chandra} COSMOS-Legacy survey data. After removing the contribution from X-ray binaries and hot gas to the stacked X-ray emission, we still find an X-ray excess in the five redshift bins that can be explained by nuclear accreting BHs. This X-ray excess is more significant for $z<0.5$. At higher redshifts, these active galactic nuclei could suffer mild obscuration, as indicated by the analysis of their hardness ratios. The average nuclear X-ray luminosities in the soft band are in the range 10$^{39}-10^{40}$ erg s$^{-1}$. Assuming that the sources accrete at $\geq$ 1\% the Eddington rate, their BH masses would be $\leq$ 10$^{5}$ M$_{\odot}$, thus in the intermediate-mass BH regime, but their mass would be smaller than the one predicted by the BH-stellar mass relation. If instead the sources follow the correlation between BH mass and stellar mass, they would have sub-Eddington accreting rates of $\sim$ 10$^{-3}$ and BH masses 1-9 $\times$ 10$^{5}$ M$_{\odot}$. We thus conclude that a population of intermediate-mass BHs exists in dwarf starburst galaxies, at least up to $z$=1.5, though their detection beyond the local Universe is challenging due to their low luminosity and mild obscuration unless deep surveys are employed.

The COSMOS-Legacy survey is a 4.6 Ms Chandra program that has imaged 2.2 deg$^2$ of the COSMOS field with an effective exposure of $\simeq$160 ks over the central 1.5 deg$^2$ and of $\simeq$80 ks in the remaining area. The survey is the combination of 56 new observations, obtained as an X-ray Visionary Project, with the previous C-COSMOS survey. We describe the reduction and analysis of the new observations and the properties of 2273 point sources detected above a spurious probability of 2$\times 10^{-5}$. We also present the updated properties of the C-COSMOS sources detected in the new data. The whole survey includes 4016 point sources (3814, 2920 and 2440 in the full, soft and hard band). The limiting depths are 2.2 $\times$ 10$^{-16}$, 1.5 $\times$ 10$^{-15}$ and 8.9$\times$ 10$^{-16}$ ${\rm erg~cm}^{-2}~{\rm s}^{-1}$ in the 0.5-2, 2-10 and 0.5-10 keV bands, respectively. The observed fraction of obscured AGN with column density $> 10^{22}$ cm$^{-2}$ from the hardness ratio (HR) is $\sim$50$^{+17}_{-16}$%. Given the large sample, we compute source number counts in the hard and soft bands, significantly reducing the uncertainties of 5-10%. For the first time, we compute number counts for obscured (HR$>$-0.2) and unobscured (HR$<$-0.2) sources and find significant differences between the two populations in the soft band. Due to the un-precedent large exposure, COSMOS-Legacy area is 3 times larger than surveys at similar depth and its depth is 3 times fainter than surveys covering similar area. The area-flux region occupied by COSMOS-Legacy is likely to remain unsurpassed for years to come.

We present the catalog of optical and infrared counterparts of the Chandra COSMOS-Legacy Survey, a 4.6 Ms Chandra program on the 2.2 square degrees of the COSMOS field, combination of 56 new overlapping observations obtained in Cycle 14 with the previous C-COSMOS survey. In this Paper we report the i, K, and 3.6 micron identifications of the 2273 X-ray point sources detected in the new Cycle 14 observations. We use the likelihood ratio technique to derive the association of optical/infrared (IR) counterparts for 97% of the X-ray sources. We also update the information for the 1743 sources detected in C-COSMOS, using new K and 3.6 micron information not available when the C-COSMOS analysis was performed. The final catalog contains 4016 X-ray sources, 97% of which have an optical/IR counterpart and a photometric redshift, while 54% of the sources have a spectroscopic redshift. The full catalog, including spectroscopic and photometric redshifts and optical and X-ray properties described here in detail, is available online. We study several X-ray to optical (X/O) properties: with our large statistics we put better constraints on the X/O flux ratio locus, finding a shift towards faint optical magnitudes in both soft and hard X-ray band. We confirm the existence of a correlation between X/O and the the 2-10 keV luminosity for Type 2 sources. We extend to low luminosities the analysis of the correlation between the fraction of obscured AGN and the hard band luminosity, finding a different behavior between the optically and X-ray classified obscured fraction.

We release the next installment of the Stripe 82 X-ray survey point-source catalog, which currently covers 31.3 deg$^2$ of the Sloan Digital Sky Survey (SDSS) Stripe 82 Legacy field. In total, 6181 unique X-ray sources are significantly detected with {\it XMM-Newton} ($>5\sigma$) and {\it Chandra} ($>4.5\sigma$). This catalog release includes data from {\it XMM-Newton} cycle AO 13, which approximately doubled the Stripe 82X survey area. The flux limits of the Stripe 82X survey are $8.7\times10^{-16}$ erg s$^{-1}$ cm$^{-2}$, $4.7\times10^{-15}$ erg s$^{-1}$ cm$^{-2}$, and $2.1\times10^{-15}$ erg s$^{-1}$ cm$^{-2}$ in the soft (0.5-2 keV), hard (2-10 keV), and full bands (0.5-10 keV), respectively, with approximate half-area survey flux limits of $5.4\times10^{-15}$ erg s$^{-1}$ cm$^{-2}$, $2.9\times10^{-14}$ erg s$^{-1}$ cm$^{-2}$, and $1.7\times10^{-14}$ erg s$^{-1}$ cm$^{-2}$. We matched the X-ray source lists to available multi-wavelength catalogs, including updated matches to the previous release of the Stripe 82X survey; 88\% of the sample is matched to a multi-wavelength counterpart. Due to the wide area of Stripe 82X and rich ancillary multi-wavelength data, including coadded SDSS photometry, mid-infrared {\it WISE} coverage, near-infrared coverage from UKIDSS and VHS, ultraviolet coverage from {\it GALEX}, radio coverage from FIRST, and far-infrared coverage from {\it Herschel}, as well as existing $\sim$30\% optical spectroscopic completeness, we are beginning to uncover rare objects, such as obscured high-luminosity AGN at high-redshift. The Stripe 82X point source catalog is a valuable dataset for constraining how this population grows and evolves, as well as for studying how they interact with the galaxies in which they live.

Heavily obscured, Compton Thick (CT, NH>10^24 cm^-2) AGN may represent an important phase in AGN/galaxy co-evolution and are expected to provide a significant contribution to the cosmic X-ray background (CXB). Through direct X-ray spectra analysis, we selected 39 heavily obscured AGN (NH>3x10^23 cm^-2) in the 2 deg^2 XMM-COSMOS survey. After selecting CT AGN based on the fit of a simple absorbed two power law model to the XMM data, the presence of CT AGN was confirmed in 80% of the sources using deeper Chandra data and more complex models. The final sample of CT AGN comprises 10 sources spanning a large range of redshift and luminosity. We collected the multi-wavelength information available for all these sources, in order to study the distribution of SMBH and host properties, such as BH mass (M_BH), Eddington ratio (\lambda_Edd), stellar mass (M*), specific star formation rate (sSFR) in comparison with a sample of unobscured AGN. We find that highly obscured sources tend to have significantly smaller M_BH and higher \lambda_edd with respect to unobscured ones, while a weaker evolution in M* is observed. The sSFR of highly obscured sources is consistent with the one observed in the main sequence of star forming galaxies, at all redshift. We also present optical spectra, spectral energy distribution (SED) and morphology for the sample of 10 CT AGN: all the available optical spectra are dominated by the stellar component of the host galaxy, and a highly obscured torus component is needed in the SED of the CT sources. Exploiting the high resolution Hubble-ACS images available, we conclude that these highly obscured sources have a significantly larger merger fraction with respect to other X-ray selected samples of AGN. Finally we discuss implications in the context of AGN/galaxy co-evolutionary models, and compare our results with the predictions of CXB synthesis models.

We present a detailed spectral analysis of new XMM-Newton data of the source CXOC J100043.1+020637, also known as CID-42, detected in the COSMOS survey at z = 0.359. Previous works suggested that CID-42 is a candidate recoiling supermassive black holes showing also an inverted P-Cygni profile in the X- ray spectra at ~6 keV (rest) with an iron emission line plus a redshifted absorption line (detected at 3sigma in previous XMM-Newton and Chandra observations). Detailed analysis of the absorption line suggested the presence of ionized material inflowing into the black hole at high velocity. In the new long XMM-Newton observation, while the overall spectral shape remains constant, the continuum 2-10 keV flux decreased of ~20% with respect to previous observation and the absorption line is undetected. The upper limit on the intensity of the absorption line is EW<162 keV. Extensive Monte Carlo simulations show that the non detection of the line is solely due to variation in the properties of the inflowing material, in agreement with the transient nature of these features, and that the intensity of the line is lower than the previously measured with a probability of 98.8%. In the scenario of CID-42 as recoiling SMBH, the absorption line can be interpreted as due to inflow of gas with variable density and located in the proximity of the SMBH and recoiling with it. New monitoring observations will be requested to further characterize this line.