Search

We present a series of five new broadband X-ray observations of the ultraluminous X-ray source Holmberg IX X-1, performed by $XMM$-$Newton$ and $NuSTAR$ in coordination. The first three of these show high soft X-ray fluxes but a near total collapse of the high-energy ($\gtrsim$15 keV) emission, previously seen to be surprisingly stable across all prior broadband observations of the source. The latter two show a recovery in hard X-rays, remarkably once again respecting the same stable high-energy flux exhibited by all of the archival observations. We also present a joint analysis of all broadband observations of Holmberg IX X-1 to date (encompassing 11 epochs in total) in order to investigate whether it shows the same luminosity-temperature behaviour as NGC 1313 X-1 (which also shows a stable high-energy flux), whereby the hotter disc component in the spectrum exhibits two distinct, positively-correlated tracks in the luminosity-temperature plane. Holmberg IX X-1 may show similar behaviour, but the results depend on whether the highest energy emission is assumed to be an up-scattering corona or an accretion column. The strongest evidence for this behaviour is found in the former case, while in the latter the new 'soft' epochs appear distinct from the other high-flux epochs. We discuss possible explanations for these new 'soft' spectra in the context of the expected structure of super-Eddington accretion flows around black holes and neutron stars, and highlight a potentially interesting analogy with the recent destruction and re-creation of the corona seen in the AGN 1ES 1927+654., Comment: 15 pages, 5 figures, submitted for publication in MNRAS (comments welcome)

Cygnus X-1 is a high mass X-ray binary where accretion onto the black hole is mediated by the stellar wind from the blue supergiant companion star HDE 226868. Depending on the position of the black hole along the orbit, X-ray observations can probe different layers of the stellar wind. Deeper wind layers can be investigated at superior conjunction (i.e. null orbital phases). We aim at characterising the stellar wind in the Cyg X-1/HDE 226868 system analysing one passage at superior conjunction covered by XMM-Newton during the CHOCBOX campaign via modelling of colour-colour diagrams. Since X-ray absorption is energy-dependent, colour indices provide information on the parameters of the stellar wind, such as the column density $N_{H,w}$ and the covering factor $f_c$. We fitted colour-colour diagrams with models that include both a continuum and a stellar wind component. We used the KDE method to infer the unknown probability distribution of the data points in the colour-colour diagram, and selected the model corresponding to the highest likelihood. In order to study the temporal evolution of the wind around superior conjunction, we extracted and fitted time-resolved colour-colour diagrams. We found that the model that best describes the shape of the colour-colour diagram of Cyg X-1 at superior conjunction requires the wind to be partially ionised. The shape of the colour-colour diagram strongly varies during the analysed observation, as due to concurrent changes of the mean $N_{H,w}$ and the $f_c$ of the wind. Our results suggest the existence of a linear scaling between the rapid variability amplitude of $N_{H,w}$ (on time scales between 10 s and 11 ks) and its long term variations (on time scales 11>ks). Using the inferred best-fit values, we estimated the stellar mass loss rate to be $\sim 7\times10^{-6} {\rm M_{\odot}yr^{-1}}$ and the clumps to have a mass of $\sim10^{17}$ g., Comment: Accepted for publication in A&A

We used the Imaging X-ray Polarimetry Explorer (IXPE) satellite to measure, for the first time, the 2-8 keV polarization of NGC 1068. We pointed IXPE for a net exposure time of 1.15 Ms on the target, in addition to two ~ 10 ks each Chandra snapshots in order to account for the potential impact of several ultraluminous X-ray source (ULXs) within IXPE's field-of-view. We measured a 2 - 8 keV polarization degree of 12.4% +/- 3.6% and an electric vector polarization angle of 101{\deg} +/- 8{\deg} at 68% confidence level. If we exclude the spectral region containing the bright Fe K lines and other soft X-ray lines where depolarization occurs, the polarization fraction rises up to 21.3% +/- 6.7% in the 3.5 - 6.0 keV band, with a similar polarization angle. The observed polarization angle is found to be perpendicular to the parsec scale radio jet. Using a combined Chandra and IXPE analysis plus multi-wavelength constraints, we estimated that the circumnuclear "torus" may sustain a half-opening angle of 50{\deg} - 55{\deg} (from the vertical axis of the system). Thanks to IXPE, we have measured the X-ray polarization of NGC 1068 and found comparable results, both in terms of polarization angle orientation with respect to the radio-jet and torus half-opening angle, to the X-ray polarimetric measurement achieved for the other archetypal Compton-thick AGN : the Circinus galaxy. Probing the geometric arrangement of parsec-scale matter in extragalactic object is now feasible thanks to X-ray polarimetry., Comment: 9 pages, 5 figures, 2 tables, submited to A&A

One of the most important objectives of solar physics is the physical understanding of the solar atmosphere, the structure of which is also described in terms of the density (N) and temperature (T) distributions of the atmospheric matter. Several multi-frequency analyses show that the characteristics of these distributions are still debated, especially for the outer coronal emission. We aim to constrain the T and N distributions of the solar atmosphere through observations in the centimetric radio domain. We employ single-dish observations from two of the INAF radio telescopes at the K-band frequencies (18 - 26 GHz). We investigate the origin of the significant brightness temperature ($T_B$) level that we detected up to the upper corona ($\sim 800$ Mm of altitude with respect to the photospheric solar surface). To probe the physical origin of the atmospheric emission and to constrain instrumental biases, we reproduced the solar signal by convolving specific 2D antenna beam models. The analysis of the solar atmosphere is performed by adopting a physical model that assumes the thermal bremsstrahlung as the emission mechanism, with specific T and N distributions. The modelled $T_B$ profiles are compared with those observed by averaging solar maps obtained during the minimum of solar activity (2018 - 2020). The T and N distributions are compatible (within $25\%$ of uncertainty) with the model up to $\sim 60$ Mm and $\sim 100$ Mm of altitude, respectively. The analysis of the role of the antenna beam pattern on our solar maps proves the physical nature of the atmospheric emission in our images up to the coronal tails seen in our $T_B$ profiles. The challenging analysis of the coronal radio emission at higher altitudes, together with the data from satellite instruments will require further multi-frequency measurements., Comment: 12 pages, 8 figures, 5 tables, accepted by A&A; v1

The Sun is an extraordinary workbench, from which several fundamental astronomical parameters can be measured with high precision. Among these parameters, the solar radius $R_{\odot}$ plays an important role in several aspects, such as in evolutionary models. Despite the efforts in obtaining accurate measurements of $R_{\odot}$, the subject is still debated and measurements are puzzling and/or lacking in many frequency ranges. We aimed to determine the mean, equatorial, and polar radii of the Sun ($R_c$, $R_{eq}$, and $R_{pol}$) in the frequency range 18.1 - 26.1 GHz. We employed single-dish observations from the newly-appointed Medicina "Gavril Grueff" Radio Telescope and the Sardinia Radio Telescope (SRT) throughout 5 years, from 2018 to mid-2023, in the framework of the SunDish project for solar monitoring. Two methods to calculate the radius at radio frequencies are considered and compared. To assess the quality of our radius determinations, we also analysed the possible degrading effects of the antenna beam pattern on our solar maps, using two 2D-models. We carried out a correlation analysis with the evolution of the solar cycle through the calculation of Pearson's correlation coefficient $\rho$. We obtained several values for the solar radius - ranging between 959 and 994 arcsec - and $\rho$, with typical errors of a few arcsec. Our $R_{\odot}$ measurements, consistent with values reported in literature, suggest a weak prolatness of the solar limb ($R_{eq}$ > $R_{pol}$), although $R_{eq}$ and $R_{pol}$ are statistically compatible within 3$\sigma$ errors. The correlation analysis using the solar images from Grueff shows (1) a positive correlation between the solar activity and the temporal variation of $R_c$ (and $R_{eq}$) at all observing frequencies, and (2) a weak anti-correlation between the temporal variation of $R_{pol}$ and the solar activity at 25.8 GHz., Comment: 18 pages, 12 figures, 6 tables, accepted by A&A; v1

The hard X-ray emission from magnetars and other isolated neutron stars remains under-explored. An instrument with higher sensitivity to hard X-rays is critical to understanding the physics of neutron star magnetospheres and also the relationship between magnetars and Fast Radio Bursts (FRBs). High sensitivity to hard X-rays is required to determine the number of magnetars with hard X-ray tails, and to track transient non-thermal emission from these sources for years post-outburst. This sensitivity would also enable previously impossible studies of the faint non-thermal emission from middle-aged rotation-powered pulsars (RPPs), and detailed phase-resolved spectroscopic studies of younger, bright RPPs. The High Energy X-ray Probe (HEX-P) is a probe-class mission concept that will combine high spatial resolution X-ray imaging ($<5$ arcsec half-power diameter (HPD) at 0.2--25 keV) and broad spectral coverage (0.2--80 keV) with a sensitivity superior to current facilities (including XMM-Newton and NuSTAR). HEX-P has the required timing resolution to perform follow-up observations of sources identified by other facilities and positively identify candidate pulsating neutron stars. Here we discuss how HEX-P is ideally suited to address important questions about the physics of magnetars and other isolated neutron stars., Comment: 22 pages, 8 figures, submitted to Frontiers

We present the results of the first X-ray polarimetric observation of the low-mass X-ray binary 4U 1957+115, performed with the Imaging X-ray Polarimetry Explorer in May 2023. The binary system has been in a high-soft spectral state since its discovery and is thought to host a black hole. The $\sim$571 ks observation reveals a linear polarisation degree of $1.9\% \pm 0.6\%$ and a polarisation angle of $-41^\circ.8 \pm 7^\circ.9$ in the 2-8 keV energy range. Spectral modelling is consistent with the dominant contribution coming from the standard accretion disc, while polarimetric data suggest a significant role of returning radiation: photons that are bent by strong gravity effects and forced to return to the disc surface, where they can be reflected before eventually reaching the observer. In this setting, we find that models with a black hole spin lower than 0.96 and an inclination lower than $50^\circ$ are disfavoured., Comment: 12 pages, 10 figures, 2 tables, accepted for publication in A&A

We report on a comprehensive analysis of simultaneous X-ray polarimetric and spectral data of the bright atoll source GX 9+9 with the Imaging X-ray Polarimetry Explorer (IXPE) and NuSTAR. The source is significantly polarized in the 4--8 keV band, with a degree of $2.2\% \pm 0.5\%$ (uncertainty at the 68% confidence level). The NuSTAR broad-band spectrum clearly shows an iron line, and is well described by a model including thermal disk emission, a Comptonized component, and reflection. From a spectro-polarimetric fit, we obtain an upper limit to the polarization degree of the disk of 4% (at 99% confidence level), while the contribution of Comptonized and reflected radiation cannot be conclusively separated. However, the polarization is consistent with resulting from a combination of Comptonization in a boundary or spreading layer, plus reflection off the disc, which gives a significant contribution in any realistic scenario., Comment: 10 pages, accepted for publication in A&A

Aims. Fast Radio Bursts are bright radio transients whose origin has not yet explained. The search for a multi-wavelength counterpart of those events can put a tight constrain on the emission mechanism and the progenitor source. Methods. We conducted a multi-wavelength observational campaign on FRB 20180916B between October 2020 and August 2021 during eight activity cycles of the source. Observations were led in the radio band by the SRT both at 336 MHz and 1547 MHz and the uGMRT at 400 MHz. Simultaneous observations have been conducted by the optical telescopes Asiago (Galileo and Copernico), CMO SAI MSU, CAHA 2.2m, RTT-150 and TNG, and X/Gamma-ray detectors on board the AGILE, Insight-HXMT, INTEGRAL and Swift satellites. Results. We present the detection of 14 new bursts detected with the SRT at 336 MHz and seven new bursts with the uGMRT from this source. We provide the deepest prompt upper limits in the optical band fro FRB 20180916B to date. In fact, the TNG/SiFAP2 observation simultaneous to a burst detection by uGMRT gives an upper limit E_optical / E_radio < 1.3 x 10^2. Another burst detected by the SRT at 336 MHz was also co-observed by Insight-HMXT. The non-detection in the X-rays yields an upper limit (1-30 keV band) of E_X-ray / E_radio in the range of (0.9-1.3) x 10^7, depending on which model is considered for the X-ray emission., Comment: A&A accepted

We present an X-ray spectro-polarimetric analysis of the bright Seyfert galaxy IC 4329A. The Imaging X-ray Polarimetry Explorer (IXPE) observed the source for ~500 ks, supported by XMM-Newton (~60 ks) and NuSTAR (~80 ks) exposures. We detect polarisation in the 2-8 keV band with 2.97 sigma confidence. We report a polarisation degree of $3.3\pm1.1$ per cent and a polarisation angle of $78\pm10$ degrees (errors are 1 sigma confidence). The X-ray polarisation is consistent with being aligned with the radio jet, albeit partially due to large uncertainties on the radio position angle. We jointly fit the spectra from the three observatories to constrain the presence of a relativistic reflection component. From this, we obtain constraints on the inclination angle to the inner disc (< 39 degrees at 99 per cent confidence) and the disc inner radius (< 11 gravitational radii at 99 per cent confidence), although we note that modelling systematics in practice add to the quoted statistical error. Our spectro-polarimetric modelling indicates that the 2-8 keV polarisation is consistent with being dominated by emission directly observed from the X-ray corona, but the polarisation of the reflection component is completely unconstrained. Our constraints on viewer inclination and polarisation degree tentatively favour more asymmetric, possibly out-flowing, coronal geometries that produce more highly polarised emission, but the coronal geometry is unconstrained at the 3 sigma level., Comment: Accepted for publication in MNRAS, 13 pages, 9 figures

We report on the second observation of the radio-quiet active galactic nucleus (AGN) MCG-05-23-16 performed with the Imaging X-ray Polarimetry Explorer (IXPE). The observation started on 2022 November 6 for a net observing time of 640 ks, and was partly simultaneous with NuSTAR (86 ks). After combining these data with those obtained in the first IXPE pointing on May 2022 (simultaneous with XMM-Newton and NuSTAR) we find a 2-8 keV polarization degree $\Pi$ = 1.6 $\pm$ 0.7 (at 68 per cent confidence level), which corresponds to an upper limit $\Pi$ = 3.2 per cent (at 99 per cent confidence level). We then compare the polarization results with Monte Carlo simulations obtained with the MONK code, with which different coronal geometries have been explored (spherical lamppost, conical, slab and wedge). Furthermore, the allowed range of inclination angles is found for each geometry. If the best fit inclination value from a spectroscopic analysis is considered, a cone-shaped corona along the disc axis is disfavoured., Comment: 9 pages, 8 figures, 3 table. Submitted to MNRAS

X Persei is a persistent low-luminosity X-ray pulsar of period of $\sim$835 s in a Be binary system. The field strength at the neutron star surface is not known precisely, but indirect signs indicate a magnetic field above $10^{13}$ G, which makes the object one of the most magnetized known X-ray pulsars. Here we present the results of observations X Persei performed with the Imaging X-ray Polarimetry Explorer (IXPE). The X-ray polarization signal was found to be strongly dependent on the spin phase of the pulsar. The energy-averaged polarization degree in 3-8 keV band varied from several to $\sim$20 per cent over the pulse with a positive correlation with the pulsed X-ray flux. The polarization angle shows significant variation and makes two complete revolutions during the pulse period resulting in nearly nil pulse-phase averaged polarization. Applying the rotating vector model to the IXPE data we obtain the estimates for the rotation axis inclination and its position angle on the sky as well as for the magnetic obliquity. The derived inclination is close to the orbital inclination reported earlier for X Persei. The polarimetric data imply a large angle between the rotation and magnetic dipole axes, which is similar to the result reported recently for the X-ray pulsar GRO J1008$-$57. After eliminating the effect of polarization angle rotation over the pulsar phase using the best-fitting rotating vector model, the strong dependence of the polarization degree with energy was discovered with its value increasing from 0% at $\sim$2 keV to 30% at 8 keV., Comment: 11 pages, 7 figures, submitted to MNRAS

We present an X-ray spectro-polarimetric analysis of the bright Seyfert galaxy NGC4151. The source has been observed with the Imaging X-ray Polarimetry Explorer (IXPE) for 700 ks, complemented with simultaneous XMM-Newton (50 ks) and NuSTAR (100 ks) pointings. A polarization degree ${\Pi} = 4.9 {\pm} 1.1 \%$ and angle ${\Psi}= 86{\deg} {\pm} 7{\deg}$ east of north ($68\%$ confidence level) are measured in the 2-8 keV energy range. The spectro-polarimetric analysis shows that the polarization could be entirely due to reflection. Given the low reflection flux in the IXPE band, this requires however a reflection with a very large ($> 38 \%$) polarization degree. Assuming more reasonable values, a polarization degree of the hot corona ranging from ${\sim}4$ to ${\sim}8\%$ is found. The observed polarization degree excludes a spherical lamppost geometry for the corona, suggesting instead a slab-like geometry, possibly a wedge, as determined via Monte Carlo simulations. This is further confirmed by the X-ray polarization angle, which coincides with the direction of the extended radio emission in this source, supposed to match the disc axis. NGC4151 is the first AGN with an X-ray polarization measure for the corona, illustrating the capabilities of X-ray polarimetry and IXPE in unveiling its geometry., Comment: 9 pages, 4 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

We report spectro-polarimetric results of an observational campaign of the bright neutron star low-mass X-ray binary Cyg X-2 simultaneously observed by IXPE, NICER and INTEGRAL. Consistently with previous results, the broad-band spectrum is characterized by a lower-energy component, attributed to the accretion disc with $kT_{\rm in} \approx$ 1 keV, plus unsaturated Comptonization in thermal plasma with temperature $kT_{\rm e} = 3$ keV and optical depth $\tau \approx 4$, assuming a slab geometry. We measure the polarization degree in the 2-8 keV band $P=1.8 \pm 0.3$ per cent and polarization angle $\phi = 140^{\circ} \pm 4^{\circ}$, consistent with the previous X-ray polarimetric measurements by OSO-8 as well as with the direction of the radio jet which was earlier observed from the source. While polarization of the disc spectral component is poorly constrained with the IXPE data, the Comptonized emission has a polarization degree $P =4.0 \pm 0.7$ per cent and a polarization angle aligned with the radio jet. Our results strongly favour a spreading layer at the neutron star surface as the main source of the polarization signal. However, we cannot exclude a significant contribution from reflection off the accretion disc, as indicated by the presence of the iron fluorescence line., Comment: 10 pages, 7 figures, accepted for publication in MNRAS

We report on the Imaging X-ray Polarimetry Explorer (IXPE) observation of the closest and X-ray brightest Compton-thick active galactic nucleus (AGN), the Circinus galaxy. We find the source to be significantly polarized in the 2--6 keV band. From previous studies, the X-ray spectrum is known to be dominated by reflection components, both neutral (torus) and ionized (ionization cones). Our analysis indicates that the polarization degree is $28 \pm 7$ per cent (at 68 per cent confidence level) for the neutral reflector, with a polarization angle of $18{\deg} \pm 5{\deg}$, roughly perpendicular to the radio jet. The polarization of the ionized reflection is unconstrained. A comparison with Monte Carlo simulations of the polarization expected from the torus shows that the neutral reflector is consistent with being an equatorial torus with a half-opening angle of 45{\deg}-55{\deg}. This is the first X-ray polarization detection in a Seyfert galaxy, demonstrating the power of X-ray polarimetry in probing the geometry of the circumnuclear regions of AGNs, and confirming the basic predictions of standard Unification Models., Comment: 9 pages, accepted for publication in MNRAS, authors list updated

We report on the first observation of a radio-quiet Active Galactic Nucleus (AGN) using polarized X-rays: the Seyfert 1.9 galaxy MCG-05-23-16. This source was pointed with the Imaging X-ray Polarimetry Explorer (IXPE) starting on May 14, 2022 for a net observing time of 486 ks, simultaneously with XMM-Newton (58 ks) and NuSTAR (83 ks). A polarization degree smaller than $\Pi<4.7\%$ (at the 99% c.l.) is derived in the 2-8 keV energy range, where emission is dominated by the primary component ascribed to the hot corona. The broad-band spectrum, inferred from a simultaneous fit to the IXPE, NuSTAR, and XMM-Newton data, is well reproduced by a power law with photon index $\Gamma=1.85\pm0.01$ and a high-energy cutoff $E_{\rm C}=120\pm15$ keV. A comparison with Monte Carlo simulations shows that a lamp-post and a conical geometry of the corona are consistent with the observed upper limit, a slab geometry is allowed only if the inclination angle of the system is less than 50$^{\circ}$., Comment: 7 pages, 3 figures, 1 table. Submitted to MNRAS Letters

We present a new solar radio imaging system implemented through the upgrade of the large single-dish telescopes of the Italian National Institute for Astrophysics (INAF), not originally conceived for solar observations. During the development and early science phase of the project (2018-2020), we obtained about 170 maps of the entire solar disk in the 18-26 GHz band, filling the observational gap in the field of solar imaging at these frequencies. These solar images have typical resolutions in the 0.7-2 arcmin range and a brightness temperature sensitivity <10 K. Accurate calibration adopting the Supernova Remnant Cas A as a flux reference, provided typical errors <3% for the estimation of the quiet-Sun level components and for active regions flux measurements. As a first early science result of the project, we present a catalog of radio continuum solar imaging observations with Medicina 32-m and SRT 64-m radio telescopes including the multi-wavelength identification of active regions, their brightness and spectral characterization. The interpretation of the observed emission as thermal bremsstrahlung components combined with gyro-magnetic variable emission pave the way to the use of our system for long-term monitoring of the Sun. We also discuss useful outcomes both for solar physics (e.g. study of the chromospheric network dynamics) and space weather applications (e.g. flare precursors studies)., Comment: 43 pages, 11 figures, 6 tables, accepted for publication in Solar Physics

StrayCats, the catalog of NuSTAR stray light observations, contains data from bright X-ray sources that fall within crowded source regions. These observations offer unique additional data with which to monitor sources like X-ray binaries that show variable timing behavior. In this work, we present a timing analysis of stray light data of the high mass X-ray binary SMC X-1, the first scientific analysis of a single source from the StrayCats project. We describe the process of screening stray light data for scientific analysis, verify the orbital ephemeris, and create both time and energy resolved pulse profiles. We find that the orbital ephemeris of SMC X-1 is unchanged and confirm a long-term spin up rate of $\dot{\nu}=(2.52\pm0.03)\times10^{-11}$ Hz s$^{-1}$. We also note that the shape of SMC X-1's pulse profile, while remaining double-peaked, varies significantly with time and only slightly with energy., Comment: 10 pages, 8 figures

Due to be launched in late 2021, the Imaging X-Ray Polarimetry Explorer (IXPE) is a NASA Small Explorer mission designed to perform polarization measurements in the 2-8 keV band, complemented with imaging, spectroscopy and timing capabilities. At the heart of the focal plane is a set of three polarization-sensitive Gas Pixel Detectors (GPD), each based on a custom ASIC acting as a charge-collecting anode. In this paper we shall review the design, manufacturing, and test of the IXPE focal-plane detectors, with particular emphasis on the connection between the science drivers, the performance metrics and the operational aspects. We shall present a thorough characterization of the GPDs in terms of effective noise, trigger efficiency, dead time, uniformity of response, and spectral and polarimetric performance. In addition, we shall discuss in detail a number of instrumental effects that are relevant for high-level science analysis -- particularly as far as the response to unpolarized radiation and the stability in time are concerned., Comment: To be published in Astroparticle Physics

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

Because many of our X-ray telescopes are optimized towards observing faint sources, observations of bright sources like X-ray binaries in outburst are often affected by instrumental biases. These effects include dead time and photon pile-up, which can dramatically change the statistical inference of physical parameters from these observations. While dead time is difficult to take into account in a statistically consistent manner, simulating dead time-affected data is often straightforward. This structure makes the issue of inferring physical properties from dead time-affected observations fall into a class of problems common across many scientific disciplines. There is a growing number of methods to address them under the name of Simulation-Based Inference (SBI), aided by new developments in density estimation and statistical machine learning. In this paper, we introduce SBI as a principled way to infer variability properties from dead time-affected light curves. We use Sequential Neural Posterior Estimation to estimate the posterior probability for variability properties. We show that this method can recover variability parameters on simulated data even when dead time is variable, and present results of an application of this approach to NuSTAR observations of the galactic black hole X-ray binary GRS 1915+105., Comment: 21 pages, 24 figures; accepted for publication in MNRAS; code available at https://github.com/dhuppenkothen/DeadTimeSBI and data products at http://doi.org/10.5281/zenodo.4670161

We report on the discovery of a new, transient ultraluminous X-ray source (ULX) in the galaxy NGC 7090. This new ULX, which we refer to as NGC 7090 ULX3, was discovered via monitoring with $Swift$ during 2019-20, and to date has exhibited a peak luminosity of $L_{\rm{X}} \sim 6 \times 10^{39}$ erg s$^{-1}$. Archival searches show that, prior to its recent transition into the ULX regime, ULX3 appeared to exhibit a fairly stable luminosity of $L_{\rm{X}} \sim 10^{38}$ erg s$^{-1}$. Such strong long-timescale variability may be reminiscent of the small population of known ULX pulsars, although deep follow-up observations with $XMM$-$Newton$ and $NuSTAR$ do not reveal any robust X-ray pulsation signals. Pulsations similar to those seen from known ULX pulsars cannot be completely excluded, however, as the limit on the pulsed fraction of any signal that remains undetected in these data is $\lesssim$20\%. The broadband spectrum from these observations is well modelled with a simple thin disc model, consistent with sub-Eddington accretion, which may instead imply a moderately large black hole accretor ($M_{\rm{BH}} \sim 40 ~ M_{\odot}$). Similarly, though, more complex models consistent with the super-Eddington spectra seen in other ULXs (and the known ULX pulsars) cannot be excluded given the limited signal-to-noise of the available broadband data. The nature of the accretor powering this new ULX therefore remains uncertain., Comment: 11 pages, 7 figures; accepted for publication in MNRAS

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

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

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

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

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

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

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

The main characteristics in the radio continuum spectra of Supernova Remnants (SNRs) result from simple synchrotron emission. In addition, electron acceleration mechanisms can shape the spectra in specific ways, especially at high radio frequencies. These features are connected to the age and the peculiar conditions of the local interstellar medium interacting with the SNR. Whereas the bulk radio emission is expected at up to $20-50$ GHz, sensitive high-resolution images of SNRs above 10 GHz are lacking and are not easily achievable, especially in the confused regions of the Galactic Plane. In the framework of the early science observations with the Sardinia Radio Telescope in February-March 2016, we obtained high-resolution images of SNRs Tycho, W44 and IC443 that provided accurate integrated flux density measurements at 21.4 GHz: 8.8 $\pm$ 0.9 Jy for Tycho, 25 $\pm$ 3 Jy for W44 and 66 $\pm$ 7 Jy for IC443. We coupled the SRT measurements with radio data available in the literature in order to characterise the integrated and spatially-resolved spectra of these SNRs, and to find significant frequency- and region-dependent spectral slope variations. For the first time, we provide direct evidence of a spectral break in the radio spectral energy distribution of W44 at an exponential cutoff frequency of 15 $\pm$ 2 GHz. This result constrains the maximum energy of the accelerated electrons in the range $6-13$ GeV, in agreement with predictions indirectly derived from AGILE and \textit{Fermi}-LAT gamma-ray observations. With regard to IC443, our results confirm the noticeable presence of a bump in the integrated spectrum around $20-70$ GHz that could result from a spinning dust emission mechanism., Comment: 12 pages, 9 figures

We report on a new NuSTAR observation and on the ongoing Swift XRT monitoring campaign of the peculiar source 1E 161348-5055, located at the centre of the supernova remnant RCW 103, which is recovering from its last outburst in June 2016. The X-ray spectrum at the epoch of the NuSTAR observation can be described by either two absorbed blackbodies ($kT_{BB_1}$ ~ 0.5 keV, $kT_{BB_2}$ ~ 1.2 keV) or an absorbed blackbody plus a power law ($kT_{BB_1}$ ~ 0.6 keV, $\Gamma$ ~ 3.9). The observed flux was ~ 9 $\times$ 10$^{-12}$ erg s$^{-1}$ cm$^{-2}$, ~ 3 times lower than what observed at the outburst onset, but about one order of magnitude higher than the historical quiescent level. A periodic modulation was detected at the known 6.67 hr periodicity. The spectral decomposition and evolution along the outburst decay are consistent with 1E 161348-5055 being a magnetar, the slowest ever detected., Comment: Accepted for publication in MNRAS

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 Astropy project supports and fosters the development of open-source and openly-developed Python packages that provide commonly-needed functionality to the astronomical community. A key element of the Astropy project is the core package Astropy, which serves as the foundation for more specialized projects and packages. In this article, we provide an overview of the organization of the Astropy project and summarize key features in the core package as of the recent major release, version 2.0. We then describe the project infrastructure designed to facilitate and support development for a broader ecosystem of inter-operable packages. We conclude with a future outlook of planned new features and directions for the broader Astropy project., Comment: Minor changes to author list and title. Comments and feedback welcome through the paper source repository: https://github.com/astropy/astropy-v2.0-paper For more information about Astropy, see http://www.astropy.org/

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

In recent years, the cross spectrum has received considerable attention as a means of characterising the variability of astronomical sources as a function of wavelength. While much has been written about the statistics of time and phase lags, the cospectrum has only recently been understood as means of mitigating instrumental effects dependent on temporal frequency in astronomical detectors, as well as a method of characterizing the coherent variability in two wavelength ranges on different time scales. In this paper, we lay out the statistical foundations of the cospectrum, starting with the simplest case of detecting a periodic signal in the presence of white noise. This case is especially relevant for detecting faint X-ray pulsars in detectors heavily affected by instrumental effects, including NuSTAR, Astrosat and IXPE. We show that the statistical distributions of both single and averaged cospectra differ considerably from those for standard periodograms. While a single cospectrum follows a Laplace distribution exactly, averaged cospectra are approximated by a Gaussian distribution only for more than ~30 averaged segments, dependent on the number of trials. We provide an instructive example of a quasi-periodic oscillation in NuSTAR and show that applying standard periodogram statistics leads to underestimated tail probabilities for period detection. We also demonstrate the application of these distributions to a NuSTAR observation of the X-ray pulsar Hercules X-1., Comment: 12 pages, 8 figures, submitted to ApJ; reproducible code available at https://github.com/dhuppenkothen/cospectra-paper ; Comments welcome

In September 2016, the microquasar Cygnus X-3 underwent a giant radio flare, which was monitored for 6 days with the Medicina Radio Astronomical Station and the Sardinia Radio Telescope. Long observations were performed in order to follow the evolution of the flare on a hourly scale, covering six frequency ranges from 1.5 GHz to 25.6 GHz. The radio emission reached a maximum of 13.2 +/- 0.7 Jy at 7.2 GHz and 10 +/- 1 Jy at 18.6 GHz. Rapid flux variations were observed at high radio frequencies at the peak of the flare, together with rapid evolution of the spectral index: alpha steepened from 0.3 to 0.6 within 5 hours. This is the first time that such fast variations are observed, giving support to the evolution from optically thick to optically thin plasmons in expansion moving outward from the core. Based on the Italian network (Noto, Medicina and SRT) and extended to the European antennas (Torun, Yebes, Onsala), VLBI observations were triggered at 22 GHz on five different occasions, four times prior to the giant flare, and once during its decay phase. Flux variations of 2-hour duration were recorded during the first session. They correspond to a mini-flare that occurred close to the core ten days before the onset of the giant flare. From the latest VLBI observation we infer that four days after the flare peak the jet emission was extended over 30 mas., Comment: 13 pages, 9 figures. Accepted for publication in MNRAS

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

Observations of supernova remnants (SNRs) are a powerful tool for investigating the later stages of stellar evolution, the properties of the ambient interstellar medium, and the physics of particle acceleration and shocks. For a fraction of SNRs, multi-wavelength coverage from radio to ultra high-energies has been provided, constraining their contributions to the production of Galactic cosmic rays. Although radio emission is the most common identifier of SNRs and a prime probe for refining models, high-resolution images at frequencies above 5 GHz are surprisingly lacking, even for bright and well-known SNRs such as IC443 and W44. In the frameworks of the Astronomical Validation and Early Science Program with the 64-m single-dish Sardinia Radio Telescope, we provided, for the first time, single-dish deep imaging at 7 GHz of the IC443 and W44 complexes coupled with spatially-resolved spectra in the 1.5-7 GHz frequency range. Our images were obtained through on-the-fly mapping techniques, providing antenna beam oversampling and resulting in accurate continuum flux density measurements. The integrated flux densities associated with IC443 are S_1.5GHz = 134 +/- 4 Jy and S_7GHz = 67 +/- 3 Jy. For W44, we measured total flux densities of S_1.5GHz = 214 +/- 6 Jy and S_7GHz = 94 +/- 4 Jy. Spectral index maps provide evidence of a wide physical parameter scatter among different SNR regions: a flat spectrum is observed from the brightest SNR regions at the shock, while steeper spectral indices (up to 0.7) are observed in fainter cooling regions, disentangling in this way different populations and spectra of radio/gamma-ray-emitting electrons in these SNRs., Comment: 13 pages, 9 figures, accepted for publication to MNRAS on 18 May 2017

[Abridged] The Sardinia Radio Telescope (SRT) is the new 64-m dish operated by INAF (Italy). Its active surface will allow us to observe at frequencies of up to 116 GHz. At the moment, three receivers, one per focal position, have been installed and tested. The SRT was officially opened in October 2013, upon completion of its technical commissioning phase. In this paper, we provide an overview of the main science drivers for the SRT, describe the main outcomes from the scientific commissioning of the telescope, and discuss a set of observations demonstrating the SRT's scientific capabilities. One of the main objectives of scientific commissioning was the identification of deficiencies in the instrumentation and/or in the telescope sub-systems for further optimization. As a result, the overall telescope performance has been significantly improved. As part of the scientific commissioning activities, different observing modes were tested and validated, and first astronomical observations were carried out to demonstrate the science capabilities of the SRT. In addition, we developed astronomer-oriented software tools, to support future observers on-site. The astronomical validation activities were prioritized based on technical readiness and scientific impact. The highest priority was to make the SRT available for joint observations as part of European networks. As a result, the SRT started to participate (in shared-risk mode) in EVN (European VLBI Network) and LEAP (Large European Array for Pulsars) observing sessions in early 2014. The validation of single-dish operations for the suite of SRT first light receivers and backends continued in the following years, and was concluded with the first call for shared-risk/early-science observations issued at the end of 2015., Comment: 26 pages, 18 figures, accepted for publication in Section 13 'Astronomical instrumentation' of Astronomy & Astrophysics

We present $\emph{NuSTAR}$ observations of neutron star (NS) low-mass X-ray binaries: 4U 1636-53, GX 17+2, and 4U 1705-44. We observed 4U 1636-53 in the hard state, with an Eddington fraction, $F_{\mathrm{Edd}}$, of 0.01; GX 17+2 and 4U 1705-44 were in the soft state with fractions of 0.57 and 0.10, respectively. Each spectrum shows evidence for a relativistically broadened Fe K$_{\alpha}$ line. Through accretion disk reflection modeling, we constrain the radius of the inner disk in 4U 1636-53 to be $R_{in}=1.03\pm0.03$ ISCO (innermost stable circular orbit) assuming a dimensionless spin parameter $a_{*}=cJ/GM^{2}=0.0$, and $R_{in}=1.08\pm0.06$ ISCO for $a_{*}=0.3$ (errors quoted at 1 $\sigma$). This value proves to be model independent. For $a_{*}=0.3$ and $M=1.4\ M_{\odot}$, for example, $1.08\pm0.06$ ISCO translates to a physical radius of $R=10.8\pm0.6$ km, and the neutron star would have to be smaller than this radius (other outcomes are possible for allowed spin parameters and masses). For GX 17+2, $R_{in}=1.00-1.04$ ISCO for $a_{*}=0.0$ and $R_{in}=1.03-1.30$ ISCO for $a_{*}=0.3$. For $a_{*}=0.3$ and $M=1.4\ M_{\odot}$, $R_{in}=1.03-1.30$ ISCO translates to $R=10.3-13.0$ km. The inner accretion disk in 4U 1705-44 may be truncated just above the stellar surface, perhaps by a boundary layer or magnetosphere; reflection models give a radius of 1.46-1.64 ISCO for $a_{*}=0.0$ and 1.69-1.93 ISCO for $a_{*}=0.3$. We discuss the implications that our results may have on the equation of state of ultradense, cold matter and our understanding of the innermost accretion flow onto neutron stars with low surface magnetic fields, and systematic errors related to the reflection models and spacetime metric around less idealized neutron stars., Comment: Accepted for publication in ApJ, 12 pages, 5 Figures, 6 Tables