IGR J17591$-$2342 is a new accreting millisecond X-ray pulsar (AMXP) that was recently discovered in outburst in 2018. Early observations revealed that the source's radio emission is brighter than that of any other known neutron star low-mass X-ray binary (NS-LMXB) at comparable X-ray luminosity, and assuming its likely $\gtrsim 6$ kpc distance. It is comparably radio bright to black hole LMXBs at similar X-ray luminosities. In this work, we present the results of our extensive radio and X-ray monitoring campaign of the 2018 outburst of IGR J17591$-$2342. In total we collected 10 quasi-simultaneous radio (VLA, ATCA) and X-ray (Swift-XRT) observations, which make IGR J17591$-$2342 one of the best-sampled NS-LMXBs. We use these to fit a power-law correlation index $\beta = 0.37^{+0.42}_{-0.40}$ between observed radio and X-ray luminosities ( $L_\mathrm{R}\propto L_\mathrm{X}^{\beta}$). However, our monitoring revealed a large scatter in IGR J17591$-$2342's radio luminosity (at a similar X-ray luminosity, $L_\mathrm{X} \sim 10^{36}$ erg s$^{-1}$, and spectral state), with $L_\mathrm{R} \sim 4 \times 10^{29}$ erg s$^{-1}$ during the first three reported observations, and up to a factor of 4 lower $L_\mathrm{R}$ during later radio observations. Nonetheless, the average radio luminosity of IGR J17591$-$2342 is still one of the highest among NS-LMXBs, and we discuss possible reasons for the wide range of radio luminosities observed in such systems during outburst. We found no evidence for radio pulsations from IGR J17591$-$2342 in our Green Bank Telescope observations performed shortly after the source returned to quiescence. Nonetheless, we cannot rule out that IGR J17591$-$2342 becomes a radio millisecond pulsar during quiescence.

Aql X-1 is one of the best-studied neutron star low-mass X-ray binaries. It was previously targeted using quasi-simultaneous radio and X-ray observations during at least 7 different accretion outbursts. Such observations allow us to probe the interplay between accretion inflow (X-ray) and jet outflow (radio). Thus far, these combined observations have only covered one order of magnitude in radio and X-ray luminosity range; this means that any potential radio - X-ray luminosity correlation, $L_R \propto L_X^{\beta}$, is not well constrained ($\beta \approx$ 0.4-0.9, based on various studies) or understood. Here we present quasi-simultaneous Very Large Array and Swift-XRT observations of Aql X-1's 2016 outburst, with which we probe one order of magnitude fainter in radio and X-ray luminosity compared to previous studies ($6 \times 10^{34} < L_X < 3 \times 10^{35}$ erg s$^{-1}$, i.e., the intermediate to low-luminosity regime between outburst peak and quiescence). The resulting radio non-detections indicate that Aql X-1's radio emission decays more rapidly at low X-ray luminosities than previously assumed - at least during the 2016 outburst. Assuming similar behaviour between outbursts, and combining all available data, this can be modelled as a steep $\beta=1.17^{+0.30}_{-0.21}$ power-law index or as a sharp radio cut-off at $L_X \lesssim 5 \times 10^{35}$ erg s$^{-1}$ (given our deep radio upper limits at X-ray luminosities below this value). We discuss these results in the context of other similar studies.

IGR J17591$-$2342 is a 527-Hz accreting millisecond X-ray pulsar that was discovered in outburst in 2018 August. In this paper, we present quasi-simultaneous radio and X-ray monitoring of this source during the early part of the outburst. IGR J17591$-$2342 is highly absorbed in X-rays, with an equivalent hydrogen absorption along the line of sight, $N_{\rm H}$, of $\approx$4.4$\times$10$^{22}$\,cm$^{-2}$, where the Galactic column density is expected to be $\approx$1 -- 2$\times$10$^{22}$\,cm$^{-2}$. The high absorption suggests that the source is either relatively distant ($>$6\,kpc), or the X-ray emission is strongly absorbed by material local to the system. Radio emission detected by the Australia Telescope Compact Array shows that, for a given X-ray luminosity and for distances greater than 3\,kpc, this source was exceptionally radio loud when compared to other accreting neutron stars in outburst ($L_{\rm X} > 10^{33}$\,erg\,s$^{-1}$). For most reasonable distances, IGR J17591$-$2342 appeared as radio luminous as actively accreting, stellar-mass black hole X-ray binaries.

Einstein's theory of gravity, general relativity, has passed stringent tests in laboratories, elsewhere in the Solar Sytem, and in pulsar binaries. Nevertheless it is known to be incompatible with quantum mechanics and must differ from the true behaviour of matter in strong fields and at small spatial scales. A key aspect of general relativity to test is the strong equivalence principle (SEP), which states that all freely falling objects, regardless of how strong their gravity, experience the same acceleration in the same gravitational field. Essentially all alternatives to general relativity violate this principle at some level. Previous direct tests of the SEP are limited by the weak gravity of the bodies in the Earth-Moon-Sun system or by the weak gravitational pull of the Galaxy on pulsar-white dwarf binaries. PSR~J0337+1715 is a hierarchical stellar triple system, where the inner binary consists of a millisecond radio pulsar in a $1.6$-day orbit with a white dwarf. This inner binary is in a $327$-day orbit with another white dwarf. In this system, the pulsar and the inner companion fall toward the outer companion with an acceleration about $10^8$ times greater than that produced by falling in the Galactic potential, and the pulsar's gravitational binding energy is roughly $10\%$ of its mass. Here we report that in spite of the pulsar's strong gravity, the accelerations experienced by it and the inner white dwarf differ by a fraction of no more than $2.6\times 10^{-6}$ ($95\%$ confidence level). We can roughly compare this to other SEP tests by using the strong-field Nordtvedt parameter $\hat\eta_N$. Our limit on $\hat\eta_N$ is a factor of ten smaller than that obtained from (weak-field) Solar-System SEP tests and a factor of almost a thousand smaller than that obtained from other strong-field SEP tests.

We report the discovery of PSR J0952$-$0607, a 707-Hz binary millisecond pulsar which is now the fastest-spinning neutron star known in the Galactic field (i.e., outside of a globular cluster). PSR J0952$-$0607 was found using LOFAR at a central observing frequency of 135 MHz, well below the 300 MHz to 3 GHz frequencies typically used in pulsar searches. The discovery is part of an ongoing LOFAR survey targeting unassociated Fermi Large Area Telescope $\gamma$-ray sources. PSR J0952$-$0607 is in a 6.42-hr orbit around a very low-mass companion ($M_\mathrm{c}\gtrsim0.02$ M$_\odot$) and we identify a strongly variable optical source, modulated at the orbital period of the pulsar, as the binary companion. The light curve of the companion varies by 1.6 mag from $r^\prime=22.2$ at maximum to $r^\prime>23.8$, indicating that it is irradiated by the pulsar wind. Swift observations place a 3-$\sigma$ upper limit on the $0.3-10$ keV X-ray luminosity of $L_X < 1.1 \times 10^{31}$ erg s$^{-1}$ (using the 0.97 kpc distance inferred from the dispersion measure). Though no eclipses of the radio pulsar are observed, the properties of the system classify it as a black widow binary. The radio pulsed spectrum of PSR J0952$-$0607, as determined through flux density measurements at 150 and 350 MHz, is extremely steep with $\alpha\sim-3$ (where $S \propto \nu^{\alpha}$). We discuss the growing evidence that the fastest-spinning radio pulsars have exceptionally steep radio spectra, as well as the prospects for finding more sources like PSR J0952$-$0607.

We present quasi-simultaneous radio (VLA) and X-ray ($Swift$) observations of the neutron star low-mass X-ray binary (NS-LMXB) 1RXS J180408.9$-$342058 (J1804) during its 2015 outburst. We found that the radio jet of J1804 was bright ($232 \pm 4 \mu$Jy at $10$ GHz) during the initial hard X-ray state, before being quenched by more than an order of magnitude during the soft X-ray state ($19 \pm 4 \mu$Jy). The source then was undetected in radio (< $13 \mu$Jy) as it faded to quiescence. In NS-LMXBs, possible jet quenching has been observed in only three sources and the J1804 jet quenching we show here is the deepest and clearest example to date. Radio observations when the source was fading towards quiescence ($L_X = 10^{34-35}$ erg s$^{-1}$) show that J1804 must follow a steep track in the radio/X-ray luminosity plane with $\beta > 0.7$ (where $L_R \propto L_X^{\beta}$). Few other sources have been studied in this faint regime, but a steep track is inconsistent with the suggested behaviour for the recently identified class of transitional millisecond pulsars. J1804 also shows fainter radio emission at $L_X < 10^{35}$ erg s$^{-1}$ than what is typically observed for accreting millisecond pulsars. This suggests that J1804 is likely not an accreting X-ray or transitional millisecond pulsar.

We report on unusually very hard spectral states in three confirmed neutron-star low-mass X-ray binaries (1RXS J180408.9-342058, EXO 1745-248, and IGR J18245-2452) at a luminosity between ~ 10^{36-37} erg s^{-1}. When fitting the Swift X-ray spectra (0.5 - 10 keV) in those states with an absorbed power-law model, we found photon indices of \Gamma ~ 1, significantly lower than the \Gamma = 1.5 - 2.0 typically seen when such systems are in their so called hard state. For individual sources very hard spectra were already previously identified but here we show for the first time that likely our sources were in a distinct spectral state (i.e., different from the hard state) when they exhibited such very hard spectra. It is unclear how such very hard spectra can be formed; if the emission mechanism is similar to that operating in their hard states (i.e., up-scattering of soft photons due to hot electrons) then the electrons should have higher temperatures or a higher optical depth in the very hard state compared to those observed in the hard state. By using our obtained \Gamma as a tracer for the spectral evolution with luminosity, we have compared our results with those obtained by Wijnands et al. (2015). We confirm their general results in that also our sample of sources follow the same track as the other neutron star systems, although we do not find that the accreting millisecond pulsars are systematically harder than the non-pulsating systems.

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.

1RXS J180408.9-342058 is a transient neutron star low-mass X-ray binary that exhibited a bright accretion outburst in 2015. We present Nustar, Swift, and Chandra observations obtained around the peak of this outburst. The source was in a soft X-ray spectral state and displayed an X-ray luminosity of Lx~(2-3)E37 (D/5.8 kpc)^2 erg cm-2 (0.5-10 keV). The Nustar data reveal a broad Fe-K emission line that we model as relativistically broadened reflection to constrain the accretion geometry. We found that the accretion disk is viewed at an inclination of i~27-35 degrees and extended close to the neutron star, down to Rin~5-7.5 gravitational radii (~11-17 km). This inner disk radius suggests that the neutron star magnetic field strength is B<2E8 G. We find a narrow absorption line in the Chandra/HEG data at an energy of ~7.64 keV with a significance of ~4.8 sigma. This feature could correspond to blue-shifted Fe xxvi and arise from an accretion disk wind, which would imply an outflow velocity of v~0.086c (~25800 km s-1). However, this would be extreme for an X-ray binary and it is unclear if a disk wind should be visible at the low inclination angle that we infer from our reflection analysis. Finally, we discuss how the X-ray and optical properties of 1RXS J180408.9-342058 are consistent with a relatively small (Porb<3 hr) binary orbit.