Search

We investigate the population of high-redshift ($3\leq z < 6$) AGN selected in the two deepest X-ray surveys, the 7 Ms \textit{Chandra} Deep Field-South and 2 Ms \textit{Chandra} Deep Field-North. Their outstanding sensitivity and spectral characterization of faint sources allow us to focus on the sub-$L_*$ regime (log$L_{\mathrm{X}}\lesssim44$), poorly sampled by previous works using shallower data, and the obscured population. Taking fully into account the individual photometric-redshift probability distribution functions, the final sample consists of $\approx102$ X-ray selected AGN at $3\leq z < 6$. The fraction of AGN obscured by column densities log$N_{\mathrm{H}}>23$ is $\sim0.6-0.8$, once incompleteness effects are taken into account, with no strong dependence on redshift or luminosity. We derived the high-redshift AGN number counts down to $F_{\mathrm{0.5-2\,keV}}=7\times10^{-18}\,\mathrm{erg\,cm^{-2}\,s^{-1}}$, extending previous results to fainter fluxes, especially at $z>4$. We put the tightest constraints to date on the low-luminosity end of AGN luminosity function at high redshift. The space-density, in particular, declines at $z>3$ at all luminosities, with only a marginally steeper slope for low-luminosity AGN. By comparing the evolution of the AGN and galaxy densities, we suggest that such a decline at high luminosities is mainly driven by the underlying galaxy population, while at low luminosities there are hints of an intrinsic evolution of the parameters driving nuclear activity. Also, the black-hole accretion rate density and star-formation rate density, which are usually found to evolve similarly at $z\lesssim3$, appear to diverge at higher redshifts., Comment: Accepted 2017 September 15 (MNRAS). 31 Pages, 25 figures

We study the X-ray variability properties of distant AGNs in the Chandra Deep Field-South region over 17 years, up to $z\sim 4$, and compare them with those predicted by models based on local samples. We use the results of Monte Carlo simulations to account for the biases introduced by the discontinuous sampling and the low-count regime. We confirm that variability is an ubiquitous property of AGNs, with no clear dependence on the density of the environment. The variability properties of high-z AGNs, over different temporal timescales, are most consistent with a Power Spectral Density (PSD) described by a broken (or bending) power-law, similar to nearby AGNs. We confirm the presence of an anti-correlation between luminosity and variability, resulting from the dependence of variability on BH mass and accretion rate. We explore different models, finding that our acceptable solutions predict that BH mass influences the value of the PSD break frequency, while the Eddington ratio $\lambda_{Edd}$ affects the PSD break frequency and, possibly, the PSD amplitude as well. We derive the evolution of the average $\lambda_{Edd}$ as a function of redshift, finding results in agreement with measurements based on different estimators. The large statistical uncertainties make our results consistent with a constant Eddington ratio, although one of our models suggest a possible increase of $\lambda_{Edd}$ with lookback time up to $z\sim 2-3$. We conclude that variability is a viable mean to trace the accretion history of supermassive BHs, whose usefulness will increase with future, wide-field/large effective area X-ray missions., Comment: 15 pages, 9 Figures, 2 Tables, in press on MNRAS

We investigate the dependence of black-hole accretion rate (BHAR) on host-galaxy star formation rate (SFR) and stellar mass ($M_*$) in the CANDELS/GOODS-South field in the redshift range of $0.5\leq z < 2.0$. Our sample consists of $\approx 18000$ galaxies, allowing us to probe galaxies with $0.1 \lesssim \mathrm{SFR} \lesssim 100\ M_\odot\ \mathrm{yr}^{-1}$ and/or $10^8 \lesssim M_* \lesssim 10^{11}\ M_{\odot}$. We use sample-mean BHAR to approximate long-term average BHAR. Our sample-mean BHARs are derived from the $Chandra$ Deep Field-South 7 Ms observations, while the SFRs and $M_*$ have been estimated by the CANDELS team through SED fitting. The average BHAR is correlated positively with both SFR and $M_*$, and the BHAR-SFR and BHAR-$M_*$ relations can both be described acceptably by linear models with a slope of unity. However, BHAR appears to be correlated more strongly with $M_*$ than SFR. This result indicates that $M_*$ is the primary host-galaxy property related to black-hole growth, and the apparent BHAR-SFR relation is largely a secondary effect due to the star-forming main sequence. Among our sources, massive galaxies ($M_* \gtrsim 10^{10} M_{\odot}$) have significantly higher BHAR/SFR ratios than less-massive galaxies, indicating the former have higher black-hole fueling efficiency and/or higher SMBH occupation fraction than the latter. Our results can naturally explain the observed proportionality between $M_{\rm BH}$ and $M_*$ for local giant ellipticals, and suggest their $M_{\rm BH}/M_*$ is higher than that of local star-forming galaxies. Among local star-forming galaxies, massive systems might have higher $M_{\rm BH}/M_*$ compared to dwarfs., Comment: 18 pages, 11 figures, 3 tables; accepted for publication in ApJ

We present $Chandra$ observations of seven broad absorption line (BAL) quasars at $z=0.863$-2.516 with redshifted BAL troughs (RSBALs). Five of our seven targets were detected by $Chandra$ in 4-13 ks exposures with ACIS-S. The $\alpha_{\rm ox}$ values, $\Delta\alpha_{\rm ox}$ values, and spectral energy distributions of our targets demonstrate they are all X-ray weak relative to expectations for non-BAL quasars, and the degree of X-ray weakness is consistent with that of appropriately-matched BAL quasars generally. Furthermore, our five detected targets show evidence for hard X-ray spectral shapes with a stacked effective power-law photon index of $\Gamma_{\rm eff}=0.5^{+0.5}_{-0.4}$. These findings support the presence of heavy X-ray absorption ($N_{\rm H}\approx 2 \times 10^{23}$ cm$^{-2}$) in RSBAL quasars, likely by the shielding gas found to be common in BAL quasars more generally. We use these X-ray measurements to assess models for the nature of RSBAL quasars, finding that a rotationally-dominated outflow model is favored while an infall model also remains plausible with some stipulations. The X-ray data disfavor a binary quasar model for RSBAL quasars in general., Comment: 11 pages, 5 figures, and 3 tables

We report on the detection of a remarkable new fast high-energy transient found in the Chandra Deep Field-South, robustly associated with a faint ($m_{\rm R}=27.5$ mag, $z_{\rm ph}$$\sim$2.2) host in the CANDELS survey. The X-ray event is comprised of 115$^{+12}_{-11}$ net 0.3-7.0 keV counts, with a light curve characterised by a $\approx$100 s rise time, a peak 0.3-10 keV flux of $\approx$5$\times$10$^{-12}$ erg s$^{-1}$ cm$^{-2}$, and a power-law decay time slope of $-1.53\pm0.27$. The average spectral slope is $\Gamma=1.43^{+0.23}_{-0.13}$, with no clear spectral variations. The \hbox{X-ray} and multi-wavelength properties effectively rule out the vast majority of previously observed high-energy transients. A few theoretical possibilities remain: an "orphan" X-ray afterglow from an off-axis short-duration Gamma-ray Burst (GRB) with weak optical emission; a low-luminosity GRB at high redshift with no prompt emission below $\sim$20 keV rest-frame; or a highly beamed Tidal Disruption Event (TDE) involving an intermediate-mass black hole and a white dwarf with little variability. However, none of the above scenarios can completely explain all observed properties. Although large uncertainties exist, the implied rate of such events is comparable to those of orphan and low-luminosity GRBs as well as rare TDEs, implying the discovery of an untapped regime for a known transient class, or a new type of variable phenomena whose nature remains to be determined., Comment: 18 pages, 7 figures, 4 tables. Accepted for publication in MNRAS

We present a sample of 10 low-mass active galactic nuclei (AGNs) selected from the 40-month NuSTAR serendipitous survey. The sample is selected to have robust NuSTAR detections at $3 - 24$~keV, to be at $z < 0.3$, and to have optical r-band magnitudes at least 0.5~mag fainter than an $L_\star$ galaxy at its redshift. The median values of absolute magnitude, stellar mass and 2--10 X-ray luminosity of our sample are $\langle M_r\rangle = -20.03$, $\langle M_\star\rangle = 4.6\times10^{9}M_\odot$, and $\langle L_{2-10\mathrm{keV}}\rangle = 3.1\times10^{42}$ erg s$^{-1}$, respectively. Five objects have detectable broad H$\alpha$ emission in their optical spectra, indicating black-hole masses of $(1.1-10.4)\times 10^6 M_\odot$. We find that $30^{+17}_{-10}\%$ of the galaxies in our sample do not show AGN-like optical narrow emission lines, and one of the ten galaxies in our sample, J115851+4243.2, shows evidence for heavy X-ray absorption. This result implies that a non-negligible fraction of low-mass galaxies might harbor accreting massive black holes that are missed by optical spectroscopic surveys and $<10$ keV X-ray surveys. The mid-IR colors of our sample also indicate these optically normal low-mass AGNs cannot be efficiently identified with typical AGN selection criteria based on WISE colors. While the hard ($>10$ keV) X-ray selected low-mass AGN sample size is still limited, our results show that sensitive NuSTAR observations are capable of probing faint hard X-ray emission originating from the nuclei of low-mass galaxies out to moderate redshift ($z<0.3$), thus providing a critical step in understanding AGN demographics in low-mass galaxies., Comment: Accepted for publication in ApJ. 27 pages, 6 tables, 10 figures, including one online figure set available at https://ctjchen.github.io/ctc_nustardwarf_figset.pdf (v2: updated acknowledgement.)

We explore the dependence of the incidence of moderate-luminosity ($L_{X} = 10^{41.9-43.7}$ erg s$^{-1}$) AGNs and the distribution of their accretion rates on host color at 0.5 < z < 2.5, using deep X-ray data in GOODS fields. We use extinction-corrected rest-frame U-V colors to divide both AGN hosts and non-AGN galaxies into red sequence (quiescent), green valley (transition), and blue cloud (star-forming) populations. We find that both the AGN fraction at fixed stellar mass and its evolution with redshift are dependent on host colors. Most notably, red galaxies have the lowest AGN fraction (~5\%) at z~1 yet with most rapid redshift evolution, increasing by a factor of 5 (~24\%) at z~2. Green galaxies exhibit the highest AGN fraction across all redshifts, which is most pronounced at z~2 with more than half of them hosting an AGN at $M_{*} > 10^{10.6} M_{\odot}$. Together with the high AGN fraction in red galaxies at z~2, this indicates that X-ray AGNs could be important in both transforming blue galaxies into red ones and subsequently maintaining their quiescence at high redshift. Furthermore, consistent with low-redshift studies, we find that the probability of hosting an AGN in the total galaxy population can be characterized by a universal Eddington ratio ($p(\lambda_{Edd}) \sim \lambda_{Edd}^{-0.4}$) and a moderate redshift evolution. Yet consistent with their different AGN fractions, different populations appear to also have different $p(\lambda_{Edd})$ with red galaxies exhibiting more rapid redshift evolution than green and blue ones. Evidence for a steeper power law of $p(\lambda_{Edd})$ in red galaxies is also presented, though larger samples are needed to confirm. These results suggest that the AGN accretion or the growth of supermassive black holes is related to their host properties, and may also influence their hosts in a different mode dependent on the host color., Comment: A&A, in press

We present X-ray source catalogs for the $\approx7$ Ms exposure of the Chandra Deep Field-South (CDF-S), which covers a total area of 484.2 arcmin$^2$. Utilizing WAVDETECT for initial source detection and ACIS Extract for photometric extraction and significance assessment, we create a main source catalog containing 1008 sources that are detected in up to three X-ray bands: 0.5-7.0 keV, 0.5-2.0 keV, and 2-7 keV. A supplementary source catalog is also provided including 47 lower-significance sources that have bright ($K_s\le23$) near-infrared counterparts. We identify multiwavelength counterparts for 992 (98.4%) of the main-catalog sources, and we collect redshifts for 986 of these sources, including 653 spectroscopic redshifts and 333 photometric redshifts. Based on the X-ray and multiwavelength properties, we identify 711 active galactic nuclei (AGNs) from the main-catalog sources. Compared to the previous $\approx4$ Ms CDF-S catalogs, 291 of the main-catalog sources are new detections. We have achieved unprecedented X-ray sensitivity with average flux limits over the central $\approx1$ arcmin$^2$ region of $\approx1.9\times10^{-17}$, $6.4\times10^{-18}$, and $2.7\times10^{-17}$ erg cm$^{-2}$ s$^{-1}$ in the three X-ray bands, respectively. We provide cumulative number-count measurements observing, for the first time, that normal galaxies start to dominate the X-ray source population at the faintest 0.5-2.0 keV flux levels. The highest X-ray source density reaches $\approx50\,500$ deg$^{-2}$, and $47\%\pm4\%$ of these sources are AGNs ($\approx23\,900$ deg$^{-2}$)., Comment: 31 pages, 31 figures, 8 tables, ApJS in press, minor text updates; full-resolution PDF version and data products available at http://www.astro.psu.edu/users/niel/cdfs/cdfs-chandra.html

We perform long-term ($\approx 15$ yr, observed-frame) X-ray variability analyses of the 68 brightest radio-quiet active galactic nuclei (AGNs) in the 6 Ms $Chandra$ Deep Field-South (CDF-S) survey; the majority are in the redshift range of $0.6-3.1$, providing access to penetrating rest-frame X-rays up to $\approx 10-30$ keV. Twenty-four of the 68 sources are optical spectral type I AGNs, and the rest (44) are type II AGNs. The time scales probed in this work are among the longest for X-ray variability studies of distant AGNs. Photometric analyses reveal widespread photon-flux variability: $90\%$ of AGNs are variable above a 95% confidence level, including many X-ray obscured AGNs and several optically classified type II quasars. We characterize the intrinsic X-ray luminosity ($L_{\rm{X}}$) and absorption ($N_{\rm{H}}$) variability via spectral fitting. Most (74%) sources show $L_{\rm{X}}$ variability; the variability amplitudes are generally smaller for quasars. A Compton-thick candidate AGN shows variability of its high-energy X-ray flux, indicating the size of reflecting material to be $\lesssim 0.3$ pc. $L_{\rm{X}}$ variability is also detected in a broad absorption line (BAL) quasar. The $N_{\rm{H}}$ variability amplitude for our sample appears to rise as time separation increases. About 16% of sources show $N_{\rm{H}}$ variability. One source transitions from an X-ray unobscured to obscured state while its optical classification remains type I; this behavior indicates the X-ray eclipsing material is not large enough to obscure the whole broad-line region., Comment: 21 pages, 17 figures, accepted for publication in ApJ

eXTP is a science mission designed to study the state of matter under extreme conditions of density, gravity and magnetism. Primary targets include isolated and binary neutron stars, strong magnetic field systems like magnetars, and stellar-mass and supermassive black holes. The mission carries a unique and unprecedented suite of state-of-the-art scientific instruments enabling for the first time ever the simultaneous spectral-timing-polarimetry studies of cosmic sources in the energy range from 0.5-30 keV (and beyond). Key elements of the payload are: the Spectroscopic Focusing Array (SFA) - a set of 11 X-ray optics for a total effective area of about 0.9 m^2 and 0.6 m^2 at 2 keV and 6 keV respectively, equipped with Silicon Drift Detectors offering <180 eV spectral resolution; the Large Area Detector (LAD) - a deployable set of 640 Silicon Drift Detectors, for a total effective area of about 3.4 m^2, between 6 and 10 keV, and spectral resolution <250 eV; the Polarimetry Focusing Array (PFA) - a set of 2 X-ray telescope, for a total effective area of 250 cm^2 at 2 keV, equipped with imaging gas pixel photoelectric polarimeters; the Wide Field Monitor (WFM) - a set of 3 coded mask wide field units, equipped with position-sensitive Silicon Drift Detectors, each covering a 90 degrees x 90 degrees FoV. The eXTP international consortium includes mostly major institutions of the Chinese Academy of Sciences and Universities in China, as well as major institutions in several European countries and the United States. The predecessor of eXTP, the XTP mission concept, has been selected and funded as one of the so-called background missions in the Strategic Priority Space Science Program of the Chinese Academy of Sciences since 2011. The strong European participation has significantly enhanced the scientific capabilities of eXTP. The planned launch date of the mission is earlier than 2025., Comment: 16 pages, 16 figures. Oral talk presented at SPIE Astronomical Telescopes and Instrumentation, June 26 to July 1, 2016, Edingurgh, UK

We present measurements of the evolution of normal-galaxy X-ray emission from $z \approx$ 0-7 using local galaxies and galaxy samples in the 6 Ms Chandra Deep Field-South (CDF-S) survey. The majority of the CDF-S galaxies are observed at rest-frame energies above 2 keV, where the emission is expected to be dominated by X-ray binary (XRB) populations; however, hot gas is expected to provide small contributions to the observed- frame < 1 keV emission at $z < 1$. We show that a single scaling relation between X-ray luminosity ($L_{\rm X}$) and star-formation rate (SFR) is insufficient for characterizing the average X-ray emission at all redshifts. We establish that scaling relations involving not only SFR, but also stellar mass ($M_\star$) and redshift, provide significantly improved characterizations of the average X-ray emission from normal galaxy populations at $z \approx$ 0-7. We further provide the first empirical constraints on the redshift evolution of X-ray emission from both low-mass XRB (LMXB) and high-mass XRB (HMXB) populations and their scalings with $M_\star$ and SFR, respectively. We find $L_{\rm 2-10 keV}({\rm LMXB})/M_\star \propto (1 + z)^{2-3}$ and $L_{\rm 2-10 keV}({\rm HMXB})$/SFR $\propto (1 + z)$, and show that these relations are consistent with XRB population-synthesis model predictions, which attribute the increase in LMXB and HMXB scaling relations with redshift as being due to declining host galaxy stellar ages and metallicities, respectively. We discuss how emission from XRBs could provide an important source of heating to the intergalactic medium in the early Universe, exceeding that of active galactic nuclei., Comment: Accepted for publication in ApJ (23 pages, 13 figures, 5 tables)

PG1247+267 is one of the most luminous known quasars at $z\sim2$ and is a strongly super-Eddington accreting SMBH candidate. We obtained NuSTAR data of this intriguing source in December 2014 with the aim of studying its high-energy emission, leveraging the broad band covered by the new NuSTAR and the archival XMM-Newton data. Several measurements are in agreement with the super-Eddington scenario for PG1247+267: the soft power law ($\Gamma=2.3\pm0.1$); the weak ionized Fe emission line and a hint of the presence of outflowing ionized gas surrounding the SMBH. The presence of an extreme reflection component is instead at odds with the high accretion rate proposed for this quasar. This can be explained with three different scenarios; all of them are in good agreement with the existing data, but imply very different conclusions: i) a variable primary power law observed in a low state, superimposed on a reflection component echoing a past, higher flux state; ii) a power law continuum obscured by an ionized, Compton thick, partial covering absorber; and iii) a relativistic disk reflector in a lamp-post geometry, with low coronal height and high BH spin. The first model is able to explain the high reflection component in terms of variability. The second does not require any reflection to reproduce the hard emission, while a rather low high-energy cutoff of $\sim100$ keV is detected for the first time in such a high redshift source. The third model require a face-on geometry, which may affect the SMBH mass and Eddington ratio measurements. Deeper X-ray broad-band data are required in order to distinguish between these possibilities., Comment: 12 pages, 10 figure. Accepted for publication in A&A

We present improved point-source catalogs for the 2 Ms Chandra Deep Field-North (CDF-N) and the 250 ks Extended Chandra Deep Field-South (E-CDF-S), implementing a number of recent improvements in Chandra source-cataloging methodology. For the CDF-N/E-CDF-S, we provide a main catalog that contains 683/1003 X-ray sources detected with wavdetect at a false-positive probability threshold of $10^{-5}$ that also satisfy a binomial-probability source-selection criterion of $P<0.004$/$P<0.002$. Such an approach maximizes the number of reliable sources detected: a total of 196/275 main-catalog sources are new compared to the Alexander et al. (2003) CDF-N/Lehmer et al. (2005) E-CDF-S main catalogs. We also provide CDF-N/E-CDF-S supplementary catalogs that consist of 72/56 sources detected at the same wavdetect threshold and having $P$ of $0.004-0.1$/$0.002-0.1$ and $K_s\le22.9/K_s\le22.3$ mag counterparts. For all $\approx1800$ CDF-N and E-CDF-S sources, including the $\approx500$ newly detected ones (these being generally fainter and more obscured), we determine X-ray source positions utilizing centroid and matched-filter techniques; we also provide multiwavelength identifications, apparent magnitudes of counterparts, spectroscopic and/or photometric redshifts, basic source classifications, and estimates of observed AGN and galaxy source densities around respective field centers. Simulations show that both the CDF-N and E-CDF-S main catalogs are highly reliable and reasonably complete. Background and sensitivity analyses indicate that the on-axis mean flux limits reached represent a factor of $\approx1.5-2.0$ improvement over the previous CDF-N and E-CDF-S limits. We make our data products publicly available., Comment: 45 pages, 57 figures, 14 tables, ApJS in press; Data products and images available at http://www2.astro.psu.edu/users/niel/hdf/hdf-chandra.html and http://www2.astro.psu.edu/users/niel/ecdfs/ecdfs-chandra.html arXiv admin note: text overlap with arXiv:1105.5643

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

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

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., Comment: 20 pages, published in the Astrophysical Journal

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., Comment: 14 pages, accepted for publication in the Astrophysical Journal

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., Comment: 14 pages, 10 figures, accepted for publication in A&A

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

We present new NuSTAR and Chandra observations of NGC 3393, a galaxy reported to host the smallest separation dual AGN resolved in the X-rays. While past results suggested a 150 pc separation dual AGN, three times deeper Chandra imaging, combined with adaptive optics and radio imaging suggest a single, heavily obscured, radio-bright AGN. Using VLA and VLBA data, we find an AGN with a two-sided jet rather than a dual AGN and that the hard X-ray, UV, optical, NIR, and radio emission are all from a single point source with a radius <0.2". We find that the previously reported dual AGN is most likely a spurious detection resulting from the low number of X-ray counts (<160) at 6-7 keV and Gaussian smoothing of the data on scales much smaller than the PSF (0.25" vs. 0.80" FWHM). We show that statistical noise in a single Chandra PSF generates spurious dual peaks of the same separation (0.55$\pm$0.07" vs. 0.6") and flux ratio (39$\pm$9% vs. 32% of counts) as the purported dual AGN. With NuSTAR, we measure a Compton-thick source (NH=$2.2\pm0.4\times10^{24}$ cm$^{-2}$) with a large torus half-opening angle, {\theta}=79 which we postulate results from feedback from strong radio jets. This AGN shows a 2-10 keV intrinsic to observed flux ratio of 150. Using simulations, we find that even the deepest Chandra observations would severely underestimate the intrinsic luminosity of NGC 3393 above z>0.2, but would detect an unobscured AGN of this luminosity out to high redshift (z=5)., Comment: Accepted for publication in ApJ. 15 Figures and 4 tables

We present the results of recent Chandra, XMM-Newton, and Hubble Space Telescope observations of the radio-loud (RL), broad absorption line (BAL) quasar PG 1004+130. We compare our new observations to archival X-ray and UV data, creating the most comprehensive, high signal-to-noise, multi-epoch, spectral monitoring campaign of a RL BAL quasar to date. We probe for variability of the X-ray absorption, the UV BAL, and the X-ray jet, on month-year timescales. The X-ray absorber has a low column density of $N_{H}=8\times10^{20}-4\times10^{21}$ cm$^{-2}$ when it is assumed to be fully covering the X-ray emitting region, and its properties do not vary significantly between the 4 observations. This suggests the observed absorption is not related to the typical "shielding gas" commonly invoked in BAL quasar models, but is likely due to material further from the central black hole. In contrast, the CIV BAL shows strong variability. The equivalent width (EW) in 2014 is EW=11.24$\pm$0.56 \AA, showing a fractional increase of $\Delta EW / \langle EW \rangle$=1.16$\pm$0.11 from the 2003 observation, 3183 days earlier in the rest-frame. This places PG 1004+130 among the most highly variable BAL quasars. By combining Chandra observations we create an exposure 2.5 times deeper than studied previously, with which to investigate the nature of the X-ray jet and extended diffuse X-ray emission. An X-ray knot, likely with a synchrotron origin, is detected in the radio jet ~8 arcsec (30 kpc) from the central X-ray source with a spatial extent of ~4 arcsec (15 kpc). No similar X-ray counterpart to the counterjet is detected. Asymmetric, non-thermal diffuse X-ray emission, likely due to inverse Compton scattering of Cosmic Microwave Background photons, is also detected., Comment: 15 pages, 7 figures, 3 tables. Accepted for publication in ApJ

We present the infrared (IR) and X-ray properties of a sample of 33 mid-IR luminous quasars ($\nu$L(6 micron)>6x10$^{44}$ erg/s) at redshift z~1-3, identified through detailed spectral energy distribution analyses of distant star-forming galaxies, using the deepest IR data from Spitzer and Herschel in the GOODS-Herschel fields. The aim is to constrain the fraction of obscured, and Compton-thick (CT, N$_H$>1.5x10$^{24}$ cm$^{-2}$) quasars at the peak era of nuclear and star-formation activities. Despite being very bright in the mid-IR band, ~30% of these quasars are not detected in the extremely deep 2 Ms and 4 Ms Chandra X-ray data available in these fields. X-ray spectral analysis of the detected sources reveals that the majority (~67%) are obscured by column densities N$_H$>10$^{22}$ cm$^{-2}$; this fraction reaches ~80% when including the X-ray undetected sources (9 out of 33), which are likely to be the most heavily-obscured, CT quasars. We constrain the fraction of CT quasars in our sample to be ~24-48%, and their space density to be $\Phi$=(6.7$\pm$2.2)x10$^{-6}$ Mpc$^{-3}$. From the investigation of the quasar host galaxies in terms of star-formation rates (SFRs) and morphological distortions, as a sign of galaxy mergers/interactions, we do not find any direct relation between SFRs and quasar luminosity or X-ray obscuration. On the other hand, there is tentative evidence that the most heavily-obscured quasars have, on average, more disturbed morphologies than the unobscured/moderately-obscured quasar hosts, which preferentially live in undisturbed systems. However, the fraction of quasars with disturbed morphology amongst the whole sample is ~40%, suggesting that galaxy mergers are not the main fuelling mechanism of quasars at z~2., Comment: 21 pages, 10 figures, 4 tables. Changes made to address referee's comments (final version). Accepted for publication in MNRAS

We present an X-ray and multiwavelength study of 33 weak emission-line quasars (WLQs) and 18 quasars that are analogs of the extreme WLQ, PHL 1811, at z ~ 0.5-2.9. New Chandra 1.5-9.5 ks exploratory observations were obtained for 32 objects while the others have archival X-ray observations. Significant fractions of these luminous type 1 quasars are distinctly X-ray weak compared to typical quasars, including 16 (48%) of the WLQs and 17 (94%) of the PHL 1811 analogs with average X-ray weakness factors of 17 and 39, respectively. We measure a relatively hard ($\Gamma=1.16_{-0.32}^{+0.37}$) effective power-law photon index for a stack of the X-ray weak subsample, suggesting X-ray absorption, and spectral analysis of one PHL 1811 analog, J1521+5202, also indicates significant intrinsic X-ray absorption. We compare composite SDSS spectra for the X-ray weak and X-ray normal populations and find several optical-UV tracers of X-ray weakness; e.g., Fe II rest-frame equivalent width and relative color. We describe how orientation effects under our previously proposed "shielding-gas" scenario can likely unify the X-ray weak and X-ray normal populations. We suggest that the shielding gas may naturally be understood as a geometrically thick inner accretion disk that shields the broad line region from the ionizing continuum. If WLQs and PHL 1811 analogs have very high Eddington ratios, the inner disk could be significantly puffed up (e.g., a slim disk). Shielding of the broad emission-line region by a geometrically thick disk may have a significant role in setting the broad distributions of C IV rest-frame equivalent width and blueshift for quasars more generally., Comment: 26 pages, 18 figures. ApJ in press

We measure the location and evolutionary vectors of 69 Herschel-detected broad-line active galactic nuclei (BLAGNs) in the M_BH-M_* plane. BLAGNs are selected from the COSMOS and CDF-S fields, and span the redshift range 0.2< z<2.1. Black-hole masses are calculated using archival spectroscopy and single-epoch virial mass estimators, and galaxy total stellar masses are calculated by fitting the spectral energy distribution (subtracting the BLAGN component). The mass-growth rates of both the black hole and galaxy are calculated using Chandra/XMM-Newton X-ray and Herschel far-infrared data, reliable measures of the BLAGN accretion and galaxy star formation rates, respectively. We use Monte Carlo simulations to account for biases in our sample, due to both selection limits and the steep slope of the massive end of the galaxy stellar-mass distribution. We find our sample is consistent with no evolution in the M_BH-M_* relation from z~2 to z~0. BLAGNs and their host galaxies which lie off the black hole mass galaxy total stellar mass relation tend to have evolutionary vectors anti-correlated with their mass ratios: that is, galaxies with over-massive (under-massive) black holes tend to have a low (high) ratio of the specific accretion rate to the specific star formation rate. We also use the measured growth rates to estimate the preferred AGN duty cycle for our galaxies to evolve most consistently onto the local M_BH-M_Bul relation. Under reasonable assumptions of exponentially declining star formation histories, the data suggest a non-evolving (no more than a factor of a few) BLAGN duty cycle among star-forming galaxies of 10% (1sigma range of 1-42% at z<1 and 2-34% at z>1)., Comment: 26 pages, 13 figures, 2 tables, ApJ accepted

Ultra-deep observations of ECDF-S with Chandra and XMM-Newton enable a search for extended X-ray emission down to an unprecedented flux of $2\times10^{-16}$ ergs s$^{-1}$ cm$^{-2}$. We present the search for the extended emission on spatial scales of 32$^{\prime\prime}$ in both Chandra and XMM data, covering 0.3 square degrees and model the extended emission on scales of arcminutes. We present a catalog of 46 spectroscopically identified groups, reaching a redshift of 1.6. We show that the statistical properties of ECDF-S, such as logN-logS and X-ray luminosity function are broadly consistent with LCDM, with the exception that dn/dz/d$\Omega$ test reveals that a redshift range of $0.2

We derive photometric redshifts (\zp) for sources in the entire ($\sim0.4$ deg$^2$) Hawaii-Hubble Deep Field-North (\hdfn) field with the EAzY code, based on point spread function-matched photometry of 15 broad bands from the ultraviolet (\bandu~band) to mid-infrared (IRAC 4.5 $\mu$m). Our catalog consists of a total of 131,678 sources. We evaluate the \zp~quality by comparing \zp~with spectroscopic redshifts (\zs) when available, and find a value of normalized median absolute deviation \sigm$=$0.029 and an outlier fraction of 5.5\% (outliers are defined as sources having $\rm |\zp - \zs|/(1+\zs) > 0.15$) for non-X-ray sources. More specifically, we obtain \sigm$=0.024$ with 2.7\% outliers for sources brighter than $R=23$~mag, \sigm$=0.035$ with 7.4\% outliers for sources fainter than $R=23$~mag, \sigm$=$0.026 with 3.9\% outliers for sources having $z<1$, and \sigm$=$0.034 with 9.0\% outliers for sources having $z>1$. Our \zp\ quality shows an overall improvement over an earlier \zp\ work that focused only on the central \hdfn\ area. We also classify each object as star or galaxy through template spectral energy distribution fitting and complementary morphological parametrization, resulting in 4959 stars and 126,719 galaxies. Furthermore, we match our catalog with the 2~Ms {\it Chandra} Deep Field-North main \xray~catalog. For the 462 matched non-stellar \xray~sources (281 having \zs), we improve their \zp~quality by adding three additional AGN templates, achieving \sigm$=0.035$ and an outlier fraction of 12.5\%. We make our catalog publicly available presenting both photometry and \zp, and provide guidance on how to make use of our catalog., Comment: 20 pages, 14 figures, 6 tables, ApJS in press; photo-z catalog available upon request (gxy909@psu.edu)

We present X-ray spectral analyses for three Seyfert 2 active galactic nuclei, NGC 424, NGC 1320, and IC 2560, observed by NuSTAR in the 3-79 keV band. The high quality hard X-ray spectra allow detailed modeling of the Compton reflection component for the first time in these sources. Using quasi-simultaneous NuSTAR and Swift/XRT data, as well as archival XMM-Newton data, we find that all three nuclei are obscured by Compton-thick material with column densities in excess of ~5 x $10^{24}$ cm$^{-2}$, and that their X-ray spectra above 3 keV are dominated by reflection of the intrinsic continuum on Compton-thick material. Due to the very high obscuration, absorbed intrinsic continuum components are not formally required by the data in any of the sources. We constrain the intrinsic photon indices and the column density of the reflecting medium through the shape of the reflection spectra. Using archival multi-wavelength data we recover the intrinsic X-ray luminosities consistent with the broadband spectral energy distributions. Our results are consistent with the reflecting medium being an edge-on clumpy torus with a relatively large global covering factor and overall reflection efficiency of the order of 1%. Given the unambiguous confirmation of the Compton-thick nature of the sources, we investigate whether similar sources are likely to be missed by commonly used selection criteria for Compton-thick AGN, and explore the possibility of finding their high-redshift counterparts., Comment: Accepted for publication in ApJ

We report NuSTAR observations of a sample of six X-ray weak broad absorption line (BAL) quasars. These targets, at z=0.148-1.223, are among the optically brightest and most luminous BAL quasars known at z<1.3. However, their rest-frame 2 keV luminosities are 14 to >330 times weaker than expected for typical quasars. Our results from a pilot NuSTAR study of two low-redshift BAL quasars, a Chandra stacking analysis of a sample of high-redshift BAL quasars, and a NuSTAR spectral analysis of the local BAL quasar Mrk 231 have already suggested the existence of intrinsically X-ray weak BAL quasars, i.e., quasars not emitting X-rays at the level expected from their optical/UV emission. The aim of the current program is to extend the search for such extraordinary objects. Three of the six new targets are weakly detected by NuSTAR with <45 counts in the 3-24 keV band, and the other three are not detected. The hard X-ray (8-24 keV) weakness observed by NuSTAR requires Compton-thick absorption if these objects have nominal underlying X-ray emission. However, a soft stacked effective photon index ({\Gamma}~1.8) for this sample disfavors Compton-thick absorption in general. The uniform hard X-ray weakness observed by NuSTAR for this and the pilot samples selected with <10 keV weakness also suggests that the X-ray weakness is intrinsic in at least some of the targets. We conclude that the NuSTAR observations have likely discovered a significant population (>33%) of intrinsically X-ray weak objects among the BAL quasars with significantly weak <10 keV emission. We suggest that intrinsically X-ray weak quasars might be preferentially observed as BAL quasars., Comment: 11 pages, 4 figures. ApJ in press

We present the results of a new analysis of the X-ray selected Active Galactic Nuclei (AGN) population in the vicinity of 135 of the most massive galaxy clusters in the redshift range of 0.2 < z < 0.9 observed with Chandra. With a sample of more than 11,000 X-ray point sources, we are able to measure, for the first time, evidence for evolution in the cluster AGN population beyond the expected evolution of field AGN. Our analysis shows that overall number density of cluster AGN scales with the cluster mass as $\sim M_{500}^{-1.2}$. There is no evidence for the overall number density of cluster member X-ray AGN depending on the cluster redshift in a manner different than field AGN, nor there is any evidence that the spatial distribution of cluster AGN (given in units of the cluster overdensity radius r_500) strongly depends on the cluster mass or redshift. The $M^{-1.2 \pm 0.7}$ scaling relation we measure is consistent with theoretical predictions of the galaxy merger rate in clusters, which is expected to scale with the cluster velocity dispersion, $\sigma$, as $ \sim \sigma^{-3}$ or $\sim M^{-1}$. This consistency suggests that AGN in clusters may be predominantly triggered by galaxy mergers, a result that is further corroborated by visual inspection of Hubble images for 23 spectroscopically confirmed cluster member AGN in our sample. A merger-driven scenario for the triggering of X-ray AGN is not strongly favored by studies of field galaxies, however, suggesting that different mechanisms may be primarily responsible for the triggering of cluster and field X-ray AGN., Comment: 21 Pages, 8 figures, 5 tables. Submitted to MNRAS. Comments are welcome, and please request Steven Ehlert for higher resolution figures

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., Comment: ApJ, in press

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., Comment: Accepted for publication in ApJ

We present a spectral analysis of four coordinated NuSTAR+XMM-Newton observations of the Seyfert galaxy NGC 1365. These exhibit an extreme level of spectral variability, which is primarily due to variable line-of-sight absorption, revealing relatively unobscured states in this source for the first time. Despite the diverse range of absorption states, each of the observations displays the same characteristic signatures of relativistic reflection from the inner accretion disk. Through time-resolved spectroscopy we find that the strength of the relativistic iron line and the Compton reflection hump relative to the intrinsic continuum are well correlated, as expected if they are two aspects of the same broadband reflection spectrum. We apply self-consistent disk reflection models to these time-resolved spectra in order to constrain the inner disk parameters, allowing for variable, partially covering absorption to account for the vastly different absorption states observed. Each of the four observations is treated independently to test the consistency of the results obtained for the black hole spin and the disk inclination, which should not vary on observable timescales. We find both the spin and the inclination determined from the reflection spectrum to be consistent, confirming NGC 1365 hosts a rapidly rotating black hole; in all cases the dimensionless spin parameter is constrained to be a* > 0.97 (at 90% statistical confidence or better)., Comment: 13 pages, 6 figures, accepted for publication in ApJ

We report on a NuSTAR and XMM-Newton program that has observed a sample of three extremely luminous, heavily obscured WISE-selected AGN at z~2 in a broad X-ray band (0.1 - 79 keV). The parent sample, selected to be faint or undetected in the WISE 3.4um (W1) and 4.6um (W2) bands but bright at 12um (W3) and 22um (W4), are extremely rare, with only ~1000 so-called W1W2-dropouts across the extragalactic sky. Optical spectroscopy reveals typical redshifts of z~2 for this population, implying rest-frame mid-IR luminosities of L(6um)~6e46 erg/s and bolometric luminosities that can exceed L(bol)~1e14 L(sun). The corresponding intrinsic, unobscured hard X-ray luminosities are L(2-10)~4e45 erg/s for typical quasar templates. These are amongst the most luminous AGN known, though the optical spectra rarely show evidence of a broad-line region and the selection criteria imply heavy obscuration even at rest-frame 1.5um. We designed our X-ray observations to obtain robust detections for gas column densities N(H)<1e24 /cm2. In fact, the sources prove to be fainter than these predictions. Two of the sources were observed by both NuSTAR and XMM-Newton, with neither being detected by NuSTAR and one being faintly detected by XMM-Newton. A third source was observed only with XMM-Newton, yielding a faint detection. The X-ray data require gas column densities N(H)>1e24 /cm2, implying the sources are extremely obscured, consistent with Compton-thick, luminous quasars. The discovery of a significant population of heavily obscured, extremely luminous AGN does not conform to the standard paradigm of a receding torus, in which more luminous quasars are less likely to be obscured. If a larger sample conforms with this finding, then this suggests an additional source of obscuration for these extreme sources., Comment: 7 pages; submitted to ApJ

We report NuSTAR observations of NuSTAR J033202-2746.8, a heavily obscured, radio-loud quasar detected in the Extended Chandra Deep Field-South, the deepest layer of the NuSTAR extragalactic survey (~400 ks, at its deepest). NuSTAR J033202-2746.8 is reliably detected by NuSTAR only at E>8 keV and has a very flat spectral slope in the NuSTAR energy band (Gamma=0.55^{+0.62}_{-0.64}; 3-30 keV). Combining the NuSTAR data with extremely deep observations by Chandra and XMM-Newton (4 Ms and 3 Ms, respectively), we constrain the broad-band X-ray spectrum of NuSTAR J033202-2746.8, indicating that this source is a heavily obscured quasar (N_H=5.6^{+0.9}_{-0.8}x10^23 cm^-2) with luminosity L_{10-40 keV}~6.4x10^44 erg s^-1. Although existing optical and near-infrared (near-IR) data, as well as follow-up spectroscopy with the Keck and VLT telescopes, failed to provide a secure redshift identification for NuSTAR J033202-2746.8, we reliably constrain the redshift z=2.00+/-0.04 from the X-ray spectral features (primarily from the iron K edge). The NuSTAR spectrum shows a significant reflection component (R=0.55^{+0.44}_{-0.37}), which was not constrained by previous analyses of Chandra and XMM-Newton data alone. The measured reflection fraction is higher than the R~0 typically observed in bright radio-loud quasars such as NuSTAR J033202-2746.8, which has L_{1.4 GHz}~10^27 W Hz^-1. Constraining the spectral shape of AGN, including bright quasars, is very important for understanding the AGN population, and can have a strong impact on the modeling of the X-ray background. Our results show the importance of NuSTAR in investigating the broad-band spectral properties of quasars out to high redshift., Comment: 13 pages, 5 figures; accepted for publication in ApJ

We present high-energy (3--30 keV) {\it NuSTAR} observations of the nearest quasar, the ultraluminous infrared galaxy (ULIRG) Markarian 231 (Mrk 231), supplemented with new and simultaneous low-energy (0.5--8 keV) data from {\it Chandra}. The source was detected, though at much fainter levels than previously reported, likely due to contamination in the large apertures of previous non-focusing hard X-ray telescopes. The full band (0.5--30 keV) X-ray spectrum suggests the active galactic nucleus (AGN) in Mrk 231 is absorbed by a patchy and Compton-thin (N$_{\rm H} \sim1.2^{+0.3}_{-0.3}\times10^{23}$ cm$^{-2}$) column. The intrinsic X-ray luminosity (L$_{\rm 0.5-30 keV}\sim1.0\times10^{43}$ erg s$^{1}$) is extremely weak relative to the bolometric luminosity where the 2--10 keV to bolometric luminosity ratio is $\sim$0.03% compared to the typical values of 2--15%. Additionally, Mrk 231 has a low X-ray-to-optical power law slope ($\alpha_{\rm OX}\sim-1.7$). It is a local example of a low-ionization broad absorption line (LoBAL) quasar that is intrinsically X-ray weak. The weak ionizing continuum may explain the lack of mid-infrared [O IV], [Ne V], and [Ne VI] fine-structure emission lines which are present in sources with otherwise similar AGN properties. We argue that the intrinsic X-ray weakness may be a result of the super-Eddington accretion occurring in the nucleus of this ULIRG, and may also be naturally related to the powerful wind event seen in Mrk 231, a merger remnant escaping from its dusty cocoon., Comment: 15 pages in ApJ style including 6 figures and 1 table. Accepted for publication in ApJ

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., Comment: 10 pages, 6 figures, 3 tables. Accepted for publication in ApJ

We report significant anisotropies in the projected two-dimensional (2D) spatial distributions of Globular Clusters (GCs) of the giant Virgo elliptical galaxy NGC4649 (M60). Similar features are found in the 2D distribution of low-mass X-ray binaries (LMXBs), both associated with GCs and in the stellar field. Deviations from azimuthal symmetry suggest an arc-like excess of GCs extending north at 4-15 kpc galactocentric radii in the eastern side of major axis of NGC4649. This feature is more prominent for red GCs, but still persists in the 2D distribution of blue GCs. High and low luminosity GCs also show some segregation along this arc, with high-luminosity GCs preferentially located in the southern end and low-luminosity GCs in the northern section of the arc. GC-LMXBs follow the anisotropy of red-GCs, where most of them reside; however, a significant overdensity of (high-luminosity) field LMXBs is present to the south of the GC arc. These results suggest that NGC4649 has experienced mergers and/or multiple accretions of less massive satellite galaxies during its evolution, of which the GCs in the arc may be the fossil remnant. We speculate that the observed anisotropy in the field LMXB spatial distribution indicates that these X-ray binaries may be the remnants of a star formation event connected with the merger, or maybe be ejected from the parent red GCs, if the bulk motion of these clusters is significantly affected by dynamical friction. We also detect a luminosity enhancement in the X-ray source population of the companion spiral galaxy NGC4647. We suggest that these may be younger high mass X-ray binaries formed as a result of the tidal interaction of this galaxy with NGC4649., Comment: 15 pages, 11 color figures, accepted for publication in ApJ