Search

Compact binary millisecond pulsars (MSPs), with orbital periods $\lesssim1$d, are key to understanding binary evolution involving massive neutron stars (NSs). Due to ablation of the companion by the rapidly spinning pulsar, these systems are also known as spiders, categorized into two main branches: redbacks (RBs; companion mass in the range of 0.1 to 0.5 M$_{\odot}$) and black widows (BWs; companion mass $\lesssim$ 0.1 M$_{\odot}$). We present models of low- and intermediate-mass X-ray binaries (XRBs) and compare them with observations of Galactic spiders (including the presence/absence of hydrogen lines in their optical spectra). We constrain and quantify the interaction between pulsar and companion. Using MESA, we create the allowed initial parameter space, modeling the irradiation of the companion by the pulsar wind. For the first time in MESA, we include the detailed evolution of the pulsar spin. Efficient mass accretion onto the NS ($70\%$ of mass transferred) with an X-ray irradiated disc, followed by strong irradiation of the companion can explain most properties of the observed spiders. Our RBs continue to the BW regime, connecting the two branches of spiders. Our models explain the lack of hydrogen in some observed BWs with ultra-light companions. The mass required to spin-up a NS to sub-milliseconds is high enough to collapse it into a black hole. Finally, we discuss the formation of RB-like spiders with giant companions and orbital periods of several days (huntsmen) showing that they are unlikely to produce super-massive NSs (maximum accreted mass $\lesssim$0.5 M$_{\odot}$). Cannibalistic MSP binary formation depends heavily on the interplay between accretion onto the pulsar and pulsar wind irradiation. Our work supports earlier claims that RBs evolve into BWs. We also show that the fastest-spinning pulsars may collapse before reaching sub-millisecond spin periods., Comment: 17 pages, 13 figures, submitted to A&A

We present GTC-OSIRIS phase-resolved optical spectroscopy of three compact binary MSPs, or `spiders': PSR J1048+2339, PSR J1810+1744, and (for the first time) PSR J1908+2105. For the companion in each system, the temperature is traced throughout its orbit, and radial velocities are measured. The radial velocities are found to vary with the absorption features used when measuring them, resulting in a lower radial velocity curve semi-amplitude measured from the day side of two of the systems when compared to the night: for J1048 ($K_\mathrm{day} = 344 \pm 4$ km s$^{-1}$, $K_\mathrm{night} = 372 \pm 3$ km s$^{-1}$) and, tentatively, for J1810 ($K_\mathrm{day} = 448 \pm 19$ km s$^{-1}$, $K_\mathrm{night} = 491 \pm 32$ km s$^{-1}$). With existing inclination constraints, this gives the neutron star (NS) and companion masses $M_\mathrm{NS} = 1.50 - 2.04$ $M_\odot$ and $M_2 = 0.32 - 0.40$ $M_\odot$ for J1048, and $M_\mathrm{NS} > 1.7$ $M_\odot$ and $M_2 = 0.05 - 0.08$ $M_\odot$ for J1810. For J1908, we find an upper limit of $K_2 < 32$ km s$^{-1}$, which constrains its mass ratio $q = M_2 / M_\mathrm{NS} > 0.55$ and inclination $i < 6.0^\circ$, revealing the previously misunderstood system to be the highest mass ratio, lowest inclination redback yet. This raises questions for the origins of its substantial radio eclipses. Additionally, we find evidence of asymmetric heating in J1048 and J1810, and signs of metal enrichment in J1908. We also explore the impact of inclination on spectroscopic temperatures, and demonstrate that the temperature measured at quadrature ($\phi = 0.25, 0.75$) is essentially independent of inclination, and thus can provide additional constraints on photometric modelling., Comment: Submitted to MNRAS. 18 pages, 19 figures

The companion to PSR J1622-0315, one of the most compact known redback millisecond pulsars, shows extremely low irradiation despite its short orbital period. We model this system to determine the binary parameters, combining optical observations from NTT in 2017 and NOT in 2022 with the binary modeling code ICARUS. We find a best-fit neutron star mass of $2.3 \pm 0.4\,\text{M}_\odot $, and a companion mass of $0.15 \pm 0.02\,\text{M}_\odot$. We detect for the first time low-level irradiation from asymmetry in the minima as well as a change in the asymmetry of the maxima of its light curves over five years. Using star spot models, we find better fits than those from symmetric direct heating models, with consistent orbital parameters. We discuss an alternative scenario where the changing asymmetry is produced by a variable intrabinary shock. In summary, we find that PSR J1622-0315 combines low irradiation with variable light curve asymmetry, and a relatively high neutron star mass., Comment: 18 pages, 11 figures

Radio emission from pulsars can be used to map out their distances through dispersion measure (DM), which quantifies the amount of radio pulse dispersion. However, this method relies on accurately modelling the free electron density in the line of sight. Here, we present a detailed study of the multiwavelength emission from PSR J1720$-$0534, a black widow compact binary millisecond pulsar discovered in 2021, which the latest electron density model of the Galaxy (Yao et al. 2017) places at only 191 pc. We obtained and analysed deep multiwavelength observations in the $\gamma$-ray (Fermi-Large Area Telescope, 2008-2022), optical (Las Cumbres Observatory, 2.7 h), near-infrared (Nordic Optical Telescope, 3.5 h), and X-ray (Swift-X-Ray Telescope, 10 ks) bands. We found no significant detection of $\gamma$-ray, optical, near-infrared, or X-ray counterparts around the radio-timing position of PSR J1720$-$0534, which we thus nickname 'the invisible black widow'. Employing the most constraining near-infrared limit ($J>23.4$ mag), we established a lower limit on the source distance, $d>1.1$ kpc, assuming conservative properties for the black widow companion star. This distance lower limit differs drastically (by a factor of more than 5) from the Yao et al. DM distance estimate. We attribute this difference to the inclusion in the Yao et al. model of a large and dense component towards the North Polar Spur. Considering our results and recent parallax distances to other pulsars in this direction, we argue that such a local and large component in the electron density model of the Galaxy is unnecessary., Comment: 12 pages, 11 figures, 3 tables. This article has been accepted for publication in MNRAS published by Oxford University Press on behalf of the Royal Astronomical Society. Typos corrected, acknowledgements added

In this paper, we study compact binary millisecond pulsars with low- and very low-mass companion stars (spiders) in the Galactic field, using data from the latest Gaia data release (DR3). We infer the parallax distances of the optical counterparts to spiders, which we use to estimate optical and X-ray luminosities. We compare the parallax distances to those derived from radio pulse dispersion measures and find that they have systematically larger values, by 40% on average. We also test the correlation between X-ray and spin-down luminosities, finding that most redbacks have a spin-down to X-ray luminosity conversion efficiency of $\sim$0.1%, indicating a contribution from the intrabinary shock. On the other hand, most black widows have an efficiency of $\sim$0.01%, similar to the majority of the pulsar population. Finally, we find that the bolometric optical luminosity significantly correlates with the orbital period, with a large scatter due to different irradiated stellar temperatures and binary properties. We interpret this correlation as the effect of the increasing size of the Roche Lobe radius with the orbital period. With this newly found correlation, an estimate of the optical magnitude can be obtained from the orbital period and a distance estimate., Comment: 24 pages, 17 figures, 9 tables. This article has been accepted for publication in MNRAS published by Oxford University Press on behalf of the Royal Astronomical Society

We present the first multi-band optical light curves of PSR J1622-0315, among the most compact known redback binary millisecond pulsars, with an orbital period Porb=3.9 h. We find a flux modulation with two maxima per orbital cycle and a peak-to-peak amplitude of about 0.3 mag, which we attribute to the ellipsoidal shape of the tidally distorted companion star. The optical colours imply a late-F to early-G spectral type companion and do not show any detectable temperature changes along the orbit. This suggests that the irradiation of the star's inner face by the pulsar wind is unexpectedly missing despite its short orbital period. To interpret these results, we introduce a new parameter fsd, defined as the ratio between the pulsar wind flux intercepted by the companion star and the companion intrinsic flux. This flux ratio fsd, which depends on the spin-down luminosity of the pulsar, the base temperature of the companion and the orbital period, can be used to quantify the effect of the pulsar wind on the companion star and turns out to be the most important factor in determining whether the companion is irradiated or not. We find that the transition between these two regimes occurs at fsd=2-4 and that the value for PSR J1622-0315 is fsd=0.7, placing it firmly in the non-irradiated regime.

Many astrophysical and terrestrial scenarios involving magnetic fields can be approached in axial geometry. Although the smoothed particle hydrodynamics (SPH) technique has been successfully extended to magneto-hydrodynamics (MHD), a well-verified, axisymmetric MHD scheme based on such technique does not exist yet. In this work we fill that gap in the scientific literature and propose and check a novel axisymmetric MHD hydrodynamic code, that can be applied to physical problems which display the adequate geometry. We show that the hydrodynamic code built following these axisymmetric hypothesis is able to produce similar results than standard 3D-SPMHD codes with equivalent resolution but with much lesser computational load., Comment: 19 pages, 13 figures, accepted for publication in MNRAS. arXiv admin note: text overlap with arXiv:2206.05324

We present the results of a deep study of the neutron star (NS) population in the globular cluster M28 (NGC 6626), using the full 330-ks 2002-2015 ACIS dataset from the Chandra X-ray Observatory and coordinated radio observations taken with the Green Bank Telescope (GBT) in 2015. We investigate the X-ray luminosity (Lx), spectrum, and orbital modulation of the 7 known compact binary millisecond pulsars (MSPs) in the cluster. We report two simultaneous detections of the redback PSR J1824-2452I (M28I) and its X-ray counterpart. We discover a double-peaked X-ray orbital flux modulation in M28I during its pulsar state, centered around pulsar inferior conjunction. We analyze the spectrum of the quiescent neutron star low-mass X-ray binary to constrain its mass and radius. Using both hydrogen and helium NS atmosphere models, we find a NS radius of R = 9.5-11.5 km and R = 13.5 - 16.7 km, respectively, for a neutron star mass of 1.4 Msun. We also search for long-term variability in the 46 brightest X-ray sources and report the discovery of six new variable low luminosity X-ray sources in M28., Comment: 26 pages, Accepted for publication in The Astrophysical Journal

Transitional millisecond pulsars (tMSPs) have emerged in the last decade as a unique class of neutron stars at the crossroads between accretion- and rotation-powered phenomena. In their (sub-luminous) accretion disk state, with X-ray luminosities of order $10^{33}-10^{34}$ erg s$^{-1}$, they switch rapidly between two distinct X-ray modes: the disk-high (DH) and disk-low (DL) states. We present a systematic XMM-Newton and Chandra analysis of the aperiodic X-ray variability of all three currently known tMSPs, with a main focus on their disk state and separating DH and DL modes. We report the discovery of flat-topped broadband noise in the DH state of two of them, with break frequencies of 2.8 mHz (PSR J1023+0038) and 0.86 mHz (M28-I). We argue that the lowest frequency variability is similar to that seen in disk-accreting X-ray binaries in the hard state, at typical luminosities at least 2 orders of magnitude higher than tMSPs. We find strong variability in the DH state around 1 Hz, not typical of hard state X-ray binaries, with fractional rms amplitudes close to 30%. We discuss our results and use them to constrain the properties of the accretion disk, assuming that the X-ray variability is produced by fluctuations in mass accretion rate, and that the break frequency corresponds to the viscous timescale at the inner edge of the disk. In this context, we find that the newly found break frequencies are broadly consistent with a disk truncated close to the light cylinder with $\dot{M}\simeq10^{13}-5\times10^{14}$ g s$^{-1}$ and a viscosity parameter $\alpha\gtrsim$0.2., Comment: MNRAS, accepted. 10 pages, 7 figures, 2 tables

A new population of neutron stars has emerged during the last decade: compact binary millisecond pulsars (CBMSPs). Because these pulsars and their companion stars are in tight orbits with typical separations of $10^{11}$ cm, their winds interact strongly forming an intrabinary shock. Electron-positron pairs reaccelerated at the shock can reach energies of about 10 TeV, which makes this new population a potential source of GeV-TeV cosmic ray positrons. We present an analytical model for the fluxes and spectra of positrons from intrabinary shocks of CBMSPs. We find that the minimum energy $E_{\min}$ of the pairs that enter the shock is critical to quantify the energy spectrum with which positrons are injected into the interstellar medium. We measure for the first time the Galactic scale height of CBMSPs, $z_e=0.4\pm0.1$ kpc, after correcting for an observational bias against finding them close to the Galactic plane. From this, we estimate a local density of 5-9 kpc$^{-3}$ and an extrapolated total of 2-7 thousand CBMSPs in the Galaxy. We then propagate the pairs in the isotropic diffusion approximation and find that the positron flux from the total population is about two times higher than that from the 52 currently known systems. For $E_{\min}$ between 1 and 50 GeV, our model predicts only a minor contribution from CBMSPs to the diffuse positron flux at 100 GeV observed at Earth. We also quantify the effects of anisotropic transport due to the ordered Galactic magnetic field, which can change the diffuse flux from nearby sources drastically. Finally, we find that a single "hidden" CBMSP close to the Galactic plane can yield a positron flux comparable to the AMS-02 measurements at 600 GeV if its line-of-sight to Earth is along the ordered Galactic field lines, while its combined electron and positron flux at higher energies would be close to the measurements of CALET, DAMPE and Fermi-LAT., Comment: Matches version published in JCAP, after minor corrections. Added acknowledgements to the PHAROS COST Action (CA16214) and to the Spanish Ministry of Science, Innovation and Universities (grant CAS19/00228). 31 pages, 10 figures, 5 tables

The maximum mass of a neutron star has important implications across multiple research fields, including astrophysics, nuclear physics and gravitational wave astronomy. Compact binary millisecond pulsars (with orbital periods shorter than about a day) are a rapidly-growing pulsar population, and provide a good opportunity to search for the most massive neutron stars. Applying a new method to measure the velocity of both sides of the companion star, we previously found that the compact binary millisecond pulsar PSR J2215+5135 hosts one of the most massive neutron stars known to date, with a mass of 2.27$\pm$0.16 M$_\odot$ (Linares, Shahbaz & Casares, 2018). We reexamine the properties of the 0.33 M$_\odot$ companion star, heated by the pulsar, and argue that irradiation in this "redback" binary is extreme yet stable, symmetric and not necessarily produced by an extended source. We also review the neutron star mass distribution in light of this and more recent discoveries. We compile a list of all (nine) systems with published evidence for super-massive neutron stars, with masses above 2 M$_\odot$. We find that four of them are compact binary millisecond pulsars (one black widow, two redbacks and one redback candidate). This shows that compact binary millisecond pulsars are key to constraining the maximum mass of a neutron star., Comment: Invited review article, submitted for publication in the proceedings of the XIII Multifrequency Behaviour of High Energy Cosmic Sources Workshop (June 3-8, 2019, Palermo, Italy); 12 pages, 1 table, 4 figures

In this White Paper we present the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to Observatory Science targets. These include flaring stars, supernova remnants, accreting white dwarfs, low and high mass X-ray binaries, radio quiet and radio loud active galactic nuclei, tidal disruption events, and gamma-ray bursts. eXTP will be excellently suited to study one common aspect of these objects: their often transient nature. Developed by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s., Comment: Accepted for publication on Sci. China Phys. Mech. Astron. (2019)

In this White Paper we present the potential of the Enhanced X-ray Timing and Polarimetry (eXTP) mission for determining the nature of dense matter; neutron star cores host an extreme density regime which cannot be replicated in a terrestrial laboratory. The tightest statistical constraints on the dense matter equation of state will come from pulse profile modelling of accretion-powered pulsars, burst oscillation sources, and rotation-powered pulsars. Additional constraints will derive from spin measurements, burst spectra, and properties of the accretion flows in the vicinity of the neutron star. Under development by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s., Comment: Accepted for publication on Sci. China Phys. Mech. Astron. (2019)

New millisecond pulsars (MSPs) in compact binaries provide a good opportunity to search for the most massive neutron stars. Their main-sequence companion stars are often strongly irradiated by the pulsar, displacing the effective center of light from their barycenter and making mass measurements uncertain. We present a series of optical spectroscopic and photometric observations of PSR J2215+5135, a "redback" binary MSP in a 4.14 hr orbit, and measure a drastic temperature contrast between the dark/cold (T$_\mathrm{N}$=5660$^{+260}_{-380}$ K) and bright/hot (T$_\mathrm{D}$=8080$^{+470}_{-280}$ K) sides of the companion star. We find that the radial velocities depend systematically on the atmospheric absorption lines used to measure them. Namely, the semi-amplitude of the radial velocity curve of J2215 measured with magnesium triplet lines is systematically higher than that measured with hydrogen Balmer lines, by 10%. We interpret this as a consequence of strong irradiation, whereby metallic lines dominate the dark side of the companion (which moves faster) and Balmer lines trace its bright (slower) side. Further, using a physical model of an irradiated star to fit simultaneously the two-species radial velocity curves and the three-band light curves, we find a center-of-mass velocity of K$_2$=412.3$\pm$5.0 km s$^{-1}$ and an orbital inclination i=63.9$^\circ$$^{+2.4}_{-2.7}$. Our model is able to reproduce the observed fluxes and velocities without invoking irradiation by an extended source. We measure masses of M$_1$=2.27$^{+0.17}_{-0.15}$ M$_\odot$ and M$_2$=0.33$^{+0.03}_{-0.02}$ M$_\odot$ for the neutron star and the companion star, respectively. If confirmed, such a massive pulsar would rule out some of the proposed equations of state for the neutron star interior., Comment: Published in The Astrophysical Journal on May 23d, 2018. 14 pages, 13 figures, 2 tables. See animated video at https://youtu.be/l41y_GqJq7c

We present the discovery of the variable optical and X-ray counterparts to the radio millisecond pulsar (MSP) PSR J1306-40, recently discovered by Keane et al. We find that both the optical and X-ray fluxes are modulated with the same period, which allows us to measure for the first time the orbital period P$_{\rm orb}$=1.09716[6] d. The optical properties are consistent with a main sequence companion with spectral type G to mid K and, together with the X-ray luminosity (8.8$\times$10$^{31}$ erg s$^{-1}$ in the 0.5-10 keV band, for a distance of 1.2 kpc), confirm the redback classification of this pulsar. Our results establish the binary nature of PSR J1306-40, which has the longest P$_{\rm orb}$ among all known compact binary MSPs in the Galactic disk. We briefly discuss these findings in the context of irradiation and intrabinary shock emission in compact binary MSPs., Comment: 5 pages, 3 figures, MNRAS, accepted (minor changes after referee report). In memory of our esteemed colleague Enrique Garc\'ia-Berro

We have monitored the transient neutron star low-mass X-ray binary 1RXS J180408.9-342058 in quiescence after its ~4.5 month outburst in 2015. The source has been observed using Swift and XMM-Newton. Its X-ray spectra were dominated by a thermal component. The thermal evolution showed a gradual X-ray luminosity decay from ~18 x 10^32 to ~4 x 10^32 (D/5.8 kpc)^2 erg s^{-1} between ~8 and ~379 d in quiescence, and the inferred neutron star surface temperature (for an observer at infinity; using a neutron star atmosphere model) decreased from ~100 to ~71 eV. This can be interpreted as cooling of an accretion-heated neutron star crust. Modelling the observed temperature curve (using NSCOOL) indicated that the source required ~1.9 MeV per accreted nucleon of shallow heating in addition to the standard deep crustal heating to explain its thermal evolution. Alternatively, the decay could also be modelled without the presence of deep crustal heating, only having a shallow heat source (again ~1.9 MeV per accreted nucleon was required). However, the XMM-Newton data statistically required an additional power-law component. This component contributed ~30 per cent of the total unabsorbed flux in 0.5-10 keV energy range. The physical origin of this component is unknown. One possibility is that it arises from low-level accretion. The presence of this component in the spectrum complicates our cooling crust interpretation because it might indicate that the smooth luminosity and temperature decay curves we observed may not be due to crust cooling but due to some other process., Comment: 10 pages, 4 figures, 3 tables. Accepted in MNRAS

We present the discovery of a variable optical counterpart to the unidentified gamma-ray source 3FGL J0212.1+5320, and argue this is a new compact binary millisecond pulsar (MSP) candidate. We show 3FGL J0212.1+5320 hosts a semi-detached binary with a 0.86955$\pm$0.00015 d orbital period and a F6-type companion star at an estimated distance of D=1.1$\pm$0.2 kpc, with a radial velocity curve semi-amplitude K$_2$=214.1$\pm$5.0 km s$^{-1}$ and a projected rotational velocity of Vsin(i)=73.2$\pm$1.6 km s$^{-1}$. We find a hard X-ray source at the same location with a 0.5$-$10 keV luminosity L$_\mathrm{X}$=2.6$\times$10$^{32}$ (D/1.1 kpc)$^2$ erg s$^{-1}$, which strengthens the MSP identification. Our results imply a mass ratio q=M$_2$/M$_1$=0.26$^{+0.02}_{-0.03}$ if the companion star fills its Roche lobe, and q$\gtrsim$0.26 in any case. This classifies 3FGL J0212.1+5320 as a "redback" binary MSP; if its MSP nature is confirmed, this will be the brightest compact binary MSP in the optical band (r'$\simeq$14.3 mag) and will have the longest orbital period among Galactic field systems (nearly 21 hr). Based on the light curve peak-to-peak amplitude ($\Delta$r=0.19 mag), we further suggest that the orbital inclination is high and the putative pulsar mass is close to canonical (M$_1$$\simeq$1.3$-$1.6 M$_\odot$). Finally, we discuss the lack of heating signatures and asymmetric optical light curves in the context of other redback MSPs., Comment: Accepted for publication in MNRAS; 9 pages, 8 figures, minor changes after referee suggestions

Tracking the spectral evolution of transiently accreting neutron stars between outburst and quiescence probes relatively poorly understood accretion regimes. Such studies are challenging because they require frequent monitoring of sources with luminosities below the thresholds of current all-sky X-ray monitors. We present the analysis of over 30 observations of the neutron star low-mass X-ray binary SAX J1750.8-2900 taken across four years with the X-ray telescope aboard Swift. We find spectral softening with decreasing luminosity both on long ($\sim$1 year) and short ($\sim$days to week) timescales. As the luminosity decreases from $4\times10^{36}$ erg s$^{-1}$ to $ \sim1\times10^{35} $ erg s$^{-1}$ (0.5-10 keV), the power law photon index increases from from 1.4 to 2.9. Although not statistically required, our spectral fits allow an additional soft component that displays a decreasing temperature as the luminosity decreases from $4 \times 10^{36} $ to $6 \times 10^{34}$ erg s$^{-1}$. Spectral softening exhibited by SAX J1750.8-2900 is consistent both with accretion emission whose spectral shape steepens with decreasing luminosity and also with being dominated by a changing soft component, possibly associated with accretion onto the neutron star surface, as the luminosity declines., Comment: 11 pages, 9 figures, 3 Tables. Resubmitted to The Astrophysical Journal after incorporating comments from the referee

We present Chandra, XMM-Newton, and Swift observations of the quiescent neutron star in the transient low-mass X-ray binary MAXI J0556-332. Observations of the source made during outburst (with the Rossi X-ray Timing Explorer) reveal tracks in its X-ray color--color and hardness--intensity diagrams that closely resemble those of the neutron-star Z sources, suggesting that MAXI J0556-332 had near- or super-Eddington luminosities for a large part of its ~16 month outburst. A comparison of these diagrams with those of other Z sources suggests a source distance of 46+/-15 kpc. Fits to the quiescent spectra of MAXI J0556-332 with a neutron-star atmosphere model (with or without a power-law component) result in distance estimates of 45+/-3 kpc, for a neutron-star radius of 10 km and a mass of 1.4 Msun. The spectra show the effective surface temperature of the neutron star decreasing monotonically over the first ~500 days of quiescence, except for two observations that were likely affected by enhanced low-level accretion. The temperatures we obtain for the fits that include a power-law (kT_eff=184-308 eV) are much higher than those seen for any other neutron star heated by accretion, while the inferred cooling (e-folding) timescale (~200 days) is similar to other sources. Fits without a power-law yield higher temperatures (kT_eff=190-336 eV) and a shorter e-folding time (~160 days). Our results suggest that the heating of the neutron-star crust in MAXI J0556-332 was considerably more efficient than for other systems, possibly indicating additional or more efficient shallow heat sources in its crust., Comment: Accepted for publication: 2014, ApJ, 795, 131 (13 pages)

Using the High Resolution Camera onboard the Chandra X-ray Observatory, we have measured an accurate position for the bright persistent neutron-star X-ray binary and atoll source GX 3+1. At a location that is consistent with this new position we have discovered the near-infrared (NIR) counterpart to GX 3+1 in images taken with the PANIC and FourStar cameras on the Magellan Baade Telescope. The identification of this K_s=15.8+-0.1 mag star as the counterpart is based on the presence of a Br-gamma emission line in a NIR spectrum taken with the FIRE spectrograph on the Baade Telescope. The absolute magnitude derived from the best available distance estimate to GX 3+1 indicates that the mass donor in the system is not a late-type giant. We find that the NIR light in GX 3+1 is likely dominated by the contribution from a heated outer accretion disk. This is similar to what has been found for the NIR flux from the brighter class of Z sources, but unlike the behavior of atolls fainter (Lx ~ 1e36 to 1e37 erg/s) than GX 3+1, where optically-thin synchrotron emission from a jet probably dominates the NIR flux., Comment: Accepted for publication in ApJ

In this work, the analysis of multi-epoch (1995-2010) X-ray observations of the Seyfert 1 galaxy H0557-385 is presented. The wealth of data presented in this analysis show that the source exhibits dramatic spectral variability, from a typical unabsorbed Seyfert 1 type spectrum to a Compton-thin absorbed state, on time scales of ~5 years. This extreme change in spectral shape can be attributed to variations in the column density and covering fraction of a neutral absorbing medium attenuating the emission from the central continuum source. Evidence for Compton reflection of the intrinsic nuclear emission is present in each of the spectra, though this feature is most prominent in the low-state spectra, where the associated Fe emission line complex is clearly visible. In addition to the variable absorbing medium, a warm absorber component has been detected in each spectral state. Optical spectroscopy concurrent with the 2010 XMM-Newton observation campaign have detected the presence of broad optical emission lines during an X-ray absorption event.From the analysis of both X-ray and optical spectroscopic data, it has been inferred that the X-ray spectral variability is a result of obscuration of the central emission region by a clumpy absorber covering >80 per cent of the source with an average column density of NH ~7x10^{23} cm^{-2}, and which is located outside the broad line region at a distance from the central source consistent with the dust sublimation radius of the AGN., Comment: Accepted for publication on MNRAS, 15 pages, 7 figures

Compact binary millisecond pulsars with main-sequence donors, often referred to as "redbacks", constitute the long-sought link between low-mass X-ray binaries and millisecond radio pulsars, and offer a unique probe of the interaction between pulsar winds and accretion flows. We present a systematic study of eight nearby redbacks, using more than 100 observations obtained with Swift's X-ray Telescope. We distinguish between three main states: pulsar, disk and outburst states. We find X-ray mode switching in the disk state of PSR J1023+0038 and XSS J12270-4859, similar to what was found in the other redback which showed evidence for accretion: rapid, recurrent changes in X-ray luminosity (0.5-10 keV, L$_\mathrm{X}$), between [6-9]$\times$10$^{32}$ erg s$^{-1}$ (disk-passive state) and [3-5]$\times$10$^{33}$ erg s$^{-1}$ (disk-active state). This strongly suggests that mode switching $-$which has not been observed in quiescent low-mass X-ray binaries$-$ is universal among redback millisecond pulsars in the disk state. We briefly explore the implications for accretion disk truncation, and find that the inferred magnetospheric radius in the disk state of PSR J1023+0038 and XSS J12270-4859 lies outside the light cylinder. Finally, we note that all three redbacks which have developed accretion disks have relatively high L$_\mathrm{X}$ in the pulsar state ($>$10$^{32}$ erg s$^{-1}$)., Comment: Accepted for publication in The Astrophysical Journal. Minor revision after referee's comments, including a brief discussion of the different variability seen in SAX J1808.4-3658. 10 pages, 5 figures, 1 table

Both the broad iron (Fe) line and the frequency of the kilohertz quasi-periodic oscillations (kHz QPOs) in neutron star low-mass X-ray binaries (LMXBs) can potentially provide independent measures of the inner radius of the accretion disc. We use XMM-Newton and simultaneous Rossi X-ray Timing Explorer observations of the LMXB 4U 1636-53 to test this hypothesis. We study the properties of the Fe-K emission line as a function of the spectral state of the source and the frequency of the kHz QPOs. We find that the inner radius of the accretion disc deduced from the frequency of the upper kHz QPO varies as a function of the position of the source in the colour-colour diagram, in accordance with previous work and with the standard scenario of accretion disc geometry. On the contrary, the inner disc radius deduced from the profile of the Fe line is not correlated with the spectral state of the source. The values of the inner radius inferred from kHz QPOs and Fe lines, in four observations, do not lead to a consistent value of the neutron star mass, regardless of the model used to fit the Fe line. Our results suggest that either the kHz QPO or the standard relativistic Fe line interpretation does not apply for this system. Furthermore, the simultaneous detection of kHz QPOs and broad Fe lines is difficult to reconcile with models in which the broadening of the Fe line is due to the reprocessing of photons in an outflowing wind., Comment: 12 pages, 6 figures, Accepted for publication in MNRAS

We report and study the outburst of a new transient X-ray binary (XRB) in Terzan 5, the third detected in this globular cluster, Swift J174805.3-244637 or Terzan 5 X-3. We find clear spectral hardening in Swift/XRT data during the outburst rise to the hard state, thanks to our early coverage (starting at L_X ~ 4x10^{34} ergs/s) of the outburst. This hardening appears to be due to the decline in relative strength of a soft thermal component from the surface of the neutron star (NS) during the rise. We identify a {Type I X-ray burst} in Swift/XRT data with a long (16 s) decay time, indicative of {hydrogen burning on the surface of the} NS. We use Swift/BAT, Maxi/GSC, Chandra/ACIS, and Swift/XRT data to study the spectral changes during the outburst, identifying a clear hard-to-soft state transition. We use a Chandra/ACIS observation during outburst to identify the transient's position. Seven archival Chandra/ACIS observations show evidence for variations in Terzan 5 X-3's non-thermal component, but not the thermal component, during quiescence. The inferred long-term time-averaged mass accretion rate, from the quiescent thermal luminosity, suggests that if this outburst is typical and only slow cooling processes are active in the neutron star core, such outbursts should recur every ~10 years., Comment: ApJ in press

The X-ray transient IGR J18245-2452 in the globular cluster M28 contains the first neutron star (NS) seen to switch between rotation-powered and accretion-powered pulsations. We analyse its 2013 March-April 25d-long outburst as observed by Swift, which had a peak bolometric luminosity of ~6% of the Eddington limit (L$_{E}$), and give detailed properties of the thermonuclear burst observed on 2013 April 7. We also present a detailed analysis of new and archival Chandra data, which we use to study quiescent emission from IGR J18245-2452 between 2002 and 2013. Together, these observations cover almost five orders of magnitude in X-ray luminosity (L$_X$, 0.5-10 keV). The Swift spectrum softens during the outburst decay (photon index $\Gamma$ from 1.3 above L$_X$/L$_{E}$=10$^{-2}$ to ~2.5 at L$_X$/L$_{E}$=10$^{-4}$), similar to other NS and black hole (BH) transients. At even lower luminosities, deep Chandra observations reveal hard ($\Gamma$=1-1.5), purely non-thermal and highly variable X-ray emission in quiescence. We therefore find evidence for a spectral transition at L$_X$/L$_{E}$~10$^{-4}$, where the X-ray spectral softening observed during the outburst decline turns into hardening as the source goes to quiescence. Furthermore, we find a striking variability pattern in the 2008 Chandra light curves: rapid switches between a high-L$_X$ "active" state (L$_X\simeq$3.9x10$^{33}$ erg/s) and a low-L$_X$ "passive" state (L$_X\simeq$5.6x10$^{32}$ erg/s), with no detectable spectral change. We put our results in the context of low luminosity accretion flows around compact objects and X-ray emission from millisecond radio pulsars. Finally, we discuss possible origins for the observed mode switches in quiescence, and explore a scenario where they are caused by fast transitions between the magnetospheric accretion and pulsar wind shock emission regimes., Comment: Accepted for publication in MNRAS. 12 pages, 6 figures, 2 tables

We report the discovery of kHz quasi periodic oscillations (QPOs) in three Rossi X-ray Timing Explorer observations of the low mass X-ray binary (LMXB) XTE J1701-407. In one of the observations we detect a kHz QPO with a characteristic frequency of 1153 +/- 5 Hz, while in the other two observations we detect twin QPOs at characteristic frequencies of 740 +/- 5 Hz, 1112 +/- 17 Hz and 740 +/- 11 Hz, 1098 +/- 5 Hz. All detections happen when XTE J1701-407 was in its high intensity soft state, and their single trial significance are in the 3.1-7.5 sigma range. The frequency difference in the centroid frequencies of the twin kHz QPOs (385 +/- 13 Hz) is one of the largest seen till date. The 3-30 keV fractional rms amplitude of the upper kHz QPO varies between ~18 % and ~30 %. XTE J1701-407, with a persistent luminosity close to 1 % of the Eddington limit, is among the small group of low luminosity kHz QPO sources and has the highest rms for the upper kHz QPO detected in any source. The X-ray spectral and variability characteristics of this source indicate its atoll source nature., Comment: 10 pages, 9 figures, accepted for publication in MNRAS

We analysed the X-ray spectra of six observations, simultaneously taken with XMM-Newton and Rossi X-ray Timing Explorer (RXTE), of the neutron star low-mass X-ray binary 4U 1636-53. The observations cover several states of the source, and therefore a large range of inferred mass accretion rate. These six observations show a broad emission line in the spectrum at around 6.5 keV, likely due to iron. We fitted this line with a set of phenomenological models of a relativistically broadened line, plus a model that accounts for relativistically smeared and ionised reflection from the accretion disc. The latter model includes the incident emission from both the neutron-star surface or boundary layer and the corona that is responsible for the high-energy emission in these systems. From the fits with the reflection model we found that in four out of the six observations the main contribution to the reflected spectrum comes from the neutron-star surface or boundary layer, whereas in the other two observations the main contribution to the reflected spectrum comes from the corona. We found that the relative contribution of these two components is not correlated to the state of the source. From the phenomenological models we found that the iron line profile is better described by a symmetric, albeit broad, profile. The width of the line cannot be explained only by Compton broadening, and we therefore explored the case of relativistic broadening. We further found that the direct emission from the disc, boundary layer, and corona generally evolved in a manner consistent with the standard accretion disc model, with the disc and boundary layer becoming hotter and the disc moving inwards as the source changed from the hard in to the soft state. The iron line, however, did not appear to follow the same trend., Comment: 20 pages, 8 figures, Accepted for publication in MNRAS

We report the discovery of burst oscillations from the intermittent accretion-powered millisecond pulsar (AMP) HETE J1900.1-2455, with a frequency approximately 1 Hz below the known spin frequency. The burst oscillation properties are far more similar to those of the non-AMPs and Aql X-1 (an intermittent AMP with a far lower duty cycle), than those of the AMPs SAX J1808.4-3658 and XTE J1814-338. We discuss the implications for models of the burst oscillation and intermittency mechanisms., Comment: 4 pages, 4 figures, submitted to ApJ Letters

XMM-Newton observed the accreting millisecond pulsar SAX J1808.4-3658 during its 2008 outburst. We present timing and spectral analyses of this observation, in particular the first pulse profile study below 2 keV, and the high-resolution spectral analysis of this source during the outburst. Combined spectral and pulse profile analyses suggest the presence of a strong unpulsed source below 2 keV that strongly reduces the pulsed fraction and a hard pulsed component that generates markedly double peaked profiles at higher energies. We also studied the high-resolution grating spectrum of SAX J1808.4-3658, and found several absorption edges and Oxygen absorption lines with whom we infer, in a model independent way, the interstellar column densities of several elements in the direction of SAX J1808.4-3658., Comment: MNRAS Letters, accepted

At present, ten years after they were first discovered, ten accreting millisecond pulsars are known. I present a study of the aperiodic X-ray variability in three of these systems, which led to the discovery of simultaneous kHz quasi periodic oscillations in XTE J1807-294 and extremely strong broadband noise at unusually low variability frequencies in IGR J00291+5934. Furthermore, we classified SWIFT J1756.9-2508 as an atoll source and measured in its 2007 outburst spectral and variability properties typical of the extreme island state. I also give detailed estimates of the total fluence during the studied outbursts., Comment: To be published in the proceedings of the workshop 'A Decade of Accreting Millisecond X-ray pulsars' (Amsterdam 14-18 April 2008; Eds. Wijnands et al.)

XTE J1701-407 is a new transient X-ray source discovered on June 8th, 2008. More than one month later it showed a rare type of thermonuclear explosion: a long type I X-ray burst. We report herein the results of our study of the spectral and flux evolution during this burst, as well as the analysis of the outburst in which it took place. We find an upper limit on the distance to the source of 6.1 kpc by considering the maximum luminosity reached by the burst. We measure a total fluence of 3.5*10^{-6} erg/cm^2 throughout the ~20 minutes burst duration and a fluence of 2.6*10^{-3} erg/cm^2 during the first two months of the outburst. We show that the flux decay is best fitted by a power law (index ~1.6) along the tail of the burst. Finally, we discuss the implications of the long burst properties, and the presence of a second and shorter burst detected by Swift ten days later, for the composition of the accreted material and the heating of the burning layer., Comment: MNRAS-Letters, accepted. Minor changes according to referee's report. 5 pages, 3 figures

SWIFT J1756.9-2508 is one of the few accreting millisecond pulsars (AMPs) discovered to date. We report here the results of our analysis of its aperiodic X-ray variability, as measured with the Rossi X-ray Timing Explorer during the 2007 outburst of the source. We detect strong (~35%) flat-topped broadband noise throughout the outburst with low characteristic frequencies (~0.1 Hz). This makes SWIFT J1756.9-2508 similar to the rest of AMPs and to other low luminosity accreting neutron stars when they are in their hard states, and enables us to classify this AMP as an atoll source in the extreme island state. We also find a hard tail in its energy spectrum extending up to 100 keV, fully consistent with such source and state classification., Comment: Accepted for publication in The Astrophysical Journal. Scheduled for the ApJ May 10, 2008, v678n 2 issue. 6 pages, 5 figures. Uses \emulateapj

We report on the X-ray aperiodic timing analysis of two accreting millisecond pulsars: XTE J1807-294 and IGR J00291+5934. On the one hand, we discovered in XTE J1807-294 seven pairs of simultaneous kilohertz quasi-periodic oscillations (kHz QPOs) separated in frequency by nearly the spin frequency of the neutron star. This confirms the suspected dichotomy in the frequency separation of kHz QPOs: sometimes once and sometimes half the spin frequency. On the other hand, we found an extreme behavior in the power spectra of IGR J00291+5934: very strong variability at very low frequencies. Namely, the fractional amplitude of the variability was ~50%, the highest value found so far in a neutron star system., Comment: To appear in the proceedings of the conference `The Multicoloured Landscape of Compact Objects and their Explosive Origins' (Cefalu, Sicily, June 2006), AIP

IGR J00291+5934 is one of the seven accreting millisecond pulsars (AMPs) discovered so far. We report on the aperiodic timing and color analysis of its X-ray flux, using all the RXTE observations of the 2004 outburst. Flat-top noise and two harmonically related quasi-periodic oscillations, all of them at very low frequencies (0.01-0.1 Hz), were present in the power spectra during most of the outburst as well as a very high fractional variability (~50%). These properties are atypical not only for AMPs but also for neutron star low-mass X-ray binaries (LMXBs) in general. There are instead some remarkable similarities with the variability observed in black hole systems, reinforcing the connections between these two types of LMXB, as well as some interesting differences. We note finally that the results of this paper are difficult to reconcile with interpretations where any break frequency of power density spectra scales inversely with the mass of the central object at an accuracy sufficient to distinguish between the masses of neutron stars and black holes in LMXBs., Comment: Accepted by The Astrophysical Journal. Scheduled for ApJ May 1, 2007, v660n1 issue. Minor changes according to referee's comments

We present Chandra and XMM-Newton X-ray observations that monitor the neutron star cooling of the quasi-persistent neutron star X-ray transients KS 1731-260 and MXB 1659-29 for approximately 4 years after these sources returned to quiescence from prolonged outbursts. In both sources the outbursts were long enough to significantly heat the neutron star crust out of thermal equilibrium with the core. We analyse the X-ray spectra by fitting absorbed neutron star atmosphere models to the observations. The results of our analysis strengthen the preliminary findings of Wijnands et al. that in both sources the neutron star crust cools down very rapidly suggesting it has a high heat conductivity and that the neutron star core requires enhanced core cooling processes. Importantly, we now detect the flattening of the cooling in both sources as the crust returns to thermal equilibrium with the core. We measure the thermal equilbrium flux and temperature in both sources by fitting a curve that decays exponentially to a constant level. The cooling curves cannot be fit with just a simple exponential decay without the constant offset. We find the constant bolometric flux and effective temperature components to be (9.2 +/- 0.9)E-14 ergs/cm^2/s and 70.0 +/- 1.6 eV in KS 1731-260 and (1.7 +/- 0.3)E-14 ergs/cm^2/s and 51.6 +/- 1.4 eV in MXB 1659-29. In addition, we find that the crust of KS 1731-260 cools faster than that of MXB 1659-29 by a factor of ~2, likely due to different crustal properties. This is the first time that the cooling of a neutron star crust into thermal equilibrium with the core has been observed in such detail., Comment: 11 pages, 3 figures, 4 tables, accepted for publication in MNRAS

We report the discovery of twin kHz QPOs in the X ray flux of XTE J1807--294, the fourth accreting millisecond pulsar (AMP). This is the second AMP exhibiting twin kHz QPOs. In contrast to the first case, SAX J1808.4--3658, the frequency separation $\Delta\nu$ between the kHz QPOs is consistent with the pulse frequency (190.6 Hz), not with half that value, confirming for the first time from pulsation measurements the inference, based on burst oscillations, that 'slow rotators' (spin frequency less than 400 Hz) have $\Delta\nu$ approximately equal to the spin frequency. While the QPOs move in frequency together over a range of more than 200 Hz, $\Delta\nu$ remains constant with an average value of 205$\pm$6 Hz. Variability components were found in the 5--130 Hz range similar to those seen in other LMXBs. The correlations between the QPO and noise frequencies are also similar to those in other sources, but shifted by a factor of 1.59 in kHz QPO frequencies, similar to the factor 1.45 shift found for SAX J1808.4--3658. Our results argue in favour of a spin-related formation mechanism for twin kHz QPOs and against a spin-related cause of the shift in the frequency correlations., Comment: Accepted for publication in ApJ (to appear in 1 December 2005, v634, 2 issue)