Search

JWST, despite not being designed to observe astrophysical phenomena that vary on rapid time scales, can be an unparalleled tool for such studies. If timing systematics can be controlled, JWST will be able to open up the sub-second infrared timescale regime. Rapid time-domain studies, such as lag measurements in accreting compact objects and Solar System stellar occultations, require both precise inter-frame timing and knowing when a time series begins to an absolute accuracy significantly below 1s. In this work we present two long-duration observations of the deeply eclipsing double white dwarf system ZTF J153932.16+502738.8, which we use as a natural timing calibrator to measure the absolute timing accuracy of JWST's clock. From our two epochs, we measure an average clock accuracy of $0.12\pm0.06$s, implying that JWST can be used for sub-second time-resolution studies down to the $\sim100$ms level, a factor $\sim5$ improvement upon the pre-launch clock accuracy requirement. We also find an asymmetric eclipse profile in the F322W2 band, which we suggest has a physical origin., Comment: 16 pages, 13 figures, accepted for publication in AJ

QS Vir is a low-accretion rate cataclysmic variable (CV), or pre-CV, as the M dwarf companion is just filling its Roche lobe. We recently identified radio emission from QS Vir in the Very Large Array Sky Survey, at a flux of ~1 mJy. The origin of radio emission from CVs is not fully understood, with evidence for synchrotron emission from jets and other coherent plasma emission processes, such as electron cyclotron maser emission (ECME) or plasma radiation. Our aim is to constrain the radio emission mechanism for QS Vir, through spectroscopic, polarisation, and time variability measurements, all while checking for correlated X-ray variations. We took 3 epochs of new observations with the VLA in S, C, and X bands, with full Stokes polarisation information, complemented by near-simultaneous Swift/XRT X-ray data. Radio spectra are extracted to search for emission features characteristic of coherent plasma emission processes (e.g. high circular polarisation and narrow-band emission). We fit the X-ray spectra with absorbed power-laws, finding no strong X-ray variability. QS Vir showed a nearly flat radio spectrum, with fluxes of 0.4-0.6 mJy in all bands. Swift/XRT showed L_X ~ 5x10^29 erg/s in all observations. We identified strong, variable circular polarisation, ranging from 33+/-3% in S band in the last observation, to <11% in the middle observation in all bands. Linear polarisation was not detected, with upper limits of at most 15%. Intriguingly, the S-band spectra show circularly polarised spectral bumps (width ~0.5 GHz) that rise and decay within <5 minutes. We suggest that the radio emission from QS Vir consists of two components: a relatively constant, low-polarisation flat-spectrum component and a band-limited, rapidly variable, and strongly circularly polarised component. This latter coherent component may be associated with ECME or plasma radiation.

We present mid-infrared (MIR) spectral-timing measurements of the prototypical Galactic microquasar GRS 1915+105. The source was observed with the Mid-Infrared Instrument (MIRI) onboard JWST in June 2023 at a MIR luminosity L(MIR)~10^{36} erg/s exceeding past IR levels by about a factor of 10. By contrast, the X-ray flux is much fainter than the historical average, in the source's now-persistent 'obscured' state. The MIRI low-resolution spectrum shows a plethora of emission lines, the strongest of which are consistent with recombination in the hydrogen Pfund (Pf) series and higher. Low amplitude (~1%) but highly significant peak-to-peak photometric variability is found on timescales of ~1,000 s. The brightest Pf(6-5) emission line lags the continuum. Though difficult to constrain accurately, this lag is commensurate with light-travel timescales across the outer accretion disc or with expected recombination timescales inferred from emission line diagnostics. Using the emission line as a bolometric indicator suggests a moderate (~5-30% Eddington) intrinsic accretion rate. Multiwavelength monitoring shows that JWST caught the source close in-time to unprecedentedly bright MIR and radio long-term flaring. Assuming a thermal bremsstrahlung origin for the MIRI continuum suggests an unsustainably high mass-loss rate during this time unless the wind remains bound, though other possible origins cannot be ruled out. PAH features previously detected with Spitzer are now less clear in the MIRI data, arguing for possible destruction of dust in the interim. These results provide a preview of new parameter space for exploring MIR spectral-timing in XRBs and other variable cosmic sources on rapid timescales., Comment: Dedicated to the memory of our colleague, Tomaso Belloni. Submitted 2024 June 21; Comments welcome

We present the results of a high angular resolution (1.1") and sensitivity (maximum of ~0.1 mJy) radio survey at 1-2 GHz in the Galactic Bulge. This complements the X-ray Chandra Galactic Bulge Survey, and investigates the full radio source population in this dense Galactic region. Radio counterparts to sources at other wavelengths can aid in classification, as there are relatively few types of objects that are reasonably detectable in radio at kiloparsec distances, and even fewer that are detected in both X-rays and radio. This survey covers about 3 square degrees of the Galactic Bulge Survey area (spanning the Galactic coordinate range of -3deg < l < +3deg and +1.6deg < b < +2.1deg as a first look into this region of the Galaxy with this combination of frequency, resolution, and sensitivity. Spectral indices within the observed band of 1-2 GHz were calculated for each source to assist in determining its emission mechanism. We find 1617 unique sources in the survey, 25 of which are radio counterparts to X-ray sources, and about 100 of which are steep-spectrum (alpha <~ -1.4) point sources that are viable pulsar candidates. Four radio sources are of particular interest: a compact binary; an infrared transient with an inverted radio spectrum; a potential transitional millisecond pulsar candidate; and a very steep spectrum radio source with an X-ray and bright infrared counterpart. We discuss other notable sources, including possible radio transients, potential new planetary nebulae, and active galactic nuclei., Comment: 20 pages, 13 figures. Accepted for publication in MNRAS

The Snake is a remarkable Galactic center radio filament with a morphology characterized by two kinks along its $\sim 20'$ extent. The major and minor kinks are located where the filament is most distorted from a linear magnetized structure running perpendicular to the Galactic plane. We present {\em Chandra}, VLA, and MeerKAT data and report the detection of an X-ray and radio source at the location of the major kink. High-resolution radio images of the major kink reveal a compact source with a steep spectrum with spectral index alpha ~ -2.7 surrounded by extended emission. The radio luminosity and steep spectrum of the compact source are consistent with a pulsar. We also show flattening of the spectrum and enhanced synchrotron emissivity away from the position of the major kink along the Snake, which suggests injection of relativistic particles along the Snake. We argue that the major kink is created by a fast-moving (~500-1000 km/s), object punching into the Snake, distorting its magnetic structure, and producing X-ray emission. X-ray emission pinpoints an active acceleration site where the interaction is taking place. A secondary kink is argued to be induced by the impact of the high-velocity object producing the major kink., Comment: 12 pages, 8 figures, MNRAS (in press)

The conditions under which accreting neutron stars launch radio-emitting jets and/or outflows are still poorly understood. The ultracompact X-ray binary X1850--087, located in the globular cluster NGC 6712, is a persistent atoll-type X-ray source that has previously shown unusual radio continuum variability. Here we present the results of a pilot radio monitoring program of X1850--087 undertaken with the Karl G. Jansky Very Large Array, with simultaneous or quasi-simultaneous Swift/XRT data obtained at each epoch. The binary is clearly detected in the radio in two of the six new epochs. When combined with previous data, these results suggest that X1850--087 shows radio emission at a slightly elevated hard state X-ray luminosity of L_X > 2x10^36 erg/s, but no radio emission in its baseline hard state L_X ~10^36 erg/s. No clear X-ray spectral changes are associated with this factor of > 10 radio variability. At all detected epochs X1850--087 has a flat-to-inverted radio spectral index, more consistent with the partially absorbed optically thick synchrotron of a compact jet rather than the evolving optically thick to thin emission associated with transient expanding synchrotron-emitting ejecta. If the radio emission in X1850--087 is indeed due to a compact jet, then it is plausibly being launched and quenched in the hard state on timescales as short as a few days. Future radio monitoring of X1850--087 could help elucidate the conditions under which compact jets are produced around hard state accreting neutron stars., Comment: 9 pages, 2 figures, 2 tables, accepted for publication in ApJ

We use the Very Long Baseline Array to conduct high precision astrometry of a sample of 33 compact, flat spectrum, variable radio sources in the direction of the Galactic plane (Becker et al. 2010). Although Becker et al. (2010) ruled out a few potential scenarios for the origin of the radio emission, the study could not rule out that these sources were black hole X-ray binaries (BHXBs). Most known BHXBs are first detected by X-ray or optical emission when they go into an outburst, leaving the larger quiescent BHXB population undiscovered. In this paper, we attempt to identify any Galactic sources amongst the Becker et al. (2010) sample by measuring their proper motions as a first step to finding quiescent BHXB candidates. Amongst the 33 targets, we could measure the proper motion of six sources. We find that G32.7193$-$0.6477 is a Galactic source and are able to constrain the parallax of this source with a 3$\sigma$ significance. We found three strong Galactic candidates, G32.5898$-$0.4468, G29.1075$-$0.1546, and G31.1494$-$0.1727, based purely on their proper motions, and suggest that G29.1075$-$0.1546, is also likely Galactic. We detected two resolved targets for multiple epochs (G30.1038+0.3984 and G29.7161$-$0.3178). We find six targets are only detected in one epoch and have an extended structure. We cross-match our VLBA detections with the currently available optical, infrared and X-ray surveys, and did not find any potential matches. We did not detect 19 targets in any VLBA epochs and suggest that this could be due to limited $uv$-coverage, drastic radio variability or faint, extended nature of the sources., Comment: 18 pages, 7 figures, Accepted in MNRAS

We observed the Galactic black hole X-ray binary V4641 Sgr with the high resolution transmission gratings on Chandra during the source's 2020 outburst. Over two epochs of Chandra gratings observations, we see numerous highly ionized metal lines, superimposed on a hot, disc-dominated X-ray continuum. The measured inner disc temperatures and luminosities imply an unfeasibly small inner disc radius, such that we suggest that the central engine of V4641 Sgr is obscured, and we are viewing scattered X-rays. We find that the emission lines in the Chandra spectra cannot be constrained by a single photoionized model, instead finding that two separate photoionized model components are required, one to reproduce the iron lines and a second for the other metals. We compare the observed X-ray spectra of V4641 Sgr to optical studies during previous outbursts of the source, suggesting that the lines originate in an accretion disc wind, potentially with a spherical geometry., Comment: 14 pages, 11 figures, accepted for publication in MNRAS

We present a detailed time-resolved photometric study of the ultra-compact X-ray binary candidate 4U 1812-12. The multicolor light curves obtained with HiPERCAM on the 10.4-m Gran Telescopio Canarias show an aprox 114 min modulation similar to a superhump. Under this interpretation, this period should lie very close to the orbital period of the system. Contrary to what its other observational properties suggest (namely, persistent dim luminosity, low optical-to-X-ray flux ratio and lack of hydrogen features in the optical spectrum), this implies that 4U1812-12 is most likely not an ultra-compact X-ray binary, which are usually defined as systems with orbital periods lower than 80 min. We discuss the nature of the system, showing that a scenario in which 4U 1812-12 is the progenitor of an ultra-compact X-ray binary may reconcile all the observables., Comment: Accepted for publication in ApJ Letters

The radiative counterpart of the supermassive black hole at the Galactic Centre, Sagittarius A*, displays flaring emission in the X-ray band atop a steady, quiescent level. Flares are also observed in the near-infrared band. The physical process producing the flares is not fully understood and it is unclear if the flaring rate varies, although some recent works suggest it has reached unprecedented variability in recent years. Using over a decade of regular X-ray monitoring of Neil Gehrels Swift Observatory, we studied the variations in count rate of Sgr A* on time scales of years. We decomposed the X-ray emission into quiescent and flaring emission, modelled as a constant and power law process, respectively. We found that the complete, multi-year dataset cannot be described by a stationary distribution of flare fluxes, while individual years follow this model better. In three of the ten studied years, the data is consistent with a purely Poissonian quiescent distribution, while for five years only an upper limit of the flare flux distribution parameter could be determined. We find that these possible changes cannot be explained fully by the different number of observations per year. Combined, these results are instead consistent with a changing flaring rate of Sgr A*, appearing more active between 2006-2007 and 2017-2019, than between 2008-2012. Finally, we discuss this result in the context of flare models and the passing of gaseous objects, and discuss the extra statistical steps taken, for instance to deal with the background in the Swift observations., Comment: 14 pages, 10 figures, Accepted for publication in MNRAS

Accreting neutron stars in low-mass X-ray binaries show outflows -- and sometimes jets -- in the general manner of accreting black holes. However, the quantitative link between the accretion flow (traced by X-rays) and outflows and/or jets (traced by radio emission) is much less well-understood for neutron stars than for black holes, other than the general observation that neutron stars are fainter in the radio at a given X-ray luminosity. We use data from the deep MAVERIC radio continuum survey of Galactic globular clusters for a systematic radio and X-ray study of six luminous (L_X > 10^34 erg/s) persistent neutron star X-ray binaries in our survey, as well as two other transient systems also captured by our data. We find that these neutron star X-ray binaries show an even larger range in radio luminosity than previously observed. In particular, in quiescence at L_X ~ 3x10^34 erg/s, the confirmed neutron star binary GRS 1747--312 in Terzan 6 sits near the upper envelope of the black hole radio/X-ray correlation, and the persistently accreting neutron star systems AC 211 (in M15) and X1850--087 (in NGC 6712) show unusual radio variability and luminous radio emission. We interpret AC 211 as an obscured "Z source" that is accreting at close to the Eddington limit, while the properties of X1850--087 are difficult to explain, and motivate future coordinated radio and X-ray observations. Overall, our results show that neutron stars do not follow a single relation between inflow and outflow, and confirm that their accretion dynamics are more complex than for black holes., Comment: 27 pages, 9 figures, 4 tables, accepted for publication in ApJ

ASASSN-21au is an ultracompact accreting white dwarf binary (AM CVn) with a period of approx. 58 min. Using multiwavelength observations of the system, we discovered a dichotomy in the behavior of outbursts in AM CVns. The binary showed an initial brightness increase which lasted for at least 82 days, followed by an additional increase which lasted 2 weeks. Afterwards ASASSN-21au went into superoutburst with a total duration of 19 days, showing an amplitude with respect to quiescence of approx. 7.5 mags in g, with a precursor and an echo outburst. A correlation between X-rays, UV and optical was identified for the first time in an AM CVn during this stage. The color evolution of ASASSN-21au indicates that during the superoutburst the dominant component was the accretion disk. The short duration, large amplitude and color evolution of the superoutburst agree with expectations from the disk instability model. These characteristics are opposite to the ones observed in SDSS J080710+485259 and SDSS J113732+405458, which have periods of approx. 53 min and 60 min, respectively. The initially slow brightness increase in the light curve of ASASSN-21au and the behavior after the superoutburst favors a scenario in which changes in the mass-transfer rate led to disk instabilities, while the outburst mechanism of SDSS J080710+485259 and SDSS J113732+405458 has been attributed to enhanced mass-transfer alone. Further observations are needed to understand the origin of this dichotomy., Comment: Published in ApJ in February 2022

We report new radio observations of a sample of thirty-six neutron star (NS) X-ray binaries, more than doubling the sample in the literature observed at current-day sensitivities. These sources include thirteen weakly-magnetised ($B<10^{10}$ G) and twenty-three strongly-magnetised ($B\geq10^{10}$ G) NSs. Sixteen of the latter category reside in high-mass X-ray binaries, of which only two systems were radio-detected previously. We detect four weakly and nine strongly-magnetised NSs; the latter are systematically radio fainter than the former and do not exceed $L_R \approx 3\times10^{28}$ erg/s. In turn, we confirm the earlier finding that the weakly-magnetized NSs are typically radio fainter than accreting stellar-mass black holes. While an unambiguous identification of the origin of radio emission in high-mass X-ray binaries is challenging, we find that in all but two detected sources (Vela X-1 and 4U 1700-37) the radio emission appears more likely attributable to a jet than the donor star wind. The strongly-magnetised NS sample does not reveal a global correlation between X-ray and radio luminosity, which may be a result of sensitivity limits. Furthermore, we discuss the effect of NS spin and magnetic field on radio luminosity and jet power in our sample. No current model can account for all observed properties, necessitating the development and refinement of NS jet models to include magnetic field strengths up to $10^{13}$ G. Finally, we discuss jet quenching in soft states of NS low-mass X-ray binaries, the radio non-detections of all observed very-faint X-ray binaries in our sample, and future radio campaigns of accreting NSs., Comment: Accepted for publication in MNRAS; 31 pages, 8 figures, 8 tables (including online supplementary materials)

Very faint X-ray transients (VFXTs) are X-ray transients with peak X-ray luminosities ($L_X$) of $L_X\lesssim10^{36}$ erg/s, which are not well-understood. We carried out a survey of 16 square degrees of the Galactic Bulge with the Swift Observatory, using short (60 s) exposures, and returning every 2 weeks for 19 epochs in 2017-18 (with a gap from November 2017 to February 2018, when the Bulge was in sun-constraint). Our main goal was to detect and study VFXT behaviour in the Galactic Bulge across various classes of X-ray sources. In this work, we explain the observing strategy of the survey, compare our results with the expected number of source detections per class, and discuss the constraints from our survey on the Galactic VFXT population. We detected 91 X-ray sources, 25 of which have clearly varied by a factor of at least 10. 45 of these X-ray sources have known counterparts: 17 chromospherically active stars, 12 X-ray binaries, 5 cataclysmic variables (and 4 candidates), 3 symbiotic systems, 2 radio pulsars, 1 AGN, and a young star cluster. The other 46 are of previously undetermined nature. We utilize X-ray hardness ratios, searches for optical/infrared counterparts in published catalogs, and flux ratios from quiescence to outburst to constrain the nature of the unknown sources. Of these 46, 7 are newly discovered hard transients, which are likely VFXT X-ray binaries. Furthermore, we find strong new evidence for a symbiotic nature of 4 sources in our full sample, and new evidence for accretion power in 6 X-ray sources with optical counterparts. Our findings indicate that a large subset of VXFTs is likely made up of symbiotic systems., Comment: 23 pages, 9 figures, 7 tables, MNRAS in press. The machine-readable catalog and the complete figure-set of light curves included in the arxiv source (accessible via "other formats")

The hard X-ray spectrum of magnetic cataclysmic variables can be modelled to provide a measurement of white dwarf mass. This method is complementary to radial velocity measurements, which depend on the (typically rather uncertain) binary inclination. Here we present results from a Legacy Survey of 19 magnetic cataclysmic variables with NuSTAR. We fit accretion column models to their 20-78 keV spectra and derive the white dwarf masses, finding a weighted average $\bar{M}_{\rm WD}=0.77\pm0.02$ $M_{\odot}$, with a standard deviation $\sigma=0.10$ $M_{\odot}$, when we include the masses derived from previous NuSTAR observations of seven additional magnetic cataclysmic variables. We find that the mass distribution of accreting magnetic white dwarfs is consistent with that of white dwarfs in non-magnetic cataclysmic variables. Both peak at a higher mass than the distributions of isolated white dwarfs and post-common-envelope binaries. We speculate as to why this might be the case, proposing that consequential angular momentum losses may play a role in accreting magnetic white dwarfs and/or that our knowledge of how the white dwarf mass changes over accretion-nova cycles may also be incomplete., Comment: 14 pages. 5 figures in main paper, 1 figure in appendix. Accepted for publication in MNRAS

Swift J1858.6-0814 is a transient neutron star X-ray binary discovered in October 2018. Multi-wavelength follow-up observations across the electromagnetic spectrum revealed many interesting properties, such as erratic flaring on minute timescales and evidence for wind outflows at both X-ray and optical wavelengths, strong and variable local absorption, and an anomalously hard X-ray spectrum. Here, we report on a detailed radio observing campaign consisting of one observation at 5.5/9 GHz with the Australia Telescope Compact Array, and nine observations at 4.5/7.5 GHz with the Karl G. Jansky Very Large Array. A radio counterpart with a flat to inverted radio spectrum is detected in all observations, consistent with a compact jet being launched from the system. Swift J1858.6-0814 is highly variable at radio wavelengths in most observations, showing significant variability when imaged on 3-to-5-minute timescales and changing up to factors of 8 within 20 minutes. The periods of brightest radio emission are not associated with steep radio spectra, implying they do not originate from the launching of discrete ejecta. We find that the radio variability is similarly unlikely to have a geometric origin, be due to scintillation, or be causally related to the observed X-ray flaring. Instead, we find that it is consistent with being driven by variations in the accretion flow propagating down the compact jet. We compare the radio properties of Swift J1858.6-0814 with those of Eddington-limited X-ray binaries with similar X-ray and optical characteristics, but fail to find a match in radio variability, spectrum, and luminosity., Comment: 14 pages, 8 figures, 3 tables; accepted for publication in MNRAS

X-ray transients, such as accreting neutron stars, periodically undergo outbursts, thought to be caused by a thermal-viscous instability in the accretion disk. Usually outbursts of accreting neutron stars are identified when the accretion disk has undergone an instability, and the persistent X-ray flux has risen to a threshold detectable by all sky monitors on X-ray space observatories. Here we present the earliest known combined optical, UV, and X-ray monitoring observations of the outburst onset of an accreting neutron star low mass X-ray binary system. We observed a significant, continuing increase in the optical i'-band magnitude starting on July 25, 12 days before the first X-ray detection with Swift/XRT and NICER (August 6), during the onset of the 2019 outburst of SAX J1808.4-3658. We also observed a 4 day optical to X-ray rise delay, and a 2 day UV to X-ray delay, at the onset of the outburst. We present the multiwavelength observations that were obtained, discussing the theory of outbursts in X-ray transients, including the disk instability model, and the implications of the delay. This work is an important confirmation of the delay in optical to X-ray emission during the onset of outbursts in low mass X-ray binaries, which has only previously been measured with less sensitive all sky monitors. We find observational evidence that the outburst is triggered by ionisation of hydrogen in the disk., Comment: 12 pages, 5 figures, accepted by MNRAS

The dynamical production of low-mass X-ray binaries and brighter cataclysmic variables (CVs) in dense globular clusters is well-established. We investigate how the X-ray emissivity of fainter X-ray binaries (principally CVs and coronally active binaries) varies between different environments. We compile calculations (largely from the literature) of the X-ray emissivity of old stellar populations, including open and globular clusters and several galaxies. We investigate three literature claims of unusual X-ray sources in low-density stellar populations. We show that a suggested quiescent neutron star in the open cluster NGC 6819 is a foreground M dwarf. We show that the suggested diffuse X-ray emission from an old nova shell in the globular cluster NGC 6366 is actually a background galaxy cluster. And we show that a suggested population of quiescent X-ray binaries in the Sculptor Dwarf Galaxy is mostly (perhaps entirely) background galaxies. We find that above densities of $10^4$ M$_{\odot}$/pc$^3$, the X-ray emissivity of globular clusters increases, due to dynamical production of X-ray emitting systems. Below this density, globular clusters have lower X-ray emissivity than the other populations, and we do not see a strong dependence of X-ray emissivity due to density effects. We find significant correlations between X-ray emissivity and binary fraction, metallicity, and density. Sampling these fits via bootstrap techniques gives less significant correlations, but confirms the effect of metallicity on low-density populations, and that of density on the full globular cluster sample., Comment: 28 pages, 23 figures. MNRAS, in press

The nature of very faint X-ray transients (VFXTs) - transient X-ray sources that peak at luminosities $L_X\lesssim10^{36} {\rm erg s^{-1}}$ - is poorly understood. The faint and often short-lived outbursts make characterising VFXTs and their multi-wavelength counterparts difficult. In 2017 April we initiated the Swift Bulge Survey, a shallow X-ray survey of $\sim$16 square degrees around the Galactic centre with the Neil Gehrels Swift Observatory. The survey has been designed to detect new and known VFXTs, with follow-up programmes arranged to study their multi-wavelength counterparts. Here we detail the optical and near-infrared follow-up of four sources detected in the first year of the Swift Bulge Survey. The known neutron star binary IGR J17445-2747 has a K4III donor, indicating a potential symbiotic X-ray binary nature and the first such source to show X-ray bursts. We also find one nearby M-dwarf (1SXPS J174215.0-291453) and one system without a clear near-IR counterpart (Swift J175233.9-290952). Finally, 3XMM J174417.2-293944 has a subgiant donor, an 8.7 d orbital period, and a likely white dwarf accretor; we argue that this is the first detection of a white dwarf accreting from a gravitationally focused wind. A key finding of our follow-up campaign is that binaries containing (sub)giant stars may make a substantial contribution to the VFXT population., Comment: 18 pages, 11 figures, accepted for publication in MNRAS Replacement: Some of the source names were missing a digit from the declination, rectified this in latest version

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., Comment: 15 pages, 11 figures, 3 tables, accepted for publication in MNRAS

Accretion at super-Eddington rates is expected to be accompanied by strong outflows. Such outflows are observed in Galactic X-ray binaries and extragalactic Ultra-luminous X-ray sources (ULXs). However, due to their large source distances, ULX outflows are challenging to detect and study in detail. Galactic neutron stars accreting from a Be-star companion at super-Eddington rates show many similarities to ULX pulsars, and therefore offer an alternative approach to study outflows in this accretion regime. Here, we present Chandra high-resolution spectroscopy of such a super-Eddington accreting neutron star, Swift J0243.6+6124, to search for wind outflow signatures during the peak of its 2017/2018 giant outburst. We detect narrow emission features at rest from Ne, Mg, S, Si, and Fe. In addition, we detect a collection of absorption features which can be identified in two ways: either as all Fe transitions at rest (with a possible contribution from Mg), or a combination of three blue-shifted Ne and Mg lines at $\sim 0.22c$, while the remaining lines are at rest. The second scenario would imply an outflow with a velocity similar to those seen in ULXs, including the ULX pulsar NGC 300 ULX-1. This result would also imply that Swift J0243.6+6124 launches both a jet, detected in radio and reported previously, and an ultra-fast wind outflow simultaneously at super-Eddington accretion rates., Comment: 14 pages, 4 figures, 3 tables. Accepted for publication in MNRAS

We report optical and infrared observations of the X-ray source CXOGBS J174623.5-310550. This Galactic object was identified as a potential quiescent low-mass X-ray binary accreting from an M-type donor on the basis of optical spectroscopy and the broad Halpha emission line. The analysis of X-shooter spectroscopy covering 3 consecutive nights supports an M2/3-type spectral classification. Neither radial velocity variations nor rotational broadening is detected in the photospheric lines. No periodic variability is found in I- and r'-band light curves. We derive r' = 20.8, I = 19.2 and Ks = 16.6 for the optical and infrared counterparts with the M-type star contributing 90% to the I-band light. We estimate its distance to be 1.3-1.8 kpc. The lack of radial velocity variations implies that the M-type star is not the donor star in the X-ray binary. This could be an interloper or the outer body in a hierarchical triple. We constrain the accreting binary to be a < 2.2 hr orbital period eclipsing cataclysmic variable or a low-mass X-ray binary lying in the foreground of the Galactic Bulge., Comment: (9 pages, 5 figures, accepted for publication in MNRAS)

We report on the evolution of the X-ray emission of the accreting neutron star (NS) low mass X-ray binary (LMXB), MXB 1659-298, during its most recent outburst in 2015-2017. We detected 60 absorption lines during the soft state (of which 21 at more than 3 $\sigma$), that disappeared in the hard state (e.g., the Fe xxv and Fe xxvi lines). The absorbing plasma is at rest, likely part of the accretion disc atmosphere. The bulk of the absorption features can be reproduced by a high column density ($log(N_H/cm^{-2})\sim23.5$) of highly ionised ($log(\xi/erg~cm~s^{-1})\sim3.8$) plasma. Its disappearance during the hard state is likely the consequence of a thermal photo-ionisation instability. MXB 1659-298's continuum emission can be described by the sum of an absorbed disk black body and its Comptonised emission, plus a black body component. The observed spectral evolution with state is in line with that typically observed in atoll and stellar mass black hole LMXB. The presence of a relativistic Fe K$\alpha$ disk-line is required during the soft state. We also tentatively detect the Fe xxii doublet, whose ratio suggests an electron density of the absorber of $n_e>10^{13} cm^{-3}$, hence, the absorber is likely located at $<7\times10^4 r_g$ from the illuminating source, well inside the Compton and outer disc radii. MXB 1659-298 is the third well monitored atoll LMXB showcasing intense Fe xxv and Fe xxvi absorption during the soft state that disappears during the hard state., Comment: MNRAS in press

We present a population study of low- and intermediate-mass X-ray binaries (LMXBs) with neutron star accretors, performed using the detailed 1D stellar evolution code MESA. We identify all plausible Roche-lobe overflowing binaries at the start of mass transfer, and compare our theoretical mass transfer tracks to the population of well-studied Milky Way LMXBs. The mass transferring evolution depends on the accepted magnetic braking (MB) law for angular momentum loss. The most common MB prescription ("Skumanich MB") originated from observations of the time-dependence of rotational braking of Sun-type stars, where the angular momentum loss rate depends on the donor mass $M_d$, donor radius $R_d$, and rotation rate $\Omega$, $\dot{J} \propto M_d R_d^{\gamma} \Omega^3$. The functional form of the Skumanich MB can be also obtained theoretically assuming a radial magnetic field, isotropic isothermal winds, and boosting of the magnetic field by rotation. Here we show that this simple form of the Skumanich MB law gives mass transfer rates an order of magnitude too weak to explain most observed persistent LMXBs. This failure suggests that the standard Skumanich MB law should not be employed to interpret Galactic, or extragalactic, LMXB populations, with either detailed stellar codes or rapid binary population synthesis codes. We investigate modifications for the MB law, and find that including a scaling of the magnetic field strength with the convective turnover time, and a scaling of MB with the wind mass loss rate, can reproduce persistent LMXBs, and does a better job at reproducing transient LMXBs., Comment: 20 pages, 10 figures, accepted to MNRAS

The X-ray light-curves of the recurring outbursts observed in low-mass X-ray binaries provide strong test beds for constraining (still) poorly understood disc-accretion processes. These light-curves act as a powerful diagnostic to probe the physics behind the mechanisms driving mass inflow and outflow in these binary systems. We have thus developed an innovative methodology, combining a foundation of Bayesian statistics, observed X-ray light-curves, and accretion disc theory. With this methodology, we characterize the angular-momentum (and mass) transport processes in an accretion disc, as well as the properties of the X-ray irradiation-heating that regulates the decay from outburst maximum in low-mass X-ray transients. We recently applied our methodology to the Galactic black-hole low-mass X-ray binary population, deriving from their lightcurves the first-ever quantitative measurements of the $\alpha$-viscosity parameter in these systems \citep{tetarenko2018}. In this paper, we continue the study of these binaries, using Bayesian methods to investigate the X-ray irradiation of their discs during outbursts of strong accretion. We find that the predictions of the disc-instability model, assuming a source of X-ray irradiation proportional to the central accretion rate throughout outburst, do not adequately describe the later stages of BH-LMXB outburst light-curves. We postulate that the complex and varied light-curve morphology observed across the population is evidence for irradiation that varies in time and space within the disc, throughout individual transient outbursts. Lastly, we demonstrate the robustness of our methodology, by accurately reproducing the synthetic model light-curves computed from numerical codes built to simulate accretion flows in binary systems., Comment: 16 pages, 7 Figures, Accepted to MNRAS

We present a new Chandra observation (performed in July 2016) of the neutron star X-ray transient IGR J17480-2446, located in the globular cluster Terzan 5. We study the continued cooling of the neutron star crust in this system that was heated during the 2010 outburst of the source. This new observation was performed two years after the last observation of IGR J17480-2446, hence, significantly extending the cooling baseline. We reanalysed all available Chandra observations of the source (but excluding observations during which one of the known transients in Terzan 5 was in outburst) and fitted the obtained cooling curve with our cooling code NSCool, which allows for much improved modelling than what was previously performed for the source. The data and our fit models indicate that the crust was still cooling ~5.5 years after the outburst ended. The neutron star crust has likely not reached crust-core thermal equilibrium yet, and further cooling is predicted (which can be confirmed with additional Chandra observations in >5 years). Intriguingly, we find indications that the thermal conductivity might be relatively low in part of the crust compared to what has been inferred for other crust-cooling sources and tentatively suggest that this layer might be located around the neutron drip. The reason for this difference is unclear, but might be related to the fact that IGR J17480-2446 harbours a relatively slowly rotating neutron star (with a spin of 11 Hz) that has a relatively strong inferred surface magnetic field ($10^{9-10}$ Gauss) compared to what is known or typically assumed for other cooling sources., Comment: 17 pages, 10 figures, 4 tables, accepted for publication in MNRAS

X-ray spectra of quiescent low-mass X-ray binaries containing neutron stars can be fit with atmosphere models to constrain the mass and the radius. Mass-radius constraints can be used to place limits on the equation of state of dense matter. We perform fits to the X-ray spectrum of a quiescent neutron star in the globular cluster M13, utilizing data from ROSAT, Chandra and XMM-Newton, and constrain the mass-radius relation. Assuming an atmosphere composed of hydrogen and a 1.4${\rm M}_{\odot}$ neutron star, we find the radius to be $R_{\rm NS}=12.2^{+1.5}_{-1.1}$ km, a significant improvement in precision over previous measurements. Incorporating an uncertainty on the distance to M13 relaxes the radius constraints slightly and we find $R_{\rm NS}=12.3^{+1.9}_{-1.7}$ km (for a 1.4${\rm M}_{\odot}$ neutron star with a hydrogen atmosphere), which is still an improvement in precision over previous measurements, some of which do not consider distance uncertainty. We also discuss how the composition of the atmosphere affects the derived radius, finding that a helium atmosphere implies a significantly larger radius., Comment: 7 pages, 5 figures, accepted for publication in MNRAS

To confirm the nature of the donor star in the ultra-compact X-ray binary candidate 47 Tuc X9, we obtained optical spectra (3,000$-$10,000 {\AA}) with the Hubble Space Telescope / Space Telescope Imaging Spectrograph. We find no strong emission or absorption features in the spectrum of X9. In particular, we place $3\sigma$ upper limits on the H$\alpha$ and HeII $\lambda 4686$ emission line equivalent widths $-$EW$_{\mathrm{H\alpha}} \lesssim 14$ {\AA} and $-$EW$_{\mathrm{HeII}} \lesssim 9$ {\AA}, respectively. This is much lower than seen for typical X-ray binaries at a similar X-ray luminosity (which, for $L_{\mathrm{2-10 keV}} \approx 10^{33}-10^{34}$ erg s$^{-1}$ is typically $-$EW$_{\mathrm{H\alpha}} \sim 50$ {\AA}). This supports our previous suggestion (by Bahramian et al.) of an H-poor donor in X9. We perform timing analysis on archival far-ultraviolet, $V$ and $I$-band data to search for periodicities. In the optical bands we recover the seven-day superorbital period initially discovered in X-rays, but we do not recover the orbital period. In the far-ultraviolet we find evidence for a 27.2 min period (shorter than the 28.2 min period seen in X-rays). We find that either a neutron star or black hole could explain the observed properties of X9. We also perform binary evolution calculations, showing that the formation of an initial black hole / He-star binary early in the life of a globular cluster could evolve into a present-day system such as X9 (should the compact object in this system indeed be a black hole) via mass-transfer driven by gravitational wave radiation., Comment: 20 pages, 12 figures. Accepted for publication in MNRAS

The X-ray spectra of intermediate polars can be modelled to give a direct measurement of white dwarf mass. Here we fit accretion column models to NuSTAR spectra of three intermediate polars; V709 Cas, NY Lup and V1223 Sgr in order to determine their masses. From fits to 3-78 keV spectra, we find masses of $M_{\rm WD}=0.88^{+0.05}_{-0.04}M_{\odot}$, $1.16^{+0.04}_{-0.02}M_{\odot}$ and $0.75\pm0.02M_{\odot}$ for V709 Cas, NY Lup and V1223 Sgr, respectively. Our measurements are generally in agreement with those determined by previous surveys of intermediate polars, but with typically a factor $\sim2$ smaller uncertainties. This work paves the way for an approved NuSTAR Legacy Survey of white dwarf masses in intermediate polars., Comment: 9 pages, 2 figures, accepted for publication in MNRAS

Recurring outbursts associated with matter flowing onto compact stellar remnants (black-holes, neutron stars, white dwarfs) in close binary systems, provide strong test beds for constraining the poorly understood accretion process. The efficiency of angular momentum (and thus mass) transport in accretion discs, which has traditionally been encoded in the $\alpha$-viscosity parameter, shapes the light-curves of these outbursts. Numerical simulations of the magneto-rotational instability that is believed to be the physical mechanism behind this transport find values of $\alpha \sim 0.1-0.2$ as required from observations of accreting white dwarfs. Equivalent $\alpha$-viscosity parameters have never been estimated in discs around neutron stars or black holes. Here we report the results of an analysis of archival X-ray light-curves of twenty-one black hole X-ray binary outbursts. Applying a Bayesian approach for a model of accretion allows us to determine corresponding $\alpha$-viscosity parameters, directly from the light curves, to be $\alpha \sim$0.2--1. This result may be interpreted either as a strong intrinsic rate of angular momentum transport in the disc, which can only be sustained by the magneto-rotational instability if a large-scale magnetic field threads the disc, or as a direct indication that mass is being lost from the disc through substantial mass outflows strongly shaping the X-ray binary outburst. Furthermore, the lack of correlation between our estimates of $\alpha$-viscosity and accretion state implies that such outflows can remove a significant fraction of disc mass in all black hole X-ray binary accretion states, favouring magnetically-driven winds over thermally-driven winds that require specific radiative conditions., Comment: 28 pages, 4 Figures, Nature in press

We present Karl G. Jansky Very Large Array radio frequency observations of the new accreting millisecond X-ray pulsar (AMXP), IGR J16597$-$3704, located in the globular cluster NGC 6256. With these data, we detect a radio counterpart to IGR J16597$-$3704, and determine an improved source position. Pairing our radio observations with quasi-simultaneous Swift/XRT X-ray observations, we place IGR J16597$-$3704 on the radio -- X-ray luminosity plane, where we find that IGR J16597$-$3704 is one of the more radio-quiet neutron star low-mass X-ray binaries known to date. We discuss the mechanisms that may govern radio luminosity (and in turn jet production and evolution) in AMXPs. Further, we use our derived radio position to search for a counterpart in archival Hubble Space Telescope and Chandra X-ray Observatory data, and estimate an upper limit on the X-ray luminosity of IGR J16597$-$3704 during quiescence., Comment: 12 pages, 4 Figures, Accepted to ApJ

Multi-decade observing campaigns of the globular clusters 47 Tucanae and M15 have led to an outstanding number of discoveries. Here, we report on the latest results of the long-term observations of the pulsars in these two clusters. For most of the pulsars in 47 Tucanae we have measured, among other things, their higher-order spin period derivatives, which have in turn provided stringent constraints on the physical parameters of the cluster, such as its distance and gravitational potential. For M15, we have studied the relativistic spin precession effect in PSR B2127+11C. We have used full-Stokes observations to model the precession effect, and to constrain the system geometry. We find that the visible beam of the pulsar is swiftly moving away from our line of sight and may very soon become undetectable. On the other hand, we expect to see the opposite emission beam sometime between 2041 and 2053., Comment: 4 pages, 1 figure, accepted for publication in the proceedings of IAU Symposium 337 - Pulsar Astrophysics: The Next Fifty Years

EXO 1745-248 is a transient neutron-star low-mass X-ray binary that resides in the globular cluster Terzan 5. We studied the transient during its quiescent state using 18 Chandra observations of the cluster acquired between 2003 and 2016. We found an extremely variable source, with a luminosity variation in the 0.5-10 keV energy range of $\sim3$ orders of magnitude (between $3\times10^{31}$ erg s$^{-1}$ and $2\times10^{34}$ erg s$^{-1}$) on timescales from years down to only a few days. Using an absorbed power-law model to fit its quiescent spectra, we obtained a typical photon index of $\sim1.4$, indicating that the source is even harder than during outburst and much harder than typical quiescent neutron stars if their quiescent X-ray spectra are also described by a single power-law model. This indicates that EXO 1745-248 is very hard throughout the entire observed X-ray luminosity range. At the highest luminosity, the spectrum fits better when an additional (soft) component is added to the model. All these quiescent properties are likely related to strong variability in the low-level accretion rate in the system. However, its extreme variable behavior is strikingly different from the one observed for other neutron star transients that are thought to still accrete in quiescence. We compare our results to these systems. We also discuss similarities and differences between our target and the transitional millisecond pulsar IGR J18245-2452, which also has hard spectra and strong variability during quiescence., Comment: Accepted for publication in MNRAS

We analyze observations of eight quiescent low-mass X-ray binaries in globular clusters and combine them to determine the neutron star mass-radius curve and the equation of state of dense matter. We determine the effect that several uncertainties may have on our results, including uncertainties in the distance, the atmosphere composition, the neutron star maximum mass, the neutron star mass distribution, the possible presence of a hotspot on the neutron star surface, and the prior choice for the equation of state of dense matter. We find that the radius of a 1.4 solar mass neutron star is most likely from 10 to 14 km and that tighter constraints are only possible with stronger assumptions about the nature of the neutron stars, the systematics of the observations, or the nature of dense matter. Strong phase transitions are preferred over other models and interpretations of the data with a Bayes factor of 8 or more, and in this case, the radius is likely smaller than 12 km. However, radii larger than 12 km are preferred if the neutron stars have uneven temperature distributions., Comment: 16 pages, 12 figures

We present near-infrared (NIR) imaging observations of three transient neutron star X-ray binaries, SAX J1753.5-2349, SAX J1806.5-2215 and AX J1754.2-2754. All three sources are members of the class of `very faint' X-ray transients which exhibit X-ray luminosities $L_X\lesssim10^{36}$ erg s$^{-1}$. The nature of this class of sources is still poorly understood. We detect NIR counterparts for all three systems and perform multi-band photometry for both SAX J1753.5-2349 and SAX J1806.5-2215, including narrow-band Br$_{\gamma}$ photometry for SAX J1806.5-2215. We find that SAX J1753.5-2349 is significantly redder than the field population, indicating that there may be absorption intrinsic to the system, or perhaps a jet is contributing to the infrared emission. SAX J1806.5-2215 appears to exhibit absorption in Br$_{\gamma}$, providing evidence for hydrogen in the system. Our observations of AX J1754.2--2754 represent the first detection of a NIR counterpart for this system. We find that none of the measured magnitudes are consistent with the expected quiescent magnitudes of these systems. Assuming that the infrared radiation is dominated by either the disc or the companion star, the observed magnitudes argue against an ultracompact nature for all three systems., Comment: 10 pages, 10 figures, accepted for publication in MNRAS

This paper is the second in a series where we report the results of the long-term timing of the millisecond pulsars (MSPs) in 47 Tucanae with the Parkes 64-m radio telescope. We obtain improved timing parameters that provide additional information for studies of the cluster dynamics: a) the pulsar proper motions yield an estimate of the proper motion of the cluster as a whole ($\mu_{\alpha}\, = \, 5.00\, \pm \, 0.14\, \rm mas \, yr^{-1}$, $\mu_{\delta}\, = \, -2.84\, \pm \, 0.12\, \rm mas \, yr^{-1}$) and the motion of the pulsars relative to each other. b) We measure the second spin-period derivatives caused by the change of the pulsar line-of-sight accelerations; 47 Tuc H, U and possibly J are being affected by nearby objects. c) For ten binary systems we now measure changes in the orbital period caused by their acceleration in the gravitational field of the cluster. From all these measurements, we derive a cluster distance no smaller than $\sim\,$4.69 kpc and show that the characteristics of these MSPs are very similar to their counterparts in the Galactic disk. We find no evidence in favour of an intermediate mass black hole at the centre of the cluster. Finally, we describe the orbital behaviour of the four "black widow" systems. Two of them, 47 Tuc J and O, exhibit orbital variability similar to that observed in other such systems, while for 47 Tuc I and R the orbits seem to be remarkably stable. It appears, therefore, that not all "black widows" have unpredictable orbital behaviour., Comment: 21 pages in journal format, 9 figures, 4 tables, accepted for publication in MNRAS, several clarifications made and typos fixed

We present 22 new (+3 confirmed) cataclysmic variables (CVs) in the non core-collapsed globular cluster 47 Tucanae (47 Tuc). The total number of CVs in the cluster is now 43, the largest sample in any globular cluster so far. For the identifications we used near-ultraviolet (NUV) and optical images from the Hubble Space Telescope, in combination with X-ray results from the Chandra X-ray Observatory. This allowed us to build the deepest NUV CV luminosity function of the cluster to date. We found that the CVs in 47 Tuc are more concentrated towards the cluster center than the main sequence turnoff stars. We compared our results to the CV populations of the core-collapsed globular clusters NGC 6397 and NGC 6752. We found that 47 Tuc has fewer bright CVs per unit mass than those two other clusters. That suggests that dynamical interactions in core-collapsed clusters play a major role creating new CVs. In 47 Tuc, the CV population is probably dominated by primordial and old dynamically formed systems. We estimated that the CVs in 47 Tuc have total masses of approx. 1.4 M_sun. We also found that the X-ray luminosity function of the CVs in the three clusters is bimodal. Additionally, we discuss a possible double degenerate system and an intriguing/unclassified object. Finally, we present four systems that could be millisecond pulsar companions given their X-ray and NUV/optical colors. For one of them we present very strong evidence for being an ablated companion. The other three could be CO- or He-WDs., Comment: Published on MNRAS. 31 Pages, 23 Figures, 5 Tables. Minor changes with respect to previous arXiv version

For the past couple of decades, the Parkes radio telescope has been regularly observing the millisecond pulsars in 47 Tucanae (47 Tuc). This long-term timing program was designed to address a wide range of scientific issues related to these pulsars and the globular cluster where they are located. In this paper, the first of a series, we address one of these objectives: the characterization of four previously known binary pulsars for which no precise orbital parameters were known, namely 47 Tuc P, V, W and X (pulsars 47 Tuc R and Y are discussed elsewhere). We determined the previously unknown orbital parameters of 47 Tuc V and X and greatly improved those of 47 Tuc P and W. For pulsars W and X we obtained, for the first time, full coherent timing solutions across the whole data span, which allowed a much more detailed characterization of these systems. 47 Tuc W, a well-known tight eclipsing binary pulsar, exhibits a large orbital period variability, as expected for a system of its class. 47 Tuc X turns out to be in a wide, extremely circular, 10.9-day long binary orbit and its position is ~3.8 arcmin away from the cluster center, more than three times the distance of any other pulsar in 47 Tuc. These characteristics make 47 Tuc X a very different object with respect to the other pulsars of the cluster., Comment: Accepted for publication by MNRAS, 18 pages, 11 figures

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., Comment: 10 pages, 3 tables, 5 figures, accepted for publication in MNRAS

Variable emission from Sgr~A*, the luminous counterpart to the super-massive black hole at the center of our Galaxy, arises from the innermost portions of the accretion flow. Better characterization of the variability is important for constraining models of the low-luminosity accretion mode powering Sgr~A*, and could further our ability to use variable emission as a probe of the strong gravitational potential in the vicinity of the $4\times10^{6}\mathrm{M}_{\odot}$ black hole. We use the \textit{Herschel} Spectral and Photometric Imaging Receiver (SPIRE) to monitor Sgr~A* at wavelengths that are difficult or impossible to observe from the ground. We find highly significant variations at 0.25, 0.35, and 0.5 mm, with temporal structure that is highly correlated across these wavelengths. While the variations correspond to $<$1% changes in the total intensity in the \textit{Herschel} beam containing Sgr~A*, comparison to independent, simultaneous observations at 0.85 mm strongly supports the reality of the variations. The lowest point in the light curves, $\sim$0.5 Jy below the time-averaged flux density, places a lower bound on the emission of Sgr~A* at 0.25 mm, the first such constraint on the THz portion of the SED. The variability on few hour timescales in the SPIRE light curves is similar to that seen in historical 1.3 mm data, where the longest time series is available, but the distribution of variations in the sub-mm do not show a tail of large-amplitude variations seen at 1.3 mm. Simultaneous X-ray photometry from XMM-Newton shows no significant variation within our observing period, which may explain the lack of very large variations if X-ray and submillimeter flares are correlated., Comment: Accepted for publication in ApJ

The observed relation between the X-ray and radio properties of low-luminosity accreting black holes has enabled the identification of multiple candidate black hole X-ray binaries (BHXBs) in globular clusters. Here we report an identification of the radio source VLA J213002.08+120904 (aka M15 S2), recently reported in Kirsten et al. 2014, as a BHXB candidate. They showed that the parallax of this flat-spectrum variable radio source indicates a 2.2$^{+0.5}_{-0.3}$ kpc distance, which identifies it as lying in the foreground of the globular cluster M15. We determine the radio characteristics of this source, and place a deep limit on the X-ray luminosity of $\sim4\times10^{29}$ erg s$^{-1}$. Furthermore, we astrometrically identify a faint red stellar counterpart in archival Hubble images, with colors consistent with a foreground star; at 2.2 kpc its inferred mass is 0.1-0.2 $M_{\odot}$. We rule out that this object is a pulsar, neutron star X-ray binary, cataclysmic variable, or planetary nebula, concluding that VLA J213002.08+120904 is the first accreting black hole X-ray binary candidate discovered in quiescence outside a globular cluster. Given the relatively small area over which parallax studies of radio sources have been performed, this discovery suggests a much larger population of quiescent BHXBs in our Galaxy, $2.6\times10^4-1.7\times10^8$ BHXBs at $3\sigma$ confidence, than has been previously estimated ($\sim10^2-10^4$) through population synthesis., Comment: 13 pages, 6 Figures, Accepted to ApJ

We present the results of VLA, ATCA, and Swift XRT observations of the 2015 outburst of the transient neutron star X-ray binary (NSXB), EXO 1745$-$248, located in the globular cluster Terzan 5. Combining (near-) simultaneous radio and X-ray measurements we measure a correlation between the radio and X-ray luminosities of $L_R\propto L_X^\beta$ with $\beta=1.68^{+0.10}_{-0.09}$, linking the accretion flow (probed by X-ray luminosity) and the compact jet (probed by radio luminosity). While such a relationship has been studied in multiple black hole X-ray binaries (BHXBs), this work marks only the third NSXB with such a measurement. Constraints on this relationship in NSXBs are strongly needed, as comparing this correlation between different classes of XB systems is key in understanding the properties that affect the jet production process in accreting objects. Our best fit disc-jet coupling index for EXO 1745$-$248 is consistent with the measured correlation in NSXB 4U 1728$-$34 ($\beta=1.5\pm 0.2$) but inconsistent with the correlation we fit using the most recent measurements from the literature of NSXB Aql X-1 ($\beta=0.76^{+0.14}_{-0.15}$). While a similar disc-jet coupling index appears to hold across multiple BHXBs in the hard accretion state, this does not appear to be the case with the three NSXBs measured so far. Additionally, the normalization of the EXO 1745$-$248 correlation is lower than the other two NSXBs, making it one of the most radio faint XBs ever detected in the hard state. We also report the detection of a type-I X-ray burst during this outburst, where the decay timescale is consistent with hydrogen burning., Comment: 11 pages, 4 figures, Accepted to MNRAS

With the advent of more sensitive all-sky instruments, the transient Universe is being probed in greater depth than ever before. Taking advantage of available resources, we have established a comprehensive database of black hole (and black hole candidate) X-ray binary (BHXB) activity between 1996 and 2015 as revealed by all-sky instruments, scanning surveys, and select narrow-field X-ray instruments aboard the INTErnational Gamma-Ray Astrophysics Laboratory (INTEGRAL), Monitor of All-Sky X-ray Image (MAXI), Rossi X-ray Timing Explorer (RXTE), and Swift telescopes; the Whole-sky Alberta Time-resolved Comprehensive black-Hole Database Of the Galaxy or WATCHDOG. Over the past two decades, we have detected 132 transient outbursts, tracked and classified behavior occurring in 47 transient and 10 persistently accreting BHs, and performed a statistical study on a number of outburst properties across the Galactic population. We find that outbursts undergone by BHXBs that do not reach the thermally dominant accretion state make up a substantial fraction ($\sim$ 40%) of the Galactic transient BHXB outburst sample over the past $\sim20$ years. Our findings suggest that this "hard-only" behavior, observed in transient and persistently accreting BHXBs, is neither a rare nor recent phenomenon and may be indicative of an underlying physical process, relatively common among binary BHs, involving the mass-transfer rate onto the BH remaining at a low level rather than increasing as the outburst evolves. We discuss how the larger number of these "hard-only" outbursts and detected outbursts in general have significant implications for both the luminosity function and mass-transfer history of the Galactic BHXB population., Comment: Accepted to ApJS, 104 pages, 30 Figures, 17 Tables

We report the detection of steady radio emission from the known X-ray source X9 in the globular cluster 47 Tuc. With a double-peaked C IV emission line in its ultraviolet spectrum providing a clear signature of accretion, this source had been previously classified as a cataclysmic variable. In deep ATCA imaging from 2010 and 2013, we identified a steady radio source at both 5.5 and 9.0 GHz, with a radio spectral index (defined as $S_{\nu}\propto\nu^{\alpha}$) of $\alpha=-0.4\pm0.4$. Our measured flux density of $42\pm4$ microJy/beam at 5.5 GHz implies a radio luminosity ($\nu L_{\nu}$) of 5.8e27 erg/s, significantly higher than any previous radio detection of an accreting white dwarf. Transitional millisecond pulsars, which have the highest radio-to-X-ray flux ratios among accreting neutron stars (still a factor of a few below accreting black holes at the same X-ray luminosity), show distinctly different patterns of X-ray and radio variability than X9. When combined with archival X-ray measurements, our radio detection places 47 Tuc X9 very close to the radio/X-ray correlation for accreting black holes, and we explore the possibility that this source is instead a quiescent stellar-mass black hole X-ray binary. The nature of the donor star is uncertain; although the luminosity of the optical counterpart is consistent with a low-mass main sequence donor star, the mass transfer rate required to produce the high quiescent X-ray luminosity of 1e33 erg/s suggests the system may instead be ultracompact, with an orbital period of order 25 minutes. This is the fourth quiescent black hole candidate discovered to date in a Galactic globular cluster, and the only one with a confirmed accretion signature from its optical/ultraviolet spectrum., Comment: 15 pages, 6 figures, accepted for publication in MNRAS

We report the discovery of the likely white dwarf companions to radio millisecond pulsars 47 Tuc Q and 47 Tuc S in the globular cluster 47 Tucanae. These blue stars were found in near-ultraviolet images from the Hubble Space Telescope for which we derived accurate absolute astrometry, and are located at positions consistent with the radio coordinates to within 0.016 arcsec (0.2sigma). We present near-ultraviolet and optical colours for the previously identified companion to millisecond pulsar 47 Tuc U, and we unambiguously confirm the tentative prior identifications of the optical counterparts to 47 Tuc T and 47 Tuc Y. For the latter, we present its radio-timing solution for the first time. We find that all five near-ultraviolet counterparts have U300-B390 colours that are consistent with He white dwarf cooling models for masses ~0.16-0.3 Msun and cooling ages within ~0.1-6 Gyr. The Ha-R625 colours of 47 Tuc U and 47 Tuc T indicate the presence of a strong Ha absorption line, as expected for white dwarfs with an H envelope., Comment: Accepted for publication on MNRAS. 12 pages, 5 figures, 3 tables

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., Comment: 11 pages, 6 figures, 5 tables, MNRAS in press. Minor (textual) changes to match published version

We measure the duty cycles for an existing sample of well observed, nearby dwarf novae using data from AAVSO, and present a quantitative empirical relation between the duty cycle of dwarf novae outbursts and the X-ray luminosity of the system in quiescence. We have found that $\log DC=0.63(\pm0.21)\times(\log L_{X}({\rm erg\,s^{-1}})-31.3)-0.95(\pm0.1)$, where DC stands for duty cycle. We note that there is intrinsic scatter in this relation greater than what is expected from purely statistical errors. Using the dwarf nova X-ray luminosity functions from \citet{Pretorius12} and \citet{Byckling10}, we compare this relation to the number of dwarf novae in the Galactic Bulge Survey which were identified through optical outbursts during an 8-day long monitoring campaign. We find a specific frequency of X-ray bright ($L_{X}>10^{31}\,{\rm erg\,s^{-1}}$) Cataclysmic Variables undergoing Dwarf Novae outbursts in the direction of the Galactic Bulge of $6.6\pm4.7\times10^{-5}\,M_{\odot}^{-1}$. Such a specific frequency would give a Solar neighborhood space density of long period CVs of $\rho=5.6\pm3.9\times10^{-6}\,$pc$^{-3}$. We advocate the use of specific frequency in future work, given that projects like LSST will detect DNe well outside the distance range over which $\rho\approx{\textrm const}$., Comment: 9 pagers, 4 figures Accepted for publication in MNRAS

We present Gemini spectroscopy for 21 candidate optical counterparts to X-ray sources discovered in the Galactic Bulge Survey (GBS). For the majority of the 21 sources, the optical spectroscopy establishes that they are indeed the likely counterparts. One of the criteria we used for the identification was the presence of an Ha emission line. The spectra of several sources revealed an Ha emission line only after careful subtraction of the F or G stellar spectral absorption lines. In a sub-class of three of these sources the residual Halpha emission line is broad (> 400 km/s) which suggests that it is formed in an accretion disk, whereas in other cases the line width is such that we currently cannot determine whether the line emission is formed in an active star/binary or in an accretion disk. GBS source CX377 shows this hidden accretion behaviour most dramatically. The previously-identified broad Ha emission of this source is not present in its Gemini spectra taken about 1 year later. However, broad emission is revealed after subtracting an F6 template star spectrum. The Gemini spectra of three sources (CX446, CX1004, and CXB2) as well as the presence of possible eclipses in light curves of these sources suggest that these sources are accreting binaries viewed under a high inclination., Comment: MNRAS published; author affiliations and references updated

Very faint X-ray binaries (VFXBs), defined as having peak luminosities Lx of 10^34-10^36 erg/s, have been uncovered in significant numbers, but remain poorly understood. We analyse three published outburst light curves of two transient VFXBs using the exponential and linear decay formalism of King and Ritter (1998). The decay timescales and brink luminosities suggest orbital periods of order 1 hour. We review various estimates of VFXB properties, and compare these with suggested explanations of the nature of VFXBs. We suggest that: 1) VFXB outbursts showing linear decays might be explained as partial drainings of the disc of "normal" X-ray transients, and many VFXB outbursts may belong to this category; 2) VFXB outbursts showing exponential decays are best explained by old, short-period systems involving mass transfer from a low-mass white dwarf or brown dwarf; 3) persistent (or quasi-persistent) VFXBs, which maintain an Lx of 10^34-10^35 erg/s for years, may be explained by magnetospheric choking of the accretion flow in a propeller effect, permitting a small portion of the flow to accrete onto the neutron star's surface. We thus predict that (quasi-)persistent VFXBs may also be transitional millisecond pulsars, turning on as millisecond radio pulsars when their Lx drops below 10^32 erg/s., Comment: 12 pages, 4 figures. MNRAS, in press