Search

Grazing-incidence X-ray optics have revolutionized X-ray astrophysics. The ability to concentrate flux to a tiny detection region provides a dramatic reduction in background and a consequent very large improvement in sensitivity. The X-ray optics also permit use of small-format, high-performance focal plane detectors and, of course, especially for high-angular-resolution optics, provide a wealth of imaging data from extended sources. This review, follows the use of X-ray optics from the first rocket-borne instruments in the 1960s through to the Observatories flying today and being developed for future use. It also includes a brief overview of the challenges of fabricating X-ray optics and the various technologies that have been used to date, Comment: 45 pages, 21 figures. Invited chapter for the 'Handbook of X-ray and Gamma-ray Astrophysics (Eds C.Bambini and A Santangelo, Springer Singapore, expected publication in 2022)

We present the first results of a search for transient hard X-ray (HXR) emission in the quiet solar corona with the \textit{Nuclear Spectroscopic Telescope Array} (\textit{NuSTAR}) satellite. While \textit{NuSTAR} was designed as an astrophysics mission, it can observe the Sun above 2~keV with unprecedented sensitivity due to its pioneering use of focusing optics. \textit{NuSTAR} first observed quiet Sun regions on 2014 November 1, although out-of-view active regions contributed a notable amount of background in the form of single-bounce (unfocused) X-rays. We conducted a search for quiet Sun transient brightenings on time scales of 100 s and set upper limits on emission in two energy bands. We set 2.5--4~keV limits on brightenings with time scales of 100 s, expressed as the temperature T and emission measure EM of a thermal plasma. We also set 10--20~keV limits on brightenings with time scales of 30, 60, and 100 s, expressed as model-independent photon fluxes. The limits in both bands are well below previous HXR microflare detections, though not low enough to detect events of equivalent T and EM as quiet Sun brightenings seen in soft X-ray observations. We expect future observations during solar minimum to increase the \textit{NuSTAR} sensitivity by over two orders of magnitude due to higher instrument livetime and reduced solar background., Comment: 11 pages, 7 figures; accepted for publication in The Astrophysical Journal

The Nuclear Spectroscopic Telescope Array (NuSTAR), launched in June 2012, flies two conical approximation Wolter-I mirrors at the end of a 10.15m mast. The optics are coated with multilayers of Pt/C and W/Si that operate from 3--80 keV. Since the optical path is not shrouded, aperture stops are used to limit the field of view from background and sources outside the field of view. However, there is still a sliver of sky (~1.0--4.0 degrees) where photons may bypass the optics altogether and fall directly on the detector array. We term these photons Stray-light. Additionally, there are also photons that do not undergo the focused double reflections in the optics and we term these Ghost Rays. We present detailed analysis and characterization of these two components and discuss how they impact observations. Finally, we discuss how they could have been prevented and should be in future observatories., Comment: Published in Journal of Astronomical Telescopes, Instruments, and Systems. Open Access. http://dx.doi.org/10.1117/1.JATIS.3.4.044003

We present a hard X-ray observation of the TeV gamma-ray binary candidate HESS J1832-093 coincident with supernova remnant (SNR) G22.7-0.2 using the Nuclear Spectroscopic Telescope Array (NuSTAR). Non-thermal X-ray emission from XMMU J183245-0921539, the X-ray source associated with HESS J1832-093, is detected up to ~30 keV and is well-described by an absorbed power-law model with the best-fit photon index $\Gamma = 1.5\pm0.1$. A re-analysis of archival Chandra and XMM-Newton data finds that the long-term X-ray flux increase of XMMU J183245-0921539 is $50^{+40}_{-20}$% (90% C.L.), much less than previously reported. A search for a pulsar spin period or binary orbit modulation yields no significant signal to a pulse fraction limit of fp < 19% in the range 4 ms < P < 40 ks. No red noise is detected in the FFT power spectrum to suggest active accretion from a binary system. While further evidence is required, we argue that the X-ray and gamma-ray properties of XMMU J183245-0921539 are most consistent with a non-accreting binary generating synchrotron X- rays from particle acceleration in the shock formed as a result of the pulsar and stellar wind collision. We also report on three nearby hard X-ray sources, one of which may be associated with diffuse emission from a fast-moving supernova fragment interacting with a dense molecular cloud., Comment: ApJ in press. 9 pages

We report on two NuSTAR observations of the HMXB A 0535+26 taken toward the end of its normal 2015 outburst at very low $3-50$ keV luminosities of ${\sim}1.4\times10^{36}$ erg/s and ${\sim}5\times10^{35}$ erg/s which are complemented by 9 Swift observations. The data clearly confirm indications seen in earlier data that the source's spectral shape softens as it becomes fainter. The smooth, exponential rollover at high energies present in the first observation evolves to a much more abrupt steepening of the spectrum at $20-30$ keV. The continuum evolution can be well described with emission from a magnetized accretion column, modeled using the compmag model modified by an additional Gaussian emission component for the fainter observation. Between the two observations, the optical depth changes from $0.75\pm0.04$ to $0.56^{+0.01}_{-0.04}$, the electron temperature remains constant, and there is an indication that the column decreases in radius. Since the energy resolved pulse profiles remain virtually unchanged in shape between the two observations, the emission properties of the accretion column, however, reflect the same accretion regime. This conclusion is also confirmed by our result that the energy of the cyclotron resonant scattering feature (CRSF) at ${\sim}45$ keV is independent of the luminosity, implying that the magnetic field in the region in which the observed radiation is produced is the same in both observations. Finally, we also constrain the evolution of the continuum parameters with rotational phase of the neutron star. The width of the CRSF could only be constrained for the brighter observation. Based on Monte-Carlo simulations of CRSF formation in single accretion columns, its pulse phase dependence supports a simplified fan beam emission pattern. The evolution of the CRSF width is very similar to that of the CRSF depth, which is in disagreement with expectations., Comment: 14 pages, 11 figures, 3 tables, accepted for publication in A&A

Understanding the origin of the flaring activity from the Galactic center supermassive black hole, Sagittarius A*, is a major scientific goal of the NuSTAR Galactic plane survey campaign. We report on the data obtained between July 2012 and April 2015, including 27 observations on Sgr A* with a total exposure of ~ 1 Ms. We found a total of ten X-ray flares detected in the NuSTAR observation window, with luminosities in the range of $L_{3-79~keV}$~$(0.2-4.0) \times 10^{35}~erg~s^{-1}$. With this largest hard X-ray Sgr A* flare dataset to date, we studied the flare spectral properties. Seven flares are detected above 5{\sigma} significance, showing a range of photon indices ({\Gamma} ~ 2.0-2.8) with typical uncertainties of +/-0.5 (90% confidence level). We found no significant spectral hardening for brighter flares as indicated by a smaller sample. The accumulation of all the flare spectra in 1-79 keV can be well fit with an absorbed power-law model with {\Gamma}=2.2+/-0.1, and does not require the existence of a spectral break. The lack of variation in X-ray spectral index with luminosity would point to a single mechanism for the flares and is consistent with the synchrotron scenario. Lastly, we present the quiescent state spectrum of Sgr A*, and derived an upper limit on the quiescent luminosity of Sgr A* above 10 keV to be $L_{Xq, 10-79 keV}$ < $(2.9{\pm}0.2) \times 10^{34}~erg~s^{-1}$., Comment: 11 pages, 4 figures, Accepted for publication in ApJ

We present a catalog of hard X-ray sources in a square-degree region surveyed by NuSTAR in the direction of the Norma spiral arm. This survey has a total exposure time of 1.7 Ms, and typical and maximum exposure depths of 50 ks and 1 Ms, respectively. In the area of deepest coverage, sensitivity limits of $5\times10^{-14}$ and $4\times10^{-14}$ erg s$^{-1}$ cm$^{-2}$ in the 3-10 and 10-20 keV bands, respectively, are reached. Twenty-eight sources are firmly detected and ten are detected with low significance; eight of the 38 sources are expected to be active galactic nuclei. The three brightest sources were previously identified as a low-mass X-ray binary, high-mass X-ray binary, and pulsar wind nebula. Based on their X-ray properties and multi-wavelength counterparts, we identify the likely nature of the other sources as two colliding wind binaries, three pulsar wind nebulae, a black hole binary, and a plurality of cataclysmic variables (CVs). The CV candidates in the Norma region have plasma temperatures of $\approx$10-20 keV, consistent with the Galactic Ridge X-ray emission spectrum but lower than temperatures of CVs near the Galactic Center. This temperature difference may indicate that the Norma region has a lower fraction of intermediate polars relative to other types of CVs compared to the Galactic Center. The NuSTAR log$N$-log$S$ distribution in the 10-20 keV band is consistent with the distribution measured by Chandra at 2-10 keV if the average source spectrum is assumed to be a thermal model with $kT\approx15$~keV, as observed for the CV candidates., Comment: 42 pages, 12 figures, 11 tables

We present observations of the occulted active region AR12222 during the third {\em NuSTAR} solar campaign on 2014 December 11, with concurrent {\em SDO/}AIA and {\em FOXSI-2} sounding rocket observations. The active region produced a medium size solar flare one day before the observations, at $\sim18$UT on 2014 December 10, with the post-flare loops still visible at the time of {\em NuSTAR} observations. The time evolution of the source emission in the {\em SDO/}AIA $335\textrm{\AA}$ channel reveals the characteristics of an extreme-ultraviolet late phase event, caused by the continuous formation of new post-flare loops that arch higher and higher in the solar corona. The spectral fitting of {\em NuSTAR} observations yields an isothermal source, with temperature $3.8-4.6$ MK, emission measure $0.3-1.8 \times 10^{46}\textrm{ cm}^{-3}$, and density estimated at $2.5-6.0 \times 10^8 \textrm{ cm}^{-3}$. The observed AIA fluxes are consistent with the derived {\em NuSTAR} temperature range, favoring temperature values in the range $4.0-4.3$ MK. By examining the post-flare loops' cooling times and energy content, we estimate that at least 12 sets of post-flare loops were formed and subsequently cooled between the onset of the flare and {\em NuSTAR} observations, with their total thermal energy content an order of magnitude larger than the energy content at flare peak time. This indicates that the standard approach of using only the flare peak time to derive the total thermal energy content of a flare can lead to a large underestimation of its value., Comment: Published in ApJ, 8 pages

The distribution of elements produced in the inner-most layers of a supernova explosion is a key diagnostic for studying the collapse of massive stars. Here we present the results of a 2.4 Ms \textit{NuSTAR} observing campaign aimed at studying the supernova remnant Cassiopeia A (Cas A). We perform spatially-resolved spectroscopic analyses of the $^{44}$Ti ejecta which we use to determine the Doppler shift and thus the three-dimensional (3D) velocities of the $^{44}$Ti ejecta. We find an initial $^{44}$Ti mass of 1.54 $\pm$ 0.21 $\times 10^{-4}$ M$_{\odot}$ which has a present day average momentum direction of 340$^{\circ}$ $\pm$ 15$^{\circ}$ projected on to the plane of the sky (measured clockwise from Celestial North) and tilted by 58$^{\circ}$ $\pm$ 20$^{\circ}$ into the plane of the sky away from the observer, roughly opposite to the inferred direction of motion of the central compact object. We find some $^{44}$Ti ejecta that are clearly interior to the reverse shock and some that are clearly exterior to the reverse shock. Where we observe $^{44}$Ti ejecta exterior to the reverse shock we also see shock-heated iron; however, there are regions where we see iron but do not observe $^{44}$Ti. This suggests that the local conditions of the supernova shock during explosive nucleosynthesis varied enough to suppress the production of $^{44}$Ti in some regions by at least a factor of two, even in regions that are assumed to be the result of processes like $\alpha$-rich freezeout that should produce both iron and titanium., Comment: 14 pages, 13 figures (6 3D animations in the online journal). Accepted for publication in ApJ

We present results from the the first campaign of dedicated solar observations undertaken by the \textit{Nuclear Spectroscopic Telescope ARray} ({\em NuSTAR}) hard X-ray telescope. Designed as an astrophysics mission, {\em NuSTAR} nonetheless has the capability of directly imaging the Sun at hard X-ray energies ($>$3~keV) with an increase in sensitivity of at least two magnitude compared to current non-focusing telescopes. In this paper we describe the scientific areas where \textit{NuSTAR} will make major improvements on existing solar measurements. We report on the techniques used to observe the Sun with \textit{NuSTAR}, their limitations and complications, and the procedures developed to optimize solar data quality derived from our experience with the initial solar observations. These first observations are briefly described, including the measurement of the Fe K-shell lines in a decaying X-class flare, hard X-ray emission from high in the solar corona, and full-disk hard X-ray images of the Sun., Comment: 11 pages, accepted to ApJ

Recently, unresolved hard (20-40 keV) X-ray emission has been discovered within the central 10 pc of the Galaxy, possibly indicating a large population of intermediate polars (IPs). Chandra and XMM-Newton measurements in the surrounding ~50 pc imply a much lighter population of IPs with $\langle M_{\rm WD} \rangle \approx 0.5 M_\odot$. Here we use broad-band NuSTAR observations of two IPs: TV Columbae, which has a fairly typical but widely varying reported mass of $M_{\rm WD} \approx 0.5-1.0 M_\odot$, and IGR J17303-0601, with a heavy reported mass of $M_{\rm WD} \approx 1.0-1.2 M_\odot$. We investigate how varying spectral models and observed energy ranges influence estimated white dwarf mass. Observations of the inner 10 pc can be accounted for by IPs with $\langle M_{\rm WD} \rangle \approx 0.9 M_\odot$, consistent with that of the CV population in general, and the X-ray observed field IPs in particular. The lower mass derived by Chandra and XMM-Newton appears to be an artifact of narrow energy band fitting. To explain the (unresolved) CHXE by IPs requires an X-ray (2-8 keV) luminosity function (XLF) extending down to at least $5\times10^{31}$ erg/s. The CHXE XLF, if extended to the surrounding ~50 pc observed by Chandra and XMM-Newton, requires at least ~20-40% of the $\sim$9000 point sources are IPs. If the XLF extends just a factor of a few lower in luminosity, then the vast majority of these sources are IPs. This is in contrast to recent observations of the Galactic ridge, where the bulk of the 2-8 keV emission is ascribed to dwarf novae., Comment: 17 pages, accepted for publication in Astrophysical Journal

We present the first survey results of hard X-ray point sources in the Galactic Center (GC) region by NuSTAR. We have discovered 70 hard (3-79 keV) X-ray point sources in a 0.6 deg^2 region around Sgr A* with a total exposure of 1.7 Ms, and 7 sources in the Sgr B2 field with 300 ks. We identify clear Chandra counterparts for 58 NuSTAR sources and assign candidate counterparts for the remaining 19. The NuSTAR survey reaches X-ray luminosities of ~4 x and ~8 x 10^32 erg s^-1 at the GC (8 kpc) in the 3-10 and 10-40 keV bands, respectively. The source list includes three persistent luminous X-ray binaries and the likely run-away pulsar called the Cannonball. New source-detection significance maps reveal a cluster of hard (>10 keV) X-ray sources near the Sgr A diffuse complex with no clear soft X-ray counterparts. The severe extinction observed in the Chandra spectra indicates that all the NuSTAR sources are in the central bulge or are of extragalactic origin. Spectral analysis of relatively bright NuSTAR sources suggests that magnetic cataclysmic variables constitute a large fraction (>40-60%). Both spectral analysis and logN-logS distributions of the NuSTAR sources indicate that the X-ray spectra of the NuSTAR sources should have kT > 20 keV on average for a single temperature thermal plasma model or an average photon index of Gamma = 1.5 - 2 for a power-law model. These findings suggest that the GC X-ray source population may contain a larger fraction of X-ray binaries with high plasma temperatures than the field population., Comment: 32 pages, 15 figures, 7 tables, accepted for publication in ApJ

We report on a simultaneous NuSTAR and Swift observation of the neutron star low-mass X-ray binary 4U 1728-34. We identified and removed four Type I X-ray bursts during the observation in order to study the persistent emission. The continuum spectrum is hard and well described by a black body with $kT=$ 1.5 keV and a cutoff power law with $\Gamma=$ 1.5 and a cutoff temperature of 25 keV. Residuals between 6 and 8 keV provide strong evidence of a broad Fe K$\alpha$ line. By modeling the spectrum with a relativistically blurred reflection model, we find an upper limit for the inner disk radius of $R_{\rm in}\leq2 R_{\rm ISCO}$. Consequently we find that $R_{\rm NS}\leq23$ km, assuming $M=1.4{\mbox{$\rm\,M_{\mathord\odot}$}}$ and $a=0.15$. We also find an upper limit on the magnetic field of $B\leq2\times10^8$ G., Comment: 9 pages, 8 figures

The high-mass X-ray binary and accreting X-ray pulsar IGR J16393-4643 was observed by NuSTAR in the 3-79 keV energy band for a net exposure time of 50 ks. We present the results of this observation which enabled the discovery of a cyclotron resonant scattering feature with a centroid energy of 29.3(+1.1/-1.3) keV. This allowed us to measure the magnetic field strength of the neutron star for the first time: B = (2.5+/-0.1)e12 G. The known pulsation period is now observed at 904.0+/-0.1 s. Since 2006, the neutron star has undergone a long-term spin-up trend at a rate of P' = -2e-8 s/s (-0.6 s per year, or a frequency derivative of nu' = 3e-14 Hz/s ). In the power density spectrum, a break appears at the pulse frequency which separates the zero slope at low frequency from the steeper slope at high frequency. This addition of angular momentum to the neutron star could be due to the accretion of a quasi-spherical wind, or it could be caused by the transient appearance of a prograde accretion disk that is nearly in corotation with the neutron star whose magnetospheric radius is around 2e8 cm., Comment: Accepted for publication in the Astrophysical Journal, 7 pages, 8 figures, 2 tables

We report on the results of NuSTAR and XMM-Newton observations of the persistent X-ray source 1E1743.1-2843, located in the Galactic Center region. The source was observed between September and October 2012 by NuSTAR and XMM-Newton, providing almost simultaneous observations in the hard and soft X-ray bands. The high X-ray luminosity points to the presence of an accreting compact object. We analyze the possibilities of this accreting compact object being either a neutron star (NS) or a black hole, and conclude that the joint XMM-Newton and NuSTAR spectrum from 0.3 to 40 $\mathrm{keV}$ fits to a black body spectrum with $kT\sim1.8~\mathrm{keV}$ emitted from a hot spot or an equatorial strip on a neutron star surface. This spectrum is thermally Comptonized by electrons with $kT_{e}\sim4.6~\mathrm{keV}$. Accepting this neutron star hypothesis, we probe the Low Mass (LMXB) or High Mass (HMXB) X-ray Binary nature of the source. While the lack of Type-I bursts can be explained in the LMXB scenario, the absence of pulsations in the 2 mHz - 49 Hz frequency range, the lack of eclipses and of an IR companion, and the lack of a $K_{\alpha}$ line from neutral or moderately ionized iron strongly disfavor interpreting this source as a HMXB. We therefore conclude that 1E1743.1-2843 is most likely a NS-LMXB located beyond the Galactic Center. There is weak statistical evidence for a soft X-ray excess possibly indicating thermal emission from an accretion disk. However, the disk normalization remains unconstrained due to the high hydrogen column density ($N_{H}\sim1.6\times10^{23}~\mathrm{cm^{-2}}$)., Comment: Accepted for publication on ApJ on March 8th 2016

We present the results of a joint Swift-NuSTAR spectroscopy campaign on M31. We focus on the five brightest globular cluster X-ray sources in our fields. Two of these had previously been argued to be black hole candidates on the basis of apparent hard-state spectra at luminosities above those for which neutron stars are in hard states. We show that these two sources are likely to be Z-sources (i.e. low magnetic field neutron stars accreting near their Eddington limits), or perhaps bright atoll sources (low magnetic field neutron stars which are just a bit fainter than this level) on the basis of simultaneous Swift and NuSTAR spectra which cover a broader range of energies. These new observations reveal spectral curvature above 6-8 keV that would be hard to detect without the broader energy coverage the NuSTAR data provide relative to Chandra and XMM-Newton. We show that the other three sources are also likely to be bright neutron star X-ray binaries, rather than black hole X-ray binaries. We discuss why it should already have been realized that it was unlikely that these objects were black holes on the basis of their being persistent sources, and we re-examine past work which suggested that tidal capture products would be persistently bright X-ray emitters. We discuss how this problem is likely due to neglecting disk winds in older work that predict which systems will be persistent and which will be transient., Comment: 12 pages, 3 figures, accepted to MNRAS

We present NuSTAR and Swift observations of the neutron star Aquila X-1 during the peak of its July 2014 outburst. The spectrum is soft with strong evidence for a broad Fe K\alpha line. Modeled with a relativistically broadened reflection model, we find that the inner disk is truncated with an inner radius of 15+/-3 R_G. The disk is likely truncated by either the boundary layer and/or a magnetic field. Associating the truncated inner disk with pressure from a magnetic field gives an upper limit of B<5+/-2x10^8G. Although the radius is truncated far from the stellar surface, material is still reaching the neutron star surface as evidenced by the X-ray burst present in the t NuSTAR observation., Comment: ApJL Accepted, 8 pages, 4 figures

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

We present a hard X-ray NuSTAR observation of PSR J0437-4715, the nearest millisecond pulsar. The known pulsations at the apparent pulse period ~5.76 ms are detected at energies up to 20 keV. We measure a photon index $\Gamma= 1.65\pm0.24$ (90% confidence) for the power law fit to the non-thermal emission. It had been shown that spectral models with two or three thermal components fit the XMM-Newton spectrum of PSR J0437-4715, depending on the slope of the power-law component, and the amount of absorption of soft X-rays. The new constraint on the high-energy emission provided by NuSTAR removes ambiguities regarding the thermal components of the emission below 3 keV. We performed a simultaneous spectral analysis of the XMM-Newton and NuSTAR data to confirm that three thermal components and a power law are required to fit the 0.3-20 keV emission of PSR J0437-4715. Adding a ROSAT-PSPC spectrum further confirmed this result and allowed us to better constrain the temperatures of the three thermal components. A phase-resolved analysis of the NuSTAR data revealed no significant change in the photon index of the high-energy emission. This NuSTAR observation provides further impetus for future observations with the NICER mission (Neutron Star Interior Composition Explorer) whose sensitivity will provide much stricter constraints on the equation of state of nuclear matter by combining model fits to the pulsar's phase-folded lightcurve with the pulsar's well-defined mass and distance from radio timing observations., Comment: Accepted for publication in MNRAS. 12 pages, 8 figures, 3 table

This paper describes the Polarization Spectroscopic Telescope Array (PolSTAR), a mission proposed to NASA's 2014 Small Explorer (SMEX) announcement of opportunity. PolSTAR measures the linear polarization of 3-50 keV (requirement; goal: 2.5-70 keV) X-rays probing the behavior of matter, radiation and the very fabric of spacetime under the extreme conditions close to the event horizons of black holes, as well as in and around magnetars and neutron stars. The PolSTAR design is based on the technology developed for the Nuclear Spectroscopic Telescope Array (NuSTAR) mission launched in June 2012. In particular, it uses the same X-ray optics, extendable telescope boom, optical bench, and CdZnTe detectors as NuSTAR. The mission has the sensitivity to measure ~1% linear polarization fractions for X-ray sources with fluxes down to ~5 mCrab. This paper describes the PolSTAR design as well as the science drivers and the potential science return., Comment: (Astroparticle Physics in press, 34 pages, 23 figures, 6 tables)

We present the first sub-arcminute images of the Galactic Center above 10 keV, obtained with NuSTAR. NuSTAR resolves the hard X-ray source IGR J17456-2901 into non-thermal X-ray filaments, molecular clouds, point sources and a previously unknown central component of hard X-ray emission (CHXE). NuSTAR detects four non-thermal X-ray filaments, extending the detection of their power-law spectra with $\Gamma\sim1.3$-$2.3$ up to ~50 keV. A morphological and spectral study of the filaments suggests that their origin may be heterogeneous, where previous studies suggested a common origin in young pulsar wind nebulae (PWNe). NuSTAR detects non-thermal X-ray continuum emission spatially correlated with the 6.4 keV Fe K$\alpha$ fluorescence line emission associated with two Sgr A molecular clouds: MC1 and the Bridge. Broad-band X-ray spectral analysis with a Monte-Carlo based X-ray reflection model self-consistently determined their intrinsic column density ($\sim10^{23}$ cm$^{-2}$), primary X-ray spectra (power-laws with $\Gamma\sim2$) and set a lower limit of the X-ray luminosity of Sgr A* flare illuminating the Sgr A clouds to $L_X \stackrel{>}{\sim} 10^{38}$ erg s$^{-1}$. Above ~20 keV, hard X-ray emission in the central 10 pc region around Sgr A* consists of the candidate PWN G359.95-0.04 and the CHXE, possibly resulting from an unresolved population of massive CVs with white dwarf masses $M_{\rm WD} \sim 0.9 M_{\odot}$. Spectral energy distribution analysis suggests that G359.95-0.04 is likely the hard X-ray counterpart of the ultra-high gamma-ray source HESS J1745-290, strongly favoring a leptonic origin of the GC TeV emission., Comment: 27 pages. Accepted for publication in the Astrophysical Journal

We present a NuSTAR, Chandra, and XMM--Newton survey of nine of the nearest ultraluminous infrared galaxies (ULIRGs). The unprecedented sensitivity of NuSTAR at energies above 10 keV enables spectral modeling with far better precision than was previously possible. Six of the nine sources observed were detected sufficiently well by NuSTAR to model in detail their broadband X-ray spectra, and recover the levels of obscuration and intrinsic X-ray luminosities. Only one source (IRAS 13120--5453) has a spectrum consistent with a Compton--thick AGN, but we cannot rule out that a second source (Arp 220) harbors an extremely highly obscured AGN as well. Variability in column density (reduction by a factor of a few compared to older observations) is seen in IRAS 05189--2524 and Mrk 273, altering the classification of these border-line sources from Compton-thick to Compton-thin. The ULIRGs in our sample have surprisingly low observed fluxes in high energy (>10 keV) X-rays, especially compared to their bolometric luminosities. They have lower ratios of unabsorbed 2--10 keV to bolometric luminosity, and unabsorbed 2--10 keV to mid-IR [O IV] line luminosity than do Seyfert 1 galaxies. We identify IRAS 08572+3915 as another candidate intrinsically X-ray weak source, similar to Mrk 231. We speculate that the X-ray weakness of IRAS 08572+3915 is related to its powerful outflow observed at other wavelengths., Comment: 17 pages, 9 figures, 4 tables; Accepted for publication in the Astrophysical Journal

Galactic Center (GC) molecular cloud Sgr B2 is the best manifestation of an X-ray reflection nebula (XRN) reprocessing a past giant outburst from the supermassive black hole Sgr A*. Alternatively, Sgr B2 could be illuminated by low-energy cosmic ray electrons (LECRe) or protons (LECRp). In 2013, NuSTAR for the first time resolved Sgr B2 hard X-ray emission on sub-arcminute scales. Two prominent features are detected above 10 keV - a newly emerging cloud G0.66-0.13 and the central 90" radius region containing two compact cores Sgr B2(M) and Sgr B2(N) surrounded by diffuse emission. It is inconclusive whether the remaining level of Sgr B2 emission is still decreasing or has reached a constant background level. A decreasing Fe K$\alpha$ emission can be best explained by XRN while a constant background emission can be best explained by LECRp. In the XRN scenario, the 3-79 keV Sgr B2 spectrum can well constrain the past Sgr A* outburst, resulting in an outburst spectrum with a peak luminosity of $L_{3-79\rm~keV} \sim 5\times10^{38} \rm~erg~s^{-1}$ derived from the maximum Compton-scattered continuum and the Fe K$\alpha$ emission consistently. The XRN scenario is preferred by the fast variability of G0.66-0.13, which could be a molecular clump located in the Sgr B2 envelope reflecting the same Sgr A* outburst. In the LECRp scenario, we derived the required CR ion power $dW/dt=(1-4)\times10^{39}\rm~erg~s^{-1}$ and the CR ionization rate $\zeta_{H}=(6-10)\times 10^{-15}\rm~H^{-1}~s^{-1}$. The Sgr B2 background level X-ray emission will be a powerful tool to constrain GC CR population., Comment: 17 pages, 6 figures, submitted to ApJ

We present spectral and timing analysis of NuSTAR observations of the accreting X-ray pulsar 2RXP J130159.6-635806. The source was serendipitously observed during a campaign focused on the gamma-ray binary PSR B1259-63 and was later targeted for a dedicated observation. The spectrum has a typical shape for accreting X-ray pulsars, consisting of a simple power law with an exponential cutoff starting at ~7 keV with a folding energy of E_fold=~18 keV. There is also an indication of the presence of a 6.4 keV iron line in the spectrum at the ~3 sigma significance level. NuSTAR measurements of the pulsation period reveal that the pulsar has undergone a strong and steady spin-up for the last 20 years. The pulsed fraction is estimated to be ~80%, and is constant with energy up to 40 keV. The power density spectrum shows a break towards higher frequencies relative to the current spin period. This, together with steady persistent luminosity, points to a long-term mass accretion rate high enough to bring the pulsar out of spin equilibrium., Comment: 13 pages, 7 figures, accepted for publication in The Astrophysical Journal

We present results from multi-wavelength simultaneous X-ray and radio observations of the black hole X-ray binary V404 Cyg in quiescence. Our coverage with NuSTAR provides the very first opportunity to study the X-ray spectrum of V404 Cyg at energies above 10 keV. The unabsorbed broad-band (0.3--30 keV) quiescent luminosity of the source is 8.9$\times$10$^{32}$ erg s$^{-1}$ for a distance of 2.4 kpc. The source shows clear variability on short time scales (an hour to a couple of hours) in radio, soft X-ray and hard X-ray bands in the form of multiple flares. The broad-band X-ray spectra obtained from XMM-Newton and NuSTAR can be characterized with a power-law model having photon index $\Gamma$=2.12$\pm$0.07 (90\% confidence errors); however, residuals at high energies indicate spectral curvature significant at a 3$\sigma$ confidence level with e-folding energy of the cutoff to be 20$^{+20}_{-7}$ keV. Such curvature can be explained using synchrotron emission from the base of a jet outflow. Radio observations using the VLA reveal that the spectral index evolves on very fast time-scales (as short as 10 min.), switching between optically thick and thin synchrotron emission, possibly due to instabilities in the compact jet or stochastic instabilities in accretion rate. We explore different scenarios to explain this very fast variability., Comment: 16 pages, 11 figures and 3 tables, Accepted for publication in The Astrophysical Journal

We present NuSTAR observations of the powerful radio galaxy Cygnus A, focusing on the central absorbed active galactic nucleus (AGN). Cygnus A is embedded in a cool-core galaxy cluster, and hence we also examine archival XMM-Newton data to facilitate the decomposition of the spectrum into the AGN and intracluster medium (ICM) components. NuSTAR gives a source-dominated spectrum of the AGN out to >70keV. In gross terms, the NuSTAR spectrum of the AGN has the form of a power law (Gamma~1.6-1.7) absorbed by a neutral column density of N_H~1.6x10^23 cm^-2. However, we also detect curvature in the hard (>10keV) spectrum resulting from reflection by Compton-thick matter out of our line-of-sight to the X-ray source. Compton reflection, possibly from the outer accretion disk or obscuring torus, is required even permitting a high-energy cutoff in the continuum source; the limit on the cutoff energy is E_cut>111keV (90% confidence). Interestingly, the absorbed power-law plus reflection model leaves residuals suggesting the absorption/emission from a fast (15,000-26,000km/s), high column-density (N_W>3x10^23 cm^-2), highly ionized (xi~2,500 erg cm/s) wind. A second, even faster ionized wind component is also suggested by these data. We show that the ionized wind likely carries a significant mass and momentum flux, and may carry sufficient kinetic energy to exercise feedback on the host galaxy. If confirmed, the simultaneous presence of a strong wind and powerful jets in Cygnus A demonstrates that feedback from radio-jets and sub-relativistic winds are not mutually exclusive phases of AGN activity but can occur simultaneously., Comment: 13 pages; accepted for publication in The Astrophysical Journal

We present results from a 244\,ks \textit{NuSTAR} observation of 3C\,273 obtained during a cross-calibration campaign with the \textit{Chandra}, \textit{INTEGRAL}, \textit{Suzaku}, \textit{Swift}, and \textit{XMM-Newton} observatories. We show that the spectrum, when fit with a power-law model using data from all observatories except \textit{INTEGRAL} over the 1--78\,keV band, leaves significant residuals in the \textit{NuSTAR} data between 30--78\,keV. The \nustar\ 3--78\,keV spectrum is well-described by an exponentially cutoff power-law ($\Gamma = 1.646 \pm 0.006$, E$_\mathrm{cutoff} = 202_{-34}^{+51}$\,keV) with a weak reflection component from cold, dense material. There is also evidence for a weak ($EW = 23 \pm 11$ eV) neutral iron line. We interpret these features as arising from coronal emission plus reflection off an accretion disk or distant material. Beyond 80\,keV \textit{INTEGRAL} data show clear excess flux relative to an extrapolation of the AGN model fit to \nustar. This high-energy power-law is consistent with the presence of a beamed jet, which begins to dominate over emission from the inner accretion flow at 30-40 keV. Modeling the jet locally (in the \textit{NuSTAR} + \textit{INTEGRAL} band) as a power-law, we find the coronal component is fit by $\Gamma_\mathrm{AGN} = 1.638 \pm 0.045$, $E_\mathrm{cutoff} = 47 \pm 15$\,keV, and jet photon index by $\Gamma_\mathrm{jet} = 1.05 \pm 0.4$. We also consider \textit{Fermi}/LAT observations of 3C\,273 and here the broad-band spectrum of the jet can be described by a log-parabolic model, peaking at $\sim 2$\,MeV. Finally, we investigate the spectral variability in the \textit{NuSTAR} band and find an inverse correlation between flux and $\Gamma$., Comment: Accepted in ApJ

We present temporal and spectral analysis of simultaneous 0.5-79 keV Swift-XRT and NuSTAR observations of the magnetar 4U 0142+61. The pulse profile changes significantly with photon energy between 3 and 35 keV. The pulse fraction increases with energy, reaching a value of ~20%, similar to that observed in 1E 1841-045 and much lower than the ~80% pulse fraction observed in 1E 2259+586. We do not detect the 55-ks phase modulation reported in previous Suzaku-HXD observations. The phase-averaged spectrum of 4U 0142+61 above 20 keV is dominated by a hard power law with a photon index, $\Gamma$ ~ 0.65, and the spectrum below 20 keV can be described by two blackbodies, a blackbody plus a soft power law, or by a Comptonized blackbody model. We study the full phase-resolved spectra using the electron-positron outflow model of Beloborodov (2013). Our results are consistent with the parameters of the active j-bundle derived from INTEGRAL data by Hascoet et al. (2014). We find that a significant degeneracy appears in the inferred parameters if the footprint of the j-bundle is allowed to be a thin ring instead of a polar cap. The degeneracy is reduced when the footprint is required to be the hot spot inferred from the soft X-ray data., Comment: 14 pages, 8 figures, 4 tables. Accepted for publication in ApJ

The black hole candidate XTE J1908+094 went into outburst for the first time since 2003 in October 2013. We report on an observation with the Nuclear Spectroscopic Telescope Array (NuSTAR) and monitoring observations with Swift during the outburst. NuSTAR caught the source in the soft state: the spectra show a broad relativistic iron line, and the light curves reveal a ~40 ks flare with the count rate peaking about 40% above the non-flare level and with significant spectral variation. A model combining a multi-temperature thermal component, a power-law, and a reflection component with an iron line provides a good description of the NuSTAR spectrum. Although relativistic broadening of the iron line is observed, it is not possible to constrain the black hole spin with these data. The variability of the power-law component, which can also be modeled as a Comptonization component, is responsible for the flux and spectral change during the flare, suggesting that changes in the corona (or possibly continued jet activity) are the likely cause of the flare., Comment: 9 pages, 6 figures, 3 tables, accepted for publication in ApJ

We report on a 350-ks NuSTAR observation of the magnetar 1E 1841-045 taken in 2013 September. During the observation, NuSTAR detected six bursts of short duration, with $T_{90}<1$ s. An elevated level of emission tail is detected after the brightest burst, persisting for $\sim$1 ks. The emission showed a power-law decay with a temporal index of 0.5 before returning to the persistent emission level. The long observation also provided detailed phase-resolved spectra of the persistent X-ray emission of the source. By comparing the persistent spectrum with that previously reported, we find that the source hard-band emission has been stable over approximately 10 years. The persistent hard X-ray emission is well fitted by a coronal outflow model, where $e^{+/-}$ pairs in the magnetosphere upscatter thermal X-rays. Our fit of phase-resolved spectra allowed us to estimate the angle between the rotational and magnetic dipole axes of the magnetar, $\alpha_{mag}=0.25$, the twisted magnetic flux, $2.5\times10^{26}\rm \ G\ cm^2$, and the power released in the twisted magnetosphere, $L_j=6\times10^{36}\rm \ erg\ s^{-1}$. Assuming this model for the hard X-ray spectrum, the soft X-ray component is well fit by a two-blackbody model, with the hotter blackbody consistent with the footprint of the twisted magnetic field lines on the star. We also report on the 3-year Swift monitoring observations obtained since 2011 July. The soft X-ray spectrum remained stable during this period, and the timing behavior was noisy, with large timing residuals., Comment: 18 pages, 12 figures. Accepted for publication in ApJ

We report on NuSTAR, XMM-Newton and Swift observations of the gamma-ray binary 1FGL J1018.6-5856. We measure the orbital period to be 16.544+/-0.008 days using Swift data spanning 1900 days. The orbital period is different from the 2011 gamma-ray measurement which was used in the previous X-ray study of An et al. (2013) using ~400 days of Swift data, but is consistent with a new gamma-ray solution reported in 2014. The light curve folded on the new period is qualitatively similar to that reported previously, having a spike at phase 0 and broad sinusoidal modulation. The X-ray flux enhancement at phase 0 occurs more regularly in time than was previously suggested. A spiky structure at this phase seems to be a persistent feature, although there is some variability. Furthermore, we find that the source flux clearly correlates with the spectral hardness throughout all orbital phases, and that the broadband X-ray spectra measured with NuSTAR, XMM-Newton, and Swift are well fit with an unbroken power-law model. This spectrum suggests that the system may not be accretion-powered., Comment: 8 pages, 4 figures. Accepted for publication in ApJ

We report results from deep observations (~750 ks) of Tycho's supernova remnant (SNR) with NuSTAR. Using these data, we produce narrow-band images over several energy bands to identify the regions producing the hardest X-rays and to search for radioactive decay line emission from 44Ti. We find that the hardest (>10 keV) X-rays are concentrated in the southwest of Tycho, where recent Chandra observations have revealed high emissivity "stripes" associated with particles accelerated to the knee of the cosmic-ray spectrum. We do not find evidence of 44Ti, and we set limits on its presence and distribution within the SNR. These limits correspond to a upper-limit 44Ti mass of M44 < 2.4x10^-4 M_sun for a distance of 2.3 kpc. We perform spatially resolved spectroscopic analysis of sixty-six regions across Tycho. We map the best-fit rolloff frequency of the hard X-ray spectra, and we compare these results to measurements of the shock expansion and ambient density. We find that the highest energy electrons are accelerated at the lowest densities and in the fastest shocks, with a steep dependence of the roll-off frequency with shock velocity. Such a dependence is predicted by models where the maximum energy of accelerated electrons is limited by the age of the SNR rather than by synchrotron losses, but this scenario requires far lower magnetic field strengths than those derived from observations in Tycho. One way to reconcile these discrepant findings is through shock obliquity effects, and future observational work is necessary to explore the role of obliquity in the particle acceleration process., Comment: 12 pages, 12 figures, ApJ in press

We present the calibration of the \textit{Nuclear Spectroscopic Telescope Array} (\nustar) X-ray satellite. We used the Crab as the primary effective area calibrator and constructed a piece-wise linear spline function to modify the vignetting response. The achieved residuals for all off-axis angles and energies, compared to the assumed spectrum, are typically better than $\pm 2$\% up to 40\,keV and 5--10\,\% above due to limited counting statistics. An empirical adjustment to the theoretical 2D point spread function (PSF) was found using several strong point sources, and no increase of the PSF half power diameter (HPD) has been observed since the beginning of the mission. We report on the detector gain calibration, good to 60\,eV for all grades, and discuss the timing capabilities of the observatory, which has an absolute timing of $\pm$ 3\,ms. Finally we present cross-calibration results from two campaigns between all the major concurrent X-ray observatories (\textit{Chandra}, \textit{Swift}, \textit{Suzaku} and \textit{XMM-Newton}), conducted in 2012 and 2013 on the sources 3C\,273 and PKS\,2155-304, and show that the differences in measured flux is within $\sim$10\% for all instruments with respect to \nustar.

We present broadband (3 -- 78 keV) NuSTAR X-ray imaging and spectroscopy of the Crab nebula and pulsar. We show that while the phase-averaged and spatially integrated nebula + pulsar spectrum is a power-law in this energy band, spatially resolved spectroscopy of the nebula finds a break at $\sim$9 keV in the spectral photon index of the torus structure with a steepening characterized by $\Delta\Gamma\sim0.25$. We also confirm a previously reported steepening in the pulsed spectrum, and quantify it with a broken power-law with break energy at $\sim$12 keV and $\Delta\Gamma\sim0.27$. We present spectral maps of the inner 100\as\ of the remnant and measure the size of the nebula as a function of energy in seven bands. These results find that the rate of shrinkage with energy of the torus size can be fitted by a power-law with an index of $\gamma = 0.094\pm 0.018$, consistent with the predictions of Kennel and Coroniti (1984). The change in size is more rapid in the NW direction, coinciding with the counter-jet where we find the index to be a factor of two larger. NuSTAR observed the Crab during the latter part of a $\gamma$-ray flare, but found no increase in flux in the 3 - 78 keV energy band.

We present deep ($>$2.4 Ms) observations of the Cassiopeia A supernova remnant with {\it NuSTAR}, which operates in the 3--79 keV bandpass and is the first instrument capable of spatially resolving the remnant above 15 keV. We find that the emission is not entirely dominated by the forward shock nor by a smooth "bright ring" at the reverse shock. Instead we find that the $>$15 keV emission is dominated by knots near the center of the remnant and dimmer filaments near the remnant's outer rim. These regions are fit with unbroken power-laws in the 15--50 keV bandpass, though the central knots have a steeper ($\Gamma \sim -3.35$) spectrum than the outer filaments ($\Gamma \sim -3.06$). We argue this difference implies that the central knots are located in the 3-D interior of the remnant rather than at the outer rim of the remnant and seen in the center due to projection effects. The morphology of $>$15 keV emission does not follow that of the radio emission nor that of the low energy ($<$12 keV) X-rays, leaving the origin of the $>$15 keV emission as an open mystery. Even at the forward shock front we find less steepening of the spectrum than expected from an exponentially cut off electron distribution with a single cutoff energy. Finally, we find that the GeV emission is not associated with the bright features in the {\it NuSTAR} band while the TeV emission may be, suggesting that both hadronic and leptonic emission mechanisms may be at work., Comment: 12 pages, 11 figures, accepted for publication in ApJ

NuSTAR observed G1.9+0.3, the youngest known supernova remnant in the Milky Way, for 350 ks and detected emission up to $\sim$30 keV. The remnant's X-ray morphology does not change significantly across the energy range from 3 to 20 keV. A combined fit between NuSTAR and CHANDRA shows that the spectrum steepens with energy. The spectral shape can be well fitted with synchrotron emission from a power-law electron energy distribution with an exponential cutoff with no additional features. It can also be described by a purely phenomenological model such as a broken power-law or a power-law with an exponential cutoff, though these descriptions lack physical motivation. Using a fixed radio flux at 1 GHz of 1.17 Jy for the synchrotron model, we get a column density of N$_{\rm H}$ = $(7.23\pm0.07) \times 10^{22}$ cm$^{-2}$, a spectral index of $\alpha=0.633\pm0.003$, and a roll-off frequency of $\nu_{\rm rolloff}=(3.07\pm0.18) \times 10^{17}$ Hz. This can be explained by particle acceleration, to a maximum energy set by the finite remnant age, in a magnetic field of about 10 $\mu$G, for which our roll-off implies a maximum energy of about 100 TeV for both electrons and ions. Much higher magnetic-field strengths would produce an electron spectrum that was cut off by radiative losses, giving a much higher roll-off frequency that is independent of magnetic-field strength. In this case, ions could be accelerated to much higher energies. A search for $^{44}$Ti emission in the 67.9 keV line results in an upper limit of $1.5 \times 10^{-5}$ $\,\mathrm{ph}\,\mathrm{cm}^{-2}\,\mathrm{s}^{-1}$ assuming a line width of 4.0 keV (1 sigma)., Comment: 9 pages, 6 figures, accepted ApJ

We report on two NuSTAR observations of GRS 1741.9-2853, a faint neutron star low mass X-ray binary burster located 10' away from the Galactic center. NuSTAR detected the source serendipitously as it was emerging from quiescence: its luminosity was $6\times 10^{34}$ erg~s$^{-1}$ on 2013 July 31, and $5\times 10^{35}$ erg~s$^{-1}$ in a second observation on 2013 August 3. A bright, 800-s long, H-triggered mixed H/He thermonuclear Type I burst with mild photospheric radius expansion (PRE) was present during the second observation. Assuming that the luminosity during the PRE was at the Eddington level, a H mass fraction $X=0.7$ in the atmosphere, and a neutron star mass $M=1.4 M_{\odot}$, we determine a new lower limit on the distance for this source of $6.3 \pm 0.5$ kpc. Combining with previous upper limits, this places GRS 1741.9-2853 at a distance of 7 kpc. Energy independent (achromatic) variability is observed during the cooling of the neutron star, which could result from the disturbance of the inner accretion disk by the burst. The large dynamic range of this burst reveals a long power-law decay tail. We also detect, at a 95.6\% confidence level (1.7 $\sigma$), a narrow absorption line at $5.46\pm0.10$ keV during the PRE phase of the burst, reminiscent of the detection by Waki et al. (1984). We propose that the line, if real, is formed in the wind above the photosphere of the neutron star by a resonant K$\alpha$ transition from H-like Cr gravitationally redshifted by a factor $1+z=1.09$, corresponding to a radius range of 29.0 -- 41.4 km for a mass range of 1.4 -- 2.0 $M_{\odot}$., Comment: Accepted for publication in The Astrophysical Journal

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

We present spectral and timing analysis of NuSTAR observations of RX J0520.5$-$6932 in the 3-79 keV band collected during its outburst in January 2014. The target was observed on two epochs and we report the detection of a cyclotron resonant scattering feature with central energies of $E_\mathrm{CRSF} = 31.3_{-0.7}^{+0.8}$ keV and $31.5_{-0.6}^{+0.7}$ keV during the two observations, respectively, corresponding to a magnetic field of $B \approx 2 \times10^{12}$ G. The 3-79 keV luminosity of the system during the two epochs assuming a nominal distance of 50 kpc was $3.667\pm0.007\times 10^{38}\,\mathrm{erg\,s^{-1}}$ and $3.983\pm0.007\times10^{38}\,\mathrm{erg\,s^{-1}}$. Both values are much higher than the critical luminosity of $\approx1.5\times10^{37}\,\mathrm{erg\,s^{-1}}$ above which a radiation dominated shock front may be expected. This adds a new object to the sparse set of three systems that have a cyclotron line observed at luminosities in excess of $10^{38}\,\mathrm{erg\,s^{-1}}$. A broad ($\sigma\approx0.45$ keV) Fe emission line is observed in the spectrum at a central energy of $6.58_{-0.05}^{+0.05}$ keV in both epochs. The pulse profile of the pulsar was observed to be highly asymmetric with a sharply rising and slowly falling profile of the primary peak. We also observed minor variations in the cyclotron line energy and width as a function of the rotation phase.% As in observations of other cyclotron absorption line sources, there is a small ($\Delta\phi\lesssim0.1$) phase difference between the peak of the cyclotron energy variation and the peak of the flux variation., Comment: 9 pages, 4 figures, 6 tables. Accepted to ApJ

Timing of high-count rate sources with the NuSTAR Small Explorer Mission requires specialized analysis techniques. NuSTAR was primarily designed for spectroscopic observations of sources with relatively low count-rates rather than for timing analysis of bright objects. The instrumental dead time per event is relatively long (~2.5 msec), and varies by a few percent event-to-event. The most obvious effect is a distortion of the white noise level in the power density spectrum (PDS) that cannot be modeled easily with the standard techniques due to the variable nature of the dead time. In this paper, we show that it is possible to exploit the presence of two completely independent focal planes and use the cross power density spectrum to obtain a good proxy of the white noise-subtracted PDS. Thereafter, one can use a Monte Carlo approach to estimate the remaining effects of dead time, namely a frequency-dependent modulation of the variance and a frequency-independent drop of the sensitivity to variability. In this way, most of the standard timing analysis can be performed, albeit with a sacrifice in signal to noise relative to what would be achieved using more standard techniques. We apply this technique to NuSTAR observations of the black hole binaries GX 339-4, Cyg X-1 and GRS 1915+105., Comment: 13 pages, 8 figures, submitted to ApJ

We report on the first hard X-ray detection of the Geminga pulsar above 10 keV using a 150 ks observation with the NuSTAR observatory. The double-peaked pulse profile of non-thermal emission seen in the soft X-ray band persists at higher energies. Broadband phase-integrated spectra over the 0.2-20 keV band with NuSTAR and archival XMM-Newton data do not fit to a conventional two-component model of a blackbody plus power-law, but instead exhibit spectral hardening above ~5 keV. We find two spectral models fit the data well: (1) a blackbody (kT1 ~ 42 eV) with a broken power-law (Gamma1 ~ 2.0, Gamma2 ~ 1.4 and Ebreak ~ 3.4 keV), and (2) two blackbody components (kT1 ~ 44 eV and kT2 ~ 195 eV) with a power-law component (Gamma ~ 1.7). In both cases, the extrapolation of the Rayleigh-Jeans tail of the thermal component is consistent with the UV data, while the non-thermal component overpredicts the near-infrared data, requiring a spectral flattening at E ~ 0.01 - 1 keV. While strong phase variation of the power-law index is present below ~5 keV, our phase-resolved spectroscopy with NuSTAR indicates that another hard non-thermal component with Gamma ~ 1.3 emerges above ~5 keV. The spectral hardening in non-thermal X-ray emission as well as spectral flattening between the optical and X-ray bands argue against the conjecture that a single power-law may account for multi-wavelength non-thermal spectra of middle-aged pulsars., Comment: Accepted to ApJ

We report on new broad band spectral and temporal observations of the magnetar 1E 2259+586, which is located in the supernova remnant CTB 109. Our data were obtained simultaneously with the Nuclear Spectroscopic Telescope Array (NuSTAR) and Swift, and cover the energy range from 0.5-79 keV. We present pulse profiles in various energy bands and compare them to previous RXTE results. The NuSTAR data show pulsations above 20 keV for the first time and we report evidence that one of the pulses in the double-peaked pulse profile shifts position with energy. The pulsed fraction of the magnetar is shown to increase strongly with energy. Our spectral analysis reveals that the soft X-ray spectrum is well characterized by an absorbed double-blackbody or blackbody plus power-law model in agreement with previous reports. Our new hard X-ray data, however, suggests that an additional component, such as a power-law, is needed to describe the NuSTAR and Swift spectrum. We also fit the data with the recently developed coronal outflow model by Beloborodov for hard X-ray emission from magnetars. The outflow from a ring on the magnetar surface is statistically preferred over outflow from a polar cap., Comment: 37 pages, 9 figures, corresponding author, vogel9@llnl.gov

We present the first images of the pulsar wind nebula (PWN) MSH 15-52 in the hard X-ray band (>8 keV), as measured with the Nuclear Spectroscopic Telescope Array (NuSTAR). Overall, the morphology of the PWN as measured by NuSTAR in the 3-7 keV band is similar to that seen in Chandra high-resolution imaging. However, the spatial extent decreases with energy, which we attribute to synchrotron energy losses as the particles move away from the shock. The hard-band maps show a relative deficit of counts in the northern region towards the RCW 89 thermal remnant, with significant asymmetry. We find that the integrated PWN spectra measured with NuSTAR and Chandra suggest that there is a spectral break at 6 keV which may be explained by a break in the synchrotron-emitting electron distribution at ~200 TeV and/or imperfect cross calibration. We also measure spatially resolved spectra, showing that the spectrum of the PWN softens away from the central pulsar B1509-58, and that there exists a roughly sinusoidal variation of spectral hardness in the azimuthal direction. We discuss the results using particle flow models. We find non-monotonic structure in the variation with distance of spectral hardness within 50" of the pulsar moving in the jet direction, which may imply particle and magnetic-field compression by magnetic hoop stress as previously suggested for this source. We also present 2-D maps of spectral parameters and find an interesting shell-like structure in the NH map. We discuss possible origins of the shell-like structure and their implications., Comment: 15 pages, 9 figures, accepted for publication in ApJ

We present a broadband (~0.5-79 keV) spectral and temporal analysis of multiple NuSTAR observations combined with archival Suzaku and Chandra data of NGC4945, the brightest extragalactic source at 100 keV. We observe hard X-ray (> 10 keV) flux and spectral variability, with flux variations of a factor 2 on timescales of 20 ksec. A variable primary continuum dominates the high energy spectrum (>10 keV) in all the states, while the reflected/scattered flux which dominates at E< 10 keV stays approximately constant. From modelling the complex reflection/transmission spectrum we derive a Compton depth along the line of sight of tau_Thomson ~ 2.9, and a global covering factor for the circumnuclear gas of ~ 0.15. This agrees with the constraints derived from the high energy variability, which implies that most of the high energy flux is transmitted, rather that Compton-scattered. This demonstrates the effectiveness of spectral analysis in constraining the geometric properties of the circumnuclear gas, and validates similar methods used for analyzing the spectra of other bright, Compton-thick AGN. The lower limits on the e-folding energy are between 200-300 keV, consistent with previous BeppoSAX, Suzaku and Swift BAT observations. The accretion rate, estimated from the X-ray luminosity and assuming a bolometric correction typical of type 2 AGN, is in the range ~0.1-0.3 lambda_Edd depending on the flux state. The substantial observed X-ray luminosity variability of NGC4945 implies that large errors can arise from using single-epoch X-ray data to derive L/L_Edd values for obscured AGNs., Comment: Accepted for publication in ApJ

We present NuSTAR spectral and timing studies of the Supergiant Fast X-ray Transient (SFXT) IGR J17544-2619. The spectrum is well-described by a ~1 keV blackbody and a hard continuum component, as expected from an accreting X-ray pulsar. We detect a cyclotron line at 17 keV, confirming that the compact object in IGR J17544-2619 is indeed a neutron star. This is the first measurement of the magnetic field in a SFXT. The inferred magnetic field strength, B = (1.45 +/- 0.03) * 10^12 G * (1+z) is typical of neutron stars in X-ray binaries, and rules out a magnetar nature for the compact object. We do not find any significant pulsations in the source on time scales of 1-2000 s., Comment: MNRAS Accepted. 8 pages, 4 figures, 3 tables

Results are presented for an initial survey of the Norma Arm gathered with the focusing hard X-ray telescope NuSTAR. The survey covers 0.2 deg$^2$ of sky area in the 3-79 keV range with a minimum and maximum raw depth of 15 ks and 135 ks, respectively. Besides a bright black-hole X-ray binary in outburst (4U 1630-47) and a new X-ray transient (NuSTAR J163433-473841), NuSTAR locates three sources from the Chandra survey of this region whose spectra are extended above 10 keV for the first time: CXOU J163329.5-473332, CXOU J163350.9-474638, and CXOU J163355.1-473804. Imaging, timing, and spectral data from a broad X-ray range (0.3-79 keV) are analyzed and interpreted with the aim of classifying these objects. CXOU J163329.5-473332 is either a cataclysmic variable or a faint low-mass X-ray binary. CXOU J163350.9-474638 varies in intensity on year-long timescales, and with no multi-wavelength counterpart, it could be a distant X-ray binary or possibly a magnetar. CXOU J163355.1-473804 features a helium-like iron line at 6.7 keV and is classified as a nearby cataclysmic variable. Additional surveys are planned for the Norma Arm and Galactic Center, and those NuSTAR observations will benefit from the lessons learned during this pilot study., Comment: Accepted for publication in the Astrophysical Journal, 10 pages, 4 figures, 3 tables

We present results of the point spread function (PSF) calibration of the hard X-ray optics of the Nuclear Spectroscopic Telescope Array (NuSTAR). Immediately post-launch, NuSTAR has observed bright point sources such as Cyg X-1, Vela X-1, and Her X-1 for the PSF calibration. We use the point source observations taken at several off-axis angles together with a ray-trace model to characterize the in-orbit angular response, and find that the ray-trace model alone does not fit the observed event distributions and applying empirical corrections to the ray-trace model improves the fit significantly. We describe the corrections applied to the ray-trace model and show that the uncertainties in the enclosed energy fraction (EEF) of the new PSF model is < 3% for extraction apertures of R > 60" with no significant energy dependence. We also show that the PSF of the NuSTAR optics has been stable over a period of ~300 days during its in-orbit operation., Comment: 10 pages, 6 figures. Presented at the SPIE conference Astronomical Telescopes + Instrumentation 2014

We report NuSTAR observations of the millisecond pulsar - low mass X-ray binary (LMXB) transition system PSR J1023+0038 from June and October 2013, before and after the formation of an accretion disk around the neutron star. Between June 10-12, a few days to two weeks before the radio disappearance of the pulsar, the 3-79 keV X-ray spectrum was well fit by a simple power law with a photon index of Gamma=1.17 +/-0.08 (at 90% confidence) with a 3-79 keV luminosity of 7.4+/-0.4 x 10^32 erg/s. Significant orbital modulation was observed with a modulation fraction of 36+/-10%. During the October 19-21 observation, the spectrum is described by a softer power law (Gamma=1.66+/-0.06) with an average luminosity of 5.8+/-0.2 x 10^33 erg/s and a peak luminosity of ~1.2 x 10^34 erg/s observed during a flare. No significant orbital modulation was detected. The spectral observations are consistent with previous and current multi-wavelength observations and show the hard X-ray power law extending to 79 keV without a spectral break. Sharp edged, flat bottomed `dips' are observed with widths between 30-1000 s and ingress and egress time-scales of 30-60 s. No change in hardness ratio was observed during the dips. Consecutive dip separations are log-normal in distribution with a typical separation of approximately 400 s. These dips are distinct from dipping activity observed in LMXBs. We compare and contrast these dips to observations of dips and state changes in the similar transition systems PSR J1824-2452I and XSS J1227.0-4859 and discuss possible interpretations based on the transitions in the inner disk., Comment: 13 pages, 17 figures, accepted for publication in ApJ