Search

We report the discovery of millihertz quasi-periodic oscillations (mHz QPOs) from the neutron star (NS) low-mass X-ray binary 4U 1730--22 using the Neutron Star Interior Composition Explorer (NICER). After being inactive for almost 50 years, 4U 1730--22 went into outburst twice between June and August 2021, and between February and July 2022. We analyse all the NICER observations of this source, and detect mHz QPOs with a significance > $4\sigma$ in 35 observations. The QPO frequency of the full data set ranged between ~4.5 and ~8.1 mHz with an average fractional rms amplitude of the order of ~2%. The X-ray colour analysis strongly suggests that 4U 1730--22 was in a soft spectral state during the QPO detections. Our findings are consistent with those reported for other sources where the mHz QPOs have been interpreted as the result of a special mode of He burning on the NS surface called marginally stable nuclear burning (MSNB). We conclude that the mHz QPOs reported in this work are also associated with the MSNB, making 4U 1730--22 the eighth source that shows this phenomenology. We discuss our findings in the context of the heat flux from the NS crust., Comment: 9 pages, 6 figures. Accepted for publication in MNRAS

IGR J17407-2808 is an enigmatic and poorly studied X-ray binary that was recently observed quasi-simultaneously with NuSTAR and XMM-Newton. In this paper we report the results of this observational campaign. During the first 60 ks of observation, the source was caught in a relatively low emission state, characterised by a modest variability and an average flux of ~8.3E-13 erg/cm^2/s (4-60 keV). Afterwards, IGR J17407-2808 entered a significantly more active emission state that persisted for the remaining ~40 ks of the NuSTAR observation. During this state, IGR J17407-2808 displayed several fast X-ray flares, featuring durations of ~1-100 s and profiles with either single or multiple peaks. The source flux in the flaring state reached values as high as 2E-9 erg/cm^2/s (4-60 keV), leading to a measured dynamic range during the NuSTAR and XMM-Newton campaign of >~ 10^3. We also analysed available archival photometric near-infrared data of IGR J17407-2808 to improve the constraints available so far on the the nature of the donor star hosted in this system. Our analysis shows that the donor star can be either a rare K or M-type sub-subgiant or an K type main sequence star, or sub-giant star. Our findings support the classification of IGR J17407-2808 as a low-mass X-ray binary. We discuss the source X-ray behaviour as recorded by NuSTAR and XMM-Newton in view of this revised classification., Accepted for publication in Astronomy & Astrophysics

We report comprehensive spectral and temporal properties of the Be/X-ray binary pulsar SMC X-2 using X-ray observations during the 2015 and 2022 outbursts. The pulse profile of the pulsar is unique and strongly luminosity dependent. It evolves from a broad-humped into a double-peaked profile above luminosity 3$\times$10$^{38}$ ergs s$^{-1}$. The pulse fraction of the pulsar is found to be a linear function of luminosity as well as energy. We also studied the spectral evolution of the source during the latest 2022 outburst with NICER. The observed photon index shows a negative and positive correlation below and above the critical luminosity, respectively, suggesting evidence of spectral transition from the sub-critical to super-critical regime. The broadband spectroscopy of four sets of NuSTAR and XRT/NICER data from both outbursts can be described using a cutoff power-law model with a blackbody component. In addition to the 6.4 keV iron fluorescence line, an absorption-like feature is clearly detected in the spectra. The cyclotron line energy observed during the 2015 outburst is below 29.5 keV, however latest estimates in the 2022 outburst suggest a value of 31.5 keV. Moreover, an increase of 3.4 keV is detected in the cyclotron line energy at equal levels of luminosity observed in 2022 with respect to 2015. The observed cyclotron line energy variation is explored in terms of accretion induced screening mechanism or geometrical variation in line forming region., 11 Pages, 11 figures, Accepted for Publication in MNRAS

When accreting X-ray pulsars (XRPs) undergo bright X-ray outbursts, their luminosity-dependent spectral and timing features can be analysed in detail. The XRP GRO J1750-27 recently underwent one of such episodes, during which it was observed with $NuSTAR$ and monitored with $NICER$. Such a data set is rarely available, as it samples the outburst over more than a month at a luminosity that is always exceeding ${\sim}5\times10^{37}\,$erg/s. This value is larger than the typical critical luminosity value, where a radiative shock is formed above the neutron star's surface. Our data analysis of the joint spectra returns a highly ($N_H\sim(5-8)\times10^{22}\,$cm$^{-2}$) absorbed spectrum showing a K$\alpha$ iron line, a soft blackbody component likely originating from the inner edge of the accretion disk, and confirms the discovery of one of the deepest cyclotron lines, at a centroid energy of ${\sim}44\,$keV corresponding to a magnetic field strength of $4.7\times10^{12}\,$G. This value is independently supported by the best-fit physical model for spectral formation in accreting XRPs which, in agreement with recent findings, favours a distance of $14$ kpc and also reflects a bulk-Comptonization dominated accretion flow. Contrary to theoretical expectations and observational evidence from other similar sources, the pulse profiles as observed by $NICER$ through the outburst raise, peak and decay remain remarkably steady. The $NICER$ spectrum, including the iron K$\alpha$ line best-fit parameters, also remain almost unchanged at all probed outburst stages, similar to the pulsed fraction behaviour. We argue that all these phenomena are linked and interpret them as resulting from a saturation effect of the accretion column's emission, which occurs in the high-luminosity regime., Comment: Accepted on A&A -- proof version

We present the pulsar timing analysis of the accreting millisecond X-ray pulsar SWIFT J1749.4-2807 monitored by NICER and XMM-Newton during its latest outburst after almost eleven years of quiescence. From the coherent timing analysis of the pulse profiles, we updated the orbital ephemerides of the system. Large phase jumps of the fundamental frequency phase of the signal are visible during the outburst, consistent with what was observed during the previous outburst. Moreover, we report on the marginally significant evidence for non-zero eccentricity ($e\simeq 4\times 10^{-5}$) obtained independently from the analysis of both the 2021 and 2010 outbursts and we discuss possible compatible scenarios. Long-term orbital evolution of SWIFT J1749.4-2807 suggests a fast expansion of both the NS projected semi-major axis $(x)$, and the orbital period $(P_{\rm orb})$, at a rate of $\dot{x}\simeq 2.6\times 10^{-13}\,\text{lt-s}\,\text{s}^{-1}$ and $\dot{P}_{\rm orb}\simeq 4 \times 10^{-10}\,\text{s}\,\text{s}^{-1}$, respectively. SWIFT J1749.4-2807 is the only accreting millisecond X-ray pulsar, so far, from which the orbital period derivative has been directly measured from appreciable changes on the observed orbital period. Finally, no significant secular deceleration of the spin frequency of the compact object is detected, which allowed us to set a constraint on the magnetic field strength at the polar caps of $B_{PC}, 13 pages, 3 figures, 1 table. Accepted for publication in MNRAS

Swift J1858.6-0814 is a recently discovered X-ray binary notable for extremely strong variability (by factors $>100$ in soft X-rays) in its discovery state. We present the detection of five thermonuclear (Type I) X-ray bursts from Swift J1858.6-0814, implying that the compact object in the system is a neutron star. Some of the bursts show photospheric radius expansion, so their peak flux can be used to estimate the distance to the system. The peak luminosity, and hence distance, can depend on several system parameters; for the most likely values, a high inclination and a helium atmosphere, $D=12.8_{-0.6}^{+0.8}$ kpc, although systematic effects allow a conservative range of $9-18$ kpc. Before one burst, we detect a QPO at $9.6\pm0.5$ mHz with a fractional rms amplitude of $2.2\pm0.2$% ($0.5-10$ keV), likely due to marginally stable burning of helium; similar oscillations may be present before the other bursts but the light curves are not long enough to allow their detection. We also search for burst oscillations but do not detect any, with an upper limit in the best case of 15% fractional amplitude (over $1-8$ keV). Finally, we discuss the implications of the neutron star accretor and this distance on other inferences which have been made about the system. In particular, we find that Swift J1858.6-0814 was observed at super-Eddington luminosities at least during bright flares during the variable stage of its outburst., 14 pages, 8 figures, 2 tables, MNRAS accepted

We report the discovery of PSR~J2251$-$3711, a radio pulsar with a spin period of 12.1 seconds, the second longest currently known. Its timing parameters imply a characteristic age of 15 Myr, a surface magnetic field of $1.3 \times 10^{13}$~G and a spin-down luminosity of $2.9 \times 10^{29}~\mathrm{erg~s}^{-1}$. Its dispersion measure of 12.12(1)~$\mathrm{pc}~\mathrm{cm}^{-3}$ leads to distance estimates of 0.5 and 1.3 kpc according to the NE2001 and YMW16 Galactic free electron density models, respectively. Some of its single pulses show an uninterrupted 180 degree sweep of the phase-resolved polarization position angle, with an S-shape reminiscent of the rotating vector model prediction. However, the fact that this sweep occurs at different phases from one pulse to another is remarkable and without straightforward explanation. Although PSR~J2251$-$3711 lies in the region of the $P-\dot{P}$ parameter space occupied by the X-ray Isolated Neutron Stars (XINS), there is no evidence for an X-ray counterpart in our Swift XRT observation; this places a 99\%-confidence upper bound on its unabsorbed bolometric thermal luminosity of $1.1 \times 10^{31}~(d / 1~\mathrm{kpc})^2~\mathrm{erg/s}$ for an assumed temperature of 85 eV, where $d$ is the distance to the pulsar. Further observations are needed to determine whether it is a rotation-powered pulsar with a true age of at least several Myr, or a much younger object such as an XINS or a recently cooled magnetar. Extreme specimens like PSR J2251$-$3711 help bridge populations in the so-called neutron star zoo in an attempt to understand their origins and evolution., Accepted for publication in MNRAS, 13 pages, 8 figures

Because of its complex clumpy wind, prominent cyclotron resonant scattering features, intrinsic variability, and convenient physical parameters (close distance, high inclination, and small orbital separation), which facilitate the observation and analysis of the system, Vela X-1 is one of the key systems for understanding 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 timescale using a partially covered power law continuum with a Fermi-Dirac cutoff to model the continuum and local absorption. We are able to 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 were not able to 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., This work has been partially funded by the Bundesministerium für Wirtschaft und Energie under the Deutsches Zentrum für Luft- und Raumfahrt Grants 50 OR 1915. S.M.N. acknowledges funding under project RTI2018-096686-B-C21 funded by MCIN/AEI/10.13039/501100011033 and by ‘ERDF A way of making Europe’, and by the Unidad de Excelencia María de Maeztu, ref. MDM-2017-0765. E.S.L. acknowledges support by DFG grant 1830Wi1860/11-1 and RFBR grant 18-502-12025. K.P. acknowledges support by NASA under award number 80GSFC17M0002.

Context. The Magellanic Clouds are our nearest star-forming galaxies. While the population of high-mass X-ray binaries (HMXBs) in the Small Magellanic Cloud (SMC) is relatively well studied, our knowledge about the Large Magellanic Cloud (LMC) is far from complete given its large angular extent and insufficient coverage with X-ray observations. Aims. We conducted a search for new HMXBs in the LMC using data from eROSITA, the soft X-ray instrument on board the Spektrum-Roentgen-Gamma (SRG) satellite. Methods. After confirming the nature of eRASSU J052914.9-662446 as a hard X-ray source positionally coincident with an early type star, we followed it up with optical spectroscopic observations from South African Large Telescope (SALT) and a dedicated NuSTAR observation. Results. We study the broadband timing and spectral behaviour of the newly discovered HMXB eRASSU J052914.9-662446 through eROSITA, Swift and NuSTAR data in X-rays and the Optical Gravitational Lensing Experiment (OGLE) and SALT RSS data at optical wavelength. We report on detection of the spin period at 1412 s and suggest an orbital period of the system of ~151 days, and thereby establish eRASSU J052914.9-662446 as an accreting pulsar. Further, through optical spectroscopic observations and the existence of H alpha emission the source is identified as a Be X-ray binary pulsar in the LMC. We also investigate the variability of the source in the optical and X-ray regime over the last decades and provide estimates on the possible magnetic field strength of the neutron star., Comment: 11 pages, 13 figures, 3 tables, accepted for publication in A&A

We present the detection of 51 thermonuclear X-ray bursts observed from 4U 1636-536 by the Neutron Star Interior Composition Explorer (NICER) over the course of a three year monitoring campaign. We performed time resolved spectroscopy for 40 of these bursts and showed the existence of a strong soft excess in all the burst spectra. The excess emission can be characterized by the use of a scaling factor (f_a method) to the persistent emission of the source, which is attributed to the increased mass accretion rate on to the neutron star due to Poynting-Robertson drag. The soft excess emission can also be characterized by the use of a model taking into account the reflection of the burst emission off of the accretion disk. We also present time resolved spectral analysis of 5 X-ray bursts simultaneously observed by NICER and AstroSat, which confirm the main results with even greater precision. Finally, we present evidence for Compton cooling using 7 X-ray bursts observed contemporaneously with \nustar, by means of a correlated decrease in the hard X-ray lightcurve of 4U 1636-536 as the bursts start., Comment: Accepted for publication in the Astrophysical Journal

We present the discovery of eclipses in the X-ray light curves of the X-ray binary Swift J1858.6–0814. From these, we find an orbital period of P=76841.3+1.3−1.4 s (≈21.3 h) and an eclipse duration of tec=4098+17−18 s (≈1.14 h). We also find several absorption dips during the pre-eclipse phase. From the eclipse duration to orbital period ratio, the inclination of the binary orbit is constrained to i > 70°. The most likely range for the companion mass suggests that the inclination is likely to be closer to this value than 90. The eclipses are also consistent with earlier data, in which strong variability (‘flares’) and the long orbital period prevent clear detection of the period or eclipses. We also find that the bright flares occurred preferentially in the post-eclipse phase of the orbit, likely due to increased thickness at the disc-accretion stream interface preventing flares being visible during the pre-eclipse phase. This supports the notion that variable obscuration is responsible for the unusually strong variability in Swift J1858.6–0814., Facultad de Ciencias Astronómicas y Geofísicas, Instituto Argentino de Radioastronomía

MAXI J1807+132 is a low-mass X-ray binary (LMXB) first detected in outburst in 2017. Observations during the 2017 outburst did not allow for an unambiguous identification of the nature of the compact object. MAXI J1807+132 was detected in outburst again in 2019 and was monitored regularly with NICER. In this paper we report on five days of observations during which we detected three thermonuclear (Type-I) X-ray bursts, identifying the system as a neutron star LMXB. Time-resolved spectroscopy of the three Type-I bursts revealed typical characteristics expected for these phenomena. All three Type-I bursts show slow rises and long decays, indicative of mixed H/He fuel. We find no strong evidence that any of the Type-I bursts reached the Eddington Luminosity; however, under the assumption that the brightest X-ray burst underwent photospheric radius expansion, we estimate a, 9 pages, 4 figures

We report on NICER observations of the Magnetar SGR~1935+2154, covering its 2020 burst storm and long-term persistent emission evolution up to $\sim90$ days post outburst. During the first 1120~seconds taken on April 28 00:40:58 UTC we detect over 217 bursts, corresponding to a burst rate of $>0.2$ bursts s$^{-1}$. Three hours later the rate is at 0.008 bursts s$^{-1}$, remaining at a comparatively low level thereafter. The $T_{90}$ burst duration distribution peaks at 840~ms; the distribution of waiting times to the next burst is fit with a log-normal with an average of 2.1 s. The 1-10 keV burst spectra are well fit by a blackbody, with an average temperature and area of $kT=1.7$ keV and $R^2=53$ km$^2$. The differential burst fluence distribution over $\sim3$ orders of magnitude is well modeled with a power-law form $dN/dF\propto F^{-1.5\pm0.1}$. The source persistent emission pulse profile is double-peaked hours after the burst storm. We find that the bursts peak arrival times follow a uniform distribution in pulse phase, though the fast radio burst associated with the source aligns in phase with the brighter peak. We measure the source spin-down from heavy-cadence observations covering days 21 to 39 post-outburst, $\dot��=-3.72(3)\times10^{-12}$ Hz s$^{-1}$; a factor 2.7 larger than the value measured after the 2014 outburst. Finally, the persistent emission flux and blackbody temperature decrease rapidly in the early stages of the outburst, reaching quiescence 40 days later, while the size of the emitting area remains unchanged., Accepted for publication in the Astrophysical Journal Letters

In November 2019, MAXI detected an X-ray outburst from the known Be X-ray binary system RX J0209.6-7427 located in the outer wing of the Small Magellanic Cloud. We followed the outburst of the system with NICER which led to the discovery of X-ray pulsations with a period of 9.3 s. We analyzed simultaneous X-ray data obtained with NuSTAR and NICER allowing us to characterize the spectrum and provide an accurate estimate of its bolometric luminosity. During the outburst the maximum broadband X-ray luminosity of the system reached $1-2\times10^{39}$ erg/s, thus exceeding by about one order of magnitude the Eddington limit for a typical 1.4 $M_{\odot}$ mass neutron star (NS). Monitoring observations with Fermi/GBM and NICER allowed us to study the spin evolution of the NS and compare it with standard accretion torque models. We found that the NS magnetic field should be of the order of $3\times10^{12}$ G. We conclude that RX J0209.6-7427 exhibited one of the brightest outbursts observed from a Be X-ray binary pulsar in the Magellanic Clouds, reaching similar luminosity level to the 2016 outburst of SMC X-3. Despite the super-Eddington luminosity of RX J0209.6-7427, the NS appears to have only a moderate magnetic field strength., 10 pages, 8 figures, accepted for publication in MNRAS

We report on the phase-coherent timing analysis of the accreting millisecond X-ray pulsar IGR J17591-2342, using Neutron Star Interior Composition Explorer (NICER) data taken during the outburst of the source between 2018 August 15 and 2018 October 17. We obtain an updated orbital solution of the binary system. We investigate the evolution of the neutron star spin frequency during the outburst, reporting a refined estimate of the spin frequency and the first estimate of the spin frequency derivative ($\dot{\nu} \sim -7\times 10^{-14}$ Hz s$^{-1}$), confirmed independently from the modelling of the fundamental frequency and its first harmonic. We further investigate the evolution of the X-ray pulse phases adopting a physical model that accounts for the accretion material torque as well as the magnetic threading of the accretion disc in regions where the Keplerian velocity is slower than the magnetosphere velocity. From this analysis we estimate the neutron star magnetic field $B_{eq} = 2.8(3)\times10^{8}$ G. Finally, we investigate the pulse profile dependence on energy finding that the observed behaviour of the pulse fractional amplitude and lags as a function of energy are compatible with a thermal Comptonisation of the soft photons emitted from the neutron star caps., Comment: Submitted to MNRAS. 10 pages, 6 figures, 2 tables

The accretion-powered X-ray pulsar GX 301-2 was observed with the balloon-borne X-Calibur hard X-ray polarimeter during late December 2018, with contiguous observations by the NICER X-ray telescope, the Swift X-ray Telescope and Burst Alert Telescope, and the Fermi Gamma-ray Burst Monitor spanning several months. The observations detected the pulsar in a rare apastron flaring state coinciding with a significant spin-up of the pulsar discovered with the Fermi GBM. The X-Calibur, NICER, and Swift observations reveal a pulse profile strongly dominated by one main peak, and the NICER and Swift data show strong variation of the profile from pulse to pulse. The X-Calibur observations constrain for the first time the linear polarization of the 15-35 keV emission from a highly magnetized accreting neutron star, indicating a polarization degree of (27+38-27)% (90% confidence limit) averaged over all pulse phases. We discuss the spin-up and the X-ray spectral and polarimetric results in the context of theoretical predictions. We conclude with a discussion of the scientific potential of future observations of highly magnetized neutron stars with the more sensitive follow-up mission XL-Calibur., Comment: Accepted for publication in the Astrophysical Journal. 20 pages, 19 figures, 4 tables

PSR J1813-1749 is one of the most energetic rotation-powered pulsars known, producing a pulsar wind nebula (PWN) and gamma-ray and TeV emission, but whose spin period is only measurable in X-ray. We present analysis of two Chandra datasets that are separated by more than ten years and recent NICER data. The long baseline of the Chandra data allows us to derive a pulsar proper motion mu_R.A.=-(0.067"+/-0.010") yr^-1 and mu_decl.=-(0.014"+/-0.007") yr^-1 and velocity v_perp~900-1600 km/s (assuming a distance d=3-5 kpc), although we cannot exclude a contribution to the change in measured pulsar position due to a change in brightness structure of the PWN very near the pulsar. We model the PWN and pulsar spectra using an absorbed power law and obtain best-fit absorption NH=(13.1+/-0.9)x10^22 cm^-2, photon index Gamma=1.5+/-0.1, and 0.3-10 keV luminosity Lx~5.4x10^34 erg/s (d/5 kpc)^2 for the PWN and Gamma=1.2+/-0.1 and Lx~9.3x10^33 erg/s (d/5 kpc)^2 for PSR J1813-1749. These values do not change between the 2006 and 2016 observations. We use NICER observations from 2019 to obtain a timing model of PSR J1813-1749, with spin frequency nu=22.35 Hz and spin frequency time derivative nudot=(-6.428+/-0.003)x10^-11 Hz/s. We also fit nu measurements from 2009-2012 and our 2019 value and find a long-term spin-down rate nudot=(-6.3445+/-0.0004)x10^-11 Hz/s. We speculate that the difference in spin-down rates is due to glitch activity or emission mode switching., 9 pages, 5 figures; accepted for publication in MNRAS

We review more than 10 years of continuous monitoring of accreting X-ray pulsars with the all-sky Gamma-ray Burst Monitor (GBM) aboard the Fermi Gamma-ray Space Telescope. Our work includes data from the start of GBM operations in August 2008, through to November 2019. Pulsations from 39 accreting pulsars are observed over an energy range of $10-50\,$keV by GBM. The GBM Accreting Pulsars Program (GAPP) performs data reduction and analysis for each accreting pulsar and makes histories of the pulse frequency and pulsed flux publicly available. We examine in detail the spin histories, outbursts and torque behaviors of the persistent and transient X-ray pulsars observed by GBM. The spin period evolution of each source is analyzed in the context of disk-accretion and quasi-spherical settling accretion driven torque models. Long-term pulse frequency histories are also analyzed over the GBM mission lifetime and compared to those available from the previous BATSE all-sky monitoring mission, revealing previously unnoticed episodes in some of the analyzed sources (such as a torque reversal in 2S 1845-024). We obtain new, or update known, orbital solutions for three sources. Our results demonstrate the capabilities of GBM as an excellent instrument for monitoring accreting X-ray pulsars and its important scientific contribution to this field., Comment: Published online by ApJ. Proofs edition + more typos fixed

The capability of the Fermi Gamma-ray Burst Monitor (GBM) to localize gamma-ray bursts (GRBs) is evaluated for two different automated algorithms: the GBM Team's RoboBA algorithm and the independently developed BALROG algorithm. Through a systematic study utilizing over 500 GRBs with known locations from instruments like Swift and the Fermi LAT, we directly compare the effectiveness of, and accurately estimate the systematic uncertainty for, both algorithms. We show simple adjustments to the GBM Team's RoboBA, in operation since early 2016, yields significant improvement in the systematic uncertainty, removing the long tail identified in the systematic, and improves the overall accuracy. The systematic uncertainty for the updated RoboBA localizations is $1.8^\circ$ for 52% of GRBs and $4.1^\circ$ for the remaining 48%. Both from public reporting by BALROG and our systematic study, we find the systematic uncertainty of $1-2^\circ$ quoted in GCN circulars for bright GRBs localized by BALROG is an underestimate of the true magnitude of the systematic, which we find to be $2.7^\circ$ for 74% of GRBs and $33^\circ$ for the remaining 26%. We show that, once the systematic uncertainty is considered, the RoboBA 90% localization confidence regions can be more than an order of magnitude smaller in area than those produced by BALROG., Comment: Accepted for publication in ApJ. Added Appendix showing localization comparisons for each GRB in the near-realtime public reporting sample

The Neutron Star Interior Composition Explorer (NICER) is an X-ray astrophysics payload on the International Space Station. It enables unprecedented high-precision timing of millisecond pulsars without the pulse broadening and delays due to dispersion and scattering within the interstellar medium that plague radio timing. We present initial timing results from a year of data on the millisecond pulsars PSR B1937+21 and PSR J0218+4232, and nine months of data on PSR B1821-24. NICER time-of-arrival uncertainties for the three pulsars are consistent with theoretical lower bounds and simulations based on their pulse shape templates and average source and background photon count rates. To estimate timing stability, we use the $\sigma_z$ measure, which is based on the average of the cubic coefficients of polynomial fits to subsets of timing residuals. So far we are achieving timing stabilities $\sigma_z \approx 3 \times 10^{-14}$ for PSR B1937+21 and on the order of $10^{-12}$ for PSRs B1821$-$24 and J0218+4232. Within the span of our \textit{NICER} data we do not yet see the characteristic break point in the slope of $\sigma_z$; detection of such a break would indicate that further improvement in the cumulative root-mean-square (RMS) timing residual is limited by timing noise. We see this break point in our comparison radio data sets for PSR B1821-24 and PSR B1937+21 on time scales of $> 2$ years., Comment: 30 pages, 15 figures; accepted to ApJ

Accreting X-ray pulsars (XRPs) undergo luminous X-ray outbursts during which the spectral and timing behavior of the neutron star can be studied in detail. We analyze a $NuSTAR$ observation of the XRP XTE J1858+034 during its outburst in 2019. The spectrum is fit with a phenomenological, a semi-empirical and a physical spectral model. A candidate cyclotron line is found at $48\,$keV, implying a magnetic field of $5.4\times10^{\rm 12}\,$G at the site of emission. This is also supported by the physical best-fit model. We propose an orbital period of about $81$ days based on the visual inspection of the X-ray outbursts recurrence time. Based on $Fermi$ Gamma-ray Burst Monitor data, the standard disk accretion-torque theory allowed us to infer a distance of $10.9\pm1.0\,$kpc. Pulse profiles are single-peaked and show a pulsed fraction that is strongly energy-dependent at least up to $40$ keV., Comment: Accepted for publication in ApJ; 9 pages, 2 tables, 7 figures

AX J0049.4-7323 is a Be/X-ray binary that shows an unusual and poorly understood optical variability that consists of periodic and bright optical outbursts, simultaneous with X-ray outbursts, characterised by a highly asymmetric profile. The periodicity of the outbursts is thought to correspond to the orbital period of the neutron star. To understand the behaviour shown by this source, we performed the first multi-wavelength monitoring campaign during the periastron passage of December 2017. The monitoring lasted for ~37 days and consisted of X-ray, near-ultraviolet, and optical data from the Neil Gehrels Swift Observatory, the optical I band from the OGLE survey, and spectroscopic observations of the H-alpha line performed with the 3.9 m Anglo-Australian Telescope. These observations revealed AX J0049.4-7323 during an anomalous outburst having remarkably different properties compared to the previous ones. In the I band, it showed a longer rise timescale (~60 days instead of 1-5 days) and a longer decay timescale. At the peak of the outburst, it showed a sudden increase in luminosity in the I band, corresponding to the onset of the X-ray outburst. The monitoring of the H-alpha emission line showed a fast and highly variable profile composed of three peaks with variable reciprocal brightness. We interpreted these results as a circumstellar disc warped by tidal interactions with the neutron star in a high eccentricity orbit during its periastron passage. The fast jump in optical luminosity at the peak of the outburst and the previous asymmetric outbursts might be caused by the reprocessing of the X-ray photons in the circumstellar disc or the tidal displacement of a large amount of material from the circumstellar disc or the outer layers of the donor star during the periastron passage of the neutron star, which led to an increase in size of the region emitting in the I band., Accepted for publication in Astronomy & Astrophysics

Be/X-ray binaries are the most populous class of High Mass X-ray Binaries. Their X-ray duty cycle is tightly related to the optical companion wind activity, which in turn can be studied through optical spectroscopical dedicated observations. We study optical spectral features of the Be circumstellar disk to test their long-term variability and their relation with the X-ray activity. Special attention has been given to the H$\alpha$ emission line, that is one of the best tracers of the disk conditions. We obtained optical broadband medium resolution spectra from a dedicated campaign with the Anglo-Australian Telescope and the Southern African Large Telescope in 2014-2015. Data span over one entire binary orbit, and cover both X-ray quiescent and moderately active periods. We used Balmer emission lines to follow the evolution of the circumstellar disk. We observe prominent spectral features, like double-peaked H$\alpha$ and H$\beta$ emission lines. The H$\alpha$ V/R ratio significantly changes over a time scale of about one year. Our observations are consistent with a system observed at a large inclination angle ($i\gtrsim60^{\circ}$). The derived circumstellar disk size shows that the disk evolves from a configuration that prevents accretion onto the neutron star, to one that allows only moderate accretion. This is in agreement with the contemporary observed X-ray activity. Our results are interpreted within the context of inefficient tidal truncation of the circumstellar disk, as expected for this source's binary configuration. We derived the H$\beta$-emitting region size, which results about half of the corresponding H$\alpha$-emitting disk, and constrain the luminosity class of V850 Cen as III-V, consistent with the previously proposed class.}, Comment: 7 pages, 5 figures, Accepted for publication on A&A

We present evidence for three spin-down glitches (or `anti-glitches') in the ultraluminous accreting X-ray pulsar NGC 300 ULX-1 in timing observations made with the Neutron Star Interior Composition Explorer (NICER). Our timing analysis reveals three sudden spin-down events of magnitudes $\Delta\nu = -23, -30,$ and $-43 \,\mu$Hz (fractional amplitudes $\Delta\nu/\nu = -4.4, -5.5,$ and $-7.7 \times 10^{-4}$). We determined fully phase-coherent timing solutions through the first two glitches, giving us high confidence in their detection, while the third candidate glitch is somewhat less secure. These are larger in magnitude (and opposite in sign) than any known radio pulsar glitch. This may be caused by the prolonged rapid spin-up of the pulsar causing a sudden transfer of angular momentum between the superfluid and non-superfluid components of the star. We find no evidence for profile or spectral changes at the epochs of the glitches, supporting the conclusion that these are due to the same process as in normal pulsar glitches, but in reverse., Comment: ApJ, accepted

In July 2015, the high-mass X-ray binary V 0332+53 underwent a giant outburst, a decade after the previous one. This system hosts a strongly magnetized neutron star. During the 2004--2005 outburst, an anti-correlation between the centroid energy of its fundamental cyclotron resonance scattering features (CRSFs) and the X-ray luminosity was observed. The long ($\approx 100$ d) and bright ($L_{\rm x} \approx 10^{38}$ erg s$^{-1}$) 2015 outburst offered the opportunity to study during another outburst the unique properties of the fundamental CRSF and its dependence on the X-ray luminosity. The source was observed by the INTEGRAL satellite for $\sim 330$ ks. We exploit the spectral resolution at high energies of the SPectrometer on INTEGRAL (SPI) and the Joint European X-ray Monitors to characterize its spectral properties, focusing in particular on the CRSF-luminosity dependence. We complement the data of the 2015 outburst with those collected by SPI in 2004--2005 and left unpublished so far. We find a highly significant anti-correlation of the centroid energy of the fundamental CRSF and the $3-100$ keV luminosity of $E_1 \propto -0.095(8)L_{37}$ keV. This trend is observed for both outbursts. We confirm the correlation between the width of the fundamental CRSF and the X-ray luminosity previously found in the JEM-X and IBIS dataset of the 2004--2005 outburst. By exploiting the RXTE/ASM and Swift/BAT monitoring data we also report on the detection of a $\sim 34$ d modulation superimposed on the mean profiles and roughly consistent with the orbital period of the pulsar. We discuss possible interpretations of such variability., Accepted for publication on A&A

We present detection and analysis of faint X-ray sources in the Chandra deep field south (CDFS) using the 4 Msec Chandra observation and adopting a new detection algorithm, based on a targeted search at the position of known high-z galaxies. This optimized technique results in the identification of 54 z>3 AGNs, 29 of which are new detections. Applying stringent completeness criteria, we derive AGN luminosity functions in the redshift bins 3-4, 4-5 and >5.8 and for 42.753 (18+17-10%). Their optical counterparts are not strongly reddened and we thus conclude that the size of the X-ray absorber is likely smaller than the dust sublimation radius. We finally report the discovery of a highly star-forming galaxy at z=3.47. If confirmed, this would be one of the farthest objects in which stellar sources are detected in X-rays., Comment: A&A in press

We report the detection of the 6.4 keV Iron K{\alpha} emission line in two infrared-luminous, massive, star-forming BzK galaxies at z = 2.578 and z = 2.90 in the CDF-S. The Chandra 4 Ms spectra of BzK4892 and BzK8608 show a reflection dominated continuum with strong Iron lines, with rest-frame equivalent widths EW-2.3 keV and 1.2 keV, respectively, demonstrating Compton thick obscuration of the central AGN. For BzK8608 the line identification closely matches the existing photometric redshift derived from the stellar emission. We use the observed luminosities of the Iron K{\alpha} line, of the rest-frame mid-IR continuum and of the UV rest-frame narrow emission lines to infer intrinsic L(2-10 keV)>1e44 erg/s, about 1.0-2.5 dex larger than the observed ones, hence confirming the presence of an absorber with N(H) > 1e24 cm-2. The two BzK galaxies have stellar masses of 5e10 M\odot and, based on VLA 1.4 GHz and submm 870{\mu}m observations, they appear to host vigorous starburst activity with SFR - 300-700 M\odot/yr that is also optically thick. We estimate that the AGN might also conceivably account for an important fraction of the bolometric far-IR emission of the galaxies. The implied volume density of Compton thick (CT) AGN with L(2-10 keV)>1e44 erg/s is in agreement with predictions from X-ray background synthesis models. These sources provide one of the first clearcut observations of the long-sought phase of simultaneous, heavily obscured quasar and star formation activity, predicted by models of massive galaxy evolution at high redshifts., Comment: ApJ Letters, in press. 6 pages, 3 figures