Globular clusters host a variety of lower-luminosity ($L_X<10^{35}$ erg s$^{-1}$) X-ray sources, including accreting neutron stars and black holes, millisecond pulsars, cataclysmic variables, and chromospherically active binaries. In this paper, we provide a comprehensive catalog of more than 1100 X-ray sources in 38 Galactic globular clusters observed by the Chandra X-ray Observatory's ACIS detector. The targets are selected to complement the MAVERIC survey's deep radio continuum maps of Galactic globular clusters. We perform photometry and spectral analysis for each source, determine a best-fit model, and assess the possibility of it being a foreground/background source based on its spectral properties and location in the cluster. We also provide basic assessments of variability. We discuss the distribution of X-ray binaries in globular clusters, their X-ray luminosity function, and carefully analyze systems with $L_X > 10^{33}$ erg s$^{-1}$. Among these moderately bright systems, we discover a new source in NGC 6539 that may be a candidate accreting stellar-mass black hole or a transitional millisecond pulsar. We show that quiescent neutron star LMXBs in globular clusters may spend ~2% of their lifetimes as transitional millisecond pulsars in their active ($L_X>10^{33}$ erg s$^{-1}$) state. Finally, we identify a substantial under-abundance of bright ($L_X>10^{33}$ erg s$^{-1}$) intermediate polars in globular clusters compared to the Galactic field, in contrast with the literature of the past two decades.

We have identified seven ultraluminous X-ray sources (ULXs) which are coincident with globular cluster candidates (GCs) associated with M87. ULXs in the old GC environment represent a new population of ULXs, and ones likely to be black holes. In this study we perform detailed X-ray spectroscopic followup to seven GC ULXs across a wealth of archival Chandra observations and long time baseline of 16 years. This study brings the total known sample of GC ULXs to 17. Two of these sources show variability in their X-ray luminosity of an order of magnitude over many years, and one of these sources shows intra-observational variability on the scale of hours. While the X-ray spectra of the majority of globular cluster ULXs are best fit by single component models, one of the sources studied in this paper is the second GC ULX to be best fit by a two component model. We compare this new sample of GC ULXs to the previously studied sample, and compare the X-ray and optical properties counterparts across the samples. We find that the clusters that host ULXs in M87 have metallicities from g-z=1.01 to g-z=1.70. The best-fit power-law indices of the X-ray spectra range from Gamma=1.37-2.21, and the best fit inner black-body disk temperatures range from kT=0.56-1.90 keV.

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.

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.

Using a 16.2 hr radio observation by the Australia Telescope Compact Array (ATCA) and archival Chandra data, we found $>5\sigma$ radio counterparts to 4 known and 3 new X-ray sources within the half-light radius ($r_\mathrm{h}$) of the Galactic globular cluster NGC 6397. The previously suggested millisecond pulsar (MSP) candidate, U18, is a steep-spectrum ($S_\nu \propto \nu^\alpha$; $\alpha=-2.0^{+0.4}_{-0.5}$) radio source with a 5.5 GHz flux density of $54.7\pm 4.3~\mathrm{\mu Jy}$. We argue that U18 is most likely a "hidden" MSP that is continuously hidden by plasma shocked at the collision between the winds from the pulsar and companion star. The nondetection of radio pulsations so far is probably the result of enhanced scattering in this shocked wind. On the other hand, we observed 5.5 GHz flux of the known MSP PSR J1740-5340 (U12) to decrease by a factor of $>2.8$ during epochs of 1.4 GHz eclipse, indicating that the radio flux is absorbed in its shocked wind. If U18 is indeed a pulsar whose pulsations are scattered, we note the contrast with U12's flux decrease in eclipse, which argues for two different eclipse mechanisms at the same radio frequency. In addition to U12 and U18, we also found radio associations for 5 other Chandra X-ray sources, four of which are likely background galaxies. The last, U97, which shows strong H$\alpha$ variability, is mysterious; it may be either a quiescent black hole low-mass X-ray binary, or something more unusual.

Using the Very Long Baseline Array and the European Very Long Baseline Interferometry Network we have made a precise measurement of the radio parallax of the black hole X-ray binary MAXI\,J1820+070, providing a model-independent distance to the source. Our parallax measurement of ($0.348\pm0.033$) mas for MAXI J1820+070 translates to a distance of ($2.96\pm0.33$) kpc. This distance implies that the source reached ($15\pm3)\%$ of the Eddington luminosity at the peak of its outburst. Further, we use this distance to refine previous estimates of the jet inclination angle, jet velocity and the mass of the black hole in MAXI J1820+070 to be ($63\pm3)^{\circ}$, ($0.89\pm0.09)c$ and ($9.2\pm1.3) M_{\odot}$, respectively.

(Abridged from Executive Summary) This white paper focuses on the interdisciplinary fields of astrostatistics and astroinformatics, in which modern statistical and computational methods are applied to and developed for astronomical data. Astrostatistics and astroinformatics have grown dramatically in the past ten years, with international organizations, societies, conferences, workshops, and summer schools becoming the norm. Canada's formal role in astrostatistics and astroinformatics has been relatively limited, but there is a great opportunity and necessity for growth in this area. We conducted a survey of astronomers in Canada to gain information on the training mechanisms through which we learn statistical methods and to identify areas for improvement. In general, the results of our survey indicate that while astronomers see statistical methods as critically important for their research, they lack focused training in this area and wish they had received more formal training during all stages of education and professional development. These findings inform our recommendations for the LRP2020 on how to increase interdisciplinary connections between astronomy and statistics at the institutional, national, and international levels over the next ten years. We recommend specific, actionable ways to increase these connections, and discuss how interdisciplinary work can benefit not only research but also astronomy's role in training Highly Qualified Personnel (HQP) in Canada.

RZ2109 is the first of several extragalactic globular clusters shown to host an ultraluminous X-ray source. RZ2109 is particularly notable because optical spectroscopy shows it has broad, luminous [OIII] 4959,5007 emission, while also having no detectable hydrogen emission. The X-ray and optical characteristics of the source in RZ2109 make it a good candidate for being a stellar mass black hole accreting from a white dwarf donor (i.e. an ultracompact black hole X-ray binary). In this paper we present optical spectroscopic monitoring of the [OIII]5007 emission line from 2007 to 2018. We find that the flux of the emission line is significantly lower in recent observations from 2016-2018 than it was in earlier observations in 2007-2011. We also explore the behaviour of the emission line shape over time. Both the core and the wings of the emission line decline over time, with some evidence that the core declines more rapidly than the wings. However, the most recent observations (in 2019) unexpectedly show the emission line core re-brightening.

We use Very Long Baseline Interferometry to measure the proper motions of three black hole X-ray binaries (BHXBs). Using these results together with data from the literature and Gaia-DR2 to collate the best available constraints on proper motion, parallax, distance and systemic radial velocity of 16 BHXBs, we determined their three dimensional Galactocentric orbits. We extended this analysis to estimate the probability distribution for the potential kick velocity (PKV) a BHXB system could have received on formation. Constraining the kicks imparted to BHXBs provides insight into the birth mechanism of black holes (BHs). Kicks also have a significant effect on BH-BH merger rates, merger sites, and binary evolution, and can be responsible for spin-orbit misalignment in BH binary systems. $75\%$ of our systems have potential kicks $>70\,\rm{km~s^{-1}}$. This suggests that strong kicks and hence spin-orbit misalignment might be common among BHXBs, in agreement with the observed quasi-periodic X-ray variability in their power density spectra. We used a Bayesian hierarchical methodology to analyse the PKV distribution of the BHXB population, and suggest that a unimodal Gaussian model with a mean of $107\pm16\,\rm{km~s^{-1}}$ is a statistically favourable fit. Such relatively high PKVs would also reduce the number of BHs likely to be retained in globular clusters. We found no significant correlation between the BH mass and PKV, suggesting a lack of correlation between BH mass and the BH birth mechanism. Our Python code allows the estimation of the PKV for any system with sufficient observational constraints.

MAXI J1621-501 is the first Swift/XRT Deep Galactic Plane Survey transient that was followed up with a multitude of space missions (NuSTAR, Swift, Chandra, NICER, INTEGRAL, and MAXI) and ground-based observatories (Gemini, IRSF, and ATCA). The source was discovered with MAXI on 2017 October 19 as a new, unidentified transient. Further observations with NuSTAR revealed 2 Type I X-ray bursts, identifying MAXI J1621-501 as a Low Mass X-ray Binary (LMXB) with a neutron star primary. Overall, 24 Type I bursts were detected from the source during a 15 month period. At energies below 10 keV, the source spectrum was best fit with three components: an absorbed blackbody with kT = 2.3 keV, a cutoff power law with index $\Gamma{}$ = 0.7, and an emission line centered on 6.3 keV. Timing analysis of the X-ray persistent emission and burst data has not revealed coherent pulsations from the source or an orbital period. We identified, however, a super-orbital period $\sim{}$78 days in the source X-ray light curve. This period agrees very well with the theoretically predicted radiative precession period of $\sim{}$82 days. Thus, MAXI J1621-501 joins a small group of sources characterized with super-orbital periods.

We report on X-ray and radio observations of the ultra-compact X-ray binary 4U 1543-624 taken in August 2017 during an enhanced accretion episode. We obtained NICER monitoring of the source over a $\sim10$ day period during which target-of-opportunity observations were also conducted with Swift, INTEGRAL, and ATCA. Emission lines were measured in the NICER X-ray spectrum at $\sim0.64$ keV and $\sim6.4$ keV that correspond to O and Fe, respectively. By modeling these line components, we are able to track changes in the accretion disk throughout this period. The innermost accretion flow appears to move inwards from hundreds of gravitational radii ($R_{g}=GM/c^{2}$) at the beginning of the outburst to $<8.7$ $R_{g}$ at peak intensity. We do not detect the source in radio, but are able to place a $3\sigma$ upper limit on the flux density at $27$ $\mu$Jy beam$^{-1}$. Comparing the radio and X-ray luminosities, we find that the source lies significantly away from the range typical of black holes in the ${L}_{{r}}$-${L}_{{x}}$ plane, suggesting a neutron star (NS) primary. This adds to the evidence that NSs do not follow a single track in the ${L}_{{r}}$-${L}_{{x}}$ plane, limiting its use in distinguishing between different classes of NSs based on radio and X-ray observations alone.

With the Australian Square Kilometre Array Pathfinder (ASKAP) we monitored the black hole candidate X-ray binary MAXI J1535--571 over seven epochs from 21 September to 2 October 2017. Using ASKAP observations, we studied the HI absorption spectrum from gas clouds along the line-of-sight and thereby constrained the distance to the source. The maximum negative radial velocities measured from the HI absorption spectra for MAXI J1535--571 and an extragalactic source in the same field of view are $-69\pm4$ km s$^{-1}$ and $-89\pm4$ km s$^{-1}$, respectively. This rules out the far kinematic distance ($9.3^{+0.5}_{-0.6}$ kpc), giving a most likely distance of $4.1^{+0.6}_{-0.5}$ kpc, with a strong upper limit of the tangent point at $6.7^{+0.1}_{-0.2}$ kpc. At our preferred distance, the peak unabsorbed luminosity of MAXI J1535--571 was $>78$ per cent of the Eddington luminosity, and shows that the soft-to-hard spectral state transition occurred at the very low luminosity of 1.2 -- 3.4 $\times$ 10$^{-5}$ times the Eddington luminosity. Finally, this study highlights the capabilities of new wide-field radio telescopes to probe Galactic transient outbursts, by allowing us to observe both a target source and a background comparison source in a single telescope pointing.

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.

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.

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.

We present a detailed study of the 2017 eruption of the classical nova ASASSN-17pf (LMCN 2017-11a), which is located in the Large Magellanic Cloud, including data from AAVSO, ASAS-SN, SALT, SMARTS, SOAR, and the Neil Gehrels \textit{Swift} Observatory. The optical light-curve is characterized by multiple maxima (flares) on top of a slowly evolving light-curve (with a decline time, $t_2>$ 100 d). The maxima correlate with the appearance of new absorption line systems in the optical spectra characterized by increasing radial velocities. We suggest that this is evidence of multiple episodes of mass-ejection with increasing expansion velocities. The line profiles in the optical spectra indicate very low expansion velocities (FWHM $\sim$ 190 km s$^{-1}$), making this nova one of the slowest expanding ever observed, consistent with the slowly evolving light-curve. The evolution of the colors and spectral energy distribution show evidence of decreasing temperatures and increasing effective radii for the pseudo-photosphere during each maximum. The optical and infrared light-curves are consistent with dust formation 125 days post-discovery. We speculate that novae showing several optical maxima have multiple mass-ejection episodes leading to shocks that may drive $\gamma$-ray emission and dust formation.

We report on observations of black hole Swift J1357.2-0933, a member of the modest population of very faint X-ray transients. This source has previously shown intense dips in the optical lightcurve, a phenomena that has been linked to the existence of a "unique toroidal structure" in the inner region of the disc, seen at a high inclination. Our observations, carried out by the Neil Gehrels Swift and NuSTAR X-ray observatories, do not show the presence of intense dips in the optical light curves. We find that the X-ray light curves do not show any features that would straightforwardly support an edge-on configuration or high inclination configuration of the orbit. This is similar to what was seen in the X-ray observations of the source during its 2011 outburst. Moreover, the broadband spectra were well described with an absorbed power-law model without any signatures of the cut-off at energies above 10 keV, or any reflection from the disc or the putative torus. Thus, the X-ray data do not support the unique obscuring torus scenario proposed for J1357. We also performed a multi-wavelength study using the data of X-ray telescope and Ultraviolet/Optical Telescope aboard Swift, taken during the 4.5 months duration of the 2017 outburst. This is consistent with what was previously inferred for this source. We found a correlation between the simultaneous X-ray and ultraviolet/optical data and our study suggests that most of the reprocessed flux must be coming out in the ultraviolet.

A number of ultraluminous X-ray sources (ULXs) are physically associated with extragalactic globular clusters (GCs). We undertake a systematic X-ray analysis of eight of the brightest of these sources. We fit the spectra of the GC ULXs to single power law and single disk models. We find that the data never require that any of the sources change between a disk and a power law across successive observations. The GC ULXs best fit by a single disk show a bimodal distribution: they either have temperatures well below 0.5 keV, or variable temperatures ranging above 0.5 keV up to 2~keV. The GC ULXs with low kT have significant changes in luminosity but show little or no change in kT. By contrast, the sources with higher kT either change in both kT and $L_X$ together, or show no significant change in either parameter. Notably, the X-ray characteristics may be related to the optical properties of these ULXs, with the two lowest kT sources showing optical emission lines.

We present an unprecedented, deep study of the primordial low-mass X-ray binary population in an isolated, lower-metallicity environment. We perform followup observations of previously-identified X-ray binary candidates in the Sculptor Dwarf Galaxy by combining a second Chandra observation with Spitzer and Gemini photometry, as well as Gemini spectroscopy of selected targets. Of the original nine bright X-ray sources identified, we are able to classify all but one as quasars, active galactic nuclei, or background galaxies. We further discover four new X-ray sources in the second-epoch Chandra observation. Three of these new sources are background sources and one is a foreground flaring star. We have found that Sculptor is effectively devoid of X-ray sources above a few 1e34 erg/s. If Sculptor is able to retain primordial binaries at a similar rate to globular clusters, this implies that bright X-ray binaries observed in globular clusters in the present epoch are all formed dynamically.

We present the first optical spectroscopy of five confirmed (or strong candidate) redback millisecond pulsar binaries, obtaining complete radial velocity curves for each companion star. The properties of these millisecond pulsar binaries with low-mass, hydrogen-rich companions are discussed in the context of the 14 confirmed and 10 candidate field redbacks. We find that the neutron stars in redbacks have a median mass of 1.78 +/- 0.09 M_sun with a dispersion of sigma = 0.21 +/- 0.09. Neutron stars with masses in excess of 2 M_sun are consistent with, but not firmly demanded by, current observations. Redback companions have median masses of 0.36 +/- 0.04 M_sun with a scatter of sigma = 0.15 +/- 0.04, and a tail possibly extending up to 0.7-0.9 M_sun. Candidate redbacks tend to have higher companion masses than confirmed redbacks, suggesting a possible selection bias against the detection of radio pulsations in these more massive candidate systems. The distribution of companion masses between redbacks and the less massive black widows continues to be strongly bimodal, which is an important constraint on evolutionary models for these systems. Among redbacks, the median efficiency of converting the pulsar spindown energy to gamma-ray luminosity is ~10%.

We improve the dynamical black hole (BH) mass estimates in three nearby low-mass early-type galaxies--NGC 205, NGC 5102, and NGC 5206. We use new \hst/STIS spectroscopy to fit the star formation histories of the nuclei in these galaxies, and use these measurements to create local color--mass-to-light ratio (\ml) relations. We then create new mass models from \hst~imaging and combined with adaptive optics kinematics, we use Jeans dynamical models to constrain their BH masses. The masses of the central BHs in NGC 5102 and NGC 5206 are both below one million solar masses and are consistent with our previous estimates, $9.12_{-1.53}^{+1.84}\times10^5$\Msun~and $6.31_{-2.74}^{+1.06}\times10^5$\Msun~(3$\sigma$ errors), respectively. However, for NGC 205, the improved models suggest the presence of a BH for the first time, with a best-fit mass of $6.8_{-6.7}^{+95.6}\times10^3$\Msun~(3$\sigma$ errors). This is the least massive central BH mass in a galaxy detected using any method. We discuss the possible systematic errors of this measurement in detail. Using this BH mass, the existing upper limits of both X-ray, and radio emissions in the nucleus of NGC 205 suggest an accretion rate $\lesssim$$10^{-5}$ of the Eddington rate. We also discuss the color--\mleff~relations in our nuclei and find that the slopes of these vary significantly between nuclei. Nuclei with significant young stellar populations have steeper color--\mleff~relations than some previously published galaxy color--\mleff~relations.

We present the discovery of ASASSN-18ey (MAXI J1820+070), a new black hole low-mass X-ray binary discovered by the All-Sky Automated Survey for SuperNovae (ASAS-SN). A week after ASAS-SN discovered ASASSN-18ey as an optical transient, it was detected as an X-ray transient by MAXI/GCS. Here, we analyze ASAS-SN and Asteroid Terrestrial-impact Last Alert System (ATLAS) pre-outburst optical light curves, finding evidence of intrinsic variability for several years prior to the outburst. While there was no long-term rise leading to outburst, as has been seen in several other systems, the start of the outburst in the optical preceded that in the X-rays by $7.20\pm0.97~\rm days$. We analyze the spectroscopic evolution of ASASSN-18ey from pre-maximum to $> 100~\rm days$ post-maximum. The spectra of ASASSN-18ey exhibit broad, asymmetric, double-peaked H$\alpha$ emission. The Bowen blend ($\lambda \approx 4650$\AA) in the post-maximum spectra shows highly variable double-peaked profiles, likely arising from irradiation of the companion by the accretion disk, typical of low-mass X-ray binaries. The optical and X-ray luminosities of ASASSN-18ey are consistent with black hole low-mass X-ray binaries, both in outburst and quiescence.

We present multi-frequency monitoring observations of the black hole X-ray binary V404 Cygni throughout its June 2015 outburst. Our data set includes radio and mm/sub-mm photometry, taken with the Karl G. Jansky Very Large Array, Arc-Minute MicroKelvin Imager Large Array, Sub-millimeter Array, James Clerk Maxwell Telescope, and the Northern Extended Millimetre Array, combined with publicly available infrared, optical, UV, and X-ray measurements. With these data, we report detailed diagnostics of the spectral and variability properties of the jet emission observed during different stages of this outburst. These diagnostics show that emission from discrete jet ejecta dominated the jet emission during the brightest stages of the outburst. We find that the ejecta became fainter, slower, less frequent, and less energetic, before the emission transitioned (over 1--2 days) to being dominated by a compact jet, as the outburst decayed toward quiescence. While the broad-band spectrum of this compact jet showed very little evolution throughout the outburst decay (with the optically thick to thin synchrotron jet spectral break residing in the near-infrared/optical bands; $\sim2-5\times10^{14}$ Hz), the emission still remained intermittently variable at mm/sub-mm frequencies. Additionally, we present a comparison between the radio jet emission throughout the 2015 and previous 1989 outbursts, confirming that the radio emission in the 2015 outburst decayed significantly faster than in 1989. Lastly, we detail our sub-mm observations taken during the December 2015 mini-outburst of V404 Cygni, which demonstrate that, similar to the main outburst, the source was likely launching jet ejecta during this short period of renewed activity.

Monitoring the cooling of neutron-star crusts heated during accretion outbursts allows us to infer the physics of the dense matter present in the crust. We examine the crust cooling evolution of the low-mass X-ray binary MXB 1659-29 up to ~505 days after the end of its 2015 outburst (hereafter outburst II) and compare it with what we observed after its previous 1999 outburst (hereafter outburst I) using data obtained from the Swift, XMM-Newton, and Chandra observatories. The observed effective surface temperature of the neutron star in MXB 1659-29 dropped from ~92 eV to ~56 eV from ~12 days to ~505 days after the end of outburst II. The most recently performed observation after outburst II suggests that the crust is close to returning to thermal equilibrium with the core. We model the crust heating and cooling for both its outbursts collectively to understand the effect of parameters that may change for every outburst (e.g., the average accretion rate, the length of outburst, the envelope composition of the neutron star at the end of the outburst) and those which can be assumed to remain the same during these two outbursts (e.g., the neutron star mass, its radius). Our modelling indicates that all parameters were consistent between the two outbursts with no need for any significant changes. In particular, the strength and the depth of the shallow heating mechanism at work (in the crust) were inferred to be the same during both outbursts, contrary to what has been found when modelling the cooling curves after multiple outburst of another source, MAXI J0556-332. This difference in source behaviour is not understood. We discuss our results in the context of our current understanding of cooling of accretion-heated neutron-star crusts, and in particular with respect to the unexplained shallow heating mechanism.

The Fermi $\gamma$-ray source 1FGL J1417.7--4407 (J1417) is a compact X-ray binary with a neutron star primary and a red giant companion in a $\sim$5.4 day orbit. This initial conclusion, based on optical and X-ray data, was confirmed when a 2.66 ms radio pulsar was found at the same location (and with the same orbital properties) as the optical/X-ray source. However, these initial studies found conflicting evidence about the accretion state and other properties of the binary. We present new optical, radio, and X-ray observations of J1417 that allow us to better understand this unusual system. We show that one of the main pieces of evidence previously put forward for an accretion disk---the complex morphology of the persistent H$\alpha$ emission line---can be better explained by the presence of a strong, magnetically driven stellar wind from the secondary and its interaction with the pulsar wind. The radio spectral index derived from VLA/ATCA observations is broadly consistent with that expected from a millisecond pulsar, further disfavoring an accretion disk scenario. X-ray observations show evidence for a double-peaked orbital light curve, similar to that observed in some redback millisecond pulsar binaries and likely due to an intrabinary shock. Refined optical light curve fitting gives a distance of 3.1$\pm$0.6 kpc, confirmed by a Gaia DR2 parallax measurement. At this distance the X-ray luminosity of J1417 is (1.0$^{+0.4}_{-0.3}$) $\times 10^{33}$ erg s$^{-1}$, which is more luminous than all known redback systems in the rotational-powered pulsar state, perhaps due to the wind from the giant companion. The unusual phenomenology of this system and its differing evolutionary path from redback millisecond pulsar binaries points to a new eclipsing pulsar "spider" subclass that is a possible progenitor of normal field millisecond pulsar binaries.

Transitional millisecond pulsars are accreting millisecond pulsars that switch between accreting X-ray binary and millisecond radio pulsar states. Only a handful of these objects have been identified so far. Terzan 5 CX1 is a variable hard X-ray source in the globular cluster Terzan 5. In this paper, we identify a radio counterpart to CX1 in deep Very Large Array radio continuum data. Chandra observations over the last fourteen years indicate that CX1 shows two brightness states: in 2003 and 2016 the source was the brightest X-ray source in the cluster (at L$_X \sim 10^{33}$ erg s$^{-1}$), while in many intermediate observations, its luminosity was almost an order of magnitude lower. We analyze all available X-ray data of CX1, showing that the two states are consistent with the spectral and variability properties observed for the X-ray active and radio pulsar states of known transitional millisecond pulsars. Finally, we discuss the prospects for the detection of CX1 as a radio pulsar in existing timing data.

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.

We present the results of long-term monitoring of the X-ray emission from the ultraluminous X-ray source XMMUJ122939.9+075333 in the extragalactic globular cluster RZ2109. The combination of the high X-ray luminosity, short term X-ray variability, X-ray spectrum, and optical emission suggest that this system is likely an accreting black hole in a globular cluster. To study the long-term behavior of the X-ray emission from this source, we analyze both new and archival Chandra and XMM-Newton observations, covering 16 years from 2000 to 2016. For all of these observations, we fit extracted spectra of RZ2109 with xspec models. The spectra are all dominated by a soft component, which is very soft with typical fit temperatures of T $\simeq$ 0.15 keV. The resulting X-ray fluxes show strong variability on short and long timescales. We also find that the X-ray spectrum often shows no significant change even with luminosity changes as large as a factor of five.

We examine the internal properties of the most massive ultracompact dwarf galaxy (UCD), M59-UCD3, by combining adaptive optics assisted near-IR integral field spectroscopy from Gemini/NIFS, and Hubble Space Telescope (HST) imaging. We use the multi-band HST imaging to create a mass model that suggests and accounts for the presence of multiple stellar populations and structural components. We combine these mass models with kinematics measurements from Gemini/NIFS to find a best-fit stellar mass-to-light ratio ($M/L$) and black hole (BH) mass using Jeans Anisotropic Models (JAM), axisymmetric Schwarzschild models, and triaxial Schwarzschild models. The best fit parameters in the JAM and axisymmetric Schwarzschild models have black holes between 2.5 and 5.9 million solar masses. The triaxial Schwarzschild models point toward a similar BH mass, but show a minimum $\chi^2$ at a BH mass of $\sim 0$. Models with a BH in all three techniques provide better fits to the central $V_{rms}$ profiles, and thus we estimate the BH mass to be $4.2^{+2.1}_{-1.7} \times 10^{6}$ M$_\odot$ (estimated 1$\sigma$ uncertainties). We also present deep radio imaging of M59-UCD3 and two other UCDs in Virgo with dynamical BH mass measurements, and compare these to X-ray measurements to check for consistency with the fundamental plane of BH accretion. We detect faint radio emission in M59cO, but find only upper limits for M60-UCD1 and M59-UCD3 despite X-ray detections in both these sources. The BH mass and nuclear light profile of M59-UCD3 suggests it is the tidally stripped remnant of a $\sim$10$^{9-10}$ M$_\odot$ galaxy.

We report on the results of the multiwavelength campaign carried out after the discovery of the INTEGRAL transient IGR J17329-2731. The optical data collected with the SOAR telescope allowed us to identify the donor star in this system as a late M giant at a distance of 2.7$^{+3.4}_{-1.2}$ kpc. The data collected quasi-simultaneously with XMM-Newton and NuSTAR showed the presence of a modulation with a period of 6680$\pm$3 s in the X-ray light curves of the source. This unveils that the compact object hosted in this system is a slowly rotating neutron star. The broadband X-ray spectrum showed the presence of a strong absorption ($\gg$10$^{23}$ cm$^{-2}$) and prominent emission lines at 6.4 keV, and 7.1 keV. These features are usually found in wind-fed systems, in which the emission lines result from the fluorescence of the X-rays from the accreting compact object on the surrounding stellar wind. The presence of a strong absorption line around $\sim$21 keV in the NuSTAR spectrum suggests a cyclotron origin, thus allowing us to estimate the neutron star magnetic field as $\sim$2.4$\times$10$^{12}$ G. All evidence thus suggests IGR J17329-2731 is a symbiotic X-ray binary. As no X-ray emission was ever observed from the location of IGR J17329-2731 by INTEGRAL (or other X-ray facilities) during the past 15 yr in orbit and considering that symbiotic X-ray binaries are known to be variable but persistent X-ray sources, we concluded that INTEGRAL caught the first detectable X-ray emission from IGR J17329-2731 when the source shined as a symbiotic X-ray binary. The Swift/XRT monitoring performed up to $\sim$3 months after the discovery of the source, showed that it maintained a relatively stable X-ray flux and spectral properties.

We present the discovery and characterization of a radio-bright binary in the Galactic globular cluster M10. First identified in deep radio continuum data from the Karl G. Jansky Very Large Array, M10-VLA1 has a flux density of $27\pm4$ $\mu$Jy at 7.4 GHz and a flat to inverted radio spectrum. Chandra imaging shows an X-ray source with $L_X \approx 10^{31}$ erg s$^{-1}$ matching the location of the radio source. This places M10-VLA1 within the scatter of the radio--X-ray luminosity correlation for quiescent stellar-mass black holes, and a black hole X-ray binary is a viable explanation for this system. The radio and X-ray properties of the source disfavor, though do not rule out, identification as an accreting neutron star or white dwarf system. Optical imaging from the Hubble Space Telescope and spectroscopy from the SOAR telescope show the system has an orbital period of 3.339 d and an unusual "red straggler" component: an evolved star found redward of M10's red giant branch. These data also show UV/optical variability and double-peaked H$\alpha$ emission characteristic of an accretion disk. However, SOAR spectroscopic monitoring reveals that the velocity semi-amplitude of the red straggler is low. We conclude that M10-VLA1 is most likely either a quiescent black hole X-ray binary with a rather face-on ($ i < $ 4$^{\circ}$) orientation or an unusual flaring RS CVn-type active binary, and discuss future observations that could distinguish between these possibilities.

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.

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.

We report the discovery of X-ray pulsations at 105.2 Hz (9.5 ms) from the transient X-ray binary IGR J16597-3704 using NuSTAR and Swift. The source was discovered by INTEGRAL in the globular cluster NGC 6256 at a distance of 9.1 kpc. The X-ray pulsations show a clear Doppler modulation implying an orbital period of ~46 minutes and a projected semi-major axis of ~5 lt-ms, which makes IGR J16597-3704 an ultra-compact X-ray binary system. We estimated a minimum companion mass of 0.0065 solar masses, assuming a neutron star mass of 1.4 solar masses, and an inclination angle of <75 degrees (suggested by the absence of eclipses or dips in its light-curve). The broad-band energy spectrum of the source is well described by a disk blackbody component (kT ~1.4 keV) plus a comptonised power-law with photon index ~2.3 and an electron temperature of ~30 keV. Radio pulsations from the source were searched for with the Parkes observatory and not detected.

We present optical and ultraviolet spectra of the first electromagnetic counterpart to a gravitational wave (GW) source, the binary neutron star merger GW170817. Spectra were obtained nightly between 1.5 and 9.5 days post-merger, using the SOAR and Magellan telescopes; the UV spectrum was obtained with the \textit{Hubble Space Telescope} at 5.5 days. Our data reveal a rapidly-fading blue component ($T\approx5500$ K at 1.5 days) that quickly reddens; spectra later than $\gtrsim 4.5$ days peak beyond the optical regime. The spectra are mostly featureless, although we identify a possible weak emission line at $\sim 7900$ \AA\ at $t\lesssim 4.5$ days. The colours, rapid evolution and featureless spectrum are consistent with a "blue" kilonova from polar ejecta comprised mainly of light $r$-process nuclei with atomic mass number $A\lesssim 140$. This indicates a sight-line within $\theta_{\rm obs}\lesssim 45^{\circ}$ of the orbital axis. Comparison to models suggests $\sim0.03$ M$_\odot$ of blue ejecta, with a velocity of $\sim 0.3c$. The required lanthanide fraction is $\sim 10^{-4}$, but this drops to $<10^{-5}$ in the outermost ejecta. The large velocities point to a dynamical origin, rather than a disk wind, for this blue component, suggesting that both binary constituents are neutron stars (as opposed to a binary consisting of a neutron star and a black hole). For dynamical ejecta, the high mass favors a small neutron star radius of $\lesssim 12$ km. This mass also supports the idea that neutron star mergers are a major contributor to $r$-process nucleosynthesis.

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.

47 Tuc X9 is a low mass X-ray binary (LMXB) in the globular cluster 47 Tucanae, and was previously thought to be a cataclysmic variable. However, Miller-Jones et al. (2015) recently identified a radio counterpart to X9 (inferring a radio/X-ray luminosity ratio consistent with black hole LMXBs), and suggested that the donor star might be a white dwarf. We report simultaneous observations of X9 performed by Chandra, Nustar and Australia Telescope Compact Array. We find a clear 28.18$\pm$0.02 min periodic modulation in the Chandra data, which we identify as the orbital period, confirming this system as an ultracompact X-ray binary. Our X-ray spectral fitting provides evidence for photoionized gas having a high oxygen abundance in this system, which indicates a C/O white dwarf donor. We also identify reflection features in the hard X-ray spectrum, making X9 the faintest LMXB to show X-ray reflection. We detect a $\sim$ 6.8 day modulation in the X-ray brightness by a factor of 10, in archival Chandra, Swift, and Rosat data. The simultaneous radio/X-ray flux ratio is consistent with either a black hole primary or a neutron star primary, if the neutron star is a transitional millisecond pulsar. Considering the measured orbital period (with other evidence of a white dwarf donor), and the lack of transitional millisecond pulsar features in the X-ray light curve, we suggest that this could be the first ultracompact black hole X-ray binary identified in our Galaxy.

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.

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.

We search the literature for reports on the spectral properties of neutron-star low-mass X-ray binaries when they have accretion luminosities between 1E34 and 1E36 ergs/s. We found that in this luminosity range the photon index (obtained from fitting a simple absorbed power-law in the 0.5-10 keV range) increases with decreasing 0.5-10 keV X-ray luminosity (i.e., the spectrum softens). Such behaviour has been reported before for individual sources, but here we demonstrate that very likely most (if not all) neutron-star systems behave in a similar manner and possibly even follow a universal relation. When comparing the neutron-star systems with black-hole systems, it is clear that most black-hole binaries have significantly harder spectra at luminosities of 1E34 - 1E35 erg/s. Despite a limited number of data points, there are indications that these spectral differences also extend to the 1E35 - 1E36 erg/s range. This observed difference between the neutron-star binaries and black-hole ones suggests that the spectral properties (between 0.5-10 keV) at 1E34 - 1E35 erg/s can be used to tentatively determine the nature of the accretor in unclassified X-ray binaries. We discuss our results in the context of properties of the accretion flow at low luminosities and we suggest that the observed spectral differences likely arise from the neutron-star surface becoming dominantly visible in the X-ray spectra. We also suggest that both the thermal component and the non-thermal component might be caused by low-level accretion onto the neutron-star surface for luminosities below a few times 1E34 erg/s.

In quiescent low-mass X-ray binaries (qLMXBs) containing neutron stars, the origin of the thermal X-ray component may be either release of heat from the core of the neutron star, or continuing low-level accretion. In general, heat from the core should be stable on timescales $<10^4$ years, while continuing accretion may produce variations on a range of timescales. While some quiescent neutron stars (e.g. Cen X-4, Aql X-1) have shown variations in their thermal components on a range of timescales, several others, particularly those in globular clusters with no detectable nonthermal hard X-rays (fit with a powerlaw), have shown no measurable variations. Here, we constrain the spectral variations of 12 low mass X-ray binaries in 3 globular clusters over $\sim10$ years. We find no evidence of variations in 10 cases, with limits on temperature variations below 11% for the 7 qLMXBs without powerlaw components, and limits on variations below 20% for 3 other qLMXBs that do show non-thermal emission. However, in 2 qLMXBs showing powerlaw components in their spectra (NGC 6440 CX 1 & Terzan 5 CX 12) we find marginal evidence for a 10% decline in temperature, suggesting the presence of continuing low-level accretion. This work adds to the evidence that the thermal X-ray component in quiescent neutron stars without powerlaw components can be explained by heat deposited in the core during outbursts. Finally, we also investigate the correlation between hydrogen column density (N$_H$) and optical extinction (A$_V$) using our sample and current models of interstellar X-ray absorption, finding $N_H ({\rm cm}^{-2}) = (2.81\pm0.13)\times10^{21} A_V$.

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.

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.

We report observations using the Swift/XRT, NuSTAR, and Chandra X-ray telescopes of the transient X-ray source CXOGC J174540.0-290005, during its 2013 outburst. Due to its location in the field of multiple observing campaigns targeting Sgr A*, this is one of the best-studied outbursts of a very faint X-ray binary (VFXB; peak $L_X<10^{36}$ erg/s) yet recorded, with detections in 173 ks of X-ray observations over 50 days. VFXBs are of particular interest, due to their unusually low outburst luminosities and time-averaged mass transfer rates, which are hard to explain within standard accretion physics and binary evolution. The 2013 outburst of CXOGC J174540.0-290005 peaked at $L_X$(2-10 keV)=$5.0\times10^{35}$ erg/s, and all data above $10^{34}$ ergs/s were well-fit by an absorbed power-law of photon index $\sim1.7$, extending from 2 keV out to $\sim$70 keV. We discuss the implications of these observations for the accretion state of CXOGC J174540.0-290005.

We use K-band spectroscopy of the counterpart to the rapidly variable X-ray transient XMMU J174445.5-295044 to identify it as a new symbiotic X-ray binary. XMMU J174445.5-295044 has shown a hard X-ray spectrum (we verify its association with an Integral/IBIS 18-40 keV detection in 2013 using a short Swift/XRT observation), high and varying N$_H$, and rapid flares on timescales down to minutes, suggesting wind accretion onto a compact star. We observed its near-infrared counterpart using the Near-infrared Integral Field Spectrograph (NIFS) at Gemini-North, and classify the companion as ~ M2 III. We infer a distance of $3.1^{+1.8}_{-1.1}$ kpc (conservative 1-sigma errors), and therefore calculate that the observed X-ray luminosity (2-10 keV) has reached to at least 4$\times10^{34}$ erg/s. We therefore conclude that the source is a symbiotic X-ray binary containing a neutron star (or, less likely, black hole) accreting from the wind of a giant.

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

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

The high stellar densities in the cores of globular clusters cause significant stellar interactions. These stellar interactions can produce close binary mass-transferring systems involving compact objects and their progeny, such as X-ray binaries and radio millisecond pulsars. Comparing the numbers of these systems and interaction rates in different clusters drives our understanding of how cluster parameters affect the production of close binaries. In this paper we estimate stellar encounter rates (Gamma) for 124 Galactic globular clusters based on observational data as opposed to the methods previously employed, which assumed "King-model" profiles for all clusters. By deprojecting cluster surface brightness profiles to estimate luminosity density profiles, we treat "King-model" and "core-collapsed" clusters in the same way. In addition, we use Monte-Carlo simulations to investigate the effects of uncertainties in various observational parameters (distance, reddening, surface brightness) on Gamma, producing the first catalog of GC stellar encounter rates with estimated errors. Comparing our results with published observations of likely products of stellar interactions (numbers of X-ray binaries, numbers of radio millisecond pulsars, and gamma-ray luminosity) we find both clear correlations and some differences with published results.

We observed the globular cluster NGC 6652 with Chandra for 47.5 ks, detecting six known X-ray sources, as well as five previously undetected X-ray sources. Source A (XB 1832-330) is a well-known bright low-mass X-ray binary (LXMB). The second brightest source, B, has a spectrum that fits well to either a power-law model (Gamma ~ 1.3) or an absorbed hot gas emission model (kT ~ 34 keV). Its unabsorbed 0.5-10 keV luminosity (L_X = 1.6+-0.1*10^34 erg/s) is suggestive of a neutron star primary; however, Source B exhibits unusual variability for a LMXB, varying by over an order of magnitude on timescales of ~ 100 s. Source C's spectrum contains a strong low-energy component below 1 keV. Its spectrum is well fit to a simplified magnetic cataclysmic variable (CV) model, thus the soft component may be explained by a hot polar cap of a magnetic CV. Source D has an average L_X (0.5-10 keV) ~ 9*10^32 erg/s, and its spectrum is well fit to a neutron star atmosphere model. This is indicative of a quiescent neutron star LXMB, suggesting Source D may be the third known LMXB in NGC 6652. Source E has L_X (0.5-10 keV) ~ 3*10^32 erg/s, while Source F has L_X (0.5-10 keV) ~ 1*10^32 erg/s. Their relatively hard X-ray spectra are well-fit by power-law or plasma emission models. Five newly detected fainter sources have luminosities between 1-5*10^31 erg/s. NGC 6652 has an unusually flat X-ray luminosity function compared to other globular clusters, which may be connected to its extremely high central density.