We discovered 2.8 s pulsations in the X-ray emission of the ultraluminous X-ray source (ULX) M51 ULX-7 within the UNSEeN project, which was designed to hunt for new pulsating ULXs (PULXs) with XMM-Newton. The pulse shape is sinusoidal and large variations of its amplitude were observed even within single exposures (pulsed fraction from less than 5% to 20%). M51 ULX-7 is a variable source, generally observed at an X-ray luminosity between $10^{39}$ and $10^{40}$ erg s$^{-1}$, located in the outskirts of the spiral galaxy M51a at a distance of 8.6 Mpc. According to our analysis, the X-ray pulsar orbits in a 2-d binary with a projected semi-major axis $a_\mathrm{X} \sin i \simeq$ 28 lt-s. For a neutron star (NS) of 1.4 $M_{\odot}$, this implies a lower limit on the companion mass of 8 $M_{\odot}$, placing the system hosting M51 ULX-7 in the high-mass X-ray binary class. The barycentric pulse period decreased by $\simeq$0.4 ms in the 31 d spanned by our May -- June 2018 observations, corresponding to a spin-up rate $\dot{P} \simeq -1.5\times10^{-10}\text{s s}^{-1}$. In an archival 2005 XMM-Newton exposure, we measured a spin period of $\sim$3.3 s, indicating a secular spin-up of $\dot{P}_{\mathrm{sec}}\simeq -10^{-9}\text{ s s}^{-1}$, a value in the range of other known PULXs. Our findings suggest that the system consists of an OB giant and a moderately magnetic (dipole field component in the range $10^{12}$ G $\lesssim B_{\mathrm{dip}}\lesssim 10^{13}$G) accreting NS with weakly beamed emission ($1/12\lesssim b\lesssim1/4$).

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

Ultraluminous X-ray sources (ULXs) are a class of accreting compact objects with X-ray luminosities above 1e39 erg/s. The ULX population counts several hundreds objects but only a minor fraction is well studied. Here we present a detailed analysis of all ULXs hosted in the galaxy NGC 7456. It was observed in X-rays only once in the past (in 2005) by XMM-Newton, but the observation was short and strongly affected by high background. In 2018, we obtained a new, deeper (~90 ks) XMM-Newton observation that allowed us to perform a detailed characterization of the ULXs hosted in the galaxy. ULX-1 and ULX-2, the two brightest objects (Lx~(6-10)e39 erg/s), have spectra that can be described by a two-thermal component model as often found in ULXs. ULX-1 shows also one order of magnitude in flux variability on short-term timescales (hundreds to thousand ks). The other sources (ULX-3 and ULX-4) show flux changes of at least an order of magnitude, and these objects may be candidate transient ULXs although longer X-ray monitoring or further studies are required to ascribe them to the ULX population. In addition, we found a previously undetected source that might be a new candidate ULX (labelled as ULX-5) with a luminosity of ~1e39 erg/s and hard power-law spectral shape, whose nature is still unclear and for which a background Active Galactic Nucleus cannot be excluded. We discuss the properties of all the ULXs in NGC 7456 within the framework of super-Eddington accretion onto stellar mass compact objects. Although no pulsations were detected, we cannot exclude that the sources host neutron stars.

Ultraluminous X-ray Sources (ULXs) provide a unique opportunities to probe the geometry and energetics of super-Eddington accretion. The radiative processes involved in super-Eddington accretion are not well understood, and so studying correlated variability between different energy bands can provide insights into the causal connection between different emitting regions. We present a spectral-timing analysis of NGC 1313 X-1 from a recent XMM-Newton campaign. The spectra can be decomposed into two thermal-like components, the hotter of which may originate from the inner accretion disc, and the cooler from an optically thick outflow. We find correlated variability between hard (2-10 keV) and soft (0.3-2 keV) bands on kilosecond timescales, and find a soft lag of ~150 seconds. The covariance spectrum suggests that emission contributing to the lags is largely associated with the hotter of the two thermal-like components, likely originating from the inner accretion flow. This is only the third ULX to exhibit soft lags. The lags range over three orders of magnitude in amplitude, but all three are ~5 to ~20 percent of the corresponding characteristic variability timescales. If these soft lags can be understood in the context of a unified picture of ULXs, then lag timescales may provide constraints on the density and extent of radiatively-driven outflows.

M82 X-2 is the first pulsating ultraluminous X-ray source (PULX)to be identified. Since the discovery in 2014, NuSTAR has observed the M82 field 15 times throughout 2015 and 2016. In this paper, we report the results of pulsation searches in all these datasets, and find only one new detection. This new detection allows us to refine the orbital period of the source and measure an average spin down rate between 2014 and 2016 of 5x10^-11 Hz/s, which is in contrast to the strong spin up seen during the 2014 observations and represent the first detection of spin down in a PULX system. Thanks to the improved orbital solution allowed by this new detection, we are also able to detect pulsations in additional segments of the original 2014 dataset. We find a glitch superimposed on the very strong and variable spin-up already reported, the first positive glitch identified in a PULX system. We discuss the new findings in the context of current leading models for PULXs.

We report NuSTAR and Chandra observations of two X-ray transients, SWIFT J174540.7$-$290015 (T15) and SWIFT J174540.2$-$290037 (T37), which were discovered by the Neil Gehrels Swift Observatory in 2016 within $r\sim1$ pc of Sgr A*. NuSTAR detected bright X-ray outbursts from T15 and T37, likely in the soft and hard states, with 3-79~keV luminosities of $8\times10^{36}$ and $3\times10^{37}$ erg/s, respectively. No X-ray outbursts have previously been detected from the two transients and our Chandra ACIS analysis puts an upper limit of $L_X \lesssim 2 \times10^{31}$ erg/s on their quiescent 2-8 keV luminosities. No pulsations, significant QPOs, or type I X-ray bursts were detected in the NuSTAR data. While T15 exhibited no significant red noise, the T37 power density spectra are well characterized by three Lorentzian components. The declining variability of T37 above $\nu \sim 10$ Hz is typical of black hole (BH) transients in the hard state. NuSTAR spectra of both transients exhibit a thermal disk blackbody, X-ray reflection with broadened Fe atomic features, and a continuum component well described by Comptonization models. Their X-ray reflection spectra are most consistent with high BH spin ($a_{*} \gtrsim 0.9$) and large disk density ($n_e\sim10^{21}$ cm$^{-3}$). Based on the best-fit ionization parameters and disk densities, we found that X-ray reflection occurred near the inner disk radius, which was derived from the relativistic broadening and thermal disk component. These X-ray characteristics suggest the outbursting BH-LMXB scenario for both transients and yield the first BH spin measurements from X-ray transients in the central 100 pc region.

This paper describes the design and implementation of Stingray, a library in Python built to perform time series analysis and related tasks on astronomical light curves. Its core functionality comprises a range of Fourier analysis techniques commonly used in spectral-timing analysis, as well as extensions for analyzing pulsar data, simulating data sets, and statistical modeling. Its modular build allows for easy extensions and incorporation of its methods into data analysis workflows and pipelines. We aim for the library to be a platform for the implementation of future spectral-timing techniques. Here, we describe the overall vision and framework, core functionality, extensions, and connections to high-level command-line and graphical interfaces. The code is well-tested, with a test coverage of currently 95%, and is accompanied by extensive API documentation and a set of step-by-step tutorials.

The X-ray Polarization Probe (XPP) is a second generation X-ray polarimeter following up on the Imaging X-ray Polarimetry Explorer (IXPE). The XPP will offer true broadband polarimetery over the wide 0.2-60 keV bandpass in addition to imaging polarimetry from 2-8 keV. The extended energy bandpass and improvements in sensitivity will enable the simultaneous measurement of the polarization of several emission components. These measurements will give qualitatively new information about how compact objects work, and will probe fundamental physics, i.e. strong-field quantum electrodynamics and strong gravity.

We report the detection of weak pulsations from the archetypal ultraluminous X-ray source (ULX) NGC 1313 X-2. Acceleration searches reveal sinusoidal pulsations in segments of two out of six new deep observations of this object, with a period of $\sim$ 1.5 s and a pulsed fraction of $\sim$ 5%. We use Monte Carlo simulations to demonstrate that the individual detections are unlikely to originate in false Poisson noise detections given their very close frequencies; their strong similarity to other pulsations detected from ULXs also argues they are real. The presence of a large bubble nebula surrounding NGC 1313 X-2 implies an age of order 1 Myr for the accreting phase of the ULX, which implies that the neutron star's magnetic field has not been suppressed over time by accreted material, nor has the neutron star collapsed into a black hole, despite an average energy output into the nebula two orders of magnitude above Eddington. This argues that most of the accreted material has been expelled over the lifetime of the ULX, favouring physical models including strong winds and/or jets for neutron star ULXs.

SS433 is the only Galactic binary system known to accrete at highly super-critical rates, analogous to tidal disruption events, and needed to explain the mass of some high redshift quasars. Probing the inner regions of SS433 in the X-rays is crucial to understanding this system, and super-critical accretion in general, but has not yet been possible due to obscuration. NuSTAR observed SS433 in the hard X-ray band across multiple phases of its super-orbital precession period. Spectral-timing tools have allowed us to confirm that the hard X-ray emission from the inner regions is scattered towards us by the walls of the wind-cone. By comparing to numerical models, we determine an intrinsic X-ray luminosity of >= 3x10^37 erg/s and that, if viewed face on, the apparent luminosity would be > 1x10^39 erg/s, confirming its long-suspected nature as an ultraluminous X-ray source (ULX). A lag due to absorption by Fe XXV/XXVI in outflowing material travelling at least 0.14-0.29c matches absorption lines seen in ULXs and - in future - will allow us to map a super-critical outflow for the first time.

Pulsar studies in the recent years have shown, more than others, to have benefited from a multi-wavelength approach. The INAF - Astronomical Observatory in Cagliari (INAF-OAC) is a growing facility with a young group devoted to pulsar and fast transients studies across the electromagnetic spectrum. Taking advantage of this expertise we have worked to provide a suite of multi-wavelength software and databases for the observations of pulsars and compact Galactic objects at the Sardinia Radio Telescope (SRT). In turn, radio pulsar observations at SRT will be made available, in a processed format, to gamma-ray searches using AGILE and Fermi gamma-ray satellite and, in a near future, they will be complementary to polarimetric X-ray observations with IXPE.

Over the last decades X-ray spectroscopy has proven to be a powerful tool for the estimation of black hole spin and several other key parameters in dozens of AGN and black hole X-ray binaries. In this White Paper, we discuss the observational and theoretical challenges expected in the exploration, discovery, and study of astrophysical black holes in the next decade. We focus on the case of accreting black holes and their electromagnetic signatures, with particular emphasis on the measurement of the relativistic reflection component in their X-ray spectra.

With the recent discoveries of massive and highly luminous quasars at high redshifts ($z\sim7$; e.g. Mortlock et al. 2011), the question of how black holes (BHs) grow in the early Universe has been cast in a new light. In order to grow BHs with $M_{\rm BH} > 10^9$ M$_{\odot}$ by less than a billion years after the Big Bang, mass accretion onto the low-mass seed BHs needs to have been very rapid (Volonteri & Rees, 2005). Indeed, for any stellar remnant seed, the rate required would need to exceed the Eddington limit. This is the point at which the outward force produced by radiation pressure is equal to the gravitational attraction experienced by the in-falling matter. In principle, this implies that there is a maximum luminosity an object of mass $M$ can emit; assuming spherical accretion and that the opacity is dominated by Thompson scattering, this Eddington luminosity is $L_{\rm{E}} = 1.38 \times 10^{38} (M/M_{\odot})$ erg s$^{-1}$. In reality, it is known that this limit can be violated, due to non-spherical geometry or various kinds of instabilities. Nevertheless, the Eddington limit remains an important reference point, and many of the details of how accretion proceeds above this limit remain unclear. Understanding how this so-called super-Eddington accretion occurs is of clear cosmological importance, since it potentially governs the growth of the first supermassive black holes (SMBHs) and the impact this growth would have had on their host galaxies (`feedback') and the epoch of reionization, as well as improving our understanding of accretion physics more generally.

We report on NuSTAR observations of transient pulsations in the neutron star X-ray binary SMC X-1. The transition from non-pulsing to pulsing states was not accompanied by a large change in flux. Instead, both pulsing and non-pulsing states were observed in a single observation during the low-flux super-orbital state. During the high-state, we measure a pulse period of $P = 0.70117(9)\,\mathrm{s}$ at $T_{ref} = 56145\,\mathrm{MJD}$. Spectral analysis during non-pulsing and pulsing states reveals that the observations can be consistently modeled by an absorbed power law with a phenomenological cutoff resembling a Fermi-Dirac distribution, or by a partially obscured cutoff power law. The shapes of the underlying continua show little variability between epochs, while the covering fraction and column density vary between super-orbital states. The strength of pulsations also varies, leading us to infer that the absence and reemergence of pulsations are related to changing obscuration, such as by a warped accretion disk. SMC X-1 is accreting near or above its Eddington limit, reaching an unabsorbed X-ray luminosity of $L_{\rm X}({\rm 2-10~keV}) \approx 5 \times 10^{38}\, {\rm erg}\, {\rm s}^{-1}$. This suggests that SMC X-1 may be a useful local analog to ultraluminous X-ray pulsars (ULXPs), which likewise exhibit strong variability in their pulsed fractions, as well as flux variability on similar timescales. In particular, the gradual pulse turn-on which has been observed in M82 X-2 is similar to the behavior we observe in SMC X-1. Thus we propose that pulse fraction variability of ULXPs may also be due to variable obscuration.

The reported discovery of a cyclotron resonance scattering feature (CRSF) in the spectrum of M51 ULX-8 may provide an important clue as to the nature of the magnetic field in those ultraluminous X-ray sources hosting neutron stars. In this paper we present the covariance (linearly correlated variability) spectrum of M51 ULX-8 on long (> 2000s) timescales. This allows us to unambiguously decompose the spectrum which requires multiple components in order to explain the broad-band emission. Having a clearer picture of the spectral decomposition leads to various tests for the dipole field strength of the neutron star which can be extended to other ULXs when certain criteria are met. In the case of M51 ULX-8, we rule out a very strong (~10^15 G) dipole solution with either a sub- or super-critical disc. Instead, our tests indicate an upper limit on the dipole field of ~10^12 G, and a classical super-critical inflow, similar to that inferred in other ULXs found to harbour neutron stars, although we do not rule out the presence of an additional, strong (~10^15 G) multipole field falling off steeply with distance from the neutron star.

We present the Spectroscopic Time-Resolving Observatory for Broadband Energy X-rays (STROBE-X), a probe-class mission concept selected for study by NASA. It combines huge collecting area, high throughput, broad energy coverage, and excellent spectral and temporal resolution in a single facility. STROBE-X offers an enormous increase in sensitivity for X-ray spectral timing, extending these techniques to extragalactic targets for the first time. It is also an agile mission capable of rapid response to transient events, making it an essential X-ray partner facility in the era of time-domain, multi-wavelength, and multi-messenger astronomy. Optimized for study of the most extreme conditions found in the Universe, its key science objectives include: (1) Robustly measuring mass and spin and mapping inner accretion flows across the black hole mass spectrum, from compact stars to intermediate-mass objects to active galactic nuclei. (2) Mapping out the full mass-radius relation of neutron stars using an ensemble of nearly two dozen rotation-powered pulsars and accreting neutron stars, and hence measuring the equation of state for ultradense matter over a much wider range of densities than explored by NICER. (3) Identifying and studying X-ray counterparts (in the post-Swift era) for multiwavelength and multi-messenger transients in the dynamic sky through cross-correlation with gravitational wave interferometers, neutrino observatories, and high-cadence time-domain surveys in other electromagnetic bands. (4) Continuously surveying the dynamic X-ray sky with a large duty cycle and high time resolution to characterize the behavior of X-ray sources over an unprecedentedly vast range of time scales. STROBE-X's formidable capabilities will also enable a broad portfolio of additional science.

Ultraluminous X-ray (ULX) pulsars are a new class of object powered by apparent super-critical accretion onto magnetized neutron stars. Three sources in this class identified so far; M82 X-2, NGC 5907 ULX-1 and NGC 7793 P13, have been found to have two properties in common; $\sim1$-s spin periods, and for NGC 5907 ULX-1 and NGC 7793 P13 periodic X-ray flux modulations on timescales of $\sim60-80$ days. M82 X-2 resides in a crowded field that includes the ULX M82 X-1 separated from X-2 by 5", as well as other bright point sources. A 60-day modulation has been observed from the region but the origin has been difficult to identify; both M82 X-1 and X-2 have been suggested as the source. In this paper we present the analysis of a systematic monitoring campaign by Chandra, the only X-ray telescope capable of resolving the crowded field. From a simple Lomb-Scargle periodogram analysis and a more sophisticated Gaussian Process analysis we find that only X-2 exhibits a periodic signal around 60 days supporting previous claims that it is the origin. We also construct a phase-averaged flux profile of the modulations from higher cadence Swift/XRT data and find that the flux variations in the Chandra data are fully consistent with the flux profile. Since the orbit of the neutron star and its companion is known to be 2.5 days, the $\sim60$-day period must be super-orbital in origin. The flux of the modulations varies by a factor of $\sim$100 from minimum to maximum, with no evidence for spectral variations, making the origin difficult to explain.

After two major outbursts in 2006 and 2011, on 2017 May 16 the magnetar CXOU J164710.2-455216, hosted within the massive star cluster Westerlund I, emitted a short (20 ms) burst, which marked the onset of a new active phase. We started a long-term monitoring campaign with Swift (45 observations), Chandra (5 observations) and NuSTAR (4 observations) from the activation until 2018 April. During the campaign, Swift BAT registered the occurrence of multiple bursts, accompanied by two other enhancements of the X-ray persistent flux. The long time span covered by our observations allowed us to study the spectral as well as the timing evolution of the source. After 11 months since the 2017 May outburst onset, the observed flux was about 15 times higher than its historical minimum level and a factor of 3 higher than the level reached after the 2006 outburst. This suggests that the crust has not fully relaxed to the quiescent level, or that the source quiescent level has changed following the multiple outburst activities in the past 10 years or so. This is another case of multiple outbursts from the same source on a yearly time scale, a somehow recently discovered behaviour in magnetars.

Two \textit{NuSTAR} observations of the luminous X-ray pulsar LMC X-4 in October and November 2015 captured several bright accretion flares from this source, which has a long history of stable pulse and superorbital behavior. We present a timing analysis of these data in which we detect a rapid pulse turn-on in association with the accretion flares, during which the source reaches super-Eddington luminosities. Pulsations, which are normally seen from this source, are found to only occur for approximately one hour before and during the bright flares. Beyond one hour before and after the flares, we find pulsations to be weak or nonexistent, with fractional RMS amplitudes of less than 0.05. At the onset of the flare, the pulse profiles exhibit a phase shift of 0.25 cycles that could be associated with a change in the emission geometry. This increase in pulse strength occurring well before the flare cannot be explained by the propeller effect, and potentially offers a connection between the magnetic properties of pulsars that accrete close to their Eddington limits and ULX pulsars.

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 more than 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 65000 km/s (0.22c) and a high ionisation factor with a log value of 3.9. This is the first direct evidence for a UFO in a neutron star ULX and also the first time that this its evidence in a ULX spectrum is seen in both soft and hard X-ray data simultaneously. 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 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 X-ray sources (ULXs) are bright X-ray sources in nearby galaxies not associated with the central supermassive black hole. Their luminosities imply they are powered by either an extreme accretion rate onto a compact stellar remnant, or an intermediate mass ($\sim100-10^5$ M$_{\odot}$) black hole. The recent detection of coherent pulsations coming from three bright ULXs demonstrates that some of these sources are powered by accretion onto a neutron star, implying accretion rates significantly in excess of the Eddington limit, a high degree of geometric beaming, or both. The physical challenges associated with the high implied accretion rates can be mitigated if the neutron star surface field is high - in the magnetar regime ($10^{14}$ G), since this suppresses the electron scattering cross section, reducing the radiation pressure that chokes off accretion for high luminosities. One of the few ways to determine surface magnetic fields is through the detection of cyclotron resonance scattering features (CRSFs) produced by the transition of charged particles between quantized Landau levels. To date, CRSFs have only been detected in Galactic accreting pulsars. Here we present the detection at 3.8-{\sigma} significance of a strong absorption line at a rest-frame energy of 4.5 keV in the Chandra spectrum of a ULX in M51. We find that this feature is likely to be a CRSF produced by the strong magnetic field of a neutron star. Assuming scattering off electrons, the magnetic field strength is implied to be $\sim10^{11}$ G, however the line is narrower than any electron CRSFs previously observed, and assuming thermal broadening, the implied temperature is significantly cooler than the accretion disk or column. The line shape is, however, consistent with a proton resonance scattering feature, implying that the neutron star has a magnetic field near the surface of B$\sim10^{15}$ G.

Dead time affects many of the instruments used in X-ray astronomy, by producing a strong distortion in power density spectra. This can make it difficult to model the aperiodic variability of the source or look for quasi-periodic oscillations. Whereas in some instruments a simple a priori correction for dead-time-affected power spectra is possible, this is not the case for others such as NuSTAR, where the dead time is non-constant and long (~2.5 ms). Bachetti et al. 2015 suggested the cospectrum obtained from light curves of independent detectors within the same instrument as a possible way out, but this solution has always only been a partial one: the measured rms was still affected by dead time, because the width of the power distribution of the cospectrum was modulated by dead time in a frequency-dependent way. In this Letter we suggest a new, powerful method to normalize cospectra and, with some caveats, even power density spectra. Our approach uses the difference of the Fourier amplitudes from two independent detectors to characterize and filter out the effect of dead time. This method is crucially important for the accurate modelling of periodograms derived from instruments affected by dead time on board current missions like NuSTAR and ASTROSAT, but also future missions such as IXPE

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.

We present a spectral and timing study of the NuSTAR and Swift observations of the black hole candidate IGR J17091-3624 in the hard state during its outburst in 2016. Disk reflection is detected in each of the NuSTAR spectra taken in three epochs. Fitting with relativistic reflection models reveals that the accretion disk is truncated during all epochs with $R_{\rm in}>10~r_{\rm g}$, with the data favoring a low disk inclination of $\sim 30^{\circ}-40^{\circ}$. The steepening of the continuum spectra between epochs is accompanied by a decrease in the high energy cut-off: the electron temperature $kT_{\rm e}$ drops from $\sim 64$ keV to $\sim 26$ keV, changing systematically with the source flux. We detect type-C QPOs in the power spectra with frequency varying between 0.131 Hz and 0.327 Hz. In addition, a secondary peak is found in the power spectra centered at about 2.3 times the QPO frequency during all three epochs. The nature of this secondary frequency is uncertain, however a non-harmonic origin is favored. We investigate the evolution of the timing and spectral properties during the rising phase of the outburst and discuss their physical implications.

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 limit of $B \lesssim 6 \times 10^{12}$ G for the dipolar component of the central neutron star's magnetic field. 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.

Relativistic plasma jets are observed in many accreting black holes. According to theory, coiled magnetic fields close to the black hole accelerate and collimate the plasma, leading to a jet being launched. Isolating emission from this acceleration and collimation zone is key to measuring its size and understanding jet formation physics. But this is challenging because emission from the jet base cannot be easily disentangled from other accreting components. Here, we show that rapid optical flux variations from a Galactic black-hole binary are delayed with respect to X-rays radiated from close to the black hole by ~0.1 seconds, and that this delayed signal appears together with a brightening radio jet. The origin of these sub-second optical variations has hitherto been controversial. Not only does our work strongly support a jet origin for the optical variations, it also sets a characteristic elevation of <~10$^3$ Schwarzschild radii for the main inner optical emission zone above the black hole, constraining both internal shock and magnetohydrodynamic models. Similarities with blazars suggest that jet structure and launching physics could potentially be unified under mass-invariant models. Two of the best-studied jetted black hole binaries show very similar optical lags, so this size scale may be a defining feature of such systems.

Aims: In light of recent discoveries of pulsating ULXs and recently introduced models, that propose neutron stars (NSs) as the central engines of ULXs, we revisit the (XMM-Newton and NuSTAR) spectra of eighteen well known ULXs, in search of indications that favor this new hypothesis. Results: We confirm the, previously noted, presence of the low energy (<6 keV) spectral rollover and argue that it could be interpreted as thermal emission. The spectra are well described by a double thermal model consisting of a "hot" (>1 keV) and a "cool" (<0.7 keV) multicolor black body (MCB). Under the assumption that the "cool" MCD emission originates in a disk truncated at the NS magnetosphere, we find that all ULXs in our sample are consistent with accretion onto a highly magnetized (B>10^12G) NS. We note a strong correlation between the strength of the magnetic field, the temperature of the "hot" thermal component and the total unabsorbed luminosity. Examination of the NuSTAR data supported this interpretation and confirmed the presence of a weak, high energy (>15 keV) tail, most likely the result of modification of the MCB emission by inverse Compton scattering. We also note that the apparent high energy tail, may simply be the result of mis-modeling of MCB emission with an atypical temperature (T) vs radius (r) gradient, using a standard MCD model with a fixed gradient of T~r^-0.75. Conclusions: We have offered a new physical interpretation for the dual-thermal spectra of ULXs. We find that the best-fit derived parameters, are in excellent agreement with recent predictions that favor super-critically accreting NSs as the engines of a large fraction of ULXs. Nevertheless, the considerable degeneracy between models and the lack of unequivocal evidence cannot rule out other equally plausible interpretations. Our results warrant further exploration of this new framework.

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.

IGR J18214-1318, a Galactic source discovered by the International Gamma-Ray Astrophysics Laboratory, is a high-mass X-ray binary (HMXB) with a supergiant O-type stellar donor. We report on the XMM-Newton and NuSTAR observations that were undertaken to determine the nature of the compact object in this system. This source exhibits high levels of aperiodic variability, but no periodic pulsations are detected with a 90% confidence upper limit of 2% fractional rms between 0.00003-88 Hz, a frequency range that includes the typical pulse periods of neutron stars (NSs) in HMXBs (0.1-10$^3$ s). Although the lack of pulsations prevents us from definitively identifying the compact object in IGR J18214-1318, the presence of an exponential cutoff with e-folding energy $\lesssim30$ keV in its 0.3-79 keV spectrum strongly suggests that the compact object is an NS. The X-ray spectrum also shows a Fe K$\alpha$ emission line and a soft excess, which can be accounted for by either a partial-covering absorber with $N_{\mathrm{H}}\approx10^{23}$ cm$^{-2}$ which could be due to the inhomogeneous supergiant wind, or a blackbody component with $kT=1.74^{+0.04}_{-0.05}$ keV and $R_{BB}\approx0.3$ km, which may originate from NS hot spots. Although neither explanation for the soft excess can be excluded, the former is more consistent with the properties observed in other supergiant HMXBs. We compare IGR J18214-1318 to other HMXBs that lack pulsations or have long pulsation periods beyond the range covered by our observations.

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.

Neutron stars are thought to be born rapidly rotating and then exhibit a phase of a rotation-powered pulsations as they slow down to 1-10 s periods. The significant population of millisecond pulsars observed in our Galaxy is explained by the recycling concept: during an epoch of accretion from a donor star in a binary system, the neutron star is spun up to millisecond periods. However, only a few pulsars are observed during this recycling process, with relatively high rotational frequencies. Here we report the detection of an X-ray pulsar with $P_{\rm spin} = 1.20$ s in the globular cluster B091D in the Andromeda galaxy, the slowest pulsar ever found in a globular cluster. This bright (up-to 30% of the Eddington luminosity) spinning-up pulsar, persistent over the 12 years of observations, must have started accreting less than 1 Myr ago and has not yet had time to accelerate to hundreds of Hz. The neutron star in this unique wide binary with an orbital period $P_{\rm orb} = 30.5$ h in a 12 Gyr old, metal rich star cluster, accretes from a low mass, slightly evolved post-main sequence companion. We argue that we are witnessing a binary formed at relatively recent epoch by getting a $\sim 0.8M_\odot$ star in a dynamical interaction -- a viable scenario in a massive dense globular cluster like B091D with high global and specific stellar encounter rates. This intensively accreting non-recycled X-ray pulsar provides therefore a long-sought missing piece in the standard pulsar recycling picture.

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).

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.

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.

We present broad-band, multi-epoch X-ray spectroscopy of the pulsating ultra-luminous X-ray source (ULX) in NGC 5907. Simultaneous XMM-Newton and NuSTAR data from 2014 are best described by a multi-color black-body model with a temperature gradient as a function of accretion disk radius significantly flatter than expected for a standard thin accretion disk (T(r) ~ r^{-p}, with p=0.608^{+0.014}_{-0.012}). Additionally, we detect a hard power-law tail at energies above 10 keV, which we interpret as being due to Comptonization. We compare this observation to archival XMM-Newton, Chandra, and NuSTAR data from 2003, 2012, and 2013, and investigate possible spectral changes as a function of phase over the 78d super-orbital period of this source. We find that observations taken around phases 0.3-0.4 show very similar temperature profiles, even though the observed flux varies significantly, while one observation taken around phase 0 has a significantly steeper profile. We discuss these findings in light of the recent discovery that the compact object is a neutron star and show that precession of the accretion disk or the neutron star can self-consistently explain most observed phenomena.

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.

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).

Supernova Remnants (SNRs) exhibit spectra featured by synchrotron radio emission arising from the relativistic electrons, and high-energy emission from both leptonic (Bremsstrahlung and Inverse Compton) and hadronic processes (${\pi}^0$ mesons decay) which are a direct signature of cosmic rays acceleration. Thanks to radio single-dish imaging observations obtained in three frequency bands (1.6, 7, 22 GHz) with the Sardinia Radio Telescope (www.srt.inaf.it), we can model different SNR regions separately. Indeed, in order to disentangle interesting and peculiar hadron contributions in the high-energy spectra (gamma-ray band) and better constrain SNRs as cosmic rays emitters, it is crucial to fully constrain lepton contributions first through radio-observed parameters. In particular, the Bremsstrahlung and Inverse Compton bumps observed in gamma-rays are bounded to synchrotron spectral slope and cut-off in the radio domain. Since these parameters vary for different SNR regions and electron populations, spatially-resolved radio spectra are then required for accurate multi-wavelength modelling.

The 6.67 hr periodicity and the variable X-ray flux of the central compact object (CCO) at the center of the SNR RCW 103, named 1E 161348-5055, have been always difficult to interpret within the standard scenarios of an isolated neutron star or a binary system. On 2016 June 22, the Burst Alert Telescope (BAT) onboard Swift detected a magnetar-like short X-ray burst from the direction of 1E 161348-5055, also coincident with a large long-term X-ray outburst. Here we report on Chandra, NuSTAR, and Swift (BAT and XRT) observations of this peculiar source during its 2016 outburst peak. In particular, we study the properties of this magnetar-like burst, we discover a hard X-ray tail in the CCO spectrum during outburst, and we study its long-term outburst history (from 1999 to July 2016). We find the emission properties of 1E 161348-5055 consistent with it being a magnetar. However in this scenario, the 6.67 hr periodicity can only be interpreted as the rotation period of this strongly magnetized neutron star, which therefore represents the slowest pulsar ever detected, by orders of magnitude. We briefly discuss the viable slow-down scenarios, favoring a picture involving a period of fall-back accretion after the supernova explosion, similarly to what is invoked (although in a different regime) to explain the "anti-magnetar" scenario for other CCOs.

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

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

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.

The recent discovery by Bachetti et al. (2014) of a pulsar in M82 that can reach luminosities of up to 10^40 ergs s^-1, a factor of ~100 the Eddington luminosity for a 1.4 Msol compact object, poses a challenge for accretion physics. In order to better understand the nature of this source and its duty cycle, and in the light of several physical models that have been subsequently published, we conduct a spectral and temporal analysis of the 0.5-8 keV X-ray emission from this source from 15 years of Chandra observations. We fit the Chandra spectra of the pulsar with a power-law model and a disk black body model, subjected to interstellar absorption in M82. We carefully assess for the effect of pile-up in our observations, where 4/19 observations have a pile-up fraction >10%, which we account for during spectral modeling with a convolution model. When fitted with a power-law model, the average photon index when the source is at high luminosity (L_X>10^39 ergs s^-1) is Gamma=1.33+/-0.15. For the disk black body model, the average temperature is T=3.24+/-0.65 keV, consistent with other luminous X-ray pulsars. We also investigated the inclusion of a soft excess component and spectral break, finding that the spectra are also consistent with these features common to luminous X-ray pulsars. In addition, we present spectral analysis from NuSTAR over the 3-50 keV range where we have isolated the pulsed component. We find that the pulsed emission in this band is best fit by a power-law with a high-energy cut-off, where Gamma=0.6+/-0.3 and E_C=14^{+5}_{-3} keV. While the pulsar has previously been identified as a transient, we find from our longer-baseline study that it has been remarkably active over the 15-year period, where for 9/19 (47%) observations that we analyzed, the pulsar appears to be emitting at a luminosity in excess of 10^39 ergs s^-1, greater than 10 times its Eddington limit.

Multi-messenger astronomy is becoming the key to understanding the Universe from a comprehensive perspective. In most cases, the data and the technology are already in place, therefore it is important to provide an easily-accessible package that combines datasets from multiple telescopes at different wavelengths. In order to achieve this, we are working to produce a data analysis pipeline that allows the data reduction from different instruments without needing detailed knowledge of each observation. Ideally, the specifics of each observation are automatically dealt with, while the necessary information on how to handle the data in each case is provided by a tutorial that is included in the program. We first focus our project on the study of pulsars and their wind nebulae (PWNe) at radio and gamma-ray frequencies. In this way, we aim to combine time-domain and imaging datasets at two extremes of the electromagnetic spectrum. In addition, the emission has the same non-thermal origin in pulsars at radio and gamma-ray frequencies, and the population of electrons is believed to be the same at these energies in PWNe. The final goal of the project will be to unveil the properties of these objects by tracking their behaviour using all of the available multi-wavelength data.

We present a new method to identify luminous off-nuclear X-ray sources in the outskirts of galaxies from large public redshift surveys, distinguishing them from foreground and background interlopers. Using the 3XMM-DR5 catalog of X-ray sources and the SDSS DR12 spectroscopic sample of galaxies, with the help of this off-nuclear cross-matching technique, we selected 98 sources with inferred X-ray luminosities in the range $10^{41} < L_{\rm X} < 10^{44}\,{\rm erg\,s}^{-1}$, compatible with hyperluminous X-ray objects (HLX). To validate the method, we verify that it allowed us to recover known HLX candidates such as ESO 243$-$49 HLX$-$1 and M82 X$-$1. From a statistical study, we conservatively estimate that up to $71 \pm 11$ of these sources may be fore- or background sources, statistically leaving at least 16 that are likely to be HLXs, thus providing support for the existence of the HLX population. We identify two good HLX candidates and using other publicly available datasets, in particular the VLA FIRST in radio, UKIDSS in the near-infrared, GALEX in the ultra-violet and CFHT Megacam archive in the optical, we present evidence that these objects are unlikely to be foreground or background X-ray objects of conventional types, e.g. active galactic nuclei, BL Lac objects, Galactic X-ray binaries or nearby stars. However, additional dedicated X-ray and optical observations are needed to confirm their association with the assumed host galaxies and thus secure their HLX classification.

The extreme extragalactic sources known as Ultraluminous X-ray Sources (ULX) represent a unique testing environment for compact objects population studies and the accretion process. Their nature has long been disputed. Their luminosity, well above the Eddington luminosity for a stellar-mass black hole, can imply the presence of an intermediate-mass black hole or a stellar black hole accreting above the Eddington limit. Both these interpretations are important to understand better the accretion process and the evolution of massive black holes. The last few years have seen a dramatic improvement of our knowledge of these sources. In particular, the super-Eddington interpretation for the bulk of the ULX population has gained a strong consensus. Nonetheless, exceptions to this general trend do exist, and in particular one ULX was shown to be a neutron star, and another was shown to be a very likely IMBH candidate. In this paper, I will review the progress done in the last few years.

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

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.

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.

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

We present results for two Ultraluminous X-ray Sources (ULXs), IC 342 X-1 and IC 342 X-2, using two epochs of XMM-Newton and NuSTAR observations separated by $\sim$7 days. We observe little spectral or flux variability above 1 keV between epochs, with unabsorbed 0.3--30 keV luminosities being $1.04^{+0.08}_{-0.06} \times 10^{40}$ erg s$^{-1}$ for IC 342 X-1 and $7.40\pm0.20 \times 10^{39}$ erg s$^{-1}$ for IC 342 X-2, so that both were observed in a similar, luminous state. Both sources have a high absorbing column in excess of the Galactic value. Neither source has a spectrum consistent with a black hole binary in low/hard state, and both ULXs exhibit strong curvature in their broadband X-ray spectra. This curvature rules out models that invoke a simple reflection-dominated spectrum with a broadened iron line and no cutoff in the illuminating power-law continuum. X-ray spectrum of IC 342 X-1 can be characterized by a soft disk-like black body component at low energies and a cool, optically thick Comptonization continuum at high energies, but unique physical interpretation of the spectral components remains challenging. The broadband spectrum of IC 342 X-2 can be fit by either a hot (3.8 keV) accretion disk, or a Comptonized continuum with no indication of a seed photon population. Although the seed photon component may be masked by soft excess emission unlikely to be associated with the binary system, combined with the high absorption column, it is more plausible that the broadband X-ray emission arises from a simple thin blackbody disk component. Secure identification of the origin of the spectral components in these sources will likely require broadband spectral variability studies.

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}$.

We report on an observation of the Galactic black hole candidate GRS 1739-278 during its 2014 outburst, obtained with NuSTAR. The source was captured at the peak of a rising "low/hard" state, at a flux of ~0.3 Crab. A broad, skewed iron line and disk reflection spectrum are revealed. Fits to the sensitive NuSTAR spectra with a number of relativistically blurred disk reflection models yield strong geometrical constraints on the disk and hard X-ray "corona". Two models that explicitly assume a "lamppost" corona find its base to have a vertical height above the black hole of h = 5 (+7, -2) GM/c^2 and h = 18 +/-4 GM/c^2 (90% confidence errors); models that do not assume a "lamppost" return emissivity profiles that are broadly consistent with coronae of this size. Given that X-ray microlensing studies of quasars and reverberation lags in Seyferts find similarly compact coronae, observations may now signal that compact coronae are fundamental across the black hole mass scale. All of the models fit to GRS 1739-278 find that the accretion disk extends very close to the black hole - the least stringent constraint is r = 5 (+3,-4) GM/c^2. Only two of the models deliver meaningful spin constraints, but a = 0.8 +/-0.2 is consistent with all of the fits. Overall, the data provide especially compelling evidence of an association between compact hard X-ray coronae and the base of relativistic radio jets in black holes.

The low-mass X-ray binary Cen X-4 is the brightest and closest (<1.2 kpc) quiescent neutron star transient. Previous 0.5-10 keV X-ray observations of Cen X-4 in quiescence identified two spectral components: soft thermal emission from the neutron star atmosphere and a hard power-law tail of unknown origin. We report here on a simultaneous observation of Cen X-4 with NuSTAR (3-79 keV) and XMM-Newton (0.3-10 keV) in 2013 January, providing the first sensitive hard X-ray spectrum of a quiescent neutron star transient. The 0.3-79 keV luminosity was 1.1 x 10^(33) erg/s (for D=1kpc), with around 60 percent in the thermal component. We clearly detect a cutoff of the hard spectral tail above 10 keV, the first time such a feature has been detected in this source class. We show that thermal Comptonization and synchrotron shock origins for the hard X-ray emission are ruled out on physical grounds. However, the hard X-ray spectrum is well fit by a thermal bremsstrahlung model with an 18 keV electron temperature, which can be understood as arising either in a hot layer above the neutron star atmosphere or in a radiatively-inefficient accretion flow (RIAF). The power-law cutoff energy may be set by the degree of Compton cooling of the bremsstrahlung electrons by thermal seed photons from the neutron star surface. Lower thermal luminosities should lead to higher (possibly undetectable) cutoff energies. We compare Cen~X-4's behavior with PSR J1023+0038, IGR J18245-2452, and XSS J12270-4859, which have shown transitions between LMXB and radio pulsar modes at a similar X-ray luminosity.

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

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.

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.

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.

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

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

The X-ray spectra of the most extreme ultra-luminous X-ray sources -- those with L > 1 E+40 erg/s -- remain something of a mystery. Spectral roll-over in the 5-10 keV band was originally detected in in the deepest XMM-Newton observations of the brightest sources; this is confirmed in subsequent NuSTAR spectra. This emission can be modeled via Comptonization, but with low electron temperatures (kT_e ~ 2 keV) and high optical depths (tau ~ 10) that pose numerous difficulties. Moreover, evidence of cooler thermal emission that can be fit with thin disk models persists, even in fits to joint XMM-Newton and NuSTAR observations. Using NGC 1313 X-1 as a test case, we show that a patchy disk with a multiple temperature profile may provide an excellent description of such spectra. In principle, a number of patches within a cool disk might emit over a range of temperatures, but the data only require a two-temperature profile plus standard Comptonization, or three distinct blackbody components. A mechanism such as the photon bubble instability may naturally give rise to a patchy disk profile, and could give rise to super-Eddington luminosities. It is possible, then, that a patchy disk (rather than a disk with a standard single-temperature profile) might be a hallmark of accretion disks close to or above the Eddington limit. We discuss further tests of this picture, and potential implications for sources such as narrow-line Seyfert-1 galaxies (NLSy1s) and other low-mass active galactic nuclei (AGN).

We present a spectral analysis of three simultaneous NuSTAR and Swift/XRT observations of the transient Be-neutron star binary KS 1947+300 taken during its outburst in 2013/2014. These broad-band observations were supported by Swift/XRT monitoring snap-shots every 3 days, which we use to study the evolution of the spectrum over the outburst. We find strong changes of the power-law photon index, which shows a weak trend of softening with increasing X-ray flux. The neutron star shows very strong pulsations with a period of P ~ 18.8 s. The 0.8-79 keV broad-band spectrum can be described by a power-law with an exponential cutoff and a black-body component at low energies. During the second observation we detect a cyclotron resonant scattering feature at 12.5 keV, which is absent in the phase-averaged spectra of observations 1 and 3. Pulse phase-resolved spectroscopy reveals that the strength of the feature changes strongly with pulse phase and is most prominent during the broad minimum of the pulse profile. At the same phases the line also becomes visible in the first and third observation at the same energy. This discovery implies that KS 1947+300 has a magnetic field strength of B ~ 1.1e12 (1+z)G, which is at the lower end of known cyclotron line sources.

The Nuclear Spectroscopic Telescope Array (NuSTAR) is the first focusing hard X-ray mission in orbit and operates in the 3-79 keV range. NuSTAR's sensitivity is roughly two orders of magnitude better than previous missions in this energy band thanks to its superb angular resolution. Since its launch in 2012 June, NuSTAR has performed excellently and observed many interesting sources including four magnetars, two rotation-powered pulsars and the cataclysmic variable AE Aquarii. NuSTAR also discovered 3.76-s pulsations from the transient source SGR J1745-29 recently found by Swift very close to the Galactic Center, clearly identifying the source as a transient magnetar. For magnetar 1E 1841-045, we show that the spectrum is well fit by an absorbed blackbody plus broken power-law model with a hard power-law photon index of ~1.3. This is consistent with previous results by INTEGRAL and RXTE. We also find an interesting double-peaked pulse profile in the 25-35 keV band. For AE Aquarii, we show that the spectrum can be described by a multi-temperature thermal model or a thermal plus non-thermal model; a multi-temperature thermal model without a non-thermal component cannot be ruled out. Furthermore, we do not see a spiky pulse profile in the hard X-ray band, as previously reported based on Suzaku observations. For other magnetars and rotation-powered pulsars observed with NuSTAR, data analysis results will be soon available.

We present the results of NuSTAR and XMM-Newton observations of the two ultraluminous X-ray sources (ULX) NGC 1313 X-1 and X-2. The combined spectral bandpass of the two satellites enables us to produce the first spectrum of X-1 between 0.3 and 30 keV, while X-2 is not significantly detected by NuSTAR above 10 keV. The NuSTAR data demonstrate that X-1 has a clear cutoff above 10 keV, whose presence was only marginally detectable with previous X-ray observations. This cutoff rules out the interpretation of X-1 as a black hole in a standard low/hard state, and it is deeper than predicted for the downturn of a broadened iron line in a reflection-dominated regime. The cutoff differs from the prediction of a single-temperature Comptonization model. Further, a cold disk-like black body component at ~0.3 keV is required by the data, confirming previous measurements by XMM-Newton only. We observe a spectral transition in X-2, from a state with high luminosity and strong variability to a lower-luminosity state with no detectable variability, and we link this behavior to a transition from a super-Eddington to a sub-Eddington regime.

We present NuSTAR observations of Vela X-1, a persistent, yet highly variable, neutron star high-mass X-ray binary (HMXB). Two observations were taken at similar orbital phases but separated by nearly a year. They show very different 3-79 keV flux levels as well as strong variability during each observation, covering almost one order of magnitude in flux. These observations allow, for the first time ever, investigations on kilo-second time-scales of how the centroid energies of cyclotron resonant scattering features (CRSFs) depend on flux for a persistent HMXB. We find that the line energy of the harmonic CRSF is correlated with flux, as expected in the sub-critical accretion regime. We argue that Vela X-1 has a very narrow accretion column with a radius of around 0.4 km that sustains a Coulomb interaction dominated shock at the observed luminosities of Lx ~ 3x10^36 erg/s. Besides the prominent harmonic line at 55 keV the fundamental line around 25 keV is clearly detected. We find that the strengths of the two CRSFs are anti-correlated, which we explain by photon spawning. This anti-correlation is a possible explanation for the debate about the existence of the fundamental line. The ratio of the line energies is variable with time and deviates significantly from 2.0, also a possible consequence of photon spawning, which changes the shape of the line. During the second observation, Vela X-1 showed a short off-state in which the power-law softened and a cut-off was no longer measurable. It is likely that the source switched to a different accretion regime at these low mass accretion rates, explaining the drastic change in spectral shape.

We report new spectral and temporal observations of the magnetar 1E 1841-045 in the Kes 73 supernova remnant obtained with the Nuclear Spectroscopic Telescope Array (NuSTAR). Combined with new Swift and archival XMM-Newton and Chandra observations, the phase-averaged spectrum is well characterized by a blackbody plus double power-law model, in agreement with previous, multi-mission X-ray results. However, we are unable to reproduce the spectral results reported using Suzaku observations. The pulsed fraction of the source is found to increase with photon energy. The measured rms pulsed fraction is ~12% and ~17% at ~20 keV and ~50 keV, respectively. We detect a new feature in the 24--35 keV band pulse profile that is uniquely double-peaked. This feature may be associated with a possible absorption or emission feature in the phase-resolved spectrum. We fit the X-ray data using the recently developed electron-positron outflow model of Beloborodov (2013) for the hard X-ray emission from magnetars. This produces a satisfactory fit allowing a constraint on the angle between the rotation and magnetic axes of the neutron star of ~20 degrees and on the angle between the rotation axis and line-of-sight of ~50 degrees. In this model, the soft X-ray component is inconsistent with a single blackbody; adding a second blackbody or a power-law component fits the data. The two-blackbody interpretation suggests a hot spot of temperature kT~0.9 keV occupying ~1% of the stellar surface.

We review recent axisymmetric and three-dimensional (3D) magnetohydrodynamic (MHD) numerical simulations of magnetospheric accretion, plasma-field interaction and outflows from the disk-magnetosphere boundary.

The microquasar 1E1740.7-2942, also known as the "Great Annihilator", was observed by NuSTAR in the Summer of 2012. We have analyzed in detail two observations taken ~2 weeks apart, for which we measure hard and smooth spectra typical of the low/hard state. A few weeks later the source flux declined significantly. Nearly simultaneous coverage by INTEGRAL is available from its Galactic Center monitoring campaign lasting ~2.5 months. These data probe the hard state spectrum from 1E1740.7-2942 before the flux decline. We find good agreement between the spectra taken with IBIS/ISGRI and NuSTAR, with the measurements being compatible with a change in flux with no spectral variability. We present a detailed analysis of the NuSTAR spectral and timing data and upper limits for reflection of the high energy emission. We show that the high energy spectrum of this X-ray binary is well described by thermal Comptonization.

We report on an observation of the neutron star low-mass X-ray binary Serpens X-1, made with NuSTAR. The extraordinary sensitivity afforded by NuSTAR facilitated the detection of a clear, robust, relativistic Fe K emission line from the inner disk. A relativistic profile is required over a single Gaussian line from any charge state of Fe at the 5-sigma level of confidence, and any two Gaussians of equal width at the same confidence. The Compton back-scattering "hump" peaking in the 10-20 keV band is detected for the first time in a neutron star X-ray binary. Fits with relativistically-blurred disk reflection models suggest that the disk likely extends close to the innermost stable circular orbit (ISCO) or stellar surface. The best-fit blurred reflection models constrain the gravitational redshift from the stellar surface to be z > 0.16. The data are broadly compatible with the disk extending to the ISCO; in that case, z > 0.22 and R < 12.6 km (assuming M = 1.4 Msun and a=0, where a = cJ/GM^2). If the star is as large or larger than its ISCO, or if the effective reflecting disk leaks across the ISCO to the surface, the redshift constraints become measurements. We discuss our results in the context of efforts to measure fundamental properties of neutron stars, and models for accretion onto compact objects.

The black hole binary Cygnus X-1 was observed in late-2012 with the Nuclear Spectroscopic Telescope Array (NuSTAR) and Suzaku, providing spectral coverage over the ~1-300 keV range. The source was in the soft state with a multi-temperature blackbody, power-law, and reflection components along with absorption from highly ionized material in the system. The high throughput of NuSTAR allows for a very high quality measurement of the complex iron line region as well as the rest of the reflection component. The iron line is clearly broadened and is well-described by a relativistic blurring model, providing an opportunity to constrain the black hole spin. Although the spin constraint depends somewhat on which continuum model is used, we obtain a*>0.83 for all models that provide a good description of the spectrum. However, none of our spectral fits give a disk inclination that is consistent with the most recently reported binary values for Cyg X-1. This may indicate that there is a >13 degree misalignment between the orbital plane and the inner accretion disk (i.e., a warped accretion disk) or that there is missing physics in the spectral models.

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

Her X-1, one of the brightest and best studied X-ray binaries, shows a cyclotron resonant scattering feature (CRSF) near 37 keV. This makes it an ideal target for detailed study with the Nuclear Spectroscopic Telescope Array (NuSTAR), taking advantage of its excellent hard X-ray spectral resolution. We observed Her X-1 three times, coordinated with Suzaku, during one of the high flux intervals of its 35d super-orbital period. This paper focuses on the shape and evolution of the hard X-ray spectrum. The broad-band spectra can be fitted with a powerlaw with a high-energy cutoff, an iron line, and a CRSF. We find that the CRSF has a very smooth and symmetric shape, in all observations and at all pulse-phases. We compare the residuals of a line with a Gaussian optical depth profile to a Lorentzian optical depth profile and find no significant differences, strongly constraining the very smooth shape of the line. Even though the line energy changes dramatically with pulse phase, we find that its smooth shape does not. Additionally, our data show that the continuum is only changing marginally between the three observations. These changes can be explained with varying amounts of Thomson scattering in the hot corona of the accretion disk. The average, luminosity-corrected CRSF energy is lower than in past observations and follows a secular decline. The excellent data quality of NuSTAR provides the best constraint on the CRSF energy to date.

We observed persistent high-frequency oscillations of the boundary layer near an accreting, weakly-magnetized star in global 3D MHD simulations. The tilted dipole magnetic field is not strong enough to open a gap between the star and the disk. Instead, it forms a highly-wrapped azimuthal field near the surface of the star which slows down rotation of the disk matter, while a small tilt of the field excites oscillations of the boundary layer with a frequency below the Keplerian frequency. This mechanism may be responsible for the high-frequency oscillations in accreting neutron stars, white dwarfs and classical T Tauri stars.

We report on the results of spectral fits made to a NuSTAR observation of the black hole GRS 1915+105 in a "plateau" state. This state is of special interest because it is similar to the "low/hard" state seen in other black holes, especially in that compact, steady jets are launched in this phase. The 3-79 keV bandpass of NuSTAR, and its ability to obtain moderate-resolution spectra free from distortions such as photon pile-up, are extremely well suited to studies of disk reflection in X-ray binaries. In only 15 ks of net exposure, an extraordinarily sensitive spectrum of GRS 1915+105 was measured across the full bandpass. Ionized reflection from a disk around a rapidly-spinning black hole is clearly required to fit the spectra; even hybrid Comptonization models including ionized reflection from a disk around a Schwarzschild black hole proved inadequate. A spin parameter of a = 0.98 +/- 0.01 (1-sigma statistical error) is measured via the best-fit model; low spins are ruled out at a high level of confidence. This result suggests that jets can be launched from a disk extending to the innermost stable circular orbit. A very steep inner disk emissivity profile is also measured, consistent with models of compact coronae above Kerr black holes. These results support an emerging association between the hard X-ray corona and the base of the relativistic jet.

The NuSTAR hard X-ray telescope observed the transient Be/X-ray binary GS 0834-430 during its 2012 outburst - the first active state of this system observed in the past 19 years. We performed timing and spectral analysis, and measured the X-ray spectrum between 3-79keV with high statistical significance. We find the phase-averaged spectrum to be consistent with that observed in many other magnetized accreting pulsars. We fail to detect cyclotron resonance scattering features that would allow us to constrain the pulsar's magnetic field in either phase-averaged or phase-resolved spectra. Timing analysis shows a clearly detected pulse period of ~12.29s in all energy bands. The pulse profiles show a strong, energy-dependent hard phase lag of up to 0.3 cycles in phase, or about 4s. Such dramatic energy-dependent lags in the pulse profile have never before been reported in high-mass X-ray binary (HMXB) pulsars. Previously reported lags have been significantly smaller in phase and restricted to low-energies (E<10keV). We investigate the possible mechanisms that might produce this energy-dependent pulse phase shift. We find the most likely explanation for this effect to be a complex beam geometry.

The Nuclear Spectroscopic Telescope Array (NuSTAR) mission, launched on 13 June 2012, is the first focusing high-energy X-ray telescope in orbit. NuSTAR operates in the band from 3 -- 79 keV, extending the sensitivity of focusing far beyond the ~10 keV high-energy cutoff achieved by all previous X-ray satellites. The inherently low-background associated with concentrating the X-ray light enables NuSTAR to probe the hard X-ray sky with a more than one-hundred-fold improvement in sensitivity over the collimated or coded-mask instruments that have operated in this bandpass. Using its unprecedented combination of sensitivity, spatial and spectral resolution, NuSTAR will pursue five primary scientific objectives, and will also undertake a broad program of targeted observations. The observatory consists of two co-aligned grazing-incidence X-ray telescopes pointed at celestial targets by a three-axis stabilized spacecraft. Deployed into a 600 km, near-circular, 6degree inclination orbit, the Observatory has now completed commissioning, and is performing consistent with pre-launch expectations. NuSTAR is now executing its primary science mission, and with an expected orbit lifetime of ten years, we anticipate proposing a guest investigator program, to begin in Fall 2014.

We investigate the quality factor and RMS amplitude of the lower kHz QPOs from XTE J1701-462, a unique X-ray source which was observed in both the so-called Z and atoll states. Correcting for the frequency drift of the QPO, we show that, as in all sources for which such a correction can be applied, the quality factor and RMS amplitude drops sharply above above a critical frequency. For XTE J1701-462 this frequency is estimated to be ~800 Hz, where the quality factor reaches a maximum of ~200 (e.g. a value consistent with the one observed from more classical systems, such as 4U~1636-536). Such a drop has been interpreted as the signature of the innermost stable circular orbit, and that interpretation is consistent with the observations we report here. The kHz QPOs in the Z state are much less coherent and lower amplitude than they are in the atoll state. We argue that the change of the QPO properties between the two source states is related to the change of the scale height of the accretion disk; a prediction of the toy model proposed by barret et al. (2007). As a by-product of our analysis, we also increased the significance of the upper kHz QPO detected in the atoll phase up to 4.8 sigma (single trial significance), and show that the frequency separation (266.5+/-13.1 Hz) is comparable with the one measured from simultaneous twin QPOs the Z phase.

We present recent results of 3D magnetohydrodynamic simulations of neutron stars with small misalignment angles, as regards the features in lightcurves produced by regular movements of the hot spots during accretion onto the star. In particular, we show that the variation of position of the hot spot created by the infalling matter, as observed in 3D simulations, can produce high frequency Quasi Periodic Oscillations with frequencies associated with the inner zone of the disk. Previously reported simulations showed that the usual assumption of a fixed hot spot near the polar region is valid only for misalignment angles relatively large. Otherwise, two phenomena challenge the assumption: one is the presence of Rayleigh-Taylor instabilities at the disk-magnetospheric boundary, which produce tongues of accreting matter that can reach the star almost anywhere between the equator and the polar region; the other one is the motion of the hot spot around the magnetic pole during stable accretion. In this paper we start by showing that both phenomena are capable of producing short-term oscillations in the lightcurves. We then use Monte Carlo techniques to produce model lightcurves based on the features of the movements observed, and we show that the main features of kHz QPOs can be reproduced. Finally, we show the behavior of the frequencies of the moving spots as the mass accretion rate changes, and propose a mechanism for the production of double QPO peaks.

3D MHD simulation of accretion onto neutron stars have shown in the last few years that the footprint (hotspot) of the accretion flow changes with time. Two different kinds of accretion, namely the funnel flow and the equatorial accretion produced by instabilities at the inner disk, produce different kinds of motion of the hotspot. The funnel flow produces hotspots that move around the magnetic pole, while instabilities produce other hotspots that appear randomly and move along the equator or slightly above. The angular velocities of the two hotspots are different, the equatorial one being higher and both close to the Keplerian velocity in the inner region. Modeling of the lightcurves of these hotspots with Monte Carlo simulations show that the signatures produced in power specra by them, if observed, are QPOs plus low frequency components. Their frequencies, general behavior and features describe correctly most of the properties of kHz QPOs, if we assume the funnel flow hotspots as the origin of the lower kHz QPO and instabilities as the origin of the upper kHz QPO.

{Timing analysis of Accretion-powered Millisecond Pulsars (AMPs) is a powerful tool to probe the physics of compact objects. The recently discovered \newigrj is the 12 discovered out of the 13 AMPs known. The Rossi XTE satellite provided an extensive coverage of the 25 days-long observation of the source outburst.} {Our goal is to investigate the complex interaction between the neutron star magnetic field and the accretion disk, determining the angular momentum exchange between them. The presence of a millisecond coherent flux modulation allows us to investigate such interaction from the study of pulse arrival times. In order to separate the neutron star proper spin frequency variations from other effects, a precise set of orbital ephemeris is mandatory.} {Using timing techniques, we analysed the pulse phase delays fitting differential corrections to the orbital parameters. To remove the effects of pulse phase fluctuations we applied the timing technique already successfully applied to the case of an another AMP, XTE J1807-294.} {We report a precise set of orbital ephemeris. We demonstrate that the companion star is a main sequence star. We find pulse phase delays fluctuations on the first harmonic with a characteristic amplitude of about 0.05, similar to what also observed in the case of the AMP XTE J1814-338. For the second time an AMP shows a third harmonic detected during the entire outburst. The first harmonic phase delays show a puzzling behaviour, while the second harmonic phase delays show a clear spin-up. Also the third harmonic shows a spin-up, although not highly significant (3$\sigma$ c.l.). The presence of a fourth harmonic is also reported. In the hypothesis that the second harmonic is a good tracer of the spin frequency of the neutron star, we find a mean spin frequency derivative for this source of \np{1.65(18)}{-13} Hz s$^{-1}$.} (continue ...)

3D Magnetohydrodynamic simulations show that when matter accretes onto neutron stars, in particular if the misalignment angle is small, it does not constantly fall at a fixed spot. Instead, the location at which matter reaches the star moves. These moving hot spots can be produced both during stable accretion, where matter falls near the magnetic poles of the star, and unstable accretion, characterized by the presence of several tongues of matter which fall on the star near the equator, due to Rayleigh-Taylor instabilities. Precise modeling with Monte Carlo simulations shows that those movements could be observed as high frequency Quasi Periodic Oscillations. We performed a number of new simulation runs with a much wider set of parameters, focusing on neutron stars with a small misalignment angle. In most cases we observe oscillations whose frequency is correlated with the mass accretion rate $\dot{M}$. Moreover, in some cases double QPOs appear, each of them showing the same correlation with $\dot{M}$.