Source author record

D. Stern

D. Stern appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

175works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

175 published item(s)

preprint2025arXiv

Unveiling the nature of the Einstein Probe transient EP 241021a

We present a multi-wavelength analysis of the fast X-ray transient EP 241021a, discovered by the Wide-field X-ray Telescope aboard the \emph{Einstein Probe} satellite on 2024 October 21. The event was not detected in gamma-rays. Follow-up observations from $\sim$1.5 to 100 days post-trigger were obtained across X-ray, UV, optical, near-infrared, and radio bands with ground- and space-based facilities. The redshift is constrained to $z = 0.7485$ from prominent optical spectral features. The optical light curve shows complex evolution: an initial $\sim t^{-0.7}$ decay, followed by a rapid re-brightening peaking at day 7.7 with $\sim t^{-1.7}$ decay, and a third phase peaking near day 19 with $\sim t^{-1.3}$ decay. The spectral energy distribution (SED) and its temporal evolution are consistent with a mix of non-thermal and thermal components. Early optical-to-X-ray spectral indices agree with optically thin synchrotron emission, while steepening of the optical SED after $\sim$20 days indicates either a shift in emission mechanism or the emergence of an additional component. Although broad-lined absorption features are absent, comparisons with type Ic-BL supernovae suggest a SN contribution at late times, suggesting a collapsar origin for EP 241021a. The likely SN in EP 241021a appears to require an additional energy source beyond $^{56}$Ni decay. These results support the view that some fast X-ray transients detected by the \emph{Einstein Probe} arise from massive stellar explosions.

preprint2022arXiv

A repeating fast radio burst associated with a persistent radio source

The dispersive sweep of fast radio bursts (FRBs) has been used to probe the ionized baryon content of the intergalactic medium, which is assumed to dominate the total extragalactic dispersion. While the host galaxy contributions to dispersion measure (DM) appear to be small for most FRBs, in at least one case there is evidence for an extreme magneto-ionic local environment and a compact persistent radio source. Here we report the detection and localization of the repeating FRB 20190520B, which is co-located with a compact, persistent radio source and associated with a dwarf host galaxy of high specific star formation rate at a redshift $z=0.241\pm0.001$. The estimated host galaxy DM $\approx 903^{+72}_{-111}$ pc cm$^{-3}$, nearly an order of magnitude higher than the average of FRB host galaxies, far exceeds the DM contribution of the intergalactic medium. Caution is thus warranted in inferring redshifts for FRBs without accurate host galaxy identifications.

preprint2022arXiv

BASS XXXI: Outflow scaling relations in low redshift X-ray AGN host galaxies with MUSE

Ionised gas kinematics provide crucial evidence of the impact that active galactic nuclei (AGN) have in regulating star formation in their host galaxies. Although the presence of outflows in AGN host galaxies has been firmly established, the calculation of outflow properties such as mass outflow rates and kinetic energy remains challenging. We present the [OIII]5007 ionised gas outflow properties of 22 z$<$0.1 X-ray AGN, derived from the BAT AGN Spectroscopic Survey using MUSE/VLT. With an average spatial resolution of 1" (0.1-1.2 kpc), the observations resolve the ionised gas clouds down to sub-kiloparsec scales. Resolved maps show that the [OIII] velocity dispersion is, on average, higher in regions ionised by the AGN, compared to star formation. We calculate the instantaneous outflow rates in individual MUSE spaxels by constructing resolved mass outflow rate maps, incorporating variable outflow density and velocity. We compare the instantaneous values with time-averaged outflow rates by placing mock fibres and slits on the MUSE field-of-view, a method often used in the literature. The instantaneous outflow rates (0.2-275 $M_{\odot}$ yr$^{-1}$) tend to be 2 orders of magnitude higher than the time-averaged outflow rates (0.001-40 $M_{\odot}$ yr$^{-1}$). The outflow rates correlate with the AGN bolometric luminosity ($L_{\rm bol}\sim$ 10$^{42.71}$-10$^{45.62}$ erg/s) but we find no correlations with black hole mass (10$^{6.1}$-10$^{8.9}$ M$_{\odot}$), Eddington ratio (0.002-1.1) and radio luminosity (10$^{21}$-10$^{26}$ W/Hz). We find the median coupling between the kinetic energy and $L_{\rm bol}$ to be 1%, consistent with the theoretical predictions for an AGN-driven outflow.

preprint2022arXiv

Euclid Preparation: XIV. The Complete Calibration of the Color-Redshift Relation (C3R2) Survey: Data Release 3

The Complete Calibration of the Color-Redshift Relation (C3R2) survey is obtaining spectroscopic redshifts in order to map the relation between galaxy color and redshift to a depth of i ~ 24.5 (AB). The primary goal is to enable sufficiently accurate photometric redshifts for Stage IV dark energy projects, particularly Euclid and the Roman Space Telescope, which are designed to constrain cosmological parameters through weak lensing. We present 676 new high-confidence spectroscopic redshifts obtained by the C3R2 survey in the 2017B-2019B semesters using the DEIMOS, LRIS, and MOSFIRE multi-object spectrographs on the Keck telescopes. Combined with the 4454 redshifts previously published by this project, the C3R2 survey has now obtained and published 5130 high-quality galaxy spectra and redshifts. If we restrict consideration to only the 0.2 < z(phot) < 2.6 range of interest for the Euclid cosmological goals, then with the current data release C3R2 has increased the spectroscopic redshift coverage of the Euclid color space from 51% (as reported by Masters et al. 2015) to the current 91%. Once completed and combined with extensive data collected by other spectroscopic surveys, C3R2 should provide the spectroscopic calibration set needed to enable photometric redshifts to meet the cosmology requirements for Euclid, and make significant headway toward solving the problem for Roman.

preprint2022arXiv

Euclid preparation: XX. The Complete Calibration of the Color-Redshift Relation survey: LBT observations and data release

The Complete Calibration of the Color-Redshift Relation survey (C3R2) is a spectroscopic programme designed to empirically calibrate the galaxy color-redshift relation to the Euclid depth (I_E=24.5), a key ingredient for the success of Stage IV dark energy projects based on weak lensing cosmology. A spectroscopic calibration sample as representative as possible of the galaxies in the Euclid weak lensing sample is being collected, selecting galaxies from a self-organizing map (SOM) representation of the galaxy color space. Here, we present the results of a near-infrared H- and K-bands spectroscopic campaign carried out using the LUCI instruments at the LBT. For a total of 251 galaxies, we present new highly-reliable redshifts in the 1.3<= z <=1.7 and 2<= z<=2.7 ranges. The newly-determined redshifts populate 49 SOM cells which previously contained no spectroscopic measurements and almost double the occupation numbers of an additional 153 SOM cells. A final optical ground-based observational effort is needed to calibrate the missing cells in particular in the redshift range 1.7<= z<=2.7 that lack spectroscopic calibration. In the end, Euclid itself will deliver telluric-free NIR spectra that can complete the calibration.

preprint2022arXiv

Extreme X-ray Reflection in the Nucleus of the Seyfert Galaxy NGC 5033

NGC 5033 is an intriguing Seyfert galaxy because its sub-classification may change with time, and because optical and sub-mm observations find that the massive black hole does not sit at the dynamical center of the galaxy, pointing to a past merger. We obtained a new optical spectrum of NGC 5033 using the 200'' Hale telescope at Palomar that clearly reveals a broad H$β$ line (FWHM$=5400\pm 300~{\rm km}~{\rm s}^{-1}$). This signals a clear view of the optical broad line region (BLR) and requires Seyfert-1.5 designation. Some spectra obtained in the past suggest a Seyfert-1.9 classification, potentially signaling a variable or "changing-look" geometry. Our analysis of a 2019 Chandra spectrum of the massive black hole reveals very little obscuration, also suggesting a clean view of the central engine. However, the narrow Fe~K$α$ emission line is measured to have an equivalent with of EW$=460^{+100}_{-90}$~eV. This value is extremely high compared to typical values in unobscured AGN. Indeed, the line is persistently strong in NGC 5033: the line equivalent width in a 2002 XMM-Newton snapshot is EW$=250^{+40}_{-40}$~eV, similar to the EW$=290^{+100}_{-100}$~eV equivalent width measured using ASCA in 1999. These results can likely be explained through a combination of an unusually high covering factor for reflection, and fluxes that are seen out of phase owing to light travel times. We examine the possibility that NGC 5033 may strengthen evidence for the X-ray Baldwin effect.

preprint2022arXiv

Identification of an X-ray Pulsar in the BeXRB system IGR J18219$-$1347

We report on observations of the candidate Be/X-ray binary IGR J18219$-$1347 with \textit{Swift}/XRT, \textit{NuSTAR}, and \textit{NICER} during Type-I outbursts in March and June 2020. Our timing analysis revealed the spin period of a neutron star with $P_\textrm{spin}=52.46$ s. This periodicity, combined with the known orbital period of $72.4$ d, indicates that the system is a BeXRB. Furthermore, by comparing the infrared counterpart's spectral energy distribution to known BeXRBs, we confirm this classification and set a distance of approximately $10-15$ kpc for the source. The source's broadband X-ray spectrum ($1.5-50$ keV) is described by an absorbed power-law with photon index $Γ$\,$\sim$\,$0.5$ and cutoff energy at $\sim$\,$13$ keV.

preprint2022arXiv

NuSTAR Observations of Intrinsically X-ray Weak Quasar Candidates: An Obscuration-Only Scenario

We utilize recent NuSTAR observations (co-added depth $\approx55$-120 ks) of PG $1001+054$, PG $1254+047$, and PHL 1811 to constrain their hard X-ray ($\gtrsim5$ keV) weakness and spectral shapes, and thus to investigate the nature of their extreme X-ray weakness. These quasars showed very weak soft X-ray emission, and they were proposed to be intrinsically X-ray weak, with the X-ray coronae producing weak continuum emission relative to their optical/UV emission. However, the new observations suggest an alternative explanation. The NuSTAR 3-24 keV spectral shapes for PG $1001+054$ and PHL 1811 are likely flat (effective power-law photon indices $Γ_{\rm eff}=1.0^{+0.5}_{-0.6}$ and $Γ_{\rm eff}=1.4^{+0.8}_{-0.7}$, respectively), while the shape is nominal for PG $1254+047$ ($Γ_{\rm eff}=1.8\pm0.3$). PG $1001+054$ and PHL 1811 are significantly weak at hard X-ray energies (by factors of $\approx26$-74 at rest-frame 8 keV) compared to the expectations from their optical/UV emission, while PG $1254+047$ is only hard X-ray weak by a factor of $\approx3$. We suggest that X-ray obscuration is present in all three quasars. We propose that, as an alternative to the intrinsic X-ray weakness + X-ray obscuration scenario, the soft and hard X-ray weakness of these quasars can be uniformly explained under an obscuration-only scenario. This model provides adequate descriptions of the multi-epoch soft and hard X-ray data of these quasars, with variable column density and leaked fraction of the partial-covering absorber. We suggest that the absorber is the clumpy dust-free wind launched from the accretion disk. These quasars probably have super-Eddington accretion rates that drive powerful and high-density winds.

preprint2021arXiv

COSMOS2020: A panchromatic view of the Universe to $z\sim10$ from two complementary catalogs

The Cosmic Evolution Survey (COSMOS) has become a cornerstone of extragalactic astronomy. Since the last public catalog in 2015, a wealth of new imaging and spectroscopic data has been collected in the COSMOS field. This paper describes the collection, processing, and analysis of this new imaging data to produce a new reference photometric redshift catalog. Source detection and multi-wavelength photometry is performed for 1.7 million sources across the $2\,\mathrm{deg}^{2}$ of the COSMOS field, $\sim$966,000 of which are measured with all available broad-band data using both traditional aperture photometric methods and a new profile-fitting photometric extraction tool, The Farmer, which we have developed. A detailed comparison of the two resulting photometric catalogs is presented. Photometric redshifts are computed for all sources in each catalog utilizing two independent photometric redshift codes. Finally, a comparison is made between the performance of the photometric methodologies and of the redshift codes to demonstrate an exceptional degree of self-consistency in the resulting photometric redshifts. The $i<21$ sources have sub-percent photometric redshift accuracy and even the faintest sources at $25<i<27$ reach a precision of $5\,\%$. Finally, these results are discussed in the context of previous, current, and future surveys in the COSMOS field. Compared to COSMOS2015, reaches the same photometric redshift precision at almost one magnitude deeper. Both photometric catalogs and their photometric redshift solutions and physical parameters will be made available through the usual astronomical archive systems (ESO Phase 3, IPAC IRSA, and CDS).

preprint2021arXiv

Euclid preparation: XVIII. Cosmic Dawn Survey. Spitzer observations of the Euclid deep fields and calibration fields

We present a new infrared survey covering the three Euclid deep fields and four other Euclid calibration fields using Spitzer's Infrared Array Camera (IRAC). We have combined these new observations with all relevant IRAC archival data of these fields in order to produce the deepest possible mosaics of these regions. In total, these observations represent nearly 11% of the total Spitzer mission time. The resulting mosaics cover a total of approximately 71.5deg$^2$ in the 3.6 and 4.5um bands, and approximately 21.8deg$^2$ in the 5.8 and 8um bands. They reach at least 24 AB magnitude (measured to sigma, in a 2.5 arcsec aperture) in the 3.6um band and up to ~ 5 mag deeper in the deepest regions. The astrometry is tied to the Gaia astrometric reference system, and the typical astrometric uncertainty for sources with 16<[3.6]<19 is <0.15 arcsec. The photometric calibration is in excellent agreement with previous WISE measurements. We have extracted source number counts from the 3.6um band mosaics and they are in excellent agreement with previous measurements. Given that the Spitzer Space Telescope has now been decommissioned these mosaics are likely to be the definitive reduction of these IRAC data. This survey therefore represents an essential first step in assembling multi-wavelength data on the Euclid deep fields which are set to become some of the premier fields for extragalactic astronomy in the 2020s.

preprint2021arXiv

Mapping the intracluster medium temperature in the era of NIKA2 and MUSTANG-2

We present preliminary results from an on-going program that aims at mapping the intracluster medium (ICM) temperature of high redshift galaxy clusters from the MaDCoWS sample using a joint analysis of shallow X-ray data obtained by $Chandra$ and high angular resolution Sunyaev-Zel'dovich (SZ) observations realized with the NIKA2 and MUSTANG-2 cameras. We also present preliminary results from an on-going Open Time program within the NIKA2 collaboration that aims at mapping the ICM temperature of a galaxy cluster at $z=0.45$ from the resolved detection of the relativistic corrections to the SZ spectrum. These studies demonstrate how high angular resolution SZ observations will play a major role in the coming decade to push the investigation of ICM dynamics and non-gravitational processes to high redshift before the next generation X-ray observatories come into play.

preprint2021arXiv

Starfall: A heavy rain of stars in 'turning on' AGN

As active galactic nuclei (AGN) `turn on', some stars end up embedded in accretion disks around supermassive black holes (SMBHs) on retrograde orbits. Such stars experience strong headwinds, aerodynamic drag, ablation and orbital evolution on short timescales. Loss of orbital angular momentum in the first $\sim 0.1$~Myr of an AGN leads to a heavy rain of stars (`starfall') into the inner disk and onto the SMBH. A large AGN loss cone ($θ_{\rm AGN,lc}$) can result from binary scatterings in the inner disk and yield tidal disruption events (TDEs). Signatures of starfall include optical/UV flares that rise in luminosity over time, particularly in the inner disk. If the SMBH mass is $M_{\rm SMBH} \ge 10^{8}M_{\odot}$, flares truncate abruptly and the star is swallowed. If $M_{\rm SMBH}<10^{8}M_{\odot}$, and if the infalling orbit lies within $θ_{\rm AGN,lc}$, the flare is followed by a TDE which can be prograde or retrograde relative to the AGN inner disk. Retrograde AGN TDEs are over-luminous and short-lived as in-plane ejecta collide with the inner disk and a lower AGN state follows. Prograde AGN TDEs add angular momentum to inner disk gas and so start off looking like regular TDEs but are followed by an AGN high state. Searches for such flare signatures test models of AGN `turn on', SMBH mass, as well as disk properties and the embedded population.

preprint2021arXiv

The black hole masses of extremely luminous radio-WISE selected galaxies

We present near-IR photometry and spectroscopy of 30 extremely luminous radio and mid-IR selected galaxies. With bolometric luminosities exceeding $\sim10^{13}$ $\rm{L_{\odot}}$ and redshifts ranging from $z = 0.880-2.853$, we use VLT instruments X-shooter and ISAAC to investigate this unique population of galaxies. Broad multi-component emission lines are detected in 18 galaxies and we measure the near-IR lines $\rm{H\rmβ}$, $\text{[OIII]}\rmλ\rmλ4959,5007$ and $\rm{H\rmα}$ in six, 15 and 13 galaxies respectively, with 10 $\rm{Lyα}$ and five CIV lines additionally detected in the UVB arm. We use the broad $\text{[OIII]}\rmλ5007$ emission lines as a proxy for the bolometric AGN luminosity, and derive lower limits to supermassive black hole masses of $10^{7.9}$-$10^{9.4}$ $\text{M}_{\odot}$ with expectations of corresponding host masses of $10^{10.4}$-$10^{12.0}$ $\text{M}_{\odot}$. We measure $\rmλ_{Edd}$ > 1 for eight of these sources at a $2σ$ significance. Near-IR photometry and SED fitting are used to compare stellar masses directly. We detect both Balmer lines in five galaxies and use these to infer a mean visual extinction of $A_{V}$ = 2.68 mag. Due to non-detections and uncertainties in our $\rm{H\rmβ}$ emission line measurements, we simulate a broad $\rm{H\rmβ}$ line of FWHM = 1480 $\rm{kms^{-1}}$ to estimate extinction for all sources with measured $\rm{H\rmα}$ emission. We then use this to infer a mean $A_{V}=3.62$ mag, demonstrating the highly-obscured nature of these galaxies, with the consequence of increasing our estimates of black-hole masses by an 0.5 orders of magnitude in the most extreme and obscured cases.

preprint2021arXiv

TOI-2257 b: A highly eccentric long-period sub-Neptune transiting a nearby M dwarf

Thanks to the relative ease of finding and characterizing small planets around M dwarf stars, these objects have become cornerstones in the field of exoplanet studies. The current paucity of planets in long-period orbits around M dwarfs make such objects particularly compelling as they provide clues about the formation and evolution of these systems. In this study, we present the discovery of TOI-2257 b (TIC 198485881), a long-period (35 d) sub-Neptune orbiting an M3 star at 57.8pc. Its transit depth is about 0.4%, large enough to be detected with medium-size, ground-based telescopes. The long transit duration suggests the planet is in a highly eccentric orbit ($e \sim 0.5$), which would make it the most eccentric planet that is known to be transiting an M-dwarf star. We combined TESS and ground-based data obtained with the 1.0-m SAINT-EX, 0.60-m TRAPPIST-North and 1.2-m FLWO telescopes to find a planetary size of 2.2 $R_{\oplus}$ and an orbital period of 35.19 days. In addition, we make use of archival data, high-resolution imaging, and vetting packages to support our planetary interpretation. With its long period and high eccentricity, TOI-2257 b falls in a novel slice of parameter space. Despite the planet's low equilibrium temperature ($\sim$ 256 K), its host star's small size ($R_* = 0.311 \pm{0.015}$) and relative infrared brightness (K$_{mag}$ = 10.7) make it a suitable candidate for atmospheric exploration via transmission spectroscopy.

preprint2020arXiv

$NuSTAR$ Observations of Four Nearby X-ray Faint AGN: Low Luminosity or Heavy Obscuration?

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

preprint2020arXiv

Candidate Electromagnetic Counterpart to the Binary Black Hole Merger Gravitational Wave Event S190521g

We report the first plausible optical electromagnetic (EM) counterpart to a (candidate) binary black hole (BBH) merger. Detected by the Zwicky Transient Facility (ZTF), the EM flare is consistent with expectations for a kicked BBH merger in the accretion disk of an active galactic nucleus (AGN), and is unlikely ($<O(0.01\%$)) due to intrinsic variability of this source. The lack of color evolution implies that it is not a supernovae and instead is strongly suggestive of a constant temperature shock. Other false-positive events, such as microlensing or a tidal disruption event, are ruled out or constrained to be $<O(0.1\%$). If the flare is associated with S190521g, we find plausible values of: total mass $ M_{\rm BBH} \sim 100 M_{\odot}$, kick velocity $v_k \sim 200\, {\rm km}\, {\rm s}^{-1}$ at $θ\sim 60^{\circ}$ in a disk with aspect ratio $H/a \sim 0.01$ (i.e., disk height $H$ at radius $a$) and gas density $ρ\sim 10^{-10}\, {\rm g}\, {\rm cm}^{-3}$. The merger could have occurred at a disk migration trap ($a \sim 700\, r_{g}$; $r_g \equiv G M_{\rm SMBH} / c^2$, where $M_{\rm SMBH}$ is the mass of the AGN supermassive black hole). The combination of parameters implies a significant spin for at least one of the black holes in S190521g. The timing of our spectroscopy prevents useful constraints on broad-line asymmetry due to an off-center flare. We predict a repeat flare in this source due to a re-encountering with the disk in $\sim 1.6\, {\rm yr}\, (M_{\rm SMBH}/10^{8}M_{\odot})\, (a/10^{3}r_{g})^{3/2}$.

preprint2020arXiv

Cold Molecular Gas and Free-Free Emission from Hot, Dust-Obscured Galaxies at $z$~3

We report on observations of redshifted CO(1-0) line emission and observed-frame $\rm\sim$ 30GHz radio continuum emission from five ultra-luminous, mid-IR selected hot, Dust-Obscured Galaxies (Hot DOGs) at $z\rm\gtrsim$ 3 using the Karl G. Jansky Very Large Array. We detect CO(1-0) line emission in all five Hot DOGs, with one of them at high signal to noise. We analyse FIR-radio spectral energy distributions, including dust, free-free and synchrotron emission for the galaxies. We find that most of the 115 GHz rest-frame continuum is mostly due to synchrotron or free-free emission, with only a potentially small contribution from thermal emission. We see a deficit in the rest-frame 115 GHz continuum emission compared to dusty star-forming galaxies (DSFGs) and sub-millimetre galaxies (SMGs) at high redshift, suggesting that Hot DOGs do not have similar cold gas reserves compared with star-forming galaxies. One target, W2305-0039, is detected in the FIRST 1.4 GHz survey, and is likely to possess compact radio jets. We compare to the FIR-radio correlation, and find that at least half of the Hot DOGs in our sample are radio-quiet with respect to normal galaxies. These findings suggest that Hot DOGs have comparably less cold molecular gas than star-forming galaxies at lower, $z\rm\sim$ 2 redshifts, and are dominated by powerful, yet radio-quiet AGN.

preprint2020arXiv

Constraining X-ray reflection in the low-luminosity AGN NGC 3718 using NuSTAR and XMM--Newton

One distinctive feature of low-luminosity active galactic nuclei (LLAGN) is the relatively weak reflection features they may display in the X-ray spectrum, which can result from the disappearance of the torus with decreasing accretion rates. Some material, however, must surround the active nucleus, i.e., the accretion flow itself and, possibly, a flattened-out or thinned torus. In this work, we study whether reflection is indeed absent or undetectable due to its intrinsically weak features together with the low statistics inherent to LLAGN. Here we focus on NGC 3718 ($L/L_{\rm Edd}\sim10^{-5}$) combining observations from XMM--Newton and the deepest to date NuSTAR (0.5--79 keV) spectrum of a LLAGN, to constrain potential reflectors, and analyze how the fitted coronal parameters depend on the reflection model. We test models representing both an accretion disc (Relxill) and a torus-like (MYTorus and Borus) neutral reflector. From a statistical point of view, reflection is not required, but its inclusion allows to place strong constraints on the geometry and physical features of the surroundings: both neutral reflectors (torus) tested should be Compton thin ($N_H<10^{23.2}$cm$^{-2}$) and preferentially cover a large fraction of the sky. If the reflected light instead arises from an ionized reflector, a highly ionized case is preferred. These models produce an intrinsic power-law spectral index in the range [1.81--1.87], where the torus models result in steeper slopes. The cut-off energy of the power-law emission also changes with the inclusion of reflection models, resulting in constrained values for the disc reflectors and unconstrained values for torus reflectors.

preprint2020arXiv

Euclid preparation: VIII. The Complete Calibration of the Colour-Redshift Relation survey: VLT/KMOS observations and data release

The Complete Calibration of the Colour-Redshift Relation survey (C3R2) is a spectroscopic effort involving ESO and Keck facilities designed to empirically calibrate the galaxy colour-redshift relation - P(z|C) to the Euclid depth (i_AB=24.5) and is intimately linked to upcoming Stage IV dark energy missions based on weak lensing cosmology. The aim is to build a spectroscopic calibration sample that is as representative as possible of the galaxies of the Euclid weak lensing sample. In order to minimise the number of spectroscopic observations to fill the gaps in current knowledge of the P(z|C), self-organising map (SOM) representations of the galaxy colour space have been constructed. Here we present the first results of an ESO@ VLT Large Programme approved in the context of C3R2, which makes use of the two VLT optical and near-infrared multi-object spectrographs, FORS2 and KMOS. This paper focuses on high-quality spectroscopic redshifts of high-z galaxies observed with the KMOS spectrograph in the H- and K-bands. A total of 424 highly-reliable z are measured in the 1.3<=z<=2.5 range, with total success rates of 60.7% in the H-band and 32.8% in the K-band. The newly determined z fill 55% of high and 35% of lower priority empty SOM grid cells. We measured Halpha fluxes in a 1."2 radius aperture from the spectra of the spectroscopically confirmed galaxies and converted them into star formation rates. In addition, we performed an SED fitting analysis on the same sample in order to derive stellar masses, E(B-V), total magnitudes, and SFRs. We combine the results obtained from the spectra with those derived via SED fitting, and we show that the spectroscopic failures come from either weakly star-forming galaxies (at z<1.7, i.e. in the H-band) or low S/N spectra (in the K-band) of z>2 galaxies.

preprint2020arXiv

Hot Dust-Obscured Galaxies with Excess Blue Light

Hot Dust-Obscured Galaxies (Hot DOGs) are among the most luminous galaxies in the Universe. Powered by highly obscured, possibly Compton-thick, active galactic nuclei (AGNs), Hot DOGs are characterized by SEDs that are very red in the mid-IR yet dominated by the host galaxy stellar emission in the UV and optical. An earlier study identified a sub-sample of Hot DOGs with significantly enhanced UV emission. One target, W0204-0506, was studied in detail and, based on Chandra observations, it was concluded that the enhanced emission was most likely due to either extreme unobscured star-formation (${\rm SFR}>1000~M_{\odot}~\rm yr^{-1}$) or to light from the highly obscured AGN scattered by gas or dust into our line of sight. Here, we present a follow-up study of W0204-0506 as well as two more Hot DOGs with excess UV emission. For the two new objects we obtained Chandra/ACIS-S observations, and for all three targets we obtained HST/WFC3 F555W and F160W imaging. We conclude that the excess UV emission is primarily dominated by light from the central highly obscured, hyper-luminous AGN that has been scattered into our line of sight. We cannot rule out, however, that star-formation may significantly contribute to the UV excess of W0204-0506.

preprint2020arXiv

Local AGN Survey (LASr): I. Galaxy sample, infrared colour selection and predictions for AGN within 100 Mpc

In order to answer some of the major open questions in the fields of supermassive black hole (SMBH) and galaxy evolution, a complete census of SMBH growth, i.e., active galactic nuclei (AGN), is required. Thanks to deep all-sky surveys, such as those by the Wide-field Infrared Survey Explorer (WISE) and the Spectrum-Roentgen-Gamma (SRG) missions, this task is now becoming feasible in the nearby Universe. Here, we present a new survey, the Local AGN Survey (LASr), with the goal of identifying AGN unbiased against obscuration and determining the intrinsic Compton-thick (CT) fraction. First, we construct the most complete all-sky sample of galaxies within 100 Mpc from astronomical databases (90% completeness for log M_star/M_sun ~9.4), four times deeper than the current local galaxy reference, the Two Micron All-Sky Survey Redshift Survey (2MRS), which turns out to miss 20% of known luminous AGN. These 49k galaxies serve as parent sample for LASr, called LASr-GPS. It contains 4.3k already known AGN, >=82% of these are estimated to have L(12um) < 10^42.3 erg/s, i.e., are low-luminosity AGN. As a first method for identifying Seyfert-like AGN, we use WISE-based infrared colours, finding 221 galaxies at L(12um) >= 10^42.3 erg/s to host an AGN at 90% reliability, This includes 61 new AGN candidates and implies and optical type 2 fraction of 50 to 71%. We quantify the efficiency of this technique and estimate the total number of AGN with L_int(2-10 keV) >= 10^42 erg/s in the volume to be 362^+145_-116 (8.6^+3.5_-2.8 10^-5 Mpc^-3). X-ray brightness estimates indicate the CT fraction to be 40-55% to explain the Swift non-detections of the infrared selected objects. One third of the AGN within 100 Mpc remain to be identified and we discuss the prospects for the eROSITA all-sky survey.

preprint2020arXiv

The Unusual Broadband X-ray Spectral Variability of NGC 1313 X-1 seen with XMM$-$Newton, Chandra and NuSTAR

We present results from the major coordinated X-ray observing program on the ULX NGC 1313 X-1 performed in 2017, combining $XMM$-$Newton$, $Chandra$ and $NuSTAR$, focusing on the evolution of the broadband ($\sim$0.3-30.0 keV) continuum emission. Clear and unusual spectral variability is observed, but this is markedly suppressed above $\sim$10-15 keV, qualitatively similar to the ULX Holmberg IX X-1. We model the multi-epoch data with two-component accretion disc models designed to approximate super-Eddington accretion, allowing for both a black hole and a neutron star accretor. With regards to the hotter disc component, the data trace out two distinct tracks in the luminosity-temperature plane, with larger emitting radii and lower temperatures seen at higher observed fluxes. Despite this apparent anti-correlation, each of these tracks individually shows a positive luminosity-temperature relation. Both are broadly consistent with $L\propto{T}^{4}$, as expected for blackbody emission with a constant area, and also with $L\propto{T}^{2}$, as may be expected for an advection-dominated disc around a black hole. We consider a variety of possibilities for this unusual behaviour. Scenarios in which the innermost flow is suddenly blocked from view by outer regions of the super-Eddington disc/wind can explain the luminosity-temperature behaviour, but are difficult to reconcile with the lack of strong variability at higher energies, assuming this emission arises from the most compact regions. Instead, we may be seeing evidence for further radial stratification of the accretion flow than is included in the simple models considered, with a combination of winds and advection resulting in the suppressed high-energy variability.

preprint2019arXiv

BAT AGN Spectroscopic Survey -- XIX: Type 1 versus Type 2 AGN dichotomy from the point of view of ionized outflows

We present a detailed study of ionized outflows in a large sample of ~650 hard X-ray detected AGN. Using optical spectroscopy from the BAT AGN Spectroscopic Survey (BASS) we are able to reveal the faint wings of the [OIII] emission lines associated with outflows covering, for the first time, an unexplored range of low AGN bolometric luminosity at low redshift (z~0.05). We test if and how the incidence and velocity of ionized outflow is related to AGN physical parameters: black hole mass, gas column density, Eddington Ratio, [OIII], X-ray, and bolometric luminosities. We find a higher occurrence of ionized outflows in type 1.9 (55%) and type 1 AGN (46%) with respect to type 2 AGN (24%). While outflows in type 2 AGN are evenly balanced between blue and red velocity offsets with respect to the [OIII] narrow component, they are almost exclusively blueshifted in type 1 and type 1.9 AGN. We observe a significant dependence between the outflow occurrence and accretion rate, which becomes relevant at high Eddington ratios (> -1.7). We interpret such behaviour in the framework of covering factor-Eddington ratio dependence. We don't find strong trends of the outflow maximum velocity with AGN physical parameters, as an increase with bolometric luminosity can be only identified when including samples of AGN at high luminosity and high redshift taken from literature.

preprint2019arXiv

Mapping the gas thermodynamic properties of the massive cluster merger MOO J1142$+$1527 at z = 1.2

We present the results of the analysis of the very massive cluster MOO J1142$+$1527 at a redshift $z = 1.2$ based on high angular resolution NIKA2 Sunyaev-Zel'dovich (SZ) and $Chandra$ X-ray data. This multi-wavelength analysis enables us to estimate the shape of the temperature profile with unprecedented precision at this redshift and to obtain a map of the gas entropy distribution averaged along the line of sight. The comparison between the cluster morphological properties observed in the NIKA2 and $Chandra$ maps together with the analysis of the entropy map allows us to conclude that MOO J1142$+$1527 is an on-going merger hosting a cool-core at the position of the X-ray peak. This work demonstrates how the addition of spatially-resolved SZ observations to low signal-to-noise X-ray data can bring valuable insights on the intracluster medium thermodynamic properties at $z>1$.

preprint2019arXiv

The nature of the broadband X-ray variability in the dwarf Seyfert galaxy NGC 4395

We present a flux-resolved X-ray analysis of the dwarf Seyfert 1.8 galaxy NGC 4395, based on three archival $XMM-Newton$ and one archival $NuSTAR$ observations. The source is known to harbor a low mass black hole ($\sim 10^4- {\rm a~ few~}\times 10^{5}~\rm M_\odot$) and shows strong variability in the full X-ray range during these observations. We model the flux-resolved spectra of the source assuming three absorbing layers: neutral, mildly ionized, and highly ionized ($N_{\rm H} \sim 1.6\times 10^{22}-3.4 \times 10^{23}~\rm cm^{-2}$, $\sim 0.8-7.8 \times 10^{22}~\rm cm^{-2}$, and $ 3.8 \times 10^{22}~\rm cm^{-2}$, respectively. The source also shows intrinsic variability by a factor of $\sim 3$, on short timescales, due to changes in the nuclear flux, assumed to be a power law ($Γ= 1.6-1.67$). Our results show a positive correlation between the intrinsic flux and the absorbers' ionization parameter. The covering fraction of the neutral absorber varies during the first $XMM-Newton$ observation, which could explain the pronounced soft X-ray variability. However, the source remains fully covered by this layer during the other two observations, largely suppressing the soft X-ray variability. This suggests an inhomogeneous and layered structure in the broad line region. We also find a difference in the characteristic timescale of the power spectra between different energy ranges and observations. We finally show simulated spectra with $XRISM$, $Athena$, and $eXTP$, which will allow us to characterize the different absorbers, study their dynamics, and will help us identify their locations and sizes.

preprint2019arXiv

Unveiling the merger dynamics of the most massive MaDCoWS cluster at $z = 1.2$ from a multi-wavelength mapping of its intracluster medium properties

The characterization of the Intra-Cluster Medium (ICM) properties of high-redshift galaxy clusters is fundamental to our understanding of large-scale structure formation processes. We present the results of a multi-wavelength analysis of the very massive cluster MOO J1142$+$1527 at a redshift $z = 1.2$ discovered as part of the Massive and Distant Clusters of WISE Survey (MaDCoWS). This analysis is based on high angular resolution $Chandra$ X-ray and NIKA2 Sunyaev-Zel'dovich (SZ) data. Although the X-ray data have only about 1700 counts, we are able to determine the ICM thermodynamic radial profiles, namely temperature, entropy, and hydrostatic mass. These have been obtained with unprecedented precision at this redshift and up to $0.7R_{500}$, thanks to the combination of high-resolution X-ray and SZ data. The comparison between the galaxy distribution mapped in infrared by $Spitzer$ and the morphological properties of the ICM derived from the combined analysis of the $Chandra$ and NIKA2 data leads us to the conclusion that the cluster is an on-going merger. We measure the hydrostatic mass profile of the cluster in four angular sectors centered on the large-scale X-ray centroid. This allows us to estimate a systematic uncertainty on the cluster total mass that characterizes both the impact of the observed deviations from spherical symmetry and of the core dynamics on the mass profile. We further combine the X-ray and SZ data at the pixel level to obtain maps of the temperature and entropy distributions averaged along the line of sight. We find a relatively low entropy core at the position of the X-ray peak and high temperature regions located on its south and west sides. The increase in ICM temperature at the location of the SZ peak is expected given the merger dynamics. (abridged)

preprint2019arXiv

VLT/SINFONI study of black hole growth in high redshift radio-loud quasars from the CARLA survey

We present VLT/SINFONI observations of 35 quasars at 2.1 < z < 3.2, the majority of which were selected from the Clusters Around Radio-Loud AGN (CARLA) survey. CARLA quasars have large CIV-based black hole (BH) masses (M(BH) > 10^9 Msun) and powerful radio emission (P(500MHz) > 27.5 W/Hz). We estimate Ha-based M(BH), finding a scatter of 0.35 dex compared to CIV. We evaluate several recipes for correcting CIV-based masses, which reduce the scatter to 0.24 dex. The radio power of the radio-loud quasars is at most weakly correlated with the interconnected quantities Ha-width, L(5100A) and M(BH), suggesting that it is governed by different physical processes. However, we do find a strong inverse correlation between CIV blueshift and radio power linked to higher Eddington ratios and L(5100A). Under standard assumptions, the BH growth time is longer than the cosmic age for many CARLA quasars, suggesting that they must have experienced more efficient growth in the past. If these BHs were growing from seeds since the epoch of reionization, it is possible that they grew at the Eddington limit like the quasars at z ~ 6-7, and then continued to grow at the reduced rates observed until z ~ 2-3. Finally, we study the relation between M(BH) and environment, finding a weak positive correlation between M(BH) and galaxy density measured by CARLA.

preprint2016arXiv

A 78 Day X-ray Period Detected from NGC 5907 ULX1 by Swift

We report the detection of a $78.1\pm0.5$ day period in the X-ray lightcurve of the extreme ultraluminous X-ray source NGC 5907 ULX1 ($L_{\rm{X,peak}}\sim5\times10^{40}$ erg s$^{-1}$), discovered during an extensive monitoring program with Swift. These periodic variations are strong, with the observed flux changing by a factor of $\sim$3-4 between the peaks and the troughs of the cycle; our simulations suggest that the observed periodicity is detected comfortably in excess of 3$σ$ significance. We discuss possible origins for this X-ray period, but conclude that at the current time we cannot robustly distinguish between orbital and super-orbital variations.

preprint2016arXiv

A growth-rate indicator for Compton-thick active galactic nuclei

Due to their heavily obscured central engines, the growth rate of Compton-thick (CT) active galactic nuclei (AGN) is difficult to measure. A statistically significant correlation between the Eddington ratio, λ$_{Edd}$, and the X-ray power-law index, Γ, observed in unobscured AGN offers an estimate of their growth rate from X-ray spectroscopy (albeit with large scatter). However, since X-rays undergo reprocessing by Compton scattering and photoelectric absorption when the line-of-sight to the central engine is heavily obscured, the recovery of the intrinsic Γ is challenging. Here we study a sample of local, predominantly Compton-thick megamaser AGN, where the black hole mass, and thus Eddington luminosity, are well known. We compile results on X-ray spectral fitting of these sources with sensitive high-energy (E> 10 keV) NuSTAR data, where X-ray torus models which take into account the reprocessing effects have been used to recover the intrinsic Γ values and X-ray luminosities, L$_X$. With a simple bolometric correction to L$_X$ to calculate λ$_{Edd}$, we find a statistically significant correlation between Γ and λ$_{Edd}$ (p = 0.007). A linear fit to the data yields Γ = (0.41$\pm$0.18)log$_{10}$λ$_{Edd}$+(2.38$\pm$ 0.20), which is statistically consistent with results for unobscured AGN. This result implies that torus modeling successfully recovers the intrinsic AGN parameters. Since the megamasers have low-mass black holes (M$_{BH}\approx10^6-10^7$ M$_{sol}$) and are highly inclined, our results extend the Γ-λ$_{Edd}$ relationship to lower masses and argue against strong orientation effects in the corona, in support of AGN unification. Finally this result supports the use of Γ as a growth-rate indicator for accreting black holes, even for Compton-thick AGN.

preprint2016arXiv

A Hard X-ray Study of the Normal Star-Forming Galaxy M83 with NuSTAR

We present results from sensitive, multi-epoch NuSTAR observations of the late-type star-forming galaxy M83 (d=4.6 Mpc), which is the first investigation to spatially resolve the hard (E>10 keV) X-ray emission of this galaxy. The nuclear region and ~ 20 off-nuclear point sources, including a previously discovered ultraluminous X-ray (ULX) source, are detected in our NuSTAR observations. The X-ray hardnesses and luminosities of the majority of the point sources are consistent with hard X-ray sources resolved in the starburst galaxy NGC 253. We infer that the hard X-ray emission is most likely dominated by intermediate accretion state black hole binaries and neutron star low-mass X-ray binaries (Z-sources). We construct the X-ray binary luminosity function (XLF) in the NuSTAR band for an extragalactic environment for the first time. The M83 XLF has a steeper XLF than the X-ray binary XLF in NGC 253, consistent with previous measurements by Chandra at softer X-ray energies. The NuSTAR integrated galaxy spectrum of M83 drops quickly above 10 keV, which is also seen in the starburst galaxies NGC253, NGC 3310 and NGC 3256. The NuSTAR observations constrain any AGN to be either highly obscured or to have an extremely low luminosity of $_{\sim}^<$10$^{38}$ erg/s (10-30 keV), implying it is emitting at a very low Eddington ratio. An X-ray point source consistent with the location of the nuclear star cluster with an X-ray luminosity of a few times 10$^{38}$ erg/s may be a low-luminosity AGN but is more consistent with being an X-ray binary.

preprint2016arXiv

A High Braking Index for a Pulsar

We present a phase-coherent timing solution for PSR J1640-4631, a young 206 ms pulsar using X-ray timing observations taken with NuSTAR. Over this timing campaign, we have measured the braking index of PSR J1640-4631 to be n = 3.15+/-0.03. Using a series of simulations, we argue that this unusually high braking index is not due to timing noise, but is intrinsic to the pulsar's spin-down. We cannot, however, rule out contamination due to an unseen glitch recovery, although the recovery timescale would have to be longer than most yet observed. If this braking index is eventually proven to be stable, it demonstrates that pulsar braking indices greater than 3 are allowed in nature, hence other physical mechanisms such as mass or magnetic quadrupoles are important in pulsar spin-down. We also present a 3-sigma upper limit on the pulsed flux at 1.4 GHz of 0.018 mJy.

preprint2016arXiv

A New Population of Compton-Thick AGN Identified Using the Spectral Curvature Above 10 keV

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.

preprint2016arXiv

An Iron K Component to the Ultrafast Outflow in NGC 1313 X-1

We present the detection of an absorpton feature at $E=8.77^{+0.05}_{-0.06}$ keV in the combined X-ray spectrum of the ultraluminous X-ray source NGC 1313 X-1 observed with XMM-Newton and NuSTAR, significant at the 3$σ$ level. If associated with blueshifted ionized iron, the implied outflow velocity is ~0.2$c$ for Fe XXVI, or ~0.25$c$ for Fe XXV. These velocities are similar to the ultrafast outflow seen in absorption recently discovered in this source at lower energies by XMM-Newton, and we therefore conclude that this is an iron component to the same outflow. Photoionization modeling marginally prefers the Fe XXV solution, but in either case the outflow properties appear to be extreme, potentially supporting a super-Eddington hypothesis for NGC 1313 X-1.

preprint2016arXiv

Broadband X-ray spectral analysis of the Seyfert 1 galaxy GRS 1734-292

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

preprint2016arXiv

Coronal properties of the luminous radio-quiet quasar QSO B2202$-$209

We present an analysis of the joint $XMM$-$Newton$ and $NuSTAR$ observations of the radio-quiet quasar QSO B2202$-$209. Using an optical observation from the Hale Telescope at the Palomar Observatory, we revise the redshift of the source from the previously reported $z=1.77$ to $z=0.532$, and we estimate the mass of the central black hole, $\log (M_{\rm BH}/M_{\odot}) = 9.08 \pm 0.18$. The X-ray spectrum of this source can be well described by a power-law of photon index $Γ= 1.82 \pm 0.05$ with $E_{\rm cut} = 152_{-54}^{+103}\,{\rm keV}$, in the rest frame of the source. Assuming a Comptonisation model, we estimate the coronal temperature to be $kT_{\rm e}=42\pm 3 \,{\rm keV}$ and $kT_{\rm e}= 56 \pm 3\,{\rm keV}$ for a spherical and a slab geometry, respectively. The coronal properties are comparable to the ones derived for local AGN, despite a difference of around one order of magnitude in black hole mass and X-ray luminosity ($L_{2-10} = 1.93\times 10^{45}\,{\rm erg\,s^{-1}}$). The quasar is X-ray loud, with an unusually flat observed optical-to-X-ray spectral slope $α_{\rm OX} = 1.00 \pm 0.02$, and has an exceptionally strong optical [O III] line. Assuming that both the X-ray emission and the [O III] line are isotropic, these two extreme properties can be explained by a nearly edge-on disk, leading to a reduction in the observed UV continuum light.

preprint2016arXiv

Discovery of coherent pulsations from the Ultraluminous X-ray Source NGC 7793 P13

We report the detection of coherent pulsations from the ultraluminous X-ray source NGC 7793 P13. The ~0.42s nearly sinusoidal pulsations were initially discovered in broadband X-ray observations using XMM-Newton and NuSTAR taken in 2016. We subsequently also found pulsations in archival XMM-Newton data taken in 2013 and 2014. The significant (>>5 sigma) detection of coherent pulsations demonstrates that the compact object in P13 is a neutron star with an observed peak luminosity of ~1e40 erg/s (assuming isotropy), well above the Eddington limit for a 1.4 M_sun accretor. This makes P13 the second ultraluminous X-ray source known to be powered by an accreting neutron star. The pulse period varies between epochs, with a slow but persistent spin up over the 2013-2016 period. This spin-up indicates a magnetic field of B ~ 1.5e12 G, typical of many accreting pulsars. The most likely explanation for the extreme luminosity is a high degree of beaming, however this is difficult to reconcile with the sinusoidal pulse profile.

preprint2016arXiv

Disk-Wind Connection During the Heartbeats of GRS 1915+105

Disk and wind signatures are seen in the soft state of Galactic black holes, while the jet is seen in the hard state. Here we study the disk-wind connection in the $ρ$ class of variability in GRS 1915+105 using a joint NuSTAR-Chandra observation. The source shows 50 sec limit cycle oscillations. By including new information provided by the reflection spectrum, and using phase-resolved spectroscopy, we find that the change in the inner disk inferred from the blackbody emission is not matched by reflection measurements. The latter is almost constant, independent of the continuum model. The two radii are comparable only if the disk temperature color correction factor changes, an effect that could be due to the changing opacity of the disk caused by changes in metal abundances. The disk inclination is similar to that inferred from the jet axis, and oscillates by ~10 deg. The simultaneous Chandra data show the presence of two wind components with velocities between 500-5000 km/s, and possibly two more with velocities reaching 20,000 km/s (~0.06 c). The column densities are ~5e22 cm$^{-2}$. An upper limit to the wind response time of 2 sec is measured, implying a launch radius of <6e10 cm. The changes in wind velocity and absorbed flux require the geometry of the wind to change during the oscillations, constraining the wind to be launched from a distance of 290 - 1300 rg from the black hole. Both datasets support fundamental model predictions in which a bulge originates in the inner disk and moves outward as the instability progresses.

preprint2016arXiv

Extremes of the jet-accretion power relation of blazars, as explored by NuSTAR

Hard X-ray observations are crucial to study the non-thermal jet emission from high-redshift, powerful blazars. We observed two bright z>2 flat spectrum radio quasars (FSRQs) in hard X-rays to explore the details of their relativistic jets and their possible variability. S5 0014+81 (at z=3.366) and B0222+185 (at z=2.690) have been observed twice by the Nuclear Spectroscopic Telescope Array (NuSTAR) simultaneously with Swift/XRT, showing different variability behaviours. We found that NuSTAR is instrumental to explore the variability of powerful high-redshift blazars, even when no gamma-ray emission is detected. The two sources have proven to have respectively the most luminous accretion disk and the most powerful jet among known blazars. Thanks to these properties, they are located at the extreme end of the jet-accretion disk relation previously found for gamma-ray detected blazars, to which they are consistent.

preprint2016arXiv

First NuSTAR observations of the BL Lac - type blazar PKS~2155-304: constraints on the jet content and distribution of radiating particles

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.

preprint2016arXiv

Hard X-ray Emission from Sh2-104: A NuSTAR search for Gamma-ray Counterparts

We present NuSTAR hard X-ray observations of Sh 2-104, a compact HII region containing several young massive stellar clusters (YMSCs). We have detected distinct hard X-ray sources coincident with localized VERITAS TeV emission recently resolved from the giant gamma-ray complex MGRO J2019+37 in the Cygnus region. Faint, diffuse X-ray emission coincident with the eastern YMSC in Sh2-104 is likely the result of colliding winds of component stars. Just outside the radio shell of Sh 2-104 lies 3XMM J201744.7+365045 and a nearby nebula NuSTAR J201744.3+364812, whose properties are most consistent with extragalactic objects. The combined XMM-Newton and NuSTAR spectrum of 3XMM J201744.7+365045 is well-fit to an absorbed power-law model with NH = (3.1 +/- 1.0)E22 cm^-2 and photon index Gamma = 2.1 +/- 0.1. Based on possible long-term flux variation and the lack of detected pulsations (< 43% modulation), this object is likely a background AGN rather than a Galactic pulsar. The spectrum of the NuSTAR nebula shows evidence of an emission line at E = 5.6 keV suggesting an optically obscured galaxy cluster at z = 0.19 +/- 0.02 (d = 800 Mpc) and Lx = 1.2E44 erg/s. Follow-up Chandra observations of Sh 2-104 will help identify the nature of the X-ray sources and their relation to MGRO J2019+37. We also show that the putative VERITAS gamma-ray excess south of Sh 2-104 is most likely associated with the newly discovered Fermi pulsar PSR J2017+3625 and not the HII region.

preprint2016arXiv

Hot Dust Obscured Galaxies with Excess Blue Light: Dual AGN or Single AGN Under Extreme Conditions?

Hot Dust-Obscured Galaxies (Hot DOGs) are a population of hyper-luminous infrared galaxies identified by the WISE mission from their very red mid-IR colors, and characterized by hot dust temperatures ($T>60~\rm K$). Several studies have shown clear evidence that the IR emission in these objects is powered by a highly dust-obscured AGN that shows close to Compton-thick absorption at X-ray wavelengths. Thanks to the high AGN obscuration, the host galaxy is easily observable, and has UV/optical colors usually consistent with those of a normal galaxy. Here we discuss a sub-population of 8 Hot DOGs that show enhanced rest-frame UV/optical emission. We discuss three scenarios that might explain the excess UV emission: (i) unobscured light leaked from the AGN by reflection over the dust or by partial coverage of the accretion disk; (ii) a second unobscured AGN in the system; or (iii) a luminous young starburst. X-ray observations can help discriminate between these scenarios. We study in detail the blue excess Hot DOG WISE J020446.13-050640.8, which was serendipitously observed by Chandra/ACIS-I for 174.5 ks. The X-ray spectrum is consistent with a single, hyper-luminous, highly absorbed AGN, and is strongly inconsistent with the presence of a secondary unobscured AGN. Based on this, we argue that the excess blue emission in this object is most likely either due to reflection or a co-eval starburst. We favor the reflection scenario as the unobscured star-formation rate needed to power the UV/optical emission would be $\gtrsim 1000~\rm M_{\odot}~\rm yr^{-1}$. Deep polarimetry observations could confirm the reflection hypothesis.

preprint2016arXiv

IC 3639 - A new bona fide Compton thick AGN unveiled by NuSTAR

We analyse high-quality NuSTAR observations of the local (z = 0.011) Seyfert 2 active galactic nucleus (AGN) IC 3639, in conjunction with archival Suzaku and Chandra data. This provides the first broadband X-ray spectral analysis of the source, spanning nearly two decades in energy (0.5-30 keV). Previous X-ray observations of the source below 10 keV indicated strong reflection/obscuration on the basis of a pronounced iron fluorescence line at 6.4 keV. The hard X-ray energy coverage of NuSTAR, together with self-consistent toroidal reprocessing models, enables direct broadband constraints on the obscuring column density of the source. We find the source to be heavily Compton-thick (CTK) with an obscuring column in excess of $3.6\times10^{24}$ cm$^{-2}$, unconstrained at the upper end. We further find an intrinsic 2-10 keV luminosity of $\textrm{log}_{10}(L_{\textrm{2-10 keV}} \textrm{[erg s}^{-1}]) = 43.4^{+0.6}_{-1.1}$ to 90% confidence, almost 400 times the observed flux, and consistent with various multi-wavelength diagnostics. Such a high intrinsic to observed flux ratio in addition to an Fe-K$α$ fluorescence line equivalent width exceeding 2 keV is extreme amongst known bona fide CTK AGN, which we suggest are both due to the high level of obscuration present around IC 3639. Our study demonstrates that broadband spectroscopic modelling with NuSTAR enables large corrections for obscuration to be carried out robustly, and emphasises the need for improved modelling of AGN tori showing intense iron fluorescence.

preprint2016arXiv

IC 751: a new changing-look AGN discovered by NuSTAR

We present the results of five NuSTAR observations of the type 2 active galactic nucleus (AGN) in IC 751, three of which were performed simultaneously with XMM-Newton or Swift/XRT. We find that the nuclear X-ray source underwent a clear transition from a Compton-thick ($N_{\rm\,H}\simeq 2\times 10^{24}\rm\,cm^{-2}$) to a Compton-thin ($N_{\rm\,H}\simeq 4\times 10^{23}\rm\,cm^{-2}$) state on timescales of $\lesssim 3$ months, which makes IC 751 the first changing-look AGN discovered by NuSTAR. Changes of the line-of-sight column density at a $\sim2σ$ level are also found on a time-scale of $\sim 48$ hours ($ΔN_{\rm\,H}\sim 10^{23}\rm\,cm^{-2}$). From the lack of spectral variability on timescales of $\sim 100$ ks we infer that the varying absorber is located beyond the emission-weighted average radius of the broad-line region, and could therefore be related either to the external part of the broad-line region or a clumpy molecular torus. By adopting a physical torus X-ray spectral model, we are able to disentangle the column density of the non-varying absorber ($N_{\rm\,H}\sim 3.8\times 10^{23}\rm\,cm^{-2}$) from that of the varying clouds [$N_{\rm\,H}\sim(1-150)\times10^{22}\rm\,cm^{-2}$], and to constrain that of the material responsible for the reprocessed X-ray radiation ($N_{\rm\,H} \sim 6 \times 10^{24}\rm\,cm^{-2}$). We find evidence of significant intrinsic X-ray variability, with the flux varying by a factor of five on timescales of a few months in the 2-10 and 10-50 keV band.

preprint2016arXiv

IGR J18293-1213 is an eclipsing Cataclysmic Variable

Studying the population of faint hard X-ray sources along the plane of the Galaxy is challenging because of high-extinction and crowding, which make the identification of individual sources more difficult. IGR J18293-1213 is part of the population of persistent sources which have been discovered by the INTEGRAL satellite. We report on NuSTAR and Swift/XRT observations of this source, performed on 2015 September 11. We detected three eclipsing intervals in the NuSTAR light curve, allowing us to constrain the duration of these eclipses, $Δt = 30.8^{+6.3}_{-0.0}$ min, and the orbital period of the system, $T = 6.92\pm0.01$ hr. Even though we only report an upper limit on the amplitude of a putative spin modulation, the orbital period and the hard thermal Bremsstrahlung spectrum of IGR J18293-1213 provide strong evidence that this source is a magnetic Cataclysmic Variable (CV). Our NuSTAR and Swift/XRT joint spectral analysis places strong constraints on the white dwarf mass $M_{\rm wd} = 0.78^{+0.10}_{-0.09}$ M$_\odot$. Assuming that the mass to radius ratio of the companion star $M_* / R_* = 1$ (solar units) and using $T$, $Δt$ and $M_{\rm wd}$, we derived the mass of the companion star $M_* = 0.82\pm0.01$ M$_\odot$, the orbital separation of the binary system $a=2.14\pm0.04$ R$_\odot$, and its orbital inclination compared to the line of sight $i=(72.2^{+2.4}_{-0.0})\pm1.0^\circ$.

preprint2016arXiv

Keck/MOSFIRE spectroscopy of five ULX counterparts

We present H-band spectra of the candidate counterparts of five ULXs (two in NGC 925, two in NGC 4136, and Holmberg II X-1) obtained with Keck/MOSFIRE. The candidate counterparts of two ULXs (J022721+333500 in NGC 925 and J120922+295559 in NGC 4136) have spectra consistent with (M-type) red supergiants (RSGs). We obtained two epochs of spectroscopy of the candidate counterpart to J022721+333500, separated by 10 months, but discovered no radial velocity variations with a 2-$σ$ upper limit of 40 km/s. If the RSG is the donor star of the ULX, the most likely options are that either the system is seen at low inclination (< 40$^\circ$), or the black hole mass is less than 100 M$_\odot$, unless the orbital period is longer than 6 years, in which case the obtained limit is not constraining. The spectrum of the counterpart to J120922+295559 shows emission lines on top of its stellar spectrum, and the remaining three counterparts do not show absorption lines that can be associated with the atmosphere of a star; their spectra are instead dominated by emission lines. Those counterparts with RSG spectra may be used in the future to search for radial velocity variations, and, if those are present, determine dynamical constraints on the mass of the accretor.

preprint2016arXiv

NuSTAR and Swift observations of the very high state in GX 339-4: Weighing the black hole with X-rays

We present results from spectral fitting of the very high state of GX~339-4 with NuSTAR and Swift. We use relativistic reflection modelling to measure the spin of the black hole and inclination of the inner disk, and find a spin of $a=0.95^{+0.02}_{-0.08}$ and inclination of $30$°$\pm1$ (statistical errors). These values agree well with previous results from reflection modelling. With the exceptional sensitivity of NuSTAR at the high-energy side of the disk spectrum, we are able to constrain multiple physical parameters simultaneously using continuum fitting. By using the constraints from reflection as input for the continuum fitting method, we invert the conventional fitting procedure to estimate the mass and distance of GX~339-4 using just the X-ray spectrum, finding a mass of $9.0^{+1.6}_{-1.2}M_\odot$ and distance of $8.4\pm0.9$ kpc (statistical errors).

preprint2016arXiv

NuSTAR and XMM-Newton Observations of the Hard X-Ray Spectrum of Centaurus A

We present simultaneous XMM-Newton and NuSTAR observations spanning 3-78 keV of the nearest radio galaxy, Centaurus A (Cen A). The accretion geometry around the central engine in Cen A is still debated, and we investigate possible configurations using detailed X-ray spectral modeling. NuSTAR imaged the central region of Cen A with sub-arcminute resolution at X-ray energies above 10 keV for the first time, but finds no evidence for an extended source or other off-nuclear point-sources. The XMM-Newton and NuSTAR spectra agree well and can be described with an absorbed power-law with a photon index Γ = 1.815 +/- 0.005 and a fluorescent Fe Kα line in good agreement with literature values. The spectrum does not require a high-energy exponential rollover, with a constraint of E_fold > 1 MeV. A thermal Comptonization continuum describes the data well, with parameters that agree with values measured by INTEGRAL, in particular an electron temperature kTe between ~100-300 keV, seed photon input temperatures between 5-50 eV. We do not find evidence for reflection or a broad iron line and put stringent upper limits of R < 0.01 on the reflection fraction and accretion disk illumination. We use archival Chandra data to estimate the contribution from diffuse emission, extra-nuclear point-sources, and the outer X-ray jet to the observed NuSTAR and XMM-Newton X-ray spectra and find the contribution to be negligible. We discuss different scenarios for the physical origin of the observed hard X-ray spectrum, and conclude that the inner disk is replaced by an advection-dominated accretion flow or that the X-rays are dominated by synchrotron self-Compton emission from the inner regions of the radio jet or a combination thereof.

preprint2016arXiv

NuSTAR observations of water megamaser AGN

Aims. Study the connection between the masing disk and obscuring torus in Seyfert 2 galaxies. Methods. We present a uniform X-ray spectral analysis of the high energy properties of 14 nearby megamaser Active Galactic Nuclei observed by NuSTAR. We use a simple analytical model to localize the maser disk and understand its connection with the torus by combining NuSTAR spectral parameters with available physical quantities from VLBI mapping. Results. Most of the sources analyzed are heavily obscured, showing a column density in excess of $\sim 10^{23}$ cm$^{-2}$. In particular, $79\%$ are Compton-thick ($N_{\rm H} > 1.5 \times 10^{24}$ cm$^{-2}$). Using column densities measured by NuSTAR, with the assumption that the torus is the extension of the maser disk, and further assuming a reasonable density profile, the torus dimensions can be predicted. They are found to be consistent with mid-IR interferometry parsec-scale observations of Circinus and NGC 1068. In this picture, the maser disk is intimately connected to the inner part of the torus. It is probably made of a large number of molecular clouds connecting the torus and the outer part of the accretion disk, giving rise to a thin disk rotating in most cases in Keplerian or sub-Keplerian motion. This toy model explains the established close connection between water megamaser emission and nuclear obscuration as a geometric effect.

preprint2016arXiv

NuSTAR reveals the extreme properties of the super-Eddington accreting SMBH in PG 1247+267

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 ($Γ=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.

preprint2016arXiv

NuSTAR unveils a heavily obscured low-luminosity Active Galactic Nucleus in the Luminous Infrared Galaxy NGC 6286

We report the detection of a heavily obscured Active Galactic Nucleus (AGN) in the luminous infrared galaxy (LIRG) NGC 6286, identified in a 17.5 ks NuSTAR observation. The source is in an early merging stage, and was targeted as part of our ongoing NuSTAR campaign observing local luminous and ultra-luminous infrared galaxies in different merger stages. NGC 6286 is clearly detected above 10 keV and, by including the quasi-simultaneous Swift/XRT and archival XMM-Newton and Chandra data, we find that the source is heavily obscured [$N_{\rm\,H}\simeq (0.95-1.32)\times 10^{24}\rm\,cm^{-2}$], with a column density consistent with being Compton-thick [CT, $\log (N_{\rm\,H}/\rm cm^{-2})\geq 24$]. The AGN in NGC 6286 has a low absorption-corrected luminosity ($L_{2-10\rm\,keV}\sim 3-20\times 10^{41}\rm\,erg\,s^{-1}$) and contributes $\lesssim$1\% to the energetics of the system. Because of its low-luminosity, previous observations carried out in the soft X-ray band ($<10$ keV) and in the infrared did not notice the presence of a buried AGN. NGC 6286 has multi-wavelength characteristics typical of objects with the same infrared luminosity and in the same merger stage, which might imply that there is a significant population of obscured low-luminosity AGN in LIRGs that can only be detected by sensitive hard X-ray observations.

preprint2016arXiv

NuSTAR, Swift, and GROND observations of the flaring MeV blazar: PMN J0641$-$0320

MeV blazars are a sub--population of the blazar family, exhibiting larger--than--average jet powers, accretion luminosities and black hole masses. Because of their extremely hard X--ray continua, these objects are best studied in the X-ray domain. Here, we report on the discovery by the $Fermi$ Large Area Telescope and subsequent follow-up observations with $NuSTAR$, $Swift$ and GROND of a new member of the MeV blazar family: PMN J0641$-$0320. Our optical spectroscopy provides confirmation that this is a flat--spectrum radio quasar located at a redshift of $z=1.196$. Its very hard $NuSTAR$ spectrum (power--law photon index of $\sim$1 up to $\sim$80 keV) indicates that the emission is produced via inverse Compton scattering off photons coming from outside the jet.The overall spectral energy distribution of PMN J0641$-$0320 is typical of powerful blazars and by reproducing it with a simple one-zone leptonic emission model we find the emission region to be located either inside the broad line region or within the dusty torus.

preprint2016arXiv

On the black hole mass of the gamma-ray emitting narrow-line Seyfert 1 galaxy 1H 0323+342

Narrow-line Seyfert 1 galaxies have been identified by the Fermi Gamma-Ray Space Telescope as a rare class of gamma-ray emitting active galactic nuclei (AGN). The lowest-redshift candidate among them is the source 1H 0323+342. Here we present quasi-simultaneous Gemini near-infrared and Keck optical spectroscopy for it, from which we derive a black hole mass based on both the broad Balmer and Paschen emission lines. We supplement these observations with a NuSTAR X-ray spectrum taken about two years earlier, from which we constrain the black hole mass based on the short timescale spectral variability. Our multiwavelength observations suggest a black hole mass of ~2x10^7 solar masses, which agrees well with previous estimates. We build the spectral energy distribution and show that it is dominated by the thermal and reprocessed emission from the accretion disc rather than the non-thermal jet component. A detailed spectral fitting with the energy-conserving accretion disc model of Done et al. constrains the Eddington ratio to L/L_Edd ~ 0.5 for a (non-rotating) Schwarzschild black hole and to L/L_Edd ~ 1 for a Kerr black hole with dimensionless spin of a*=0.8. Higher spin values and so higher Eddington ratios are excluded, since they would strongly overpredict the observed soft X-ray flux.

preprint2016arXiv

Optical and near-infrared observations of SN 2013dx associated with GRB 130702A

We present optical and near-infrared light curves and optical spectra of SN 2013dx, associated with the nearby (redshift 0.145) gamma-ray burst GRB 130702A. The prompt isotropic gamma-ray energy released from GRB 130702A is measured to be $E_{γ,iso}=6.4_{-1.0}^{+1.3}\times10^{50}$erg (1keV-10MeV in the rest frame), placing it intermediate between low-luminosity GRBs like GRB 980425/SN 1998bw and the broader cosmological population. We compare the observed $g'r'i'z'$ light curves of SN 2013dx to a SN 1998bw template, finding that SN 2013dx evolves ~20% faster (steeper rise time), with a comparable peak luminosity. Spectroscopically, SN 2013dx resembles other broad-lined Type Ic supernovae, both associated with (SN 2006aj and SN 1998bw) and lacking (SN 1997ef, SN 2007I, and SN 2010ah) gamma-ray emission, with photospheric velocities around peak of ~21,000km s$^{-1}$. We construct a quasi-bolometric ($g'r'i'z'yJ$) light curve for SN 2013dx, only the fifth GRB-associated SN with extensive NIR coverage and the third with a bolometric light curve extending beyond $Δt>40$d. Together with the measured photospheric velocity, we derive basic explosion parameters using simple analytic models. We infer a $^{56}$Ni mass of $M_{\mathrm{Ni}}=0.37\pm0.01$M$_{\odot}$, an ejecta mass of $M_{\mathrm{ej}}=3.1\pm0.1$M$_{\odot}$, and a kinetic energy of $E_{\mathrm{K}}=(8.2\pm0.43)\times10^{51}$ erg (statistical uncertainties only), consistent with previous GRB-associated SNe. When considering the ensemble population of GRB-associated SNe, we find no correlation between the mass of synthesized $^{56}$Ni and high-energy properties, despite clear predictions from numerical simulations that $M_{\mathrm{Ni}}$ should correlate with the degree of asymmetry. On the other hand, $M_{\mathrm{Ni}}$ clearly correlates with the kinetic energy of the supernova ejecta across a wide range of core-collapse events.

preprint2016arXiv

Spectro-timing study of GX 339-4 in a hard intermediate state

We present an analysis of NuSTAR observations of a hard intermediate state of the transient black hole GX 339-4 taken in January 2015. As the source softened significantly over the course of the 1.3 d-long observation we split the data into 21 sub-sets and find that the spectrum of all of them can be well described by a power-law continuum with an additional relativistically blurred reflection component. The photon index increases from ~1.69 to ~1.77 over the course of the observation. The accretion disk is truncated at around 9 gravitational radii in all spectra. We also perform timing analysis on the same 21 individual data sets, and find a strong type-C quasi-periodic oscillation (QPO), which increase in frequency from ~0.68 to ~1.05 Hz with time. The frequency change is well correlated with the softening of the spectrum. We discuss possible scenarios for the production of the QPO and calculate predicted inner radii in the relativistic precession model as well as the global disk mode oscillations model. We find discrepancies with respect to the observed values in both models unless we allow for a black hole mass of ~100 M_sun , which is highly unlikely. We discuss possible systematic uncertainties, in particular with the measurement of the inner accretion disk radius in the relativistic reflection model. We conclude that the combination of observed QPO frequencies and inner accretion disk radii, as obtained from spectral fitting, is difficult to reconcile with current models.

preprint2016arXiv

The Chandra COSMOS Legacy survey: overview and point source catalog

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

preprint2016arXiv

The first X-ray imaging spectroscopy of quiescent solar active regions with NuSTAR

We present the first observations of quiescent active regions (ARs) using NuSTAR, a focusing hard X-ray telescope capable of studying faint solar emission from high temperature and non-thermal sources. We analyze the first directly imaged and spectrally resolved X-rays above 2~keV from non-flaring ARs, observed near the west limb on 2014 November 1. The NuSTAR X-ray images match bright features seen in extreme ultraviolet and soft X-rays. The NuSTAR imaging spectroscopy is consistent with isothermal emission of temperatures $3.1-4.4$~MK and emission measures $1-8\times 10^{46}$~cm$^{-3}$. We do not observe emission above 5~MK but our short effective exposure times restrict the spectral dynamic range. With few counts above 6~keV, we can place constraints on the presence of an additional hotter component between 5 and 12~MK of $\sim 10^{46}$cm$^{-3}$ and $\sim 10^{43}$ cm$^{-3}$, respectively, at least an order of magnitude stricter than previous limits. With longer duration observations and a weakening solar cycle (resulting in an increased livetime), future NuSTAR observations will have sensitivity to a wider range of temperatures as well as possible non-thermal emission.

preprint2016arXiv

The nature of the torus in the heavily obscured AGN Markarian 3: an X-ray study

In this paper we report the results of an X-ray monitoring campaign on the heavily obscured Seyfert galaxy Markarian 3 carried out between the fall of 2014 and the spring of 2015 with NuSTAR, Suzaku and XMM-Newton. The hard X-ray spectrum of Markarian 3 is variable on all the time scales probed by our campaign, down to a few days. The observed continuum variability is due to an intrinsically variable primary continuum seen in transmission through a large, but still Compton-thin column density (N_H~0.8-1.1$\times$10$^{24}$ cm$^{-2}$). If arranged in a spherical-toroidal geometry, the Compton scattering matter has an opening angle ~66 degrees and is seen at a grazing angle through its upper rim (inclination angle ~70 degrees). We report a possible occultation event during the 2014 campaign. If the torus is constituted by a system of clouds sharing the same column density, this event allows us to constrain their number (17$\pm$5) and individual column density, [~(4.9$\pm$1.5)$\times$10$^{22}$ cm$^{-2}$]. The comparison of IR and X-ray spectroscopic results with state-of-the art "torus" models suggests that at least two thirds of the X-ray obscuring gas volume might be located within the dust sublimation radius. We report also the discovery of an ionized absorber, characterised by variable resonant absorption lines due to He- and H-like iron. This discovery lends support to the idea that moderate column density absorbers could be due to clouds evaporated at the outer surface of the torus, possibly accelerated by the radiation pressure due to the central AGN emission leaking through the patchy absorber.

preprint2016arXiv

The NuSTAR Extragalactic Surveys: The Number Counts of Active Galactic Nuclei and the Resolved Fraction of the Cosmic X-ray Background

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.

preprint2016arXiv

The Pan-STARRS1 distant z>5.6 quasar survey: more than 100 quasars within the first Gyr of the universe

Luminous quasars at z>5.6 can be studied in detail with the current generation of telescopes and provide us with unique information on the first gigayear of the universe. Thus far these studies have been statistically limited by the number of quasars known at these redshifts. Such quasars are rare and therefore wide-field surveys are required to identify them and multiwavelength data are needed to separate them efficiently from their main contaminants, the far more numerous cool dwarfs. In this paper, we update and extend the selection for z~6 quasars presented in Banados et al. (2014) using the Pan-STARRS1 (PS1) survey. We present the PS1 distant quasar sample, which currently consists of 124 quasars in the redshift range 5.6<z<6.7 that satisfy our selection criteria. Seventy-seven of these quasars have been discovered with PS1, and 63 of them are newly identified in this paper. We present composite spectra of the PS1 distant quasar sample. This sample spans a factor of ~20 in luminosity and shows a variety of emission line properties. The number of quasars at z>5.6 presented in this work almost double the quasars previously known at these redshifts, marking a transition phase from studies of individual sources to statistical studies of the high-redshift quasar population, which was impossible with earlier, smaller samples.

preprint2016arXiv

The Soft State of Cygnus X-1 Observed with NuSTAR: A Variable Corona and a Stable Inner Disk

We present a multi-epoch hard X-ray analysis of Cygnus X-1 in its soft state based on four observations with NuSTAR. Despite the basic similarity of the observed spectra, there is clear spectral variability between epochs. To investigate this variability, we construct a model incorporating both the standard disk-corona continuum and relativistic reflection from the accretion disk, based on prior NuSTAR work on Cygnus X-1, and apply this model to each epoch independently. We find excellent consistency for the black hole spin, and the iron abundance of the accretion disk, which are expected to remain constant on observational timescales. In particular, we confirm that Cygnus X-1 hosts a rapidly rotating black hole, 0.93<a*<0.96, in broad agreement with the majority of prior studies of the relativistic disk reflection and constraints on the spin obtained through studies of the thermal accretion disk continuum. Our work also confirms the apparent misalignment between the inner disk and the orbital plane of the binary system reported previously, finding the magnitude of this warp to be ~10-15deg. This level of misalignment does not significantly change (and may even improve) the agreement between our reflection results and the thermal continuum results regarding the black hole spin. The spectral variability observed by NuSTAR is dominated by the primary continuum, implying variability in the temperature of the scattering electron plasma. Finally, we consistently observe absorption from ionized iron at ~6.7 keV, which varies in strength as a function of orbital phase in a manner consistent with the absorbing material being an ionized phase of the focused stellar wind from the supergiant companion star.

preprint2015arXiv

A Hard X-Ray Study of Ultraluminous X-ray Source NGC 5204 X-1 with NuSTAR and XMM-Newton

We present the results from coordinated X-ray observations of the ultraluminous X-ray source NGC 5204 X-1 performed by NuSTAR and XMM-Newton in early 2013. These observations provide the first detection of NGC 5204 X-1 above 10 keV, extending the broadband coverage to 0.3-20 keV. The observations were carried out in two epochs separated by approximately 10 days, and showed little spectral variation, with an observed luminosity of Lx = (4.95+/-0.11)e39 erg/s. The broadband spectrum confirms the presence of a clear spectral downturn above 10 keV seen in some previous observations. This cutoff is inconsistent with the standard low/hard state seen in Galactic black hole binaries, as would be expected from an intermediate mass black hole accreting at significantly sub-Eddington rates given the observed luminosity. The continuum is apparently dominated by two optically thick thermal-like components, potentially accompanied by a faint high energy tail. The broadband spectrum is likely associated with an accretion disk that differs from a standard Shakura & Sunyaev thin disk.

preprint2015arXiv

A mature galaxy cluster at z=1.58 around the radio galaxy 7C1753+6311

We report on the discovery of a z=1.58 mature cluster around the high-redshift radio galaxy 7C1753+6311, first identified in the Clusters Around Radio-Loud AGN survey. Two-thirds of the excess galaxies within the central 1Mpc lie on a red sequence with a colour that is consistent with an average formation redshift of zf~3. We show that 80+/-6% of the red sequence galaxies in the cluster core are quiescent, while the remaining 20% are red due to dusty star formation. We demonstrate that the cluster has an enhanced quiescent galaxy fraction that is three times that of the control field. We also show that this enhancement is mass dependent: 91+/-9% of the M* >10^{10.5}Msun cluster galaxies are quiescent, compared to only 36+/-2% of field galaxies, whereas the fraction of quiescent galaxies with lower masses is the same in the cluster and field environments. The presence of a dense core and a well-formed, quiescent red sequence suggest that this is a mature cluster. This means that distant radio galaxies do not solely reside in young, uncollapsed protoclusters, rather they can be found in clusters in a wide range of evolutionary states.

preprint2015arXiv

A NuSTAR observation of the center of the Coma Cluster

We present the results of a 55ks NuSTAR observation of the core of the Coma Cluster. The global spectrum can be explained by thermal gas emission, with a conservative 90% upper limit to non-thermal inverse Compton (IC) emission of $5.1 \times 10^{-12}$ erg cm$^{-2}$ s$^{-1}$ in a 12 arcmin $\times$ 12 arcmin field of view . The brightness of the thermal component in this central region does not allow more stringent upper limits on the IC component when compared with non-imaging instruments with much larger fields of view where claims of detections have been made. Future mosaic NuSTAR observations of Coma will further address this issue. The temperature map shows a relatively uniform temperature distribution with a gradient from the hot northwest side to the cooler southeast, in agreement with previous measurements. The temperature determination is robust given the flat effective area and low background in the 3-20 keV band, making NuSTAR an ideal instrument to measure high temperatures in the intracluster medium.

preprint2015arXiv

A Soft X-Ray Spectral Episode for the Clocked Burster, GS 1826-24 as Measured by Swift and NuSTAR

We report on NuSTAR and Swift observations of a soft state of the neutron star low-mass X-ray binary GS 1826-24, commonly known as the "clocked" burster. The transition to the soft state was recorded in 2014 June through an increase of the 2-20 keV source intensity measured by MAXI, simultaneous with a decrease of the 15-50 keV intensity measured by Swift/BAT. The episode lasted approximately two months, after which the source returned to its usual hard state. We analyze the broad-band spectrum measured by Swift/XRT and NuSTAR, and estimate the accretion rate during the soft episode to be about 13% of Eddington, within the range of previous observations. However, the best fit spectral model, adopting the double Comptonization used previously, exhibits significantly softer components. We detect seven type-I X-ray bursts, all significantly weaker (and with shorter rise and decay times) than observed previously. The burst profiles and recurrence times vary significantly, ruling out the regular bursts that are typical for this source. One burst exhibited photospheric radius expansion, and we estimate the source distance at about (5.7 / xi_b^1/2) kpc, where xi_b parameterizes the possible anisotropy of the burst emission. Interpreting the soft state as a transition from an optically thin inner flow to an optically thick flow passing through a boundary layer, as is commonly observed in similar systems, is contradicted by the lower optical depth measured for the double Comptonization model we find for this soft state. The effect of a change in disk geometry on the burst behavior remains unclear.

preprint2015arXiv

An Extragalactic Spectroscopic Survey of the SSA22 Field

We present VLT VIMOS, Keck DEIMOS and Keck LRIS multi-object spectra of 367 sources in the field of the z ~ 3.09 protocluster SSA22. Sources are spectroscopically classified via template matching, allowing new identifications for 206 extragalactic sources, including 36 z > 2 Lyman-break galaxies (LBGs) and Lyman-α emitters (LAEs), 8 protocluster members, and 94 X-ray sources from the ~ 400 ks Chandra deep survey of SSA22. Additionally, in the area covered by our study, we have increased by ~ 4, 13, and 6 times the number of reliable redshifts of sources at 1.0 < z < 2.0, at z > 3.4, and with X-Ray emission, respectively. We compare our results with past spectroscopic surveys of SSA22 to investigate the completeness of the LBGs and the X-Ray properties of the new spectroscopically-classified sources in the SSA22 field.

preprint2015arXiv

Black hole feedback in the luminous quasar PDS 456

The evolution of galaxies is connected to the growth of supermassive black holes in their centers. During the quasar phase, a huge luminosity is released as matter falls onto the black hole, and radiation-driven winds can transfer most of this energy back to the host galaxy. Over five different epochs, we detected the signatures of a nearly spherical stream of highly ionized gas in the broadband X-ray spectra of the luminous quasar PDS 456. This persistent wind is expelled at relativistic speeds from the inner accretion disk, and its wide aperture suggests an effective coupling with the ambient gas. The outflow's kinetic power larger than 10^46 ergs per second is enough to provide the feedback required by models of black hole and host galaxy co-evolution.

preprint2015arXiv

Broad Iron Emission from Gravitationally Lensed Quasars Observed by Chandra

Recent work has demonstrated the potential of gravitationally lensed quasars to extend measurements of black hole spin out to high-redshift with the current generation of X-ray observatories. Here we present an analysis of a large sample of 27 lensed quasars in the redshift range 1.0<z<4.5 observed with Chandra, utilizing over 1.6 Ms of total observing time, focusing on the rest-frame iron K emission from these sources. Although the X-ray signal-to-noise (S/N) currently available does not permit the detection of iron emission from the inner accretion disk in individual cases in our sample, we find significant structure in the stacked residuals. In addition to the narrow core, seen almost ubiquitously in local AGN, we find evidence for an additional underlying broad component from the inner accretion disk, with a clear red wing to the emission profile. Based on simulations, we find the detection of this broader component to be significant at greater than the 3-sigma level. This implies that iron emission from the inner disk is relatively common in the population of lensed quasars, and in turn further demonstrates that, with additional observations, this population represents an opportunity to significantly extend the sample of AGN spin measurements out to high-redshift.

preprint2015arXiv

Broadband Observations of the Compton-thick Nucleus of NGC 3393

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, θ=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).

preprint2015arXiv

Constraining the radio-loud fraction of quasars at z>5.5

Radio-loud Active Galactic Nuclei at z~2-4 are typically located in dense environments and their host galaxies are among the most massive systems at those redshifts, providing key insights for galaxy evolution. Finding radio-loud quasars at the highest accessible redshifts (z~6) is important to study their properties and environments at even earlier cosmic time. They would also serve as background sources for radio surveys intended to study the intergalactic medium beyond the epoch of reionization in HI 21 cm absorption. Currently, only five radio-loud ($R=f_{ν,5{\rm GHz}}/f_{ν,4400Å}>10$) quasars are known at z~6. In this paper we search for 5.5 < z < 7.2 quasars by cross-matching the optical Pan-STARRS1 and radio FIRST surveys. The radio information allows identification of quasars missed by typical color-based selections. While we find no good 6.4 < z <7.2 quasar candidates at the sensitivities of these surveys, we discover two new radio-loud quasars at z~6. Furthermore, we identify two additional z~6 radio-loud quasars which were not previously known to be radio-loud, nearly doubling the current z~6 sample. We show the importance of having infrared photometry for z>5.5 quasars to robustly classify them as radio-quiet or radio-loud. Based on this, we reclassify the quasar J0203+0012 (z=5.72), previously considered radio-loud, to be radio-quiet. Using the available data in the literature, we constrain the radio-loud fraction of quasars at z~6, using the Kaplan--Meier estimator, to be $8.1^{+5.0}_{-3.2}\%$. This result is consistent with there being no evolution of the radio-loud fraction with redshift, in contrast to what has been suggested by some studies at lower redshifts.

preprint2015arXiv

Dark matter line emission constraints from NuSTAR observations of the Bullet Cluster

Line emission from dark matter is well motivated for some candidates e.g. sterile neutrinos. We present the first search for dark matter line emission in the 3-80keV range in a pointed observation of the Bullet Cluster with NuSTAR. We do not detect any significant line emission and instead we derive upper limits (95% CL) on the flux, and interpret these constraints in the context of sterile neutrinos and more generic dark matter candidates. NuSTAR does not have the sensitivity to constrain the recently claimed line detection at 3.5keV, but improves on the constraints for energies of 10-25keV.

preprint2015arXiv

Determining the covering factor of Compton-thick active galactic nuclei with NuSTAR

The covering factor of Compton-thick obscuring material associated with the torus in active galactic nuclei (AGN) is at present best understood through the fraction of sources exhibiting Compton-thick absorption along the line of sight ($N_{H}>1.5\times10^{24}$ cm$^{-2}$) in the X-ray band, which reveals the average covering factor. Determining this Compton-thick fraction is difficult however, due to the extreme obscuration. With its spectral coverage at hard X-rays ($>$10 keV), NuSTAR is sensitive to the AGN covering factor since Compton scattering of X-rays off optically thick material dominates at these energies. We present a spectral analysis of 10 AGN observed with NuSTAR where the obscuring medium is optically thick to Compton scattering, so called Compton-thick (CT) AGN. We use the torus models of Brightman & Nandra which predict the X-ray spectrum from reprocessing in a torus and include the torus opening angle as a free parameter and aim to determine the covering factor of the Compton-thick gas in these sources individually. Across the sample we find mild to heavy Compton-thick columns, with $N_{H}$ measured from $10^{24}-10^{26}$ cm$^{-2}$, and a wide range of covering factors, where individual measurements range from 0.2-0.9. We find that the covering factor, $f_{c}$, is a strongly decreasing function of the intrinsic 2-10 keV luminosity, $L_{X}$, where $f_{c}=(-0.41\pm0.13)$log$_{10}$($L_{X}$/erg s$^{-1}$)$+18.31\pm5.33$, across more than two orders of magnitude in $L_{X}$ (10$^{41.5}-10^{44}$ erg s$^{-1}$). The covering factors measured here agree well with the obscured fraction as a function of $L_{X}$ as determined by studies of local AGN with $L_{X}>10^{42.5}$ erg s$^{-1}$.

preprint2015arXiv

Distorted cyclotron line profile in Cep X-4 as observed by NuSTAR

We present spectral analysis of NuSTAR and Swift observations of Cep X-4 during its outburst in 2014. We observed the source once during the peak of the outburst and once during the decay, finding good agreement in the spectral shape between the observations. We describe the continuum using a powerlaw with a Fermi-Dirac cutoff at high energies. Cep X-4 has a very strong cyclotron resonant scattering feature (CRSF) around 30 keV. A simple absorption-like line with a Gaussian optical depth or a pseudo-Lorentzian profile both fail to describe the shape of the CRSF accurately, leaving significant deviations at the red side of the line. We characterize this asymmetry with a second absorption feature around 19 keV. The line energy of the CRSF, which is not influenced by the addition of this feature, shows a small but significant positive luminosity dependence. With luminosities between (1-6)e36 erg/s, Cep X-4 is below the theoretical limit where such a correlation is expected. This behavior is similar to Vela X-1 and we discuss parallels between the two systems.

preprint2015arXiv

First NuSTAR Observations of Mrk 501 within a Radio to TeV Multi-Instrument Campaign

We report on simultaneous broadband observations of the TeV-emitting blazar Markarian 501 between 1 April and 10 August 2013, including the first detailed characterization of the synchrotron peak with Swift and NuSTAR. During the campaign, the nearby BL Lac object was observed in both a quiescent and an elevated state. The broadband campaign includes observations with NuSTAR, MAGIC, VERITAS, the Fermi Large Area Telescope (LAT), Swift X-ray Telescope and UV Optical Telescope, various ground-based optical instruments, including the GASP-WEBT program, as well as radio observations by OVRO, Metsähovi and the F-Gamma consortium. Some of the MAGIC observations were affected by a sand layer from the Saharan desert, and had to be corrected using event-by-event corrections derived with a LIDAR (LIght Detection And Ranging) facility. This is the first time that LIDAR information is used to produce a physics result with Cherenkov Telescope data taken during adverse atmospheric conditions, and hence sets a precedent for the current and future ground-based gamma-ray instruments. The NuSTAR instrument provides unprecedented sensitivity in hard X-rays, showing the source to display a spectral energy distribution between 3 and 79 keV consistent with a log-parabolic spectrum and hard X-ray variability on hour timescales. None (of the four extended NuSTAR observations) shows evidence of the onset of inverse-Compton emission at hard X-ray energies. We apply a single-zone equilibrium synchrotron self-Compton model to five simultaneous broadband spectral energy distributions. We find that the synchrotron self-Compton model can reproduce the observed broadband states through a decrease in the magnetic field strength coinciding with an increase in the luminosity and hardness of the relativistic leptons responsible for the high-energy emission.

preprint2015arXiv

NuSTAR and multifrequency study of the two high-redshift blazars S5 0836+710 and PKS 2149-306

The most powerful blazars are the flat spectrum radio quasars whose emission is dominated by a Compton component peaking between a few hundred keV and a few hundred MeV. We selected two bright blazars, PKS 2149-306 at redshift z=2.345 and S5 0836+710 at z=2.172, in order to observe them in the hard X-ray band with the NuSTAR satellite. In this band the Compton component is rapidly rising almost up to the peak of the emission. Simultaneous soft-X-rays and UV-optical observations were performed with the Swift satellite, while near-infrared (NIR) data were obtained with the REM telescope. To study their variability, we repeated these observations for both sources on a timescale of a few months. While no fast variability was detected during a single observation, both sources were found to be variable in the X-ray band, up to 50%, between the two observations, with larger variability at higher energies. No variability was detected in the optical/NIR band. These data together with Fermi-LAT, WISE and other literature data are then used to study the overall spectral energy distributions (SEDs) of these blazars. Although the jet non-thermal emission dominates the SED, it leaves the UV band unhidden, allowing us to detect the thermal emission of the disc and to estimate the mass of the black hole. The non-thermal emission is well reproduced by a one-zone leptonic model. The non-thermal radiative processes are synchrotron, self-Compton and external Compton using seed photons from both the broad-line region (BLR) and the torus. We find that our data are better reproduced if we assume that the location of the dissipation region of the jet, R_diss, is in-between the torus, (at R_torus), and the BLR (R_torus>R_diss>R_BLR). The observed variability is explained by changing a minimum number of model parameters by a very small amount.

preprint2015arXiv

NuSTAR and Suzaku observations of the hard state in Cygnus X-1: locating the inner accretion disk

We present simultaneous Nuclear Spectroscopic Telescope Array (NuSTAR ) and Suzaku observations of the X-ray binary Cygnus X-1 in the hard state. This is the first time this state has been observed in Cyg X-1 with NuSTAR, which enables us to study the reflection and broad-band spectra in unprecedented detail. We confirm that the iron line cannot be fit with a combination of narrow lines and absorption features, and instead requires a relativistically blurred profile in combination with a narrow line and absorption from the companion wind. We use the reflection models of Garcia et al. (2014) to simultaneously measure the black hole spin, disk inner radius, and coronal height in a self-consistent manner. Detailed fits to the iron line profile indicate a high level of relativistic blurring, indicative of reflection from the inner accretion disk. We find a high spin, a small inner disk radius, and a low source height, and rule out truncation to greater than three gravitational radii at the 3σ confidence level. In addition, we find that the line profile has not changed greatly in the switch from soft to hard states, and that the differences are consistent with changes in the underlying reflection spectrum rather than the relativistic blurring. We find that the blurring parameters are consistent when fitting either just the iron line or the entire broad-band spectrum, which is well modelled with a Comptonized continuum plus reflection model.

preprint2015arXiv

NuSTAR and Suzaku X-ray Spectroscopy of NGC 4151: Evidence for Reflection from the Inner Accretion Disk

We present X-ray timing and spectral analyses of simultaneous 150 ks Nuclear Spectroscopic Telescope Array (NuSTAR) and Suzaku X-ray observations of the Seyfert 1.5 galaxy NGC 4151. We disentangle the continuum emission, absorption, and reflection properties of the active galactic nucleus (AGN) by applying inner accretion disk reflection and absorption-dominated models. With a time-averaged spectral analysis, we find strong evidence for relativistic reflection from the inner accretion disk. We find that relativistic emission arises from a highly ionized inner accretion disk with a steep emissivity profile, which suggests an intense, compact illuminating source. We find a preliminary, near-maximal black hole spin a>0.9 accounting for statistical and systematic modeling errors. We find a relatively moderate reflection fraction with respect to predictions for the lamp post geometry, in which the illuminating corona is modeled as a point source. Through a time-resolved spectral analysis, we find that modest coronal and inner disk reflection flux variation drives the spectral variability during the observations. We discuss various physical scenarios for the inner disk reflection model, and we find that a compact corona is consistent with the observed features.

preprint2015arXiv

NuSTAR Observations of the Compton-thick Active Galactic Nucleus and Ultraluminous X-ray Source Candidate in NGC 5643

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

preprint2015arXiv

NuSTAR Reveals Extreme Absorption in z < 0.5 Type 2 Quasars

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.

preprint2015arXiv

NuSTAR Reveals Relativistic Reflection But No Ultra-Fast Outflow In The Quasar PG 1211+143

We report on four epochs of observations of the quasar PG 1211+143 using NuSTAR. The net exposure time is 300 ks. Prior work on this source found suggestive evidence of an 'ultra-fast outflow' (or, UFO) in the Fe K band, with a velocity of approximately 0.1c. The putative flow would carry away a high mass flux and kinetic power, with broad implications for feedback and black hole-galaxy co-evolution. NuSTAR detects PG 1211+143 out to 30 keV, meaning that the continuum is well-defined both through and above the Fe K band. A characteristic relativistic disk reflection spectrum is clearly revealed, via a broad Fe K emission line and Compton back-scattering curvature. The data offer only weak constraints on the spin of the black hole. A careful search for UFO's show no significant absorption feature above 90% confidence. The limits are particularly tight when relativistic reflection is included. We discuss the statistics and the implications of these results in terms of connections between accretion onto quasars, Seyferts, and stellar-mass black holes, and feedback into their host environments.

preprint2015arXiv

NuSTAR, XMM-Newton and Suzaku Observations of the Ultraluminous X-ray Source Holmberg II X-1

We present the first broadband 0.3-25.0 kev X-ray observations of the bright ultraluminous X-ray source (ULX) Holmberg II X-1, performed by NuSTAR, XMM-Newton and Suzaku in September 2013. The NuSTAR data provide the first observations of Holmberg II X-1 above 10 keV, and reveal a very steep high-energy spectrum, similar to other ULXs observed by NuSTAR to date. These observations further demonstrate that ULXs exhibit spectral states that are not typically seen in Galactic black hole binaries. Comparison with other sources implies that Holmberg II X-1 accretes at a high fraction of its Eddington accretion rate, and possibly exceeds it. The soft X-ray spectrum (E<10 keV) appears to be dominated by two blackbody-like emission components, the hotter of which may be associated with an accretion disk. However, all simple disk models under-predict the NuSTAR data above ~10 keV and require an additional emission component at the highest energies probed, implying the NuSTAR data does not fall away with a Wien spectrum. We investigate physical origins for such an additional high-energy emission component, and favor a scenario in which the excess arises from Compton scattering in a hot corona of electrons with some properties similar to the very-high state seen in Galactic binaries. The observed broadband 0.3-25.0 keV luminosity inferred from these epochs is Lx = (8.1+/-0.1)e39 erg/s, typical for Holmberg II X-1, with the majority of the flux (~90%) emitted below 10 keV.

preprint2015arXiv

Optical spectroscopic observations of gamma-ray blazar candidates IV. Results of the 2014 follow-up campaign

The extragalactic gamma-ray sky is dominated by the emission arising from blazars, one of the most peculiar classes of radio-loud active galaxies. Since the launch of Fermi several methods were developed to search for blazars as potential counterparts of unidentified gamma-ray sources (UGSs). To confirm the nature of the selected candidates, optical spectroscopic observations are necessary. In 2013 we started a spectroscopic campaign to investigate gamma-ray blazar candidates selected according to different procedures. The main goals of our campaign are: 1) to confirm the nature of these candidates, and 2) whenever possible determine their redshifts. Optical spectroscopic observations will also permit us to verify the robustness of the proposed associations and check for the presence of possible source class contaminants to our counterpart selection. This paper reports the results of observations carried out in 2014 in the Northern hemisphere with Kitt Peak National Observatory (KPNO) and in the Southern hemisphere with the Southern Astrophysical Research (SOAR) telescopes. We also report three sources observed with the Magellan and Palomar telescopes. Our selection of blazar-like sources that could be potential counterparts of UGSs is based on their peculiar IR colors and on their combination with radio observations both at high and low frequencies (i.e., above and below ~ 1 GHz) in publicly available large radio surveys. We present the optical spectra of 27 objects. [....] We also report the case for WISE J173052.85-035247.2, candidate counterpart of the source 2FGL J1730.6-0353, which has no radio counterpart in the major radio surveys. We confirm that our selection of gamma-ray blazars candidates can successfully indentify low-energy counterparts to Fermi unassociated sources and allow us to discover new blazars.

preprint2015arXiv

Planck 2015 results. I. Overview of products and scientific results

The European Space Agency's Planck satellite, dedicated to studying the early Universe and its subsequent evolution, was launched 14~May 2009 and scanned the microwave and submillimetre sky continuously between 12~August 2009 and 23~October 2013. In February~2015, ESA and the Planck Collaboration released the second set of cosmology products based on data from the entire Planck mission, including both temperature and polarization, along with a set of scientific and technical papers and a web-based explanatory supplement. This paper gives an overview of the main characteristics of the data and the data products in the release, as well as the associated cosmological and astrophysical science results and papers. The science products include maps of the cosmic microwave background (CMB), the thermal Sunyaev-Zeldovich effect, and diffuse foregrounds in temperature and polarization, catalogues of compact Galactic and extragalactic sources (including separate catalogues of Sunyaev-Zeldovich clusters and Galactic cold clumps), and extensive simulations of signals and noise used in assessing the performance of the analysis methods and assessment of uncertainties. The likelihood code used to assess cosmological models against the Planck data are described, as well as a CMB lensing likelihood. Scientific results include cosmological parameters deriving from CMB power spectra, gravitational lensing, and cluster counts, as well as constraints on inflation, non-Gaussianity, primordial magnetic fields, dark energy, and modified gravity.

preprint2015arXiv

Planck 2015 results. XXVII. The Second Planck Catalogue of Sunyaev-Zeldovich Sources

We present the all-sky Planck catalogue of Sunyaev-Zeldovich (SZ) sources detected from the 29 month full-mission data. The catalogue (PSZ2) is the largest SZ-selected sample of galaxy clusters yet produced and the deepest all-sky catalogue of galaxy clusters. It contains 1653 detections, of which 1203 are confirmed clusters with identified counterparts in external data-sets, and is the first SZ-selected cluster survey containing > $10^3$ confirmed clusters. We present a detailed analysis of the survey selection function in terms of its completeness and statistical reliability, placing a lower limit of 83% on the purity. Using simulations, we find that the Y5R500 estimates are robust to pressure-profile variation and beam systematics, but accurate conversion to Y500 requires. the use of prior information on the cluster extent. We describe the multi-wavelength search for counterparts in ancillary data, which makes use of radio, microwave, infra-red, optical and X-ray data-sets, and which places emphasis on the robustness of the counterpart match. We discuss the physical properties of the new sample and identify a population of low-redshift X-ray under- luminous clusters revealed by SZ selection. These objects appear in optical and SZ surveys with consistent properties for their mass, but are almost absent from ROSAT X-ray selected samples.

preprint2015arXiv

Rapid Variability of Blazar 3C 279 during Flaring States in 2013-2014 with Joint Fermi-LAT, NuSTAR, Swift, and Ground-Based Multi-wavelength Observations

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 $γ$-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 $γ$-ray spectrum during one of the flares was very hard, with an index of $Γ_γ= 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_γ/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 $ΔΓ_{\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 $γ$-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 $γ$-ray and optical spectra cannot satisfactorily explain the production of X-ray emission.

preprint2015arXiv

Simultaneous NuSTAR/Chandra observations of the Bursting Pulsar GRO J1744-28 during its third reactivation

We report on a 10 ks simultaneous Chandra/HETG-NuSTAR observation of the Bursting Pulsar, GRO J1744-28, during its third detected outburst since discovery and after nearly 18 years of quiescence. The source is detected up to 60 keV with an Eddington persistent flux level. Seven bursts, followed by dips, are seen with Chandra, three of which are also detected with NuSTAR. Timing analysis reveals a slight increase in the persistent emission pulsed fraction with energy (from 10% to 15%) up to 10 keV, above which it remains constant. The 0.5-70 keV spectra of the persistent and dip emission are the same within errors, and well described by a blackbody (BB), a power-law with an exponential rolloff, a 10 keV feature, and a 6.7 keV emission feature, all modified by neutral absorption. Assuming that the BB emission originates in an accretion disc, we estimate its inner (magnetospheric) radius to be about 4x10^7 cm, which translates to a surface dipole field B~9x10^10 G. The Chandra/HETG spectrum resolves the 6.7 keV feature into (quasi-)neutral and highly ionized Fe XXV and Fe XXVI emission lines. XSTAR modeling shows these lines to also emanate from a truncated accretion disk. The burst spectra, with a peak flux more than an order of magnitude higher than Eddington, are well fit with a power-law with an exponential rolloff and a 10~keV feature, with similar fit values compared to the persistent and dip spectra. The burst spectra lack a thermal component and any Fe features. Anisotropic (beamed) burst emission would explain both the lack of the BB and any Fe components.

preprint2015arXiv

The 0.3-30 keV Spectra of Powerful Starburst Galaxies: NuSTAR and Chandra Observations of NGC 3256 and NGC 3310

We present nearly simultaneous Chandra and NuSTAR observations of two actively star-forming galaxies within 50 Mpc: NGC 3256 and NGC 3310. Both galaxies are detected by both Chandra and NuSTAR, which together provide the first-ever spectra of these two galaxies spanning 0.3-30 keV. The X-ray emission from both galaxies is spatially resolved by Chandra; we find that hot gas dominates the E < 1-3 keV emission while ultraluminous X-ray sources (ULXs) dominate at E > 1-3 keV. The NuSTAR galaxy-wide spectra of both galaxies follow steep power-law distributions with Gamma ~ 2.6 at E > 5-7 keV, similar to the spectra of bright individual ULXs and other galaxies that have been studied by NuSTAR. We find that both NGC 3256 and NGC 3310 have X-ray detected sources coincident with nuclear regions; however, the steep NuSTAR spectra of both galaxies restricts these sources to be either low luminosity AGN or non-AGN in nature (e.g., ULXs or crowded X-ray sources that reach L2-10 keV ~ 10^40 erg/s cannot be ruled out). Combining our constraints on the 0.3-30 keV spectra of NGC 3256 and NGC 3310 with equivalent measurements for nearby star-forming galaxies M83 and NGC 253, we analyze the SFR-normalized spectra of these starburst galaxies. The spectra of all four galaxies show sharply declining power-law slopes above 3-6 keV due to ULX populations. Our observations therefore constrain the average spectra of luminous accreting binaries (i.e., ULXs). This result is similar to the super-Eddington accreting ULXs that have been studied individually in a targeted NuSTAR ULX program. We also find that NGC 3310 exhibits a factor of ~3-10 elevation of X-ray emission over the other star-forming galaxies. We argue that the excess is most likely explained by the relatively low metallicity of the young stellar population in NGC 3310.

preprint2015arXiv

The complex accretion geometry of GX 339-4 as seen by NuSTAR and Swift

We present spectral analysis of five NuSTAR and Swift observations of GX 339-4 taken during a failed outburst in summer 2013. These observations cover Eddington luminosity fractions in the range ~0.9-6%. Throughout this outburst, GX 339-4 stayed in the hard state, and all five observations show similar X-ray spectra with a hard power-law with a photon index near 1.6 and significant contribution from reflection. Using simple reflection models we find unrealistically high iron abundances. Allowing for different photon indices for the continuum incident on the reflector relative to the underlying observed continuum results in a statistically better fit and reduced iron abundances. With a photon index around 1.3, the input power-law on the reflector is significantly harder than that which is directly observed. We study the influence of different emissivity profiles and geometries and consistently find an improvement when using separate photon indices. The inferred inner accretion disk radius is strongly model dependent, but we do not find evidence for a truncation radius larger than 100 r_g in any model. The data do not allow independent spin constraints but the results are consistent with the literature (i.e., a>0). Our best-fit models indicate an inclination angle in the range 40-60 degrees, consistent with limits on the orbital inclination but higher than reported in the literature using standard reflection models. The iron line around 6.4 keV is clearly broadened, and we detect a superimposed narrow core as well. This core originates from a fluorescence region outside the influence of the strong gravity of the black hole and we discuss possible geometries.

preprint2015arXiv

The Compton hump and variable blue wing in the extreme low-flux NuSTAR observations of 1H0707-495

The Narrow-line Seyfert I galaxy, 1H0707-495, has been well observed in the 0.3-10 keV band, revealing a dramatic drop in flux in the iron K alpha band, a strong soft excess, and short timescale reverberation lags associated with these spectral features. In this paper, we present the first results of a deep 250 ks NuSTAR observation of 1H0707-495, which includes the first sensitive observations above 10 keV. Even though the NuSTAR observations caught the source in an extreme low-flux state, the Compton hump is still significantly detected. NuSTAR, with its high effective area above 7 keV, clearly detects the drop in flux in the iron K alpha band, and by comparing these observations with archival XMM-Newton observations, we find that the energy of this drop increases with increasing flux. We discuss possible explanations for this, the most likely of which is that the drop in flux is the blue wing of the relativistically broadened iron K alpha emission line. When the flux is low, the coronal source height is low, thus enhancing the most gravitationally redshifted emission.

preprint2015arXiv

The formation history of massive cluster galaxies as revealed by CARLA

We use a sample of 37 of the densest clusters and protoclusters across $1.3 \le z \le 3.2$ from the Clusters Around Radio-Loud AGN (CARLA) survey to study the formation of massive cluster galaxies. We use optical $i'$-band and infrared 3.6$μ$m and 4.5$μ$m images to statistically select sources within these protoclusters and measure their median observed colours; $\langle i'-[3.6] \rangle$. We find the abundance of massive galaxies within the protoclusters increases with decreasing redshift, suggesting these objects may form an evolutionary sequence, with the lower redshift clusters in the sample having similar properties to the descendants of the high redshift protoclusters. We find that the protocluster galaxies have an approximately unevolving observed-frame $i'-[3.6]$ colour across the examined redshift range. We compare the evolution of the $\langle i'-[3.6] \rangle$ colour of massive cluster galaxies with simplistic galaxy formation models. Taking the full cluster population into account, we show that the formation of stars within the majority of massive cluster galaxies occurs over at least 2Gyr, and peaks at $z \sim 2$-3. From the median $i'-[3.6]$ colours we cannot determine the star formation histories of individual galaxies, but their star formation must have been rapidly terminated to produce the observed red colours. Finally, we show that massive galaxies at $z>2$ must have assembled within 0.5Gyr of them forming a significant fraction of their stars. This means that few massive galaxies in $z>2$ protoclusters could have formed via dry mergers.

preprint2015arXiv

The hard X-ray emission of the luminous infrared galaxy NGC 6240 as observed by NuSTAR

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.

preprint2015arXiv

The identification of z-dropouts in Pan-STARRS1: three quasars at 6.5<z<6.7

Luminous distant quasars are unique probes of the high redshift intergalactic medium (IGM) and of the growth of massive galaxies and black holes in the early universe. Absorption due to neutral Hydrogen in the IGM makes quasars beyond a redshift of z~6.5 very faint in the optical $z$-band, thus locating quasars at higher redshifts require large surveys that are sensitive above 1 micron. We report the discovery of three new z>6.5 quasars, corresponding to an age of the universe of <850 Myr, selected as z-band dropouts in the Pan-STARRS1 survey. This increases the number of known z>6.5 quasars from 4 to 7. The quasars have redshifts of z=6.50, 6.52, and 6.66, and include the brightest z-dropout quasar reported to date, PSO J036.5078+03.0498 with M_1450=-27.4. We obtained near-infrared spectroscopy for the quasars and from the MgII line we estimate that the central black holes have masses between 5x10^8 and 4x10^9 M_sun, and are accreting close to the Eddington limit (L_Bol/L_Edd=0.13-1.2). We investigate the ionized regions around the quasars and find near zone radii of R_NZ=1.5-5.2 proper Mpc, confirming the trend of decreasing near zone sizes with increasing redshift found for quasars at 5.7<z<6.4. By combining R_NZ of the PS1 quasars with those of 5.7<z<7.1 quasars in the literature, we derive a luminosity corrected redshift evolution of R_NZ,corrected=(7.2+/-0.2)-(6.1+/-0.7)x(z-6) Mpc. However, the large spread in R_NZ in the new quasars implies a wide range in quasar ages and/or a large variation in the neutral Hydrogen fraction along different lines of sight.

preprint2015arXiv

The Lick AGN Monitoring Project 2011: Spectroscopic Campaign and Emission-Line Light Curves

In the Spring of 2011 we carried out a 2.5 month reverberation mapping campaign using the 3 m Shane telescope at Lick Observatory, monitoring 15 low-redshift Seyfert 1 galaxies. This paper describes the observations, reductions and measurements, and data products from the spectroscopic campaign. The reduced spectra were fitted with a multicomponent model in order to isolate the contributions of various continuum and emission-line components. We present light curves of broad emission lines and the AGN continuum, and measurements of the broad H-beta line widths in mean and root-mean square (rms) spectra. For the most highly variable AGNs we also measured broad H-beta line widths and velocity centroids from the nightly spectra. In four AGNs exhibiting the highest variability amplitudes, we detect anticorrelations between broad H-beta width and luminosity, demonstrating that the broad-line region "breathes" on short timescales of days to weeks in response to continuum variations. We also find that broad H-beta velocity centroids can undergo substantial changes in response to continuum variations; in NGC 4593 the broad H-beta velocity shifted by ~250 km/s over a one-month duration. This reverberation-induced velocity shift effect is likely to contribute a significant source of confusion noise to binary black hole searches that use multi-epoch quasar spectroscopy to detect binary orbital motion. We also present results from simulations that examine biases that can occur in measurement of broad-line widths from rms spectra due to the contributions of continuum variations and photon-counting noise.

preprint2015arXiv

The Massive and Distant Clusters of WISE Survey. III. Sunyaev-Zel'dovich Masses of Galaxy Clusters at $z \sim 1$

We present CARMA 30 GHz Sunyaev-Zel'dovich (SZ) observations of five high-redshift ($z \gtrsim 1$), infrared-selected galaxy clusters discovered as part of the all-sky Massive and Distant Clusters of WISE Survey (MaDCoWS). The SZ decrements measured toward these clusters demonstrate that the MaDCoWS selection is discovering evolved, massive galaxy clusters with hot intracluster gas. Using the SZ scaling relation calibrated with South Pole Telescope clusters at similar masses and redshifts, we find these MaDCoWS clusters have masses in the range $M_{200} \approx 2-6 \times 10^{14}$ $M_\odot$. Three of these are among the most massive clusters found to date at $z\gtrsim 1$, demonstrating that MaDCoWS is sensitive to the most massive clusters to at least $z = 1.3$. The added depth of the AllWISE data release will allow all-sky infrared cluster detection to $z \approx 1.5$ and beyond.

preprint2015arXiv

The Multi-Layer Variable Absorbers in NGC 1365 Revealed by XMM-Newton and NuSTAR

Between July 2012 and February 2013, NuSTAR and XMM-Newton performed four long-look joint obser- vations of the type 1.8 Seyfert, NGC 1365. We have analyzed the variable absorption seen in these observations in order to characterize the geometry of the absorbing material. Two of the observations caught NGC 1365 in an unusually low absorption state, revealing complexity in the multi-layer absorber which had previously been hidden. We find the need for three distinct zones of neutral absorption in addition to the two zones of ionized absorption and the Compton-thick torus previously seen in this source. The most prominent absorber is likely associated with broad line region clouds with column densities of around $\sim\,$10$^{23}$ cm$^{-2}$ and a highly clumpy nature as evidenced by an occultation event in February 2013. We also find evidence of a patchy absorber with a variable column around $\sim\,10^{22}$ cm$^{-2}$ and a line of sight covering fraction of 0.3$-$0.9 which responds directly to the intrinsic source flux, possibly due to a wind geometry. A full-covering, constant absorber with a low column density of $\sim\,1 \times$ 10$^{22}$ cm$^{-2}$ is also present, though the location of this low density haze is unknown.

preprint2015arXiv

The NuSTAR Extragalactic Surveys: First Direct Measurements of the >10 keV X-Ray Luminosity Function for Active Galactic Nuclei at z>0.1

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.

preprint2015arXiv

The NuSTAR Extragalactic Surveys: Initial Results and Catalog from the Extended Chandra Deep Field South

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.

preprint2015arXiv

The NuSTAR Extragalactic Surveys: Overview and Catalog from the COSMOS Field

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.

preprint2015arXiv

The NuSTAR view of reflection and absorption in NGC 7582

NGC 7582 is a well-studied X-ray bright Seyfert 2 with moderately heavy ($N_{\text{H}}\sim10^{23}-10^{24}$~cm$^{-2}$), highly variable absorption and strong reflection spectral features. The spectral shape changed around the year 2000, dropping in observed flux and becoming much more highly absorbed. Two scenarios have been put forth to explain this spectral change: 1) the central X-ray source partially ``shut off'' around this time, decreasing in intrinsic luminosity, with a delayed decrease in reflection features due to the light-crossing time of the Compton-thick material or 2) the source became more heavily obscured, with only a portion of the power law continuum leaking through. NuSTAR observed NGC~7582 twice in 2012, two weeks apart, in order to quantify the reflection using high-quality data above 10 keV. We find that the most plausible scenario is that NGC 7582 has recently become more heavily absorbed by a patchy torus with a covering fraction of $\sim\,80-90\%$ and an equatorial column density of $\sim 3 \times10^{24}$ cm$^{-2}$. We find the need for an additional highly variable full-covering absorber with $N_{\text{H}}= 4-6 \times10^{23}$ cm$^{-2}$ in the line of sight, possibly associated with a hidden broad line region.

preprint2015arXiv

The NuSTAR X-ray spectrum of the low-luminosity AGN in NGC 7213

We present an analysis of the 3-79 keV NuSTAR spectrum of the low-luminosity active galactic nucleus NGC 7213. In agreement with past observations, we find a lower limit to the high-energy cut-off of Ec > 140 keV, no evidence for a Compton-reflected continuum, and the presence of an iron Kalpha complex, possibly produced in the broad-line region. From the application of the MYTorus model, we find that the line-emitting material is consistent with the absence of a significant Compton reflection if arising from a Compton-thin torus of gas with a column density of 5.0(+2.0,-1.6) cm^-2. We report variability of the equivalent width of the iron lines on the time-scale of years using archival observations from XMM-Newton, Chandra and Suzaku. This analysis suggests a possible contribution from dusty gas. A fit with a Comptonization model indicates the presence of a hot corona with a temperature kTe > 40 keV and an optical depth tau < 1, assuming a spherical geometry.

preprint2015arXiv

The Strikingly Uniform, Highly Turbulent Interstellar Medium of the Most Luminous Galaxy in the Universe

Observed at z = 4.601 and with L_bol = 3.5 x 10^14 Lsun, W2246-0526 is the most luminous galaxy known in the Universe, and hosts a deeply-buried active galactic nucleus (AGN)/super-massive black hole (SMBH). Discovered using the Wide-field Infrared Survey Explorer (WISE), W2246-0526 is classified as a Hot Dust Obscured Galaxy (Hot DOG), based on its luminosity and dust temperature. Here we present spatially resolved ALMA [CII]157.7um observations of W2246-0526, providing unique insight into the kinematics of its interstellar medium (ISM). The measured [CII]-to-far-infrared ratio is ~2 x 10^-4, implying ISM conditions that compare only with the most obscured, compact starbursts and AGN in the local Universe today. The spatially resolved [CII] line is strikingly uniform and very broad, 500-600 km/s wide, extending throughout the entire galaxy over about 2.5 kpc, with modest shear. Such a large, homogeneous velocity dispersion indicates a highly turbulent medium. W2246-0526 is unstable in terms of the energy and momentum that are being injected into the ISM, strongly suggesting that the gas is being blown away from the system isotropically, likely reflecting a cathartic state on its road to becoming an un-obscured quasar. W2246-0526 provides an extraordinary laboratory to study and model the properties and kinematics of gas in an extreme environment under strong feedback, at a time when the Universe was 1/10 of its current age: a system pushing the limits that can be reached during galaxy formation.

preprint2015arXiv

Wide-Field InfrarRed Survey Telescope-Astrophysics Focused Telescope Assets WFIRST-AFTA 2015 Report

This report describes the 2014 study by the Science Definition Team (SDT) of the Wide-Field Infrared Survey Telescope (WFIRST) mission. It is a space observatory that will address the most compelling scientific problems in dark energy, exoplanets and general astrophysics using a 2.4-m telescope with a wide-field infrared instrument and an optical coronagraph. The Astro2010 Decadal Survey recommended a Wide Field Infrared Survey Telescope as its top priority for a new large space mission. As conceived by the decadal survey, WFIRST would carry out a dark energy science program, a microlensing program to determine the demographics of exoplanets, and a general observing program utilizing its ultra wide field. In October 2012, NASA chartered a Science Definition Team (SDT) to produce, in collaboration with the WFIRST Study Office at GSFC and the Program Office at JPL, a Design Reference Mission (DRM) for an implementation of WFIRST using one of the 2.4-m, Hubble-quality telescope assemblies recently made available to NASA. This DRM builds on the work of the earlier WFIRST SDT, reported by Green et al. (2012) and the previous WFIRST-2.4 DRM, reported by Spergel et. (2013). The 2.4-m primary mirror enables a mission with greater sensitivity and higher angular resolution than the 1.3-m and 1.1-m designs considered previously, increasing both the science return of the primary surveys and the capabilities of WFIRST as a Guest Observer facility. The addition of an on-axis coronagraphic instrument to the baseline design enables imaging and spectroscopic studies of planets around nearby stars.

preprint2014arXiv

{\it NuSTAR} Reveals an Intrinsically X-ray Weak Broad Absorption Line Quasar in the Ultraluminous Infrared Galaxy Markarian 231

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 ($α_{\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.

preprint2014arXiv

A Focused, Hard X-ray Look at Arp 299 with NuSTAR

We report on simultaneous observations of the local starburst system Arp 299 with NuSTAR and Chandra, which provides the first resolved images of this galaxy up to energies of ~ 45 keV. Fitting the 3-40 keV spectrum reveals a column density of $N_{\rm H}$ ~ 4 x10^{24} cm^{-2}, characteristic of a Compton-thick AGN, and a 10-30 keV luminosity of 1.2x 10^{43} ergs s^{-1}. The hard X-rays detected by NuSTAR above 10 keV are centered on the western nucleus, Arp 299-B, which previous X-ray observations have shown to be the primary source of neutral Fe-K emission. Other X-ray sources, including Arp 299-A, the eastern nucleus which is also thought to harbor an AGN, as well as X-ray binaries, contribute $\lesssim 10%$ to the 10-20 keV emission from the Arp 299 system. The lack of significant emission above 10 keV other than that attributed to Arp 299-B suggests that: a) any AGN in Arp 299-A must be heavily obscured ($N_{\rm H}$ > 10^{24} cm^{-2}) or have a much lower luminosity than Arp 299-B and b) the extranuclear X-ray binaries have spectra that cut-off above ~10 keV. Such soft spectra are characteristic of ultraluminous X-ray (ULX) sources observed to date by NuSTAR.

preprint2014arXiv

An Ultraluminous X-ray Source Powered by An Accreting Neutron Star

Ultraluminous X-ray sources (ULX) are off-nuclear point sources in nearby galaxies whose X-ray luminosity exceeds the theoretical maximum for spherical infall (the Eddington limit) onto stellar-mass black holes. Their luminosity ranges from $10^{40}$ erg s$^{-1} < L_X$(0.5 - 10 keV) $<10^{40}$ erg s$^{-1}$. Since higher masses imply less extreme ratios of the luminosity to the isotropic Eddington limit theoretical models have focused on black hole rather than neutron star systems. The most challenging sources to explain are those at the luminous end ($L_X$ > $10^{40}$ erg s$^{-1}$), which require black hole masses MBH >50 solar masses and/or significant departures from the standard thin disk accretion that powers bright Galactic X-ray binaries. Here we report broadband X-ray observations of the nuclear region of the galaxy M82, which contains two bright ULXs. The observations reveal pulsations of average period 1.37 s with a 2.5-day sinusoidal modulation. The pulsations result from the rotation of a magnetized neutron star, and the modulation arises from its binary orbit. The pulsed flux alone corresponds to $L_X$(3 - 30 keV) = $4.9 \times 10^{39}$ erg s$^{-1}$. The pulsating source is spatially coincident with a variable ULX which can reach $L_X$ (0.3 - 10 keV) = $1.8 \times 10^{40}$ erg s$^{-1}$. This association implies a luminosity ~100 times the Eddington limit for a 1.4 solar mass object, or more than ten times brighter than any known accreting pulsar. This finding implies that neutron stars may not be rare in the ULX population, and it challenges physical models for the accretion of matter onto magnetized compact objects.

preprint2014arXiv

Asymmetries in core-collapse supernovae from maps of radioactive $^{44}$Ti in CassiopeiaA

Asymmetry is required by most numerical simulations of stellar core-collapse explosions, but the form it takes differs significantly among models. The spatial distribution of radioactive 44Ti, synthesized in an exploding star near the boundary between material falling back onto the collapsing core and that ejected into the surrounding medium, directly probes the explosion asymmetries. CassiopeiaA is a young, nearby, core-collapse remnant from which 44Ti emission has previously been detected but not imaged. Asymmetries in the explosion have been indirectly inferred from a high ratio of observed 44Ti emission to estimated 56Ni emission, from optical light echoes, and from jet-like features seen in the X-ray and optical ejecta. Here we report spatial maps and spectral properties of the 44Ti in Cassiopeia A. This may explain the unexpected lack of correlation between the 44Ti and iron X-ray emission, the latter being visible only in shock-heated material. The observed spatial distribution rules out symmetric explosions even with a high level of convective mixing, as well as highly asymmetric bipolar explosions resulting from a fast-rotating progenitor. Instead, these observations provide strong evidence for the development of low-mode convective instabilities in core-collapse supernovae.

preprint2014arXiv

Broadband X-ray Spectra of the Ultraluminous X-ray Source Holmberg IX X-1 observed with NuSTAR, XMM-Newton and Suzaku

We present results from the coordinated broadband X-ray observations of the extreme ultraluminous X-ray source Holmberg IX X-1 performed by NuSTAR, XMM-Newton and Suzaku in late 2012. These observations provide the first high-quality spectra of Holmberg IX X-1 above 10 keV to date, extending the X-ray coverage of this remarkable source up to ~30 keV. Broadband observations were undertaken at two epochs, between which Holmberg IX X-1 exhibited both flux and strong spectral variability, increasing in luminosity from Lx = (1.90+/-0.03)e40 erg/s to Lx = (3.35+/-0.03)e40 erg/s. Neither epoch exhibits a spectrum consistent with emission from the standard low/hard accretion state seen in Galactic black hole binaries, that would have been expected if Holmberg IX X-1 harbors a truly massive black hole accreting at substantially sub-Eddington accretion rates. The NuSTAR data confirm that the curvature observed previously in the 3-10 keV bandpass does represent a true spectral cutoff. During each epoch, the spectrum appears to be dominated by two optically thick thermal components, likely associated with an accretion disk. The spectrum also shows some evidence for a non-thermal tail at the highest energies, which may further support this scenario. The available data allow for either of the two thermal components to dominate the spectral evolution, although both scenarios require highly non-standard behavior for thermal accretion disk emission.

preprint2014arXiv

Candidate Clusters of Galaxies at $z>1.3$ Identified in the Spitzer SPT Deep Field Survey

We present 279 galaxy cluster candidates at $z > 1.3$ selected from the 94 deg$^{2}$ Spitzer South Pole Telescope Deep Field (SSDF) survey. We use a simple algorithm to select candidate high-redshift clusters of galaxies based on Spitzer/IRAC mid-infrared data combined with shallow all-sky optical data. We identify distant cluster candidates in SSDF adopting an overdensity threshold that results in a high purity (80%) cluster sample based on tests in the Spitzer Deep, Wide-Field Survey of the Boötes field. Our simple algorithm detects all three $1.4 < z \leq 1.75$ X-ray detected clusters in the Boötes field. The uniqueness of the SSDF survey resides not just in its area, one of the largest contiguous extragalactic fields observed with Spitzer, but also in its deep, multi-wavelength coverage by the South Pole Telescope (SPT), Herschel/SPIRE and XMM-Newton. This rich dataset will allow direct or stacked measurements of Sunyaev-Zel'dovich effect decrements or X-ray masses for many of the SSDF clusters presented here, and enable systematic study of the most distant clusters on an unprecedented scale. We measure the angular correlation function of our sample and find that these candidates show strong clustering. Employing the COSMOS/UltraVista photometric catalog in order to infer the redshift distribution of our cluster selection, we find that these clusters have a comoving number density $n_c = (0.7^{+6.3}_{-0.6}) \times 10^{-7} h^{3} \mathrm{Mpc}^{-3}$ and a spatial clustering correlation scale length $r_0 = (32 \pm 7) h^{-1} \rm{Mpc}$. Assuming our sample is comprised of dark matter halos above a characteristic minimum mass, $M_{\rm min}$, we derive that at $z=1.5$ these clusters reside in halos larger than $M_{\rm min} = 1.5^{+0.9}_{-0.7} \times 10^{14} h^{-1} M_{\odot}$. (abridged)

preprint2014arXiv

Coronal Properties of the Seyfert 1.9 Galaxy MCG -05-23-016 Determined from Hard X-ray Spectroscopy with NuSTAR

Measurements of the high-energy cut-off in the coronal continuum of active galactic nuclei have long been elusive for all but a small number of the brightest examples. We present a direct measurement of the cut-off energy in the nuclear continuum of the nearby Seyfert 1.9 galaxy MCG -05-23-016 with unprecedented precision. The high sensitivity of NuSTAR up to 79 keV allows us to clearly disentangle the spectral curvature of the primary continuum from that of its reflection component. Using a simple phenomenological model for the hard X-ray spectrum, we constrain the cut-off energy to $116_{-5}^{+6}$ keV with 90% confidence. Testing for more complex models and nuisance parameters that could potentially influence the measurement, we find that the cut-off is detected robustly. We further use simple Comptonized plasma models to provide independent constraints for both the kinetic temperature of the electrons in the corona and its optical depth. At the 90% confidence level, we find $kT_e=29\pm2$ keV and $τ_e=1.23\pm0.08$ assuming a slab (disk-like) geometry, and $kT_e=25\pm2$ keV and $τ_e=3.5\pm0.2$ assuming a spherical geometry. Both geometries are found to fit the data equally well and their two principal physical parameters are correlated in both cases. With the optical depth in the optically thick regime, the data are pushing the currently available theoretical models of the Comptonized plasma to the limits of their validity. Since the spectral features and variability arising from the inner accretion disk have been observed previously in MCG -05-23-016, the inferred high optical depth implies that a spherical or disk-like corona cannot be homogeneous.

preprint2014arXiv

Iron K and Compton hump reverberation in SWIFT J2127.4+5654 and NGC 1365 revealed by NuSTAR and XMM-Newton

In the past five years, a flurry of X-ray reverberation lag measurements of accreting supermassive black holes have been made using the XMM-Newton telescope in the 0.3-10 keV energy range. In this work, we use the NuSTAR telescope to extend the lag analysis up to higher energies for two Seyfert galaxies, SWIFT J2127.4+5654 and NGC 1365. X-ray reverberation lags are due to the light travel time delays between the direct continuum emission and the reprocessed emission from the inner radii of an ionised accretion disc. XMM-Newton has been particularly adept at measuring the lag associated with the broad Fe K emission line, where the gravitationally redshifted wing of the line is observed to respond before the line centroid at 6.4 keV, produced at larger radii. Now we use NuSTAR to probe the lag at higher energies, where the spectrum shows clear evidence for Compton reflection, known as the Compton 'hump'. The XMM-Newton data show Fe K lags in both SWIFT J2127.4+5654 and NGC 1365. The NuSTAR data provide independent confirmation of these Fe K lags, and also show evidence for the corresponding Compton hump lags, especially in SWIFT J2127.4+5654. These broadband lag measurements confirm that the Compton hump and Fe K lag are produced at small radii. At low-frequencies in NGC 1365, where the spectrum shows evidence for eclipsing clouds in the line of sight, we find a clear negative (not positive) lag from 2--10 keV, which can be understood as the decrease in column density from a neutral eclipsing cloud moving out of our line of sight during the observation.

preprint2014arXiv

NuSTAR and XMM-Newton Observations of Luminous, Heavily Obscured, WISE-Selected Quasars at z ~ 2

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.

preprint2014arXiv

NuSTAR and XMM-Newton observations of the extreme ultraluminous X-ray source NGC 5907 ULX1: A Vanishing Act

We present results obtained from two broadband X-ray observations of the extreme ultraluminous X-ray source (ULX) NGC5907 ULX1, known to have a peak X-ray luminosity of ~5e40 erg/s. These XMM-Newton and NuSTAR observations, separated by only ~4 days, revealed an extreme level of short-term flux variability. In the first epoch, NGC5907 ULX1 was undetected by NuSTAR, and only weakly detected (if at all) with XMM-Newton, while in the second NGC5907 ULX1 was clearly detected at high luminosity by both missions. This implies an increase in flux of ~2 orders of magnitude or more during this ~4 day window. We argue that this is likely due to a rapid rise in the mass accretion rate, rather than to a transition from an extremely obscured to an unobscured state. During the second epoch we observed the broadband 0.3-20.0 keV X-ray luminosity to be (1.55+/-0.06)e40 erg/s, similar to the majority of the archival X-ray observations. The broadband X-ray spectrum obtained from the second epoch is inconsistent with the low/hard accretion state observed in Galactic black hole binaries, but is well modeled with a simple accretion disk model incorporating the effects of photon advection. This strongly suggests that, when bright, NGC5907 ULX1 is a high-Eddington accretor.

preprint2014arXiv

NuSTAR J033202-2746.8: direct constraints on the Compton reflection in a heavily obscured quasar at z~2

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.

preprint2014arXiv

NuSTAR Observations of Heavily Obscured Quasars at z ~ 0.5

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.

preprint2014arXiv

NuSTAR Reveals the Comptonizing Corona of the Broad-Line Radio Galaxy 3C 382

Broad-line radio galaxies (BLRGs) are active galactic nuclei that produce powerful, large-scale radio jets, but appear as Seyfert 1 galaxies in their optical spectra. In the X-ray band, BLRGs also appear like Seyfert galaxies, but with flatter spectra and weaker reflection features. One explanation for these properties is that the X-ray continuum is diluted by emission from the jet. Here, we present two NuSTAR observations of the BLRG 3C 382 that show clear evidence that the continuum of this source is dominated by thermal Comptonization, as in Seyfert 1 galaxies. The two observations were separated by over a year and found 3C 382 in different states separated by a factor of 1.7 in flux. The lower flux spectrum has a photon-index of $Γ=1.68^{+0.03}_{-0.02}$, while the photon-index of the higher flux spectrum is $Γ=1.78^{+0.02}_{-0.03}$. Thermal and anisotropic Comptonization models provide an excellent fit to both spectra and show that the coronal plasma cooled from $kT_e=330\pm 30$ keV in the low flux data to $231^{+50}_{-88}$ keV in the high flux observation. This cooling behavior is typical of Comptonizing corona in Seyfert galaxies and is distinct from the variations observed in jet-dominated sources. In the high flux observation, simultaneous Swift data are leveraged to obtain a broadband spectral energy distribution and indicates that the corona intercepts $\sim 10$% of the optical and ultraviolet emitting accretion disk. 3C 382 exhibits very weak reflection features, with no detectable relativistic Fe K$α$ line, that may be best explained by an outflowing corona combined with an ionized inner accretion disk.

preprint2014arXiv

NuSTAR unveils a Compton-thick Type 2 quasar in Mrk 34

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.

preprint2014arXiv

Observations of MCG-5-23-16 with Suzaku, XMM-Newton and NuSTAR: Disk tomography and Compton hump reverberation

MCG-5-23-16 is one of the first AGN where relativistic reverberation in the iron K line originating in the vicinity of the supermassive black hole was found, based on a short XMM-Newton observation. In this work, we present the results from long X-ray observations using Suzaku, XMM-Newton and NuSTAR designed to map the emission region using X-ray reverberation. A relativistic iron line is detected in the lag spectra on three different time-scales, allowing the emission from different regions around the black hole to be separated. Using NuSTAR coverage of energies above 10 keV reveals a lag between these energies and the primary continuum, which is detected for the first time in an AGN. This lag is a result of the Compton reflection hump responding to changes in the primary source in a manner similar to the response of the relativistic iron K line.

preprint2014arXiv

Rapidly growing black holes and host galaxies in the distant Universe from the Herschel Radio Galaxy Evolution Project

We present results from a survey of 70 radio galaxies (RGs) at redshifts 1<z<5.2 using the PACS and SPIRE on-board Herschel. Combined with existing mid-IR photometry from Spitzer and observations obtained with LABOCA, the SEDs of galaxies in our sample are continuously covered across 3.6-870um. The total infrared luminosities of these RGs are such that they almost all are either ultra-or hyper-luminous infrared galaxies. We fit the infrared SEDs with a set of empirical templates which represent dust heated (1) by a variety of SB and (2) by a AGN. We find that the SEDs of RGs require the dust to be heated by both AGN and SB, but the luminosities of these two components are not strongly correlated. Assuming empirical relations and simple physical assumptions, we calculate the SFR, the black hole mass accretion rate (MdotBH), and the black hole mass (MBH) for each RG. We find that the host galaxies and their BHs are growing extremely rapidly, having SFR~100-5000 Msun/yr and MdotBH~1-100 Msun/yr. The mean sSFR of RGs at z>2.5 are higher than the sSFR of typical star-forming galaxies over the same redshift range but are similar or perhaps lower than the galaxy population for RGs at z<2.5. By comparing the sSFR and the specific black hole mass accretion rate, we conclude that BHs in radio loud AGN are already, or soon will be, overly massive compared to their host galaxies in terms of expectations from the local MBH-MGal relation. In order to ``catch up'' with the BH, the galaxies require about an order-of magnitude more time to grow in mass, at the observed SFRs, compared to the time the BH is actively accreting. However, during the current cycle of activity, we argue that this catching-up is likely to be difficult due to the short gas depletion times. Finally, we speculate on how the host galaxies might grow sufficiently in stellar mass to ultimately fall onto the local MBH-MGal relation.

preprint2014arXiv

The 2-79 keV X-ray Spectrum of the Circinus Galaxy with NuSTAR, XMM-Newton and Chandra: a Fully Compton-Thick AGN

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 $Γ= 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-$Γ$ 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.

preprint2014arXiv

The Angular Clustering of WISE-Selected AGN: Different Haloes for Obscured and Unobscured AGN

We calculate the angular correlation function for a sample of 170,000 AGN extracted from the Wide-field Infrared Survey Explorer (WISE) catalog, selected to have red mid-IR colors (W1 - W2 > 0.8) and 4.6 micron flux densities brighter than 0.14 mJy). The sample is expected to be >90% reliable at identifying AGN, and to have a mean redshift of z=1.1. In total, the angular clustering of WISE-AGN is roughly similar to that of optical AGN. We cross-match these objects with the photometric SDSS catalog and distinguish obscured sources with (r - W2) > 6 from bluer, unobscured AGN. Obscured sources present a higher clustering signal than unobscured sources. Since the host galaxy morphologies of obscured AGN are not typical red sequence elliptical galaxies and show disks in many cases, it is unlikely that the increased clustering strength of the obscured population is driven by a host galaxy segregation bias. By using relatively complete redshift distributions from the COSMOS survey, we find obscured sources at mean redshift z=0.9 have a bias of b = 2.9 \pm 0.6 and are hosted in dark matter halos with a typical mass of log(M/M_odot)~13.5. In contrast, unobscured AGN at z~1.1 have a bias of b = 1.6 \pm 0.6 and inhabit halos of log(M/M_odot)~12.4. These findings suggest that obscured AGN inhabit denser environments than unobscured AGN, and are difficult to reconcile with the simplest AGN unification models, where obscuration is driven solely by orientation.

preprint2014arXiv

The Broad-band X-ray Spectrum of IC 4329A from a Joint NuSTAR/Suzaku Observation

We have obtained a deep, simultaneous observation of the bright, nearby Seyfert galaxy IC 4329A with Suzaku and NuSTAR. Through a detailed spectral analysis, we are able to robustly separate the continuum, absorption and distant reflection components in the spectrum. The absorbing column is found to be modest at $N_H = 6 \times 10^{21}$ cm$^2$, and does not introduce any significant curvature in the Fe K band. We are able to place a strong constraint on the presence of a broadened Fe Kα line: $E = 6.46^{+0.08}_{-0.07}$ keV rest frame with $σ = 0.33^{+0.08}_{-0.07}$ keV and $EW = 34^{+8}_{-7}$ eV, though we are not able to constrain any of the parameters of a relativistic reflection model. These results highlight the range in broad Fe Kα line strengths observed in nearby, bright AGN (roughly an order of magnitude), and imply a corresponding range in the physical properties of the inner accretion disk in these sources. We have also updated our previously reported measurement of the high-energy cutoff of the hard X-ray emission using both observatories rather than just NuSTAR alone: $E_{cut} = 186 \pm 14$ keV. This high-energy cutoff acts as a proxy for the temperature of the coronal electron plasma, enabling us to further separate this parameter from the optical depth of the plasma and to update our results for these parameters as well. We derive $kT = 50^{+6}_{-3}$ keV with $τ = 2.34^{+0.16}_{-0.11}$ using a spherical geometry, $kT = 61 \pm 1$ keV with $τ = 0.68 \pm 0.02$ for a slab geometry, with both having an equivalent goodness-of-fit.

preprint2014arXiv

The hard X-ray spectrum of NGC 5506 as seen by NuSTAR

NuSTAR observed the bright Compton-thin, narrow line Seyfert 1 galaxy, NGC 5506, for about 56 ks. In agreement with past observations, the spectrum is well fit by a power law with Gamma~1.9, a distant reflection component and narrow ionized iron lines. A relativistically blurred reflection component is not required by the data. When an exponential high energy cutoff is added to the power law, a value of 720(+130,-190) keV (90% confidence level) is found. Even allowing for systematic uncertainties, we find a 3 sigma lower limit to the high-energy cutoff of 350 keV, the highest lower limit to the cutoff energy found so far in an AGN by NuSTAR.

preprint2014arXiv

The NuSTAR spectrum of Mrk 335: Extreme relativistic effects within 2 gravitational radii of the event horizon?

We present 3-50 keV NuSTAR observations of the AGN Mrk 335 in a very low flux state. The spectrum is dominated by very strong features at the energies of the iron line at 5-7 keV and Compton hump from 10-30 keV. The source is variable during the observation, with the variability concentrated at low energies, which suggesting either a relativistic reflection or a variable absorption scenario. In this work we focus on the reflection interpretation, making use of new relativistic reflection models that self consistently calculate the reflection fraction, relativistic blurring and angle-dependent reflection spectrum for different coronal heights to model the spectra. We find that the spectra can be well fit with relativistic reflection, and that the lowest flux state spectrum is described by reflection alone, suggesting the effects of extreme light-bending occurring within ~2 gravitational radii of the event horizon. The reflection fraction decreases sharply with increasing flux, consistent with a point source moving up to above 10 Rg as the source brightens. We constrain the spin parameter to greater than 0.9 at the 3 sigma confidence level. By adding a spin-dependent upper limit on the reflection fraction to our models, we demonstrate that this can be a powerful way of constraining the spin parameter, particularly in reflection dominated states. We also calculate a detailed emissivity profile for the iron line, and find that it closely matches theoretical predictions for a compact source within a few Rg of the black hole.

preprint2014arXiv

The NuSTAR View of Nearby Compton-thick AGN: The Cases of NGC 424, NGC 1320 and IC 2560

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.

preprint2014arXiv

The Seyfert 2 galaxy NGC 2110: hard X-ray emission observed by NuSTAR and variability of the iron K$α$ line

We present NuSTAR observations of the bright Seyfert 2 galaxy NGC 2110 obtained in 2012, when the source was at the highest flux level ever observed, and in 2013, when the source was at a more typical flux level. We include archival observations from other X-ray satellites, namely XMM-Newton, Suzaku, BeppoSAX, Chandra and Swift. Simultaneous NuSTAR and Swift broad band spectra (in the 3-80 keV range) indicate a cutoff energy $E_{\rm c}>210$ keV, with no detectable contribution from Compton reflection. NGC 2110 is one of the very few sources where no evidence for distant Compton thick scattering is found and, by using temporal information collected over more than a decade, we investigate variations of the iron K$α$ line on time scales of years. The Fe K$α$ line is likely the sum of two components: one constant (originating from distant Compton-thick material) and the other one variable and linearly correlated with the source flux (possibly arising from Compton-thin material much closer to the black hole).

preprint2014arXiv

The soft X-ray emission of Ark 120. XMM-Newton, NuSTAR, and the importance of taking the broad view

We present simultaneous XMM-Newton and NuSTAR observations of the `bare' Seyfert 1 galaxy, Ark 120, a system in which ionized absorption is absent. The NuSTAR hard X-ray spectral coverage allows us to constrain different models for the excess soft X-ray emission. Among phenomenological models, a cutoff power law best explains the soft X-ray emission. This model likely corresponds to Comptonization of the accretion disk seed UV photons by a population of warm electrons: using Comptonization models, a temperature of ~0.3 keV and an optical depth of ~13 are found. If the UV-to-X-ray optxagnf model is applied, the UV fluxes from the XMM-$Newton$ Optical Monitor suggest an intermediate black hole spin. Contrary to several other sources observed by NuSTAR, no high energy cutoff is detected, with a lower limit of 190 keV.

preprint2014arXiv

Timing and Flux Evolution of the Galactic Center Magnetar SGR J1745-2900

We present the X-ray timing and spectral evolution of the Galactic Center magnetar SGR J1745-2900 for the first ~4 months post-discovery using data obtained with the Nuclear Spectroscopic Telescope Array (NuSTAR)} and Swift observatories. Our timing analysis reveals a large increase in the magnetar spin-down rate by a factor of 2.60 +/- 0.07 over our data span. We further show that the change in spin evolution was likely coincident with a bright X-ray burst observed in 2013 June by Swift, and if so, there was no accompanying discontinuity in the frequency. We find that the source 3-10 keV flux has declined monotonically by a factor of ~2 over an 80-day period post-outburst accompanied by a ~20% decrease in the source's blackbody temperature, although there is evidence for both flux and kT having levelled off. We argue that the torque variations are likely to be magnetospheric in nature and will dominate over any dynamical signatures of orbital motion around Sgr A*.

preprint2014arXiv

Weak Hard X-ray Emission from Broad Absorption Line Quasars: Evidence for Intrinsic X-ray Weakness

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 (Γ~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.

preprint2014arXiv

Why z > 1 radio-loud galaxies are commonly located in proto-clusters

Distant powerful radio-loud active galactic nuclei (RLAGN) tend to reside in dense environments and are commonly found in proto-clusters at z > 1.3. We examine whether this occurs because RLAGN are hosted by massive galaxies, which preferentially reside in rich environments. We compare the environments of powerful RLAGN at 1.3 < z < 3.2 from the CARLA survey to a sample of radio-quiet galaxies matched in mass and redshift. We find the environments of RLAGN are significantly denser than those of radio-quiet galaxies, implying that not more than 50% of massive galaxies in this epoch can host powerful radio-loud jets. This is not an observational selection effect as we find no evidence to suggest it is easier to observe the radio emission when the galaxy resides in a dense environment. We therefore suggest that the dense Mpc-scale environment fosters the formation of a radio-jet from an AGN. We show that the number density of potential RLAGN host galaxies is consistent with every > 10^14 solar mass cluster having experienced powerful radio-loud feedback of duration ~60 Myr during 1.3 < z < 3.2. This feedback could heat the intracluster medium to the extent of 0.5-1 keV per gas particle, which could limit the amount of gas available for further star formation in the proto-cluster galaxies.

preprint2013arXiv

An Extremely Luminous and Variable Ultraluminous X-ray Source in the Outskirts of Circinus Observed with NuSTAR

Following a serendipitous detection with the NuSTAR observatory, we present a multi-epoch spectral and temporal analysis of an extreme ultraluminous X-ray source (ULX) located in the outskirts of the Circinus galaxy, hereafter Circinus ULX5, including coordinated follow-up observations with XMM-Newton and NuSTAR. The NuSTAR data presented here represent one of the first instances of a ULX reliably detected at hard (E>10 keV) X-rays. Circinus ULX5 is variable on long timescales by at least a factor of ~5 in flux, and was caught in a historically bright state during our 2013 observations (0.3-30.0 keV luminosity of 1.6e40 erg/s). During this epoch, the source displayed a curved 3-10 keV spectrum, broadly similar to other bright ULXs. While in many cases the interpretation of spectral components in ULXs is uncertain, the spectral and temporal properties of the all the high quality datasets currently available strongly support a simple disk-corona model reminiscent of that invoked for Galactic binaries, with the accretion disk becoming more prominent as the luminosity increases. However, although the disk temperature and luminosity are remarkably well correlated across all timescales currently probed, the observed luminosity follows T^(1.70+/-0.17), flatter than that expected for simple blackbody radiation. The spectral variability displayed by Circinus ULX5 is highly reminiscent of that observed from some known Galactic black hole binaries (BHBs) at high luminosities. This comparison would imply a black hole mass of ~90 Msun for Circinus ULX5. However, given the diverse behavior observed from Galactic BHB accretion disks, this mass estimate is still uncertain. Finally, the limits placed on any undetected iron absorption features with the 2013 dataset imply that we are not viewing the central regions of Circinus ULX5 through any extreme super-Eddington outflow. [abridged]

preprint2013arXiv

First Results from NuSTAR Observations of Mkn 421

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.

preprint2013arXiv

Galaxy Clusters around radio-loud AGN at 1.3 < z < 3.2 as seen by Spitzer

We report the first results from the Clusters Around Radio-Loud AGN (CARLA) program, a Cycle 7 and 8 Spitzer Space Telescope snapshot program to investigate the environments of a large sample of obscured and unobscured luminous radio-loud AGN at 1.2 < z < 3.2. These data, obtained for 387 fields, reach 3.6 and 4.5 um depths of [3.6] (AB) = 22.6 and [4.5] (AB) = 22.9 at the 95% completeness level, which is two to three times fainter than L* in this redshift range. By using the color cut [3.6]-[4.5] > -0.1 (AB), which efficiently selects high-redshift (z > 1.3) galaxies of all types, we identify galaxy cluster member candidates in the fields of the radio-loud AGN. The local density of these IRAC-selected sources is compared to the density of similarly selected sources in blank fields. We find that 92% of the radio-loud AGN reside in environments richer than average. The majority (55%) of the radio-loud AGN fields are found to be overdense at a > 2 σ level; 10% are overdense at a > 5 σ level. A clear rise in surface density of IRAC-selected sources towards the position of the radio-loud AGN strongly supports an association of the majority of the IRAC-selected sources with the radio-loud AGN. Our results provide solid statistical evidence that radio-loud AGN are likely beacons for finding high-redshift galaxy (proto-)clusters. We investigate how environment depends on AGN type (unobscured radio-loud quasars vs. obscured radio galaxies), radio luminosity and redshift, finding no correlation with either AGN type or radio luminosity. We find a decrease in density with redshift, consistent with galaxy evolution for this uniform, flux-limited survey. These results are consistent with expectations from the orientation-driven AGN unification model, at least for the high radio luminosity regimes considered in this sample.

preprint2013arXiv

Improved constraints on the expansion rate of the Universe up to z~1.1 from the spectroscopic evolution of cosmic chronometers

We present new improved constraints on the Hubble parameter H(z) in the redshift range 0.15 < z < 1.1, obtained from the differential spectroscopic evolution of early-type galaxies as a function of redshift. We extract a large sample of early-type galaxies (\sim11000) from several spectroscopic surveys, spanning almost 8 billion years of cosmic lookback time (0.15 < z < 1.42). We select the most massive, red elliptical galaxies, passively evolving and without signature of ongoing star formation. Those galaxies can be used as standard cosmic chronometers, as firstly proposed by Jimenez & Loeb (2002), whose differential age evolution as a function of cosmic time directly probes H(z). We analyze the 4000 Å break (D4000) as a function of redshift, use stellar population synthesis models to theoretically calibrate the dependence of the differential age evolution on the differential D4000, and estimate the Hubble parameter taking into account both statistical and systematical errors. We provide 8 new measurements of H(z) (see Tab. 4), and determine its change in H(z) to a precision of 5-12% mapping homogeneously the redshift range up to z \sim 1.1; for the first time, we place a constraint on H(z) at z \neq 0 with a precision comparable with the one achieved for the Hubble constant (about 5-6% at z \sim 0.2), and covered a redshift range (0.5 < z < 0.8) which is crucial to distinguish many different quintessence cosmologies. These measurements have been tested to best match a ΛCDM model, clearly providing a statistically robust indication that the Universe is undergoing an accelerated expansion. This method shows the potentiality to open a new avenue in constrain a variety of alternative cosmologies, especially when future surveys (e.g. Euclid) will open the possibility to extend it up to z \sim 2.

preprint2013arXiv

Measuring the Coronal Properties of IC4329A with NuSTAR

We present an analysis of a ~160 ks NuSTAR observation of the nearby bright Seyfert galaxy IC4329A. The high-quality broadband spectrum enables us to separate the effects of distant reflection from the direct coronal continuum, and to therefore accurately measure the high-energy cutoff to be $E_{cut}=178^{+74}_{-40}$ keV. The coronal emission arises from accretion disk photons Compton up-scattered by a thermal plasma, with the spectral index and cutoff being due to a combination of the finite plasma temperature and optical depth. Applying standard Comptonization models, we measure both physical properties independently using the best signal-to-noise obtained to date in an AGN over the 3-79 keV band. We derive $kT_e=37^{+7}_{-6}$ keV with $τ=1.25^{+0.20}_{-0.10}$ assuming a slab geometry for the plasma, and $kT_e=33^{+6}_{-6}$ keV with $τ=3.41^{+0.58}_{-0.38}$ for a spherical geometry, with both having an equivalent goodness-of-fit.

preprint2013arXiv

NuSTAR detection of the blazar B2 1023+25 at redshift 5.3

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.

preprint2013arXiv

SN 2010jl: Optical to hard X-ray observations reveal an explosion embedded in a ten solar mass cocoon

(Abridged) Some supernovae (SNe) may be powered by the interaction of the SN ejecta with a large amount of circumstellar matter (CSM). Here we outline a method to measure the mass of the optically thick CSM around such SNe. We present observations of SN2010jl, including the first detection of a SN using NuSTAR. The total radiated luminosity of SN2010jl is extreme, at least 9e50 erg. By modeling the visible-light data, we robustly show that the mass of the circumstellar material within ~1e16 cm of the progenitor was in excess of 10 solar masses, likely ejected tens of years prior to the SN explosion. Our modeling suggests that the shock velocity during shock breakout was ~6000 km/s, decelerating to ~2600 km/s about two years after maximum light. Our late-time NuSTAR+XMM spectra of the SN presumably provide the first direct measurement of SN shock velocity two years after the SN maximum light -- measured to be in the range of 2000 to 4500 km/s if the ions and electrons are in equilibrium, and >~2000 km/s if they are not in equilibrium. This measurement is in agreement with the shock velocity predicted by our modeling of the optical data. We also show that the mean radial density distribution of the CSM roughly follows an r^-2 law. A possible explanation for the massive CSM with a wind-like profile is that they are the result of multiple pulsational pair instability events prior to the SN explosion, separated from each other by years.

preprint2013arXiv

Spectroscopic Needs for Imaging Dark Energy Experiments: Photometric Redshift Training and Calibration

Large sets of objects with spectroscopic redshift measurements will be needed for imaging dark energy experiments to achieve their full potential, serving two goals:_training_, i.e., the use of objects with known redshift to develop and optimize photometric redshift algorithms; and_calibration_, i.e., the characterization of moments of redshift (or photo-z error) distributions. Better training makes cosmological constraints from a given experiment stronger, while highly-accurate calibration is needed for photo-z systematics not to dominate errors. In this white paper, we investigate the required scope of spectroscopic datasets which can serve both these purposes for ongoing and next-generation dark energy experiments, as well as the time required to obtain such data with instruments available in the next decade. Large time allocations on kilo-object spectrographs will be necessary, ideally augmented by infrared spectroscopy from space. Alternatively, precision calibrations could be obtained by measuring cross-correlation statistics using samples of bright objects from a large baryon acoustic oscillation experiment such as DESI. We also summarize the additional work on photometric redshift methods needed to prepare for ongoing and future dark energy experiments.

preprint2013arXiv

The Cluster and Field Galaxy AGN Fraction at z = 1 to 1.5: Evidence for a Reversal of the Local Anticorrelation Between Environment and AGN Fraction

The fraction of cluster galaxies that host luminous AGN is an important probe of AGN fueling processes, the cold ISM at the centers of galaxies, and how tightly black holes and galaxies co-evolve. We present a new measurement of the AGN fraction in a sample of 13 clusters of galaxies (M >= 10^{14} Msun) at 1<z<1.5 selected from the Spitzer/IRAC Shallow Cluster Survey, as well as the field fraction in the immediate vicinity of these clusters, and combine these data with measurements from the literature to quantify the relative evolution of cluster and field AGN from the present to z~3. We estimate that the cluster AGN fraction at 1<z<1.5 is f_A = 3.0^{+2.4}_{-1.4}% for AGN with a rest-frame, hard X-ray luminosity greater than L_{X,H} >= 10^{44} erg/s. This fraction is measured relative to all cluster galaxies more luminous than M*_{3.6}(z)+1, where M*_{3.6}(z) is the absolute magnitude of the break in the galaxy luminosity function at the cluster redshift in the IRAC 3.6um bandpass. The cluster AGN fraction is 30 times greater than the 3sigma upper limit on the value for AGN of similar luminosity at z~0.25, as well as more than an order of magnitude greater than the AGN fraction at z~0.75. AGN with L_{X,H} >= 10^{43} erg/s exhibit similarly pronounced evolution with redshift. In contrast with the local universe, where the luminous AGN fraction is higher in the field than in clusters, the X-ray and MIR-selected AGN fractions in the field and clusters are consistent at 1<z<1.5. This is evidence that the cluster AGN population has evolved more rapidly than the field population from z~1.5 to the present. This environment-dependent AGN evolution mimics the more rapid evolution of star-forming galaxies in clusters relative to the field.

preprint2013arXiv

The Era of Star Formation in Galaxy Clusters

We analyze the star formation properties of 16 infrared-selected, spectroscopically confirmed galaxy clusters at $1 < z < 1.5$ from the Spitzer/IRAC Shallow Cluster Survey (ISCS). We present new spectroscopic confirmation for six of these high-redshift clusters, five of which are at $z>1.35$. Using infrared luminosities measured with deep Spitzer/MIPS observations at 24 $μ$m, along with robust optical+IRAC photometric redshifts and SED-fitted stellar masses, we present the dust-obscured star-forming fractions, star formation rates and specific star formation rates in these clusters as functions of redshift and projected clustercentric radius. We find that $z\sim 1.4$ represents a transition redshift for the ISCS sample, with clear evidence of an unquenched era of cluster star formation at earlier times. Beyond this redshift the fraction of star-forming cluster members increases monotonically toward the cluster centers. Indeed, the specific star formation rate in the cores of these distant clusters is consistent with field values at similar redshifts, indicating that at $z>1.4$ environment-dependent quenching had not yet been established in ISCS clusters. Combining these observations with complementary studies showing a rapid increase in the AGN fraction, a stochastic star formation history, and a major merging episode at the same epoch in this cluster sample, we suggest that the starburst activity is likely merger-driven and that the subsequent quenching is due to feedback from merger-fueled AGN. The totality of the evidence suggests we are witnessing the final quenching period that brings an end to the era of star formation in galaxy clusters and initiates the era of passive evolution.

preprint2013arXiv

The NuSTAR Extragalactic Survey: A First Sensitive Look at the High-Energy Cosmic X-ray Background Population

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

preprint2013arXiv

The Spitzer-South Pole Telescope Deep Field: Survey Design and IRAC Catalogs

The Spitzer-South Pole Telescope Deep Field (SSDF) is a wide-area survey using Spitzer's Infrared Array Camera (IRAC) to cover 94 square degrees of extragalactic sky, making it the largest IRAC survey completed to date outside the Milky Way midplane. The SSDF is centered at 23:30,-55:00, in a region that combines observations spanning a broad wavelength range from numerous facilities. These include millimeter imaging from the South Pole Telescope, far-infrared observations from Herschel/SPIRE, X-ray observations from the XMM XXL survey, near-infrared observations from the VISTA Hemisphere Survey, and radio-wavelength imaging from the Australia Telescope Compact Array, in a panchromatic project designed to address major outstanding questions surrounding galaxy clusters and the baryon budget. Here we describe the Spitzer/IRAC observations of the SSDF, including the survey design, observations, processing, source extraction, and publicly available data products. In particular, we present two band-merged catalogs, one for each of the two warm IRAC selection bands. They contain roughly 5.5 and 3.7 million distinct sources, the vast majority of which are galaxies, down to the SSDF 5-sigma sensitivity limits of 19.0 and 18.2 Vega mag (7.0 and 9.4 microJy) at 3.6 and 4.5 microns, respectively.

preprint2013arXiv

Weak Hard X-ray Emission from Two Broad Absorption Line Quasars Observed with NuSTAR: Compton-thick Absorption or Intrinsic X-ray Weakness?

We present NuSTAR hard X-ray observations of two X-ray weak broad absorption line (BAL) quasars, PG 1004+130 (radio loud) and PG 1700+518 (radio quiet). Many BAL quasars appear X-ray weak, probably due to absorption by the shielding gas between the nucleus and the accretion-disk wind. The two targets are among the optically brightest BAL quasars, yet they are known to be significantly X-ray weak at rest-frame 2-10 keV (16-120 times fainter than typical quasars). We would expect to obtain ~400-600 hard X-ray (>10 keV) photons with NuSTAR, provided that these photons are not significantly absorbed (NH<1E24 cm^{-2}). However, both BAL quasars are only detected in the softer NuSTAR bands (e.g., 4-20 keV) but not in its harder bands (e.g., 20-30 keV), suggesting that either the shielding gas is highly Compton-thick or the two targets are intrinsically X-ray weak. We constrain the column densities for both to be NH~7E24 cm^{-2} if the weak hard X-ray emission is caused by obscuration from the shielding gas. We discuss a few possibilities for how PG 1004+130 could have Compton-thick shielding gas without strong Fe K? line emission; dilution from jet-linked X-ray emission is one likely explanation. We also discuss the intrinsic X-ray weakness scenario based on a coronal-quenching model relevant to the shielding gas and disk wind of BAL quasars. Motivated by our NuSTAR results, we perform a Chandra stacking analysis with the Large Bright Quasar Survey BAL quasar sample and place statistical constraints upon the fraction of intrinsically X-ray weak BAL quasars; this fraction is likely 17-40%.

preprint2012arXiv

A Measurement of the Correlation of Galaxy Surveys with CMB Lensing Convergence Maps from the South Pole Telescope

We compare cosmic microwave background lensing convergence maps derived from South Pole Telescope (SPT) data with galaxy survey data from the Blanco Cosmology Survey, the Wide-field Infrared Survey Explorer, and a new large Spitzer/IRAC field designed to overlap with the SPT survey. Using optical and infrared catalogs covering between 17 and 68 square degrees of sky, we detect correlation between the SPT convergence maps and each of the galaxy density maps at >4 sigma, with zero cross-correlation robustly ruled out in all cases. The amplitude and shape of the cross-power spectra are in good agreement with theoretical expectations and the measured galaxy bias is consistent with previous work. The detections reported here utilize a small fraction of the full 2500 square degree SPT survey data and serve as both a proof of principle of the technique and an illustration of the potential of this emerging cosmological probe.

preprint2012arXiv

Characterizing the Mid-IR Extragalactic Sky with WISE and SDSS

The Wide-field Infrared Survey Explorer (WISE) has completed its all-sky survey at 3.4-22 micron. We merge the WISE data with optical SDSS data and provide a phenomenological characterization of mid-IR, extragalactic sources. WISE is most sensitive at 3.4micron(W1) and least at 22micron(W4). The W1 band probes massive early-type galaxies out to z\gtrsim1. This is more distant than SDSS identified early-type galaxies, consistent with the fact that 28% of 3.4micron sources have faint or no r-band counterparts (r>22.2). In contrast, 92-95% of 12 and 22micron sources have SDSS optical counterparts with r<22.2. WISE 3.4micron detects 89.8% of the entire SDSS QSO catalog at SNR(W1)>7, but only 18.9% at 22micron with SNR(W4)>5. We show that WISE colors alone are effective in isolating stars (or local early-type galaxies), star-forming galaxies and strong AGN/QSOs at z<3. We highlight three major applications of WISE colors: (1) Selection of strong AGN/QSOs at z<3 using W1-W2>0.8 and W2<15.2 criteria, producing a census of this population. The surface density of these strong AGN/QSO candidates is 67.5+-0.14/deg^2. (2) Selection of dust-obscured, type-2 AGN/QSO candidates. We show that WISE W1-W2>0.8, W2<15.2 combined with r-W2>6 (Vega) colors can be used to identify type-2 AGN candidates. The fraction of these type-2 AGN candidates is 1/3rd of all WISE color-selected AGNs. (3) Selection of ULIRGs at z\sim2 with extremely red colors, r-W4>14 or well-detected 22micron sources lacking detections in the 3.4 and 4.6micron bands. The surface density of z~2 r-W4>14 ULIRGs is 0.9+-0.07/deg^2 at SNR(W4)>5 (flux(W4)>=2.5mJy), which is consistent with that inferred from smaller area Spitzer surveys. Optical spectroscopy of a small number of these high-redshift ULIRGs confirms our selection, and reveals a possible trend that optically fainter or r-W4 redder candidates are at higher redshifts.

preprint2012arXiv

IDCS J1426.5+3508: Sunyaev-Zel'dovich Measurement of a Massive IR-selected Cluster at z=1.75

We report 31 GHz CARMA observations of IDCS J1426.5+3508, an infrared-selected galaxy cluster at z = 1.75. A Sunyaev-Zel'dovich decrement is detected towards this cluster, indicating a total mass of M200 = (4.3 +/- 1.1) x 10^{14} Msun in agreement with the approximate X-ray mass of ~5 x 10^{14} Msun. IDCS J1426.5+3508 is by far the most distant cluster yet detected via the Sunyaev-Zel'dovich effect, and the most massive z >= 1.4 galaxy cluster found to date. Despite the mere ~1% probability of finding it in the 8.82 deg^2 IRAC Distant Cluster Survey, IDCS J1426.5+3508 is not completely unexpected in LCDM once the area of large, existing surveys is considered. IDCS J1426.5+3508 is, however, among the rarest, most extreme clusters ever discovered, and indeed is an evolutionary precursor to the most massive known clusters at all redshifts. We discuss how imminent, highly sensitive Sunyaev-Zel'dovich experiments will complement infrared techniques for statistical studies of the formation of the most massive galaxy clusters in the z > 1.5 Universe, including potential precursors to IDCS J1426.5+3508.

preprint2012arXiv

Jet and Torus Orientations in High Redshift Radio Galaxies

We examine the relative orientation of radio jets and dusty tori surrounding the AGN in powerful radio galaxies at z > 1. The radio core dominance R = P(20 GHz) /P(500 MHz) serves as an orientation indicator, measuring the ratio between the anisotropic Doppler-beamed core extended core emission and the isotropic lobe emission. Assuming a fixed cylindrical geometry for the hot, dusty torus, we derive its inclination i by fitting optically-thick radiative transfer models to spectral energy distributions obtained with the Spitzer Space Telescope. We find a highly significant anti-correlation (p < 0.0001) between R and i in our sample of 35 type 2 AGN combined with a sample of 18 z = 1 3CR sources containing both type 1 and 2 AGN. This analysis provides observational evidence both for the Unified scheme of AGN and for the common assumption that radio jets are in general perpendicular to the plane of the torus. The use of inclinations derived from mid-infrared photometry breaks several degeneracies which have been problematic in earlier analyses. We illustrate this by deriving the core Lorentz factor Gamma from the R-i anti-correlation, finding Gamma > 1.3.

preprint2012arXiv

Origin of the 12um Emission Across Galaxy Populations from WISE and SDSS Surveys

We cross-matched Wide-field Infrared Survey Explorer (WISE) sources brighter than 1 mJy at 12um with the Sloan Digital Sky Survey (SDSS) galaxy spectroscopic catalog to produce a sample of ~10^5 galaxies at <z>=0.08, the largest of its kind. This sample is dominated (70%) by star-forming (SF) galaxies from the blue sequence, with total IR luminosities in the range ~10^8-10^12 L_sun. We identify which stellar populations are responsible for most of the 12um emission. We find that most (~80%) of the 12um emission in SF galaxies is produced by stellar populations younger than 0.6 Gyr. In contrast, the 12um emission in weak AGN (L[OIII]<10^7 L_sun) is produced by older stars, with ages of ~1-3 Gyr. We find that L_[12um] linearly correlates with stellar mass for SF galaxies. At fixed 12um luminosity, weak AGN deviate toward higher masses since they tend to be hosted by massive, early-type galaxies with older stellar populations. Star-forming galaxies and weak AGN follow different L_[12um]-SFR (star formation rate) relations, with weak AGN showing excess 12um emission at low SFR (~0.02-1 M_sun/yr). This is likely due to dust grains heated by older stars. While the specific star formation rate (SSFR) of SF galaxies is nearly constant, the SSFR of weak AGN decreases by ~3 orders of magnitude, reflecting the very different star formation efficiencies between SF galaxies and massive, early-type galaxies. Stronger type II AGN in our sample (L_[OIII]>10^7 L_sun), act as an extension of massive SF galaxies, connecting the SF and weak AGN sequences. This suggests a picture where galaxies form stars normally until an AGN (possibly after a starburst episode) starts to gradually quench the SF activity. We also find that 4.6-12um color is a useful first-order indicator of SF activity in a galaxy when no other data are available.

preprint2012arXiv

Rapid Coeval Black Hole and Host Galaxy Growth in MRC 1138-262: The Hungry Spider

We present a detailed study of the infrared spectral energy distribution of the high-redshift radio galaxy MRC 1138-26 at z = 2.156, also known as the Spiderweb Galaxy. By combining photometry from Spitzer, Herschel and LABOCA we fit the rest-frame 5-300 um emission using a two component, starburst and active galactic nucleus (AGN), model. The total infrared (8 - 1000 um) luminosity of this galaxy is (1.97+/-0.28)x10^13 Lsun with (1.17+/-0.27) and (0.79+/-0.09)x10^13 Lsun due to the AGN and starburst components respectively. The high derived AGN accretion rate of \sim20% Eddington, and the measured star formation rate (SFR) of 1390pm150 Msun/yr, suggest that this massive system is in a special phase of rapid central black hole and host galaxy growth, likely caused by a gas rich merger in a dense environment. The accretion rate is sufficient to power both the jets and the previously observed large outflow. The high SFR and strong outflow suggest this galaxy could potentially exhaust its fuel for stellar growth in a few tens of Myr, although the likely merger of the radio galaxy with nearby satellites suggest bursts of star formation may recur again on time scales of several hundreds of Myr. The age of the radio lobes implies the jet started after the current burst of star formation, and therefore we are possibly witnessing the transition from a merger-induced starburst phase to a radio-loud AGN phase. We also note tentative evidence for [CII]158um emission. This paper marks the first results from the Herschel Galaxy Evolution Project (Project HeRGE), a systematic study of the evolutionary state of 71 high redshift, 1 < z < 5.2, radio galaxies.

preprint2012arXiv

The Herschel view of the environment of the radio galaxy 4C+41.17 at z = 3.8

We present Herschel observations at 70, 160, 250, 350 and 500 micron of the environment of the radio galaxy 4C+41.17 at z = 3.792. About 65% of the extracted sources are securely identified with mid-IR sources observed with the Spitzer Space Telescope at 3.6, 4.5, 5.8, 8 and 24 micron. We derive simple photometric redshifts, also including existing 850 micron and 1200 micron data, using templates of AGN, starburst-dominated systems and evolved stellar populations. We find that most of the Herschel sources are foreground to the radio galaxy and therefore do not belong to a structure associated with 4C+41.17. We do, however, find that the SED of the closest (~ 25" offset) source to the radio galaxy is fully consistent with being at the same redshift as 4C+41.17. We show that finding such a bright source that close to the radio galaxy at the same redshift is a very unlikely event, making the environment of 4C+41.17 a special case. We demonstrate that multi-wavelength data, in particular on the Rayleigh-Jeans side of the spectral energy distribution, allow us to confirm or rule out the presence of protocluster candidates that were previously selected by single wavelength data sets.

preprint2012arXiv

The Intense Starburst HDF850.1 in a Galaxy Overdensity at z=5.2 in the Hubble Deep Field

The Hubble Deep Field (HDF) is a region in the sky that provides one of the deepest multi-wavelength views of the distant universe and has led to the detection of thousands of galaxies seen throughout cosmic time. An early map of the HDF at a wavelength of 850 microns that is sensitive to dust emission powered by star formation revealed the brightest source in the field, dubbed HDF850.1. For more than a decade, this source remained elusive and, despite significant efforts, no counterpart at shorter wavelengths, and thus no redshift, size or mass, could be identified. Here we report, using a millimeter wave molecular line scan, an unambiguous redshift determination for HDF850.1 of z=5.183. This places HDF850.1 in a galaxy overdensity at z~5.2 in the HDF, corresponding to a cosmic age of only 1.1 Gyr after the Big Bang. This redshift is significantly higher than earlier estimates and higher than most of the >100 sub-millimeter bright galaxies identified to date. The source has a star formation rate of 850 M_sun/yr and is spatially resolved on scales of 5 kpc, with an implied dynamical mass of ~1.3x10^11 M_sun, a significant fraction of which is present in the form of molecular gas. Despite our accurate redshift and position, a counterpart arising from starlight remains elusive.

preprint2012arXiv

Using the Bright Ultra-Hard XMM-Newton Survey to define an IR selection of luminous AGN based on WISE colours

We present a highly complete and reliable mid-infrared (MIR) colour selection of luminous AGN candidates using the 3.4, 4.6, and 12 um bands of the WISE survey. The MIR colour wedge was defined using the wide-angle Bright Ultra-Hard XMM-Newton Survey (BUXS), one of the largest complete flux-limited samples of bright (f(4.5-10 keV)>6x10^{-14} erg cm^-2 s^-1) "ultra-hard" (4.5-10 keV) X-ray selected AGN to date. BUXS includes 258 objects detected over a total sky area of 44.43 deg^2 of which 251 are spectroscopically identified and classified, with 145 being type-1 AGN and 106 type-2 AGN. Our technique is designed to select objects with red MIR power-law spectral energy distributions (SED) in the three shortest bands of WISE and properly accounts for the errors in the photometry and deviations of the MIR SEDs from a pure power-law. The completeness of the MIR selection is a strong function of luminosity. At L(2-10 keV)>10^{44} erg s^-1, where the AGN is expected to dominate the MIR emission, 97.1_{-4.8}^{+2.2}% and 76.5_{-18.4}^{+13.3}% of the BUXS type-1 and type-2 AGN meet the selection. Our technique shows one of the highest reliability and efficiency of detection of the X-ray selected luminous AGN population with WISE amongst those in the literature. In the area covered by the BUXS survey our selection identifies 2755 AGN candidates detected with SNR>5 in the three shorter wavelength bands of WISE with 38.5% having a detection at 2-10 keV X-ray energies. We also analyzed the possibility of including the 22um WISE band to select AGN candidates, but neither the completeness nor the reliability of the selection improves. This is likely due to both the significantly shallower depth at 22um compared with the first three bands of WISE and star-formation contributing to the 22um emission at the WISE 22um sensitivity.

preprint2012arXiv

Wide-Field InfraRed Survey Telescope (WFIRST) Final Report

In December 2010, NASA created a Science Definition Team (SDT) for WFIRST, the Wide Field Infra-Red Survey Telescope, recommended by the Astro 2010 Decadal Survey as the highest priority for a large space mission. The SDT was chartered to work with the WFIRST Project Office at GSFC and the Program Office at JPL to produce a Design Reference Mission (DRM) for WFIRST. Part of the original charge was to produce an interim design reference mission by mid-2011. That document was delivered to NASA and widely circulated within the astronomical community. In late 2011 the Astrophysics Division augmented its original charge, asking for two design reference missions. The first of these, DRM1, was to be a finalized version of the interim DRM, reducing overall mission costs where possible. The second of these, DRM2, was to identify and eliminate capabilities that overlapped with those of NASA's James Webb Space Telescope (henceforth JWST), ESA's Euclid mission, and the NSF's ground-based Large Synoptic Survey Telescope (henceforth LSST), and again to reduce overall mission cost, while staying faithful to NWNH. This report presents both DRM1 and DRM2.

preprint2011arXiv

AGES: The AGN and Galaxy Evolution Survey

The AGN and Galaxy Evolution Survey (AGES) is a redshift survey covering, in its standard fields, 7.7 square degrees of the Bootes field of the NOAO Deep Wide-Field Survey (NDWFS). The final sample consists of 23745 redshifts. There are well-defined galaxy samples in ten bands (the Bw, R, I, J, K, IRAC 3.6, 4.5, 5.8 and 8.0 micron and MIPS 24 micron bands) to a limiting magnitude of I<20 mag for spectroscopy. For these galaxies, we obtained 18163 redshifts from a sample of 35200 galaxies, where random sparse sampling was used to define statistically complete sub-samples in all ten photometric bands. The median galaxy redshift is 0.31, and 90% of the redshifts are in the range 0.085<z<0.66. AGN were selected as radio, X-ray, IRAC mid-IR and MIPS 24 micron sources to fainter limiting magnitudes (I<22.5 mag for point sources). Redshifts were obtained for 4764 quasars and galaxies with AGN signatures, with 2926, 1718, 605, 119 and 13 above redshifts of 0.5, 1, 2, 3 and 4, respectively. We detail all the AGES selection procedures and present the complete spectroscopic redshift catalogs, spectra, and spectral energy distribution decompositions. The photometric redshift estimates are for all sources in the AGES samples.

preprint2011arXiv

Black Hole Mass Estimates Based on CIV are Consistent with Those Based on the Balmer Lines

Using a sample of high-redshift lensed quasars from the CASTLES project with observed-frame ultraviolet or optical and near-infrared spectra, we have searched for possible biases between supermassive black hole (BH) mass estimates based on the CIV, Halpha and Hbeta broad emission lines. Our sample is based upon that of Greene, Peng & Ludwig, expanded with new near-IR spectroscopic observations, consistently analyzed high S/N optical spectra, and consistent continuum luminosity estimates at 5100A. We find that BH mass estimates based on the FWHM of CIV show a systematic offset with respect to those obtained from the line dispersion, sigma_l, of the same emission line, but not with those obtained from the FWHM of Halpha and Hbeta. The magnitude of the offset depends on the treatment of the HeII and FeII emission blended with CIV, but there is little scatter for any fixed measurement prescription. While we otherwise find no systematic offsets between CIV and Balmer line mass estimates, we do find that the residuals between them are strongly correlated with the ratio of the UV and optical continuum luminosities. Removing this dependency reduces the scatter between the UV- and optical-based BH mass estimates by a factor of approximately 2, from roughly 0.35 to 0.18 dex. The dispersion is smallest when comparing the CIV sigma_l mass estimate, after removing the offset from the FWHM estimates, and either Balmer line mass estimate. The correlation with the continuum slope is likely due to a combination of reddening, host contamination and object-dependent SED shapes. When we add additional heterogeneous measurements from the literature, the results are unchanged.

preprint2011arXiv

The Hubble Space Telescope Cluster Supernova Survey: V. Improving the Dark Energy Constraints Above z>1 and Building an Early-Type-Hosted Supernova Sample

We present ACS, NICMOS, and Keck AO-assisted photometry of 20 Type Ia supernovae SNe Ia from the HST Cluster Supernova Survey. The SNe Ia were discovered over the redshift interval 0.623 < z < 1.415. Fourteen of these SNe Ia pass our strict selection cuts and are used in combination with the world's sample of SNe Ia to derive the best current constraints on dark energy. Ten of our new SNe Ia are beyond redshift $z=1$, thereby nearly doubling the statistical weight of HST-discovered SNe Ia beyond this redshift. Our detailed analysis corrects for the recently identified correlation between SN Ia luminosity and host galaxy mass and corrects the NICMOS zeropoint at the count rates appropriate for very distant SNe Ia. Adding these supernovae improves the best combined constraint on the dark energy density ρ_{DE}(z) at redshifts 1.0 < z < 1.6 by 18% (including systematic errors). For a LambdaCDM universe, we find Ω_Λ= 0.724 +0.015/-0.016 (68% CL including systematic errors). For a flat wCDM model, we measure a constant dark energy equation-of-state parameter w = -0.985 +0.071/-0.077 (68% CL). Curvature is constrained to ~0.7% in the owCDM model and to ~2% in a model in which dark energy is allowed to vary with parameters w_0 and w_a. Tightening further the constraints on the time evolution of dark energy will require several improvements, including high-quality multi-passband photometry of a sample of several dozen z>1 SNe Ia. We describe how such a sample could be efficiently obtained by targeting cluster fields with WFC3 on HST.

preprint2011arXiv

The Lick AGN Monitoring Project 2011: Reverberation Mapping of Markarian 50

The Lick AGN Monitoring Project 2011 observing campaign was carried out over the course of 11 weeks in Spring 2011. Here we present the first results from this program, a measurement of the broad-line reverberation lag in the Seyfert 1 galaxy Mrk 50. Combining our data with supplemental observations obtained prior to the start of the main observing campaign, our dataset covers a total duration of 4.5 months. During this time, Mrk 50 was highly variable, exhibiting a maximum variability amplitude of a factor of 4 in the U-band continuum and a factor of 2 in the H-beta line. Using standard cross-correlation techniques, we find that H-beta and H-gamma lag the V-band continuum by tau_cen = 10.64(-0.93,+0.82) and 8.43(-1.28,+1.30) days, respectively, while the lag of He II 4686 is unresolved. The H-beta line exhibits a symmetric velocity-resolved reverberation signature with shorter lags in the high-velocity wings than in the line core, consistent with an origin in a broad-line region dominated by orbital motion rather than infall or outflow. Assuming a virial normalization factor of f=5.25, the virial estimate of the black hole mass is (3.2+-0.5)*10^7 solar masses. These observations demonstrate that Mrk 50 is among the most promising nearby active galaxies for detailed investigations of broad-line region structure and dynamics.

preprint2011arXiv

The Mid-IR and X-ray Selected QSO Luminosity Function

We present the J-band luminosity function of 1838 mid-infrared and X-ray selected AGNs in the redshift range 0<z<5.85. These luminosity functions are constructed by combining the deep multi-wavelength broad-band observations from the UV to the mid-IR of the NDWFS Bootes field with the X-ray observations of the XBootes survey and the spectroscopic observations of the same field by AGES. Our sample is primarily composed of IRAC-selected AGNs, targeted using modifications of the Stern et al.(2005) criteria, complemented by MIPS 24 microns and X-ray selected AGNs to alleviate the biases of IRAC mid-IR selection against z~4.5 quasars and AGNs faint with respect to their hosts. This sample provides an accurate link between low and high redshift AGN luminosity functions and does not suffer from the usual incompleteness of optical samples at z~3. We find that the space density of the brightest quasars strongly decreases from z=3 to z=0, while the space density of faint quasars is at least flat, and possibly increasing, over the same redshift range. At z>3 we observe a decrease in the space density of quasars of all brightnesses. We model the luminosity function by a double power-law and find that its evolution cannot be described by either pure luminosity or pure density evolution, but must be a combination of both. Our best-fit model has bright and faint power-law indices consistent with the low redshift measurements based on the 2QZ and 2SLAQ surveys and it generally agrees with the number of bright quasars predicted by other LFs at all redshifts. If we construct the QSO luminosity function using only the IRAC-selected AGNs, we find that the biases inherent to this selection method significantly modify the behavior of phi*(z) only for z<1 and have no significant impact upon the characteristic magnitude M*_J(z).

preprint2011arXiv

Three QSOs acting as strong gravitational lenses

We report the discovery of three new cases of QSOs acting as strong gravitational lenses on background emission line galaxies: SDSS J0827+5224 (zQSO = 0.293, zs = 0.412), SDSS J0919+2720 (zQSO = 0.209, zs = 0.558), SDSS J1005+4016 (zQSO = 0.230, zs = 0.441). The selection was carried out using a sample of 22,298 SDSS spectra displaying at least four emission lines at a redshift beyond that of the foreground QSO. The lensing nature is confirmed from Keck imaging and spectroscopy, as well as from HST/WFC3 imaging in the F475W and F814W filters. Two of the QSOs have face-on spiral host galaxies and the third is a QSO+galaxy pair. The velocity dispersion of the host galaxies, inferred from simple lens modeling, is between σ_v = 210 and 285 km/s, making these host galaxies comparable in mass with the SLACS sample of early-type strong lenses.

preprint2010arXiv

Imaging the molecular gas in a submm galaxy at z = 4.05: cold mode accretion or a major merger?

We present a high resolution (down to 0.18"), multi-transition imaging study of the molecular gas in the z = 4.05 submillimeter galaxy GN20. GN20 is one of the most luminous starburst galaxy known at z > 4, and is a member of a rich proto-cluster of galaxies at z = 4.05 in GOODS-North. We have observed the CO 1-0 and 2-1 emission with the VLA, the CO 6-5 emission with the PdBI Interferometer, and the 5-4 emission with CARMA. The H_2 mass derived from the CO 1-0 emission is 1.3 \times 10^{11} (α/0.8) Mo. High resolution imaging of CO 2-1 shows emission distributed over a large area, appearing as partial ring, or disk, of ~ 10kpc diameter. The integrated CO excitation is higher than found in the inner disk of the Milky Way, but lower than that seen in high redshift quasar host galaxies and low redshift starburst nuclei. The VLA CO 2-1 image at 0.2" resolution shows resolved, clumpy structure, with a few brighter clumps with intrinsic sizes ~ 2 kpc. The velocity field determined from the CO 6-5 emission is consistent with a rotating disk with a rotation velocity of ~ 570 km s^{-1} (using an inclination angle of 45^o), from which we derive a dynamical mass of 3 \times 10^{11} \msun within about 4 kpc radius. The star formation distribution, as derived from imaging of the radio synchrotron and dust continuum, is on a similar scale as the molecular gas distribution. The molecular gas and star formation are offset by ~ 1" from the HST I-band emission, implying that the regions of most intense star formation are highly dust-obscured on a scale of ~ 10 kpc. The large spatial extent and ordered rotation of this object suggests that this is not a major merger, but rather a clumpy disk accreting gas rapidly in minor mergers or smoothly from the proto-intracluster medium. ABSTRACT TRUNCATED

preprint2010arXiv

Low Resolution Spectral Templates For AGNs and Galaxies From 0.03 -- 30 microns

We present a set of low resolution empirical SED templates for AGNs and galaxies in the wavelength range from 0.03 to 30 microns based on the multi-wavelength photometric observations of the NOAO Deep-Wide Field Survey Bootes field and the spectroscopic observations of the AGN and Galaxy Evolution Survey. Our training sample is comprised of 14448 galaxies in the redshift range 0<~z<~1 and 5347 likely AGNs in the range 0<~z<~5.58. We use our templates to determine photometric redshifts for galaxies and AGNs. While they are relatively accurate for galaxies, their accuracies for AGNs are a strong function of the luminosity ratio between the AGN and galaxy components. Somewhat surprisingly, the relative luminosities of the AGN and its host are well determined even when the photometric redshift is significantly in error. We also use our templates to study the mid-IR AGN selection criteria developed by Stern et al.(2005) and Lacy et al.(2004). We find that the Stern et al.(2005) criteria suffers from significant incompleteness when there is a strong host galaxy component and at z =~ 4.5, when the broad Halpha emission line is redshifted into the [3.6] band, but that it is little contaminated by low and intermediate redshift galaxies. The Lacy et al.(2004) criterion is not affected by incompleteness at z =~ 4.5 and is somewhat less affected by strong galaxy host components, but is heavily contaminated by low redshift star forming galaxies. Finally, we use our templates to predict the color-color distribution of sources in the upcoming WISE mission and define a color criterion to select AGNs analogous to those developed for IRAC photometry. We estimate that in between 640,000 and 1,700,000 AGNs will be identified by these criteria, but will have serious completeness problems for z >~ 3.4.

preprint2010arXiv

SDWFS-MT-1: A Self-Obscured Luminous Supernova at z~0.2

We report the discovery of a six-month-long mid-infrared transient, SDWFS-MT-1 (aka SN 2007va), in the Spitzer Deep, Wide-Field Survey of the NOAO Deep Wide-Field Survey Bootes field. The transient, located in a z=0.19 low luminosity (M_[4.5]~-18.6 mag, L/L_MilkyWay~0.01) metal-poor (12+log(O/H)~7.8) irregular galaxy, peaked at a mid-infrared absolute magnitude of M_[4.5]~-24.2 in the 4.5 micron Spitzer/IRAC band and emitted a total energy of at least 10^51 ergs. The optical emission was likely fainter than the mid-infrared, although our constraints on the optical emission are poor because the transient peaked when the source was "behind" the Sun. The Spitzer data are consistent with emission by a modified black body with a temperature of ~1350 K. We rule out a number of scenarios for the origin of the transient such as a Galactic star, AGN activity, GRB, tidal disruption of a star by a black hole and gravitational lensing. The most plausible scenario is a supernova exploding inside a massive, optically thick circumstellar medium, composed of multiple shells of previously ejected material. If the proposed scenario is correct, then a significant fraction (~10%) of the most luminous supernova may be self-enshrouded by dust not only before but also after the supernova occurs. The spectral energy distribution of the progenitor of such a supernova would be a slightly cooler version of eta Carina, peaking at 20-30 microns.

preprint2010arXiv

The Great Observatories Origins Deep Survey: Constraints on the Lyman Continuum Escape Fraction Distribution of Lyman--Break Galaxies at 3.4<z<4.5

We use ultra-deep ultraviolet VLT/VIMOS intermediate-band and VLT/FORS1 narrow-band imaging in the GOODS Southern field to derive limits on the distribution of the escape fraction (f_esc) of ionizing radiation for L >~ L*(z=3) Lyman Break Galaxies (LBGs) at redshift 3.4--4.5. Only one LBG, at redshift z=3.795, is detected in its Lyman continuum (LyC; S/N~5.5), the highest redshift galaxy currently known with a direct detection. Its ultraviolet morphology is quite compact (R_eff=0.8, kpc physical). Three out of seven AGN are also detected in their LyC, including one at redshift z=3.951 and z850 = 26.1. From stacked data (LBGs) we set an upper limit to the average f_esc in the range 5%--20%, depending on the how the data are selected (e.g., by magnitude and/or redshift). We undertake extensive Monte Carlo simulations that take into account intergalactic attenuation, stellar population synthesis models, dust extinction and photometric noise in order to explore the moments of the distribution of the escaping radiation. Various distributions (exponential, log-normal and Gaussian) are explored. We find that the median f_esc is lower than ~6% with an 84% percentile limit not larger than 20%. If this result remains valid for fainter LBGs down to current observational limits, then the LBG population might be not sufficient to account for the entire photoionization budget at the redshifts considered here, with the exact details dependent upon the assumed ionizing background and QSO contribution thereto. It is possible that f_esc depends on the UV luminosity of the galaxies, with fainter galaxies having higher f_esc, and estimates of f_esc from a sample of faint LBG from the HUDF (i775<28.5) are in broad quantitative agreement with such a scenario.

preprint2010arXiv

Very High Gas Fractions and Extended Gas Reservoirs in z=1.5 Disk Galaxies

We present evidence for very high gas fractions and extended molecular gas reservoirs in normal, near-infrared selected (BzK) galaxies at z~1.5, based on multi-configuration CO[2-1] observations obtained at the IRAM PdBI. Six of the six galaxies observed were securely detected. High resolution observations resolve the CO emission in four of them, implying sizes of order of 6-11 kpc and suggesting the presence of rotation. The UV morphologies are consistent with clumpy, unstable disks, and the UV sizes are consistent with the CO sizes. The star formation efficiencies are homogeneously low and similar to local spirals - the resulting gas depletion times are ~0.5 Gyr, much higher than what is seen in high-z submm galaxies and quasars. The CO luminosities can be predicted to within 0.15 dex from the star formation rates and stellar masses, implying a tight correlation of the gas mass with these quantities. We use dynamical models of clumpy disk galaxies to derive dynamical masses. These models are able to reproduce the peculiar spectral line shapes of the CO emission. After accounting for the stellar and dark matter masses we derive gas masses of 0.4-1.2x10^11 Msun. The conversion factor is very high: alpha_CO=3.6+-0.8, consistent with the Galaxy but four times higher than that of local ultra-luminous IR galaxies. The gas accounts for an impressive 50-65% of the baryons within the galaxies' half light radii. We are witnessing truly gas-dominated galaxies at z~1.5, a finding that explains the high specific SFRs observed for z>1 galaxies. The BzK galaxies can be viewed as scaled-up versions of local disk galaxies, with low efficiency star formation taking place inside extended, low excitation gas disks. They are markedly different than local ULIRGs and high-z submm galaxies, which have more excited and compact gas.

preprint2010arXiv

X-Ray Emission from Two Infrared-Selected Galaxy Clusters at z>1.4 in the IRAC Shallow Cluster Survey

We report the X-ray detection of two z>1.4 infrared-selected galaxy clusters from the IRAC Shallow Cluster Survey (ISCS). We present new data from the Hubble Space Telescope and the W. M. Keck Observatory that spectroscopically confirm cluster ISCS J1432.4+3250 at z=1.49, the most distant of 18 confirmed z>1 clusters in the ISCS to date. We also present new spectroscopy for ISCS J1438.1+3414, previously reported at z = 1.41, and measure its dynamical mass. Clusters ISCS J1432.4+3250 and ISCS J1438.1+3414 are detected in 36ks and 143ks Chandra exposures at significances of 5.2 sigma and 9.7 sigma, from which we measure total masses of log(M_{200,Lx}/Msun) = 14.4 +/- 0.2 and 14.35^{+0.14}_{-0.11}, respectively. The consistency of the X-ray and dynamical properties of these high redshift clusters further demonstrates that the ISCS is robustly detecting massive clusters to at least z = 1.5.

preprint2009arXiv

An Intensive HST Survey for z>1 Supernovae by Targeting Galaxy Clusters

We present a new survey strategy to discover and study high redshift Type Ia supernovae (SNe Ia) using the Hubble Space Telescope (HST). By targeting massive galaxy clusters at 0.9<z<1.5, we obtain a twofold improvement in the efficiency of finding SNe compared to an HST field survey and a factor of three improvement in the total yield of SN detections in relatively dust-free red-sequence galaxies. In total, sixteen SNe were discovered at z>0.95, nine of which were in galaxy clusters. This strategy provides a SN sample that can be used to decouple the effects of host galaxy extinction and intrinsic color in high redshift SNe, thereby reducing one of the largest systematic uncertainties in SN cosmology.

preprint2009arXiv

Mid-Infrared Galaxy Luminosity Functions from the AGN and Galaxy Evolution Survey

We present galaxy luminosity functions at 3.6, 4.5, 5.8, and 8.0 micron measured by combining photometry from the IRAC Shallow Survey with redshifts from the AGN and Galaxy Evolution Survey of the NOAO Deep Wide-Field Survey Bootes field. The well-defined IRAC samples contain 3800-5800 galaxies for the 3.6-8.0 micron bands with spectroscopic redshifts and z < 0.6. We obtained relatively complete luminosity functions in the local redshift bin of z < 0.2 for all four IRAC channels that are well fit by Schechter functions. We found significant evolution in the luminosity functions for all four IRAC channels that can be fit as an evolution in M* with redshift, ΔM* = Qz. While we measured Q=1.2\pm0.4 and 1.1\pm0.4 in the 3.6 and 4.5 micron bands consistent with the predictions from a passively evolving population, we obtained Q=1.8\pm1.1 in the 8.0 micron band consistent with other evolving star formation rate estimates. We compared our LFs with the predictions of semi-analytical galaxy formation and found the best agreement at 3.6 and 4.5 micron, rough agreement at 8.0 micron, and a large mismatch at 5.8 micron. These models also predicted a comparable Q value to our luminosity functions at 8.0 micron, but predicted smaller values at 3.6 and 4.5 micron. We also measured the luminosity functions separately for early and late-type galaxies. While the luminosity functions of late-type galaxies resemble those for the total population, the luminosity functions of early-type galaxies in the 3.6 and 4.5 micron bands indicate deviations from the passive evolution model, especially from the measured flat luminosity density evolution. Combining our estimates with other measurements in the literature, we found (53\pm18)% of the present stellar mass of early-type galaxies has been assembled at z=0.7.

preprint2009arXiv

Spectroscopic Observations of Lyman-Break Galaxies at Redshift ~ 4, 5 and 6 in the GOODS-South Field

We report on observations of Lyman-break galaxies (LBGs) selected from the Great Observatories Origins Deep Survey at mean redshift z~4, 5 and 6, obtained with FORS2 spectrograph at the ESO-VLT.This program has yielded spectroscopic identifications for 114 galaxies of which 51 are at z~4, 31 at z~5, and 32 at z~6. The adopted selection criteria are effective, identifying galaxies at the expected redshift with minimal foreground contamination. Once selection effects are properly accounted for, the rest-frame UV spectra of the higher-redshift LBGs appear to be similar to their counterparts at z~3. As at z~3, LBGs at z~4 and z~5 show Lya both in emission and in absorption; when in absorption, strong interstellar lines and red UV spectra are observed, a fact also measured at z~2 and 3. At z~6, sensitivity issues bias our sample towards galaxies with Lya in emission; nevertheless, these spectra appear to be similar to their lower-redshift counterparts. As in other studies at similar redshifts, we find clear evidence that brighter LBGs tend to have weaker Lya emission lines. At fixed rest-frame UV luminosity, the equivalent width of the Lya emission line is larger at higher redshifts. At all redshifts where the measurements can be reliably made, the redshift of the Lya emission line turns out to be larger than that of the interstellar absorption lines,with a median velocity difference DeltaV~400 km/s at z~4 and 5, consistent with results at lower redshifts. This shows that powerful, large-scale winds are common at high redshift.In general,there is no strong correlation between the morphology of the UV light and the spectroscopic properties. However, galaxies with deep interstellar absorption lines and strong Lya absorption appear to be more diffuse than galaxies with Lya in emission.(abridged)

preprint2008arXiv

The Environments of High Redshift QSOs

We present a sample of $i_{775}$-dropout candidates identified in five Hubble Advanced Camera for Surveys fields centered on Sloan Digital Sky Survey QSOs at redshift $z\sim 6$. Our fields are as deep as the Great Observatory Origins Deep Survey (GOODS) ACS images which are used as a reference field sample. We find them to be overdense in two fields, underdense in two fields, and as dense as the average density of GOODS in one field. The two excess fields show significantly different color distributions from that of GOODS at the 99% confidence level, strengthening the idea that the excess objects are indeed associated with the QSO. The distribution of $i_{775}$-dropout counts in the five fields is broader than that derived from GOODS at the 80% to 96% confidence level, depending on which selection criteria were adopted to identify $i_{775}$-dropouts; its width cannot be explained by cosmic variance alone. Thus, QSOs seem to affect their environments in complex ways. We suggest the picture where the highest redshift QSOs are located in very massive overdensities and are therefore surrounded by an overdensity of lower mass halos. Radiative feedback by the QSO can in some cases prevent halos from becoming galaxies, thereby generating in extreme cases an underdensity of galaxies. The presence of both enhancement and suppression is compatible with the expected differences between lines of sight at the end of reionization as the presence of residual diffuse neutral hydrogen would provide young galaxies with shielding from the radiative effects of the QSO.

preprint2008arXiv

Two bright submillimeter galaxies in a z=4.05 proto-cluster in GOODS-North, and accurate radio-infrared photometric redshifts

We present the serendipitous discovery of z=4.05 molecular gas CO emission lines with the IRAM Plateau de Bure interferometer coincident with GN20 and GN20.2, two luminous submillimeter galaxies (SMGs) in the Great Observatories Origins Deep Survey North field (GOODS-N). These are among the most distant submillimeter-selected galaxies reliably identified through CO emission and also some of the most luminous known. In terms of CO to bolometric luminosities, stellar mass and star formation rates (SFRs), these newly discovered z>4 SMGs are similar to z~1.5-3 SMGs studied to date. These z~4 SMGs have much higher specific SFRs than typical B-band dropout Lyman break galaxies at the same redshift. The stellar mass-SFR correlation for normal galaxies does not seem to evolve much further, between z~2 and z~4. A significant z=4.05 spectroscopic redshift spike is observed in GOODS-N, and a strong spatial overdensity of B-band dropouts and IRAC selected z>3.5 galaxies appears to be centered on the GN20 and GN20.2 galaxies. This suggests a proto-cluster structure with total mass ~10^14 Msun. Using photometry at mid-IR, submm and radio wavelengths, we show that reliable photometric redshifts (Dz/(1+z)~0.1) can be derived for SMGs over 1<z<4. This new photometric redshift technique has been used to provide a first estimate of the space density of 3.5<z<6 hyper-luminous starburst galaxies, and to show that they contribute substantially to the SFR density at early epochs. Many of these high-redshift starbursts will be within reach of Herschel. We find that the radio to mid-IR flux density ratio can be used to select z>3.5 starbursts, regardless of their submm/mm emission [abridged].

preprint2007arXiv

Galaxy Cluster Correlation Function to z ~ 1.5 in the IRAC Shallow Cluster Survey

We present the galaxy cluster autocorrelation function of 277 galaxy cluster candidates with 0.25 \le z \le 1.5 in a 7 deg^2 area of the IRAC Shallow Cluster Survey. We find strong clustering throughout our galaxy cluster sample, as expected for these massive structures. Specifically, at <z> = 0.5 we find a correlation length of r_0 = 17.40^{+3.98}_{-3.10} h^-1 Mpc, in excellent agreement with the Las Campanas Distant Cluster Survey, the only other non-local measurement. At higher redshift, <z> = 1, we find that strong clustering persists, with a correlation length of r_0=19.14^{+5.65}_{-4.56} h^-1 Mpc. A comparison with high resolution cosmological simulations indicates these are clusters with halo masses of \sim 10^{14} Msun, a result supported by estimates of dynamical mass for a subset of the sample. In a stable clustering picture, these clusters will evolve into massive (10^{15} Msun) clusters by the present day.

preprint2007arXiv

The Great Observatories Origins Deep Survey. VLT/FORS2 Spectroscopy in the GOODS-South Field: Part III

Aims. We present the full data set of the spectroscopic campaign of the ESO/GOODS program in the GOODS-South field, obtained with the FORS2 spectrograph at the ESO/VLT. Method. Objects were selected as candidates for VLT/FORS2 observations primarily based on the expectation that the detection and measurement of their spectral features would benefit from the high throughput and spectral resolution of FORS2. The reliability of the redshift estimates is assessed using the redshift-magnitude and color-redshift diagrams, and comparing the results with public data. Results. Including the third part of the spectroscopic campaign (12 masks) to the previous work (26 masks, Vanzella et al. 2005, 2006), 1715 spectra of 1225 individual targets have been analyzed. The actual spectroscopic catalog provides 887 redshift determinations. The typical redshift uncertainty is estimated to be sigma(z) ~ 0.001. Galaxies have been selected adopting different color criteria and using photometric redshifts. The resulting redshift distribution typically spans two domains: from z=0.5 to 2 and z=3.5 to 6.3. The reduced spectra and the derived redshifts are released to the community through the ESO web page http://www.eso.org/science/goods/

preprint2007arXiv

The reversal of the star formation-density relation in the distant universe

We study the relationship between the local environment of galaxies and their star formation rate (SFR) in the Great Observatories Origins Deep Survey, GOODS, at z~1, from ultradeep imaging at 24 microns with the MIPS camera onboard Spitzer. We show that the star formation-density relation observed locally was reversed at z~1: the average SFR of an individual galaxy increased with local galaxy density when the universe was less than half its present age. Hierarchical galaxy formation models (simulated lightcones from the Millennium model) predicted such a reversal to occur only at only at earlier epochs (z>2) and at a lower level. We present a remarkable structure at z~1.016, containing X-ray traced galaxy concentrations, which will eventually merge into a Virgo-like cluster. This structure illustrates how the individual SFR of galaxies increases with density at the ~1-2 Mpc scale. The SFR of z~1 galaxies is found to correlate with stellar mass suggesting that mass plays a role in the observed star formation-density trend. However the specific SFR (=SFR/M*) decreases with stellar mass while it increases with galaxy density, which implies that the environment does directly affect the star formation activity of galaxies. Major mergers do not appear to be the unique or even major cause for this effect since nearly half (46%) of the luminous infrared galaxies (LIRGs) at z~1 present the HST-ACS morphology of spirals, while only a third present a clear signature of major mergers. The remaining galaxies are divided into compact (9%) and irregular (14%) galaxies. Moreover, the specific SFR of major mergers is only marginally stronger than that of spirals. Reproducing the SFR-density relation at z ~ 1 is a new challenge for models, requiring a correct balance between mergers and in-situ star formation at early epochs.

preprint2006arXiv

Mid-Infrared Selection of Brown Dwarfs and High-Redshift Quasars

We discuss color selection of rare objects in a wide-field, multiband survey spanning from the optical to the mid-infrared. Simple color criteria simultaneously identify and distinguish two of the most sought after astrophysical sources: the coolest brown dwarfs and the most distant quasars. We present spectroscopically-confirmed examples of each class identified in the IRAC Shallow Survey of the Bootes field of the NOAO Deep Wide-Field Survey. ISS J142950.9+333012 is a T4.5 brown dwarf at a distance of approximately 42 pc, and ISS J142738.5+331242 is a radio-loud quasar at redshift z=6.12. Our selection criteria identify a total of four candidates over 8 square degrees of the Bootes field. The other two candidates are both confirmed 5.5<z<6 quasars, previously reported by Cool et al. (2006). We discuss the implications of these discoveries and conclude that there are excellent prospects for extending such searches to cooler brown dwarfs and higher redshift quasars.

preprint2005arXiv

Direct evidence for an early reionization of the Universe?

We examine the possible reionization of the intergalactic medium (IGM) by the source UDF033238.7-274839.8 (hereafter HUDF-JD2), which was discovered in deep {\it HST}/VLT/{\it Spitzer} images obtained as part of the Great Observatory Origins Deep Survey and {\it Hubble} Ultra-Deep Field projects. Mobasher et al (2005) have identified HUDF-JD2 as a massive ($\sim6\times10^{11}M_\odot$) post-starburst galaxy at redshift z$\gtrsim6.5$. We find that HUDF-JD2 may be capable of reionizing its surrounding region of the Universe, starting the process at a redshift as high as z$\approx 15 \pm5$.

preprint2005arXiv

Evidence for a Massive Post-Starburst Galaxy at z ~ 6.5

We present results from a search for high-redshift J--band ``dropout'' galaxies in the portion of the GOODS southern field that is covered by extremely deep imaging from the Hubble Ultradeep Field (HUDF).Using observations at optical, near-infrared and mid-infrared wavelengths from the Hubble and Spitzer Space Telescopes and the ESO-VLT, we search for very massive galaxies at high redshifts and find one particularly remarkable candidate. Its spectral energy distribution is consistent with a galaxy at z ~ 6.5 and a stellar mass of 6x10e11 M(sun) (for a Salpeter IMF). We interpret a prominent photometric break between the near-infrared and Spitzer bandpasses as the 3646A Balmer discontinuity. The best-fitting models have low reddening and ages of several hundred Myr, placing the formation of the bulk of the stars at z > 9. Alternative models of dusty galaxies at z ~ 2.5 are possible but provide significantly poorer fits. The object is detected with Spitzer at 24 micron. This emission originats from an obscured active nucleus or star formation. We present optical and near-infrared spectroscopy which has, thus far, failed to detect any spectral features. This helps limit the solution in which the galaxy is a starburst or active galaxy at z ~ 2.5. If the high-redshift interpretation is correct, this object would be an example of a galaxy that formed by a process strongly resembling traditional models of monolithic collapse, in a way which a very large mass of stars formed within a remarkably short period of time, at very high redshift.

preprint2004arXiv

Number Counts At 3 < lambda < 10 um From the Spitzer Space Telescope

Infrared source counts at wavelengths 3 < lambda < 10 um cover more than 10 magnitudes in source brightness, four orders of magnitude in surface density, and reach an integrated surface density of 10^5 sources/deg^2. At m<14 mag, most of the sources are Galactic stars, in agreement with models. After removal of Galactic stars, galaxy counts are consistent with what few measurements exist at nearby wavelengths. At 3.6 and 4.5 um, the galaxy counts follow the expectations of a Euclidean world model down to ~16 mag and drop below the Euclidean curve for fainter magnitudes. Counts at these wavelengths begin to show decreasing completeness around magnitude 19.5. At 5.8 and 8 um, the counts relative to a Euclidean world model show a large excess at bright magnitudes. This is probably because local galaxies emit strongly in the aromatic dust (``PAH'') features. The counts at 3.6 um resolve <50% of the Cosmic Infrared Background at that wavelength.