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 examine morphology-separated color-mass diagrams to study the quenching of star formation in $\sim 100,000$ ($z\sim0$) Sloan Digital Sky Survey (SDSS) and $\sim 20,000$ ($z\sim1$) Cosmic Assembly Near-Infrared Deep Extragalactic Legacy Survey (CANDELS) galaxies. To classify galaxies morphologically, we developed Galaxy Morphology Network (GaMorNet), a convolutional neural network that classifies galaxies according to their bulge-to-total light ratio. GaMorNet does not need a large training set of real data and can be applied to data sets with a range of signal-to-noise ratios and spatial resolutions. GaMorNet's source code as well as the trained models are made public as part of this work ( http://www.astro.yale.edu/aghosh/gamornet.html ). We first trained GaMorNet on simulations of galaxies with a bulge and a disk component and then transfer learned using $\sim25\%$ of each data set to achieve misclassification rates of $\lesssim5\%$. The misclassified sample of galaxies is dominated by small galaxies with low signal-to-noise ratios. Using the GaMorNet classifications, we find that bulge- and disk-dominated galaxies have distinct color-mass diagrams, in agreement with previous studies. For both SDSS and CANDELS galaxies, disk-dominated galaxies peak in the blue cloud, across a broad range of masses, consistent with the slow exhaustion of star-forming gas with no rapid quenching. A small population of red disks is found at high mass ($\sim14\%$ of disks at $z\sim0$ and $2\%$ of disks at $z \sim 1$). In contrast, bulge-dominated galaxies are mostly red, with much smaller numbers down toward the blue cloud, suggesting rapid quenching and fast evolution across the green valley. This inferred difference in quenching mechanism is in agreement with previous studies that used other morphology classification techniques on much smaller samples at $z\sim0$ and $z\sim1$.

Theory predicts that a supermassive black hole binary (SMBHB) could be observed as a luminous active galactic nucleus (AGN) that periodically varies on the order of its orbital timescale. In X-rays, periodic variations could be caused by mechanisms including relativistic Doppler boosting and shocks. Here we present the first systematic search for periodic AGNs using $941$ hard X-ray light curves (14-195 keV) from the first 105 months of the Swift Burst Alert Telescope (BAT) survey (2004-2013). We do not find evidence for periodic AGNs in Swift-BAT, including the previously reported SMBHB candidate MCG+11$-$11$-$032. We find that the null detection is consistent with the combination of the upper-limit binary population in AGNs in our adopted model, their expected periodic variability amplitudes, and the BAT survey characteristics. We have also investigated the detectability of SMBHBs against normal AGN X-ray variability in the context of the eROSITA survey. Under our assumptions of a binary population and the periodic signals they produce which have long periods of hundreds of days, up to $13$% true periodic binaries can be robustly distinguished from normal variable AGNs with the ideal uniform sampling. However, we demonstrate that realistic eROSITA sampling is likely to be insensitive to long-period binaries because longer observing gaps reduce their detectability. In contrast, large observing gaps do not diminish the prospect of detecting binaries of short, few-day periods, as 19% can be successfully recovered, the vast majority of which can be identified by the first half of the survey.

The clustering of active galactic nuclei (AGN) sheds light on their typical large (Mpc-scale) environments, which can constrain the growth and evolution of supermassive black holes. Here we measure the clustering of luminous X-ray-selected AGN in the Stripe 82X and XMM-XXL-North surveys around the peak epoch of black hole growth, in order to investigate the dependence of luminosity on large-scale AGN environment. We compute the auto-correlation function of AGN in two luminosity bins, $10^{43}\leq L_X<10^{44.5}$ erg s$^{-1}$ at $z\sim 0.8$ and $L_X\geq 10^{44.5}$ erg s$^{-1}$ at $z\sim 1.8$, and calculate the AGN bias taking into account the redshift distribution of the sources using three different methods. Our results show that while the less luminous sample has an inferred typical halo mass that is smaller than for the more luminous AGN, the host halo mass may be less dependent on luminosity than suggested in previous work. Focusing on the luminous sample, we calculate a typical host halo mass of $\sim 10^{13}$ M$_{\odot}~h^{-1}$, which is similar to previous measurements of moderate-luminosity X-ray AGN and significantly larger than the values found for optical quasars of similar luminosities and redshifts. We suggest that the clustering differences between different AGN selection techniques are dominated by selection biases, and not due to a dependence on AGN luminosity. We discuss the limitations of inferring AGN triggering mechanisms from halo masses derived by large-scale bias.

We have conducted 22 GHz radio imaging at 1" resolution of 100 low-redshift AGN selected at 14-195 keV by the Swift-BAT. We find a radio core detection fraction of 96%, much higher than lower-frequency radio surveys. Of the 96 radio-detected AGN, 55 have compact morphologies, 30 have morphologies consistent with nuclear star formation, and 11 have sub-kpc to kpc-scale jets. We find that the total radio power does not distinguish between nuclear star formation and jets as the origin of the radio emission. For 87 objects, we use optical spectroscopy to test whether AGN physical parameters are distinct between radio morphological types. We find that X-ray luminosities tend to be higher if the 22 GHz morphology is jet-like, but find no significant difference in other physical parameters. We find that the relationship between the X-ray and core radio luminosities is consistent with the $L_R/L_X \sim 10^{-5}$ of coronally active stars. We further find that the canonical fundamental planes of black hole activity systematically over-predict our radio luminosities, particularly for objects with star formation morphologies.

We present the highest resolution --- 15 pc (0.03'') --- ALMA $^{12}$CO(2-1) line emission and 1.3mm continuum maps, tracers of the molecular gas and dust, respectively, in the nearby merging galaxy system NGC 6240, that hosts two supermassive black holes growing simultaneously. These observations provide an excellent spatial match to existing Hubble optical and near-infrared observations of this system. A significant molecular gas mass, $\sim$9$\times$10$^9$M$_\odot$, is located in between the two nuclei, forming a clumpy stream kinematically dominated by turbulence, rather than a smooth rotating disk as previously assumed from lower resolution data. Evidence for rotation is seen in the gas surrounding the southern nucleus, but not in the northern one. Dynamical shells can be seen, likely associated with nuclear supernovae remnants. We further detect the presence of significant high velocity outflows, some of them reaching velocities $>$500 km/s, affecting a significant fraction, $\sim$11\% of the molecular gas in the nuclear region. Inside the spheres of influence of the northern and southern supermassive black holes we find molecular masses of 7.4$\times$10$^8$M$_\odot$ and 3.3$\times$10$^9$M$_\odot$, respectively. We are thus directly imaging the reservoir of gas that can accrete onto each supermassive black hole. These new ALMA maps highlight the critical need for high resolution observations of molecular gas in order to understand the feeding of supermassive black holes and its connection to galaxy evolution in the context of a major galaxy merger.

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 cross-correlation studies of gamma-ray (0.1-300 GeV), X-ray (0.2-10 keV) and optical (R-band) variability of a sample of 26 blazars during 2008-2016. The light curves are from Fermi-LAT, Swift-XRT, and the Yale-SMARTS blazar monitoring program. We stack the discrete cross-correlation functions of the blazars such that the features that are consistently present in a large fraction of the sample become more prominent in the final result. We repeat the same analysis for two subgroups, namely, low synchrotron peaked (LSP) and high synchrotron peaked (HSP) blazars. We find that, on average, the variability at multiple bands is correlated, with a time lag consistent with zero in both subgroups. We describe this correlation with a leptonic model of non-thermal emission from blazar jets. By comparing the model results with those from the actual data we find that the inter-band cross-correlations are consistent with an emission region of size nearly 0.1 pc within the broad line region for LSP blazars. We rule out large changes of magnetic field (> 0.5 Gauss) across the emission region or small values of magnetic field (e.g. 0.2 Gauss) for this population. We also find that the observed variability of the HSP blazars can be explained if the emission region is much larger than the distance to the broad line region from the central black hole.

With the launch of eROSITA (extended Roentgen Survey with an Imaging Telescope Array), successfully occurred on 2019 July 13, we are facing the challenge of computing reliable photometric redshifts for 3 million of active galactic nuclei (AGNs) over the entire sky, having available only patchy and inhomogeneous ancillary data. While we have a good understanding of the photo-z quality obtainable for AGN using spectral energy distribution (SED)-fitting technique, we tested the capability of machine learning (ML), usually reliable in computing photo-z for QSO in wide and shallow areas with rich spectroscopic samples. Using MLPQNA as example of ML, we computed photo-z for the X-ray-selected sources in Stripe 82X, using the publicly available photometric and spectroscopic catalogues. Stripe 82X is at least as deep as eROSITA will be and wide enough to include also rare and bright AGNs. In addition, the availability of ancillary data mimics what can be available in the whole sky. We found that when optical, and near- and mid-infrared data are available, ML and SED fitting perform comparably well in terms of overall accuracy, realistic redshift probability density functions, and fraction of outliers, although they are not the same for the two methods. The results could further improve if the photometry available is accurate and including morphological information. Assuming that we can gather sufficient spectroscopy to build a representative training sample, with the current photometry coverage we can obtain reliable photo-z for a large fraction of sources in the Southern hemisphere well before the spectroscopic follow-up, thus timely enabling the eROSITA science return. The photo-z catalogue is released here.

We present a multi wavelength analysis of 28 of the most luminous low-redshift narrow-line, ultra-hard X-ray selected active galactic nuclei (AGN) drawn from the 70 month Swift/BAT all-sky survey, with bolometric luminosities of log(L_bol/erg/s) > 45.25. The broad goal of our study is to determine whether these objects have any distinctive properties, potentially setting them aside from lower-luminosity obscured AGN in the local Universe. Our analysis relies on the first data release of the BAT AGN Spectroscopic Survey (BASS/DR1) and on dedicated observations with the VLT, Palomar, and Keck observatories. We find that the vast majority of our sources agree with commonly used AGN selection criteria which are based on emission line ratios and on mid-infrared colours. Our AGN are predominantly hosted in massive galaxies (9.8 < log(M_*/M_sun) < 11.7); based on visual inspection of archival optical images, they appear to be mostly ellipticals. Otherwise, they do not have distinctive properties. Their radio luminosities, determined from publicly available survey data, show a large spread of almost 4 orders of magnitude - much broader than what is found for lower X-ray luminosity obscured AGN in BASS. Moreover, our sample shows no preferred combination of black hole masses (M_BH) and/or Eddington ratio (lambda_Edd), covering 7.5 < log(M_BH/M_sun) < 10.3 and 0.01 < lambda_Edd < 1. Based on the distribution of our sources in the lambda_Edd-N_H plane, we conclude that our sample is consistent with a scenario where the amount of obscuring material along the line of sight is determined by radiation pressure exerted by the AGN on the dusty circumnuclear gas.

The recently released 105-month {\it Swift}-Burst Alert Telescope (BAT) all-sky hard X-ray survey catalog presents an opportunity to study astrophysical objects detected in the deepest look at the entire hard X-ray (14$-$195 keV) sky. Here we report the results of a multifrequency study of 146 blazars from this catalog, quadrupling the number compared to past studies, by utilizing recent data from the {\it Fermi}-Large Area Telescope (LAT), Swift-BAT, and archival measurements. In our $\gamma$-ray analysis of $\sim$10 years of the LAT data, 101 are found as $\gamma$-ray emitters, whereas, 45 remains LAT undetected. We model the broadband spectral energy distributions with a synchrotron-inverse Compton radiative model. On average, BAT detected sources host massive black holes ($M_{\rm bh}\sim10^9$ M$_{\odot}$) and luminous accretion disks ($L_{\rm d}\sim10^{46}$ erg s$^{-1}$). At high-redshifts ($z>2$), BAT blazars host more powerful jets with luminous accretion disks compared to those detected only with the {\it Fermi}-LAT. We find good agreement in the black hole masses derived from the single-epoch optical spectroscopic measurements and standard accretion disk modeling approaches. Other physical properties of BAT blazars are similar to those known for {\it Fermi}-LAT detected objects.

The origin and properties of black hole seeds that grow to produce the detected population of supermassive black holes are unconstrained at present. Despite the existence of several potentially feasible channels for the production of initial seeds in the high redshift universe, since even actively growing seeds are not directly observable at these epochs, discriminating between models remains challenging. Several new observables that encapsulate information about seeding have been proposed in recent years, and these offer exciting prospects for truly unraveling the nature of black hole seeds in the coming years. One of the key challenges for this task lies in the complexity of the problem, the required disentangling of the confounding effects of accretion physics and mergers, as mergers and accretion events over cosmic time stand to erase these initial conditions. Nevertheless, some unique signatures of seeding do survive and still exist in: local scaling relations between black holes and their galaxy hosts at low-masses; in high-redshift luminosity functions of accreting black holes; and in the total number and mass functions of gravitational wave coalescence events from merging binary black holes. One of the clearest discriminants for seed models are these high redshift gravitational wave detections of mergers from space detectable in the milliHertz range. These predicted event rates offer the most direct constraints on the properties of initial black hole seeds. Improving our theoretical understanding of black hole dynamics and accretion will also be pivotal in constraining seeding models in combination with the wide range of multi-messenger data.

We report the results of a Sloan Digital Sky Survey-IV eBOSS program to target X-ray sources and mid-infrared-selected WISE AGN candidates in a 36.8 deg$^2$ region of Stripe 82. About half this survey (15.6 deg$^2$) covers the largest contiguous portion of the Stripe 82 X-ray survey. This program represents the largest spectroscopic survey of AGN candidates selected solely by their WISE colors. We combine this sample with X-ray and WISE AGN in the field identified via other sources of spectroscopy, producing a catalog of 4847 sources that is 82% complete to $r\sim22$. Based on X-ray luminosities or WISE colors, 4730 of these sources are AGN, with a median sample redshift of $z\sim1$. About 30% of the AGN are optically obscured (i.e., lack broad lines in their optical spectra). BPT analysis, however, indicates that 50% of the WISE AGN at $z<0.5$ have emission line ratios consistent with star-forming galaxies, so whether they are buried AGN or star-forming galaxy contaminants is currently unclear. We find that 61% of X-ray AGN are not selected as MIR AGN, with 22% of X-ray AGN undetected by WISE. Most of these latter AGN have high X-ray luminosities ($L_{\rm x} > 10^{44}$ erg s$^{-1}$), indicating that MIR selection misses a sizable fraction of the highest luminosity AGN, as well as lower luminosity sources where AGN heated dust is not dominating the MIR emission. Conversely, $\sim$58% of WISE AGN are undetected by X-rays, though we do not find that they are preferentially redder than the X-ray detected WISE AGN.

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.

Aims: We probe the radiatively-efficient, hot wind feedback mode in two nearby luminous unobscured (type 1) AGN from the Close AGN Reference Survey (CARS), which show intriguing kpc-scale arc-like features of extended [OIII] ionized gas as mapped with VLT-MUSE. We aimed to detect hot gas bubbles that would indicate the existence of powerful, galaxy-scale outflows in our targets, HE 0227-0931 and HE 0351+0240, from deep (200 ks) Chandra observations. Methods: By measuring the spatial and spectral properties of the extended X-ray emission and comparing with the sub kpc-scale IFU data, we are able to constrain feedback scenarios and directly test if the ionized gas is due to a shocked wind. Results: No extended hot gas emission on kpc-scales was detected. Unless the ambient medium density is low ($n_{H}\sim~1$ cm$^{-3}$ at 100 pc), the inferred upper limits on the extended X-ray luminosities are well below what is expected from theoretical models at matching AGN luminosities. Conclusions: We conclude that the highly-ionized gas structures on kpc scales are not inflated by a hot outflow in either target, and instead are likely caused by photo-ionization of pre-existing gas streams of different origins. Our non-detections suggest that extended X-ray emission from an AGN-driven wind is not universal, and may lead to conflicts with current theoretical predictions.

We present the results from the Hitomi Soft Gamma-ray Detector (SGD) observation of the Crab nebula. The main part of SGD is a Compton camera, which in addition to being a spectrometer, is capable of measuring polarization of gamma-ray photons. The Crab nebula is one of the brightest X-ray / gamma-ray sources on the sky, and, the only source from which polarized X-ray photons have been detected. SGD observed the Crab nebula during the initial test observation phase of Hitomi. We performed the data analysis of the SGD observation, the SGD background estimation and the SGD Monte Carlo simulations, and, successfully detected polarized gamma-ray emission from the Crab nebula with only about 5 ks exposure time. The obtained polarization fraction of the phase-integrated Crab emission (sum of pulsar and nebula emissions) is (22.1 $\pm$ 10.6)% and, the polarization angle is 110.7$^o$ + 13.2 / $-$13.0$^o$ in the energy range of 60--160 keV (The errors correspond to the 1 sigma deviation). The confidence level of the polarization detection was 99.3%. The polarization angle measured by SGD is about one sigma deviation with the projected spin axis of the pulsar, 124.0$^o$ $\pm$0.1$^o$.

We present ALMA and MUSE observations of the Brightest Cluster Galaxy in Abell 2597, a nearby (z=0.0821) cool core cluster of galaxies. The data map the kinematics of a three billion solar mass filamentary nebula that spans the innermost 30 kpc of the galaxy's core. Its warm ionized and cold molecular components are both cospatial and comoving, consistent with the hypothesis that the optical nebula traces the warm envelopes of many cold molecular clouds that drift in the velocity field of the hot X-ray atmosphere. The clouds are not in dynamical equilibrium, and instead show evidence for inflow toward the central supermassive black hole, outflow along the jets it launches, and uplift by the buoyant hot bubbles those jets inflate. The entire scenario is therefore consistent with a galaxy-spanning "fountain", wherein cold gas clouds drain into the black hole accretion reservoir, powering jets and bubbles that uplift a cooling plume of low-entropy multiphase gas, which may stimulate additional cooling and accretion as part of a self-regulating feedback loop. All velocities are below the escape speed from the galaxy, and so these clouds should rain back toward the galaxy center from which they came, keeping the fountain long-lived. The data are consistent with major predictions of chaotic cold accretion, precipitation, and stimulated feedback models, and may trace processes fundamental to galaxy evolution at effectively all mass scales.

The Shocked POststarburst Galaxy Survey (SPOGS) aims to identify galaxies in the transitional phase between actively star-forming and quiescence with nebular lines that are excited from shocks rather than star formation processes. We explored the ultraviolet (UV) properties of objects with near-ultraviolet (NUV) and far-ultraviolet (FUV) photometry from archival GALEX data; 444 objects were detected in both bands, 365 in only NUV, and 24 in only FUV, for a total of 833 observed objects. We compared SPOGs to samples of Star-forming galaxies (SFs), Quiescent galaxies (Qs), classical E+A post-starburst galaxies, active galactic nuclei (AGN) host galaxies, and interacting galaxies. We found that SPOGs have a larger range in their FUV-NUV and NUV-r colors compared to most of the other samples, although all of our comparison samples occupied color space inside of the SPOGs region. Based on their UV colors, SPOGs are a heterogeneous group, possibly made up of a mixture of SFs, Qs, and/or AGN. Using Gaussian mixture models, we are able to recreate the distribution of FUV-NUV colors of SPOGs and E+A galaxies with different combinations of SFs, Qs, and AGN. We find that the UV colors of SPOGs require a >60% contribution from SFs, with either Qs or AGN representing the remaining contribution, while UV colors of E+A galaxies required a significantly lower fraction of SFs, supporting the idea that SPOGs are at an earlier point in their transition from quiescent to star-forming than E+A galaxies.

We present a spectroscopically complete sample of 147 infrared-color-selected AGN down to a 22 $\mu$m flux limit of 20 mJy over the $\sim$270 deg$^2$ of the SDSS Stripe 82 region. Most of these sources are in the QSO luminosity regime ($L_{\rm bol} \gtrsim 10^{12} L_\odot$) and are found out to $z\simeq3$. We classify the AGN into three types, finding: 57 blue, unobscured Type-1 (broad-lined) sources; 69 obscured, Type-2 (narrow-lined) sources; and 21 moderately-reddened Type-1 sources (broad-lined and $E(B-V) > 0.25$). We study a subset of this sample in X-rays and analyze their obscuration to find that our spectroscopic classifications are in broad agreement with low, moderate, and large amounts of absorption for Type-1, red Type-1 and Type-2 AGN, respectively. We also investigate how their X-ray luminosities correlate with other known bolometric luminosity indicators such as [O III] line luminosity ($L_{\rm [OIII]}$) and infrared luminosity ($L_{6 \mu{\rm m}}$). While the X-ray correlation with $L_{\rm [OIII]}$ is consistent with previous findings, the most infrared-luminous sources appear to deviate from established relations such that they are either under-luminous in X-rays or over-luminous in the infrared. Finally, we examine the luminosity function (LF) evolution of our sample, and by AGN type, in combination with the complementary, infrared-selected, AGN sample of Lacy et al. (2013), spanning over two orders of magnitude in luminosity. We find that the two obscured populations evolve differently, with reddened Type-1 AGN dominating the obscured AGN fraction ($\sim$30%) for $L_{5 \mu{\rm m}} > 10^{45}$ erg s$^{-1}$, while the fraction of Type-2 AGN with $L_{5 \mu{\rm m}} < 10^{45}$ erg s$^{-1}$ rises sharply from 40% to 80% of the overall AGN population.

We characterize the environments of local accreting supermassive black holes by measuring the clustering of AGN in the Swift/BAT Spectroscopic Survey (BASS). With 548 AGN in the redshift range 0.01<z<0.1 over the full sky from the DR1 catalog, BASS provides the largest, least biased sample of local AGN to date due to its hard X-ray selection (14-195 keV) and rich multiwavelength/ancillary data. By measuring the projected cross-correlation function between the AGN and 2MASS galaxies, and interpreting it via halo occupation distribution (HOD) and subhalo-based models, we constrain the occupation statistics of the full sample, as well as in bins of absorbing column density and black hole mass. We find that AGN tend to reside in galaxy group environments, in agreement with previous studies of AGN throughout a large range of luminosity and redshift, and that on average they occupy their dark matter halos similar to inactive galaxies of comparable stellar mass. We also find evidence that obscured AGN tend to reside in denser environments than unobscured AGN, even when samples were matched in luminosity, redshift, stellar mass, and Eddington ratio. We show that this can be explained either by significantly different halo occupation distributions or statistically different host halo assembly histories. Lastly, we see that massive black holes are slightly more likely to reside in central galaxies than black holes of smaller mass.

We present a framework to link and describe AGN variability on a wide range of timescales, from days to billions of years. In particular, we concentrate on the AGN variability features related to changes in black hole fuelling and accretion rate. In our framework, the variability features observed in different AGN at different timescales may be explained as realisations of the same underlying statistical properties. In this context, we propose a model to simulate the evolution of AGN light curves with time based on the probability density function (PDF) and power spectral density (PSD) of the Eddington ratio ($L/L_{\rm Edd}$) distribution. Motivated by general galaxy population properties, we propose that the PDF may be inspired by the $L/L_{\rm Edd}$ distribution function (ERDF), and that a single (or limited number of) ERDF+PSD set may explain all observed variability features. After outlining the framework and the model, we compile a set of variability measurements in terms of structure function (SF) and magnitude difference. We then combine the variability measurements on a SF plot ranging from days to Gyr. The proposed framework enables constraints on the underlying PSD and the ability to link AGN variability on different timescales, therefore providing new insights into AGN variability and black hole growth phenomena.

We present results from the Hitomi X-ray observation of a young composite-type supernova remnant (SNR) G21.5$-$0.9, whose emission is dominated by the pulsar wind nebula (PWN) contribution. The X-ray spectra in the 0.8-80 keV range obtained with the Soft X-ray Spectrometer (SXS), Soft X-ray Imager (SXI) and Hard X-ray Imager (HXI) show a significant break in the continuum as previously found with the NuSTAR observation. After taking into account all known emissions from the SNR other than the PWN itself, we find that the Hitomi spectra can be fitted with a broken power law with photon indices of $\Gamma_1=1.74\pm0.02$ and $\Gamma_2=2.14\pm0.01$ below and above the break at $7.1\pm0.3$ keV, which is significantly lower than the NuSTAR result ($\sim9.0$ keV). The spectral break cannot be reproduced by time-dependent particle injection one-zone spectral energy distribution models, which strongly indicates that a more complex emission model is needed, as suggested by recent theoretical models. We also search for narrow emission or absorption lines with the SXS, and perform a timing analysis of PSR J1833$-$1034 with the HXI and SGD. No significant pulsation is found from the pulsar. However, unexpectedly, narrow absorption line features are detected in the SXS data at 4.2345 keV and 9.296 keV with a significance of 3.65 $\sigma$. While the origin of these features is not understood, their mere detection opens up a new field of research and was only possible with the high resolution, sensitivity and ability to measure extended sources provided by an X-ray microcalorimeter.

We present optical and near-IR Integral Field Unit (IFU) and ALMA band 6 observations of the nearby dual Active Galactic Nuclei (AGN) Mrk 463. At a distance of 210 Mpc, and a nuclear separation of $\sim$4 kpc, Mrk 463 is an excellent laboratory to study the gas dynamics, star formation processes and supermassive black hole (SMBH) accretion in a late-stage gas-rich major galaxy merger. The IFU observations reveal a complex morphology, including tidal tails, star-forming clumps, and emission line regions. The optical data, which map the full extent of the merger, show evidence for a biconical outflow and material outflowing at $>$600 km s$^{-1}$, both associated with the Mrk 463E nucleus, together with large scale gradients likely related to the ongoing galaxy merger. We further find an emission line region $\sim$11 kpc south of Mrk 463E that is consistent with being photoionized by an AGN. Compared to the current AGN luminosity, the energy budget of the cloud implies a luminosity drop in Mrk 463E by a factor 3-20 over the last 40,000 years. The ALMA observations of $^{12}$CO(2-1) and adjacent 1mm continuum reveal the presence of $\sim$10$^{9}$M$_\odot$ in molecular gas in the system. The molecular gas shows velocity gradients of $\sim$800 km/s and $\sim$400 km/s around the Mrk 463E and 463W nuclei, respectively. We conclude that in this system the infall of $\sim$100s $M_\odot$/yr of molecular gas is in rough balance with the removal of ionized gas by a biconical outflow being fueled by a relatively small, $<$0.01% of accretion onto each SMBH.

The present paper investigates the temperature structure of the X-ray emitting plasma in the core of the Perseus cluster using the 1.8--20.0 keV data obtained with the Soft X-ray Spectrometer (SXS) onboard the Hitomi Observatory. A series of four observations were carried out, with a total effective exposure time of 338 ks and covering a central region $\sim7'$ in diameter. The SXS was operated with an energy resolution of $\sim$5 eV (full width at half maximum) at 5.9 keV. Not only fine structures of K-shell lines in He-like ions but also transitions from higher principal quantum numbers are clearly resolved from Si through Fe. This enables us to perform temperature diagnostics using the line ratios of Si, S, Ar, Ca, and Fe, and to provide the first direct measurement of the excitation temperature and ionization temperature in the Perseus cluster. The observed spectrum is roughly reproduced by a single temperature thermal plasma model in collisional ionization equilibrium, but detailed line ratio diagnostics reveal slight deviations from this approximation. In particular, the data exhibit an apparent trend of increasing ionization temperature with increasing atomic mass, as well as small differences between the ionization and excitation temperatures for Fe, the only element for which both temperatures can be measured. The best-fit two-temperature models suggest a combination of 3 and 5 keV gas, which is consistent with the idea that the observed small deviations from a single temperature approximation are due to the effects of projection of the known radial temperature gradient in the cluster core along the line of sight. Comparison with the Chandra/ACIS and the XMM-Newton/RGS results on the other hand suggests that additional lower-temperature components are present in the ICM but not detectable by Hitomi SXS given its 1.8--20 keV energy band.

The Hitomi SXS spectrum of the Perseus cluster, with $\sim$5 eV resolution in the 2-9 keV band, offers an unprecedented benchmark of the atomic modeling and database for hot collisional plasmas. It reveals both successes and challenges of the current atomic codes. The latest versions of AtomDB/APEC (3.0.8), SPEX (3.03.00), and CHIANTI (8.0) all provide reasonable fits to the broad-band spectrum, and are in close agreement on best-fit temperature, emission measure, and abundances of a few elements such as Ni. For the Fe abundance, the APEC and SPEX measurements differ by 16%, which is 17 times higher than the statistical uncertainty. This is mostly attributed to the differences in adopted collisional excitation and dielectronic recombination rates of the strongest emission lines. We further investigate and compare the sensitivity of the derived physical parameters to the astrophysical source modeling and instrumental effects. The Hitomi results show that an accurate atomic code is as important as the astrophysical modeling and instrumental calibration aspects. Substantial updates of atomic databases and targeted laboratory measurements are needed to get the current codes ready for the data from the next Hitomi-level mission.

We present Hitomi observations of N132D, a young, X-ray bright, O-rich core-collapse supernova remnant in the Large Magellanic Cloud (LMC). Despite a very short observation of only 3.7 ks, the Soft X-ray Spectrometer (SXS) easily detects the line complexes of highly ionized S K and Fe K with 16-17 counts in each. The Fe feature is measured for the first time at high spectral resolution. Based on the plausible assumption that the Fe K emission is dominated by He-like ions, we find that the material responsible for this Fe emission is highly redshifted at ~800 km/s compared to the local LMC interstellar medium (ISM), with a 90% credible interval of 50-1500 km/s if a weakly informative prior is placed on possible line broadening. This indicates (1) that the Fe emission arises from the supernova ejecta, and (2) that these ejecta are highly asymmetric, since no blue-shifted component is found. The S K velocity is consistent with the local LMC ISM, and is likely from swept-up ISM material. These results are consistent with spatial mapping that shows the He-like Fe concentrated in the interior of the remnant and the S tracing the outer shell. The results also show that even with a very small number of counts, direct velocity measurements from Doppler-shifted lines detected in extended objects like supernova remnants are now possible. Thanks to the very low SXS background of ~1 event per spectral resolution element per 100 ks, such results are obtainable during short pointed or slew observations with similar instruments. This highlights the power of high-spectral-resolution imaging observations, and demonstrates the new window that has been opened with Hitomi and will be greatly widened with future missions such as the X-ray Astronomy Recovery Mission (XARM) and Athena.

Super-massive black holes weighing up to $\sim 10^9 \, \mathrm{M_{\odot}}$ are in place by $z \sim 7$, when the age of the Universe is $\lesssim 1 \, \mathrm{Gyr}$. This implies a time crunch for their growth, since such high masses cannot be easily reached in standard accretion scenarios. Here, we explore the physical conditions that would lead to optimal growth wherein stable super-Eddington accretion would be permitted. Our analysis suggests that the preponderance of optimal conditions depends on two key parameters: the black hole mass and the host galaxy central gas density. In the high-efficiency region of this parameter space, a continuous stream of gas can accrete onto the black hole from large to small spatial scales, assuming a global isothermal profile for the host galaxy. Using analytical initial mass functions for black hole seeds, we find an enhanced probability of high-efficiency growth for seeds with initial masses $\gtrsim 10^4 \, \mathrm{M_{\odot}}$. Our picture suggests that a large population of high-$z$ lower-mass black holes that formed in the low-efficiency region, with low duty cycles and accretion rates, might remain undetectable as quasars, since we predict their bolometric luminosities to be $\lesssim 10^{41} \, \mathrm{erg \, s^{-1}}$. The presence of these sources might be revealed only via gravitational wave detections of their mergers.

We report a Hitomi observation of IGR J16318-4848, a high-mass X-ray binary system with an extremely strong absorption of N_H~10^{24} cm^{-2}. Previous X-ray studies revealed that its spectrum is dominated by strong fluorescence lines of Fe as well as continuum emission. For physical and geometrical insight into the nature of the reprocessing material, we utilize the high spectroscopic resolving power of the X-ray microcalorimeter (the soft X-ray spectrometer; SXS) and the wide-band sensitivity by the soft and hard X-ray imager (SXI and HXI) aboard Hitomi. Even though photon counts are limited due to unintended off-axis pointing, the SXS spectrum resolves Fe K{\alpha_1} and K{\alpha_2} lines and puts strong constraints on the line centroid and width. The line width corresponds to the velocity of 160^{+300}_{-70} km s^{-1}. This represents the most accurate, and smallest, width measurement of this line made so far from any X-ray binary, much less than the Doppler broadening and shift expected from speeds which are characteristic of similar systems. Combined with the K-shell edge energy measured by the SXI and HXI spectra, the ionization state of Fe is estimated to be in the range of Fe I--IV. Considering the estimated ionization parameter and the distance between the X-ray source and the absorber, the density and thickness of the materials are estimated. The extraordinarily strong absorption and the absence of a Compton shoulder component is confirmed. These characteristics suggest reprocessing materials which are distributed in a narrow solid angle or scattering primarily with warm free electrons or neutral hydrogen.

The origin of the narrow Fe-K{\alpha} fluorescence line at 6.4 keV from active galactic nuclei has long been under debate; some of the possible sites are the outer accretion disk, the broad line region, a molecular torus, or interstellar/intracluster media. In February-March 2016, we performed the first X-ray microcalorimeter spectroscopy with the Soft X-ray Spectrometer (SXS) onboard the Hitomi satellite of the Fanaroff-Riley type I radio galaxy NGC 1275 at the center of the Perseus cluster of galaxies. With the high energy resolution of ~5 eV at 6 keV achieved by Hitomi/SXS, we detected the Fe-K{\alpha} line with ~5.4 {\sigma} significance. The velocity width is constrained to be 500-1600 km s$^{-1}$ (FWHM for Gaussian models) at 90% confidence. The SXS also constrains the continuum level from the NGC 1275 nucleus up to ~20 keV, giving an equivalent width ~20 eV of the 6.4 keV line. Because the velocity width is narrower than that of broad H{\alpha} line of ~2750 km s$^{-1}$, we can exclude a large contribution to the line flux from the accretion disk and the broad line region. Furthermore, we performed pixel map analyses on the Hitomi/SXS data and image analyses on the Chandra archival data, and revealed that the Fe-K{\alpha} line comes from a region within ~1.6 kpc from the NGC 1275 core, where an active galactic nucleus emission dominates, rather than that from intracluster media. Therefore, we suggest that the source of the Fe-K{\alpha} line from NGC 1275 is likely a low-covering fraction molecular torus or a rotating molecular disk which probably extends from a pc to hundreds pc scale in the active galactic nucleus system.

Extending the earlier measurements reported in Hitomi collaboration (2016, Nature, 535, 117), we examine the atmospheric gas motions within the central 100~kpc of the Perseus cluster using observations obtained with the Hitomi satellite. After correcting for the point spread function of the telescope and using optically thin emission lines, we find that the line-of-sight velocity dispersion of the hot gas is remarkably low and mostly uniform. The velocity dispersion reaches maxima of approximately 200~km~s$^{-1}$ toward the central active galactic nucleus (AGN) and toward the AGN inflated north-western `ghost' bubble. Elsewhere within the observed region, the velocity dispersion appears constant around 100~km~s$^{-1}$. We also detect a velocity gradient with a 100~km~s$^{-1}$ amplitude across the cluster core, consistent with large-scale sloshing of the core gas. If the observed gas motions are isotropic, the kinetic pressure support is less than 10\% of the thermal pressure support in the cluster core. The well-resolved optically thin emission lines have Gaussian shapes, indicating that the turbulent driving scale is likely below 100~kpc, which is consistent with the size of the AGN jet inflated bubbles. We also report the first measurement of the ion temperature in the intracluster medium, which we find to be consistent with the electron temperature. In addition, we present a new measurement of the redshift to the brightest cluster galaxy NGC~1275.

Thanks to its high spectral resolution (~5 eV at 6 keV), the Soft X-ray Spectrometer (SXS) on board Hitomi enables us to measure the detailed structure of spatially resolved emission lines from highly ionized ions in galaxy clusters for the first time. In this series of papers, using the SXS we have measured the velocities of gas motions, metallicities and the multi-temperature structure of the gas in the core of the Perseus cluster. Here, we show that when inferring physical properties from line emissivities in systems like Perseus, the resonant scattering (RS) effect should be taken into account. In the Hitomi waveband, RS mostly affects the FeXXV He$\alpha$ line ($w$) - the strongest line in the spectrum. The flux measured by Hitomi in this line is suppressed by a factor ~1.3 in the inner ~30 kpc, compared to predictions for an optically thin plasma; the suppression decreases with the distance from the center. The $w$ line also appears slightly broader than other lines from the same ion. The observed distortions of the $w$ line flux, shape and distance dependence are all consistent with the expected effect of the resonant scattering in the Perseus core. By measuring the ratio of fluxes in optically thick ($w$) and thin (FeXXV forbidden, He$\beta$, Ly$\alpha$) lines, and comparing these ratios with predictions from Monte Carlo radiative transfer simulations, the velocities of gas motions have been obtained. The results are consistent with the direct measurements of gas velocities from line broadening described elsewhere in this series, although the systematic and statistical uncertainties remain significant. Further improvements in the predictions of line emissivities in plasma models, and deeper observations with future X-ray missions will enable RS measurements to provide powerful constraints on the amplitude and anisotropy of clusters gas motions.

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

To search for giant X-ray pulses correlated with the giant radio pulses (GRPs) from the Crab pulsar, we performed a simultaneous observation of the Crab pulsar with the X-ray satellite Hitomi in the 2 -- 300 keV band and the Kashima NICT radio observatory in the 1.4 -- 1.7 GHz band with a net exposure of about 2 ks on 25 March 2016, just before the loss of the Hitomi mission.The timing performance of the Hitomi instruments was confirmed to meet the timing requirement and about 1,000 and 100 GRPs were simultaneously observed at the main and inter-pulse phases, respectively, and we found no apparent correlation between the giant radio pulses and the X-ray emission in either the main or inter-pulse phases.All variations are within the 2 sigma fluctuations of the X-ray fluxes at the pulse peaks, and the 3 sigma upper limits of variations of main- or inter- pulse GRPs are 22\% or 80\% of the peak flux in a 0.20 phase width, respectively, in the 2 -- 300 keV band.The values become 25\% or 110\% for main or inter-pulse GRPs, respectively, when the phase width is restricted into the 0.03 phase.Among the upper limits from the Hitomi satellite, those in the 4.5-10 keV and the 70-300 keV are obtained for the first time, and those in other bands are consistent with previous reports.Numerically, the upper limits of main- and inter-pulse GRPs in the 0.20 phase width are about (2.4 and 9.3) $\times 10^{-11}$ erg cm$^{-2}$, respectively. No significant variability in pulse profiles implies that the GRPs originated from a local place within the magnetosphere and the number of photon-emitting particles temporally increases.However, the results do not statistically rule out variations correlated with the GRPs, because the possible X-ray enhancement may appear due to a $>0.02$\% brightening of the pulse-peak flux under such conditions.

We present arc-second-resolution data in the radio, IR, optical and X-ray for 4C+19.44 (=PKS 1354+195), the longest and straightest quasar jet with deep X-ray observations. We report results from radio images with half to one arc-second angular resolution at three frequencies, plus HST and Spitzer data. The Chandra data allow us to measure the X-ray spectral index in 10 distinct regions along the 18 arcsec jet and compare with the radio index. The radio and X-ray spectral indices of the jet regions are consistent with a value of $\alpha =0.80$ throughout the jet, to within 2 sigma uncertainties. The X-ray jet structure to the south extends beyond the prominent radio jet and connects to the southern radio lobe, and there is extended X-ray emission in the direction of the unseen counter jet and coincident with the northern radio lobe. This jet is remarkable since its straight appearance over a large distance allows the geometry factors to be taken as fixed along the jet. Using the model of inverse Compton scattering of the cosmic microwave background (iC/CMB) by relativistic electrons, we find that the magnetic field strengths and Doppler factors are relatively constant along the jet. If instead the X-rays are synchrotron emission, they must arise from a population of electrons distinct from the particles producing the radio synchrotron spectrum.

The Crab nebula originated from a core-collapse supernova (SN) explosion observed in 1054 A.D. When viewed as a supernova remnant (SNR), it has an anomalously low observed ejecta mass and kinetic energy for an Fe-core collapse SN. Intensive searches were made for a massive shell that solves this discrepancy, but none has been detected. An alternative idea is that the SN1054 is an electron-capture (EC) explosion with a lower explosion energy by an order of magnitude than Fe-core collapse SNe. In the X-rays, imaging searches were performed for the plasma emission from the shell in the Crab outskirts to set a stringent upper limit to the X-ray emitting mass. However, the extreme brightness of the source hampers access to its vicinity. We thus employed spectroscopic technique using the X-ray micro-calorimeter onboard the Hitomi satellite. By exploiting its superb energy resolution, we set an upper limit for emission or absorption features from yet undetected thermal plasma in the 2-12 keV range. We also re-evaluated the existing Chandra and XMM-Newton data. By assembling these results, a new upper limit was obtained for the X-ray plasma mass of <~ 1Mo for a wide range of assumed shell radius, size, and plasma temperature both in and out of the collisional equilibrium. To compare with the observation, we further performed hydrodynamic simulations of the Crab SNR for two SN models (Fe-core versus EC) under two SN environments (uniform ISM versus progenitor wind). We found that the observed mass limit can be compatible with both SN models if the SN environment has a low density of <~ 0.03 cm-3 (Fe core) or <~ 0.1 cm-3 (EC) for the uniform density, or a progenitor wind density somewhat less than that provided by a mass loss rate of 10-5 Mo yr-1 at 20 km s-1 for the wind environment.

We evaluate the optical/near-infrared (OIR) color variability of 3C 279 in both gamma-ray flaring and non-flaring states over 7-year timescales using the Small and Medium Aperture Research Telescope System (SMARTS) in Cerro Tololo, Chile and gamma-ray fluxes obtained from the Fermi Gamma-ray Space Telescope. This observing strategy differs from previous blazar color variability studies in two key ways: 1) the reported color variability is assessed across optical through near-infrared wavelengths, and 2) the color variability is assessed over timescales significantly longer than an individual flare or ground-based observing season. We highlight 3C 279 because of its complex color variability, which is difficult to reconcile with the simple "redder when brighter" behavior often associated with Flat Spectrum Radio Quasar (FSRQ) color variability. We suggest that the observed OIR color changes depend on a combination of the jet and disk emission. We parameterize this behavior in terms of a single variable, $\zeta^m_n$, representing a smooth transition from disk-dominated, to a mixed contribution, to a jet-dominated system, which provides an explanation of the long-term OIR color variability in the same blazar over time. This suggests a general scheme that could apply to OIR color variability in other blazars.

We investigate the optical and Wide-field Survey Explorer (WISE) colors of "E+A" identified post-starburst galaxies, including a deep analysis on 190 post-starbursts detected in the 2{\mu}m All Sky Survey Extended Source Catalog. The post-starburst galaxies appear in both the optical green valley and the WISE Infrared Transition Zone (IRTZ). Furthermore, we find that post-starbursts occupy a distinct region [3.4]-[4.6] vs. [4.6]-[12] WISE colors, enabling the identification of this class of transitioning galaxies through the use of broad-band photometric criteria alone. We have investigated possible causes for the WISE colors of post-starbursts by constructing a composite spectral energy distribution (SED), finding that mid-infrared (4-12{\mu}m) properties of post-starbursts are consistent with either 11.3{\mu}m polycyclic aromatic hydrocarbon emission, or Thermally Pulsating Asymptotic Giant Branch (TP-AGB) and post-AGB stars. The composite SED of extended post- starburst galaxies with 22{\mu}m emission detected with signal to noise >3 requires a hot dust component to produce their observed rising mid-infrared SED between 12 and 22{\mu}m. The composite SED of WISE 22{\mu}m non-detections (S/N<3), created by stacking 22{\mu}m images, is also flat, requiring a hot dust component. The most likely source of this mid-infrared emission of these E+A galaxies is a buried active galactic nucleus. The inferred upper limit to the Eddington ratios of post-starbursts are 1e-2 to 1e-4, with an average of 1e-3. This suggests that AGNs are not radiatively dominant in these systems. This could mean that including selections able to identify active galactic nuclei as part of a search for transitioning and post-starburst galaxies would create a more complete census of the transition pathways taken as a galaxy quenches its star formation.

High-resolution X-ray spectroscopy with Hitomi was expected to resolve the origin of the faint unidentified E=3.5 keV emission line reported in several low-resolution studies of various massive systems, such as galaxies and clusters, including the Perseus cluster. We have analyzed the Hitomi first-light observation of the Perseus cluster. The emission line expected for Perseus based on the XMM-Newton signal from the large cluster sample under the dark matter decay scenario is too faint to be detectable in the Hitomi data. However, the previously reported 3.5 keV flux from Perseus was anomalously high compared to the sample-based prediction. We find no unidentified line at the reported high flux level. Taking into account the XMM measurement uncertainties for this region, the inconsistency with Hitomi is at a 99% significance for a broad dark-matter line and at 99.7% for a narrow line from the gas. We do not find anomalously high fluxes of the nearby faint K line or the Ar satellite line that were proposed as explanations for the earlier 3.5 keV detections. We do find a hint of a broad excess near the energies of high-n transitions of Sxvi (E=3.44 keV rest-frame) -- a possible signature of charge exchange in the molecular nebula and another proposed explanation for the unidentified line. While its energy is consistent with XMM pn detections, it is unlikely to explain the MOS signal. A confirmation of this interesting feature has to wait for a more sensitive observation with a future calorimeter experiment.

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 study the significance of mergers in the quenching of star formation in galaxies at z~1 by examining their color-mass distributions for different morphology types. We perform two-dimensional light profile fits to GOODS iz images of ~5000 galaxies and X-ray selected active galactic nucleus (AGN) hosts in the CANDELS/GOODS-north and south fields in the redshift range 0.7<z<1.3. Distinguishing between bulge-dominated and disk-dominated morphologies, we find that disks and spheroids have distinct color-mass distributions, in agreement with studies at z~0. The smooth distribution across colors for the disk galaxies corresponds to a slow exhaustion of gas, with no fast quenching event. Meanwhile, blue spheroids most likely come from major mergers of star-forming disk galaxies, and the dearth of spheroids at intermediate green colors is suggestive of rapid quenching. The distribution of moderate luminosity X-ray AGN hosts is even across colors, in contrast, and we find similar numbers and distributions among the two morphology types with no apparent dependence on Eddington ratio. The high fraction of bulge-dominated galaxies that host an AGN in the blue cloud and green valley is consistent with the scenario in which the AGN is triggered after a major merger, and the host galaxy then quickly evolves into the green valley. This suggests AGN feedback may play a role in the quenching of star formation in the minority of galaxies that undergo major mergers.

We consider energy budgets and radiative history of 8 fading AGN, identified from mismatch between the ionizion of very extended (>10 kpc) ionized clouds and the luminosity of the nucleus viewed directly. All show significant fading on ~50,000-year timescales. We explore the use of minimum ionizing luminosity Q derived from photoionization balance in the brightest pixels in H-alpha at each projected radius. Tests using PG QSOs, and one target with detailed photoionization modeling, suggest that we can derive useful histories of individual AGN; the minimum ionizing luminosity is always an underestimate and subject to fine structure in the ionized material. These tests suggest that the underestimation from the upper envelope of Q values is roughly constant for a given object. These AGN show rapid drops and standstills; the common feature is a rapid drop in the last 20,000 years before our view of the nucleus. E-folding timescales are mostly thousands of years, with a few episodes as short as 400. In the limit of largely obscured AGN, we find additional evidence for fading, comparing lower limits from recombination balance and the maximum luminosities derived from from infrared fluxes. We compare these long-term light curves to simulations of AGN accretion; the strongest variations on these timespans are seen in models with strong and local feedback. Gemini integral-field optical spectroscopy shows a very limited role for outflows in these structures. While rings and loops of emission are common, their kinematic structure shows some to be in regular rotation. UGC 7342 exhibits local signatures of outflows <300 km/s, largely associated with very diffuse emission. Only in the Teacup AGN do we see outflow signatures of order 1000 km/s. Clouds around these fading AGN consist largely of tidal debris being externally illuminated but not displaced by AGN outflows. (Abridged)

We present the Chandra imaging and spectral analysis of NGC 4968, a nearby (z = 0.00986) Seyfert 2 galaxy. We discover extended ($\sim$1 kpc) X-ray emission in the soft band (0.5 - 2 keV) that is neither coincident with the narrow line region nor the extended radio emission. Based on spectral modeling, it is linked to on-going star formation ($\sim$2.6-4 M$_{\sun}$ yr$^{-1}$). The soft emission at circumnuclear scales (inner $\sim$400 pc) originates from hot gas, with kT $\sim$ 0.7 keV, while the most extended thermal emission is cooler (kT $\sim$ 0.3 keV). We refine previous measurements of the extreme Fe K$\alpha$ equivalent width in this source (EW = 2.5$^{+2.6}_{-1.0}$ keV), which suggests the central engine is completely embedded within Compton-thick levels of obscuration. Using physically motivated models fit to the Chandra spectrum, we derive a Compton-thick column density ($N_{\rm H} > 1.25\times10^{24}$ cm$^{-2}$) and an intrinsic hard (2-10 keV) X-ray luminosity of $\sim$3-8$\times 10^{42}$ erg s$^{-1}$ (depending on the presumed geometry of the obscurer), which is over two orders of magnitude larger than that observed. The large Fe K$\alpha$ EW suggests a spherical covering geometry, which could be confirmed with X-ray measurements above 10 keV. NGC 4968 is similar to other active galaxies that exhibit extreme Fe K$\alpha$ EWs (i.e., $>$2 keV) in that they also contain ongoing star formation. This work supports the idea that gas associated with nuclear star formation may increase the covering factor of the enshrouding gas and play a role in obscuring AGN.

We report the first hard X-ray observations with NuSTAR of the BL Lac type blazar PKS 2155-304, augmented with soft X-ray data from XMM-Newton and gamma-ray data from the Fermi Large Area Telescope, obtained in April 2013 when the source was in a very low flux state. A joint NuSTAR and XMM spectrum, covering the energy range 0.5 - 60 keV, is best described by a model consisting of a log-parabola component with curvature beta = 0.3(+0.2,-0.1) and a (local) photon index 3.04 +/- 0.15 at photon energy of 2 keV, and a hard power-law tail with photon index 2.2 +/- 0.4. The hard X-ray tail can be smoothly joined to the quasi-simultaneous gamma-ray spectrum by a synchrotron self-Compton component produced by an electron distribution with index p = 2.2. Assuming that the power-law electron distribution extends down to the minimum electron Lorentz factor gamma_min = 1 and that there is one proton per electron, an unrealistically high total jet power L_p of roughly 10^47 erg/s is inferred. This can be reduced by two orders of magnitude either by considering a significant presence of electron-positron pairs with lepton-to-proton ratio of at least 30, or by introducing an additional, low-energy break in the electron energy distribution at the electron Lorentz factor gamma_br1 of roughly 100. In either case, the jet composition is expected to be strongly matter-dominated.

We present a population study of the star formation history of 1244 Type 2 AGN host galaxies, compared to 6107 inactive galaxies. A Bayesian method is used to determine individual galaxy star formation histories, which are then collated to visualise the distribution for quenching and quenched galaxies within each population. We find evidence for some of the Type 2 AGN host galaxies having undergone a rapid drop in their star formation rate within the last 2 Gyr. AGN feedback is therefore important at least for this population of galaxies. This result is not seen for the quenching and quenched inactive galaxies whose star formation histories are dominated by the effects of downsizing at earlier epochs, a secondary effect for the AGN host galaxies. We show that histories of rapid quenching cannot account fully for the quenching of all the star formation in a galaxy's lifetime across the population of quenched AGN host galaxies, and that histories of slower quenching, attributed to secular (non-violent) evolution, are also key in their evolution. This is in agreement with recent results showing both merger-driven and non-merger processes are contributing to the co-evolution of galaxies and supermassive black holes. The availability of gas in the reservoirs of a galaxy, and its ability to be replenished, appear to be the key drivers behind this co-evolution.

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 CO(1-0) observations of objects within the Shocked POststarburst Galaxy Survey taken with the Institut de Radioastronomie Millimetrique (IRAM) 30m single dish and the Combined Array for Research for Millimeter Astronomy (CARMA) interferometer. Shocked Poststarburst Galaxies (SPOGs) represent a transitioning population of galaxies, with deep Balmer absorption (Hdelta>5A), consistent with an intermediate-age (A-star) stellar population, and ionized gas line ratios inconsistent with pure star formation. The CO(1-0) subsample was selected from SPOGs detected by the Wide-field Infrared Survey Explorer with 22um flux detected at a signal-to-noise (S/N)>3. Of the 52 objects observed in CO(1-0), 47 are detected with S/N>3. A large fraction (37-46%) of our CO-SPOG sample were visually classified as morphologically disrupted. The H2 masses detected were between 10^(8.7-10.8) Msuns, consistent with the gas masses found in normal galaxies, though approximately an order of magnitude larger than the range seen in poststarburst galaxies. When comparing the 22um and CO(1-0) fluxes, SPOGs diverge from the normal star-forming relation, having 22um fluxes in excess by a factor of <EMIR,SPOG>=4.91+0.42-0.39. The Na I D characteristics of CO-SPOGs show that it is likely that many of these objects host interstellar winds. Objects with the large Na I D enhancements also tend to emit in the radio, suggesting possible AGN-driving of neutral winds.

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.

Supermassive black holes in galaxy centres can grow by the accretion of gas, liberating energy that might regulate star formation on galaxy-wide scales. The nature of the gaseous fuel reservoirs that power black hole growth is nevertheless largely unconstrained by observations, and is instead routinely simplified as a smooth, spherical inflow of very hot gas. Recent theory and simulations instead predict that accretion can be dominated by a stochastic, clumpy distribution of very cold molecular clouds - a departure from the "hot mode" accretion model - although unambiguous observational support for this prediction remains elusive. Here we report observations that reveal a cold, clumpy accretion flow towards a supermassive black hole fuel reservoir in the nucleus of the Abell 2597 Brightest Cluster Galaxy (BCG), a nearby (redshift z=0.0821) giant elliptical galaxy surrounded by a dense halo of hot plasma. Under the right conditions, thermal instabilities can precipitate from this hot gas, producing a rain of cold clouds that fall toward the galaxy's centre, sustaining star formation amid a kiloparsec-scale molecular nebula that inhabits its core. The observations show that these cold clouds also fuel black hole accretion, revealing "shadows" cast by the molecular clouds as they move inward at about 300 kilometres per second towards the active supermassive black hole in the galaxy centre, which serves as a bright backlight. Corroborating evidence from prior observations of warmer atomic gas at extremely high spatial resolution, along with simple arguments based on geometry and probability, indicate that these clouds are within the innermost hundred parsecs of the black hole, and falling closer towards it.

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 a new metric that uses the spectral curvature (SC) above 10 keV to identify Compton-thick AGN in low-quality Swift BAT X-ray data. Using NuSTAR, we observe nine high SC-selected AGN. We find that high-sensitivity spectra show the majority are Compton-thick (78% or 7/9) and the remaining two are nearly Compton-thick (NH~5-8x10^23 cm^-2). We find the SC_bat and SC_nustar measurements are consistent, suggesting this technique can be applied to future telescopes. We tested the SC method on well-known Compton-thick AGN and find it is much more effective than broad band ratios (e.g. 100% using SC vs. 20% using 8-24/3-8 keV). Our results suggest that using the >10 keV emission may be the only way to identify this population since only two sources show Compton-thick levels of excess in the OIII to X-ray emission ratio (F_OIII/F_2-10 keV>1) and WISE colors do not identify most of them as AGN. Based on this small sample, we find that a higher fraction of these AGN are in the final merger stage than typical BAT AGN. Additionally, these nine obscured AGN have, on average, four times higher accretion rates than other BAT-detected AGN (Edd ratio=0.068 compared to 0.016). The robustness of SC at identifying Compton-thick AGN implies a higher fraction of nearby AGN may be Compton-thick (~22%) and the sum of black hole growth in Compton-thick AGN (Eddington ratio times population percentage), is nearly as large as mildly obscured and unobscured AGN.

We combined multi-wavelength data for blazars from the Roma-BZCAT catalog and analyzed hundreds of X-ray spectra. We present the fluxes and Spectral Energy Distributions (SEDs), in 12 frequency bands from radio to $\gamma$-rays, for a final sample of 2214 blazars. Using a model-independent statistical approach, we looked for systematic trends in the SEDs; the most significant trends involved the radio luminosities and X-ray spectral indices of the blazars. We used a Principal Component Analysis (PCA), to determine the basis vectors of the blazar SEDs and, in order to maximize the size of the sample, imputed missing fluxes using the K-nearest neighbors method. Using more than an order of magnitude more data than was available when Fossati et al. (1997, 1998) first reported trends of SED shape with blazar luminosity, we confirmed the anti-correlation between radio luminosity and synchrotron peak frequency, although with greater scatter than was seen in the smaller sample. The same trend can be seen between bolometric luminosity and synchrotron peak frequency. Finally, we used all available blazar data to determine an empirical SED description that depends only on the radio luminosity at 1.4~GHz and the redshift. We verified that this statistically significant relation was not a result of the luminosity-luminosity correlations that are natural in flux-limited samples (i.e., where the correlation is actually caused by the redshift rather than the luminosity).

We describe the first data release from the Spitzer-IRAC Equatorial Survey (SpIES); a large-area survey of 115 deg^2 in the Equatorial SDSS Stripe 82 field using Spitzer during its 'warm' mission phase. SpIES was designed to probe sufficient volume to perform measurements of quasar clustering and the luminosity function at z > 3 to test various models for "feedback" from active galactic nuclei (AGN). Additionally, the wide range of available multi-wavelength, multi-epoch ancillary data enables SpIES to identify both high-redshift (z > 5) quasars as well as obscured quasars missed by optical surveys. SpIES achieves 5{\sigma} depths of 6.13 {\mu}Jy (21.93 AB magnitude) and 5.75 {\mu}Jy (22.0 AB magnitude) at 3.6 and 4.5 microns, respectively - depths significantly fainter than WISE. We show that the SpIES survey recovers a much larger fraction of spectroscopically-confirmed quasars (98%) in Stripe 82 than are recovered by WISE (55%). This depth is especially powerful at high-redshift (z > 3.5), where SpIES recovers 94% of confirmed quasars, whereas WISE only recovers 25%. Here we define the SpIES survey parameters and describe the image processing, source extraction, and catalog production methods used to analyze the SpIES data. In addition to this survey paper, we release 234 images created by the SpIES team and three detection catalogs: a 3.6 {\mu}m-only detection catalog containing 6.1 million sources, a 4.5 {\mu}m-only detection catalog containing 6.5 million sources, and a dual-band detection catalog containing 5.4 million sources.

Some reddened quasars appear to be transitional objects in the merger-induced black hole growth/galaxy evolution paradigm, where a heavily obscured nucleus starts to be unveiled by powerful quasar winds evacuating the surrounding cocoon of dust and gas. Hard X-ray observations are able to peer through this gas and dust, revealing the properties of circumnuclear obscuration. Here, we present NuSTAR and XMM-Newton/Chandra observations of FIRST-2MASS selected red quasars F2M 0830+3759 and F2M 1227+3214. We find that though F2M 0830+3759 is moderately obscured ($N_{\rm H,Z} = 2.1\pm0.2 \times10^{22}$ cm$^{-2}$) and F2M 1227+3214 is mildly absorbed ($N_{\rm H,Z} = 3.4^{+0.8}_{-0.7}\times10^{21}$ cm$^{-2}$) along the line-of-sight, heavier global obscuration may be present in both sources, with $N_{\rm H,S} = 3.7^{+4.1}_{-2.6} \times 10^{23}$ cm$^{-2}$ and $< 5.5\times10^{23}$ cm$^{-2}$, for F2M 0830+3759 and F2M 1227+3214, respectively. F2M 0830+3759 also has an excess of soft X-ray emission below 1 keV which is well accommodated by a model where 7% of the intrinsic AGN X-ray emission is scattered into the line-of-sight. While F2M 1227+3214 has a dust-to-gas ratio ($E(B-V)$/$N_{\rm H}$) consistent with the Galactic value, the $E(B-V)$/$N_{\rm H}$ value for F2M 0830+3759 is lower than the Galactic standard, consistent with the paradigm that the dust resides on galactic scales while the X-ray reprocessing gas originates within the dust-sublimation zone of the broad-line-region. The X-ray and 6.1$\mu$m luminosities of these red quasars are consistent with the empirical relations derived for high-luminosity, unobscured quasars, extending the parameter space of obscured AGN previously observed by NuSTAR to higher luminosities.

We present deep Chandra X-ray observations of the core of IC 2497, the galaxy associated with Hanny's Voorwerp and hosting a fading AGN. We find extended soft X-ray emission from hot gas around the low intrinsic luminosity (unobscured) AGN ($L_{\rm bol} \sim 10^{42}-10^{44}$ erg s$^{-1}$). The temperature structure in the hot gas suggests the presence of a bubble or cavity around the fading AGN ($\mbox{E$_{\rm bub}$} \sim 10^{54} - 10^{55}$ erg). A possible scenario is that this bubble is inflated by the fading AGN, which after changing accretion state is now in a kinetic mode. Other possibilities are that the bubble has been inflated by the past luminous quasar ($L_{\rm bol} \sim 10^{46}$ erg s$^{-1}$), or that the temperature gradient is an indication of a shock front from a superwind driven by the AGN. We discuss the possible scenarios and the implications for the AGN-host galaxy interaction, as well as an analogy between AGN and X-ray binaries lifecycles. We conclude that the AGN could inject mechanical energy into the host galaxy at the end of its lifecycle, and thus provide a source for mechanical feedback, in a similar way as observed for X-ray binaries.

Capitalizing on the all-sky coverage of {\it WISE}, and the 35\% and 50\% sky coverage from SDSS and Pan-STARRS, respectively, we explore the efficacy of $m_{R}$ (optical) - $m_{3.4 \mu m}$ (mid-infrared), hereafter $R-W1$, as a color diagnostic to identify obscured supermassive black hole accretion in wide-area X-ray surveys. We use the $\sim$16.5 deg$^2$ Stripe 82 X-ray survey data as a test-bed to compare $R-W1$ with $R-K$, an oft-used obscured AGN selection criterion, and examine where different classes of objects lie in this parameter space. Most stars follow a well-defined path in $R-K$ vs. $R-W1$ space. We demonstrate that optically normal galaxies hosting X-ray AGN at redshifts $0.5<z<1$ can be recovered with an $R-W1>4$ color-cut, while they typically are not selected as AGN based on their $W1-W2$ colors. Additionally, different observed X-ray luminosity bins favor different regions in $R-W1$ parameter space: moderate luminosity AGN ($10^{43}$ erg s$^{-1} < L_{\rm 0.5-10 keV} < 10^{44}$ erg s$^{-1}$) tend to have red colors while the highest luminosity AGN ($L_{\rm 0.5-10 keV} > 10^{45}$ erg s$^{-1}$) have bluer colors; higher spectroscopic completeness of the Stripe 82X sample is needed to determine whether this is a selection effect or an intrinsic property. Finally, we parameterize X-ray obscuration of Stripe 82X AGN by calculating their hardness ratios (HRs) and find no clear trends between HR and optical reddening. Our results will help inform best-effort practices in following-up obscured AGN candidates in current and future wide-area, shallow X-ray surveys, including the all-sky {\it eROSITA} mission.

We present coordinated multiwavelength observations of the bright, nearby BL Lac object Mrk 421 taken in 2013 January-March, involving GASP-WEBT, Swift, NuSTAR, Fermi-LAT, MAGIC, VERITAS, and other collaborations and instruments, providing data from radio to very-high-energy (VHE) gamma-ray bands. NuSTAR yielded previously unattainable sensitivity in the 3-79 keV range, revealing that the spectrum softens when the source is dimmer until the X-ray spectral shape saturates into a steep power law with a photon index of approximately 3, with no evidence for an exponential cutoff or additional hard components up to about 80 keV. For the first time, we observed both the synchrotron and the inverse-Compton peaks of the spectral energy distribution (SED) simultaneously shifted to frequencies below the typical quiescent state by an order of magnitude. The fractional variability as a function of photon energy shows a double-bump structure which relates to the two bumps of the broadband SED. In each bump, the variability increases with energy which, in the framework of the synchrotron self-Compton model, implies that the electrons with higher energies are more variable. The measured multi-band variability, the significant X-ray-to-VHE correlation down to some of the lowest fluxes ever observed in both bands, the lack of correlation between optical/UV and X-ray flux, the low degree of polarization and its significant (random) variations, the short estimated electron cooling time, and the significantly longer variability timescale observed in the NuSTAR light curves point toward in-situ electron acceleration, and suggest that there are multiple compact regions contributing to the broadband emission of Mrk 421 during low-activity states.

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.

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

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

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

Supermassive black holes (SMBHs) and their host galaxies are generally thought to coevolve, so that the SMBH achieves up to about 0.2 to 0.5% of the host galaxy mass in the present day. The radiation emitted from the growing SMBH is expected to affect star formation throughout the host galaxy. The relevance of this scenario at early cosmic epochs is not yet established. We present spectroscopic observations of a galaxy at redshift z = 3.328, which hosts an actively accreting, extremely massive BH, in its final stages of growth. The SMBH mass is roughly one-tenth the mass of the entire host galaxy, suggesting that it has grown much more efficiently than the host, contrary to models of synchronized coevolution. The host galaxy is forming stars at an intense rate, despite the presence of a SMBH-driven gas outflow.

New images from the Hubble Space Telescope of the FRII radio galaxy Pictor A reveal a previously undiscovered tidal tail, as well as a number of jet knots coinciding with a known X-ray and radio jet. The tidal tail is approximately 5" wide (3 kpc projected), starting 18" (12 kpc) from the center of Pictor A, and extends more than 90" (60 kpc). The knots are part of a jet observed to be about 4' (160 kpc) long, extending to a bright hotspot. These images are the first optical detections of this jet, and by extracting knot flux densities through three filters we set constraints on emission models. While the radio and optical flux densities are usually explained by synchrotron emission, there are several emission mechanisms which might be used to explain the X-ray flux densities. Our data rule out Doppler boosted inverse Compton scattering as a source of the high energy emission. Instead, we find that the observed emission can be well described by synchrotron emission from electrons with a low energy index ($p\sim2$) that dominates the radio band, while a high energy index ($p\sim3$) is needed for the X-ray band and the transition occurs in the optical/infrared band. This model is consistent with a continuous electron injection scenario.

We used the Hubble Space Telescope WFC3 near-infrared camera to image the host galaxies of a sample of eleven luminous, dust-reddened quasars at z ~ 2 -- the peak epoch of black hole growth and star formation in the Universe -- to test the merger-driven picture for the co-evolution of galaxies and their nuclear black holes. The red quasars come from the FIRST+2MASS red quasar survey and a newer, deeper, UKIDSS+FIRST sample. These dust-reddened quasars are the most intrinsically luminous quasars in the Universe at all redshifts, and may represent the dust-clearing transitional phase in the merger-driven black hole growth scenario. Probing the host galaxies in rest-frame visible light, the HST images reveal that 8/10 of these quasars have actively merging hosts, while one source is reddened by an intervening lower redshift galaxy along the line-of-sight. We study the morphological properties of the quasar hosts using parametric Sersic fits as well as the non-parametric estimators (Gini coefficient, M_{20} and asymmetry). Their properties are heterogeneous but broadly consistent with the most extreme morphologies of local merging systems such as Ultraluminous Infrared galaxies. The red quasars have a luminosity range of log(L_bol) = 47.8 - 48.3 (erg/s) and the merger fraction of their AGN hosts is consistent with merger-driven models of luminous AGN activity at z=2, which supports the picture in which luminous quasars and galaxies co-evolve through major mergers that trigger both star formation and black hole growth.

We report the results of a multi-band observing campaign on the famous blazar 3C 279 conducted during a phase of increased activity from 2013 December to 2014 April, including first observations of it with NuSTAR. The $\gamma$-ray emission of the source measured by Fermi-LAT showed multiple distinct flares reaching the highest flux level measured in this object since the beginning of the Fermi mission, with $F(E > 100\,{\rm MeV})$ of $10^{-5}$ photons cm$^{-2}$ s$^{-1}$, and with a flux doubling time scale as short as 2 hours. The $\gamma$-ray spectrum during one of the flares was very hard, with an index of $\Gamma_\gamma = 1.7 \pm 0.1$, which is rarely seen in flat spectrum radio quasars. The lack of concurrent optical variability implies a very high Compton dominance parameter $L_\gamma/L_{\rm syn} > 300$. Two 1-day NuSTAR observations with accompanying Swift pointings were separated by 2 weeks, probing different levels of source activity. While the 0.5$-$70 keV X-ray spectrum obtained during the first pointing, and fitted jointly with Swift-XRT is well-described by a simple power law, the second joint observation showed an unusual spectral structure: the spectrum softens by $\Delta\Gamma_{\rm X} \simeq 0.4$ at $\sim$4 keV. Modeling the broad-band SED during this flare with the standard synchrotron plus inverse Compton model requires: (1) the location of the $\gamma$-ray emitting region is comparable with the broad line region radius, (2) a very hard electron energy distribution index $p \simeq 1$, (3) total jet power significantly exceeding the accretion disk luminosity $L_{\rm j}/L_{\rm d} \gtrsim 10$, and (4) extremely low jet magnetization with $L_{\rm B}/L_{\rm j} \lesssim 10^{-4}$. We also find that single-zone models that match the observed $\gamma$-ray and optical spectra cannot satisfactorily explain the production of X-ray emission.

We present narrow- and medium-band HST imaging, with additional supporting ground-based data, for 8 galaxies identified as hosting fading AGN. These have AGN-ionized gas projected >10 kpc from the nucleus, and significant shortfall of ionizing radiation between the distant gas and the AGN, indicating fading AGN on ~50,000-year timescales. Every system shows evidence of ongoing or past interactions; a similar sample of obscured AGN with extended ionized clouds shares this incidence of disturbances. Several systems show multiple dust lanes in different orientations, broadly fit by differentially precessing disks of accreted material ~1.5 Gyr after initial arrival. The gas has lower metallicity than the nuclei; three systems have abundances uniformly well below solar, consistent with an origin in tidally disrupted low-luminosity galaxies, while some systems have more nearly solar abundances (accompanied by such signatures as multiple Doppler components), which may suggest redistribution of gas by outflows within the host galaxies themselves. These aspects are consistent with a tidal origin for the extended gas in most systems, although the ionized gas and stellar tidal features do not always match closely. In contrast to clouds near radio-loud AGN, these are dominated by rotation, in some cases in warped disks. Outflows are important only in localized regions near some of the AGN. In UGC 7342 and UGC 11185, luminous star clusters are seen within projected ionization cones, potentially triggered by outflows. As in the discovery example Hanny's Voorwerp/IC 2497, some clouds lack a strong correlation between H-alpha surface brightness and ionization parameter, indicating unresolved fine structure. Together with thin coherent filaments spanning several kpc, persistence of these structures over their orbital lifetimes may require a role for magnetic confinement. (Abridged)

We present multi-epoch optical spectroscopy of seven southern Fermi-monitored blazars from 2008 - 2013 using the Small and Medium Aperture Research Telescope System (SMARTS), with supplemental spectroscopy and polarization data from the Steward Observatory. We find that the emission lines are much less variable than the continuum; 4 of 7 blazars had no detectable emission line variability over the 5 years. This is consistent with photoionization primarily by an accretion disk, allowing us to use the lines as a probe of disk activity. Comparing optical emission line flux with Fermi $\gamma$-ray flux and optical polarized flux, we investigate whether relativistic jet variability is related to the accretion flow. In general, we see no such dependence, suggesting the jet variability is likely caused by internal processes like turbulence or shock acceleration rather than a variable accretion rate. However, three sources showed statistically significant emission line flares in close temporal proximity to very large Fermi $\gamma$-ray flares. While we do not have sufficient emission line data to quantitatively assess their correlation with the $\gamma$-ray flux, it appears that in some cases, the jet might provide additional photoionizing flux to the broad line region, which implies some gamma-rays are produced within the broad line region, at least for these large flares.

Recent work suggests blue ellipticals form in mergers and migrate quickly from the blue cloud of star-forming galaxies to the red sequence of passively evolving galaxies, perhaps as a result of black hole feedback. Such rapid reddening of stellar populations implies that large gas reservoirs in the pre-merger star-forming pair must be depleted on short time scales. Here we present pilot observations of atomic hydrogen gas in four blue early-type galaxies that reveal increasing spatial offsets between the gas reservoirs and the stellar components of the galaxies, with advancing post-starburst age. Emission line spectra show associated nuclear activity in two of the merged galaxies, and in one case radio lobes aligned with the displaced gas reservoir. These early results suggest that a kinetic process (possibly feedback from black hole activity) is driving the quick truncation of star formation in these systems, rather than a simple exhaustion of gas supply.

We investigate early black hole (BH) growth through the methodical search for $z\gtrsim5$ AGN in the $Chandra$ Deep Field South. We base our search on the $Chandra$ 4-Ms data with flux limits of $9.1\times\ 10^{-18}$ (soft, 0.5 - 2 keV) and $5.5\times\ 10^{-17}\ \mathrm{erg}\ \mathrm{s}^{-1}\ \mathrm{cm}^{-2}$ (hard, 2 - 8 keV). At $z\sim5$ this corresponds to luminosities as low as $\sim10^{42}$ ($\sim10^{43}$) $\mathrm{erg}\ \mathrm{s}^{-1}$ in the soft (hard) band and should allow us to detect Compton-thin AGN with $M_\mathrm{BH}>10^7 M_{\odot}$ and Eddington ratios > 0.1. Our field ($0.03~\mathrm{deg}^2$) contains over 600 $z\sim5$ Lyman Break Galaxies. Based on lower redshift relations we would expect $\sim20$ of them to host AGN. After combining the $Chandra$ data with GOODS/ACS, CANDELS/WFC3 and $Spitzer$/IRAC data, the sample consists of 58 high-redshift candidates. We run a photometric redshift code, stack the GOODS/ACS data, apply colour criteria and the Lyman Break Technique and use the X-ray Hardness Ratio. We combine our tests and using additional data find that all sources are most likely at low redshift. We also find five X-ray sources without a counterpart in the optical or infrared which might be spurious detections. We conclude that our field does not contain any convincing $z\gtrsim5$ AGN. Explanations for this result include a low BH occupation fraction, a low AGN fraction, short, super-Eddington growth modes, BH growth through BH-BH mergers or in optically faint galaxies. By searching for $z\gtrsim5$ AGN we are setting the foundation for constraining early BH growth and seed formation scenarios.

SDSS J015957.64+003310.5 is an X-ray selected, $z=0.31$ AGN from the Stripe 82X survey that transitioned from a Type 1 quasar to a Type 1.9 AGN between 2000 and 2010. This is the most distant AGN, and first quasar, yet observed to have undergone such a dramatic change. We re-observed the source with the double spectrograph on the Palomar 5m telescope in July 2014 and found that the spectrum is unchanged since 2010. From fitting the optical spectra, we find that the AGN flux dropped by a factor of 6 between 2000 and 2010 while the broad H$\alpha$ emission faded and broadened. Serendipitous X-ray observations caught the source in both the bright and dim state, showing a similar 2-10 keV flux diminution as the optical while lacking signatures of obscuration. The optical and X-ray changes coincide with $g$-band magnitude variations over multiple epochs of Stripe 82 observations. We demonstrate that variable absorption, as might be expected from the simplest AGN unification paradigm, does not explain the observed photometric or spectral properties. We interpret the changing state of J0159+0033 to be caused by dimming of the AGN continuum, reducing the supply of ionizing photons available to excite gas in the immediate vicinity around the black hole. J0159+0033 provides insight into the intermittency of black hole growth in quasars, as well as an unprecedented opportunity to study quasar physics (in the bright state) and the host galaxy (in the dim state), which has been impossible to do in a single sources until now.

We report on observations of NGC1068 with NuSTAR, which provide the best constraints to date on its $>10$~keV spectral shape. We find no strong variability over the past two decades, consistent with its Compton-thick AGN classification. The combined NuSTAR, Chandra, XMM-Newton, and Swift-BAT spectral dataset offers new insights into the complex reflected emission. The critical combination of the high signal-to-noise NuSTAR data and a spatial decomposition with Chandra allow us to break several model degeneracies and greatly aid physical interpretation. When modeled as a monolithic (i.e., a single N_H) reflector, none of the common Compton-reflection models are able to match the neutral fluorescence lines and broad spectral shape of the Compton reflection. A multi-component reflector with three distinct column densities (e.g., N_H~1.5e23, 5e24, and 1e25 cm^{-2}) provides a more reasonable fit to the spectral lines and Compton hump, with near-solar Fe abundances. In this model, the higher N_H components provide the bulk of the Compton hump flux while the lower N_H component produces much of the line emission, effectively decoupling two key features of Compton reflection. We note that ~30% of the neutral Fe Kalpha line flux arises from >2" (~140 pc), implying that a significant fraction of the <10 keV reflected component arises from regions well outside of a parsec-scale torus. These results likely have ramifications for the interpretation of poorer signal-to-noise observations and/or more distant objects [Abridged].

Limited studies have been performed on the radio-loud fraction in X-ray selected type 1 AGN samples. The consistency between various radio-loudness definitions also needs to be checked. We measure the radio-loudness of the 407 type 1 AGNs in the XMM-COSMOS quasar sample using nine criteria from the literature (six defined in the rest-frame and three defined in the observed frame): $R_L=\log(L_{5GHz}/L_B)$, $q_{24}=\log(L_{24\mu m}/L_{1.4GHz})$, $R_{uv}=\log(L_{5GHz}/L_{2500\AA})$, $R_{i}=\log(L_{1.4GHz}/L_i)$, $R_X=\log(\nu L_{\nu}(5GHz)/L_X)$, $P_{5GHz}=\log(P_{5GHz}(W/Hz/Sr))$, $R_{L,obs}=\log(f_{1.4GHz}/f_B)$ (observed frame), $R_{i,obs}=\log(f_{1.4GHz}/f_i)$ (observed frame), and $q_{24, obs}=\log(f_{24\mu m}/f_{1.4GHz})$ (observed frame). Using any single criterion defined in the rest-frame, we find a low radio-loud fraction of $\lesssim 5\%$ in the XMM-COSMOS type 1 AGN sample, except for $R_{uv}$. Requiring that any two criteria agree reduces the radio-loud fraction to $\lesssim 2\%$ for about 3/4 of the cases. The low radio-loud fraction cannot be simply explained by the contribution of the host galaxy luminosity and reddening. The $P_{5GHz}=\log(P_{5GHz}(W/Hz/Sr))$ gives the smallest radio-loud fraction. Two of the three radio-loud fractions from the criteria defined in the observed frame without k-correction ($R_{L,obs}$ and $R_{i,obs}$) are much larger than the radio-loud fractions from other criteria.

We present Nustar 3-40 keV observations of the optically selected Type 2 quasar (QSO2) SDSS J1034+6001 or Mrk 34. The high-quality hard X-ray spectrum and archival XMM-Newton data can be fitted self-consistently with a reflection-dominated continuum and strong Fe Kalpha fluorescence line with equivalent-width >1 keV. Prior X-ray spectral fitting below 10 keV showed the source to be consistent with being obscured by Compton-thin column densities of gas along the line-of-sight, despite evidence for much higher columns from multiwavelength data. NuSTAR now enables a direct measurement of this column, and shows that Nh lies in the Compton-thick (CT) regime. The new data also show a high intrinsic 2-10 keV luminosity of L_{2-10}~10^{44} erg/s, in contrast to previous low-energy X-ray measurements for which L_{2-10}<~10^{43} erg/s (i.e. X-ray selection below 10 keV does not pick up this source as an intrinsically luminous obscured quasar). Both the obscuring column and the intrinsic power are about an order of magnitude (or more) larger than inferred from pre-NuSTAR X-ray spectral fitting. Mrk34 is thus a 'gold standard' CT QSO2 and is the nearest non-merging system in this class, in contrast to the other local CT quasar NGC6240 which is currently undergoing a major merger coupled with strong star-formation. For typical X-ray bolometric correction factors, the accretion luminosity of Mrk34 is high enough to potentially power the total infrared luminosity. X-ray spectral fitting also shows that thermal emission related to star-formation is unlikely to drive the observed bright soft component below ~3 keV, favoring photionization instead.

The Circinus galaxy is one of the nearest obscured AGN, making it an ideal target for detailed study. Combining archival Chandra and XMM-Newton data with new NuSTAR observations, we model the 2-79 keV spectrum to constrain the primary AGN continuum and to derive physical parameters for the obscuring material. Chandra's high angular resolution allows a separation of nuclear and off-nuclear galactic emission. In the off-nuclear diffuse emission we find signatures of strong cold reflection, including high equivalent-width neutral Fe lines. This Compton-scattered off-nuclear emission amounts to 18% of the nuclear flux in the Fe line region, but becomes comparable to the nuclear emission above 30 keV. The new analysis no longer supports a prominent transmitted AGN component in the observed band. We find that the nuclear spectrum is consistent with Compton-scattering by an optically-thick torus, where the intrinsic spectrum is a powerlaw of photon index $\Gamma = 2.2-2.4$, the torus has an equatorial column density of $N_{\rm H} = (6-10)\times10^{24}$cm$^{-2}$ and the intrinsic AGN $2-10$ keV luminosity is $(2.3-5.1)\times 10^{42}$ erg/s. These values place Circinus along the same relations as unobscured AGN in accretion rate-vs-$\Gamma$ and $L_X$-vs-$L_{IR}$ phase space. NuSTAR's high sensitivity and low background allow us to study the short time-scale variability of Circinus at X-ray energies above 10 keV for the first time. The lack of detected variability favors a Compton-thick absorber, in line with the the spectral fitting results.

Using self-consistent, physically motivated models, we investigate the X-ray obscuration in 19 Type 2 [OIII] 5007 \AA\ selected AGN, 9 of which are local Seyfert 2 galaxies and 10 of which are Type 2 quasar candidates. We derive reliable line-of-sight and global column densities for these objects, which is the first time this has been reported for an AGN sample; 4 AGN have significantly different global and line-of-sight column densities. Five sources are heavily obscured to Compton-thick. We comment on interesting sources revealed by our spectral modeling, including a candidate ``naked'' Sy2. After correcting for absorption, we find that the ratio of the rest-frame, 2-10 keV luminosity (L$_{\rm 2-10keV,in}$) to L$_{\rm [OIII]}$ is 1.54 $\pm$ 0.49 dex which is essentially identical to the mean Type 1 AGN value. The Fe K$\alpha$ luminosity is significantly correlated with L$_{\rm [OIII]}$, but with substantial scatter. Finally, we do not find a trend between L$_{\rm 2-10keV,in}$ and global or line-of-sight column density, between column density and redshift, between column density and scattering fraction or between scattering fraction and redshift.

We use SDSS+\textit{GALEX}+Galaxy Zoo data to study the quenching of star formation in low-redshift galaxies. We show that the green valley between the blue cloud of star-forming galaxies and the red sequence of quiescent galaxies in the colour-mass diagram is not a single transitional state through which most blue galaxies evolve into red galaxies. Rather, an analysis that takes morphology into account makes clear that only a small population of blue early-type galaxies move rapidly across the green valley after the morphologies are transformed from disk to spheroid and star formation is quenched rapidly. In contrast, the majority of blue star-forming galaxies have significant disks, and they retain their late-type morphologies as their star formation rates decline very slowly. We summarize a range of observations that lead to these conclusions, including UV-optical colours and halo masses, which both show a striking dependence on morphological type. We interpret these results in terms of the evolution of cosmic gas supply and gas reservoirs. We conclude that late-type galaxies are consistent with a scenario where the cosmic supply of gas is shut off, perhaps at a critical halo mass, followed by a slow exhaustion of the remaining gas over several Gyr, driven by secular and/or environmental processes. In contrast, early-type galaxies require a scenario where the gas supply and gas reservoir are destroyed virtually instantaneously, with rapid quenching accompanied by a morphological transformation from disk to spheroid. This gas reservoir destruction could be the consequence of a major merger, which in most cases transforms galaxies from disk to elliptical morphology, and mergers could play a role in inducing black hole accretion and possibly AGN feedback.

We present NuSTAR hard X-ray (3-79 keV) observations of three Type 2 quasars at z ~ 0.4-0.5, optically selected from the Sloan Digital Sky Survey (SDSS). Although the quasars show evidence for being heavily obscured Compton-thick systems on the basis of the 2-10 keV to [OIII] luminosity ratio and multiwavelength diagnostics, their X-ray absorbing column densities (N_H) are poorly known. In this analysis: (1) we study X-ray emission at >10 keV, where X-rays from the central black hole are relatively unabsorbed, in order to better constrain N_H; (2) we further characterize the physical properties of the sources through broad-band near-UV to mid-IR spectral energy distribution (SED) analyses. One of the quasars is detected with NuSTAR at >8 keV with a no-source probability of <0.1%, and its X-ray band ratio suggests near Compton-thick absorption with N_H \gtrsim 5 x 10^23 cm^-2. The other two quasars are undetected, and have low X-ray to mid-IR luminosity ratios in both the low energy (2-10 keV) and high energy (10-40 keV) X-ray regimes that are consistent with extreme, Compton-thick absorption (N_H \gtrsim 10^24 cm^-2). We find that for quasars at z ~ 0.5, NuSTAR provides a significant improvement compared to lower energy (<10 keV) Chandra and XMM-Newton observations alone, as higher column densities can now be directly constrained.

The mid-infrared to ultraviolet (0.1 -- 10 $\mu m$) spectral energy distribution (SED) shapes of 407 X-ray-selected radio-quiet type 1 AGN in the wide-field ``Cosmic Evolution Survey" (COSMOS) have been studied for signs of evolution. For a sub-sample of 200 radio-quiet quasars with black hole mass estimates and host galaxy corrections, we studied their mean SEDs as a function of a broad range of redshift, bolometric luminosity, black hole mass and Eddington ratio, and compared them with the Elvis et al. (1994, E94) type 1 AGN mean SED. We found that the mean SEDs in each bin are closely similar to each other, showing no statistical significant evidence of dependence on any of the analyzed parameters. We also measured the SED dispersion as a function of these four parameters, and found no significant dependencies. The dispersion of the XMM-COSMOS SEDs is generally larger than E94 SED dispersion in the ultraviolet, which might be due to the broader ``window function'' for COSMOS quasars, and their X-ray based selection.

Mkn 421 is a nearby active galactic nucleus dominated at all wavelengths by a very broad non-thermal continuum thought to arise from a relativistic jet seen at a small angle to the line of sight. Its spectral energy distribution peaks in the X-ray and TeV gamma-ray bands, where the energy output is dominated by cooling of high-energy electrons in the jet. In order to study the electron distribution and its evolution, we carried out a dedicated multi-wavelength campaign, including extensive observations by the recently launched highly sensitive hard X-ray telescope NuSTAR, between December 2012 and May 2013. Here we present some initial results based on NuSTAR data from January through March 2013, as well as calibration observations conducted in 2012. Although the observations cover some of the faintest hard X-ray flux states ever observed for Mkn 421, the sensitivity is high enough to resolve intra-day spectral variability. We find that in this low state the dominant flux variations are smooth on timescales of hours, with typical intra-hour variations of less than 5%. We do not find evidence for either a cutoff in the hard X-ray spectrum, or a rise towards a high-energy component, but rather that at low flux the spectrum assumes a power law shape with a photon index of approximately 3. The spectrum is found to harden with increasing brightness.

We present multi-epoch optical observations of the blazar 3C 454.3 (z = 0.859) from 2008 August through 2011 December, using the SMARTS Consortium 1.5m+RCSpectrograph and 1.3m+ANDICAM in Cerro Tololo, Chile. The spectra reveal that the broad optical emission lines Mg II, H-beta and H-gamma are far less variable than the optical or gamma-ray continuum. Although, the gamma-rays varied by a factor of 100 above the EGRET era flux, the lines generally vary by a factor of 2 or less. Smaller variations in the gamma-ray flux did not produce significant variation in any of the observed emission lines. Therefore, to first order, the ionizing flux from the disk changes only slowly during large variations of the jet. However, two exceptions in the response of the broad emission lines are reported during the largest gamma-ray flares in 2009 December and 2010 November, when significant deviations from the mean line flux in H-gamma and Mg II were observed. H-gamma showed a maximum 3-sigma and 4-sigma deviation in each flare, respectively, corresponding to a factor of 1.7 and 2.5 increase in flux. Mg II showed a 2-sigma deviation in both flares; no variation was detected in H-beta during either flare. These significant deviations from the mean line flux also coincide with 7mm core ejections reported previously (Jorstad et al. 2012). The correlation of the increased emission line flux with mm core ejections, and gamma-ray, optical and UV flares suggests that the broad line region extends beyond the gamma-emitting region during the 2009 and 2010 flares.

We have analyzed the {\it XMM-Newton} and {\it Chandra} data overlapping $\sim$16.5 deg$^2$ of Sloan Digital Sky Survey Stripe 82, including $\sim$4.6 deg$^2$ of proprietary {\it XMM-Newton} data that we present here. In total, 3362 unique X-ray sources are detected at high significance. We derive the {\it XMM-Newton} number counts and compare them with our previously reported {\it Chandra} Log$N$-Log$S$ relations and other X-ray surveys. The Stripe 82 X-ray source lists have been matched to multi-wavelength catalogs using a maximum likelihood estimator algorithm. We discovered the highest redshift ($z=5.86$) quasar yet identified in an X-ray survey. We find 2.5 times more high luminosity (L$_x \geq 10^{45}$ erg s$^{-1}$) AGN than the smaller area {\it Chandra} and {\it XMM-Newton} survey of COSMOS and 1.3 times as many identified by XBo\"otes. Comparing the high luminosity AGN we have identified with those predicted by population synthesis models, our results suggest that this AGN population is a more important component of cosmic black hole growth than previously appreciated. Approximately a third of the X-ray sources not detected in the optical are identified in the infrared, making them candidates for the elusive population of obscured high luminosity AGN in the early universe.

B2 1023+25 is an extremely radio-loud quasar at z=5.3 which was first identified as a likely high-redshift blazar candidate in the SDSS+FIRST quasar catalog. Here we use the Nuclear Spectroscopic Telescope Array (NuSTAR) to investigate its non-thermal jet emission, whose high-energy component we detected in the hard X-ray energy band. The X-ray flux is ~5.5x10^(-14) erg cm^(-2)s^(-1) (5-10keV) and the photon spectral index is Gamma_X=1.3-1.6. Modeling the full spectral energy distribution, we find that the jet is oriented close to the line of sight, with a viewing angle of ~3deg, and has significant Doppler boosting, with a large bulk Lorentz factor ~13, which confirms the identification of B2 1023+25 as a blazar. B2 1023+25 is the first object at redshift larger than 5 detected by NuSTAR, demonstrating the ability of NuSTAR to investigate the early X-ray Universe and to study extremely active supermassive black holes located at very high redshift.

We report on the first ten identifications of sources serendipitously detected by the NuSTAR to provide the first sensitive census of the cosmic X-ray background (CXB) source population at >10 keV. We find that these NuSTAR-detected sources are ~100x fainter than those previously detected at >10 keV and have a broad range in redshift and luminosity (z=0.020-2.923 and L_10-40 keV~4x10^{41}-5x10^{45} erg/s); the median redshift and luminosity are z~0.7 and L_10-40 keV~3x10^{44} erg/s, respectively. We characterize these sources on the basis of broad-band ~0.5-32 keV spectroscopy, optical spectroscopy, and broad-band ultraviolet-to-mid-infrared SED analyzes. We find that the dominant source population is quasars with L_10-40 keV>10^{44} erg/s, of which ~50% are obscured with N_H>10^{22} cm^{-2}. However, none of the ten NuSTAR sources are Compton thick (N_H>10^{24} cm^{-2}) and we place a 90% confidence upper limit on the fraction of Compton-thick quasars (L_10-40 keV>10^{44} erg/s) selected at >10 keV of ~33% over the redshift range z=0.5-1.1. We jointly fitted the rest-frame ~10-40 keV data for all of the non-beamed sources with L_10-40 keV>10^{43} erg/s to constrain the average strength of reflection; we find R<1.4 for \Gamma=1.8, broadly consistent with that found for local AGNs observed at >10 keV. We also constrain the host galaxy masses and find a median stellar mass of ~10^{11} M_sun, a factor ~5 times higher than the median stellar mass of nearby high-energy selected AGNs, which may be at least partially driven by the order of magnitude higher X-ray luminosities of the NuSTAR sources. Within the low source-statistic limitations of our study, our results suggest that the overall properties of the NuSTAR sources are broadly similar to those of nearby high-energy selected AGNs but scaled up in luminosity and mass.

We define a quasar-galaxy mixing diagram using the slopes of their spectral energy distributions (SEDs) from 1\mu m to 3000\AA\ and from 1\mu m to 3\mu m in the rest frame. The mixing diagram can easily distinguish among quasar-dominated, galaxy-dominated and reddening-dominated SED shapes. By studying the position of the 413 XMM selected Type 1 AGN in the wide-field "Cosmic Evolution Survey" (COSMOS) in the mixing diagram, we find that a combination of the Elvis et al. (1994, hereafter E94) quasar SED with various contributions from galaxy emission and some dust reddening is remarkably effective in describing the SED shape from 0.3-3\mu m for large ranges of redshift, luminosity, black hole mass and Eddington ratio of type 1 AGN. In particular, the location in the mixing diagram of the highest luminosity AGN is very close (within 1\sigma) to that of the E94 SED. The mixing diagram can also be used to estimate the host galaxy fraction and reddening in quasar. We also show examples of some outliers which might be AGN in different evolutionary stages compared to the majority of AGN in the quasar-host galaxy co-evolution cycle.

We report multi-wavelength (IR-O, UV, X-ray, gamma-ray) data of the flat spectrum radio quasar (FSRQ) PKS B1424-418 (z=1.52) gathered during an active phase in 2013 April . As for a few other cases of FSRQs reported in literature, the detection by the Large Area Telescope onboard Fermi of gamma rays with energy above 10 GeV indicates that the emission likely occurs beyond the highly opaque (tau ~10) broad line region. This conclusion is strengthened by a model fit to the spectral energy distribution, which allows us to locate the emission region even beyond the distance generally assumed for the dusty torus. The consequent large size (~1 pc) inferred for the emission region cannot account for the observed daily-scale variability of the gamma-ray flux. We discuss the possibility that short-term variability results from fast magnetic reconnection events, as proposed in recent works.

We present the redshift lower limit of z>0.6035 for the very-high-energy (VHE; E>100 GeV) emitting blazar PKS 1424+240 (PG 1424+240). This limit is inferred from Lyman beta and gamma absorption observed in the far-ultraviolet spectra from the Hubble Space Telescope/Cosmic Origins Spectrograph. No VHE-detected blazar has shown solid spectroscopic evidence of being more distant. At this distance, VHE observations by VERITAS are shown to sample historically large gamma-ray opacity values at 500 GeV, extending beyond tau=4 for low-level models of the extragalactic background light (EBL) and beyond tau=5 for high-levels. The majority of the z=0.6035 absorption-corrected VHE spectrum appears to exhibit a lower flux than an extrapolation of the contemporaneous LAT power-law fit beyond 100 GeV. However, the highest energy VERITAS point is the only point showing agreement with this extrapolation, possibly implying the overestimation of the gamma-ray opacity or the onset of an unexpected VHE spectral feature. A curved log parabola is favored when fitting the full range of gamma-ray data (0.5 to 500 GeV). While fitting the absorption-corrected VHE data alone results in a harder differential power law than that from the full range, the indices derived using three EBL models are consistent with the physically motivated limit set by Fermi acceleration processes.

We present the first results of a wide area X-ray survey within the Sloan Digital Sky Survey (SDSS) Stripe 82, a 300 deg$^2$ region of the sky with a substantial investment in multi-wavelength coverage. We analyzed archival {\it Chandra} observations that cover 7.5 deg$^2$ within Stripe 82 ("Stripe 82 ACX"), reaching 4.5$\sigma$ flux limits of 7.9$\times10^{-16}$, 3.4$\times10^{-15}$ and 1.8$\times10^{-15}$ erg s$^{-1}$ cm$^{-2}$ in the soft (0.5-2 keV), hard (2-7 keV) and full (0.5-7 keV) bands, to find 774, 239 and 1118 X-ray sources, respectively. Three hundred twenty-one sources are detected only in the full band and 9 sources are detected solely in the soft band. Utilizing data products from the {\it Chandra} Source Catalog, we construct independent Log$N$-Log$S$ relationships, detailing the number density of X-ray sources as a function of flux, which show general agreement with previous {\it Chandra} surveys. We compare the luminosity distribution of Stripe 82 ACX with the smaller, deeper CDF-S + E-CDFS surveys and with {\it Chandra}-COSMOS, illustrating the benefit of wide-area surveys in locating high luminosity AGN. We also investigate the differences and similarities of X-ray and optical selection to uncover obscured AGN in the local Universe. Finally, we estimate the population of AGN we expect to find with increased coverage of 100 deg$^2$ or 300 deg$^2$, which will provide unprecedented insight into the high redshift, high luminosity regime of black hole growth currently under-represented in X-ray surveys.

The Yale/SMARTS optical-near-IR monitoring program has followed the variations in emission of the Fermi-LAT monitored blazars in the southern sky with closely spaced observations since 2008. We report the discovery of an optical-near-IR (OIR) outburst with no accompanying gamma-rays in the blazar PKS 0208-512, one of the targets of this program. While the source undergoes three outbursts of 1 mag or more at OIR wavelengths lasting for longer than 3 months during 2008-2011, only interval 1 and 3 have corresponding bright phases in GeV energies lasting longer than 1 month. The OIR outburst during interval 2 is comparable in brightness and temporal extent to the OIR flares during intervals 1 and 3 which do have gamma-ray counterparts. Gamma-ray and OIR variability are very well-correlated in most cases in the Fermi blazars and the lack of correlation in this case is anomalous. By analyzing the gamma-ray, OIR, and supporting multi-wavelength variability data in details, we speculate that the location of the outburst in the jet during interval 2 was closer to the black hole where the jet is more compact and the magnetic field strength is higher, and the bulk Lorentz factor of the material in the jet is smaller. These result in a much lower Compton dominance and no observable gamma-ray outburst during interval 2.

Using the spectra of a sample of $\sim$28,000 nearby obscured active galaxies from Data Release 7 of the Sloan Digital Sky Survey (SDSS), we probe the connection between AGN activity and star formation over a range of radial scales in the host galaxy. We use the extinction-corrected luminosity of the [OIII] 5007 \AA\ line as a proxy of intrinsic AGN power and supermassive black hole (SMBH) accretion rate. The star formation rates (SFRs) are taken from the MPA-JHU value-added catalog and are measured through the 3$^{\prime\prime}$ SDSS aperture. We construct matched samples of galaxies covering a range in redshifts. With increasing redshift, the projected aperture size encompasses increasing amounts of the host galaxy. This allows us to trace the radial distribution of star-formation as a function of AGN luminosity. We find that the star formation becomes more centrally concentrated with increasing AGN luminosity and Eddington ratio. This implies that such circumnuclear star formation is associated with AGN activity, and that it increasingly dominates over omnipresent disk star formation at higher AGN luminosities, placing critical constraints on theoretical models that link host galaxy star formation and SMBH fueling. We parametrize this relationship and find that the star formation on radial scales $<$1.7 kpc, when including a constant disk component, has a sub-linear dependence on SMBH accretion rate: $SFR \propto \dot{M}^{0.36}$, suggesting that angular momentum transfer through the disk limits accretion efficiency rather than the supply from stellar mass loss.

The Nuclear Spectroscopic Telescope Array (NuSTAR) mission, launched on 13 June 2012, is the first focusing high-energy X-ray telescope in orbit. NuSTAR operates in the band from 3 -- 79 keV, extending the sensitivity of focusing far beyond the ~10 keV high-energy cutoff achieved by all previous X-ray satellites. The inherently low-background associated with concentrating the X-ray light enables NuSTAR to probe the hard X-ray sky with a more than one-hundred-fold improvement in sensitivity over the collimated or coded-mask instruments that have operated in this bandpass. Using its unprecedented combination of sensitivity, spatial and spectral resolution, NuSTAR will pursue five primary scientific objectives, and will also undertake a broad program of targeted observations. The observatory consists of two co-aligned grazing-incidence X-ray telescopes pointed at celestial targets by a three-axis stabilized spacecraft. Deployed into a 600 km, near-circular, 6degree inclination orbit, the Observatory has now completed commissioning, and is performing consistent with pre-launch expectations. NuSTAR is now executing its primary science mission, and with an expected orbit lifetime of ten years, we anticipate proposing a guest investigator program, to begin in Fall 2014.

We report the discovery of an anomalous flare in a bright blazar, namely, PKS 0208-512, one of the targets of the Yale/SMARTS optical-near-infrared (OIR) monitoring program of Fermi blazars. We identify three intervals during which PKS 0208-512 undergoes outbursts at OIR wavelengths lasting for longer than 3 months. Its brightness increases and then decreases again by at least 1 magnitude in these intervals. In contrast, the source undergoes bright phases in GeV energies lasting for longer than 1 month during intervals 1 and 3 only. The OIR outburst during interval 2 is comparable in brightness and temporal extent to the OIR flares during intervals 1 and 3 which do have gamma-ray counterparts. By analyzing the gamma-ray, OIR, and supporting multi-wavelength variability data in details, we speculate that the OIR outburst during interval 2 was caused by a change in the magnetic field without any change in the total number of emitting electrons or Doppler factor of the emitting region. Alternatively, it is possible that the location of the outburst in the jet during interval 2 was closer to the black hole where the jet is more compact and the bulk Lorentz factor of the material in the jet is smaller. We also discuss the complex OIR spectral behavior during these three intervals.

The growth of supermassive black holes appears to be driven by galaxy mergers, violent merger-free processes and/or `secular' processes. In order to quantify the effects of secular evolution on black hole growth, we study a sample of active galactic nuclei (AGN) in galaxies with a calm formation history free of significant mergers, a population that heretofore has been difficult to locate. Here we present an initial sample of 13 AGN in massive (M_* >~ 1e10 M_sun) bulgeless galaxies -- which lack the classical bulges believed inevitably to result from mergers -- selected from the Sloan Digital Sky Survey using visual classifications from Galaxy Zoo. Parametric morphological fitting confirms the host galaxies lack classical bulges; any contributions from pseudobulges are very small (typically < 5%). We compute black hole masses for the two broad-line objects in the sample (4.2e6 and 1.2e7 M_sun) and place lower limits on black hole masses for the remaining sample (typically M_BH >~ 1e6 M_sun), showing that significant black hole growth must be possible in the absence of mergers or violent disk instabilities. The black hole masses are systematically higher than expected from established bulge-black hole relations. However, if the mean Eddington ratio of the systems with measured black hole masses (L/L_Edd = 0.065) is typical, 10 of 13 sources are consistent with the correlation between black hole mass and total stellar mass. That pure disk galaxies and their central black holes may be consistent with a relation derived from elliptical and bulge-dominated galaxies with very different formation histories implies the details of stellar galaxy evolution and dynamics may not be fundamental to the co-evolution of galaxies and black holes.

We compute black hole masses and bolometric luminosities for 57 active galactic nuclei (AGN) in the redshift range 1.25 < z < 2.67, selected from the GOODS-South deep multi-wavelength survey field via their X-ray emission. We determine host galaxy morphological parameters by separating the galaxies from their central point sources in deep HST images, and host stellar masses and colors by multi-wavelength SED fitting. 90% of GOODS AGN at these redshifts have detected rest-frame optical nuclear point sources; bolometric luminosities range from 2e43 - 2e46 erg/s. The black holes are growing at a range of accretion rates, with at least 50% of the sample having L/L_Edd < 0.1. 70% of host galaxies have stellar masses M* > 1e10 M_sun, with a range of colors suggesting a complex star formation history. We find no evolution of AGN bolometric luminosity within the sample, and no correlation between AGN bolometric luminosity and host stellar mass, color or morphology. Fully half the sample of host galaxies is disk-dominated, with another 25% having strong disk components. Fewer than 15% of the systems appear to be at some stage of a major merger. These moderate-luminosity AGN hosts are therefore inconsistent with a dynamical history dominated by mergers strong enough to destroy disks, indicating minor mergers or secular processes dominate the co-evolution of galaxies and their central black holes at z ~ 2.

Using multiwavelength surveys of active galactic nuclei across a wide range of bolometric luminosities (10^{43}<L_{bol}(erg/s<5x10^{46}) and redshifts (0<z<3), we find a strong, redshift-independent correlation between the AGN luminosity and the fraction of host galaxies undergoing a major merger. That is, only the most luminous AGN phases are connected to major mergers, while less luminous AGN appear to be driven by secular processes. Combining this trend with AGN luminosity functions to assess the overall cosmic growth of black holes, we find that ~50% by mass is associated with major mergers, while only 10% of AGN by number, the most luminous, are connected to these violent events. Our results suggest that to reach the highest AGN luminosities -where the most massive black holes accreted the bulk of their mass - a major merger appears to be required. The luminosity dependence of the fraction of AGN triggered by major mergers can successfully explain why the observed scatter in the M-\sigma relation for elliptical galaxies is significantly lower than in spirals. The lack of a significant redshift dependence of the L_{bol}-f_{merger} relation suggests that downsizing, i.e., the general decline in AGN and star formation activity with decreasing redshift, is driven by a decline in the frequency of major mergers combined with a decrease in the availability of gas at lower redshifts.

The "Cosmic Evolution Survey" (COSMOS) enables the study of the Spectral Energy Distributions (SEDs) of Active Galactic Nuclei (AGN) because of the deep coverage and rich sampling of frequencies from X-ray to radio. Here we present a SED catalog of 413 X-ray (\xmm) selected type 1 (emission line FWHM$>2000$ km s$^{-1}$) AGN with Magellan, SDSS or VLT spectrum. The SEDs are corrected for the Galactic extinction, for broad emission line contributions, constrained variability, and for host galaxy contribution. We present the mean SED and the dispersion SEDs after the above corrections in the rest frame 1.4 GHz to 40 keV, and show examples of the variety of SEDs encountered. In the near-infrared to optical (rest frame $\sim 8\mu m$-- 4000\AA), the photometry is complete for the whole sample and the mean SED is derived from detections only. Reddening and host galaxy contamination could account for a large fraction of the observed SED variety. The SEDs are all available on-line.

The blazar AO 0235+164 (z = 0.94) has been one of the most active objects observed by Fermi Large Area Telescope (LAT) since its launch in Summer 2008. In addition to the continuous coverage by Fermi, contemporaneous observations were carried out from the radio to {\gamma} -ray bands between 2008 September and 2009 February. In this paper, we summarize the rich multi-wavelength data collected during the campaign (including F-GAMMA, GASP- WEBT, Kanata, OVRO, RXTE, SMARTS, Swift, and other instruments), examine the cross-correlation between the light curves measured in the different energy bands, and interpret the resulting spectral energy distributions in the context of well-known blazar emission models. We find that the {\gamma} -ray activity is well correlated with a series of near-IR/optical flares, accompanied by an increase in the optical polarization degree. On the other hand, the X-ray light curve shows a distinct 20 day high state of unusually soft spectrum, which does not match the extrapolation of the optical/UV synchrotron spectrum. We tentatively interpret this feature as the bulk Compton emission by cold electrons contained in the jet, which requires an accretion disk corona with an effective covering factor of 19% at a distance of 100 Rg . We model the broadband spectra with a leptonic model with external radiation dominated by the infrared emission from the dusty torus.

We present the results from a Chandra pilot study of 12 massive galaxy mergers selected from Galaxy Zoo. The sample includes major mergers down to a host galaxy mass of 10$^{11}$ $M_\odot$ that already have optical AGN signatures in at least one of the progenitors. We find that the coincidences of optically selected active nuclei with mildly obscured ($N_H \lesssim 1.1 \times 10^{22}$ cm$^{-2}$) X-ray nuclei are relatively common (8/12), but the detections are too faint ($< 40$ counts per nucleus; $f_{2-10 keV} \lesssim 1.2 \times 10^{-13}$ erg s$^{-1}$ cm$^{-2}$) to reliably separate starburst and nuclear activity as the origin of the X-ray emission. Only one merger is found to have confirmed binary X-ray nuclei, though the X-ray emission from its southern nucleus could be due solely to star formation. Thus, the occurrences of binary AGN in these mergers are rare (0-8%), unless most merger-induced active nuclei are very heavily obscured or Compton thick.

We present the time variability properties of a sample of six blazars, AO 0235+164, 3C 273, 3C 279, PKS 1510-089, PKS 2155-304, and 3C 454.3, at optical-IR as well as gamma-ray energies. These observations were carried out as a part of the Yale/SMARTS program during 2008-2010 that has followed the variations in emission of the bright Fermi-LAT-monitored blazars in the southern sky with closely-spaced observations at BVRJK bands. We find the optical/IR time variability properties of these blazars to be remarkably similar to those at the gamma-ray energies. The power spectral density (PSD) functions of the R-band variability of all six blazars are fit well by simple power-law functions with negative slope such that there is higher amplitude variability on longer timescales. No clear break is identified in the PSD of any of the sources. The average slope of the PSD of R-band variability of these blazars is similar to what was found by the Fermi team for the gamma-ray variability of a larger sample of bright blazars. This is consistent with leptonic models where the optical-IR and gamma-ray emission is generated by the same population of electrons through synchrotron and inverse-Compton processes, respectively. The prominent flares present in the optical-IR as well as the gamma-ray light curves of these blazars are predominantly symmetric, i.e., have similar rise and decay timescales, indicating that the long-term variability is dominated by the crossing time of radiation or a disturbance through the emission region rather than by the acceleration or energy-loss timescales of the radiating electrons. In the blazar 3C 454.3, which has the highest-quality light curves, the location of a large gamma-ray outburst during 2009 December is consistent with being in the jet at ~18 pc from the central engine. This poses strong constraints on the models of high energy emission in the jets of blazars.

We present multiwavelength data for twelve blazars observed from 2008-2010 as part of an ongoing optical-infrared photometric monitoring project. Sources were selected to be bright, southern (dec < 20 deg) blazars observed by the Fermi Gamma-Ray Space Telescope, with daily and weekly gamma-ray fluxes made available from the start of the Fermi mission. Light curves are presented for the twelve blazars in BVRJK at near-daily cadence. We find that optical and infrared fluxes are well correlated in all sources. Gamma-ray bright flat spectrum radio quasars (FSRQs) in our sample have optical/infrared emission correlated with gamma-rays consistent with inverse Compton-scattering models for GeV emission. In FSRQs, the variability amplitude decreases towards optical/IR wavelengths, consistent with the presence of a thermal emission component from the accretion disk varying on significantly longer timescales than the jet synchrotron emission. In BL Lac objects, variability is mainly constant across wavelengths, consistent with a weak or radiatively inefficient disk. FSRQs have redder optical-infrared colors when they are brighter, while BL Lac objects show no such trend. Several objects show complicated color-magnitude behavior: AO 0235+164 appears in two different states depending on whether it is gamma-ray bright or not. OJ 287 and 3C 279 show some hysteresis tracks in their color-magnitude diagrams. Individual flares may be achromatic or otherwise depart from the trend, suggesting different jet components becoming important at different times. We present a time-dependent spectral energy distribution of the bright FSRQ 3C 454.3 during its December 2009 flare, which is well fit by an external Compton model in the bright state, although day to day changes pose challenges to a simple one-zone model. All data from the SMARTS monitoring program are publicly available on our website.

Significant progress has been made in the last few years on understanding how supermassive black holes form and grow. In this paper, we begin by reviewing the spectral signatures of Active Galactic Nuclei (AGN) ranging from radio to hard X-ray wavelengths. We then describe the most commonly used methods to find these sources, including optical/UV, radio, infrared and X-ray emission and optical emission lines. We then describe the main observational properties of the obscured and unobscured AGN population. Finally, we summarize the cosmic history of black hole accretion, i.e., when in the history of the Universe supermassive black holes were getting most of their mass. We finish with a summary of open questions and a description of planned and future observatories that are going to help answer them.

We compute black hole masses and bolometric luminosities for 87 obscured AGN in the redshift range 0.25 < z < 1.25, selected from the GOODS deep multi-wavelength survey fields via their X-ray emission. We fit the optical images and obtain morphological parameters for the host galaxy, separating the galaxy from its central point source, thereby obtaining a four-band optical SED for each active nucleus. We calculate bolometric luminosities for these AGN by reddening a normalized mean SED of GOODS broad-line AGN to match the observed central point-source SED of each obscured AGN. This estimate of Lbol has a smaller spread than simple bolometric corrections to the X-ray luminosity or direct integration of the observed multi-wavelength SED, suggesting it is a better measure. We estimate central black hole masses from the bulge luminosities. The black hole masses span a wide range, 7 x 10^6 M_sun to 6 x 10^9 M_sun; the median black hole mass is 5 x 10^8 M_sun. The majority of these AGN have L/L_Edd < 0.01, and we detect no significant evolution of the mean Eddington ratio to z = 1.25. This implies that the bulk of black hole growth in these obscured AGN must have occurred at z > 1 and that we are observing these AGN in a slow- or no-growth state.

We show how accretion rate governs the physical properties of a sample of unobscured broad-line, narrow-line, and lineless active galactic nuclei (AGNs). We avoid the systematic errors plaguing previous studies of AGN accretion rate by using accurate accretion luminosities (L_int) from well-sampled multiwavelength SEDs from the Cosmic Evolution Survey (COSMOS), and accurate black hole masses derived from virial scaling relations (for broad-line AGNs) or host-AGN relations (for narrow-line and lineless AGNs). In general, broad emission lines are present only at the highest accretion rates (L_int/L_Edd > 0.01), and these rapidly accreting AGNs are observed as broad-line AGNs or possibly as obscured narrow-line AGNs. Narrow-line and lineless AGNs at lower specific accretion rates (L_int/L_Edd < 0.01) are unobscured and yet lack a broad line region. The disappearance of the broad emission lines is caused by an expanding radiatively inefficient accretion flow (RIAF) at the inner radius of the accretion disk. The presence of the RIAF also drives L_int/L_Edd < 10^-2 narrow-line and lineless AGNs to 10 times higher ratios of radio to optical/UV emission than L_int/L_Edd > 0.01 broad-line AGNs, since the unbound nature of the RIAF means it is easier to form a radio outflow. The IR torus signature also tends to become weaker or disappear from L_int/L_Edd < 0.01 AGNs, although there may be additional mid-IR synchrotron emission associated with the RIAF. Together these results suggest that specific accretion rate is an important physical "axis" of AGN unification, described by a simple model.

We present the rest-frame optical morphologies of active galactic nucleus (AGN) host galaxies at 1.5<z<3, using near-infrared imaging from the Hubble Space Telescope Wide Field Camera 3, the first such study of AGN host galaxies at these redshifts. The AGN are X-ray selected from the Chandra Deep Field South and have typical luminosities of 1E42 < L_X < 1E44 erg/s. Accreting black holes in this luminosity and redshift range account for a substantial fraction of the total space density and black hole mass growth over cosmic time; they thus represent an important mode of black hole growth in the universe. We find that the majority (~80%) of the host galaxies of these AGN have low Sersic indices indicative of disk-dominated light profiles, suggesting that secular processes govern a significant fraction of the cosmic growth of black holes. That is, many black holes in the present-day universe grew much of their mass in disk-dominated galaxies and not in early-type galaxies or major mergers. The properties of the AGN host galaxies are furthermore indistinguishable from their parent galaxy population and we find no strong evolution in either effective radii or morphological mix between z~2 and z~0.05.

Galaxy formation is significantly modulated by energy output from supermassive black holes at the centers of galaxies which grow in highly efficient luminous quasar phases. The timescale on which black holes transition into and out of such phases is, however, unknown. We present the first measurement of the shutdown timescale for an individual quasar using X-ray observations of the nearby galaxy IC 2497, which hosted a luminous quasar no more than 70,000 years ago that is still seen as a light echo in `Hanny's Voorwerp', but whose present-day radiative output is lower by at least 2 and more likely by over 4 orders of magnitude. This extremely rapid shutdown provides new insights into the physics of accretion in supermassive black holes, and may signal a transition of the accretion disk to a radiatively inefficient state.

We take advantage of the rich multi-wavelength data available in the Chandra Deep Field South (CDF-S), including the 4 Msec Chandra observations (the deepest X-ray data to date), in order to search for heavily-obscured low-luminosity AGN among infrared-luminous galaxies. In particular, we obtained a stacked rest-frame X-ray spectrum for samples of galaxies binned in terms of their IR luminosity or stellar mass. We detect a significant signal at E~1 to 8 keV, which we interpret as originating from a combination of emission associated with star-formation processes at low energies combined with a heavily-obscured AGN at E>5 keV. We further find that the relative strength of this AGN signal decays with decreasing IR luminosity, indicating a higher AGN fraction for more luminous IR sources. Together, these results strongly suggest the presence of a large number of obscured AGN in IR-luminous galaxies. Using samples binned in terms of stellar mass in the host galaxy, we find a significant excess at E=6-7 keV for sources with M>10^{11} Msun, consistent with a large obscured AGN population in high mass galaxies. In contrast, no strong evidence of AGN activity was found for less-massive galaxies. The integrated intensity at high energies indicates that a significant fraction of the total black hole growth, ~22%, occurs in heavily-obscured systems that are not individually detected in even the deepest X-ray observations. There are also indications that the number of low-luminosity, heavily-obscured AGN does not evolve significantly with redshift, in contrast to the strong evolution seen in higher luminosity sources.

Using new, highly accurate photometric redshifts from the MUSYC medium-band survey in the Extended Chandra Deep Field South (ECDF-S), we fit synthetic stellar population models to compare AGN host galaxies to inactive galaxies at 0.8 < z < 1.2. We find that AGN host galaxies are predominantly massive galaxies on the red sequence and in the green valley of the color-mass diagram. Because both passive and dusty galaxies can appear red in optical colors, we use rest-frame near-infrared colors to separate passively evolving stellar populations from galaxies that are reddened by dust. As with the overall galaxy population, ~25% of the `red' AGN host galaxies and ~75% of the `green' AGN host galaxies have colors consistent with young stellar populations reddened by dust. The dust-corrected rest-frame optical colors are the blue colors of star-forming galaxies, which implies that these AGN hosts are not passively aging to the red sequence. At z~1, AGN activity is roughly evenly split between two modes of black hole growth: the first in passively evolving host galaxies, which may be heating up the galaxy's gas and preventing future episodes of star formation, and the second in dust-reddened young galaxies, which may be ionizing the galaxy's interstellar medium and shutting down star formation.

We present deep optical 18-medium-band photometry from the Subaru telescope over the ~30' x 30' Extended Chandra Deep Field-South (ECDF-S), as part of the Multiwavelength Survey by Yale-Chile (MUSYC). This field has a wealth of ground- and space-based ancillary data, and contains the GOODS-South field and the Hubble Ultra Deep Field. We combine the Subaru imaging with existing UBVRIzJHK and Spitzer IRAC images to create a uniform catalog. Detecting sources in the MUSYC BVR image we find ~40,000 galaxies with R_AB<25.3, the median 5 sigma limit of the 18 medium bands. Photometric redshifts are determined using the EAZY code and compared to ~2000 spectroscopic redshifts in this field. The medium band filters provide very accurate redshifts for the (bright) subset of galaxies with spectroscopic redshifts, particularly at 0.1 < z < 1.2 and at z > 3.5. For 0.1 < z < 1.2, we find a 1 sigma scatter in \Delta z/(1+z) of 0.007, similar to results obtained with a similar filter set in the COSMOS field. As a demonstration of the data quality, we show that the red sequence and blue cloud can be cleanly identified in rest-frame color-magnitude diagrams at 0.1 < z < 1.2. We find that ~20% of the red-sequence-galaxies show evidence of dust-emission at longer rest-frame wavelengths. The reduced images, photometric catalog, and photometric redshifts are provided through the public MUSYC website.

A Chandra X-ray imaging observation of the jet in Pictor A showed a feature that appears to be a flare that faded between 2000 and 2002. The feature was not detected in a follow-up observation in 2009. The jet itself is over 150 kpc long and a kpc wide, so finding year-long variability is surprising. Assuming a synchrotron origin of the observed high-energy photons and a minimum energy condition for the outflow, the synchrotron loss time of the X-ray emitting electrons is of order 1200 yr, which is much longer than the observed variability timescale. This leads to the possibility that the variable X-ray emission arises from a very small sub-volume of the jet, characterized by magnetic field that is substantially larger than the average over the jet.

Despite observed strong correlations between central supermassive black holes (SMBHs) and star-formation in galactic nuclei, uncertainties exist in our understanding of their coupling. We present observations of the ratio of heavily-obscured to unobscured quasars as a function of cosmic epoch up to z~3, and show that a simple physical model describing mergers of massive, gas-rich galaxies matches these observations. In the context of this model, every obscured and unobscured quasar represent two distinct phases that result from a massive galaxy merger event. Much of the mass growth of the SMBH occurs during the heavily-obscured phase. These observations provide additional evidence for a causal link between gas-rich galaxy mergers, accretion onto the nuclear SMBH and coeval star formation.

Models of galaxy formation invoke the major merger of gas-rich progenitor galaxies as the trigger for significant phases of black hole growth and the associated feedback that suppresses star formation to create red spheroidal remnants. However, the observational evidence for the connection between mergers and active galactic nucleus (AGN) phases is not clear. We analyze a sample of low-mass early-type galaxies known to be in the process of migrating from the blue cloud to the red sequence via an AGN phase in the green valley. Using deeper imaging from SDSS Stripe 82, we show that the fraction of objects with major morphological disturbances is high during the early starburst phase, but declines rapidly to the background level seen in quiescent early-type galaxies by the time of substantial AGN radiation several hundred Myr after the starburst. This observation empirically links the AGN activity in low-redshift early-type galaxies to a significant merger event in the recent past. The large time delay between the merger-driven starburst and the peak of AGN activity allows for the merger features to decay to the background and hence may explain the weak link between merger features and AGN activity in the literature.

We use data from large surveys of the local Universe (SDSS+Galaxy Zoo) to show that the galaxy-black hole connection is linked to host morphology at a fundamental level. The fraction of early-type galaxies with actively growing black holes, and therefore the AGN duty cycle, declines significantly with increasing black hole mass. Late-type galaxies exhibit the opposite trend: the fraction of actively growing black holes increases with black hole mass.

We use data from the Sloan Digital Sky Survey and visual classifications of morphology from the Galaxy Zoo project to study black hole growth in the nearby Universe (z < 0.05) and to break down the AGN host galaxy population by color, stellar mass and morphology. We find that black hole growth at luminosities L_OIII >1E40 erg/s in early- and late-type galaxies is fundamentally different. AGN host galaxies as a population have a broad range of stellar masses (1E10-1E11 Msun), reside in the green valley of the color-mass diagram and their central black holes have median masses around 1E6.5 Msun. However, by comparing early- and late-type AGN host galaxies to their non-active counterparts, we find several key differences: in early-type galaxies, it is preferentially the galaxies with the least massive black holes that are growing, while late-type galaxies, it is preferentially the most massive}black holes that are growing. The duty cycle of AGN in early-type galaxies is strongly peaked in the green valley below the low-mass end (1E10 Msun) of the red sequence at stellar masses where there is a steady supply of blue cloud progenitors. The duty cycle of AGN in late-type galaxies on the other hand peaks in massive (1E11 Msun) green and red late-types which generally do not have a corresponding blue cloud population of similar mass. At high Eddington ratios (L/L_Edd > 0.1), the only population with a substantial fraction of AGN are the low-mass green valley early-type galaxies. Finally, the Milky Way likely resides in the "sweet spot" on the color-mass diagram where the AGN duty cycle of late-type galaxies is highest. We discuss the implications of these results for our understanding of the role of AGN in the evolution of galaxies

We constrain the number density and evolution of Compton-thick Active Galactic Nuclei (AGN), and their contribution to the extragalactic X-ray background. In the local Universe we use the wide area surveys from the Swift and INTEGRAL satellites, while for high redshifts we explore candidate selections based on mid-IR parameters. We present the properties of a sample of 211 heavily-obscured AGN candidates in the Extended Chandra Deep Field-South (ECDF-S) selecting objects with f24/fR>1000 and R-K>4.5. The X-ray to mid-IR ratios for these sources are significantly larger than that of star-forming galaxies and ~2 orders of magnitude smaller than for the general AGN population, suggesting column densities of NH>5x10^24 cm^-2. The space density of CT AGN at z~2 derived from these observations is ~10^-5 Mpc^{-3}, finding a strong evolution in the number of LX>10^44 erg/s sources from z=1.5 to 2.5.

We study the properties of a sample of 211 heavily-obscured Active Galactic Nucleus (AGN) candidates in the Extended Chandra Deep Field-South selecting objects with f_24/f_R>1000 and R-K>4.5. Of these, 18 were detected in X-rays and found to be obscured AGN with neutral hydrogen column densities of ~10^23 cm^-2. In the X-ray undetected sample, the following evidence suggests a large fraction of heavily-obscured (Compton Thick) AGN: (i) The stacked X-ray signal of the sample is strong, with an observed ratio of soft to hard X-ray counts consistent with a population of ~90% heavily obscured AGN combined with 10% star-forming galaxies. (ii) The X-ray to mid-IR ratios for these sources are significantly larger than that of star-forming galaxies and ~2 orders of magnitude smaller than for the general AGN population, suggesting column densities of N_H>5x10^24 cm^-2. (iii) The Spitzer near- and mid-IR colors of these sources are consistent with those of the X-ray-detected sample if the effects of dust self-absorption are considered. Spectral fitting to the rest-frame UV/optical light (dominated by the host galaxy) returns stellar masses of ~10^11 M_sun and <E(B-V)> =0.5, and reveals evidence for a significant young stellar population, indicating that these sources are experiencing considerable star-formation. This sample of heavily-obscured AGN candidates implies a space density at z~2 of ~10^-5 Mpc^-3, finding a strong evolution in the number of L_X>10^44 erg/s sources from z=1.5 to 2.5, possibly consistent with a short-lived heavily-obscured phase before an unobscured quasar is visible.

Following the recent discovery of gamma rays from the radio-loud narrow-line Seyfert 1 galaxy PMN J0948+0022 (z=0.5846), we started a multiwavelength campaign from radio to gamma rays, which was carried out between the end of March and the beginning of July 2009. The source displayed activity at all the observed wavelengths: a general decreasing trend from optical to gamma-ray frequencies was followed by an increase of radio emission after less than two months from the peak of the gamma-ray emission. The largest flux change, about a factor of about 4, occurred in the X-ray band. The smallest was at ultraviolet and near-infrared frequencies, where the rate of the detected photons dropped by a factor 1.6-1.9. At optical wavelengths, where the sampling rate was the highest, it was possible to observe day-scale variability, with flux variations up to a factor of about 3. The behavior of PMN J0948+0022 observed in this campaign and the calculated power carried out by its jet in the form of protons, electrons, radiation and magnetic field are quite similar to that of blazars, specifically of flat-spectrum radio quasars. These results confirm the idea that radio-loud narrow-line Seyfert 1 galaxies host relativistic jets with power similar to that of average blazars.

The emission mechanisms in extragalactic jets include synchrotron and various inverse-Compton processes. At low (radio through infrared) energies, it is widely agreed that synchrotron emission dominates in both low-power (FR I) and high-power (FR II and quasar) jets, because of the power-law nature of the spectra observed and high polarizations. However, at higher energies, the emission mechanism for high-power jets at kpc scales is hotly debated. Two mechanisms have been proposed: either inverse-Compton of cosmic microwave background photons or synchrotron emission from a second, high-energy population of electrons. Here we discuss optical polarimetry as a method for diagnosing the mechanism for the high-energy emission in quasar jets, as well as revealing the jet's three-dimensional energetic and magnetic field structure. We then discuss high-energy emission mechanisms for powerful jets in the light of the HST polarimetry of PKS 1136-135.

We investigate a class of rapidly growing emission line galaxies, known as "Green Peas", first noted by volunteers in the Galaxy Zoo project because of their peculiar bright green colour and small size, unresolved in SDSS imaging. Their appearance is due to very strong optical emission lines, namely [O III] 5007 A, with an unusually large equivalent width of up to ~1000 A. We discuss a well-defined sample of 251 colour-selected objects, most of which are strongly star forming, although there are some AGN interlopers including 8 newly discovered narrow Line Seyfert 1 galaxies. The star-forming Peas are low mass galaxies (M~10^8.5 - 10^10 M_sun) with high star formation rates (~10 M_sun/yr), low metallicities (log[O/H] + 12 ~ 8.7) and low reddening (E(B-V) < 0.25) and they reside in low density environments. They have some of the highest specific star formation rates (up to ~10^{-8} yr^{-1}) seen in the local Universe, yielding doubling times for their stellar mass of hundreds of Myrs. The few star-forming Peas with HST imaging appear to have several clumps of bright star-forming regions and low surface density features that may indicate recent or ongoing mergers. The Peas are similar in size, mass, luminosity and metallicity to Luminous Blue Compact Galaxies. They are also similar to high redshift UV-luminous galaxies, e.g., Lyman-break galaxies and Lyman-alpha emitters, and therefore provide a local laboratory with which to study the extreme star formation processes that occur in high-redshift galaxies. Studying starbursting galaxies as a function of redshift is essential to understanding the build up of stellar mass in the Universe.

The Chandra COSMOS Survey (C-COSMOS) is a large, 1.8 Ms, Chandra} program that has imaged the central 0.5 sq.deg of the COSMOS field (centered at 10h, +02deg) with an effective exposure of ~160ksec, and an outer 0.4sq.deg. area with an effective exposure of ~80ksec. The limiting source detection depths are 1.9e-16 erg cm(-2) s(-1) in the Soft (0.5-2 keV) band, 7.3e(-16) erg cm^-2 s^-1 in the Hard (2-10 keV) band, and 5.7e(-16) erg cm(-2) s(-1) in the Full (0.5-10 keV) band. Here we describe the strategy, design and execution of the C-COSMOS survey, and present the catalog of 1761 point sources detected at a probability of being spurious of <2e(-5) (1655 in the Full, 1340 in the Soft, and 1017 in the Hard bands). By using a grid of 36 heavily (~50%) overlapping pointing positions with the ACIS-I imager, a remarkably uniform (to 12%) exposure across the inner 0.5 sq.deg field was obtained, leading to a sharply defined lower flux limit. The widely different PSFs obtained in each exposure at each point in the field required a novel source detection method, because of the overlapping tiling strategy, which is described in a companion paper. (Puccetti et al. Paper II). This method produced reliable sources down to a 7-12 counts, as verified by the resulting logN-logS curve, with sub-arcsecond positions, enabling optical and infrared identifications of virtually all sources, as reported in a second companion paper (Civano et al. Paper III). The full catalog is described here in detail, and is available on-line.

We present a simple formalism to interpret two galaxy statistics, the UV luminosity function and two-point correlation functions for star-forming galaxies at z~4, 5, 6 in the context of LCDM cosmology. Both statistics are the result of how star formation takes place in DM halos, and thus are used to constrain how UV light depends on halo properties such as mass. The two measures were taken from the GOODS data, thus ideal for joint analysis. The two physical quantities we explore are the SF duty cycle, and the range of L_UV that a halo of mass M can have (mean and variance). The former addresses the typical duration of SF activity in halos while the latter addresses the averaged SF history and regularity of gas inflow into these systems. We explore various physical models consistent with data, and find the following: 1) the typical duration of SF observed in the data is <0.4 Gyr (1 sig), 2) the inferred scaling law between L_UV and halo mass M from the observed slope of the LFs is roughly linear at all redshifts, and 3) L_UV for a fixed halo mass decreases with time, implying that the SF efficiency (after dust extinction) is higher at earlier times. We explore several physical scenarios relating star formation to halo mass, but find that these scenarios are indistinguishable due to the limited range of halo mass probed by our data. In order to discriminate between different scenarios, we discuss constraining the bright-faint galaxy cross-correlation functions and luminosity-dependence of galaxy bias. (Abridged)

We present a new K-selected, optical-to-near-infrared photometric catalog of the Extended Chandra Deep Field South (ECDFS), making it publicly available to the astronomical community. The dataset is founded on publicly available imaging, supplemented by original zJK imaging data obtained as part of the MUltiwavelength Survey by Yale-Chile (MUSYC). The final photometric catalog consists of photometry derived from nine band U-K imaging covering the full 0.5x0.5 sq. deg. of the ECDFS, plus H band data for approximately 80% of the field. The 5sigma flux limit for point-sources is K = 22.0 (AB). This is also the nominal completeness and reliability limit of the catalog: the empirical completeness for 21.75 < K < 22.00 is 85+%. We have verified the quality of the catalog through both internal consistency checks, and comparisons to other existing and publicly available catalogs. As well as the photometric catalog, we also present catalogs of photometric redshifts and restframe photometry derived from the ten band photometry. We have collected robust spectroscopic redshift determinations from published sources for 1966 galaxies in the catalog. Based on these sources, we have achieved a (1sigma) photometric redshift accuracy of Dz/(1+z) = 0.036, with an outlier fraction of 7.8%. Most of these outliers are X-ray sources. Finally, we describe and release a utility for interpolating restframe photometry from observed SEDs, dubbed InterRest. Particularly in concert with the wealth of already publicly available data in the ECDFS, this new MUSYC catalog provides an excellent resource for studying the changing properties of the massive galaxy population at z < 2. (Abridged)

We constrain the number density and evolution of Compton-thick Active Galactic Nuclei (AGN). In the local Universe we use the wide area surveys from the Swift and INTEGRAL satellites, while for high redshifts we explore candidate selections based on a combination of X-ray and mid-IR parameters. We find a significantly lower space density of Compton-thick AGN in the local Universe than expected from published AGN population synthesis models to explain the X-ray background. This can be explained by the numerous degeneracies in the parameters of those models; we use the high-energy surveys described here to remove those degeneracies. We show that only direct observations of CT AGN can currently constrain the number of heavily-obscured supermassive black holes. At high redshift, the inclusion of IR-selected Compton-thick AGN candidates leads to a much higher space density, implying (a) a different (steeper) evolution for these sources compared to less-obscured AGN, (b) that the IR selection includes a large number of interlopers, and/or (c) that there is a large number of reflection-dominated AGN missed in the INTEGRAL and Swift observations. The contribution of CT AGN to the X-ray background is small, ~9%, with a comparable contribution to the total cosmic accretion, unless reflection-dominated CT AGN significantly outnumber transmission-dominated CT AGN, in which case their contribution can be much higher. Using estimates derived here for the accretion luminosity over cosmic time we estimate the local mass density in supermassive black holes and find a good agreement with available constraints for an accretion efficiency of ~10%. Transmission-dominated CT AGN contribute only ~8% to total black hole growth.

The COSMOS survey is a multiwavelength survey aimed to study the evolution of galaxies, AGN and the large scale structure. The XMM-COSMOS is a deep X-ray survey over the full 2 deg2 of the COSMOS area. It consists of 55 XMM-Newton pointings for a total exposure of ~1.5 Ms with an average vignetting corrected depth of 40 ks across the field of view and a sky coverage of 2.13 deg2. We present the catalogue of point-like X-ray sources detected with the EPIC CCD cameras, the logN-logS relations and the X-ray colour-colour diagrams. The analysis was performed in the 0.5-2 keV, 2-10 keV and 5-10 keV energy bands. The completeness of the catalogue as well as logN-logS have been calibrated using Monte Carlo simulations. The catalogs contains a total of 1887 unique sources detected in at least one band. The survey, that shows unprecedented homogeneity, has a flux limit of ~1.7x10-15 erg cm-2 s-1, ~9.3x10-15 erg cm-2 s-1 and ~1.3x10-14 erg cm-2 s-1 over 90% of the area (1.92 deg2) in the 0.5-2 keV, 2-10 keV and 5-10 keV energy band, respectively. Thanks to the rather homogeneous exposure over a large area, the derived logN-logS relations are very well determined over the flux range sampled by XMM-COSMOS. These relations have been compared with XRB synthesis models, which reproduce the observations with an agreement of ~10% in the 5-10 keV and 2-10 keV band, while in the 0.5-2 keV band the agreement is of the order of ~20%. The hard X-ray colors confirmed that the majority of the extragalactic sources, in a bright subsample, are actually Type I or Type II AGN. About 20% of the sources have X-ray luminosity typical of AGN (L_X >1042 erg/s) although they do not show any clear signature of nuclear activity in the optical spectrum.

The growth of supermassive black holes and their host galaxies are thought to be linked, but the precise nature of this symbiotic relationship is still poorly understood. Both observations and simulations of galaxy formation suggest that the energy input from active galactic nuclei (AGN), as the central supermassive black hole accretes material and grows, heats the interstellar material and suppresses star formation. In this Letter, we show that most host galaxies of moderate-luminosity supermassive black holes in the local universe have intermediate optical colors that imply the host galaxies are transitioning from star formation to quiescence, the first time this has been shown to be true for all AGN independent of obscuration. The intermediate colors suggest that star formation in the host galaxies ceased roughly 100 Myr ago. This result indicates that either the AGN are very long-lived, accreting for more than 1 Gyr beyond the end of star formation, or there is a ~100 Myr time delay between the shutdown of star formation and detectable black hole growth. The first explanation is unlikely given current estimates for AGN lifetimes, so low-lumiosity AGN must shut down star formation before substantial black hole accretion activity is detected. The scarcity of AGN host galaxies in the blue cloud reported here challenges scenarios where significant star formation and black hole growth are coeval. Lastly, these observations also strongly support the `Unified Model' of AGN as the host galaxy colors are independent of obscuration towards the central engine.

The blazar 3C 454.3 was revealed by the Fermi Gamma-ray Space Telescope to be in an exceptionally high flux state in July 2008. Accordingly, we performed a multi-wavelength monitoring campaign on this blazar using IR and optical observations from the SMARTS telescopes, optical, UV and X-ray data from the Swift satellite, and public-release gamma-ray data from Fermi. We find an excellent correlation between the IR, optical, UV and gamma-ray light curves, with a time lag of less than one day. The amplitude of the infrared variability is comparable to that in gamma-rays, and larger than at optical or UV wavelengths. The X-ray flux is not strongly correlated with either the gamma-rays or longer wavelength data. These variability characteristics find a natural explanation in the external Compton model, in which electrons with Lorentz factor gamma~10^(3-4) radiate synchrotron emission in the infrared-optical and also scatter accretion disk or emission line photons to gamma-ray energies, while much cooler electrons (gamma~10^(1-2)) produce X-rays by scattering synchrotron or other ambient photons.

We present the first results of our optical spectroscopy program aimed to provide redshifts and identifications for the X-ray sources in the Extended Chandra Deep Field South. A total of 339 sources were targeted using the IMACS spectrograph at the Magellan telescopes and the VIMOS spectrograph at the VLT. We measured redshifts for 186 X-ray sources, including archival data and a literature search. We find that the AGN host galaxies have on average redder rest-frame optical colors than non-active galaxies, and that they live mostly in the "green valley". The dependence of the fraction of AGN that are obscured on both luminosity and redshift is confirmed at high significance and the observed AGN space density is compared with the expectations from existing luminosity functions. These AGN show a significant difference in the mid-IR to X-ray flux ratio for obscured and unobscured AGN, which can be explained by the effects of dust self-absorption on the former. This difference is larger for lower luminosity sources, which is consistent with the dust opening angle depending on AGN luminosity.

We present the color-magnitude and color-stellar mass diagrams for galaxies with z_phot < ~2, based on a K < 22 (AB) catalog of the Extended Chandra Deep Field South (ECDFS) from the MUltiwavelength Survey by Yale-Chile (MUSYC). Our main sample of 7840 galaxies contains 1297 M_* > 10^11 M_Sol galaxies in the range 0.2 < z_phot < 1.8. We show empirically that this catalog is approximately complete for M_* > 10^11 M_Sol galaxies for z_phot < 1.8. For this mass-limited sample, we show that the locus of the red sequence color-stellar mass relation evolves as Del(u-r) ~ (-0.44+/-0.02) z_phot for z_phot < ~1.2. For z_phot > ~1.3, however, we are no longer able to reliably distinguish red and blue subpopulations from the observed color distribution; we show that this would require much deeper near infrared data. At 1.5 < z_phot <1.8, the comoving number density of M_* > 10^11 M_Sol galaxies is ~50% of the local value, with a red fraction of ~33%. Making a parametric fit to the observed evolution, we find n_tot(z) ~ (1+z_phot)^(-0.52+/-0.12(+/-0.20)). We find stronger evolution in the red fraction: f_red(z) ~ (1+z_phot)^(-1.17+/-0.18(+/-0.21)). Through a series of sensitivity analyses, we show that the most important sources of systematic error are: 1. systematic differences in the analysis of the z~0 and z>>0 samples; 2. systematic effects associated with details of the photometric redshift calculation; and 3. uncertainties in the photometric calibration. With this in mind, we show that our results based on photometric redshifts are consistent with a completely independent analysis which does not require redshift information for individual galaxies. Our results suggest that, at most, 1/5 of local red sequence galaxies with M_* >10^11 M_Sol were already in place at z ~ 2.

We present an analysis of the radio properties of large samples of Lyman Break Galaxies (LBGs) at $z \sim 3$, 4, and 5 from the COSMOS field. The median stacking analysis yields a statistical detection of the $z \sim 3$ LBGs (U-band drop-outs), with a 1.4 GHz flux density of $0.90 \pm 0.21 \mu$Jy. The stacked emission is unresolved, with a size $< 1"$, or a physical size $< 8$kpc. The total star formation rate implied by this radio luminosity is $31\pm 7$ $M_\odot$ year$^{-1}$, based on the radio-FIR correlation in low redshift star forming galaxies. The star formation rate derived from a similar analysis of the UV luminosities is 17 $M_\odot$ year$^{-1}$, without any correction for UV dust attenuation. The simplest conclusion is that the dust attenuation factor is 1.8 at UV wavelengths. However, this factor is considerably smaller than the standard attenuation factor $\sim 5$, normally assumed for LBGs. We discuss potential reasons for this discrepancy, including the possibility that the dust attenuation factor at $z \ge 3$ is smaller than at lower redshifts. Conversely, the radio luminosity for a given star formation rate may be systematically lower at very high redshift. Two possible causes for a suppressed radio luminosity are: (i) increased inverse Compton cooling of the relativistic electron population due to scattering off the increasing CMB at high redshift, or (ii) cosmic ray diffusion from systematically smaller galaxies. The radio detections of individual sources are consistent with a radio-loud AGN fraction of 0.3%. One source is identified as a very dusty, extreme starburst galaxy (a 'submm galaxy').

In order to assess the accuracy with which we can determine the morphologies of AGN host galaxies, we have simulated more than 50,000 ACS images of galaxies with z < 1.25, using image and noise properties appropriate for the GOODS survey. We test the effect of central point-source brightness on host galaxy parameter recovery with a set of simulated AGN host galaxies made by adding point sources to the centers of normal galaxies. We extend this analysis and also quantify the recovery of intrinsic morphological parameters of AGN host galaxies with a set of fully simulated inactive and AGN host galaxies. We can reliably separate good from poor fit results using a combination of reasonable error cuts, in the regime where L_{host}:L_{PS} > 1:4. We give quantitative estimates of parameter errors as a function of host-to-point-source ratio. In general, we separate host and point-source magnitudes reliably at all redshifts; point sources are well recovered more than 90% of the time, although spurious detection of central point sources can be as high as 25% for bulge-dominated sources. We find a general correlation between Sersic index and intrinsic bulge-to-total ratio, such that a host galaxy with Sersic n < 1.5 generally has at least 80% of its light from a disk component. Likewise, "bulge-dominated" galaxies with n > 4 typically derive at least 70% of their total host galaxy light from a bulge, but this number can be as low as 55%. Single-component Sersic fits to an AGN host galaxy are statistically very reliable to z < 1.25 (for ACS survey data like ours). In contrast, two-component fits involving separate bulge and disk components tend to over-estimate the bulge fraction by ~10%, with uncertainty of order 50%.

We present the mid-infrared colors of X-ray-detected AGN and explore mid-infrared selection criteria. Using a statistical matching technique, the likelihood ratio, over 900 IRAC counterparts were identified with a new MUSYC X-ray source catalog that includes ~1000 published X-ray sources in the Chandra Deep Field-South and Extended Chandra Deep Field-South. Most X-ray-selected AGN have IRAC spectral shapes consistent with power-law slopes, f_{nu} ~ nu^{alpha}, and display a wide range of colors, -2 < alpha < 2. Although X-ray sources typically fit to redder (more negative alpha) power-laws than non-X-ray detected galaxies, more than 50% do have flat or blue (galaxy-like) spectral shapes in the observed 3-8 micron band. Only a quarter of the X-ray selected AGN detected at 24 micron are well fit by featureless red power laws in the observed 3.6-24 micron, likely the subset of our sample whose infrared spectra are dominated by emission from the central AGN region. Most IRAC color-selection criteria fail to identify the majority of X-ray-selected AGN, finding only the more luminous AGN, the majority of which have broad emission lines. In deep surveys, these color-selection criteria select 10-20% of the entire galaxy population and miss many moderate luminosity AGN.

In recent years deep X-ray and infrared surveys have provided an efficient way to find accreting supermassive black holes, otherwise known as active galactic nuclei (AGN), in the young universe. Such surveys can, unlike optical surveys, find AGN obscured by high column densities of gas and dust. In those cases, deep optical data show only the host galaxy, which can then be studied in greater detail than in unobscured AGN. Some years ago the hard spectrum of the X-ray "background" suggested that most AGN were obscured. Now GOODS, MUSYC, COSMOS and other surveys have confirmed this picture and given important quantitative constraints on AGN demographics. Specifically, we show that most AGN are obscured at all redshifts and the amount of obscuration depends on both luminosity and redshift, at least out to redshift z~2, the epoch of substantial black holes and galaxy growth. Larger-area deep infrared and hard X-ray surveys will be needed to reach higher redshifts and to probe fully the co-evolution of galaxies and black holes.

We present the first clustering results of X-ray selected AGN at z~3. Using Chandra X-ray imaging and UVR optical colors from MUSYC photometry in the ECDF-S field, we selected a sample of 58 z~3 AGN candidates. From the optical data we also selected 1385 LBG at 2.8<z< 3.8 with R<25.5. We performed auto-correlation and cross-correlation analyses, and here we present results for the clustering amplitudes and dark matter halo masses of each sample. For the LBG we find a correlation length of r_0,LBG = 6.7 +/- 0.5 Mpc, implying a bias value of 3.5 +/- 0.3 and dark matter (DM) halo masses of log(Mmin/Msun) = 11.8 +/- 0.1. The AGN-LBG cross-correlation yields r_0,AGN-LBG = 8.7 +/- 1.9 Mpc, implying for AGN at 2.8<z<3.8 a bias value of 5.5 +/- 2.0 and DM halo masses of log(Mmin/Msun) = 12.6 +0.5/-0.8. Evolution of dark matter halos in the Lambda CDM cosmology implies that today these z~3 AGN are found in high mass galaxies with a typical luminosity of 7+4/-2 L*.

We studied the clustering properties and multiwavelength spectral energy distributions of a complete sample of 162 Ly Alpha-Emitting (LAE) galaxies at z=3.1 discovered in deep narrow-band MUSYC imaging of the Extended Chandra Deep Field South. LAEs were selected to have observed frame equivalent widths >80A and emission line fluxes >1.5E-17 erg/cm^2/s. Only 1% of our LAE sample appears to host AGN. The LAEs exhibit a moderate spatial correlation length of r_0=3.6+0.8-1.0 Mpc, corresponding to a bias factor b=1.7+0.3-0.4, which implies median dark matter halo masses of log10(M_med) = 10.9+0.5-0.9 M_sun. Comparing the number density of LAEs, (1.5+-0.3)E-3/Mpc^3, with the number density of these halos finds a mean halo occupation ~1-10%. The evolution of galaxy bias with redshift implies that most z=3.1 LAEs evolve into present-day galaxies with L<2.5L*, whereas other z>3 galaxy populations typically evolve into more massive galaxies. Halo merger trees show that z=0 descendants occupy halos with a wide range of masses, with a median descendant mass close to that of L*. Only 30% of LAEs have sufficient stellar mass (>~3E9 M_sun) to yield detections in deep Spitzer-IRAC imaging. A two-population SED fit to the stacked UBVRIzJK+[3.6,4.5,5.6,8.0]micron fluxes of the IRAC-undetected objects finds that the typical LAE has low stellar mass (1.0+0.6-0.4 E9 M_sun), moderate star formation rate (2+-1 M_sun/yr), a young component age of 20+30-10 Myr, and little dust (A_V<0.2). The best fit model has 20% of the mass in the young stellar component, but models without evolved stars are also allowed.

We combine deep UBVRIzJK photometry from the Multiwavelength Survey by Yale-Chile (MUSYC) with redshifts from the COMBO-17 survey to perform a large-scale study of the rest-frame ultraviolet (UV) properties of 674 high-redshift (0.5<z<1) early-type galaxies, drawn from the Extended Chandra Deep Field South (E-CDFS). Galaxy morphologies are determined through visual inspection of Hubble Space Telescope (HST) images taken from the GEMS survey. We harness the sensitivity of the UV to young (<1 Gyr old) stars to quantify the recent star formation history of early-type galaxies across a range of luminosities (-23.5 < M(V) < -18). Comparisons to simple stellar populations forming at high redshift indicate that only ~1.1 percent of early-types in this sample are consistent with purely passive ageing since z=2. Parametrising the recent star formation (RSF) in terms of the mass fraction of stars less than a Gyr old, we find that the early-type population as a whole shows a typical RSF between 5 and 13% in the redshift range 0.5<z<1. Early-types on the UV red sequence show RSF values less than 5% while the reddest early-types are virtually quiescent with RSF values of ~1%. We find compelling evidence that early-types of all luminosities form stars over the lifetime of the Universe, although the bulk of their star formation is already complete at high redshift. This tail-end of star formation is measurable and not negligible, with luminous (-23<M(V)<-20.5) early-types potentially forming 10-15% of their mass since z=1, with their less luminous (M(V)>-20.5) counterparts potentially forming 30-60 percent of their mass in the same redshift range. (abridged)

Using the Gemini Near-InfraRed Spectrograph (GNIRS), we have assembled a complete sample of 20 K-selected galaxies at 2.0<z<2.7 with high quality near-infrared spectra. As described in a previous paper, 9 of these 20 galaxies have strongly suppressed star formation and no detected emission lines. The present paper concerns the 11 galaxies with detected Halpha emission, and studies the origin of the line emission using the GNIRS spectra and follow-up observations with SINFONI on the VLT. Based on their [NII]/Halpha ratios, the spatial extent of the line emission and several other diagnostics, we infer that four of the eleven emission-line galaxies host narrow line active galactic nuclei (AGNs). The AGN host galaxies have stellar populations ranging from evolved to star-forming. Combining our sample with a UV-selected galaxy sample at the same redshift that spans a broader range in stellar mass, we find that black-hole accretion is more effective at the high-mass end of the galaxy distribution (~2.9x10^11 Msun) at z~2.3. Furthermore, by comparing our results with SDSS data, we show that the AGN activity in massive galaxies has decreased significantly between z~2.3 and z~0. AGNs with similar normalized accretion rates as those detected in our K-selected galaxies reside in less massive galaxies (~4.0x10^10 Msun) at low redshift. This is direct evidence for downsizing of AGN host galaxies. Finally, we speculate that the typical stellar mass-scale of the actively accreting AGN host galaxies, both at low and at high redshift, might be similar to the mass-scale at which star-forming galaxies seem to transform into red, passive systems.

This paper presents multiwavelength imaging and broad-band spectroscopy of the relativistic jets in the two nearby radio galaxies 3C 371 and PKS 2201+044, acquired with Chandra, HST, VLA, and Merlin. Radio polarization images are also available. The two sources stand out as "intermediate'' between FRIs and FRIIs; their cores are classified as BL Lacs, although broad and narrow optical emission lines were detected at times. The multiwavelength images show jet morphologies with the X-ray emission peaking closer to the nucleus than the longer wavelengths. The jets are resolved at all wavelengths in a direction perpendicular to the jet axis. The jets SEDs are consistent with a single spectral component from radio to X-rays, interpreted as synchrotron emission. The SEDs show a progressive softening from the inner to the outer regions of the jet, indicating that the electron break energy moves to lower energies with distance from the core. Overall, the X-ray and multiwavelength properties of the jets of 3C 371 and PKS 2201+044 appear intermediate between those of FRIs and FRIIs.

We describe the results of an extremely deep, 0.28 deg^2 survey for z = 3.1 Ly-alpha emission-line galaxies in the Extended Chandra Deep Field South. By using a narrow-band 5000 Anstrom filter and complementary broadband photometry from the MUSYC survey, we identify a statistically complete sample of 162 galaxies with monochromatic fluxes brighter than 1.5 x 10^-17 ergs cm^-2 s^-1 and observers frame equivalent widths greater than 80 Angstroms. We show that the equivalent width distribution of these objects follows an exponential with a rest-frame scale length of w_0 = 76 +/- 10 Angstroms. In addition, we show that in the emission line, the luminosity function of Ly-alpha galaxies has a faint-end power-law slope of alpha = -1.49 +/- 0.4, a bright-end cutoff of log L^* = 42.64 +/- 0.2, and a space density above our detection thresholds of 1.46 +/- 0.12 x 10^-3 h70^3 galaxies Mpc^-3. Finally, by comparing the emission-line and continuum properties of the LAEs, we show that the star-formation rates derived from Ly-alpha are ~3 times lower than those inferred from the rest-frame UV continuum. We use this offset to deduce the existence of a small amount of internal extinction within the host galaxies. This extinction, coupled with the lack of extremely-high equivalent width emitters, argues that these galaxies are not primordial Pop III objects, though they are young and relatively chemically unevolved.

We present Spitzer IRAC imaging of the large-scale jet in the quasar PKS 1136-135 at wavelengths of 3.6 and 5.8 um, combined with previous VLA, HST, and Chandra observations. We clearly detect infrared emission from the jet, resulting in the most detailed multifrequency data among the jets in lobe-dominated quasars. The spectral energy distributions of the jet knots have significant variations along the jet, like the archetypal jet in 3C 273. The infrared measurements with IRAC are consistent with the previous idea that the jet has two spectral components, namely (1) the low-energy synchrotron spectrum extending from radio to infrared, and (2) the high-energy component responsible for the X-ray flux. The optical fluxes may be a mixture of the two components. We consider three radiation models for the high-energy component: inverse Compton scattering off CMB photons by radio-emitting electrons in a highly relativistic jet, synchrotron radiation by a second distinct electron population, and synchrotron radiation by ultra high energy protons. Each hypothesis leads to important insights into and constraints on particle acceleration in the jet, as well as the basic physical properties of the jet such as bulk velocity, transporting power, and particle contents.

We present new Chandra and HST data for four gamma-ray blazars selected on the basis of radio morphology with the aim of revealing X-ray and optical emission from their jets at large scales. All the sources have been detected. Spectral Energy Distributions of the large scale jets are obtained as well as new X-ray spectra for the blazar cores. Modeling for each object the core (sub-pc scale) and large-scale (>100 kpc) jet SEDs, we derive the properties of the same jet at the two scales. The comparison of speeds and powers at different scales supports a simple scenario for the dynamics and propagation of high power relativistic jets.

The Cosmic Evolution Survey (COSMOS) was initiated with an extensive allocation (590 orbits in Cycles 12-13) using the Hubble Space Telescope (HST) for high resolution imaging. Here we review the characteristics of the HST imaging with the Advanced Camera for Surveys (ACS) and parallel observations with NICMOS and WFPC2. A square field (1.8$\sq$\deg) has been imaged with single-orbit ACS I-F814W exposures with 50% completeness for sources 0.5\arcsec in diameter at I$_{AB} $ = 26.0 mag. The ACS imaging is a key part of the COSMOS survey, providing very high sensitivity and high resolution (0.09\arcsec FWHM, 0.05\arcsec pixels) imaging and detecting 1.2 million objects to a limiting magnitude of 26.5 (AB). These images yield resolved morphologies for several hundred thousand galaxies. The small HST PSF also provides greatly enhanced sensitivity for weak lensing investigations of the dark matter distribution.

The Spitzer-COSMOS survey (S-COSMOS) is a Legacy program (Cycles 2+3) designed to carry out a uniform deep survey of the full 2 sq deg COSMOS field in all seven Spitzer bands (3.6, 4.5, 5.6, 8.0, 24.0, 70.0, 160.0 u). This paper describes the survey parameters, mapping strategy, data reduction procedures, achieved sensitivities to date, and the complete data set for future reference. We show that the observed infrared backgrounds in the S-COSMOS field are within 10% of the predicted background levels. The fluctuations in the background at 24u have been measured and do not show any significant contribution from cirrus, as expected. In addition, we report on the number of asteroid detections in the low galactic latitude COSMOS field. We use the Cycle 2 S-COSMOS data to determine preliminary number counts, and compare our results with those from previous Spitzer Legacy surveys (e.g. SWIRE, GOODS). The results from this "first analysis" confirm that the S-COSMOS survey will have sufficient sensitivity with IRAC to detect ~ L* disks and spheroids out to z ~ 3, and with MIPS to detect ultraluminous starbursts and AGN out to z ~3 at 24u and out to z ~1.5-2 at 70u and 160u.

We present deep near-infrared JHK imaging of four 10'x10' fields. The observations were carried out as part of the Multiwavelength Survey by Yale-Chile (MUSYC) with ISPI on the CTIO 4m telescope. The typical point source limiting depths are J~22.5, H~21.5, and K~21 (5sigma; Vega). The effective seeing in the final images is ~1.0". We combine these data with MUSYC UBVRIz imaging to create K-selected catalogs that are unique for their uniform size, depth, filter coverage, and image quality. We investigate the rest-frame optical colors and photometric redshifts of galaxies that are selected using common color selection techniques, including distant red galaxies (DRGs), star-forming and passive BzKs, and the rest-frame UV-selected BM, BX, and Lyman break galaxies (LBGs). These techniques are effective at isolating large samples of high redshift galaxies, but none provide complete or uniform samples across the targeted redshift ranges. The DRG and BM/BX/LBG criteria identify populations of red and blue galaxies, respectively, as they were designed to do. The star-forming BzKs have a very wide redshift distribution, a wide range of colors, and may include galaxies with very low specific star formation rates. In comparison, the passive BzKs are fewer in number, have a different distribution of K magnitudes, and have a somewhat different redshift distribution. By combining these color selection criteria, it appears possible to define a reasonably complete sample of galaxies to our flux limit over specific redshift ranges. However, the redshift dependence of both the completeness and sampled range of rest-frame colors poses an ultimate limit to the usefulness of these techniques.

[ABRIGED] We present the optical identification of a sample of 695 X-ray sources detected in the first 1.3 deg^2 of the XMM-COSMOS survey, down to a 0.5-2 keV (2-10 keV) limiting flux of ~10^-15 erg cm-2 s-1 (~5x10^-15 erg cm^-2 s-1). We were able to associate a candidate optical counterpart to ~90% (626) of the X-ray sources, while for the remaining ~10% of the sources we were not able to provide a unique optical association due to the faintness of the possible optical counterparts (I_AB>25) or to the presence of multiple optical sources within the XMM-Newton error circles. We also cross-correlated the candidate optical counterparts with the Subaru multicolor and ACS catalogs and with the Magellan/IMACS, zCOSMOS and literature spectroscopic data; the spectroscopic sample comprises 248 objects (~40% of the full sample). Our analysis reveals that for ~80% of the counterparts there is a very good agreement between the spectroscopic classification, the morphological parameters as derived from ACS data, and the optical to near infrared colors. About 20% of the sources show an apparent mismatch between the morphological and spectroscopic classifications. All the ``extended'' BL AGN lie at redshift <1.5, while the redshift distribution of the full BL AGN population peaks at z~1.5. Our analysis also suggests that the Type 2/Type 1 ratio decreases towards high luminosities, in qualitative agreement with the results from X-ray spectral analysis and the most recent modeling of the X-ray luminosity function evolution.

We present the first set of XMM-Newton EPIC observations in the 2 square degree COSMOS field. The strength of the COSMOS project is the unprecedented combination of a large solid angle and sensitivity over the whole multiwavelength spectrum. The XMM-Newton observations are very efficient in localizing and identifying active galactic nuclei (AGN) and clusters as well as groups of galaxies. One of the primary goals of the XMM-Newton Cosmos survey is to study the co-evolution of active galactic nuclei as a function of their environment in the Cosmic web. Here we present the log of observations, images and a summary of first research highlights for the first pass of 25 XMM-Newton pointings across the field. In the existing dataset we have detected 1416 new X-ray sources in the 0.5-2, 2-4.5 and 4.5-10 keV bands to an equivalent 0.5-2 keV flux limit of 7x10-16 erg cm-2 s-1. The number of sources is expected to grow to almost 2000 in the final coverage of the survey. From an X-ray color color analysis we identify a population of heavily obscured, partially leaky or reflecting absorbers, most of which are likely to be nearby, Compton-thick AGN.

The VLA-COSMOS large project is described and its scientific objective is discussed. We present a catalog of ~ 3,600 radio sources found in the 2deg^2 COSMOS field at 1.4 GHz. The observations in the VLA A and C configuration resulted in a resolution of 1.5''x1.4'' and a mean rms noise of ~ 10.5(15) uJy/beam in the central 1(2)deg^2. 80 radio sources are clearly extended consisting of multiple components, and most of them appear to be double-lobed radio galaxies. The astrometry of the catalog has been thoroughly tested and the uncertainty in the relative and absolute astrometry are 130mas and <55mas, respectively.

In order to study the evolution of the relative fraction of obscured Active Galactic Nuclei (AGN) we constructed the largest sample to date of AGN selected in hard X-rays. The full sample contains 2341 X-ray-selected AGN, roughly 4 times the largest previous samples studied in this connection. Of these, 1229 (53%) have optical counterparts for which redshifts are available; these span the redshift range $z$=0-4. The observed fraction of obscured AGN declines only slightly with redshift. Correcting for selection bias, we find that the intrinsic fraction of obscured AGN must actually increase with redshift, as (1+$z$)$^{\alpha}$, with $\alpha$$\simeq0.4\pm$0.1. This evolution is consistent with the integrated X-ray background, which provides the strongest constraints at relatively low redshift, $z\sim$1. Summing over all AGN, we estimate the bolometric AGN light to be 3.8 nW m$^{-2}$ sr$^{-1}$, or $\lesssim$8% of the total extragalactic light. Together with the observed black hole mass density in the local Universe, this implies an accretion efficiency of $\eta \sim$0.1-0.2, consistent with the values typically assumed.

We have performed infrared imaging of the jet of the quasar 3C 273 at wavelengths 3.6 and 5.8 microns with the Infrared Array Camera (IRAC) on the Spitzer Space Telescope. When combined with the radio, optical and X-ray measurements, the IRAC photometry clearly shows that the optical emission is dominated by the high-energy component of the jet, not by the radio synchrotron component, as had been assumed to date. The high-energy component may be due to a second synchrotron component or to IC scattering of ambient photons. In the former case, we argue that the acceleration of protons exceeding 10^16 eV or possibly even to 10^19 eV would be taking place in the jet. In contrast, the IC model, into which highly relativistic Doppler beaming has to be incorporated, requires very low-energy electrons (~ 1 MeV). The present polarization data in the radio and optical would favor the former interpretation in the case of the 3C 273 jet. Sensitive and detailed measurements of optical polarization are important to establish the radiation mechanism responsible for the high-energy emission. The present study offers new clues as to the controversial origin of the X-ray emission seen in many quasar jets.

We selected 40 candidate Lyman Alpha Emitting galaxies (LAEs) at z ~=3.1 with observed frame equivalent widths >150A and inferred emission line fluxes >2.5x10^-17 ergs/cm^2/s from deep narrow-band and broad-band MUSYC images of the Extended Chandra Deep Field South. Covering 992 sq. arcmin, this is the largest ``blank field'' surveyed for LAEs at z ~3, allowing an improved estimate of the space density of this population of 3+-1x10^-4 h_70^3/Mpc^3. Spectroscopic follow-up of 23 candidates yielded 18 redshifts, all at z ~=3.1. Over 80% of the LAEs are dimmer in continuum magnitude than the typical Lyman break galaxy spectroscopic limit of R= 25.5 (AB), with a median continuum magnitude R ~=27 and very blue continuum colors, (V-z) ~=0. Over 80% of the LAEs have the right UVR colors to be selected as Lyman break galaxies, but only 10% also have R<=25.5. Stacking the UBVRIzJK fluxes reveals that LAEs have stellar masses ~=5x10^8 h_70^-2 M_sun and minimal dust extinction, A_V < ~ 0.1. Inferred star formation rates are ~=6 h_70^-2 M_sun/yr, yielding a cosmic star formation rate density of 2x10^-3 h_70 M_sun/yr/Mpc^3. None of our LAE candidates show evidence for rest-frame emission line equivalent widths EW_rest>240A which might imply a non-standard IMF. One candidate is detected by Chandra, implying an AGN fraction of 2+-2% for LAE candidate samples. In summary, LAEs at z ~ 3 have rapid star formation, low stellar mass, little dust obscuration and no evidence for a substantial AGN component.

The Extended Chandra Deep Field-South (ECDFS) survey consists of 4 Chandra ACIS-I pointings and covers $\approx$ 1100 square arcminutes ($\approx$ 0.3 deg$^2$) centered on the original CDF-S field to a depth of approximately 228 ks. This is the largest Chandra survey ever conducted at such depth, and only one XMM-Newton survey reaches a lower flux limit in the hard 2.0--8.0 keV band. We detect 651 unique sources -- 587 using a conservative source detection threshold and 64 using a lower source detection threshold. These are presented as two separate catalogs. Of the 651 total sources, 561 are detected in the full 0.5--8.0 keV band, 529 in the soft 0.5--2.0 keV band, and 335 in the hard 2.0--8.0 keV band. For point sources near the aim point, the limiting fluxes are approximately $1.7 \times 10^{-16}$ $\rm{erg cm^{-2} s^{-1}}$ and $3.9 \times 10^{-16}$ $\rm{erg cm^{-2} s^{-1}}$ in the 0.5--2.0 keV and 2.0--8.0 keV bands, respectively. Using simulations, we determine the catalog completeness as a function of flux and assess uncertainties in the derived fluxes due to incomplete spectral information. We present the differential and cumulative flux distributions, which are in good agreement with the number counts from previous deep X-ray surveys and with the predictions from an AGN population synthesis model that can explain the X-ray background. In general, fainter sources have harder X-ray spectra, consistent with the hypothesis that these sources are mainly obscured AGN.

Modeling the multiwavelength emission of successive regions in the jet of the quasar PKS 1136-135 we find indication that the jet suffers deceleration near its end on a (deprojected) scale of ~400 kpc. Adopting a continuous flow approximation we discuss the possibility that the inferred deceleration from a Lorentz factor Gamma=6.5 to Gamma=2.5 is induced by entrainment of external gas. Some consequences of this scenario are discussed.

We present mid-infrared observations of AGN in the GOODS fields, performed with the Spitzer Space Telescope. These are the deepest infrared and X-ray fields to date and cover a total area of ~0.1 square degrees. AGN are selected on the basis of their hard (2-8 keV) X-ray emission. The median AGN infrared luminosity is at least 10 times larger than the median for normal galaxies with the same redshift distribution, suggesting that the infrared emission is dominated by the central nucleus. The X-ray to infrared luminosity ratios of GOODS AGN, most of which are at 0.5<z<1.5, are similar to the values obtained for AGN in the local Universe. The observed infrared flux distribution has an integral slope of ~1.5 and there are 1000 sources per square degree brighter than ~50 uJy at 3-6 microns. The counts approximately match the predictions of models based on AGN unification, in which the majority of AGN are obscured. This agreement confirms that the faintest X-ray sources, which are dominated by the host galaxy light in the optical, are obscured AGN. Using these Spitzer data, the AGN contribution to the extragalactic infrared background light is calculated by correlating the X-ray and infrared catalogues. This is likely to be a lower limit given that the most obscured AGN are missed in X-rays. We estimate the contribution of AGN missed in X-rays, using a population synthesis model, to be ~45% of the observed AGN contribution, making the AGN contribution to the infrared background at most ~2-10% in the 3-24 micron range, depending on wavelength, lower than most previous estimates. The AGN contribution to the infrared background remains roughly constant with source flux in the IRAC bands but decreases with decreasing flux in the MIPS 24 um band, where the galaxy population becomes more important.

The Great Observatories Origins Deep Survey (GOODS) combines deep HST and Spitzer imaging with the deepest Chandra/XMM observations to probe obscured AGN at higher redshifts than previous multiwavelength surveys. We present a self-consistent implementation of the AGN unification paradigm, which postulates obscured AGN wherever there are unobscured AGN, to successfully explain the infrared, optical, and X-ray number counts of X-ray sources detected in the GOODS fields. Assuming either a constant ratio of obscured to unobscured AGN of 3:1 (the local value), or a ratio that decreases with luminosity, and including Compton-thick sources, we can explain the spectral shape and normalization of the extragalactic X-ray "background" as a superposition of unresolved AGN, predominantly at z~0.5-1.5 and L_x~10^43-10^44 ergs/s. The possible dependence of the obscured to unobscured ratio with redshift is not well constrained; present data allow it to decrease or increase substantially beyond z~1.

The Extended Chandra Deep Field-South (E-CDF-S) survey consists of 4 Chandra ACIS-I pointings and covers ~1100 square arcminutes (~0.3 deg^2) centered on the original CDF-S field to a depth of approximately 228 ks (PI: Niel Brandt; Lehmer et al, 2005). This is the largest Chandra survey ever conducted at such depth. In our analysis (Virani et al, 2005), we detect 651 unique sources -- 587 using a conservative source detection limit and 64 using a lower source detection limit. Of these 651 sources, 561 are detected in the full 0.5--8.0 keV band, 529 in the soft 0.5--2.0 keV band, and 335 in the hard 2.0--8.0 keV band. For point sources near the aim point, the limiting fluxes are approximately $1.7 \times 10^{-16}$ $\rm{erg cm^{-2} s^{-1}}$ and $3.9 \times 10^{-16}$ $\rm{erg cm^{-2} s^{-1}}$ in the 0.5--2.0 keV and 2.0--8.0 keV bands, respectively. We present the optical properties of these X-ray sources, specifically the R-band magnitude distribution and a preliminary spectroscopic redshift distribution. One exciting result is the discovery of 7 new Extreme X-ray-to-Optical flux ratio objects (EXOs) found in the E-CDF-S field.

We present UBVRIz' optical images of the 0.32 square degree Extended Hubble Deep Field South. This is one of four fields comprising the MUSYC survey, which is optimized for the study of galaxies at z=3, AGN demographics, and Galactic structure. We calculate corrected aperture photometry and its uncertainties and find through tests that these provide a significant improvement upon standard techniques. Our photometric catalog of 62968 objects is complete to a total magnitude of R_AB=25. We select z=3 Lyman break galaxy (LBG) candidates from their UVR colors and find a sky surface density of 1.4/sq. arcmin and an angular correlation function w(theta) = 2.3+-1.0 theta^{-0.8}. (Abridged)

Emission from the relativistic jet located at hundreds of kpc from the core of the superluminal quasar PKS 0637-752 was detected at 3.6 and 5.8 microns with the Infrared Array Camera (IRAC) on the Spitzer Space Telescope. The unprecedented sensitivity and arcsecond resolution of IRAC allows us to explore the mid-infrared emission from kiloparsec-scale quasar jets for the first time. The mid-infrared flux from the jet knots, when combined with radio and optical fluxes, confirms a synchrotron origin of the radio-to-optical emission and constrains very well the high energy end of the nonthermal electron distribution. Assuming the X-rays are produced in the relativistically moving knots via inverse Compton scattering of cosmic microwave background (CMB) radiation, the infrared observation puts constraints on the matter content of the quasar extended jet. Specifically, pure electron-positoron pair jet models are unfavorable based on the lack of an infrared bump associated with ``bulk Comptonization'' of CMB photons by an ultrarelativistic jet.

We investigate morphological structure parameters and local environments of distant moderate-luminosity active galactic nucleus (AGN) host galaxies in the overlap between the HST/ACS observations of the Great Observatories Origins Deep Survey (GOODS) and the two Chandra Deep Fields. We compute near-neighbor counts and BViz asymmetry (A) and concentration (C) indices for ~35,500 GOODS/ACS galaxies complete to z_850 ~ 26.6, including the resolved hosts of 322 X-ray-selected AGNs. Distributions of (1) z_850 asymmetry for 130 AGN hosts at z_850 < 23 and (2) near-neighbor counts for 173 AGN hosts at z_850 < 24 are both consistent with non-AGN control samples. This implies no close connection between recent galaxy mergers and moderate-luminosity AGN activity out to appreciable look-back times (z < 1.3), approaching the epoch of peak AGN activity in the universe. The distribution of z_850 concentrations for the AGN hosts is offset by +0.5 compared to the non-AGN, a 6.4-sigma discrepancy much larger than can be explained by the possible influence of unresolved emission from the AGN or a circumnuclear starburst. The local universe association between AGN and bulge-dominated galaxies thus persists to substantial look-back time. We discuss implications in the context of the low-redshift supermassive central black hole mass correlation with host galaxy properties, including concentration.

We report stellar velocity dispersion measurements for a sample of 28 AGN host galaxies including our previous work. Using the mass-dispersion ($M_{\bullet}-\sigma$) and the fundamental plane relations, we estimate the black hole mass for a sample of 66 BL Lac objects and investigate the role of black hole mass in the energetics of BL Lac objects. The black hole mass range for different BL Lac spectral types is similar, $10^{7} < M_{\bullet} < 4 \times 10^{9}$. Neither X-ray nor radio luminosity correlates with black hole mass. Low-frequency-peaked BL Lac objects have higher Eddington ratios on average, because of either more beaming or higher intrinsic power. For the black hole mass range $3 \times 10^{7} < M_{\bullet} < 10^{9}$, the radio luminosity of BL Lac objects and flat-spectrum radio quasars spans over 4 orders of magnitude, with BL Lac objects being low-power AGNs. We also investigate the evolution of host galaxies for 39 AGNs out to $z \approx 0.5$ with measuredstellar velocity dispersions. Comparing the mass-to-light ratio evolution in the observed frame with population synthesis models, we find that single burst star formation models with $z_{form} = 1.4^{+0.9}_{-0.2} $ are consistent with the observations. From our $z_{form}=1.4$ model, we estimated the intrinsic mass-to-light ratio evolution in the Cousins $R$ band, $\Delta log (M/L)/ \Delta z = -0.502 \pm 0.08$, consistent with that of normal early type galaxies.

Active Galactic Nuclei (AGN) unification, which implies a large number of obscured AGN, can explain the optical, infrared and X-ray content of deep multiwavelength surveys. Here, we show that the same model also successfully explains the spectral shape and intensity of the X-ray background. The simplest possible unified model assumes a constant ratio of obscured to unobscured AGN, independent of redshift or luminosity. With almost no free parameters, the predicted X-ray background agrees remarkably well with observations. Both the observed properties of deep AGN samples and the X-ray background can be explained by the same model without invoking different evolution for obscured and unobscured AGN. Recent observational evidence shows that the ratio of obscured to unobscured AGN may depend on luminosity. An AGN unification model that allows for such a dependence also fits the X-ray background spectrum, and the predicted numbers counts in the 2-10 keV X-ray band are in good agreement with observations. We present predictions for the source counts in hard X-rays (10-100 keV), which indicate that future missions like the Black Hole Finder will observe thousands of heavily obscured AGN.

We present analysis of spectral energy distributions (SEDs) from mid-infrared through X-ray of a sample of 420 hard X-ray selected, z-band and Spitzer/IRAC detected active galactic nuclei (AGN) and AGN candidates from the GOODS multiwavelength survey. We fit local empirical SED templates of both normal and active galaxies to the rest-frame luminosities calculated from spectroscopic (where available) and photometric redshifts. The majority of the optically stellar-dominated (with early-type galaxy fits) sources are moderately luminous (L_X,2-10 keV ~ 10^43 erg/s) hard X-ray sources with high X-ray hardness ratios (HR > 0.2), high MIR luminosities and red MIR colors in excess of a typical stellar dominated elliptical galaxy. These sources likely harbor heavily obscured (though Compton-thin) AGN. The observed ratio of obscured to unobscured AGN has an integrated mean of ~3.4:1 but declines with increasing redshift. This effect has been explained by Treister et al. (2004) as an observational bias triggered by the lack of spectroscopic redshifts at R > 24 which are predominantly higher redshift obscured sources.

Using the known numbers and evolution of unobscured AGN, and assuming that there is a large number of corresponding obscured AGN -- as expected if AGN unification is correct -- we predict the infrared number counts of AGN. These agree well with the observed GOODS Spitzer data at 3.6, 8 and 24 $\mu$m, provided that the dust content of the obscuring torus is $\sim 2$x less than assumed. This strongly supports the AGN unification picture, and implies a large population of infrared bright obscured AGN at all redshifts and luminosities.

Using a simple unification model for AGN, with a fixed ratio of obscured to unobscured AGN of 3:1, we explain the X-ray, optical and infrared properties of the X-ray sources in the GOODS fields. That is, the GOODS data are consistent with the existence of a large population of obscured AGN out to high redshifts. About half of these are so obscured that they are missed even by hard X-ray observations, and are detected only in deep infrared observations. The previously reported decreasing trend of obscured to total AGN ratio with increasing X-ray luminosity can be explained entirely as a selection effect.

We present the first results from the Cal\'an-Yale Deep Extragalactic Research (CYDER) survey. The main goal of this survey is to study serendipitous X-ray sources detected by Chandra in an intermediate flux range ($10^{-15}-10^{-12}$ ergs s$^{-1}$) that comprises most of the X-ray background. 267 X-ray sources spread over 5 archived fields were detected. The $\log N-\log S$ distribution obtained for this sample is consistent with the results of other surveys. Deep $V$ and $I$ images were taken of these fields in order to calculate X-ray-to-optical flux ratios. Identifications and redshifts were obtained for 106 sources using optical spectroscopy from 8-m class telescopes to reach the optically faintest sources, to the same level as deeper X-ray fields like the Chandra Deep Fields, showing that the nature of sources detected depends mostly on the optical limit for spectroscopy. In general, sources optically classified as obscured Active Galactic Nuclei (AGNs) have redder optical colors than unobscured AGN. A rough correlation between $f_X/f_{\rm opt}$ and hard X-ray luminosity was found for obscured AGN confirming the prediction by existing models that in obscured AGN the optical light is completely dominated by the host galaxy. The previously claimed decrease of the obscured to unobscured AGN ratio with increasing X-ray luminosity is observed. However, this correlation can be explained as a selection effect caused by the lower optical flux of obscured AGN. Comparison between the observed $N_H$ distribution and predictions by existing models shows that the sample appears complete up to $N_H<3\times 10^{22}$ cm$^{-2}$, while for more obscured sources incompleteness plays an important role in the observed obscured to unobscured AGN ratio.

We investigate the luminosity evolution of the host galaxies of radio-loud AGN through Hubble Space Telescope imaging of 72 BL Lac objects, including new STIS imaging of nine z > 0.6 BL Lacs. With their intrinsically low accretion rates and their strongly beamed jets, BL Lacs provide a unique opportunity to probe host galaxy evolution independent of the biases and ambiguities implicit in quasar studies. We find that the host galaxies of BL Lacs evolve strongly, consistent with passive evolution from a period of active star formation in the range 0.5 <~ z <~ 2.5, and inconsistent with either passive evolution from a high formation redshift or a non-evolving population. This evolution is broadly consistent with that observed in the hosts of other radio-loud AGN, and inconsistent with the flatter luminosity evolution of quiescent early types and radio-quiet hosts. This indicates that active star formation, and hence galaxy interactions, are associated with the formation for radio-loud AGN, and that these host galaxies preferentially accrete less material after their formation epoch than galaxies without powerful radio jets. We discuss possible explanations for the link between merger history and the incidence of a radio jet.

We analysed the optical and infrared properties of X-ray sources in the Great Observatories Origins Deep Survey (GOODS), a deep multiwavelength survey covering 0.1 square degrees in two fields. The HST ACS data are well explained by a unified AGN scheme that postulates roughly 3 times as many obscured as unobscured AGN, as are the spectroscopic and photometric redshift distributions once selection effects are considered. Our model predicts infrared number counts of AGN that agree well with the preliminary Spitzer data, confirming that large numbers of obscured AGN are present in the early Universe (z>1).

H1517+656 is an unusual source, even for a BL Lac object. It is one of the most luminous BL Lacs known, with extreme emission properties at radio, optical, and X-ray frequencies. Furthermore, in our WFPC2 snapshot survey we discovered a series of peculiar arcs describing a 2.4 arcsec radius ring surrounding the source. This paper describes follow-up observations with additional WFPC2 bands and the STIS longpass filter, which have revealed this structure to be the remnants of a very recent galaxy merger. Population synthesis modelling has shown that regions of the arcs have stellar populations with age < 20 Myrs. Additionally, the circularity of the arcs indicates that the plane of the collision and hence accretion is very close to the plane of the sky. Given that BL Lac jets are thought to be aligned with the line of sight, this observation may provide a direct link between the transfer of angular momentum in an interaction and the generation of a radio source.

The deep X-ray, optical, and far-infrared fields that constitute GOODS are sensitive to obscured AGN (N_H>10^{22} cm^{-2}) at the quasar epoch (z~2-3), as well as to unobscured AGN as distant as z~7. Luminous X-ray emission is a sign of accretion onto a supermassive black hole and thus reveals all but the most heavily obscured AGN. We combine X-ray luminosity functions with appropriate spectral energy distributions for AGN to model the X-ray, optical and far-infrared flux distributions of the X-ray sources in the GOODS fields. A simple model based on the unified paradigm for AGN, with ~3 times as many obscured AGN as unobscured, successfully reproduces the z-band flux distributions measured in the deep HST ACS observations on the GOODS North and South fields. This model is also consistent with the observed spectroscopic and photometric redshift distributions once selection effects are considered. The previously reported discrepancy between observed spectroscopic redshift distributions and the predictions of population synthesis models for the X-ray background can be explained by bias against the most heavily obscured AGN generated both by X-ray observations and the identification of sources via optical spectroscopy. We predict the AGN number counts for Spitzer MIPS 24 um and IRAC 3.6-8 um observations in the GOODS fields, which will verify whether most AGN in the early Universe are obscured in the optical. Such AGN should be very bright far-infrared sources and include some obscured AGN missed even by X-ray observations.

We report on deep radio images and optical spectroscopy of two BL Lac objects that have very close compact companions. The two targets, 1ES 0502+675 and 1ES 1440+122, were selected from the HST imaging survey of 110 BL Lacs as candidates for possible gravitational lensing. The new observations clearly demonstrate that the companion objects are not secondary images of the active nuclei but, in spite of the relatively low chance projection probability, foreground Galactic stars. Gravitational lensing appears to be unimportant to the BL Lac phenomenon. We discuss the radio properties of the BL Lac objects in the context of standard beaming models, and show they are as expected for beamed FRI radio galaxies.

The Chandra discovery of bright X-ray emission from kpc-scale jets allows insight into the physical parameters of the jet flow at large scale. At the opposite extreme, extensive studies of the inner relativistic jets in Blazars with multiwavelength observations, yield comparable information on sub-parsec scales. In the framework of simple radiation models for the emission regions we compare the physical parameters of jets on these two very different scales in the only two well studied Blazars for which large-scale emission has been resolved by Chandra. Notably, we find that the relativistic Doppler factors and powers derived independently at the two scales are consistent, suggesting that the jet does not suffer severe deceleration or dissipation. Moreover the internal equipartition pressures in the inner jet and in the external X-ray bright knots scale inversely with the jet cross section as expected in the simple picture of a freely expanding jet in equipartition.

The Great Observatories Origins Deep Survey provides significant constraints on the space density of less luminous QSOs at high redshift, which is particularly important to understand the interplay between the formation of galaxies and super-massive black holes and to measure the QSO contribution to the UV ionizing background. We present the results of a search for high-z QSOs, identified in the two GOODS fields on the basis of deep imaging in the optical (with HST) and X-ray (Chandra), and discuss the allowed space density of QSOs in the early universe.

We present the results from an X-ray and optical survey of a sample of 17 radio jets in AGN performed with Chandra and HST. The sample was selected from the radio and is unbiased toward detection at shorter wavelengths, but preferentially it includes beamed sources. We find that X-ray emission is common on kpc-scales, with over half radio jets exhibiting at least one X-ray knot on the Chandra images. The distributions of the radio-to-X-ray and radio-to-optical spectral indices for the detected jets are similar to the limits for the non-detections,suggesting all bright radio jets have X-ray counterparts which will be visible in longer observations. Comparing the radio and X-ray morphologies shows that the majority of the X-ray jets have structures that closely map the radio. Analysis of the SED of the jet knots suggest the knots in which the X-ray and radio morphologies track each other produce X-rays by IC scattering of the Cosmic Microwave Background. The remaining knots produce X-rays by the synchrotron process. Spectral changes are detected along the jets, with the ratio of the X-ray-to-radio and optical-to-radio flux densities decreasing from the inner to the outer regions. This suggests the presence of an additional contribution to the X-ray flux in the jet's inner part, either from synchrotron or IC of the stellar light. Alternatively, in a pure IC/CMB scenario, the plasma decelerates as it flows from the inner to the outer regions. Finally, the X-ray spectral indices for the brightest knots are flat, indicating that the bulk of the luminosity of the jets is emitted at GeV energies, and raising the interesting possibility of future detections with GLAST.

The paradigm for AGN unification is reviewed, in terms of its optical manifestation as obscuration in the equatorial plane and its radio manifestation as relativistic beaming of the jet emission. Within this paradigm, observed AGN properties depend strongly on orientation angle. The predictions of unification are commensurate with the local numbers of AGN of various types. The outstanding question concerns a possible population of obscured AGN at high redshift, which are thought to produce the X-ray background, although few are observed directly. We describe early results from GOODS, an HST-SIRTF-Chandra multiwavelength survey, which has as one of its principal goals the search for such obscured AGN.

We present a sample of 17 high-redshift (3.5<z<5.2) QSO candidates in the 320 sq.arcmin area of the Great Observatories Origins Deep Survey, selected in the magnitude range 22.45<z_{850}<25.25 using deep imaging with the Advanced Camera for Surveys onboard the Hubble Space Telescope and the Advanced CCD Imaging Spectrometer onboard the Chandra X-ray Observatory. On the basis of seven spectroscopic and ten photometric redshifts we estimate that the final sample will contain between two and four QSOs with 4<z<5.2. A dearth of high-redshift, moderate-luminosity (M_{145}=~-23) QSOs is observed with respect to predictions based on a) the extrapolation of the z~2.7 luminosity function (LF), according to a pure luminosity evolution calibrated by the results of the Sloan Digital Sky Survey; and b) a constant universal efficiency in the formation of super-massive black holes (SMBHs) in dark-matter halos. Evidence is gathered in favor of a density evolution of the LF at high redshift and of a suppression of the formation or feeding of SMBHs in low-mass halos.

Multi-band near-infrared images of twelve BL Lacertae objects were obtained with the 2.5m telescope at the Las Campanas Observatory in order to determine the properties of their underlying host galaxies. Resolved emission was clearly detected in eight of the lowest redshift targets (up to z~0.3), and was modeled with a de Vaucouleurs r^{1/4} surface brightness law. We find that the morphologies match the elliptical galaxy profiles well, and that the BL Lac objects reside in large and luminous, but otherwise normal hosts -- consistent with previous studies done predominantly at optical wavelengths. The median absolute K-band magnitude of the galaxies in this study is -26.2, the average half-light radius is 4.2 +/- 2.3 kpc, and their average integrated R-K color is 2.7 +/- 0.3 mag. These are well within the range of values measured previously in the H-band by Kotilainen et al. and Scarpa et al. in a comparable number of targets. Taking their data together with our results, we find a best-fit K-band Kormendy relation of \mu_{e} = 4.3 log_{10} (r_{e}/kpc) + 14.2 mag arcsec^{-2}, virtually identical to that obtained for normal ellipticals. Finally, the near-infrared colors determined for five galaxies (average J-K = 0.8 +/- 0.3 mag) are the first such measurements for BL Lac hosts, and match those expected from old stellar populations at the BL Lac redshifts.

We describe a possible new class of X-ray sources that have robust detections in ultra-deep Chandra data, yet have no detections at all in our deep multi-band GOODS Hubble Space Telescope (HST) ACS images, which represent the highest quality optical imaging obtained to date on these fields. These extreme X-ray / Optical ratio sources ("EXO"s) have values of Fx/Fopt at least an order of magnitude above those generally found for other AGN, even those that are harbored by reddened hosts. We thus infer two possible scenarios: (1) if these sources lie at redshifts z < 6, then their hosts need to be exceedingly underluminous, or more reddened, compared with other known sources; (2) if these sources lie above z ~ 6-7, such that even their Lyman-alpha emission is redshifted out of the bandpass of our ACS z(850) filter, then their optical and X-ray fluxes can be accounted for in terms of relatively normal L* hosts and moderate-luminosity AGN.

The radio galaxy Centaurus A was observed by the BeppoSAX satellite five times from 1997 to 2000. From July 6 1999 to August 17 1999, the source was also simultaneously pointed by COMPTEL on-board of the satellite CGRO. Centaurus A has a complex spectrum with multiple extended components and a strongly absorbed (N_H\sim 10^{23}$ cm${-2}) nucleus well fitted by a power law (\Gamma \sim 1.8) which bends at high energies. When the BeppoSAX and COMPTEL observations are combined together, an exponential cutoff with e-folding energy \sim 1000 keV gives an adequate description of the spectral steepening. A complex feature in emission at 6-7 keV is resolved into two Fe components, one narrow cold line and an ionized line centred at 6.8 keV. Significant variations have been observed in the iron feature, with the less prominent ionized line seemingly being the only one responsible for them: its variations do not appear to correlate with the strength of the continuum. The high energy cutoff and the Fe feature suggest the presence of an accretion flow in the Centaurus A nucleus. However the absence of a significant reflection, the narrowness of the cold line as well as the lack of correlation between the continuum and 6.8 keV line variations disfavour a standard cold/ionized thin disk (at least in the inner regions). A more plausible configuration might be a hot thick optically thin accretion flow surrounded by material with different opacities. Finally, we note that high energy break observed by BeppoSAX and COMPTEL could be also reasonably explained by Inverse Compton radiation from a jet. If this is the case, a structured jet with outer slow layers surrounding a beamed inner region is necessary to explain the strong Fe feature observed by BeppoSAX.

We present the X-ray properties of a sample of 17 radio sources observed with the Chandra X-ray Observatory as part of a project aimed at studying the X-ray emission from their radio jets. In this paper, we concentrate on the X-ray properties of the unresolved cores. The sample includes 16 quasars (11 core-dominated and 5 lobe-dominated) in the redshift range z=0.30--1.96, and one low-power radio-galaxy at z=0.064. No diffuse X-ray emission is present around the cores of the quasars, except for the nearby low-power galaxy that has diffuse emission on a scale and with a luminosity consistent with other FRIs. No high-amplitude, short-term variability is detected within the relatively short Chandra exposures. However, 1510-089 shows low-amplitude flux changes with a timescale of $\sim$25 minutes. The X-ray spectra of the quasar cores are generally well described by a single power law model with Galactic absorption. However, in six quasars we find soft X-ray excess emission below 1.6 keV. Interestingly, we detect an Fe K-shell emission line, consistent with fluorescent Kalpha emission from cold Iron, in one lobe- and two core-dominated sources. The average X-ray photon index for the quasars in the sample is Gamma=1.66 and dispersion, sigma=0.23. The average spectral slope for our sample is flatter than the slope found for radio-quiet quasars and for radio-loud AGNs with larger jet orientations; this indicates that beaming affects the X-ray emission from the cores in our sample of quasars.

We have calculated inverse-Compton Doppler factors for 3C 279 using the collection of VLBI data recently published by us, and the collection of multiwavelength spectra recently published by Hartman et al. From the Doppler factor and superluminal apparent speed, we then calculate the Lorentz factor and angle to the line-of-sight of the parsec-scale relativistic jet. We model the jet components as homogeneous spheres and the VLBI core as an unresolved inhomogeneous conical jet. The conical-jet model can be made to match both the observed X-ray emission and the VLBI properties of the core with a suitable choice of Doppler factor, implying the core makes a significant contribution to the X-ray emission. The parameters of the conical models indicate the jet is particle dominated at the radii that produce significant emission, and is not in equipartition. At the inner radius of the conical jet the magnetic field is of order 0.1 G and the relativistic-particle number density is of order 10 cm^{-3}. When all components are included in the calculation, then on average the core produces about half of the X-rays, with the other half being split between the long-lived component C4 and the brightest inner-jet component. We calculate an average speed and angle to the line-of-sight for the region of the jet interior to 1 mas of v=0.992c (gamma=8) and 4 degrees, and an average speed and angle to the line-of-sight for C4 (at a distance from the core of 3 mas) of v=0.997c (gamma=13) and 2 degrees. These values imply average Doppler factors of delta=12 for the inner jet, and delta=21 for C4.

We revisit the issue of whether radio loudness in AGN is associated with central black hole mass, as has been suggested in the literature. We present new estimates of black hole mass for 295 AGN (mostly radio-quiet), calculating their radio loudnesses from published radio and optical fluxes, and combine with our previously published values, for a sample of 452 AGN for which both black hole mass and radio loudness (or upper limits thereto) are known. Among the radio-quiet AGN, there are now many black holes with mass larger than $10^{9}$ $M_{\odot}$, extending to the same high black hole masses as radio-loud AGN of the same redshifts. Over the full sample, the black hole masses of radio-loud and radio-quiet AGN span the same large range, $10^6 - 10^{10} M_\odot$. We conclude that radio loudness in AGN does not depend strongly on central black hole mass.

We have studied the host galaxies of a sample of radio-loud AGN spanning more than four decades in the energy output of the nucleus. The core sample includes 40 low-power sources (BL Lac objects) and 22 high-power sources (radio-loud quasars) spanning the redshift range z~0.15 to z~0.5, all imaged with the high spatial resolution of HST. All of the sources are found to lie in luminous elliptical galaxies, which follow the Kormendy relation for normal ellipticals. A very shallow trend is detected between nuclear brightness (corrected for beaming) and host galaxy luminosity. Black hole masses are estimated for the entire sample, using both the bulge luminosity--black hole mass and the velocity dispersion--black hole mass relations for local galaxies. The latter involves a new method, using the host galaxy morphological parameters, mu_e and r_e, to infer the velocity dispersion, sigma, via the fundamental plane correlation. Both methods indicate that the entire sample of radio-loud AGN are powered by very massive central black holes, with M_{black hole} ~ 10^8 to 10^10 M_{sun}$. Eddington ratios range from L/L_{Eddington} ~ 2 x 10^-4 to ~1, with the high- power sources having higher Eddington ratios than the low-power sources. Overall, radio-loud AGN appear to span a very large range in accretion efficiency, which is all but independent of the mass of the host galaxy.

PKS 0537-441, a bright gamma-ray emitting blazar, was observed at radio, optical, UV and X-ray frequencies during various EGRET pointings, often quasi-simultaneously. In 1995 the object was found in an intense emission state at all wavelengths. BeppoSAX observations made in 1998, non-simultaneously with exposures at other frequencies, allow us to characterize precisely the spectral shape of the high energy blazar component, which we attribute to inverse Compton scattering. The optical-to-gamma-ray spectral energy distributions at the different epochs show that the gamma-ray luminosity dominates the bolometric output. This, together with the presence of optical and UV line emission, suggests that, besides the synchrotron self-Compton mechanism, the Compton upscattering of photons external to the jet (e.g., in the broad line region) may have a significant role for high energy radiation. The multiwavelength variability can be reproduced by changes of the plasma bulk Lorentz factor. The spectrum secured by IUE in 1995 appears to be partially absorbed shortward of ~1700 A. However, this signature is not detected in the HST spectrum taken during a lower state of the source. The presence of intervening absorbers is not supported by optical imaging and spectroscopy of the field.

Black hole mass, along with mass accretion rate, is a fundamental property of active galactic nuclei. Black hole mass sets an approximate upper limit to AGN energetics via the Eddington limit. We collect and compare all AGN black hole mass estimates from the literature; these 177 masses are mostly based on the virial assumption for the broad emission lines, with the broad-line region size determined from either reverberation mapping or optical luminosity. We introduce 200 additional black hole mass estimates based on properties of the host galaxy bulges, using either the observed stellar velocity dispersion or using the fundamental plane relation to infer $\sigma$; these methods assume that AGN hosts are normal galaxies. We compare 36 cases for which black hole mass has been generated by different methods and find, for individual objects, a scatter as high as a couple of orders of magnitude. The less direct the method, the larger the discrepancy with other estimates, probably due to the large scatter in the underlying correlations assumed. Using published fluxes, we calculate bolometric luminosities for 234 AGNs and investigate the relation between black hole mass and luminosity. In contrast to other studies, we find no significant correlation of black hole mass with luminosity, other than those induced by circular reasoning in the estimation of black hole mass. The Eddington limit defines an approximate upper envelope to the distribution of luminosities, but the lower envelope depends entirely on the sample of AGN included. For any given black hole mass, there is a range in Eddington ratio of up to three orders of magnitude.

XMM observed the BL Lac PKS 2155-304 for a full orbit (about 150 ksec) on 2000 November 19-21. Preliminary results on the temporal and spectral analysis of data from the EPIC PN camera and Optical Monitor are presented. The variability amplitude depends systematically on energy, however the slopes of the structure functions of the light-curves in different bands do not appear to be significantly different. No evidence of time lags is found by cross correlating the light-curves in different bands.

We present the first results from an X-ray and optical survey of a sample of AGN radio jets with Chandra and HST. We focus here on the first six sources observed at X-rays, in four of which a bright X-ray jet was detected for the first time. In three out of four cases optical emission from the jet is also detected in our HST images. We compare the X-ray morphology with the radio as derived from improved processing of archival VLA data and we construct spectral energy distributions (SED) for the most conspicuous emission knots. In most cases the SEDs, together with the similarity of the X-ray and radio morphologies, favor an inverse Compton origin of the X-rays. The most likely origin of the seed photons is the Cosmic Microwave Background, implying the jets are still relativistic on kiloparsec scales. However, in the first knot of the PKS 1136-135 jet, X-rays are likely produced via the synchrotron process. In all four cases bulk Lorentz factors of a few are required. The radio maps of the two jets not detected by either Chandra or HST suggest that they are less beamed at large scales than the other four detected sources. Our results demonstrate that, at the sensitivity and resolution of Chandra, X-ray emission from extragalactic jets is common, yielding essential information on their physical properties.

We present the BeppoSAX observations of 3C 279 performed in 1997 January 13-23 simultaneously with gamma-ray pointings with EGRET on board CGRO and optical observations at the Torino Observatory (1997 January 11-18). ISO data close to this epoch are also available. We compare the derived Spectral Energy Distribution with those obtained at all other epochs with adequate multiwavelength coverage. Simple spectral models fitted to the multiepoch SEDs suggest that changes in intensity can be ascribed mainly to the variation of the bulk Lorentz factor of the plasma in the jet.

Three uninterrupted, long (lasting respectively 7, 10, and 10 days) ASCA observations of the well-studied TeV-bright blazars Mrk 421, Mrk 501 and PKS 2155-304 all show continuous strong X-ray flaring. Despite the relatively faint intensity states in 2 of the 3 sources, there was no identifiable quiescent period in any of the observations. Structure function analysis shows that all blazars have a characteristic time scale of ~ a day, comparable to the recurrence time and to the time scale of the stronger flares. On the other hand, examination of these flares in more detail reveals that each of the strong flares is not a smooth increase and decrease, but exhibits substructures of shorter flares having time scales of ~10 ks. We verify via simulations that in order to explain the observed structure function, these shorter flares ("shots") are unlikely to be fully random, but in some way are correlated with each other. The energy dependent cross-correlation analysis shows that inter-band lags are not universal in TeV blazars. This is important since in the past, only positive detections of lags were reported. In this work, we determine that the sign of a lag may differ from flare to flare; significant lags of both signs were detected from several flares, while no significant lag was detected from others. However, we also argue that the nature of the underlying component can affect these values. The facts that all flares are nearly symmetric and that fast variability shorter than the characteristic time scale is strongly suppressed, support the scenario where the light crossing time dominates the variability time scales of the day-scale flares.

We compare the host galaxies properties of BL Lac objects with those of radio loud quasars (RLQs) imaged by the WFPC2 on board of HST. The considered objects (z<0.5) are always well resolved and their host galaxies satisfactorily modelled by ellipticals. After homogeneous treatment of the data we find RLQs hosts are systematically more luminous (by ~0.7 mag) with respect to the hosts of BL Lacs. Using the M_BH - L_bulge relation, derived for nearby elliptical galaxies, we have evaluated the central black hole masses of our sample of active galaxies. These data are discussed in conjunction with the nuclear luminosity and the Eddington ratio.

We report the detection at X-rays of the radio/optical jet of 3C 371, from a short (10 ks) Chandra exposure in March 2000. We also present a new MERLIN observation at 1.4 GHz together with a renalysis of the archival HST WFPC2 F555W image. Despite the limited signal-to-noise ratio of the Chandra data, the X-ray morphology is clearly different from that of the radio/optical emission, with the brightest X-ray knot at 1.7" from the nucleus and little X-ray emission from the brightest radio/optical knot at 3.1". We construct the spectral energy distributions for the two emission regions at 1.7" and 3.1". Both show that the X-ray flux is below the extrapolation from the radio-to-optical continuum, suggesting moderately beamed synchrotron from an electron population with decreasing high energy cut-off as a plausible emission mechanism.

We present new BeppoSAX observations of seven BL Lacertae objects selected from the 1 Jy sample plus one additional source. The collected data cover the energy range 0.1 - 10 keV (observer's frame), reaching ~ 50 keV for one source (BL Lac). All sources characterized by a peak in their multifrequency spectra at infrared/optical energies (i.e., of the LBL type) display a relatively flat (alpha_x ~ 0.9) X-ray spectrum, which we interpret as inverse Compton emission. Four objects (2/3 of the LBL) show some evidence for a low-energy steepening which is likely due to the synchrotron tail merging into the inverse Compton component around ~ 1 - 3 keV. If this were generally the case with LBL, it would explain why the 0.1 - 2.4 keV ROSAT spectra of our sources are systematically steeper than the BeppoSAX ones (by ~ 0.5 in alpha_x). The broad-band spectral energy distributions fully confirm this picture and a synchrotron inverse Compton model allows us to derive the physical parameters (intrinsic power, magnetic field, etc.) of our sources. Combining our results with those obtained by BeppoSAX on BL Lacs covering a wide range of synchrotron peak frequencies, we confirm and clarify the dependence of the X-ray spectral index on synchrotron peak frequency originally found in ROSAT data.

We test the hypothesis that radio galaxies are a random subset of otherwise normal elliptical galaxies. Starting with the observed optical luminosity functions for elliptical galaxies, we show that the probability of an elliptical forming a radio source is a continuous, increasing function of optical luminosity, proportional to L squared. With this probability function and the luminosity function of normal elliptical galaxies as input to Monte Carlo simulations, we reproduce the observed distribution of radio galaxies in the radio-optical luminosity plane. Our results show that radio galaxies are a luminosity-biased but otherwise random sample of elliptical galaxies. This unified view of radio-loud and radio-quiet ellipticals also explains the well known difference of ~0.5 mag in average optical luminosity between FRI and FRII radio galaxies as a simple selection effect. Specifically, FRII appear preferentially in smaller galaxies because both radio and optical luminosity functions are steep, so there is a negligible probability of observing a powerful radio source associated with a bright galaxy; no intrinsic physical differences between FRI and FRII host galaxies are required.

We present a detailed X-ray study of the Broad Line Radio Galaxy 3C382, observed with the BeppoSAX satellite in a very bright state. The continuum emission is well modeled with a power law that steepens at high energies, with an e-folding energy of about 120 keV. At soft energies a clear excess of emission is detected, which can not be explained solely by the extended thermal halo seen in a ROSAT HRI image. A second, more intense soft X-ray component, possibly related to an accretion disk, is required by the data. Both a reflection component (R=0.3) and an iron line (EW \sim 50) are detected, at levels much weaker than in Seyfert galaxies, suggesting a common origin. Combining our measurements with results from the literature we find that the iron line has remained approximately constant over 9 years while the continuum varied by a factor of 5. Thus the fluorescent gas does not respond promptly to the variations of the X-ray primary source, suggesting that the reprocessing site is located away, likely at parsec distances. While the continuum shape indicates that X-rays derive from a thermal Comptonization process, the weakness of other spectral features implies that either the upper layers of the optically thick accretion disk are completely ionized or the corona above the disk is outflowing with mildly relativistic velocity.

The possibility of radio galaxies being random sample of otherwise normal elliptical galaxies is tested. Starting with the observed optical luminosity functions for elliptical galaxies, it is shown that the probability of an elliptical forming a radio source is a continuous, increasing function of optical luminosity, precisely proportional to square of the optical luminosity of the galaxy. Once the probability function is fixed, the luminosity function of normal elliptical galaxies is used as input for Monte Carlo simulations that reproduce the distribution of radio galaxies in the radio-optical luminosity plane. Our results show that radio galaxies are luminosity biased, but otherwise random sample of elliptical galaxies. This unified view of radio and non-radio ellipticals also explains the well known difference of 0.5 mag in average optical luminosity between FRI and FRII radio galaxies as a simple selection effect. Specifically, FRII appear preferentially in smaller galaxies because both radio and optical luminosity functions are steep, so there is a negligible probability of observing a powerful radio source associated with a bright galaxy; no intrinsic physical differences between FRI and FRII host galaxies is required.

We report results from Chandra observations of the X-ray jet of 3C~273 during the calibration phase in 2000 January. The zeroeth-order images and spectra from two 40-ks exposures with the HETG and LETG+ACIS-S show a complex X-ray structure. The brightest optical knots are detected and resolved in the 0.2-8 keV energy band. The X-ray morphology tracks well the optical. However, while the X-ray brightness decreases along the jet, the outer parts of the jet tend to be increasingly bright with increasing wavelength. The spectral energy distributions of four selected regions can best be explained by inverse Compton scattering of (beamed) cosmic microwave background photons. The model parameters are compatible with equipartition and a moderate Doppler factor, which is consistent with the one-sidedness of the jet. Alternative models either imply implausible physical conditions and energetics (the synchrotron self-Compton model) or are sufficiently ad hoc to be unconstrained by the present data (synchrotron radiation from a spatially or temporally distinct particle population).

As an overview of blazar research, we briefly discuss some of the themes of this conference, including the various types of blazars, blazar samples and selection effects, unified schemes, evolution, and jet physics.

We conducted a multi-frequency campaign for the TeV blazar Mrk~421 in 1998 April. The campaign started from a pronounced high amplitude flare recorded by SAX and Whipple; ASCA observation started three days later. In the X-ray data, we detected multiple flares, occuring on time scales of about one day. ASCA data clearly reveal spectral variability. The comparison of the data from ASCA, EUVE and RXTE indicates that the variability amplitudes in the low energy synchrotron component are larger at higher photon energies. In TeV Gamma-rays, large intra-day variations -- which were correlated with the X-ray flux -- were observed when results from three Cherenkov telescopes are combined. The RMS variability of TeV Gamma--rays was similar to that observed in hard X-rays, above 10 keV. The X-ray light curve reveals flares which are almost symmetric for most of cases, implying the dominant time scale is the light crossing time through the emitting region. The structure function analysis based on the continuous X-ray light curve of seven days indicates that the characteristic time scale is ~0.5 day. The analysis of ASCA light curves in various energy bands appears to show both soft (positive) and hard (negative) lags. These may not be real, as systematic effects could also produce these lags, which are all much smaller than an orbit. If the lags of both signs are real, these imply that the particle acceleration and X-ray cooling time scales are similar.

We propose that the X-ray emission detected by Chandra from the 100-kiloparsec-scale jet of PKS 0637-752 is produced through inverse Compton scattering of the Cosmic Microwave Background (CMB). We analyze the physical state of the jet and show that inverse Compton scattering from the CMB is consistent with equipartition for a moderate beaming of the emission, with Doppler factor delta~10. The power transported by the jet is then similar to that of other powerful blazars, Lj~10^{48} erg/s and the jet has low radiative efficiency. The radiative cooling times of the electrons are a few thousand years, compatible with the size of the knot. The low-energy cutoff of the electron distribution is constrained to be gamma_min~10, the first such constraint from spectral considerations. A parallel analysis for the SSC model yields far less reasonable physical conditions.

Empirical evidence suggests that elliptical galaxies hosting a radio source may not be different from normal non-radio ellipticals. To test this possibility, we use Monte Carlo simulations to reproduce the distribution of radio galaxies in the radio-optical luminosity plane. The input parameters of the simulation are the optical luminosity function (LF) of ellipticals and the radio LF of radio galaxies, linked by a function giving the probability for an elliptical to be a radio galaxy. Simulations reproduce the observations well, supporting unification of radio and non-radio ellipticals, provided that the probability of an elliptical hosting a radio source is proportional to the square of its optical luminosity. The difference of 0.5 magnitudes in average optical luminosity between FRI and FRII radio galaxies is also explained in this framework.

We present the preliminary results of a long BeppoSAX observation of the BL Lac object Mkn501 carried out in June 1999. The source was fainter than found during the BeppoSAX pointings of 1997 and 1998, but is still detected with a good signal-to-noise ratio up to ~40 keV. The X-ray spectrum in the energy range 0.1-40 keV, produced through synchrotron radiation, is steeper than in the previous years, it is clearly curved, and peaks (in nu*F_nu) at ~0.5 keV. This energy is much lower than those at which the synchrotron component was found to peak in 1997 and 1998. Some intraday variability suggests that activity of the source on small time scales accompanies the large long time scale changes of brightness and spectrum.

We report on the optical properties of a sample of 30 BL Lac host galaxies in the redshift range $0.03<z<0.2$, as derived from HST observations. All galaxies are fully resolved in the WFPC2 (F702W filter) images, allowing a quantitative analysis in two dimensions. Most and possibly all these galaxies have characteristics very similar to those of ``normal'' giant ellipticals. The luminosity, ellipticity, isophote twisting and amount of disky or boxy isophotes are consistent with those found in non-active ellipticals and in radio galaxies. In all cases the BL Lac nucleus is well centered in the main body of its host galaxy, a result that argues strongly against the microlensing hypothesis for any significant fraction of the population. A search for faint sub-structures in the host galaxies has not revealed notable signatures of tidal distortions or sub-components (faint disks, bars, X features, etc.), and with only one exception, there are no prominent dusty features in the central regions. Instead, the BL Lac host galaxies are smooth and unperturbed, suggesting that strong external gravitational interactions are not important to ongoing activity. A careful examination of the environment around the nucleus, however, shows a high incidence of close companion objects, whose nature remains unclear pending spectroscopic observations.

The HST NICMOS Camera 2 was used for H-band imaging of 12 BL Lacertae objects taken from the larger sample observed with the WFPC2 in the R band (Urry et al. 2000; Scarpa et al. 2000). Ten of the 12 BL Lacs are clearly resolved, and the detected host galaxies are large, bright ellipticals with average H-band absolute magnitude M=-26.2+-0.45 mag and effective radius 10+-5 kpc. The rest-frame integrated color of the host galaxies is on average R-H=2.3+-0.3, consistent with the value for both radio galaxies and normal, non-active elliptical galaxies, and indicating the dominant stellar population is old. The host galaxies tend to be bluer in their outer regions than in their cores, with average color gradient Delta(R-H)/Delta(log r)=-0.2 mag, again consistent with results for normal non-active elliptical galaxies. The infrared Kormendy relation, derived for the first time for BL Lac host galaxies, is m(e) = 3.8*log(R)+14.8 (where m(e) is the surface brightness at the effective radius R), fully in agreement with the relation for normal ellipticals. The close similarity between BL Lac host galaxies and normal ellipticals suggests the active nucleus has surprisingly little effect on the host galaxy. This supports a picture in which all elliptical galaxies harbor black holes which can be actively accreting for some fraction of their lifetime.

We briefly review BeppoSAX observations of X-ray bright radio-galaxies. Their X-ray spectra are quite varied, and perhaps surprisingly, any similarity between radio-loud AGN and Seyfert galaxies is the exception rather than the rule. When detected, reprocessing features (iron line and reflection) are generally weak, suggesting two possible scenarios: either: (1) non-thermal (jet?) radiation dilutes the X-ray emission from the disk in radio-loud objects, or (2) the solid angle subtended by the X-ray reprocessing material is smaller in radio-loud than in radio-quiet AGN due to different characteristics of the accretion disk itself.

We present the first observations at medium-hard X-rays with ASCA in 1998 August--November of four Flat Spectrum Radio Quasars (FSRQs), characterized by unusually steep soft X-ray spectra (photon index, \Gamma_{0.2-2.4 keV} \sim 2-2.5), as previously measured with ROSAT. Such steep X-ray slopes are similar to those observed in synchrotron-dominated BL Lacs and are unexpected in the context of the recent blazar paradigm, where sources with strong emission lines (such as FSRQs) are dominated in soft X-rays by a flat inverse Compton tail. We find that the ASCA spectra of the four FSRQs are consistent with a power law model with \Gamma_{2-10 keV} \sim 1.8, flatter than their ROSAT spectra. This indicates the onset of an inverse Compton component at energies \gtrsim 2 keV, in agreement with the blazar unification scheme. However, these objects are still anomalous within the blazar class for their steep soft X-ray continua which, together with non-simultaneous data at longer wavelengths, hint at the possibility that the synchrotron emission extends to soft X-rays. This would imply an anomalously high synchotron peak frequency for a quasar with luminous broad lines, challenging current blazar unification schemes. Alternatively, a plausible explanation for the steep optical-to-soft X-ray continua of the four FSRQs is thermal emission from the accretion disk, similar to the blazars 3C~273 and 3C~345. In the Appendix, we present fits to the SIS data in an effort to contribute to the ongoing calibration of the the time-dependence of the SIS response at low energies.

The energy budget of all known optical jets is discussed. It is found that to power the extended radio lobes of radio galaxies, the jet bulk motion on kpc scales must be relativistic, on average. Based on various constraints, a ``most probable'' region centered around bulk Lorentz factor ~7.5 and jet inclination \~20 degrees is found. Because of the consequent relativistic beaming, the rest frame magnetic field is lower and electron lifetimes longer. Combining the effects of time dilation and lower emission rate, the electron diffusion length becomes fully consistent with the deprojected jet length, without the need for reacceleration.

We report on a large HST imaging survey of BL Lac objects, at spatial resolution ~10 times better than previous ground-based surveys. We focus on data reduction and analysis, describing the procedures used to model the host galaxy surface brightness radial profiles. A total of 69 host galaxies were resolved out of 110 objects observed, including almost all sources at z < 0.5. We classify them morphologically by fitting with either an exponential disk or a de~Vaucouleurs profile; when one fit is preferred over the other, in 58 of 69 cases, it is invariably the elliptical morphology. This is a very strong result given the large number of BL Lac objects, the unprecedented spatial resolution, and the homogeneity of the data set. With the present reclassification of the host galaxy of 1418+546 as an elliptical, there remain no undisputed examples of a disk galaxy hosting a BL Lac nucleus. This implies that, at 99% confidence, fewer than 7% of BL Lacs can be in disk galaxies. The apparent magnitude of the host galaxies varies with distance as expected if the absolute magnitudes are approximately the same, with a spread of +-1 mag, out to redshift z < 0.5. At larger redshifts, only 6 of 23 BL Lacs are resolved so the present data do not constrain possible luminosity evolution of the host galaxies. The collective Hubble diagram for BL Lac host galaxies and radio galaxies strongly supports their unification.

We have used the HST WFPC2 camera to survey 132 BL Lac objects comprising seven complete radio-, X-ray-, and optically-selected samples. We obtained useful images for 110 targets spanning the redshift range 0 < z < 1.3. In two thirds of the BL Lac images, host galaxies are detected, including nearly all for z < 0.5 (58 of 63). The highest redshift host galaxy detected is in a BL Lac object at z=0.664. In 58 of the 72 resolved host galaxies, a de Vaucouleurs profile is significantly preferred, at >99% confidence, over a pure exponential disk; the two fits are comparable in the remaining 14 cases. These results limit the number of disk systems to at most 8% of BL Lacs (at 99% confidence), and are consistent with all BL~Lac host galaxies being ellipticals. The detected host galaxies are luminous ellipticals with a median absolute K-corrected magnitude of M_R= -23.7 +- 0.6 mag, at least one magnitude brighter than M* and comparable to brightest cluster galaxies. The galaxy morphologies are generally smooth and undisturbed, with small or negligible ellipticities (<0.2). There is no correlation between host galaxy and observed nuclear magnitude or estimated jet power corrected for beaming. If black hole mass is correlated linearly with bulge mass in general, this implies a large range in Eddington ratio. Present data strongly support the unification picture with FR I galaxies constituting the bulk of the parent population of BL Lac objects.

We analyze images of BL Lacertae objects obtained with the {\it Hubble Space Telescope} WFPC2 and the F814W filter. The nine objects cover a redshift range of 0.19 to 0.997. The relatively deep images are sufficient to detect galaxies at least one magnitude below M$^{*}_{I}$ (--21.4) and in most cases to three magnitudes below M$^{*}$. Galaxy enhancement over the average background is found around four out of the nine objects. Results for some cases are confirmed by ground-based imaging. In the other cases, the redshifts of the target BL Lac objects may be incorrect or they are truly isolated. These findings reinforce the idea that on average, BL Lac objects are found in regions of above average galaxy density. However, isolated objects apparently can host BL Lac nuclei too, a result that has implications for the processes that trigger/fuel the nuclear activity.

The TeV blazar PKS 2155-304 was monitored with the X-ray satellite ASCA in 1994 May, as part of a multiwavelength campaign from the radio to X-ray bands. At the beginning of the two-day continuous observation, we detected a large flare, where the 2-10 keV flux changed by a factor of 2 on a time scale of 3 x 10^4 sec. During the flare, the increase in the hard X-ray flux clearly preceded that observed in the soft X-rays, the spectral evolution tracking a `clockwise loop' in the flux versus photon index plane. Ascribing the energy-dependent variability to differential synchrotron cooling of relativistic electrons, we estimate the magnetic field B in the emission region. We tested two different methods of comparing the time series in various X-ray bands: (i) fitting the light curves to a Gaussian and searching for the time shift of the peak of the flare, and (ii) calculating the discrete correlation function. Both methods yielded a consistent solution of B ~ 0.1 G. We also found that the flare amplitude becomes larger as the photon energy increases, while the duration of the flare stays roughly constant through the ASCA energy band (0.7-7.5 keV). In the framework of the time dependent synchrotron self-Compton model in a homogeneous region, we consider a flare where the maximum Lorentz factor (Gmax) of the injected electrons increases uniformly throughout the emission volume. The temporal evolution of spectra as well as the light curves were reproduced with the physical parameters self-consistently determined from seven observables. We obtain B ~ 0.1-0.2 G and a region size R ~ 10^-2 pc, for relativistic beaming with a Doppler factor of D ~ 20-30. We discuss the significance of light travel time effects.

(Abridged) We present a detailed power density spectrum and cross-correlation analysis of the X-ray light curves of the BL Lac object PKS 2155-304, observed with BeppoSAX in 1997 and 1996, aimed at exploring the rapid variability properties and the inter-band cross correlations in the X-rays. We also perform the same analysis on the (archival) X-ray light curve obtained with ASCA in 1994.

HST observations have led to the discovery of the optical counterpart of the radio jet of PKS 2201+044, and to a detailed analysis of the optical jets of PKS 0521-365 and 3C371. At HST spatial resolution these jets are well resolved, displaying knotty morphologies. When compared with radio maps of appropriate resolution, a clear one-to-one correspondence between optical and radio structures is found, showing that all detected optical structures are indeed related to the radio synchrotron emission. Photometry of the brightest knots shows that the radio-to-optical spectral index and the derived intensity of the equipartition magnetic field are approximately constant along the jet. Thus, present observations suggest that the electron energy distribution does not change significantly all along the jet.

The gamma-ray blazar 3C 279 was monitored on a nearly daily basis with IUE, ROSAT and EGRET for three weeks between December 1992 and January 1993. During this period, the blazar was at a historical minimum at all wavelengths. Here we present the UV data obtained during the above multiwavelength campaign. A maximum UV variation of ~50% is detected, while during the same period the X-ray flux varied by no more than 13%. At the lowest UV flux level the average spectrum in the 1230-2700 A interval is unusually flat for this object (<alpha_UV>~1). The flattening could represent the lowest energy tail of the inverse Compton component responsible for the X-ray emission, or could be due to the presence of a thermal component at ~20000 K possibly associated with an accretion disk. The presence of an accretion disk in this blazar object, likely observable only in very low states and otherwise hidden by the beamed, variable synchrotron component, would be consistent with the scenario in which the seed photons for the inverse Compton mechanism producing the gamma-rays are external to the relativistic jet. We further discuss the long term correlation of the UV flux with the X-ray and gamma-ray fluxes obtained at various epochs. All UV archival data are included in the analysis. Both the X- and gamma-ray fluxes are generally well correlated with the UV flux, approximately with square root and quadratic dependences, respectively.