We present NuSTAR observations of the luminous neutron star low-mass X-ray binary (NS LMXB) and Z source GX 5-1. During our three observations made with separations of roughly two days, the source traced out an almost complete Z track. We extract spectra from the various branches and fit them with a continuum model that has been successfully applied to other Z sources. Surprisingly, and unlike most of the (luminous) NS-LMXBs observed with NuSTAR, we do not find evidence for reflection features in any of the spectra of GX 5-1. We discuss several possible explanations for the absence of reflection features. Based on a comparison with other accreting neutron star systems and given the high luminosity of GX 5-1 (~1.6-2.3 times the Eddington luminosity, for a distance of 9 kpc), we consider a highly ionized disk the most likely explanation for the absence of reflection features in GX 5-1.

The transient neutron star (NS) low-mass X-ray binary MAXI J0556$-$332 provides a rare opportunity to study NS crust heating and subsequent cooling for multiple outbursts of the same source. We examine {\it MAXI}, {\it Swift}, {\it Chandra}, and {\it XMM-Newton} data of MAXI J0556$-$332 obtained during and after three accretion outbursts of different durations and brightness. We report on new data obtained after outburst III. The source has been tracked up to $\sim$1800 d after the end of outburst I. Outburst I heated the crust strongly, but no significant reheating was observed during outburst II. Cooling from $\sim$333 eV to $\sim$146 eV was observed during the first $\sim$1200 d. Outburst III reheated the crust up to $\sim$167 eV, after which the crust cooled again to $\sim$131 eV in $\sim$350 d. We model the thermal evolution of the crust and find that this source required a different strength and depth of shallow heating during each of the three outbursts. The shallow heating released during outburst I was $\sim$17 MeV nucleon$^{-1}$ and outburst III required $\sim$0.3 MeV nucleon$^{-1}$. These cooling observations could not be explained without shallow heating. The shallow heating for outburst II was not well constrained and could vary from $\sim$0--2.2 MeV nucleon$^{-1}$, i.e., this outburst could in principle be explained without invoking shallow heating. We discuss the nature of the shallow heating and why it may occur at different strengths and depths during different outbursts.

Recent work on jets and disk winds in low-mass X-ray binaries (LMXBs) suggests that they are to a large extent mutually exclusive, with jets observed in spectrally hard states and disk winds observed in spectrally soft states. In this paper we use existing literature on jets and disk winds in the luminous neutron star (NS) LMXB GX 13+1, in combination with archival Rossi X-ray Timing Explorer data, to show that this source is likely able to produce jets and disk winds simultaneously. We find that jets and disk winds occur in the same location on the source's track in its X-ray color-color diagram. A further study of literature on other luminous LMXBs reveals that this behavior is more common, with indications for simultaneous jets and disk winds in the black hole LMXBs V404 Cyg and GRS 1915+105 and the NS LMXBs Sco X-1 and Cir X-1. For the three sources for which we have the necessary spectral information, we find that the simultaneous jets/winds all occur in their spectrally hardest states. Our findings indicate that in LMXBs with luminosities above a few tens of percent of the Eddington luminosity, jets and disk winds are not mutually exclusive, and that the presence of disk winds does not necessarily result in jet suppression.

We report on a detailed analysis of the so-called ~1 Hz quasi-periodic oscillation (QPO) in the eclipsing and dipping neutron-star low-mass X-ray binary EXO 0748-676. This type of QPO has previously been shown to have a geometric origin. Our study focuses on the evolution of the QPO as the source moves through the color-color diagram, in which it traces out an atoll-source-like track. The QPO frequency increases from ~0.4 Hz in the hard state to ~25 Hz as the source approaches the soft state. Combining power spectra based on QPO frequency reveals additional features that strongly resemble those seen in non-dipping/eclipsing atoll sources. We show that the low-frequency QPOs in atoll sources and the ~1 Hz QPO in EXO 0748-676 follow similar relations with respect to the noise components in their power spectra. We conclude that the frequencies of both types of QPOs are likely set by (the same) precession of a misaligned inner accretion disk. For high-inclination systems, like EXO 0748-676, this results in modulations of the neutron-star emission due to obscuration or scattering, while for lower-inclination systems the modulations likely arise from relativistic Doppler boosting and light-bending effects.

The X-ray transient XTE J1701-462 was the first source observed to evolve through all known subclasses of low-magnetic-field neutron star low-mass X-ray binaries (NS-LMXBs), as a result of large changes in its mass accretion rate. To investigate to what extent similar evolution is seen in other NS-LMXBs we have performed a detailed study of the color-color and hardness-intensity diagrams (CDs and HIDs) of Cyg X-2, Cir X-1, and GX 13+1 -- three luminous X-ray binaries, containing weakly magnetized neutron stars, known to exhibit strong secular changes in their CD/HID tracks. Using the full set of Rossi X-ray Timing Explorer Proportional Counter Array data collected for the sources over the 16 year duration of the mission, we show that Cyg X-2 and Cir X-1 display CD/HID evolution with close similarities to XTE J1701-462. Although GX 13+1 shows behavior that is in some ways unique, it also exhibits similarities to XTE J1701-462, and we conclude that its overall CD/HID properties strongly indicate that it should be classified as a Z source, rather than as an atoll source. We conjecture that the secular evolution of Cyg X-2, Cir X-1, and GX 13+1 -- illustrated by sequences of CD/HID tracks we construct -- arises from changes in the mass accretion rate. Our results strengthen previous suggestions that within single sources Cyg-like Z source behavior takes place at higher luminosities and mass accretion rates than Sco-like Z behavior, and lend support to the notion that the mass accretion rate is the primary physical parameter distinguishing the various NS-LMXB subclasses.

When neutron stars reside in transient X-ray binaries, their crustal layers become heated during accretion outbursts and subsequently cool in quiescence. Observing and modeling this thermal response has yielded valuable insight into the physics of neutron star crusts. However, one unresolved problem is the evidence in several neutron stars for an extra energy source, located at shallow depth in the crust, that is not accounted for by standard heating models. Its origin remains puzzling, and it is currently unclear whether this additional heating occurs in all neutron stars, and if the magnitude is always the same. Here, we report on Chandra observations that cover two years after the 2012 outburst of the transient neutron star X-ray binary Swift J174805.3-244637 in the globular cluster Terzan 5. The temperature of the neutron star was elevated during the first two months following its ~8 week accretion episode, but had decayed to the pre-outburst level within ~100 days. Interpreting this as rapid cooling of the accretion-heated crust, we model the observed temperature curve with a thermal evolution code. We find that there is no need to invoke shallow heating for this neutron star, although an extra energy release up to ~1.4 MeV/nucleon is allowed by the current data (2-sigma confidence). We also present two new data points on the crust cooling curve of the 11-Hz X-ray pulsar IGR J17480-2446 in Terzan 5, which was active in 2010. The temperature of this neutron star remained significantly above its pre-outburst level, but we detect no temperature change since the previous measurements of 2013 February. This is consistent with the slower cooling expected several years post-outburst.

A recent joint XMM-Newton/Nuclear Spectroscopic Telescope Array (NuSTAR) observation of the accreting neutron star Cen X-4 ($L_{\rm X}\sim10^{33}{\rm~erg~s}^{-1}$) revealed a hard power-law component ($\Gamma\sim1$-$1.5$) with a relatively low cut-off energy (~10 keV), suggesting bremsstrahlung emission. The physical requirements for bremsstrahlung combined with other observed properties of Cen X-4 suggest the emission comes from a boundary layer rather than the accretion flow. The accretion flow itself is thus undetected (with an upper limit of $L_{\rm flow}\lesssim0.3 L_{\rm X}$). A deep search for coherent pulsations (which would indicate a strong magnetic field) places a 6 per cent upper limit on the fractional amplitude of pulsations, suggesting the flow is not magnetically regulated. Considering the expected energy balance between the accretion flow and the boundary layer for different values of the neutron star parameters (size, magnetic field, and spin) we use the upper limit on $L_{\rm flow}$ to set an upper limit of $\varepsilon\lesssim0.3$ for the intrinsic radiative efficiency of the accretion flow for the most likely model of a fast-spinning, non-magnetic neutron star. The non-detection of the accretion flow provides the first direct evidence that this flow is indeed 'radiatively inefficient', i.e. most of the gravitational potential energy lost by the flow before it hits the star is not emitted as radiation.

The accreting millisecond X-ray pulsar SAX J1808.4--3658 shows peculiar low luminosity states known as "reflares" after the end of the main outburst. During this phase the X-ray luminosity of the source varies by up to three orders of magnitude in less than 1-2 days. The lowest X-ray luminosity observed reaches a value of ~1e32 erg/s, only a factor of a few brighter than its typical quiescent level. We investigate the 2008 and 2005 reflaring state of SAX J1808.4-3658 to determine whether there is any evidence for a change in the accretion flow with respect to the main outburst. We perform a multiwavelength photometric and spectral study of the 2005 and 2008 reflares with data collected during an observational campaign covering the near-infrared, optical, ultra-violet and X-ray band. We find that the NIR/optical/UV emission, expected to some from the outer accretion disk shows variations in luminosity which are 1--2 orders of magnitude shallower than in X-rays. The X-ray spectral state observed during the reflares does not change substantially with X-ray luminosity indicating a rather stable configuration of the accretion flow. We investigate the most likely configuration of the innermost regions of the accretion flow and we infer an accretion disk truncated at or near the co-rotation radius. We interpret these findings as due to either a strong outflow (due to a propeller effect) or a trapped disk (with limited/no outflow) in the inner regions of the accretion flow.

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.

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.

The transient neutron star low-mass X-ray binary and 11 Hz X-ray pulsar IGR J17480-2446 in the globular cluster Terzan 5 exhibited an 11-week accretion outburst in 2010. Chandra observations performed within five months after the end of the outburst revealed evidence that the crust of the neutron star became substantially heated during the accretion episode and was subsequently cooling in quiescence. This provides the rare opportunity to probe the structure and composition of the crust. Here, we report on new Chandra observations of Terzan 5 that extend the monitoring to ~2.2 yr into quiescence. We find that the thermal flux and neutron star temperature have continued to decrease, but remain significantly above the values that were measured before the 2010 accretion phase. This suggests that the crust has not thermally relaxed yet, and may continue to cool. Such behavior is difficult to explain within our current understanding of heating and cooling of transiently accreting neutron stars. Alternatively, the quiescent emission may have settled at a higher observed equilibrium level (for the same interior temperature), in which case the neutron star crust may have fully cooled.

We present new Chandra X-ray observations of the transient black hole X-ray binary MAXI J1659-152 in quiescence. These observations were made more than one year after the end of the source's 2010-2011 outburst. We detect the source at a 0.5-10 keV flux of 2.8(8)e-15 erg/cm^2/s, which corresponds to a luminosity of ~1.2e31 (d/ 6 kpc)^2 erg/s. This level, while being the lowest at which the source has been detected, is within factors of ~2 of the levels seen at the end of the initial decay of the outburst and soon after a major reflare of the source. The quiescent luminosity of MAXI J1659-152, which is the shortest-orbital-period black hole X-ray binary (~2.4 hr), is lower than that of neutron-star X-ray binaries with similar periods. However, it is higher than the quiescent luminosities found for black hole X-ray binaries with orbital periods ~2-4 times longer. This could imply that a minimum quiescent luminosity may exist for black hole X-ray binaries, around orbital periods of ~5-10 hr, as predicted by binary-evolution models for the mass transfer rate. Compared to the hard state we see a clear softening of the power-law spectrum in quiescence, from an index of 1.55(4) to an index of 2.5(4). We constrain the luminosity range in which this softening starts to fractional Eddington luminosities of (0.18-6.2)e-5 (d/ 6 kpc)^2 (M/ 8 Msun), which is consistent with the ranges inferred for other sources.

The quasi-persistent neutron star low-mass X-ray binary MXB 1659-29 went into quiescence in 2001, and we have followed its quiescent X-ray evolution since. Observations over the first 4 years showed a rapid drop in flux and temperature of the neutron star atmosphere, interpreted as cooling of the neutron star crust which had been heated during the 2.5 year outburst. However, observations taken approximately 1400 and 2400 days into quiescence were consistent with each other, suggesting the crust had reached thermal equilibrium with the core. Here we present a new Chandra observation of MXB 1659-29 taken 11 years into quiescence and 4 years since the last Chandra observation. This new observation shows an unexpected factor of ~3 drop in count rate and change in spectral shape since the last observation, which cannot be explained simply by continued cooling. Two possible scenarios are that either the neutron star temperature has remained unchanged and there has been an increase in the column density, or, alternatively the neutron star temperature has dropped precipitously and the spectrum is now dominated by a power-law component. The first scenario may be possible given that MXB 1659-29 is a near edge-on system, and an increase in column density could be due to build-up of material in, and a thickening of, a truncated accretion disk during quiescence. But, a large change in disk height may not be plausible if standard accretion disk theory holds during quiescence. Alternatively, the disk may be precessing, leading to a higher column density during this latest observation.

We present the first multi-wavelength follow-up observations of two candidate gravitational-wave (GW) transient events recorded by LIGO and Virgo in their 2009-2010 science run. The events were selected with low latency by the network of GW detectors and their candidate sky locations were observed by the Swift observatory. Image transient detection was used to analyze the collected electromagnetic data, which were found to be consistent with background. Off-line analysis of the GW data alone has also established that the selected GW events show no evidence of an astrophysical origin; one of them is consistent with background and the other one was a test, part of a "blind injection challenge". With this work we demonstrate the feasibility of rapid follow-ups of GW transients and establish the sensitivity improvement joint electromagnetic and GW observations could bring. This is a first step toward an electromagnetic follow-up program in the regime of routine detections with the advanced GW instruments expected within this decade. In that regime multi-wavelength observations will play a significant role in completing the astrophysical identification of GW sources. We present the methods and results from this first combined analysis and discuss its implications in terms of sensitivity for the present and future instruments.

Aims. A transient astrophysical event observed in both gravitational wave (GW) and electromagnetic (EM) channels would yield rich scientific rewards. A first program initiating EM follow-ups to possible transient GW events has been developed and exercised by the LIGO and Virgo community in association with several partners. In this paper, we describe and evaluate the methods used to promptly identify and localize GW event candidates and to request images of targeted sky locations. Methods. During two observing periods (Dec 17 2009 to Jan 8 2010 and Sep 2 to Oct 20 2010), a low-latency analysis pipeline was used to identify GW event candidates and to reconstruct maps of possible sky locations. A catalog of nearby galaxies and Milky Way globular clusters was used to select the most promising sky positions to be imaged, and this directional information was delivered to EM observatories with time lags of about thirty minutes. A Monte Carlo simulation has been used to evaluate the low-latency GW pipeline's ability to reconstruct source positions correctly. Results. For signals near the detection threshold, our low-latency algorithms often localized simulated GW burst signals to tens of square degrees, while neutron star/neutron star inspirals and neutron star/black hole inspirals were localized to a few hundred square degrees. Localization precision improves for moderately stronger signals. The correct sky location of signals well above threshold and originating from nearby galaxies may be observed with ~50% or better probability with a few pointings of wide-field telescopes.

We present the results of continued monitoring of the quiescent neutron star low-mass X-ray binary XTE J1701-462 with Chandra and Swift. A new Chandra observation from 2010 October extends our tracking of the neutron star surface temperature from ~800 days to ~1160 days since the end of an exceptionally luminous 19 month outburst. This observation indicates that the neutron star crust may still be slowly cooling toward thermal equilibrium with the core; another observation further into quiescence is needed to verify this. The shape of the overall cooling curve is consistent with that of a broken power law, although an exponential decay to a constant level cannot be excluded with the present data. To investigate possible low-level activity, we conducted a monitoring campaign of XTE J1701-462 with Swift during 2010 April-October. Short-term flares - presumably arising from episodic low-level accretion - were observed up to a luminosity of ~1e35 erg/s, ~20 times higher than the normal quiescent level. We conclude that flares of this magnitude are not likely to have significantly affected the equilibrium temperature of the neutron star and are probably not able to have a measurable impact on the cooling curve. However, it is possible that brighter and longer periods of low-level activity have had an appreciable effect on the equilibrium temperature.

A number of studies have revealed variability from neutron star low-mass X-ray binaries during quiescence. Such variability is not well characterised, or understood, but may be a common property that has been missed due to lack of multiple observations. One such source where variability has been observed is Aql X-1. Here, we analyse 14 Chandra and XMM-Newton observations of Aql X-1 in quiescence, covering a period of approximately 2 years. There is clear variability between the epochs, with the most striking feature being a flare-like increase in the flux by a factor of 5. Spectral fitting is inconclusive as to whether the power-law and/or thermal component is variable. We suggest that the variability and flare-like behaviour during quiescence is due to accretion at low rates which might reach the neutron star surface.

We report on an analysis of RXTE data of the transient neutron star low-mass X-ray binary (NS-LMXB) XTE J1701-462, obtained during its 2006-2007 outburst. The X-ray properties of the source changed between those of various types of NS-LMXB subclasses. At high luminosities the source switched between two types of Z source behavior and at low luminosities we observed a transition from Z source to atoll source behavior. These transitions between subclasses primarily manifest themselves as changes in the shapes of the tracks in X-ray color-color and hardness-intensity diagrams, but they are accompanied by changes in the kHz quasi-periodic oscillations, broad-band variability, burst behavior, and/or X-ray spectra. We find that the low-energy X-ray flux is a good parameter to track the gradual evolution of the tracks in color-color and hardness-intensity diagrams, allowing us to resolve the evolution of the source in greater detail than before and relate the observed properties to other NS-LMXBs. We further find that during the transition from Z to atoll, characteristic behavior known as the atoll upper banana can equivalently be described as the final stage of a weakening Z source flaring branch, thereby blurring the line between the two subclasses. Our findings strongly suggest that the wide variety in behavior observed in NS-LXMBs with different luminosities can be linked through changes in a single variable parameter, namely the mass accretion rate, without the need for additional differences in the neutron star parameters or viewing angle. We briefly discuss the implications of our findings for the spectral changes observed in NS LMXBs and suggest that, contrary to what is often assumed, the position along the color-color tracks of Z sources is not determined by the instantaneous mass accretion rate.

We present Rossi X-Ray Timing Explorer and Swift observations made during the final three weeks of the 2006-2007 outburst of the super-Eddington neutron star (NS) transient XTE J1701-462, as well as Chandra and XMM-Newton observations covering the first ~800 days of the subsequent quiescent phase. The source transitioned quickly from active accretion to quiescence, with the luminosity dropping by over 3 orders of magnitude in ~13 days. The spectra obtained during quiescence exhibit both a thermal component, presumed to originate in emission from the NS surface, and a non-thermal component of uncertain origin, which has shown large and irregular variability. We interpret the observed decay of the inferred effective surface temperature of the NS in quiescence as the cooling of the NS crust after having been heated and brought out of thermal equilibrium with the core during the outburst. The interpretation of the data is complicated by an apparent temporary increase in temperature ~220 days into quiescence, possibly due to an additional spurt of accretion. We derive an exponential decay timescale of ~120 (+30/-20) days for the inferred temperature (excluding observations affected by the temporary increase). This short timescale indicates a highly conductive NS crust. Further observations are needed to confirm whether the crust is still slowly cooling or has already reached thermal equilibrium with the core at a surface temperature of ~125 eV. The latter would imply a high equilibrium bolometric thermal luminosity of ~5x10^{33} erg/s for an assumed distance of 8.8 kpc.

We analyze data from five Chandra observations of the spiral galaxy NGC 6946 and from three Chandra observations of the irregular/spiral interacting galaxy pair NGC 4485/4490, with an emphasis on investigating the long-term variability exhibited by the source populations. We detect 90 point sources coincident with NGC 6946 down to luminosities of a few times 10^36 erg/s, and 38 sources coincident with NGC 4485/90 down to a luminosity of ~1x10^37 erg/s. Twenty-five (15) sources in NGC 6946 (NGC 4485/90) exhibit long-term (i.e., weeks to years) variability in luminosity; 11 (4) are transient candidates. The single ultraluminous X-ray source (ULX) in NGC 6946 and all but one of the eight ULXs in NGC 4485/90 exhibit long-term flux variability. Two of the ULXs in NGC 4485/90 have not been identified before as ultraluminous sources. The widespread variability in both systems is indicative of the populations being dominated by X-ray binaries, and this is supported by the X-ray colors of the sources. The distribution of colors among the sources indicates a large fraction of high-mass X-ray binaries in both systems. The shapes of the X-ray luminosity functions of the galaxies do not change significantly between observations and can be described by power laws with cumulative slopes ~0.6-0.7 (NGC 6946) and ~0.4 (NGC 4485/90).