Search

Strong stellar winds are an important feature in wind-accreting high-mass X-ray binary (HMXB) systems, providing insights into stellar evolution and their impact on surrounding environments. However, the long-term evolution and temporal variability of these winds are not fully understood. This work probes the archetypal wind-accreting HMXB Vela X-1 using MAXI observations over 14 years, focusing on orbit-to-orbit absorption variability in the 2-10 keV band. Additionally, the relation between hardness ratio trends in binary orbits and neutron star spin states is investigated. We calculate hardness ratios to track absorption variability, comparing flux changes across energy bands, as the effect of absorption on the flux is energy-dependent. Variability is analyzed by comparing hardness ratio trends across binary orbits to the MAXI long-term averaged evolution. The long-term averaged hardness ratio evolution displays a stable pattern. Yet, individual binary orbits reveal different hardness ratio evolutions between consecutive orbits with no evident periodicity. Less than half of the binary orbits align with the long-term evolution. Moreover, neutron star spin-up episodes exhibit harder-than-average hardness trends compared to spin-down episodes, although their distributions overlap considerably. The long-term averaged hardness ratio dispersion is consistent with absorption column densities reported in literature from shorter observations, suggesting that heterogeneous wind structures, including accretion wakes and wind clumps, drive observed variations. The orbit-to-orbit variability indicates that pointed X-ray observations provide limited insight into wind structure. The link between neutron star spin states and hardness trends underscores the influence of accretion on absorption, with variability tied to stellar wind density fluctuations.

In 2021, the high-mass X-ray binary EXO 2030+375 underwent a giant X-ray outburst, the first since 2006, that reached a peak flux of ${\sim}600\,\mathrm{mCrab}$ (3-50\,keV). The goal of this work is to study the spectral evolution over the course of the outburst, search for possible cyclotron resonance scattering features (CRSFs), and to associate spectral components with the emission pattern of the accretion column. We used broadband spectra taken with the Nuclear Spectroscopic Telescope Array (NuSTAR), the Neutron Star Interior Composition Explorer (NICER), and Chandra near the peak and during the decline phase of the outburst. We describe the data with established empirical continuum models and perform pulse-phase-resolved spectroscopy. We compare the spectral evolution with pulse phase using a proposed geometrical emission model. We find a significant spectral hardening toward lower luminosity, a behavior that is expected for super-critical sources. The continuum shape and evolution cannot be described by a simple power-law model with exponential cutoff; it requires additional absorption or emission components. We can confirm the presence of a narrow absorption feature at ${\sim}10\,\mathrm{keV}$ in both NuSTAR observations. The absence of harmonics puts into question the interpretation of this feature as a CRSF. The empirical spectral components cannot be directly associated with identified emission components from the accretion column., Comment: 12 pages, 8 figures, accepted in A&A

From late September 2017 to February 2018, the Be X-ray binary (BeXB) Swift J0243.6+6124 underwent an unprecedently bright giant outburst. The reported X-ray luminosities were so high that the system was classified as an Ultraluminous X-ray source (ULX). It was also the first BeXB pulsar showing radio jet emission. The source was not only bright in X-rays and radio, but also in optical and UV wavelenghts. In this paper we aim to understand the origin of the observed optical/UV fluxes simultaneous to the X-ray emission. We studied the optical/UV light curves in comparison with the X-ray fluxes along the outburst, considering the main mechanisms that can explain the optical/UV emission in X-ray binaries. Due to the tight correlation observed between the optical/UV and X-ray light curves, reprocessing of X-rays seems to be the most plausible explanation. We calculated the timescales of the light curves decays and studied the correlation indexes between the optical and X-ray emission. Finally, we built a physical model considering X-ray heating of the surface of the donor star, irradiation of the accretion disk, and emission from a viscously heated accretion disk, in order to reproduce the observed optical/UV SEDs along the outburst. We considered in our model that the Be circumstellar disk was co-planar to the orbit, and then we neglected its irradiation in the current model. As an input of the model, we used as incident X-ray luminosities those calculated from the bolometric X-ray fluxes obtained from the spectral fit of the Swift/XRT and BAT observations. We conclude that reprocessing of X-rays as X-ray heating of the Be star surface and irradiation of the accretion disk are the two main mechanisms that can reproduce the observed optical/UV emission during the 2017-2018 giant outburst of Swift J0243.6+6124., Comment: Accepted for publication in A&A on 7 December 2023; 16 pages, 12 figures, 2 tables

Accreting neutron stars (NSs) represent a unique laboratory for probing the physics of accretion in the presence of strong magnetic fields ($B\gtrsim 10^8$ G). Additionally, the matter inside the NS itself exists in an ultra-dense, cold state that cannot be reproduced in Earth-based laboratories. Hence, observational studies of these objects are a way to probe the most extreme physical regimes. Here we present an overview of the field and discuss the most important outstanding problems related to NS accretion. We show how these open questions regarding accreting NSs in both low-mass and high-mass X-ray binary systems can be addressed with the High-Energy X-ray Probe (HEX-P) via simulated data. In particular, with the broad X-ray passband and improved sensitivity afforded by a low X-ray background, HEX-P will be able to 1) distinguish between competing continuum emission models; 2) provide tighter upper limits on NS radii via reflection modeling techniques that are independent and complementary to other existing methods; 3) constrain magnetic field geometry, plasma parameters, and accretion column emission patterns by characterizing fundamental and harmonic cyclotron lines and exploring their behavior with pulse phase; 4) directly measure the surface magnetic field strength of highly magnetized NSs at the lowest accretion luminosities; as well as 5) detect cyclotron line features in extragalactic sources and probe their dependence on luminosity in the super-Eddington regime in order to distinguish between geometrical evolution and accretion-induced decay of the magnetic field. In these ways HEX-P will provide an essential new tool for exploring the physics of NSs, their magnetic fields, and the physics of extreme accretion., Comment: 37 pages, 13 figures, 4 tables, To be published in a special topical issue by Frontiers in Astronomy and Space Sciences on probe-class mission concept HEX-P

LS V +44 17 is a persistent Be/X-ray binary (BeXRB) that displayed a bright, double-peaked period of X-ray activity in late 2022/early 2023. We present a radio monitoring campaign of this outburst using the Very Large Array. Radio emission was detected, but only during the second, X-ray brightest, peak, where the radio emission followed the rise and decay of the X-ray outburst. LS V +44 17 is therefore the third neutron star BeXRB with a radio counterpart. Similar to the other two systems (Swift J0243.6+6124 and 1A 0535+262), its X-ray and radio luminosity are correlated: we measure a power law slope $\beta = 1.25^{+0.64}_{-0.30}$ and a radio luminosity of $L_R = (1.6\pm0.2)\times10^{26}$ erg/s at a $0.5-10$ keV X-ray luminosity of $2\times10^{36}$ erg/s (i.e. $\sim 1\%$ $L_{\rm Edd}$). This correlation index is slightly steeper than measured for the other two sources, while its radio luminosity is higher. We discuss the origin of the radio emission, specifically in the context of jet launching. The enhanced radio brightness compared to the other two BeXRBs is the first evidence of scatter in the giant BeXRB outburst X-ray--radio correlation, similar to the scatter observed in sub-classes of low-mass X-ray binaries. While a universal explanation for such scatter is not known, we explore several options: we conclude that the three sources do not follow proposed scalings between jet power and neutron star spin or magnetic field, and instead briefly explore the effects that ambient stellar wind density may have on BeXRB jet luminosity., Comment: Accepted for publication in MNRAS on 31 October 2023. 12 pages, 3 tables, 2 figures

Neutron stars accreting from OB supergiants are often divided between persistently and transiently accreting systems, called Supergiant X-ray Binaries (SgXBs) and Supergiant Fast X-ray Transients (SFXTs). This dichotomy in accretion behaviour is typically attributed to systematic differences in the massive stellar wind, binary orbit, or magnetic field configuration, but direct observational evidence for these hypotheses remains sparse. To investigate their stellar winds, we present the results of pilot 100-GHz observations of one SFXT and one SgXB with the Northern Extended Millimetre Array. The SFXT, IGR J18410-0535, is detected as a point source at $63.4 \pm 9.6$ $\mu$Jy, while the SgXB, IGR J18410-0535 remains undetected. Radio observations of IGR J18410-0535 imply a flat or inverted low-frequency spectrum, arguing for wind emission and against non-thermal flaring. Due to the uncertain SFXT distance, however, the observations do not necessarily imply a difference between the wind properties of the SFXT and SgXB. We compare the mm constraints with other HMXBs and isolated OB supergiants, before considering how future mm campaigns can constrain HMXB wind properties by including X-ray measurements. Specifically, we discuss caveats and future steps to successfully measure wind mass loss rates and velocities in HMXBs with coordinated mm, radio, and X-ray campaigns., Comment: Accepted for publication in MNRAS Letters. 7 pages, 1 figure

The accretion of matter onto black holes and neutron stars often leads to the launching of outflows that can greatly affect the environments surrounding the compact object. In supermassive black holes, these outflows can even be powerful enough to dictate the evolution of the entire host galaxy, and yet, to date, we do not understand how these so-called accretion disk winds are launched - whether by radiation pressure, magnetic forces, thermal irradiation, or a combination thereof. An important means of studying disk winds produced near the central compact object is through X-ray absorption line spectroscopy, which allows us to probe outflow properties along a single line of sight, but usually provides little information about the global 3D disk wind structure that is vital for understanding the launching mechanism and total wind energy budget. Here, we study Hercules X-1, a unique, nearly edge-on X-ray binary with a warped accretion disk precessing with a period of about 35 days. This disk precession results in changing sightlines towards the neutron star, through the ionized outflow. We perform time-resolved X-ray spectroscopy over the precession phase and detect a strong decrease in the wind column density by three orders of magnitude as our sightline progressively samples the wind at greater heights above the accretion disk. The wind becomes clumpier as it rises upwards and expands away from the neutron star. Modelling the warped disk shape, we create a 2D map of wind properties. This unique measurement of the vertical structure of an accretion disk wind allows direct comparisons to 3D global simulations to reveal the outflow launching mechanism., Comment: Published in Nature Astronomy on April 10, available at https://www.nature.com/articles/s41550-023-01929-7

High-Mass X-ray Binaries (HMXBs) offer a unique opportunity for the investigation of accretion onto compact objects and of wind structure in massive stars. A key source for such studies is the bright neutron star HMXB Vela X-1 whose convenient physical and orbital parameters facilitate the analysis and in particular enable studies of the wind structure in HMXBs. Here, we analyse simultaneous XMM-Newton and NuSTAR observations at $\phi_{\mathrm{orb}} \approx$ 0.36-0.52 and perform time-resolved spectral analysis down to the pulse period of the neutron star, based on our previous NuSTAR-only results. For the first time, we are able to trace the onset of the wakes in a broad 0.5-78 keV range with a high-time resolution of $\sim$283 s and compare to theoretical predictions. We observe a clear rise of the absorption column density of the stellar wind $N_{\mathrm{H,1}}$ starting at orbital phase $\sim$0.44, corresponding to the wake structure entering our line of sight towards the neutron star, together with local extrema throughout the observation possibly associated with clumps or other structures in the wind. Periods of high absorption reveal the presence of multiple fluorescent emission lines of highly ionised species, mainly in the soft X-ray band between 0.5 and 4 keV, indicating photoionisation of the wind., Comment: 18 pages, 18 figures

NGC 5907 ULX1 is the most luminous ultra-luminous X-ray pulsar (ULXP) known to date, reaching luminosities in excess of 1e41 erg/s. The pulsar is known for its fast spin-up during the on-state. Here, we present a long-term monitoring of the X-ray flux and the pulse period between 2003-2022. We find that the source was in an off- or low-state between mid-2017 to mid-2020. During this state, our pulse period monitoring shows that the source had spun down considerably. We interpret this spin-down as likely being due to the propeller effect, whereby accretion onto the neutron star surface is inhibited. Using state-of-the-art accretion and torque models, we use the spin-up and spin-down episodes to constrain the magnetic field. For the spin-up episode, we find solutions for magnetic field strengths of either around 1e12G or 1e13G, however, the strong spin-down during the off-state seems only to be consistent with a very high magnetic field, namely, >1e13G. This is the first time a strong spin-down is seen during a low flux state in a ULXP. Based on the assumption that the source entered the propeller regime, this gives us the best estimate so far for the magnetic field of NGC 5907 ULX1., Comment: 10 pages, 1 figure, accepted for publication in A&A

Despite two decades of studies, it is still not clear whether ULX spectral transitions are due to stochastic variability in the wind or variations in the accretion rate or in the source geometry. The compact object is also unknown for most ULXs. In order to place constraints onto such scenarios and on the structure of the accretion disc, we studied the temporal evolution of the spectral components of the variable source NGC 55 ULX-1. Using recent and archival data obtained with the XMM-Newton satellite, we modelled the spectra with two blackbody components which we interpret as thermal emission from the inner accretion flow and the regions around or beyond the spherization radius. The luminosity-temperature (L-T) relation of each spectral component agrees with the L proportional T^4 relationship expected from a thin disc model, which suggests that the accretion rate is close to the Eddington limit. However, there are some small deviations at the highest luminosities, possibly due to an expansion of the disc and a contribution from the wind at higher accretion rates. Assuming that such deviations are due to the crossing of the Eddington or supercritical accretion rate, we estimate a compact object mass of 6-14 Msun, favouring a stellar-mass black hole as the accretor., Comment: 12 pages, 9 figures, accepted for publication in MNRAS

Because of its complex clumpy wind, prominent cyclotron resonant scattering features, intrinsic variability and convenient physical parameters (close distance, high inclination, small orbital separation) which facilitate the observation and analysis of the system, Vela X-1 is one of the key systems to understand accretion processes in high-mass X-ray binaries on all scales. We revisit Vela X-1 with two new observations taken with NuSTAR at orbital phases ~0.68-0.78 and ~0.36-0.52 which show a plethora of variability and allow us to study the accretion geometry and stellar wind properties of the system. We follow the evolution of spectral parameters down to the pulse period time-scale using a partially covered powerlaw continuum with a Fermi-Dirac cut-off to model the continuum and local absorption. We could confirm anti-correlations between the photon index and the luminosity and, for low fluxes, between the folding energy and the luminosity, implying a change of properties in the Comptonising plasma. We could not confirm a previously seen correlation between the cyclotron line energy and the luminosity of the source in the overall observation, but we observed a drop in the cyclotron line energy following a strong flare. We see strong variability in absorption between the two observations and within one observation (for the ~0.36-0.52 orbital phases) that can be explained by the presence of a large-scale structure, such as accretion- and photoionisation wakes in the system and our variable line of sight through this structure., Comment: 18 pages, 13 figures, A&A accepted

Soft Ultra-Luminous X-ray (ULXs) sources are a subclass of the ULXs that can switch from a supersoft spectral state, where most of the luminosity is emitted below 1 keV, to a soft spectral state with significant emission above 1 keV. In a few systems, dips have been observed. The mechanism behind this state transition and the dips nature are still debated. To investigate these issues, we obtained a long XMM-Newton monitoring campaign of a member of this class, NGC 247 ULX-1. We computed the hardness-intensity diagram for the whole dataset and identified two different branches: the normal branch and the dipping branch, which we study with four and three hardness-intensity resolved spectra, respectively. All seven spectra are well described by two thermal components: a colder ($kT_{\rm bb}$ $\sim$ 0.1-0.2 keV) black-body, interpreted as emission from the photo-sphere of a radiatively-driven wind, and a hotter ($kT_{\rm disk}$ $\sim$ 0.6 keV) multicolour disk black-body, likely due to reprocessing of radiation emitted from the innermost regions. In addition, a complex pattern of emission and absorption lines has been taken into account based on previous high-resolution spectroscopic results. We studied the evolution of spectral parameters and the flux of the two thermal components along the two branches and discuss two scenarios possibly connecting the state transition and the dipping phenomenon. One is based on geometrical occultation of the emitting regions, the other invokes the onset of a propeller effect., Comment: 13 pages, accepted by MNRAS on 2021/08/16

It is thought that ultraluminous X-ray sources (ULXs) are mainly powered by super-Eddington accreting neutron stars or black holes as shown by recent discovery of X-ray pulsations and relativistic winds. This work presents a follow up study of the spectral evolution over two decades of the pulsing ULX NGC 1313 X-2, in order to understand the structure of the accretion disc. The primary objective is to determine the shape and nature of the dominant spectral components by investigating their variability with the changes in the source luminosity. We have performed a spectral analysis over the canonical 0.3-10 keV energy band of all the high signal-to-noise XMM-Newton observations, and we have tested a number of different spectral models, which should approximate super-Eddington accretion discs. The baseline model consists of two thermal blackbody components with different temperatures plus an exponential cutoff powerlaw. In particular, the hotter and brighter thermal component describes the emission from the super-Eddington inner disc and the cutoff powerlaw the contribution from the accretion column of the neutron star. Instead, the cooler component describes the emission from the outer region of the disc close to the spherisation radius and the wind. The luminosity-temperature relation for the cool component follows a negative trend, which is not consistent with L$\propto$T$^4$, as expected from a sub-Eddington thin disc of Shakura-Sunayev, nor with L$\propto$T$^2$, as expected for advection-dominated disc, but would rather agree with a wind-dominated X-ray emitting region. Instead, the (L,T) relation for the hotter component is somewhere in between the first two theoretical scenarios. Our findings agree with the super-Eddington scenario and provide further detail on the disc structure. The source spectral evolution is qualitatively similar to that seen in NGC1313 X-1 and HolmbergIX X-1., Comment: 19 pages, including 10 figures and 8 tables, accepted for pubblication in A&A

Ultra-luminous X-ray pulsars (ULXPs) provide a unique opportunity to study super-Eddington accretion. We present the results of a monitoring campaign of ULXP NGC 7793 P13. Over our four-year monitoring campaign with Swift, XMM-Newton, and NuSTAR, we measured a continuous spin-up with $\dot P$ ~ -3.8e-11 s/s. The strength of the spin-up is independent of the observed X-ray flux, indicating that despite a drop in observed flux in 2019, accretion onto the source has continued at largely similar rates. The source entered an apparent off-state in early 2020, which might have resulted in a change in the accretion geometry as no pulsations were found in observations in July and August 2020. We used the long-term monitoring to update the orbital ephemeris and the periodicities seen in both the observed optical/UV and X-ray fluxes. We find that the optical/UV period is very stable over the years, with $P_\text{UV}$ = 63.75 (+0.17, -0.12) d. The best-fit orbital period determined from our X-ray timing results is 64.86 +/- 0.19 d, which is almost a day longer than previously implied, and the X-ray flux period is 65.21+/- 0.15 d, which is slightly shorter than previously measured. The physical origin of these different flux periods is currently unknown. We study the hardness ratio to search for indications of spectral changes. We find that the hardness ratios at high energies are very stable and not directly correlated with the observed flux. At lower energies we observe a small hardening with increased flux, which might indicate increased obscuration through outflows at higher luminosities. We find that the pulsed fraction is significantly higher at low fluxes. This seems to imply that the accretion geometry already changed before the source entered the deep off-state. We discuss possible scenarios to explain this behavior, which is likely driven by a precessing accretion disk., Comment: 10 pages, 6 figures, accepted for publication in A&A

Most ultraluminous X-ray sources (ULXs) are believed to be stellar mass black holes or neutron stars accreting beyond the Eddington limit. Determining the nature of the compact object and the accretion mode from broadband spectroscopy is currently a challenge, but the observed timing properties provide insight into the compact object and details of the geometry and accretion processes. Here we report a timing analysis for an 800 ks XMM-Newton campaign on the supersoft ultraluminous X-ray source, NGC 247 ULX-1. Deep and frequent dips occur in the X-ray light curve, with the amplitude increasing with increasing energy band. Power spectra and coherence analysis reveals the dipping preferentially occurs on $\sim 5$ ks and $\sim 10$ ks timescales. The dips can be caused by either the occultation of the central X-ray source by an optically thick structure, such as warping of the accretion disc, or from obscuration by a wind launched from the accretion disc, or both. This behaviour supports the idea that supersoft ULXs are viewed close to edge-on to the accretion disc., Comment: Accepted MNRAS version following minor comments. 8 pages, 6 figures

Most ULXs are believed to be powered by super-Eddington accreting neutron stars and, perhaps, black holes. Above the Eddington rate the disc is expected to thicken and to launch powerful winds through radiation pressure. Winds have been recently discovered in several ULXs. However, it is yet unclear whether the thickening of the disc or the wind variability causes the switch between the classical soft and supersoft states observed in some ULXs. In order to understand such phenomenology and the overall super-Eddington mechanism, we undertook a large (800 ks) observing campaign with XMM-Newton to study NGC 247 ULX-1, which shifts between a supersoft and classical soft ULX state. The new observations show unambiguous evidence of a wind in the form of emission and absorption lines from highly-ionised ionic species, with the latter indicating a mildly-relativistic outflow (-0.17c) in line with the detections in other ULXs. Strong dipping activity is observed in the lightcurve and primarily during the brightest observations, which is typical among soft ULXs, and indicates a close relationship between the accretion rate and the appearance of the dips. The latter is likely due to a thickening of the disc scale-height and the wind as shown by a progressively increasing blueshift in the spectral lines., Comment: 18 pages, 16 figures, accepted for publication in MNRAS

Recent observations of X-ray pulsars at low luminosities allow, for the first time, to compare theoretical models for the emission from highly magnetized neutron star atmospheres at low mass accretion rates ($\dot{M} \lesssim 10^{15}$ g s$^{-1}$) with the broadband X-ray data. The purpose of this paper is to investigate the spectral formation in the neutron star atmosphere at low $\dot{M}$ and to conduct a parameter study of physical properties of the emitting region. We obtain the structure of the static atmosphere, assuming that Coulomb collisions are the dominant deceleration process. The upper part of the atmosphere is strongly heated by the braking plasma, reaching temperatures of 30-40 keV, while its denser isothermal interior is much cooler (~2 keV). We numerically solve the polarized radiative transfer in the atmosphere with magnetic Compton scattering, free-free processes, and non-thermal cyclotron emission due to possible collisional excitations of electrons. The strongly polarized emitted spectrum has a double-hump shape that is observed in low-luminosity X-ray pulsars. A low-energy "thermal" component is dominated by extraordinary photons that can leave the atmosphere from deeper layers due to their long mean free path at soft energies. We find that a high-energy component is formed due to resonant Comptonization in the heated non-isothermal part of the atmosphere even in the absence of collisional excitations. The latter, however, affect the ratio of the two components. A strong cyclotron line originates from the optically thin, uppermost zone. A fit of the model to NuSTAR and Swift/XRT observations of GX 304-1 provides an accurate description of the data with reasonable parameters. The model can thus reproduce the characteristic double-hump spectrum observed in low-luminosity X-ray pulsars and provides insights into spectral formation., Comment: 18 pages, 10 figures, A&A accepted

We present an analysis of the first observation of the iconic High Mass X-ray Binary \so with the \chandra High Energy Transmission Gratings during an X-ray eclipse. The goal of the observation was to study the structure/physical conditions in the clumpy stellar wind through high resolution spectroscopy. We find that: a) emission line brightness from K shell transitions, corresponding to near neutral species, directly correlates with continuum illumination. However, these lines do not greatly diminish during eclipse. This is readily explained if fluorescence K$\alpha$ emission comes from the bulk of the wind. b) The highly ionised Fexxv and Fexxvi Ly$\alpha$ diminish during eclipse. Thus, they must be produced in the vicinity of the compact object where $\log \xi >3$. c) to describe the emission line spectrum, the sum of two self consistent photo ionisation models with low ionisation ($\log \xi\sim -1$) and high ionisation ($\log \xi\sim 2.4$) is required. From their emission measures, the clump-to-interclump density ratio can be estimated to be $n_c/n_i\sim 300$. To fit the complex He-like \ion{Si}{xiii}{} profile, the plasma requires a broadening with $v_{\rm bulk}\sim 840$ km s$^{-1}$. Reproducing the observed $r\approx f$ line fluxes requires the addition of a third collisionally ionised plasma. d) Emission lines widths appear unresolved at the \textsc{hetg} gratings resolution with exception of Silicon. There is no clear radial segregation between (quasi)neutral and ionised species, consistent with cold wind clumps interspersed in a hot rarefied interclump medium.

We present an improved model for excess variance spectra describing ultra-fast outflows and successfully apply it to the luminous (L ~ 10^47 erg/s) low-redshift (z = 0.184) quasar PDS 456. The model is able to account well for the broadening of the spike-like features of these outflows in the excess variance spectrum of PDS 456, by considering two effects: a correlation between the outflow velocity and the logarithmic X-ray flux and intrinsic Doppler broadening with v_int = 10^4 km/s. The models were generated by calculating the fractional excess variance of count spectra from a Monte Carlo simulation. We find evidence that the outflow in PDS 456 is structured, i.e., that there exist two or more layers with outflow velocities 0.27-0.30 c, 0.41-0.49 c, and 0.15-0.20 c for a possible third layer, which agrees well with the literature. We discuss the prospects of generally applicable models for excess variance spectra for detecting ultra-fast outflows and investigating their structure. We provide an estimate for the strength of the correlation between the outflow velocity and the logarithmic X-ray flux and investigate its validity.

We present the spectral analysis of the LMXB GRO J1744-28 using $\sim$29 ks of NuSTAR data taken in 2017 February at a low luminosity of $3.2\times 10^{36}$ erg/s (3-50 keV). The continuum spectrum is modeled with an absorbed power-law with exponential cut-off, and an additional iron line component. We find no obvious indications for a CRSF and therefore perform a detailed cyclotron line search using statistical methods on the pulse phase-averaged as well as phase-resolved spectra. The previously detected Type II X-ray bursts are absent. Clear pulsations at a period of 2.141124(9) Hz are detected. The pulse profile shows an indication of a secondary peak, which was not seen at higher flux. The 4$\sigma$ upper limit for the strength of a CRSF in the 3-20 keV band is 0.07 keV, lower than the strength of the line found at higher luminosity. The detection of pulsations shows that the source did not enter the "propeller" regime, even though the source flux of $4.15\times 10^{-10}$ erg/cm$^{2}$/s was almost one order of magnitude below the threshold for the propeller regime claimed in previous studies on this source. The transition into the propeller regime in GRO J1744-28 must therefore be below a luminosity of $3.2\times 10^{36}$ erg/s, which implies a surface magnetic field $\lesssim 2.9\times 10^{11}$ G and mass accretion rate $\lesssim 1.7\times 10^{16}$ g/s. A change of the CRSF depth as function of luminosity is not unexpected and has been observed in other sources. This result possibly implies a change in emission geometry as function of mass accretion rate to reduce the depth of the line below our detection limit.

We summarize the results of a dedicated effort between 2012 and 2019 to follow the evolution of the cyclotron line in Her~X-1 through repeated NuSTAR observations. The previously observed nearly 20-year long decay of the cyclotron line energy has ended around 2012: from there onward the pulse phase averaged flux corrected cyclotron line energy has remained stable and constant at an average value of Ecyc= (37.44+/-0.07) keV (normalized to a flux level of 6.8 RXTE/ASM-cts/s). The flux dependence of Ecyc discovered in 2007 is now measured with high precision, giving a slope of (0.675+/-0.075) keV/(ASM-cts/s), corresponding to an increase of 6.5% of Ecyc for an increase in flux by a factor of two. We also find that all line parameters as well as the continuum parameters show a correlation with X-ray flux. While a correlation between Ecyc and X-ray flux (both positive and negative) is now known for several accreting binaries with various suggestions for the underlying physics, the phenomenon of a long-term decay has so far only been seen in Her~X-1 and Vela~X-1, with far less convincing explanations., Comment: 13 pages, 17 figures, Astronomy and Astrophysics (in press)

We present a $\sim$130 ks observation of the prototypical wind-accreting, high-mass X-ray binary Vela X-1 collected with XMM-Newton at orbital phases between 0.12 and 0.28. A strong flare took place during the observation that allows us to investigate the reaction of the clumpy stellar wind to the increased X-ray irradiation. To examine the wind's reaction to the flare, we performed both time-averaged and time-resolved analyses of the RGS spectrum and examined potential spectral changes. We focused on the high-resolution XMM-Newton RGS spectra and divided the observation into pre-flare, flare, and post-flare phases. We modeled the time-averaged and time-resolved spectra with phenomenological components and with the self-consistent photoionization models calculated via CLOUDY and XSTAR in the pre-flare phase, where strong emission lines due to resonant transitions of highly ionized ions are seen. In the spectra, we find emission lines corresponding to K-shell transitions in highly charged ions of oxygen, neon, magnesium, and silicon as well as radiative recombination continua (RRC) of oxygen. Additionally, we observe potential absorption lines of magnesium at a lower ionization stage and features identified as iron L lines. The CLOUDY and XSTAR photoionization models provide contradictory results, either pointing towards uncertainties in theory or possibly a more complex multi-phase plasma, or both. We are able to demonstrate the existence of a plethora of variable narrow features, including the firm detection of oxygen lines and RRC that RGS enables to observe in this source for the first time. We show that Vela X-1 is an ideal source for future high-resolution missions, such as XRISM and Athena., Comment: 21 pages, 11 figures, A&A accepted

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., Comment: Accepted on ApJ; 10 pages, 5 figures, 3 tables

We present a qualitative search for ultra-fast outflows (UFOs) in excess variance spectra of radio-quiet active galactic nuclei (AGN). We analyse 42 sources from the Tombesi et al. (2010) spectroscopic UFO detection sample, and an additional 22 different sources from the Kara et al. (2016) variability sample. A total of 58 sources have sufficient observational data from XMM-Newton EPIC-pn and variability for an excess variance spectrum to be calculated. We examine these spectra for peaks corresponding to variable blue-shifted H- and He-like ion absorption lines from UFOs. We find good evidence for such outflows in 28% of the AGN sample and weak evidence in a further 31%, meaning that $\sim$ 30-60% of the AGN sample hosts such UFOs. The mean and median blue-shifted velocity is found to be $\sim$ 0.14c and 0.12c, respectively. Current variability methods allow for a fast, model-independent determination of UFOs, however, further work needs to be undertaken to better characterize the statistical significance of the peaks in these spectra by more rigorous modelling. Detecting good evidence for variable UFO lines in a large number of sources also lays the groundwork for detailed analysis of the variability timescales of the absorbers. This will allow us to probe their densities and hence distances from the central super-massive black hole., Comment: 23 pages, 13 figures, accepted for publication in MNRAS

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., Comment: 44 pages, 13 figures, 15 tables; accepted for publication in A&A

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., Comment: 20 pages, 16 figures, 4 tables, accepted for publication in MNRAS

MAXI J1820+070 (optical counterpart ASASSN-18ey) is a black hole candidate discovered through its recent very bright outburst. The low extinction column and long duration at high flux allow detailed measurements of the accretion process to be made. In this work, we compare the evolution of X-ray spectral and timing properties through the initial hard state of the outburst. We show that the inner accretion disc, as measured by relativistic reflection, remains steady throughout this period of the outburst. Nevertheless, subtle spectral variability is observed, which is well explained by a change in coronal geometry. However, characteristic features of the temporal variability - low-frequency roll-over and QPO frequency - increase drastically in frequency, as the outburst proceeds. This suggests that the variability timescales are governed by coronal conditions rather than solely by the inner disc radius. We also find a strong correlation between X-ray luminosity and coronal temperature. This can be explained by electron pair production with a changing effective radius and a non-thermal electron fraction of ~20%., Comment: Re-submitted to MNRAS after reviewers comments. 14 pages, 14 figures

SN2010da/NGC 300 ULX-1 was first detected as a supernova impostor in May 2010 and was recently discovered to be a pulsating ultraluminous X-ray source. In this letter, we present VLT/X-shooter spectra of this source obtained in October 2018, covering the wavelength range 350-2300 nm. The $J$- and $H$-bands clearly show the presence of a red supergiant donor star that is best matched by a MARCS stellar atmosphere with $T_{\rm eff} = 3650 - 3900$ K and $\log(L_{\rm bol}/L_{\odot}) = 4.25\pm0.10$, which yields a stellar radius $R = 310 \pm 70 R_{\odot}$. To fit the full spectrum, two additional components are required: a blue excess that can be fitted either by a hot blackbody (T $\gtrsim 20,000$ K) or a power law (spectral index $\alpha \approx 4$) and is likely due to X-ray emission reprocessed in the outer accretion disk or the donor star; and a red excess that is well fitted by a blackbody with a temperature of $\sim 1100$ K, and is likely due to warm dust in the vicinity of SN2010da. The presence of a red supergiant in this system implies an orbital period of at least 0.8-2.1 years, assuming Roche lobe overflow. Given the large donor-to-compact object mass ratio, orbital modulations of the radial velocity of the red supergiant are likely undetectable. However, the radial velocity amplitude of the neutron star is large enough (up to 40-60 km s$^{-1}$) to potentially be measured in the future, unless the system is viewed at a very unfavorable inclination., Comment: 10 pages, 6 figures, accepted for publication in ApJ Letters

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

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., Comment: 6 pages, 4 figures, accepted for publication by MNRAS. Comments are welcomed

Vela X-1 is an eclipsing high mass X-ray binary (HMXB) consisting of a 283s accreting X-ray pulsar in a close orbit of 8.964 days around the B0.5Ib supergiant HD77581 at a distance of just 2.4 kpc. The system is considered a prototype of wind-accreting HMXB and it has been used as a baseline in different theoretical or modelling studies. We discuss the observational properties of the system and the use of the observational data as laboratory to test recent developments in modelling the accretion process in High-Mass X-ray Binaries (e.g., Sander et al. 2018; El Mellah et al. 2018), which range from detailed descriptions of the wind acceleration to modelling of the structure of the flow of matter close to the neutron star and its variations., Comment: 4 pages, 2 figures, proceedings of the 12th INTEGRAL conference "INTEGRAL looks AHEAD to Multimessenger astronomy" in Geneva (Switzerland) on 11-15 February 2019

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., Comment: 6 pages, 4 figures, accepted for publication in Astronomische Nachrichten (proceedings of the XMM-Newton Workshop 2019)

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

We report on our search for the optical counterparts of two ultraluminous X-ray pulsars with known orbital periods, M82 X-2 and NGC 5907 X-1, in new and archival HST observations, in an effort to characterize the donor stars in these systems. We detect five near-infrared sources consistent with the position of M82 X-2 that are too bright to be single stars. We also detect seven sources in the WFC3/UVIS F336W image whose photometry matches that of 10-15 M$_\odot$ stars turning off the main sequence. Such stars have densities consistent with the properties of the donor star of M82 X-2 as inferred from X-ray timing analysis, although it is also possible that the donor is a lower mass star below our detection limit or that there is a significant contribution from the accretion disc to the optical emission. We detect three candidate counterparts to NGC 5907 X-1 in the near-infrared. All of these are too bright to be the donor star of the ULX, which based on its orbital period is a red giant. The high background at the location of NGC 5907 X-1 precludes us from detecting this expected donor star. The recently discovered NGC 5907 ULX-2 also falls within the field of view of the near-infrared imaging; we detect four sources in the error circle, with photometry that matches AGB stars. The star suggested to be the counterpart of NGC 5907 ULX-2 by Pintore et al. (2018) falls outside our 2-$\sigma$ error circle., Comment: 10 pages, 3 figures, accepted for publication in ApJ

Cyclotron lines, also called cyclotron resonant scattering features (CRSF) are spectral features, generally appearing in absorption, in the X-ray spectra of objects containing highly magnetized neutron stars, allowing the direct measurement of the magnetic field strength in these objects. Cyclotron features are thought to be due to resonant scattering of photons by electrons in the strong magnetic fields. The main content of this contribution focusses on electron cyclotron lines as found in accreting X-ray binary pulsars (XRBP) with magnetic fields on the order of several 1012 Gauss. Also, possible proton cyclotron lines from single neutron stars with even stronger magnetic fields are briefly discussed.With regard to electron cyclotron lines, we present an updated list of XRBPs that show evidence of such absorption lines. The first such line was discovered in a 1976 balloon observation of the accreting binary pulsar Hercules X-1, it is considered to be the first direct measurement of the magnetic field of a neutron star. As of today (mid 2018), we list 36 XRBPs showing evidence of one ore more electron cyclotron absorption line(s). A few have been measured only once and must be confirmed (several more objects are listed as candidates). In addition to the Tables of objects, we summarize the evidence of variability of the cyclotron line as a function of various parameters (especially pulse phase, luminosity and time), and add a discussion of the different observed phenomena and associated attempts of theoretical modeling. We also discuss our understanding of the underlying physics of accretion onto highly magnetized neutron stars. For proton cyclotron lines, we present tables with seven neutron stars and discuss their nature and the physics in these objects., Comment: 32 pages, 15 figures, 8 Tables, accepted by A&A 2018

SS433 is the only Galactic binary system known to persistently accrete at highly super-critical (or hyper-critical) rates, similar to those in tidal disruption events, and likely needed to explain the rapid growth of those very high redshift quasars containing massive SMBHs. Probing the inner regions of SS433 in the X-rays is crucial to understanding this system, and super-critical accretion in general, but is highly challenging due to obscuration by the surrounding wind, driven from the accretion flow. NuSTAR observed SS433 in the hard X-ray band across multiple phases of its 162 day super-orbital precession period. Spectral-timing tools allow us to infer that the hard X-ray emission from the inner regions is likely being scattered towards us by the walls of the wind-cone. By comparing to numerical models, we determine an intrinsic X-ray luminosity of $\ge$ 2$\times$10$^{37}$ erg/s and that, if viewed face on, we would infer an apparent luminosity of $>$ 1$\times$10$^{39}$ erg/s, confirming SS433's long-suspected nature as an ultraluminous X-ray source (ULX). We present the discovery of a narrow, $\sim 100$s lag due to atomic processes occurring in outflowing material travelling at least 0.14-0.29c, which matches absorption lines seen in ULXs and -- in the future -- will allow us to map a super-critical outflow for the first time., Comment: 14 pages, 10 figures, accepted for publication in MNRAS

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., Comment: 18 pages, 10 figures, accepted for publication in A&A

Based on our new NuSTAR X-ray telescope data, we rule out any cyclotron line up to 60 keV in the spectra of the high mass X-ray binary 4U2206+54. In particular, we do not find any evidence of the previously claimed line around 30 keV, independently of the source flux, along the spin pulse. The spin period has increased significantly, since the last observation, up to $5750\pm10$ s, confirming the rapid spin down rate $\dot \nu=-1.8\times 10^{-14}$~Hz s$^{-1}$. This behaviour might be explained by the presence of a strongly magnetized neutron star ($B_{\rm s}>$ several times $10^{13}$ G) accreting from the slow wind of its main sequence O9.5 companion., Comment: Accepted for publication in MNRAS, Main Journal

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., Comment: Accepted on MNRAS, 5 pages, 2 figure, 1 table

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., Comment: 10 pages, 6 figures, accepted for publication in A&A

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., Comment: 6 pages, 3 figures. Accepted for publication in ApJL

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., Comment: 13 pages, 6 figures; accepted for publication in ApJ

Ultraluminous pulsars are a definite proof that persistent super-Eddington accretion occurs in nature. They support the scenario according to which most Ultraluminous X-ray Sources (ULXs) are super-Eddington accretors of stellar mass rather than sub-Eddington intermediate mass black holes. An important prediction of theories of supercritical accretion is the existence of powerful outflows of moderately ionized gas at mildly relativistic speeds. In practice, the spectral resolution of X-ray gratings such as RGS onboard XMM-Newton is required to resolve their observational signatures in ULXs. Using RGS, outflows have been discovered in the spectra of 3 ULXs (none of which are currently known to be pulsars). Most recently, the fourth ultraluminous pulsar was discovered in NGC 300. Here we report detection of an ultrafast outflow (UFO) in the X-ray spectrum of the object, with a significance of more than 3{\sigma}, during one of the two simultaneous observations of the source by XMM-Newton and NuSTAR in December 2016. The outflow has a projected velocity of 65000 km/s (0.22c) and a high ionisation factor with a log value of 3.9. This is the first direct evidence for a UFO in a neutron star ULX and also the first time that this its evidence in a ULX spectrum is seen in both soft and hard X-ray data simultaneously. We find no evidence of the UFO during the other observation of the object, which could be explained by either clumpy nature of the absorber or a slight change in our viewing angle of the accretion flow., Comment: 10 pages, 4 figures. Accepted to MNRAS

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., Comment: 14 pages, 9 figures, accepted to MNRAS

Bright and eclipsing, the high-mass X-ray binary Vela X-1 offers a unique opportunity to study accretion onto a neutron star from clumpy winds of O/B stars and to disentangle the complex accretion geometry of these systems. In Chandra-HETGS spectroscopy at orbital phase ~0.25, when our line of sight towards the source does not pass through the large-scale accretion structure such as the accretion wake, we observe changes in overall spectral shape on timescales of a few kiloseconds. This spectral variability is, at least in part, caused by changes in overall absorption and we show that such strongly variable absorption cannot be caused by unperturbed clumpy winds of O/B stars. We detect line features from high and low ionization species of silicon, magnesium and neon whose strengths and presence depend on the overall level of absorption. They imply a co-existence of cool and hot gas phases in the system that we interpret as a highly variable, structured accretion flow close to the compact object such as has been recently seen in simulations of wind accretion in high-mass X-ray binaries., Comment: 17 pages, 14 figures; A&A accepted

Recent observations of Her X-1 with NuSTAR and INTEGRAL in 2016 have provided evidence that the 20-year decay of the cyclotron line energy found between 1996 and 2015 has ended and that an inversion with a new increase, possibly similar to the one observed around 1990-1993, has started. We consider this a strong motivation for further observations and for enhanced efforts to significantly improve our theoretical understanding of the accretion process in binary X-ray pulsars. We speculate about the physics behind the long-term decay and its inversion, a possible cyclic behavior, and correlations with other variable observables., Comment: 5 pages, 1 figure, accepted by A&A

Based on new Chandra X-ray telescope data, we present empirical evidence of plasma Compton cooling during a flare in the non pulsating massive X-ray binary 4U1700-37. This behaviour might be explained by quasispherical accretion onto a slowly rotating magnetised neutron star. In quiescence, the neutron star in 4U1700-37 is surrounded by a hot radiatively cooling shell. Its presence is supported by the detection of mHz quasi periodic oscillations likely produced by its convection cells. The high plasma temperature and the relatively low X-ray luminosity observed during the quiescence, point to a small emitting area about 1 km, compatible with a hot spot on a NS surface. The sudden transition from a radiative to a significantly more efficient Compton cooling regime triggers an episode of enhanced accretion resulting in a flare. During the flare, the plasma temperature drops quickly. The predicted luminosity for such transitions, Lx = 3 x 10^35 erg s-1, is very close to the luminosity of 4U1700-37 during quiescence. The transition may be caused by the accretion of a clump in the stellar wind of the donor star. Thus, a magnetised NS nature of the compact object is strongly favoured., Comment: Accepted for publication in MNRAS

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., Comment: 52 pages, 15 figures, accepted, invited review (to appear in New Astronomy Reviews), v3: updated figure and text

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., Comment: 11 pages, 6 figures, accepted for publication in A&A (v2 including language edits)

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., Comment: 15 pages, 7 figures, accepted for publication in MNRAS