• The main characteristics in the radio continuum spectra of Supernova Remnants (SNRs) result from simple synchrotron emission. In addition, electron acceleration mechanisms can shape the spectra in specific ways, especially at high radio frequencies. These features are connected to the age and the peculiar conditions of the local interstellar medium interacting with the SNR. Whereas the bulk radio emission is expected at up to $20-50$ GHz, sensitive high-resolution images of SNRs above 10 GHz are lacking and are not easily achievable, especially in the confused regions of the Galactic Plane. In the framework of the early science observations with the Sardinia Radio Telescope in February-March 2016, we obtained high-resolution images of SNRs Tycho, W44 and IC443 that provided accurate integrated flux density measurements at 21.4 GHz: 8.8 $\pm$ 0.9 Jy for Tycho, 25 $\pm$ 3 Jy for W44 and 66 $\pm$ 7 Jy for IC443. We coupled the SRT measurements with radio data available in the literature in order to characterise the integrated and spatially-resolved spectra of these SNRs, and to find significant frequency- and region-dependent spectral slope variations. For the first time, we provide direct evidence of a spectral break in the radio spectral energy distribution of W44 at an exponential cutoff frequency of 15 $\pm$ 2 GHz. This result constrains the maximum energy of the accelerated electrons in the range $6-13$ GeV, in agreement with predictions indirectly derived from AGILE and \textit{Fermi}-LAT gamma-ray observations. With regard to IC443, our results confirm the noticeable presence of a bump in the integrated spectrum around $20-70$ GHz that could result from a spinning dust emission mechanism.
  • We report on a new NuSTAR observation and on the ongoing Swift XRT monitoring campaign of the peculiar source 1E 161348-5055, located at the centre of the supernova remnant RCW 103, which is recovering from its last outburst in June 2016. The X-ray spectrum at the epoch of the NuSTAR observation can be described by either two absorbed blackbodies ($kT_{BB_1}$ ~ 0.5 keV, $kT_{BB_2}$ ~ 1.2 keV) or an absorbed blackbody plus a power law ($kT_{BB_1}$ ~ 0.6 keV, $\Gamma$ ~ 3.9). The observed flux was ~ 9 $\times$ 10$^{-12}$ erg s$^{-1}$ cm$^{-2}$, ~ 3 times lower than what observed at the outburst onset, but about one order of magnitude higher than the historical quiescent level. A periodic modulation was detected at the known 6.67 hr periodicity. The spectral decomposition and evolution along the outburst decay are consistent with 1E 161348-5055 being a magnetar, the slowest ever detected.
  • Based on phase-resolved broadband spectroscopy using $XMM$-$Newton$ and $NuSTAR$, we report on a potential cyclotron resonant scattering feature at $E \sim 13$ keV in the pulsed spectrum of the recently discoverd ULX pulsar NGC 300 ULX1. If this interpretation is correct, the implied magnetic field of the central neutron star is $B \sim 10^{12}$ G (assuming scattering off electrons), similar to that estimated from the observed spin-up of the star, and also similar to known Galactic X-ray pulsars. We discuss the implications of this result for the connection between NGC 300 ULX1 and the other known ULX pulsars, particularly in light of the recent discovery of a likely proton Cyclotron line in another ULX, M51 ULX-8.
  • We present broadband X-ray analyses of a sample of bright ultraluminous X-ray sources with the goal of investigating the spectral similarity of this population to the known ULX pulsars, M82 X-2, NGC7793 P13 and NGC5907 ULX. We perform a phase-resolved analysis of the broadband XMM-Newton+NuSTAR dataset of NGC5907 ULX, finding that the pulsed emission from the accretion column in this source exhibits a similar spectral shape to that seen in both M82 X-2 and NGC7793 P13, and that this is responsible for the excess emission observed at the highest energies when the spectra are fit with accretion disk models. We then demonstrate that similar 'hard' excesses are seen in all the ULXs in the broadband sample. Finally, for the ULXs where the nature of the accretor is currently unknown, we test whether the hard excesses are all consistent with being produced by an accretion column similar to those present in M82 X-2, NGC7793 P13 and NGC5907 ULX. Based on the average shape of the pulsed emission, we find that in all cases a similar accretion column can successfully reproduce the observed data, consistent with the hypothesis that this ULX sample may be dominated by neutron star accretors. Compared to the known pulsar ULXs, our spectral fits for the remaining ULXs suggest that the non-pulsed emission from the accretion flow beyond the magnetosphere makes a stronger relative contribution than the component associated with the accretion column. If these sources do also contain neutron star accretors, this may help to explain the lack of detected pulsations.
  • Ultraluminous pulsars are a definite proof that persistent super-Eddington accretion occurs in nature. They support the scenario according to which most Ultraluminous X-ray Sources (ULXs) are super-Eddington accretors of stellar mass rather than sub-Eddington intermediate mass black holes. An important prediction of theories of supercritical accretion is the existence of powerful outflows of moderately ionized gas at mildly relativistic speeds. In practice, the spectral resolution of X-ray gratings such as RGS onboard XMM-Newton is required to resolve their observational signatures in ULXs. Using RGS, outflows have been discovered in the spectra of 3 ULXs (none of which are currently known to be pulsars). Most recently, the fourth ultraluminous pulsar was discovered in NGC 300. Here we report detection of an ultrafast outflow (UFO) in the X-ray spectrum of the object, with a significance of at least 3{\sigma}, during one of the two simultaneous observations of the source by XMM-Newton and NuSTAR in December 2016. The outflow has a projected velocity of about 72000 km/s (0.24c) and a high ionisation factor with a log value of 3.9. This is the first time a detection of an ultrafast wind in a ULX was achieved using both soft and hard X-ray data simultaneously, and it is also the first detection of a UFO in a neutron star ULX. We find no evidence of the UFO during the other observation of the object, which could be explained by either clumpy nature of the absorber or a slight change in our viewing angle of the accretion flow.
  • The Astropy Collaboration, A. M. Price-Whelan, B. M. Sipőcz, H. M. Günther, P. L. Lim, S. M. Crawford, S. Conseil, D. L. Shupe, M. W. Craig, N. Dencheva, A. Ginsburg, J. T. VanderPlas, L. D. Bradley, D. Pérez-Suárez, M. de Val-Borro, T. L. Aldcroft, K. L. Cruz, T. P. Robitaille, E. J. Tollerud, C. Ardelean, T. Babej, M. Bachetti, A. V. Bakanov, S. P. Bamford, G. Barentsen, P. Barmby, A. Baumbach, K. L. Berry, F. Biscani, M. Boquien, K. A. Bostroem, L. G. Bouma, G. B. Brammer, E. M. Bray, H. Breytenbach, H. Buddelmeijer, D. J. Burke, G. Calderone, J. L. Cano Rodríguez, M. Cara, J. V. M. Cardoso, S. Cheedella, Y. Copin, D. Crichton, D. DÁvella, C. Deil, É. Depagne, J. P. Dietrich, A. Donath, M. Droettboom, N. Earl, T. Erben, S. Fabbro, L. A. Ferreira, T. Finethy, R. T. Fox, L. H. Garrison, S. L. J. Gibbons, D. A. Goldstein, R. Gommers, J. P. Greco, P. Greenfield, A. M. Groener, F. Grollier, A. Hagen, P. Hirst, D. Homeier, A. J. Horton, G. Hosseinzadeh, L. Hu, J. S. Hunkeler, Ž. Ivezić, A. Jain, T. Jenness, G. Kanarek, S. Kendrew, N. S. Kern, W. E. Kerzendorf, A. Khvalko, J. King, D. Kirkby, A. M. Kulkarni, A. Kumar, A. Lee, D. Lenz, S. P. Littlefair, Z. Ma, D. M. Macleod, M. Mastropietro, C. McCully, S. Montagnac, B. M. Morris, M. Mueller, S. J. Mumford, D. Muna, N. A. Murphy, S. Nelson, G. H. Nguyen, J. P. Ninan, M. Nöthe, S. Ogaz, S. Oh, J. K. Parejko, N. Parley, S. Pascual, R. Patil, A. A. Patil, A. L. Plunkett, J. X. Prochaska, T. Rastogi, V. Reddy Janga, J. Sabater, P. Sakurikar, M. Seifert, L. E. Sherbert, H. Sherwood-Taylor, A. Y. Shih, J. Sick, M. T. Silbiger, S. Singanamalla, L. P. Singer, P. H. Sladen, K. A. Sooley, S. Sornarajah, O. Streicher, P. Teuben, S. W. Thomas, G. R. Tremblay, J. E. H. Turner, V. Terrón, M. H. van Kerkwijk, A. de la Vega, L. L. Watkins, B. A. Weaver, J. B. Whitmore, J. Woillez, V. Zabalza
    Jan. 16, 2018 astro-ph.IM
    The Astropy project supports and fosters the development of open-source and openly-developed Python packages that provide commonly-needed functionality to the astronomical community. A key element of the Astropy project is the core package Astropy, which serves as the foundation for more specialized projects and packages. In this article, we provide an overview of the organization of the Astropy project and summarize key features in the core package as of the recent major release, version 2.0. We then describe the project infrastructure designed to facilitate and support development for a broader ecosystem of inter-operable packages. We conclude with a future outlook of planned new features and directions for the broader Astropy project.
  • In recent years, the cross spectrum has received considerable attention as a means of characterising the variability of astronomical sources as a function of wavelength. While much has been written about the statistics of time and phase lags, the cospectrum has only recently been understood as means of mitigating instrumental effects dependent on temporal frequency in astronomical detectors, as well as a method of characterizing the coherent variability in two wavelength ranges on different time scales. In this paper, we lay out the statistical foundations of the cospectrum, starting with the simplest case of detecting a periodic signal in the presence of white noise. This case is especially relevant for detecting faint X-ray pulsars in detectors heavily affected by instrumental effects, including NuSTAR, Astrosat and IXPE. We show that the statistical distributions of both single and averaged cospectra differ considerably from those for standard periodograms. While a single cospectrum follows a Laplace distribution exactly, averaged cospectra are approximated by a Gaussian distribution only for more than ~30 averaged segments, dependent on the number of trials. We provide an instructive example of a quasi-periodic oscillation in NuSTAR and show that applying standard periodogram statistics leads to underestimated tail probabilities for period detection. We also demonstrate the application of these distributions to a NuSTAR observation of the X-ray pulsar Hercules X-1.
  • The presence of neutron stars in at least three ultraluminous X-ray sources is now firmly established and offers an unambiguous view of super-critical accretion. All three systems show long-timescale periods (60-80 days) in the X-rays and/or optical, two of which are known to be super-orbital in nature. Should the flow be classically super critical, i.e. the Eddington limit is reached locally in the disc (implying surface dipole fields that are sub-magnetar in strength), then the large scale-height flow can precess through the Lense-Thirring effect which could provide an explanation for the observed super-orbital periods. By connecting the details of the Lense-Thirring effect with the observed pulsar spin period, we are able to infer the moment-of-inertia and therefore equation-of-state of the neutron star without relying on the inclination of, or distance to the system. We apply our technique to the case of NGC 7793 P13 and demonstrate that stronger magnetic fields imply stiffer equations of state. We discuss the caveats and uncertainties, many of which can be addressed through forthcoming radiative magnetohydrodynamic (RMHD) simulations and their connection to observation.
  • In September 2016, the microquasar Cygnus X-3 underwent a giant radio flare, which was monitored for 6 days with the Medicina Radio Astronomical Station and the Sardinia Radio Telescope. Long observations were performed in order to follow the evolution of the flare on a hourly scale, covering six frequency ranges from 1.5 GHz to 25.6 GHz. The radio emission reached a maximum of 13.2 +/- 0.7 Jy at 7.2 GHz and 10 +/- 1 Jy at 18.6 GHz. Rapid flux variations were observed at high radio frequencies at the peak of the flare, together with rapid evolution of the spectral index: alpha steepened from 0.3 to 0.6 within 5 hours. This is the first time that such fast variations are observed, giving support to the evolution from optically thick to optically thin plasmons in expansion moving outward from the core. Based on the Italian network (Noto, Medicina and SRT) and extended to the European antennas (Torun, Yebes, Onsala), VLBI observations were triggered at 22 GHz on five different occasions, four times prior to the giant flare, and once during its decay phase. Flux variations of 2-hour duration were recorded during the first session. They correspond to a mini-flare that occurred close to the core ten days before the onset of the giant flare. From the latest VLBI observation we infer that four days after the flare peak the jet emission was extended over 30 mas.
  • We present a detailed, broadband X-ray spectral analysis of the ULX pulsar NGC 7793 P13, a known super-Eddington source, utilizing data from the $XMM$-$Newton$, $NuSTAR$ and $Chandra$ observatories. The broadband $XMM$-$Newton+NuSTAR$ spectrum of P13 is qualitatively similar to the rest of the ULX sample with broadband coverage, suggesting that additional ULXs in the known population may host neutron star accretors. Through time-averaged, phase-resolved and multi-epoch studies, we find that two non-pulsed thermal blackbody components with temperatures $\sim$0.5 and $\sim$1.5 keV are required to fit the data below 10 keV, in addition to a third continuum component which extends to higher energies and is associated with the pulsed emission from the accretion column. The characteristic radii of the thermal components appear to be similar, and are too large to be associated with the neutron star itself, so the need for two components likely indicates the accretion flow outside the magnetosphere is complex. We suggest a scenario in which the thick inner disc expected for super-Eddington accretion begins to form, but is terminated by the neutron star's magnetic field soon after its onset, implying a magnetic field of $B \lesssim 6 \times 10^{12}$ G for the central neutron star. Evidence of similar termination of the disc in other sources may offer a further means of identifying additional neutron star ULXs. Finally, we examine the spectrum exhibited by P13 during one of its unusual 'off' states. These data require both a hard powerlaw component, suggesting residual accretion onto the neutron star, and emission from a thermal plasma, which we argue is likely associated with the P13 system.
  • Observations of supernova remnants (SNRs) are a powerful tool for investigating the later stages of stellar evolution, the properties of the ambient interstellar medium, and the physics of particle acceleration and shocks. For a fraction of SNRs, multi-wavelength coverage from radio to ultra high-energies has been provided, constraining their contributions to the production of Galactic cosmic rays. Although radio emission is the most common identifier of SNRs and a prime probe for refining models, high-resolution images at frequencies above 5 GHz are surprisingly lacking, even for bright and well-known SNRs such as IC443 and W44. In the frameworks of the Astronomical Validation and Early Science Program with the 64-m single-dish Sardinia Radio Telescope, we provided, for the first time, single-dish deep imaging at 7 GHz of the IC443 and W44 complexes coupled with spatially-resolved spectra in the 1.5-7 GHz frequency range. Our images were obtained through on-the-fly mapping techniques, providing antenna beam oversampling and resulting in accurate continuum flux density measurements. The integrated flux densities associated with IC443 are S_1.5GHz = 134 +/- 4 Jy and S_7GHz = 67 +/- 3 Jy. For W44, we measured total flux densities of S_1.5GHz = 214 +/- 6 Jy and S_7GHz = 94 +/- 4 Jy. Spectral index maps provide evidence of a wide physical parameter scatter among different SNR regions: a flat spectrum is observed from the brightest SNR regions at the shock, while steeper spectral indices (up to 0.7) are observed in fainter cooling regions, disentangling in this way different populations and spectra of radio/gamma-ray-emitting electrons in these SNRs.
  • [Abridged] The Sardinia Radio Telescope (SRT) is the new 64-m dish operated by INAF (Italy). Its active surface will allow us to observe at frequencies of up to 116 GHz. At the moment, three receivers, one per focal position, have been installed and tested. The SRT was officially opened in October 2013, upon completion of its technical commissioning phase. In this paper, we provide an overview of the main science drivers for the SRT, describe the main outcomes from the scientific commissioning of the telescope, and discuss a set of observations demonstrating the SRT's scientific capabilities. One of the main objectives of scientific commissioning was the identification of deficiencies in the instrumentation and/or in the telescope sub-systems for further optimization. As a result, the overall telescope performance has been significantly improved. As part of the scientific commissioning activities, different observing modes were tested and validated, and first astronomical observations were carried out to demonstrate the science capabilities of the SRT. In addition, we developed astronomer-oriented software tools, to support future observers on-site. The astronomical validation activities were prioritized based on technical readiness and scientific impact. The highest priority was to make the SRT available for joint observations as part of European networks. As a result, the SRT started to participate (in shared-risk mode) in EVN (European VLBI Network) and LEAP (Large European Array for Pulsars) observing sessions in early 2014. The validation of single-dish operations for the suite of SRT first light receivers and backends continued in the following years, and was concluded with the first call for shared-risk/early-science observations issued at the end of 2015.
  • We present results from four new broadband X-ray observations of the extreme ultraluminous X-ray source Holmberg IX X-1 ($L_{\rm{X}} > 10^{40}$ erg/s), performed by $Suzaku$ and $NuSTAR$ in coordination. Combined with the archival data, we now have broadband observations of this remarkable source from six separate epochs. Two of these new observations probe lower fluxes than seen previously, allowing us to extend our knowledge of the broadband spectral variability exhibited. The spectra are well fit by two thermal blackbody components, which dominate the emission below 10 keV, as well as a steep ($\Gamma \sim 3.5$) powerlaw tail which dominates above $\sim$15 keV. Remarkably, while the 0.3-10.0 keV flux varies by a factor of $\sim$3 between all these epochs, the 15-40 keV flux varies by only $\sim$20%. Although the spectral variability is strongest in the $\sim$1-10 keV band, both of the thermal components are required to vary when all epochs are considered. We also re-visit the search for iron absorption features, leveraging the high-energy $NuSTAR$ data to improve our sensitivity to extreme velocity outflows in light of the ultra-fast outflow recently detected in NGC 1313 X-1. Iron absorption from a similar outflow along our line of sight can be ruled out in this case. We discuss these results in the context of super-Eddington accretion models that invoke a funnel-like geometry for the inner flow, and propose a scenario in which we have an almost face-on view of a funnel that expands to larger radii with increasing flux, resulting in an increasing degree of geometrical collimation for the emission from intermediate temperature regions.
  • We present first results from a series of $NuSTAR$ observations of the black hole X-ray binary V404 Cyg obtained during its summer 2015 outburst, primarily focusing on observations during the height of this outburst activity. The $NuSTAR$ data show extreme variability in both the flux and spectral properties of the source. This is partly driven by strong and variable line-of-sight absorption, similar to previous outbursts. The latter stages of this observation are dominated by strong flares, reaching luminosities close to Eddington. During these flares, the central source appears to be relatively unobscured and the data show clear evidence for a strong contribution from relativistic reflection, providing a means to probe the geometry of the innermost accretion flow. Based on the flare properties, analogy with other Galactic black hole binaries, and also the simultaneous onset of radio activity, we argue that this intense X-ray flaring is related to transient jet activity during which the ejected plasma is the primary source of illumination for the accretion disk. If this is the case, then our reflection modelling implies that these jets are launched in close proximity to the black hole (as close as a few gravitational radii), consistent with expectations for jet launching models that tap either the spin of the central black hole, or the very innermost accretion disk. Our analysis also allows us to place the first constraints on the black hole spin for this source, which we find to be $a^* > 0.92$ (99% statistical uncertainty, based on an idealized lamppost geometry).
  • In June 2015, the black hole X-ray binary (BHXRB) V404 Cygni went into outburst for the first time since 1989. Here, we present a comprehensive search for quasi-periodic oscillations (QPOs) of V404 Cygni during its recent outburst, utilizing data from six instruments on board five different X-ray missions: Swift/XRT, Fermi/GBM, Chandra/ACIS, INTEGRAL's IBIS/ISGRI and JEM-X, and NuSTAR. We report the detection of a QPO at 18 mHz simultaneously with both Fermi/GBM and Swift/XRT, another example of a rare but slowly growing new class of mHz-QPOs in BHXRBs linked to sources with a high orbital inclination. Additionally, we find a duo of QPOs in a Chandra/ACIS observation at 73 mHz and 1.03 Hz, as well as a QPO at 136 mHz in a single Swift/XRT observation that can be interpreted as standard Type-C QPOs. Aside from the detected QPOs, there is significant structure in the broadband power, with a strong feature observable in the Chandra observations between 0.1 and 1 Hz. We discuss our results in the context of current models for QPO formation.
  • We present $\emph{NuSTAR}$ observations of neutron star (NS) low-mass X-ray binaries: 4U 1636-53, GX 17+2, and 4U 1705-44. We observed 4U 1636-53 in the hard state, with an Eddington fraction, $F_{\mathrm{Edd}}$, of 0.01; GX 17+2 and 4U 1705-44 were in the soft state with fractions of 0.57 and 0.10, respectively. Each spectrum shows evidence for a relativistically broadened Fe K$_{\alpha}$ line. Through accretion disk reflection modeling, we constrain the radius of the inner disk in 4U 1636-53 to be $R_{in}=1.03\pm0.03$ ISCO (innermost stable circular orbit) assuming a dimensionless spin parameter $a_{*}=cJ/GM^{2}=0.0$, and $R_{in}=1.08\pm0.06$ ISCO for $a_{*}=0.3$ (errors quoted at 1 $\sigma$). This value proves to be model independent. For $a_{*}=0.3$ and $M=1.4\ M_{\odot}$, for example, $1.08\pm0.06$ ISCO translates to a physical radius of $R=10.8\pm0.6$ km, and the neutron star would have to be smaller than this radius (other outcomes are possible for allowed spin parameters and masses). For GX 17+2, $R_{in}=1.00-1.04$ ISCO for $a_{*}=0.0$ and $R_{in}=1.03-1.30$ ISCO for $a_{*}=0.3$. For $a_{*}=0.3$ and $M=1.4\ M_{\odot}$, $R_{in}=1.03-1.30$ ISCO translates to $R=10.3-13.0$ km. The inner accretion disk in 4U 1705-44 may be truncated just above the stellar surface, perhaps by a boundary layer or magnetosphere; reflection models give a radius of 1.46-1.64 ISCO for $a_{*}=0.0$ and 1.69-1.93 ISCO for $a_{*}=0.3$. We discuss the implications that our results may have on the equation of state of ultradense, cold matter and our understanding of the innermost accretion flow onto neutron stars with low surface magnetic fields, and systematic errors related to the reflection models and spacetime metric around less idealized neutron stars.
  • In the frame of the Astronomical Validation activities for the 64m Sardinia Radio Telescope, we performed 5-22 GHz imaging observations of the complex-morphology supernova remnants (SNRs) W44 and IC443. We adopted innovative observing and mapping techniques providing unprecedented accuracy for single-dish imaging of SNRs at these frequencies, revealing morphological details typically available only at lower frequencies through interferometry observations. High-frequency studies of SNRs in the radio range are useful to better characterize the spatially-resolved spectra and the physical parameters of different regions of the SNRs interacting with the ISM. Furthermore, synchrotron-emitting electrons in the high-frequency radio band are also responsible for the observed high-energy phenomenology as -e.g.- Inverse Compton and bremsstrahlung emission components observed in gamma-rays, to be disentangled from hadron emission contribution (providing constraints on the origin of cosmic rays).
  • We report the detection of a $78.1\pm0.5$ day period in the X-ray lightcurve of the extreme ultraluminous X-ray source NGC 5907 ULX1 ($L_{\rm{X,peak}}\sim5\times10^{40}$ erg s$^{-1}$), discovered during an extensive monitoring program with Swift. These periodic variations are strong, with the observed flux changing by a factor of $\sim$3-4 between the peaks and the troughs of the cycle; our simulations suggest that the observed periodicity is detected comfortably in excess of 3$\sigma$ significance. We discuss possible origins for this X-ray period, but conclude that at the current time we cannot robustly distinguish between orbital and super-orbital variations.
  • We present new observations of the galaxy cluster 3C 129 obtained with the Sardinia Radio Telescope in the frequency range 6000-7200 MHz, with the aim to image the large-angular-scale emission at high-frequency of the radio sources located in this cluster of galaxies. The data were acquired using the recently-commissioned ROACH2-based backend to produce full-Stokes image cubes of an area of 1 deg x 1 deg centered on the radio source 3C 129. We modeled and deconvolved the telescope beam pattern from the data. We also measured the instrumental polarization beam patterns to correct the polarization images for off-axis instrumental polarization. Total intensity images at an angular resolution of 2.9 arcmin were obtained for the tailed radio galaxy 3C 129 and for 13 more sources in the field, including 3C 129.1 at the galaxy cluster center. These data were used, in combination with literature data at lower frequencies, to derive the variation of the synchrotron spectrum of 3C 129 along the tail of the radio source. If the magnetic field is at the equipartition value, we showed that the lifetimes of radiating electrons result in a radiative age for 3C 129 of t_syn = 267 +/- 26 Myrs. Assuming a linear projected length of 488 kpc for the tail, we deduced that 3C 129 is moving supersonically with a Mach number of M=v_gal/c_s=1.47. Linearly polarized emission was clearly detected for both 3C 129 and 3C 129.1. The linear polarization measured for 3C 129 reaches levels as high as 70% in the faintest region of the source where the magnetic field is aligned with the direction of the tail.
  • We present the detection of an absorpton feature at $E=8.77^{+0.05}_{-0.06}$ keV in the combined X-ray spectrum of the ultraluminous X-ray source NGC 1313 X-1 observed with XMM-Newton and NuSTAR, significant at the 3$\sigma$ level. If associated with blueshifted ionized iron, the implied outflow velocity is ~0.2$c$ for Fe XXVI, or ~0.25$c$ for Fe XXV. These velocities are similar to the ultrafast outflow seen in absorption recently discovered in this source at lower energies by XMM-Newton, and we therefore conclude that this is an iron component to the same outflow. Photoionization modeling marginally prefers the Fe XXV solution, but in either case the outflow properties appear to be extreme, potentially supporting a super-Eddington hypothesis for NGC 1313 X-1.
  • We report on observations of the neutron star (NS) residing in the low-mass X-ray binary 1RXS J180408.9-34205 taken 2015 March by $\emph{NuSTAR}$ and $\emph{XMM-Newton}$ while the source was in the hard spectral state. We find multiple reflection features (Fe K$_{\alpha}$ detected with $\emph{NuSTAR}$; N VII, O VII, and O VIII detected in the RGS) from different ionization zones. Through joint fits using the self consistent relativistic reflection model {\sc relxill}, we determine the inner radius to be $\leq 11.1\ R_{g}$. For a 1.4 M$_{\odot}$ NS with a spin of $a_{*}=0$, this is an inner disk radius of $\leq22.2$ km. We find the inclination of the system to be between $18^{\circ}$-$29^{\circ}$. If the disk is truncated at a radius greater than the neutron star radius, it could be truncated by a boundary layer on the neutron star surface. It is also possible that the disk is truncated at the magnetospheric radius; conservative estimates would then imply $B\leq(0.3 -1.0)\times10^{9}$ G at the magnetic poles, though coherent pulsations have not been detected and the source is not identified as a pulsar.
  • We present the results from coordinated X-ray observations of the ultraluminous X-ray source NGC 5204 X-1 performed by NuSTAR and XMM-Newton in early 2013. These observations provide the first detection of NGC 5204 X-1 above 10 keV, extending the broadband coverage to 0.3-20 keV. The observations were carried out in two epochs separated by approximately 10 days, and showed little spectral variation, with an observed luminosity of Lx = (4.95+/-0.11)e39 erg/s. The broadband spectrum confirms the presence of a clear spectral downturn above 10 keV seen in some previous observations. This cutoff is inconsistent with the standard low/hard state seen in Galactic black hole binaries, as would be expected from an intermediate mass black hole accreting at significantly sub-Eddington rates given the observed luminosity. The continuum is apparently dominated by two optically thick thermal-like components, potentially accompanied by a faint high energy tail. The broadband spectrum is likely associated with an accretion disk that differs from a standard Shakura & Sunyaev thin disk.
  • 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.
  • We present spectral analysis of NuSTAR and Swift observations of Cep X-4 during its outburst in 2014. We observed the source once during the peak of the outburst and once during the decay, finding good agreement in the spectral shape between the observations. We describe the continuum using a powerlaw with a Fermi-Dirac cutoff at high energies. Cep X-4 has a very strong cyclotron resonant scattering feature (CRSF) around 30 keV. A simple absorption-like line with a Gaussian optical depth or a pseudo-Lorentzian profile both fail to describe the shape of the CRSF accurately, leaving significant deviations at the red side of the line. We characterize this asymmetry with a second absorption feature around 19 keV. The line energy of the CRSF, which is not influenced by the addition of this feature, shows a small but significant positive luminosity dependence. With luminosities between (1-6)e36 erg/s, Cep X-4 is below the theoretical limit where such a correlation is expected. This behavior is similar to Vela X-1 and we discuss parallels between the two systems.
  • We present the first broadband 0.3-25.0 kev X-ray observations of the bright ultraluminous X-ray source (ULX) Holmberg II X-1, performed by NuSTAR, XMM-Newton and Suzaku in September 2013. The NuSTAR data provide the first observations of Holmberg II X-1 above 10 keV, and reveal a very steep high-energy spectrum, similar to other ULXs observed by NuSTAR to date. These observations further demonstrate that ULXs exhibit spectral states that are not typically seen in Galactic black hole binaries. Comparison with other sources implies that Holmberg II X-1 accretes at a high fraction of its Eddington accretion rate, and possibly exceeds it. The soft X-ray spectrum (E<10 keV) appears to be dominated by two blackbody-like emission components, the hotter of which may be associated with an accretion disk. However, all simple disk models under-predict the NuSTAR data above ~10 keV and require an additional emission component at the highest energies probed, implying the NuSTAR data does not fall away with a Wien spectrum. We investigate physical origins for such an additional high-energy emission component, and favor a scenario in which the excess arises from Compton scattering in a hot corona of electrons with some properties similar to the very-high state seen in Galactic binaries. The observed broadband 0.3-25.0 keV luminosity inferred from these epochs is Lx = (8.1+/-0.1)e39 erg/s, typical for Holmberg II X-1, with the majority of the flux (~90%) emitted below 10 keV.