• ### An efficient differential quadrature method for fractional advection-diffusion equation(1611.01597)

Oct. 18, 2017 math.NA
This article studies a direct numerical approach for fractional advection-diffusion equations (ADEs). Using a set of cubic trigonometric B-splines as test functions, a differential quadrature (DQ) method is firstly proposed for the 1D and 2D time-fractional ADEs of order $(0,1]$. The weighted coefficients are determined, and with them, the original equation is transformed into a group of general ordinary differential equations (ODEs), which are discretized by an effective difference scheme or Runge-Kutta method. The stability is investigated under a mild theoretical condition. Secondly, based on a set of cubic B-splines, we develop a new Crank-Nicolson type DQ method for the 2D space-fractional ADEs without advection. The DQ approximations to fractional derivatives are introduced and the values of the fractional derivatives of B-splines are computed by deriving explicit formulas. The presented DQ methods are evaluated on five benchmark problems and the concrete simulations of the unsteady propagation of solitons and Gaussian pulse. In comparison with the existing algorithms in the open literature, numerical results finally illustrate the validity and accuracy.
• ### A deep X-ray view of the bare AGN Ark120. IV. XMM-Newton and NuSTAR spectra dominated by two temperature (warm, hot) Comptonization processes(1707.08907)

July 27, 2017 astro-ph.HE
We perform an X-ray spectral analysis of the brightest and cleanest bare AGN known so far, Ark 120, in order to determine the process(es) at work in the vicinity of the SMBH. We present spectral analysis of data from an extensive campaign observing Ark 120 in X-rays with XMM-Newton (4$\times$120 ks, 2014 March 18-24), and NuSTAR (65.5 ks, 2014 March 22). During this very deep X-ray campaign, the source was caught in a high flux state similar to the earlier 2003 XMM-Newton observation, and about twice as bright as the lower-flux observation in 2013. The spectral analysis confirms the "softer when brighter" behaviour of Ark 120. The four XMM-Newton/pn spectra are characterized by the presence of a prominent soft X-ray excess and a significant FeK$\alpha$ complex. The continuum is very similar above about 3 keV, while significant variability is present for the soft X-ray excess. We find that relativistic reflection from a constant-density, flat accretion disk cannot simultaneously produce the soft excess, broad FeK$\alpha$ complex, and hard X-ray excess. Instead, Comptonization reproduces the broadband (0.3-79 keV) continuum well, together with a contribution from a mildly relativistic disk reflection spectrum. During this 2014 observational campaign, the soft X-ray spectrum of Ark 120 below $\sim$0.5 keV was found to be dominated by Comptonization of seed photons from the disk by a warm ($kT_{\rm e}$$\sim0.5 keV), optically-thick corona (\tau$$\sim$9). Above this energy, the X-ray spectrum becomes dominated by Comptonization from electrons in a hot optically thin corona, while the broad FeK$\alpha$ line and the mild Compton hump result from reflection off the disk at several tens of gravitational radii.
• ### A Long Look At MCG-5-23-16 With NuSTAR: I- Relativistic Reflection And Coronal Properties(1701.02309)

Jan. 9, 2017 astro-ph.HE
MCG-5-23-16 was targeted in early 2015 with a half mega-seconds observing campaign using NuSTAR. Here we present the spectral analysis of these datasets along with an earlier observation and study the relativistic reflection and the primary coronal source. The data show strong reflection features in the form of both narrow and broad iron lines plus a Compton reflection hump. A cutoff energy is significantly detected in all exposures. The shape of the reflection spectrum does not change in the two years spanned by the observations, suggesting a stable geometry. A strong positive correlation is found between the cutoff energy and both the hard X-ray flux and spectral index. The measurements imply that the coronal plasma is not at the runaway electron-positron pair limit, and instead contains mostly electrons. The observed variability in the coronal properties is driven by a variable optical depth. A constant heating to cooling ratio is measured implying that there is a feedback mechanism in which a significant fraction of the photons cooling the corona are due to reprocessed hard X-rays.
• ### NuSTAR observations of WISE J1036+0449, a Galaxy at z$\sim1$ obscured by hot dust(1609.04808)

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

Oct. 27, 2016 astro-ph.GA, astro-ph.HE
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$\alpha$ 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.
• ### Disk-Wind Connection During the Heartbeats of GRS 1915+105(1610.05772)

Oct. 18, 2016 astro-ph.HE
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 $\rho$ 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.
• ### The NuSTAR Extragalactic Surveys: The Number Counts of Active Galactic Nuclei and the Resolved Fraction of the Cosmic X-ray Background(1511.04183)

Sept. 10, 2016 astro-ph.HE
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.
• ### First NuSTAR observations of the BL Lac - type blazar PKS~2155-304: constraints on the jet content and distribution of radiating particles(1609.02203)

Sept. 7, 2016 astro-ph.HE
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.
• ### The Phoenix galaxy as seen by NuSTAR(1609.00374)

Sept. 1, 2016 astro-ph.HE
Aims. We study the long-term variability of the well-known Seyfert 2 galaxy Mrk 1210 (a.k.a. UGC 4203, or the Phoenix galaxy). Methods. The source was observed by many X-ray facilities in the last 20 years. Here we present a NuSTAR observation and put the results in context of previously published observations. Results. NuSTAR observed Mrk 1210 in 2012 for 15.4 ks. The source showed Compton-thin obscuration similar to that observed by Chandra, Suzaku, BeppoSAX and XMM-Newton over the past two decades, but different from the first observation by ASCA in 1995, in which the active nucleus was caught in a low flux state - or obscured by Compton-thick matter, with a reflection-dominated spectrum. Thanks to the high-quality hard X-ray spectrum obtained with NuSTAR and exploiting the long-term spectral coverage spanning 16.9 years, we can precisely disentangle the transmission and reflection components and put constraints on both the intrinsic long-term variability and hidden nucleus scenarios. In the former case, the distance between the reflector and the source must be at least ~ 2 pc, while in the latter one the eclipsing cloud may be identified with a water maser-emitting clump.
• ### Extremes of the jet-accretion power relation of blazars, as explored by NuSTAR(1510.08849)

Aug. 22, 2016 astro-ph.HE
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.
• ### Hard X-ray Emission from Sh2-104: A NuSTAR search for Gamma-ray Counterparts(1605.03641)

May 17, 2016 astro-ph.HE
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.
• ### A Hard X-ray Study of the Normal Star-Forming Galaxy M83 with NuSTAR(1604.07441)

April 25, 2016 astro-ph.GA, astro-ph.HE
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.
• ### NuSTAR and Swift observations of the very high state in GX 339-4: Weighing the black hole with X-rays(1603.03777)

March 11, 2016 astro-ph.HE
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${\deg}$\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).
• ### The first X-ray imaging spectroscopy of quiescent solar active regions with NuSTAR(1603.01069)

March 3, 2016 astro-ph.SR
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.
• ### NuSTAR, Swift, and GROND observations of the flaring MeV blazar: PMN J0641$-$0320(1602.06446)

Feb. 20, 2016 astro-ph.CO, astro-ph.HE
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.
• ### NuSTAR observations of water megamaser AGN(1602.03185)

Feb. 9, 2016 astro-ph.GA, astro-ph.HE
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.
• ### IC 751: a new changing-look AGN discovered by NuSTAR(1602.00702)

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

Sept. 3, 2015 astro-ph.HE
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.
• ### Dark matter line emission constraints from NuSTAR observations of the Bullet Cluster(1507.01378)

July 6, 2015 astro-ph.CO
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.
• ### NuSTAR Reveals Extreme Absorption in z < 0.5 Type 2 Quasars(1506.05120)

June 16, 2015 astro-ph.GA, astro-ph.HE
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.
• ### The complex accretion geometry of GX 339-4 as seen by NuSTAR and Swift(1506.01381)

June 3, 2015 astro-ph.HE
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.
• ### NuSTAR and Suzaku observations of the hard state in Cygnus X-1: locating the inner accretion disk(1506.00007)

May 29, 2015 astro-ph.HE
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{\sigma} 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.
• ### Rapid Variability of Blazar 3C 279 during Flaring States in 2013-2014 with Joint Fermi-LAT, NuSTAR, Swift, and Ground-Based Multi-wavelength Observations(1502.04699)

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

Feb. 25, 2015 astro-ph.HE
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}$.