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.

Most ultraluminous X-ray sources (ULXs) are thought to be powered by neutron stars and black holes accreting beyond the Eddington limit. If the compact object is a black hole or a neutron star with a magnetic field $\lesssim10^{12}$ G, the accretion disc is expected to thicken and launch powerful winds driven by radiation pressure. Evidence of such winds has been found in ULXs through the high-resolution spectrometers onboard XMM-Newton, but several unknowns remain, such as the geometry and launching mechanism of these winds. In order to better understand ULX winds and their link to the accretion regime, we have undertaken a major campaign with XMM-Newton to study the ULX NGC 1313 X-1, which is known to exhibit strong emission and absorption features from a mildly-relativistic wind. The new observations show clear changes in the wind with a significantly weakened fast component (0.2c) and the rise of a new wind phase which is cooler and slower (0.06-0.08c). We also detect for the first time variability in the emission lines which indicates an origin within the accretion disc or in the wind. We describe the variability of the wind in the framework of variable super-Eddington accretion rate and discuss a possible geometry for the accretion disc.

High-mass X-ray binaries (HMXBs) are exceptional astrophysical laboratories that offer a rare glimpse into the physical processes that govern accretion on compact objects, massive-star winds, and stellar evolution. In a subset of the HMXBs, the compact objects accrete matter solely from winds of massive donor stars. These so-called wind-fed HMXBs are divided in persistent HMXBs and supergiant fast X-ray transients (SFXTs) according to their X-ray properties. While it has been suggested that this dichotomy depends on the characteristics of stellar winds, they have been poorly studied. With this investigation, we aim to remedy this situation by systematically analyzing donor stars of wind-fed HMXBs that are observable in the UV, concentrating on those with neutron star (NS) companions. We obtained Swift X-ray data, HST UV spectra, and additional optical spectra for all our targets. Our multi-wavelength approach allows us to provide stellar and wind parameters for six donor stars (four wind-fed systems and two OBe X-ray binaries). The wind properties are in line with the predictions of the line-driven wind theory. Three of the donor stars are in an advanced evolutionary stage, while for some of the stars, the abundance pattern indicates that processed material might have been accreted. When passing by the NS in its tight orbit, the donor star wind has not yet reached its terminal velocity but it is still significantly slower; its speed is comparable with the orbital velocity of the NS companion. There are no systematic differences between the two types of wind-fed HMXBs (persistent versus transients) with respect to the donor stars. For the SFXTs in our sample, the orbital eccentricity is decisive for their transient X-ray nature. Based on the orbital parameters and the further evolution of the donor stars, the investigated HMXBs will presumably form Thorne-\.Zytkow objects in the future.

Ultraluminous X-ray sources (ULXs) are extragalactic X-ray emitters located off-center of their host galaxy and with a luminosity in excess of a few ${10^{39}\text{ erg s}^{-1}}$, if emitted isotropically. The discovery of periodic modulation revealed that in some ULXs the accreting compact object is a neutron star, indicating luminosities substantially above their Eddington limit. The most extreme object in this respect is ${NGC 5907~ULX-1}$ (ULX1), with a peak luminosity that is 500 times its Eddington limit. During a Chandra observation to probe a low state of ULX1, we detected diffuse X-ray emission at the position of ULX1. Its diameter is $2.7 \pm 1.0$ arcsec and contains 25 photons, none below 0.8 keV. We interpret this extended structure as an expanding nebula powered by the wind of ULX1. Its diameter of about ${200\text{ pc}}$, characteristic energy of ${\sim 1.9\text{ keV}}$, and luminosity of ${\sim 2\times10^{38}\text{ erg s}^{-1}}$ imply a mechanical power of ${1.3\times10^{41}\text{ erg s}^{-1}}$ and an age ${\sim 7 \times 10^{4}\text{ yr}}$. This interpretation suggests that a genuinely super-Eddington regime can be sustained for time scales much longer than the spin-up time of the neutron star powering the system. As the mechanical power from a single ULX nebula can rival the injection rate of cosmic rays of an entire galaxy, ULX nebulae could be important cosmic ray accelerators.

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.

X-ray radiation from accreting compact objects is an important part of stellar feedback. The metal-poor galaxy ESO 338-4 has experienced vigorous starburst during the last 40 Myr and contains some of the most massive super star clusters in the nearby Universe. Given its starburst age and its star-formation rate, ESO 338-4 is one of the most efficient nearby manufactures of neutron stars and black holes, hence providing an excellent laboratory for feedback studies. We compared X-ray images and spectra obtained by XMM-Newton and Chandra telescopes with integral field spectroscopic VLT MUSE observations in the optical to constrain the nature of strong X-ray emitters. X-ray observations uncover three ultraluminous X-ray sources (ULXs) in ESO 338-4. The brightest among them, ESO 338~X-1, has X-ray luminosity in excess of 10^{40} erg/s. We speculate that ESO 338-4 is powered by accretion on an intermediate-mass (~300Msun) black hole. We show that X-ray radiation from ULXs and hot superbubbles strongly contributes to HeII ionization and general stellar feedback in this template starburst galaxy.

We report on Nuclear Spectroscopic Telescope Array (NuSTAR), Neil Gehrels Swift Observatory (Swift) X-ray Telescope (XRT) and Swift Burst Alert Telescope (BAT) observations of IGR J16493-4348, a wind-fed Supergiant X-ray Binary (SGXB) showing significant superorbital variability. From a discrete Fourier transform of the BAT light curve, we refine its superorbital period to be 20.058 $\pm$ 0.007 days. The BAT dynamic power spectrum and a fractional root mean square analysis both show strong variations in the amplitude of the superorbital modulation, but no observed changes in the period were found. The superorbital modulation is significantly weaker between MJD 55,700 and MJD 56,300. The joint NuSTAR and XRT observations, which were performed near the minimum and maximum of one cycle of the 20 day superorbital modulation, show that the flux increases by more than a factor of two between superorbital minimum and maximum. We find no significant changes in the 3-50 keV pulse profiles between superorbital minimum and maximum, which suggests a similar accretion regime. Modeling the pulse-phase averaged spectra we find a possible Fe K$\alpha$ emission line at 6.4 keV at superorbital maximum. The feature is not significant at superorbital minimum. While we do not observe any significant differences between the pulse-phase averaged spectral continua apart from the overall flux change, we find that the hardness ratio near the broad main peak of the pulse profile increases from superorbital minimum to maximum. This suggests the spectral shape hardens with increasing luminosity. We discuss different mechanisms that might drive the observed superorbital modulation.

We present an analysis of the spectral shape and pulse profile of the accretion-powered pulsar 4U 1626-67 observed with Suzaku and NuSTAR during a spin-up state. The pulsar, which experienced a torque reversal to spin-up in 2008, has a spin period of 7.7 s. Comparing the phase-averaged spectra obtained with Suzaku in 2010 and with NuSTAR in 2015, we find that the spectral shape changed between the two observations: the 3-10 keV flux increased by 5% while the 30-60 keV flux decreased significantly by 35%. Phase-averaged and phase-resolved spectral analysis shows that the continuum spectrum observed by NuSTAR is well described by an empirical NPEX continuum with an added broad Gaussian emission component around the spectral peak at 20 keV. Taken together with the observed Pdot value obtained from Fermi/GBM, we conclude that the spectral change between the Suzaku and NuSTAR observations was likely caused by an increase of the accretion rate. We also report the possible detection of asymmetry in the profile of the fundamental cyclotron line. Furthermore, we present a study of the energy-resolved pulse profiles using a new relativistic ray tracing code, where we perform a simultaneous fit to the pulse profiles assuming a two-column geometry with a mixed pencil- and fan-beam emission pattern. The resulting pulse profile decompositions enable us to obtain geometrical parameters of accretion columns (inclination, azimuthal and polar angles) and a fiducial set of beam patterns. This information is important to validate the theoretical predictions from radiation transfer in a strong magnetic field.

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.

Strongly magnetized, accreting neutron stars show periodic and aperiodic variability over a wide range of time scales. By obtaining spectral and timing information on these different time scales, we can have a closer look into the physics of accretion close to the neutron star and the properties of the accreted material. One of the most prominent time scales is the strong pulsation, i.e., the rotation period of the neutron star itself. Over one rotation, our view of the accretion column and the X-ray producing region changes significantly. This allows us to sample different physical conditions within the column but at the same time requires that we have viewing-angle-resolved models to properly describe them. In wind-fed high-mass X-ray binaries, the main source of aperiodic variability is the clumpy stellar wind, which leads to changes in the accretion rate (i.e., luminosity) as well as absorption column. This variability allows us to study the behavior of the accretion column as a function of luminosity, as well as to investigate the structure and physical properties of the wind, which we can compare to winds in isolated stars.

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.

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.

We present two observations of the high-mass X-ray binary GX 301-2 with NuSTAR, taken at different orbital phases and different luminosities. We find that the continuum is well described by typical phenomenological models, like a very strongly absorbed NPEX model. However, for a statistically acceptable description of the hard X-ray spectrum we require two cyclotron resonant scattering features (CRSF), one at ~35 keV and the other at ~50 keV. Even though both features strongly overlap, the good resolution and sensitivity of NuSTAR allows us to disentangle them at >=99.9% significance. This is the first time that two CRSFs are seen in GX 301-2. We find that the CRSFs are very likely independently formed, as their energies are not harmonically related and, if it were a single line, the deviation from a Gaussian shape would be very large. We compare our results to archival Suzaku data and find that our model also provides a good fit to those data. We study the behavior of the continuum as well as the CRSF parameters as function of pulse phase in seven phase bins. We find that the energy of the 35 keV CRSF varies smoothly as function of phase, between 30-38 keV. To explain this variation, we apply a simple model of the accretion column, taking the altitude of the line-forming region, the velocity of the in-falling material, and the resulting relativistic effects into account. We find that in this model the observed energy variation can be explained simply due to a variation of the projected velocity and beaming factor of the line forming region towards us.

We present a spectral and timing analysis of the newly reported Galactic X-ray transient Swift J1658.2-4242 observed by NuSTAR and Swift. The broad-band X-ray continuum is typical of a black hole binary in the bright hard state, with a photon index of $\Gamma=1.63\pm0.02$ and a low coronal temperature of $kT_{\rm e}=22\pm1$ keV, corresponding to a low spectral cutoff well constrained by NuSTAR. Spectral modeling of the relativistic disk reflection features, consisting of a broad Fe K$\alpha$ line and the Compton reflection hump, reveals that the black hole is rapidly spinning with the spin parameter of $a^{*}>0.96$, and the inner accretion disk is viewed at a high inclination angle of $i=64^{+2}_{-3}{^\circ}$ (statistical errors, 90% confidence). The high inclination is independently confirmed by dips in the light curves, which can be explained by absorbing material located near the disk plane temporarily obscuring the central region. In addition, we detect an absorption line in the NuSTAR spectra centered at $7.03^{+0.04}_{-0.03}$ keV. If associated with ionized Fe K absorption lines, this provides evidence for the presence of outflowing material in the low/hard state of a black hole binary candidate. A timing analysis shows the presence of a type-C QPO in the power spectrum, with the frequency increasing from $\sim0.14$ Hz to $\sim0.21$ Hz during the single NuSTAR exposure. Our analysis reveals that Swift J1658.2-4242 displays characteristics typical for a black hole binary that is viewed at a high inclination angle, making it a good system for studying the accretion geometry in black hole binaries.

We report on the serendipitous discovery of a new transient in NGC 5907, at a peak luminosity of 6.4x10^{39} erg/s. The source was undetected in previous 2012 Chandra observations with a 3 sigma upper limit on the luminosity of 1.5x10^{38} erg/s, implying a flux increase of a factor of >35. We analyzed three recent 60ks/50ks Chandra and 50ks XMM-Newton observations, as well as all the available Swift observations performed between August 2017/March 2018. Until the first half of October 2017, Swift observations do not show any emission from the source. The transient entered the ULX regime in less than two weeks and its outburst was still on-going at the end of February 2018. The 0.3-10 keV spectrum is consistent with a single multicolour blackbody disc (kT~1.5 keV). The source might be a ~30 solar mass black hole accreting at the Eddington limit. However, although we did not find evidence of pulsations, we cannot rule-out the possibility that this ULX hosts an accreting neutron star.

We present a timing analysis of multiple XMM-Newton and NuSTAR observations of the ultra-luminous pulsar NGC 7793 P13 spread over its 65d variability period. We use the measured pulse periods to determine the orbital ephemeris, confirm a long orbital period with P_orb = 63.9 (+0.5,-0.6) d, and find an eccentricity of e <= 0.15. The orbital signature is imprinted on top of a secular spin-up, which seems to get faster as the source becomes brighter. We also analyse data from dense monitoring of the source with Swift and find an optical photometric period of 63.9 +/- 0.5 d and an X-ray flux period of 66.8 +/- 0.4 d. The optical period is consistent with the orbital period, while the X-ray flux period is significantly longer. We discuss possible reasons for this discrepancy, which could be due to a super-orbital period caused by a precessing accretion disk or an orbital resonance. We put the orbital period of P13 into context with the orbital periods implied for two other ultra-luminous pulsars, M82 X-2 and NGC 5907 ULX and discuss possible implications for the system parameters.

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.

We report on a NuSTAR observation of the recently discovered bright black hole candidate MAXI J1535-571. NuSTAR observed the source on MJD 58003 (five days after the outburst was reported). The spectrum is characteristic of a black hole binary in the hard state. We observe clear disk reflection features, including a broad Fe K$\alpha$ line and a Compton hump peaking around 30 keV. Detailed spectral modeling reveals narrow Fe K$\alpha$ line complex centered around 6.5 keV on top of the strong relativistically broadened Fe K$\alpha$ line. The narrow component is consistent with distant reflection from moderately ionized material. The spectral continuum is well described by a combination of cool thermal disk photons and a Comptonized plasma with the electron temperature $kT_{\rm e}=19.7\pm{0.4}$ keV. An adequate fit can be achieved for the disk reflection features with a self-consistent relativistic reflection model that assumes a lamp-post geometry for the coronal illuminating source. The spectral fitting measures a black hole spin $a>0.84$, inner disk radius $R_{\rm in}<2.01~r_{\rm ISCO}$, and a lamp-post height $h=7.2^{+0.8}_{-2.0} r_{\rm g}$ (statistical errors, 90% confidence), indicating no significant disk truncation and a compact corona. Although the distance and mass of this source are not currently known, this suggests the source was likely in the brighter phases of the hard state during this NuSTAR observation.

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.

We present results from a coordinated $XMM$-$Newton$+$NuSTAR$ observation of the type 1.8 Seyfert galaxy IRAS 13197-1627. This is a highly complex source, with strong contributions from relativistic reflection from the inner accretion disk, neutral absorption and further reprocessing by more distant material, and ionised absorption from an outflow. We undertake a detailed spectral analysis combining the broadband coverage provided by $XMM$-$Newton$+$NuSTAR$ with a multi-epoch approach incorporating archival observations performed by $XMM$-$Newton$ and $Suzaku$. Our focus is on characterising the reflection from the inner accretion disk, which previous works have suggested may dominate the AGN emission, and constraining the black hole spin. Using lamppost disk reflection models, we find that the results for the inner disk are largely insensitive to assumptions regarding the geometry of the distant reprocessor and the precise form of the illuminating X-ray continuum. However, these results do depend on the treatment of the iron abundance of the distant absorber/reprocessor. The multi-epoch data favour a scenario in which the AGN is chemically homogeneous, and we find that a rapidly rotating black hole is preferred, with $a^* \geq 0.7$, but a slowly-rotating black hole is not strongly excluded. In addition to the results for the inner disk, we also find that both the neutral and ionised absorbers vary from epoch to epoch, implying that both have some degree of inhomogeneity in their structure.

Whilst astronomy as a science is historically founded on observations at optical wavelengths, studying the Universe in other bands has yielded remarkable discoveries, from pulsars in the radio, signatures of the Big Bang at submm wavelengths, through to high energy emission from accreting, gravitationally-compact objects and the discovery of gamma-ray bursts. Unsurprisingly, the result of combining multiple wavebands leads to an enormous increase in diagnostic power, but powerful insights can be lost when the sources studied vary on timescales shorter than the temporal separation between observations in different bands. In July 2015, the workshop "Paving the way to simultaneous multi-wavelength astronomy" was held as a concerted effort to address this at the Lorentz Center, Leiden. It was attended by 50 astronomers from diverse fields as well as the directors and staff of observatories and spaced-based missions. This community white paper has been written with the goal of disseminating the findings of that workshop by providing a concise review of the field of multi-wavelength astronomy covering a wide range of important source classes, the problems associated with their study and the solutions we believe need to be implemented for the future of observational astronomy. We hope that this paper will both stimulate further discussion and raise overall awareness within the community of the issues faced in a developing, important field.

The Be X-ray binary EXO 2030+375 was in an extended low luminosity state during most of 2016. We observed this state with NuSTAR and Swift, supported by INTEGRAL observations as well as optical spectroscopy with the NOT. We present a comprehensive spectral and timing analysis of these data here to study the accretion geometry and investigate a possible onset of the propeller effect. The H-alpha data show that the circumstellar disk of the Be-star is still present. We measure equivalent widths similar to values found during more active phases in the past, indicating that the low-luminosity state is not simply triggered by a smaller Be disk. The NuSTAR data, taken at a 3-78 keV luminosity of ~6.8e35 erg/s (for a distance of 7.1 kpc), are well described by standard accreting pulsar models, such as an absorbed power-law with a high-energy cutoff. We find that pulsations are still clearly visible at these luminosities, indicating that accretion is continuing despite the very low mass transfer rate. In phase-resolved spectroscopy we find a peculiar variation of the photon index from ~1.5 to ~2.5 over only about 3% of the rotational period. This variation is similar to that observed with XMM-Newton at much higher luminosities. It may be connected to the accretion column passing through our line of sight. With Swift/XRT we observe luminosities as low as 1e34 erg/s during which the data quality did not allow us to search for pulsations, but the spectrum is much softer and well described by either a blackbody or soft power-law continuum. This softer spectrum might be due to the fact that accretion has been stopped by the propeller effect and we only observe the neutron star surface cooling.

The Be X-ray binary pulsar EXO 2030+375, first detected in 1985, has shown a significant detected X-ray outburst at nearly every periastron passage of its 46-day orbit for the past ~25 years, with one low state accompanied by a torque reversal in the 1990s. In early 2015 the outbursts progressively became fainter and less regular while the monotonic spin-up flattened. At the same time a decrease in the H$\alpha$ line equivalent width was reported, indicating a change in the disk surrounding the mass donor. In order to explore the source behaviour in the poorly explored low-flux state with a possible transition to a state of centrifugal inhibition of accretion we have undertaken an observing campaign with Swift/XRT, NuSTAR and the Nordic Optical Telescope (NOT). This conference contribution reports the preliminary results obtained from our campaign.

We present results from four new broadband X-ray observations of the extreme ultraluminous X-ray source Holmberg IX X-1 ($L_{\rm{X}} > 10^{40}$ erg/s), performed by $Suzaku$ and $NuSTAR$ in coordination. Combined with the archival data, we now have broadband observations of this remarkable source from six separate epochs. Two of these new observations probe lower fluxes than seen previously, allowing us to extend our knowledge of the broadband spectral variability exhibited. The spectra are well fit by two thermal blackbody components, which dominate the emission below 10 keV, as well as a steep ($\Gamma \sim 3.5$) powerlaw tail which dominates above $\sim$15 keV. Remarkably, while the 0.3-10.0 keV flux varies by a factor of $\sim$3 between all these epochs, the 15-40 keV flux varies by only $\sim$20%. Although the spectral variability is strongest in the $\sim$1-10 keV band, both of the thermal components are required to vary when all epochs are considered. We also re-visit the search for iron absorption features, leveraging the high-energy $NuSTAR$ data to improve our sensitivity to extreme velocity outflows in light of the ultra-fast outflow recently detected in NGC 1313 X-1. Iron absorption from a similar outflow along our line of sight can be ruled out in this case. We discuss these results in the context of super-Eddington accretion models that invoke a funnel-like geometry for the inner flow, and propose a scenario in which we have an almost face-on view of a funnel that expands to larger radii with increasing flux, resulting in an increasing degree of geometrical collimation for the emission from intermediate temperature regions.

We present first results from a series of $NuSTAR$ observations of the black hole X-ray binary V404 Cyg obtained during its summer 2015 outburst, primarily focusing on observations during the height of this outburst activity. The $NuSTAR$ data show extreme variability in both the flux and spectral properties of the source. This is partly driven by strong and variable line-of-sight absorption, similar to previous outbursts. The latter stages of this observation are dominated by strong flares, reaching luminosities close to Eddington. During these flares, the central source appears to be relatively unobscured and the data show clear evidence for a strong contribution from relativistic reflection, providing a means to probe the geometry of the innermost accretion flow. Based on the flare properties, analogy with other Galactic black hole binaries, and also the simultaneous onset of radio activity, we argue that this intense X-ray flaring is related to transient jet activity during which the ejected plasma is the primary source of illumination for the accretion disk. If this is the case, then our reflection modelling implies that these jets are launched in close proximity to the black hole (as close as a few gravitational radii), consistent with expectations for jet launching models that tap either the spin of the central black hole, or the very innermost accretion disk. Our analysis also allows us to place the first constraints on the black hole spin for this source, which we find to be $a^* > 0.92$ (99% statistical uncertainty, based on an idealized lamppost geometry).

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

The RXTE observed four outbursts of the accreting X-ray binary transient source, GX 304-1 in 2010 and 2011. We present results of detailed 3-100 keV spectral analysis of 69 separate observations, and report a positive correlation between cyclotron line parameters, as well as other spectral parameters, with power law flux. The cyclotron line energy, width and depth versus flux, and thus luminosity, correlations show a flattening of the relationships with increasing luminosity, which are well described by quasi-spherical or disk accretion that yield the surface magnetic field to be ~60 keV. Since HEXTE cluster A was fixed aligned with the PCA field of view and cluster B was fixed viewing a background region 1.5 degrees off of the source direction during these observations near the end of the RXTE mission, the cluster A background was estimated from cluster B events using HEXTEBACKEST. This made possible the detection of the ~55 keV cyclotron line and an accurate measurement of the continuum. Correlations of all spectral parameters with the primary 2-10 keV power law flux reveal it to be the primary driver of the spectral shape. The accretion is found to be in the collisionless shock braking regime.

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

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 detailed spectral analysis of XMM-Newton and NuSTAR observations of the accreting transient black hole GRS 1739-278 during a very faint low hard state at ~0.02% of the Eddington luminosity (for a distance of 8.5 kpc and a mass of 10 M_sun ). The broad-band X-ray spectrum between 0.5-60 keV can be well-described by a power law continuum with an exponential cutoff. The continuum is unusually hard for such a low luminosity, with a photon index of Gamma = 1.39 +/- 0.04. We find evidence for an additional reflection component from an optically thick accretion disk at the 98% likelihood level. The reflection fraction is low with R_refl = 0.043(+0.033,-0.023). In combination with measurements of the spin and inclination parameters made with NuSTAR during a brighter hard state by Miller and co-workers, we seek to constrain the accretion disk geometry. Depending on the assumed emissivity profile of the accretion disk, we find a truncation radius of 15-35 Rg (5-12 R_ISCO ) at the 90% confidence limit. These values depend strongly on the assumptions and we discuss possible systematic uncertainties.

Electron cyclotron resonant scattering features (CRSFs) are observed as absorption-like lines in the spectra of X-ray pulsars. A significant fraction of the computing time for Monte Carlo simulations of these quantum mechanical features is spent on the calculation of the mean free path for each individual photon before scattering, since it involves a complex numerical integration over the scattering cross section and the (thermal) velocity distribution of the scattering electrons. We aim to numerically calculate interpolation tables which can be used in CRSF simulations to sample the mean free path of the scattering photon and the momentum of the scattering electron. The tables also contain all the information required for sampling the scattering electron's final spin. The tables were calculated using an adaptive Simpson integration scheme. The energy and angle grids were refined until a prescribed accuracy is reached. The tables are used by our simulation code to produce artificial CRSF spectra. The electron momenta sampled during these simulations were analyzed and justified using theoretically determined boundaries. We present a complete set of tables suited for mean free path calculations of Monte Carlo simulations of the cyclotron scattering process for conditions expected in typical X-ray pulsar accretion columns (0.01<B/B_{crit}<=0.12, where B_{crit}=4.413x10^{13} G and 3keV<=kT<15keV). The sampling of the tables is chosen such that the results have an estimated relative error of at most 1/15 for all points in the grid. The tables are available online at http://www.sternwarte.uni-erlangen.de/research/cyclo.

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 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 present the results of a NuSTAR study of the dynamically confirmed stellar-mass black hole GS 1354-645. The source was observed during its 2015 "hard" state outburst; we concentrate on spectra from two relatively bright phases. In the higher-flux observation, the broadband NuSTAR spectra reveal a clear, strong disk reflection spectrum, blurred by a degree that requires a black hole spin of a = cJ/GM^2 > 0.98 (1 sigma statistical limits only). The fits also require a high inclination: theta = 75(2) degrees. Strong "dips" are sometimes observed in the X-ray light curves of sources viewed at such an angle; these are absent, perhaps indicating that dips correspond to flared disk structures that only manifest at higher accretion rates. In the lower-flux observation, there is evidence of radial truncation of the thin accretion disk. We discuss these results in the context of spin in stellar-mass black holes, and inner accretion flow geometries at moderate accretion rates.

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

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

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

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

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.

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 report on a multiwavelength observational campaign of the black hole X-ray binary Swift J1753.5-0127 that consists of an ESO/X-shooter spectrum supported by contemporaneous Swift/XRT+UVOT and ATCA data. ISM absorption lines in the X-shooter spectrum allows us to determine E(B-V)=0.45+/-0.02 along the line-of-sight to the source. We also report detection of emission signatures of He II at 4686 angstrom, H alpha, and, for the first time, H I at 10906 angstrom and Paschen Beta. The double-peaked morphology of these four lines is typical of the chromosphere of a rotating accretion disk. Nonetheless, the paucity of disk features points towards a low level of irradiation in the system. This is confirmed through spectral energy distribution modeling and we find that the UVOT+X-shooter continuum mostly stems from the thermal emission of a viscous disk. We speculate that the absence of reprocessing is due to the compactness of an illumination-induced envelope that fails to reflect enough incoming hard X-ray photons back to the outer regions. The disk also marginally contributes to the Compton-dominated X-ray emission and is strongly truncated, with an inner radius about a thousand times larger than the black hole's gravitational radius. A near-infrared excess is present, and we associate it with synchrotron radiation from a compact jet. However, the measured X-ray flux is significantly higher than what can be explained by the optically thin synchrotron jet component. We discuss these findings in the framework of the radio quiet versus X-ray bright hypothesis, favoring the presence of a residual disk, predicted by evaporation models, that contributes to the X-ray emission without enhancing the radio flux.

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

(abridged) We report on multi-wavelength measurements of Swift J1753.5-0127 in the hard state at L=2.7e36 erg/s (assuming d=3 kpc) in 2014. The radio emission is optically thick synchrotron, presumably from a compact jet. We take advantage of the low extinction and model the near-IR to UV emission with a multi-temperature disk model. Assuming a BH mass of M_BH=5 Msun and a system inclination of 40 deg, the fits imply an inner radius for the disk of Rin/Rg>212 d_3 (5Msun/M_BH). The outer radius is R_out/R_g=90,000 d_3 (5Msun/M_BH), which corresponds to 6.6e10 d_3 cm, consistent with the expected size of the disk. The 0.5-240 keV spectrum measured by Swift/XRT, Suzaku, and NuSTAR is relatively well characterized by a power-law with a photon index of Gamma=1.722+/-0.003, but a significant improvement is seen when a second continuum component is added. Reflection is a possibility, but no iron line is detected, implying a low iron abundance. We are able to fit the entire SED with a multi-temperature disk component, a Comptonization component, and a broken power-law, representing the emission from the compact jet. The broken power-law cannot significantly contribute to the soft X-ray emission, and this may be related to why Swift J1753.5-0127 is an outlier in the radio/X-ray correlation. The broken power-law might dominate above 20 keV, which would constrain the break frequency to be between 2.4e10 Hz and 3.6e12 Hz. Although the fits to the full SED do not include significant thermal emission in the X-ray band, previous observations have consistently seen such a component, and we find that there is evidence at the 3.1-sigma level for a disk-blackbody component with a temperature of 150(+30)(-20) eV and an inner radius of 5-14 R_g. If this component is real, it might imply the presence of an inner accretion disk in addition to the strongly truncated (R_in>212 R_g) disk.

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

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

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

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

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.

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

As of 2014 August, the Nuclear Spectroscopic Telescope Array (NuSTAR) had observed ~30 X-ray binaries either as part of the planned program, as targets of opportunity, or for instrument calibration. The main science goals for the observations include probing the inner part of the accretion disk and constraining black hole spins via reflection components, providing the first observations of hard X-ray emission from quiescent Low Mass X-ray Binaries (LMXBs), measuring cyclotron lines from accreting pulsars, and studying type I X-ray bursts from neutron stars. Here, we describe the science objectives in more depth and give an overview of the NuSTAR observations that have been carried out to achieve the objectives. These include observation of four "IGR" High Mass X-ray Binaries (HMXBs) discovered by INTEGRAL. We also summarize the results that have been obtained and their implications. Among the IGR HMXBs, we focus on the discovery of a cyclotron line in the spectrum of IGR J17544-2619.

Vela X-1 is the prototype of the class of wind-fed accreting pulsars in high mass X-ray binaries hosting a supergiant donor. We have analyzed in a systematic way ten years of INTEGRAL data of Vela X-1 (22-50 keV) and we found that when outside the X-ray eclipse, the source undergoes several luminosity drops where the hard X-rays luminosity goes below 3x10^35 erg/s, becoming undetected by INTEGRAL. These drops in the X-ray flux are usually referred to as "off-states" in the literature. We have investigated the distribution of these off-states along the Vela X-1 ~8.9 d orbit, finding that their orbital occurrence displays an asymmetric distribution, with a higher probability to observe an off-state near the pre-eclipse than during the post-eclipse. This asymmetry can be explained by scattering of hard X-rays in a region of ionized wind, able to reduce the source hard X-ray brightness preferentially near eclipse ingress. We associate this ionized large-scale wind structure with the photoionization wake produced by the interaction of the supergiant wind with the X-ray emission from the neutron star. We emphasize that this observational result could be obtained thanks to the accumulation of a decade of INTEGRAL data, with observations covering the whole orbit several times, allowing us to detect an asymmetric pattern in the orbital distribution of off-states in Vela X-1.

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

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.

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.

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

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 results from the first Suzaku observation of the high-mass X-ray binary 4U 1538-522. The broad-band spectral coverage of Suzaku allows for a detailed spectral analysis, characterizing the cyclotron resonance scattering feature at $23.0 \pm 0.4$ keV and the iron K$\alpha$ line at $6.426 \pm 0.008$ keV, as well as placing limits on the strengths of the iron K$\beta$ line and the iron K edge. We track the evolution of the spectral parameters both in time and in luminosity, notably finding a significant positive correlation between cyclotron line energy and luminosity. A dip and spike in the lightcurve is shown to be associated with an order-of-magnitude increase in column density along the line of sight, as well as significant variation in the underlying continuum, implying the accretion of a overdense region of a clumpy stellar wind. We also present a phase-resolved analysis, with most spectral parameters of interest showing significant variation with phase. Notably, both the cyclotron line energy and the iron K$\alpha$ line intensity vary significantly with phase, with the iron line intensity significantly out-of-phase with the pulse profile. We discuss the implications of these findings in the context of recent work in the areas of accretion column physics and cyclotron resonance scattering feature formation.

The Circinus galaxy is one of the nearest obscured AGN, making it an ideal target for detailed study. Combining archival Chandra and XMM-Newton data with new NuSTAR observations, we model the 2-79 keV spectrum to constrain the primary AGN continuum and to derive physical parameters for the obscuring material. Chandra's high angular resolution allows a separation of nuclear and off-nuclear galactic emission. In the off-nuclear diffuse emission we find signatures of strong cold reflection, including high equivalent-width neutral Fe lines. This Compton-scattered off-nuclear emission amounts to 18% of the nuclear flux in the Fe line region, but becomes comparable to the nuclear emission above 30 keV. The new analysis no longer supports a prominent transmitted AGN component in the observed band. We find that the nuclear spectrum is consistent with Compton-scattering by an optically-thick torus, where the intrinsic spectrum is a powerlaw of photon index $\Gamma = 2.2-2.4$, the torus has an equatorial column density of $N_{\rm H} = (6-10)\times10^{24}$cm$^{-2}$ and the intrinsic AGN $2-10$ keV luminosity is $(2.3-5.1)\times 10^{42}$ erg/s. These values place Circinus along the same relations as unobscured AGN in accretion rate-vs-$\Gamma$ and $L_X$-vs-$L_{IR}$ phase space. NuSTAR's high sensitivity and low background allow us to study the short time-scale variability of Circinus at X-ray energies above 10 keV for the first time. The lack of detected variability favors a Compton-thick absorber, in line with the the spectral fitting results.

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

We present a spectral analysis of four coordinated NuSTAR+XMM-Newton observations of the Seyfert galaxy NGC 1365. These exhibit an extreme level of spectral variability, which is primarily due to variable line-of-sight absorption, revealing relatively unobscured states in this source for the first time. Despite the diverse range of absorption states, each of the observations displays the same characteristic signatures of relativistic reflection from the inner accretion disk. Through time-resolved spectroscopy we find that the strength of the relativistic iron line and the Compton reflection hump relative to the intrinsic continuum are well correlated, as expected if they are two aspects of the same broadband reflection spectrum. We apply self-consistent disk reflection models to these time-resolved spectra in order to constrain the inner disk parameters, allowing for variable, partially covering absorption to account for the vastly different absorption states observed. Each of the four observations is treated independently to test the consistency of the results obtained for the black hole spin and the disk inclination, which should not vary on observable timescales. We find both the spin and the inclination determined from the reflection spectrum to be consistent, confirming NGC 1365 hosts a rapidly rotating black hole; in all cases the dimensionless spin parameter is constrained to be a* > 0.97 (at 90% statistical confidence or better).

The symbiotic X-ray binary 4U 1954+319 is a rare system hosting a peculiar neutron star (NS) and an M-type optical companion. Its ~5.4h NS spin period is the longest among all known accretion-powered pulsars and exhibited large (~7%) fluctuations over 8 years. A spin trend transition was detected with Swift/BAT around an X-ray brightening in 2012. The source was in quiescent and bright states before and after this outburst based on 60 ks Suzaku observations in 2011 and 2012. The observed continuum is well described by a Comptonized model with the addition of a narrow 6.4 keV Fe Kalpha line during the outburst. Spectral similarities to slowly rotating pulsars in high-mass X-ray binaries, its high pulsed fraction (~60-80%), and the location in the Corbet diagram favor high B-field (>~1e+12 G) over a weak field as in low-mass X-ray binaries. The observed low X-ray luminosity (1e+33-1e+35 erg/s), probable wide orbit, and a slow stellar wind of this SyXB make quasi-spherical accretion in the subsonic settling regime a plausible model. Assuming a ~1e+13 G NS, this scheme can explain the ~5.4 h equilibrium rotation without employing the magnetar-like field (~1e+16 G) required in the disk accretion case. The time-scales of multiple irregular flares (~50 s) can also be attributed to the free-fall time from the Alfven shell for a ~1e+13 G field. A physical interpretation of SyXBs beyond the canonical binary classifications is discussed.

GRO J1008$-$57 is a high-mass X-ray binary for which several claims of a cyclotron resonance scattering feature near 80 keV have been reported. We use NuSTAR, Suzaku, and Swift data from its giant outburst of November 2012 to confirm the existence of the 80 keV feature and perform the most sensitive search to date for cyclotron scattering features at lower energies. We find evidence for a 78$^{+3}_{-2}$ keV line in the NuSTAR and Suzaku data at $> 4\sigma$ significance, confirming the detection using Suzaku alone by Yamamoto et al. (2014). A search of both the time-integrated and phase-resolved data rules out a fundamental at lower energies with optical depth larger than 5% of the 78 keV line. These results indicate that GRO J1008$-$57 has a magnetic field of $6.7\times10^{12} (1+z)$ G, the highest of known accreting pulsars.

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.

We present a broad band spectral analysis of the joint XMM-Newton and NuSTAR observational campaign of the Narrow Line Seyfert 1 SWIFT J2127.4+5654, consisting of 300 ks performed during three XMM-Newton orbits. We detect a relativistic broadened iron K$\alpha$ line originating from the innermost regions of the accretion disc surrounding the central black hole, from which we infer an intermediate spin of $a$=$0.58^{+0.11}_{-0.17}$. The intrinsic spectrum is steep ($\Gamma=2.08\pm0.01$) as commonly found in Narrow Line Seyfert 1 galaxies, while the cutoff energy (E$_{\rm c}=108^{+11}_{-10}$ keV) falls within the range observed in Broad Line Seyfert 1 Galaxies. We measure a low-frequency lag that increases steadily with energy, while at high frequencies, there is a clear lag following the shape of the broad Fe K emission line. Interestingly, the observed Fe K lag in SWIFT J2127.4+5654 is not as broad as in other sources that have maximally spinning black holes. The lag amplitude suggests a continuum-to-reprocessor distance of about $ 10-20\ r_{\mathrm{g}}$. These timing results independently support an intermediate black hole spin and a compact corona.

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

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

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.

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]

Mkn 421 is a nearby active galactic nucleus dominated at all wavelengths by a very broad non-thermal continuum thought to arise from a relativistic jet seen at a small angle to the line of sight. Its spectral energy distribution peaks in the X-ray and TeV gamma-ray bands, where the energy output is dominated by cooling of high-energy electrons in the jet. In order to study the electron distribution and its evolution, we carried out a dedicated multi-wavelength campaign, including extensive observations by the recently launched highly sensitive hard X-ray telescope NuSTAR, between December 2012 and May 2013. Here we present some initial results based on NuSTAR data from January through March 2013, as well as calibration observations conducted in 2012. Although the observations cover some of the faintest hard X-ray flux states ever observed for Mkn 421, the sensitivity is high enough to resolve intra-day spectral variability. We find that in this low state the dominant flux variations are smooth on timescales of hours, with typical intra-hour variations of less than 5%. We do not find evidence for either a cutoff in the hard X-ray spectrum, or a rise towards a high-energy component, but rather that at low flux the spectrum assumes a power law shape with a photon index of approximately 3. The spectrum is found to harden with increasing brightness.

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 present a spectral and timing analysis of INTEGRAL observations of two high mass X-ray binaries, 4U 1538-522 and 4U 1907+09. Our timing measurements for 4U 1538-522 find the pulse period to have exhibited a spin-up trend until approximately 2009, after which there is evidence for a torque reversal, with the source beginning to spin down to the most recently-measured period of 525.407 +/- 0.001 s. The most recent INTEGRAL observations of 4U 1907+09 are not found to yield statistically significant pulse periods due to the significantly lower flux from the source compared to 4U 1538-522. A spectral model consisting of a power-law continuum with an exponential cutoff and modified by two cyclotron resonance scattering features is found to fit both sources well, with the cyclotron scattering features detected at ~22 and ~49 keV for 4U 1538-522 and at ~18 and ~36 keV in 4U 1907+09. The spectral parameters of 4U 1538-522 are generally not found to vary significantly with flux, and there is little to no variation across the torque reversal. Examining our results in conjunction with previous work, we find no evidence for a correlation between cyclotron line energy and luminosity for 4U 1538-522. 4U 1907+09 shows evidence for a positive correlation between cyclotron line energy and luminosity, which would make it the fourth, and lowest-luminosity, cyclotron line source to exhibit this relationship

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.

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

We have analyzed 3 observations of the High Mass X-ray Binary A0535+26 performed by the Rossi X-ray Timing Explorer (RXTE) 3, 5, and 6 months after the last outburst in 2011 February. We detect pulsations only in the second observation. The 3-20 keV spectra can be fit equally well with either an absorbed power law or absorbed thermal bremsstrahlung model. Re-analysis of 2 earlier RXTE observations made 4 years after the 1994 outburst, original BeppoSAX observations 2 years later, re-analysis of 4 EXOSAT observations made 2 years after the last 1984 outburst, and a recent XMM-Newton observation in 2012 reveal a stacked, quiescent flux level decreasing from ~2 to <1 x 10^{-11} ergs/cm2/s over 6.5 years after outburst. Detection of pulsations during half of the quiescent observations would imply that accretion onto the magnetic poles of the neutron star continues despite the fact that the circumstellar disk may no longer be present. The accretion could come from material built-up at the corotation radius or from an isotropic stellar wind.

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 present an analysis of the neutron star High Mass X-ray Binary (HMXB) 4U 1909+07 mainly based on Suzaku data. We extend the pulse period evolution, which behaves in a random-walk like manner, indicative of direct wind accretion. Studying the spectral properties of 4U 1909+07 between 0.5 to 90 keV we find that a power-law with an exponential cutoff can describe the data well, when additionally allowing for a blackbody or a partially covering absorber at low energies. We find no evidence for a cyclotron resonant scattering feature (CRSF), a feature seen in many other neutron star HMXBs sources. By performing pulse phase resolved spectroscopy we investigate the origin of the strong energy dependence of the pulse profile, which evolves from a broad two-peak profile at low energies to a profile with a single, narrow peak at energies above 20 keV. Our data show that it is very likely that a higher folding energy in the high energy peak is responsible for this behavior. This in turn leads to the assumption that we observe the two magnetic poles and their respective accretion columns at different phases, and that these accretions column have slightly different physical conditions.

We present an in-depth study of the High Mass X-ray Binary (HMXB) GX 301-2 during its pre-periastron flare using data from the XMM-Newton satellite. The energy spectrum shows a power law continuum absorbed by a large equivalent hydrogen column on the order of 10^24 cm^2 and a prominent Fe K-alpha fluorescent emission line. Besides the Fe K-alpha line, evidence for Fe K-beta, Ni K-alpha, Ni K-beta, S K-alpha, Ar K-alpha, Ca K-alpha, and Cr K-alpha fluorescent lines is found. The observed line strengths are consistent with fluorescence in a cold absorber. This is the first time that a Cr K-alpha line is seen in emission in the X-ray spectrum of a HMXB. In addition to the modulation by the strong pulse period of ~685 sec the source is highly variable and shows different states of activity. We perform time-resolved as well as pulse-to-pulse resolved spectroscopy to investigate differences between these states of activity. We find that the fluorescent line fluxes are strongly variable and generally follow the overall flux. The N_H value is variable by a factor of 2, but not correlated to the continuum normalization. We find an interval of low flux in the light curve in which the pulsations cease almost completely, without any indication of an increasing absorption column. We investigate this dip in detail and argue that it is most likely that during the dip the accretion ceased and the afterglow of the fluorescent iron accounted for the main portion of the X-ray flux. A similar dip was found earlier in RXTE data, and we compare our findings to these results

We present the analysis of two Suzaku observations of GX 301-2 at two orbital phases after the periastron passage. Variations in the column density of the line-of-sight absorber are observed, consistent with accretion from a clumpy wind. In addition to a CRSF, multiple fluorescence emission lines were detected in both observations. The variations in the pulse profiles and the CRSF throughout the pulse phase have a signature of a magnetic dipole field. Using a simple dipole model we calculated the expected magnetic field values for different pulse phases and were able to extract a set of geometrical angles, loosely constraining the dipole geometry in the neutron star. From the variation of the CRSF width and energy, we found a geometrical solution for the dipole, making the inclination consistent with previously published values.

We present an analysis of the highly variable accreting X-ray pulsar 3A 1954+319 using 2005-2009 monitoring data obtained with INTEGRAL and Swift. This considerably extends the pulse period history and covers flaring episodes in 2005 and 2008. In 2006 the source was identified as one of only a few known symbiotic X-ray binaries (SyXBs), i.e., systems composed of a neutron star accreting from the inhomogeneous medium around an M-giant star. The extremely long pulse period of 5.3 hr is directly visible in the 2008 INTEGRAL-ISGRI outburst light curve. The pulse profile is double peaked and generally not significantly energy dependent although there is an indication of possible softening during the main pulse. During the outburst a strong spin-up of -1.8 10^(-4) hr hr^(-1) occurred. Between 2005 and 2008 a long-term spin-down trend of 2.1 10^-5 hr hr^(-1) was observed for the first time for this source. The 3-80 keV pulse peak spectrum of 3A 1954+319 during the 2008 flare could be well described by a thermal Comptonization model. We interpret the results within the framework of a recently developed quasi-spherical accretion model for SyXBs.

We present broad band analysis of the Be/X-ray transient 1A 1118-61 by Suzaku at the peak of its 3rd observed outburst in January 2009 and 2 weeks later when the source flux had decayed by an order of magnitude. The continuum was modeled with a \texttt{cutoffpl} model as well as a compTT model, with both cases requiring an additional black body component at lower energies. We confirm the detection of a cyclotron line at ~5 keV and discuss the possibility of a first harmonic at ~110 keV. Pulse profile comparisons show a change in the profile structure at lower energies, an indication for possible changes in the accretion geometry. Phase resolved spectroscopy in the outburst data show a change in the continuum throughout the pulse period. The decrease in the CRSF centroid energy also indicates that the viewing angle on the accretion column is changing throughout the pulse period.

We present the first detailed spectral and timing analysis of the High Mass X-ray Binary (HMXB) 4U 1909+07 with INTEGRAL and RXTE. 4U 1909+07 is detected in the ISGRI 20-40 keV energy band with an average countrate of 2.6 cps. The pulse period of ~604 sec is not stable, but changing erratically on timescales of years. The pulse profile is strongly energy dependent: it shows a double peaked structure at low energies, the secondary pulse decreases rapidly with increasing energy and above 20 keV only the primary pulse is visible. This evolution is consistent between PCA, HEXTE, and ISGRI. The phase averaged spectrum can be well described by the sum of a photoabsorbed power law with a cutoff at high energies and a blackbody component. To investigate the pulse profile, we performed phase resolved spectral analysis. We find that the changing spectrum can be best described with a variation of the folding energy. We rule out a correlation between the black body component and the continuum variation and discuss possible accretion geometries.

We present the first detailed spectral and timing analysis of the High Mass X-ray Binary (HMXB) 4U 1909+07 with INTEGRAL and RXTE. 4U 1909+07 is detected with an average of 2.4cps in ISGRI, but shows flares up to ~50cps. The system shows a pulse period of 605s, but we found that the period changes erratically around this value. The pulse profile is extremely energy dependent: while it shows a double peaked structure at low energies, the secondary pulse decreases rapidly with increasing energy and above 20keV only the primary pulse is visible. This evolution is consistent between PCA, HEXTE and ISGRI. We find that the phase averaged spectrum can be well fitted with a photoabsorbed power law with a cutoff at high energies and a blackbody component. To investigate the peculiar pulse profile, we performed phase resolved spectral analysis. We find that a change in the cutoff energy is required to fit the changing spectrum of the different pulse phases.