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.

Black holes in binary systems execute patterns of outburst activity where two characteristic X-ray states are associated with different behaviours observed at radio wavelengths. The hard state is associated with radio emission indicative of a continuously replenished, collimated, relativistic jet, whereas the soft state is rarely associated with radio emission, and never continuously, implying the absence of a quasi-steady jet. Here we report radio observations of the black hole transient MAXI J1820$+$070 during its 2018 outburst. As the black hole transitioned from the hard to soft state we observed an isolated radio flare, which, using high angular resolution radio observations, we connect with the launch of bi-polar relativistic ejecta. This flare occurs as the radio emission of the core jet is suppressed by a factor of over 800. We monitor the evolution of the ejecta over 200 days and to a maximum separation of 10$''$, during which period it remains detectable due to in-situ particle acceleration. Using simultaneous radio observations sensitive to different angular scales we calculate an accurate estimate of energy content of the approaching ejection. This energy estimate is far larger than that derived from state transition radio flare, suggesting a systematic underestimate of jet energetics.

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.

This white paper gives a brief summary of the time domain science that has been performed with the JCMT in recent years and highlights the opportunities for continuing work in this field over the next decade. The main focus of this document is the JCMT Transient Survey, a large program initiated in 2015 to measure the frequency and amplitude of variability events associated with protostars in nearby star-forming regions. After summarising the major accomplishments so far, an outline is given for extensions to the current survey, featuring a discussion on what will be possible with the new 850 micron camera that is expected to be installed in late 2022. We also discuss possible applications of submillimetre monitoring to active galactic nuclei, X-ray binaries, asymptotic giant branch stars, and flare stars.

MAXI J1535-571 is a Galactic black hole candidate X-ray binary that was discovered going into outburst in 2017 September. In this paper, we present comprehensive radio monitoring of this system using the Australia Telescope Compact Array (ATCA), as well as the MeerKAT radio observatory, showing the evolution of the radio jet during its outburst. Our radio observations show the early rise and subsequent quenching of the compact jet as the outburst brightened and then evolved towards the soft state. We constrain the compact jet quenching factor to be more than 3.5 orders of magnitude. We also detected and tracked (for 303 days) a discrete, relativistically-moving jet knot that was launched from the system. From the motion of the apparently superluminal knot, we constrain the jet inclination (at the time of ejection) and speed to $\leq 45^{\circ}$ and $\geq0.69$c, respectively. Extrapolating its motion back in time, our results suggest that the jet knot was ejected close in time to the transition from the hard intermediate state to soft intermediate state. The launching event also occurred contemporaneously with a short increase in X-ray count rate, a rapid drop in the strength of the X-ray variability, and a change in the type-C quasi-periodic oscillation (QPO) frequency that occurs $>$2.5 days before the first appearance of a possible type-B QPO.

Powerful relativistic jets are one of the main ways in which accreting black holes provide kinetic feedback to their surroundings. Jets launched from or redirected by the accretion flow that powers them should be affected by the dynamics of the flow, which in accreting stellar-mass black holes has shown increasing evidence for precession due to frame dragging effects that occur when the black hole spin axis is misaligned with the orbital plane of its companion star. Recently, theoretical simulations have suggested that the jets can exert an additional torque on the accretion flow, although the full interplay between the dynamics of the accretion flow and the launching of the jets is not yet understood. Here we report a rapidly changing jet orientation on a timescale of minutes to hours in the black hole X-ray binary V404 Cygni, detected with very long baseline interferometry during the peak of its 2015 outburst. We show that this can be modelled as Lense-Thirring precession of a vertically-extended slim disk that arises from the super-Eddington accretion rate. Our findings suggest that the dynamics of the precessing inner accretion disk could play a role in either directly launching or redirecting the jets within the inner few hundred gravitational radii. Similar dynamics should be expected in any strongly-accreting black hole whose spin is misaligned with the inflowing gas, both affecting the observational characteristics of the jets, and distributing the black hole feedback more uniformly over the surrounding environment.

The main outburst of the candidate black hole low-mass X-ray binary (BH LMXB) MAXI J1535-571 ended in 2018 May and was followed by at least five episodes of re-brightenings. We have monitored this re-brightening phenomenon at X-ray and radio wavelengths using the {\it Neil Gehrels Swift Observatory} and Australia Telescope Compact Array, respectively. The first two re-brightenings exhibited a high peak X-ray luminosity (implying a high mass accretion rate) and were observed to transition from the hard to the soft state. However, unlike the main outburst, these re-brightenings did not exhibit clear hysteresis. During the re-brightenings, when MAXI J1535-571 was in the hard state, we observed the brightening of a compact radio jet which was subsequently quenched when the source transitioned to a similar soft state as was observed during the main outburst. We report on the first investigation of disk-jet coupling over multiple rapidly evolving re-brightenings in a BH LMXB. We find that the accretion flow properties and the accompanying compact jet evolve on a similarly rapid time scale of ~days rather than the typical value of ~weeks as observed for most other BH LMXBs during their main outburst events.

This paper outlines the importance of understanding jets from compact binaries for the problem of understanding the broader phenomenology of jet production. Because X-ray binaries are nearby and bright, have well-measured system parameters, and vary by factors of $\sim 10^6$ on $\sim$ year timescales, they provide a unique opportunity to understand how various aspects of the jet physics change in response to changes in the accretion flow, giving the possibility of looking for trends within individual systems and testing their universality with other systems, rather than trying to interpret large samples of objects on a statistical basis.

We present simultaneous multi-band radio and X-ray observations of the black hole X-ray binary Cygnus X-1, taken with the Karl G. Jansky Very Large Array and the Nuclear Spectroscopic Telescope Array. With these data, we detect clear flux variability consistent with emission from a variable compact jet. To probe how the variability signal propagates down the jet flow, we perform detailed timing analyses of our data. We find that the radio jet emission shows no significant power at Fourier frequencies $f\gtrsim0.03$ Hz (below $\sim30$ sec timescales), and that the higher frequency radio bands (9/11 GHz) are strongly correlated over a range of timescales, displaying a roughly constant time lag with Fourier frequency of a few tens of seconds. However, in the lower frequency radio bands (2.5/3.5 GHz) we find a significant loss of coherence over the same range of timescales. Further, we detect a correlation between the X-ray/radio emission, measuring time lags between the X-ray/radio bands on the order of tens of minutes. We use these lags to solve for the compact jet speed, finding that the Cyg X-1 jet is more relativistic than usually assumed for compact jets, where $\beta=0.92^{+0.03}_{-0.06}$, ($\Gamma=2.59^{+0.79}_{-0.61}$). Lastly, we constrain how the jet size scale changes with frequency, finding a shallower relation ($\propto \nu^{-0.4}$) than predicted by simple jet models ($\propto \nu^{-1}$), and estimate a jet opening angle of $\phi\sim0.4-1.8$ degrees. With this study, we have developed observational techniques designed to overcome the challenges of radio timing analyses and created the tools needed to connect rapid radio jet variability properties to internal jet physics.

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.

We report on the results of optical, near-infrared (NIR) and mid-infrared observations of the black hole X-ray binary candidate (BHB) MAXI J1535-571 during its 2017/2018 outburst. During the first part of the outburst (MJD 58004-58012), the source shows an optical-NIR spectrum that is consistent with an optically thin synchrotron power-law from a jet. After MJD 58015, however, the source faded considerably, the drop in flux being much more evident at lower frequencies. Before the fading, we measure a de-reddened flux density of $\gtrsim$100 mJy in the mid-infrared, making MAXI J1535-571 one of the brightest mid-infrared BHBs known so far. A significant softening of the X-ray spectrum is evident contemporaneous with the infrared fade. We interpret it as due to the suppression of the jet emission, similar to the accretion-ejection coupling seen in other BHBs. However, MAXI J1535-571 did not transition smoothly to the soft state, instead showing X-ray hardness deviations, associated with infrared flaring. We also present the first mid-IR variability study of a BHB on minute timescales, with a fractional rms variability of the light curves of $\sim 15-22 \%$, which is similar to that expected from the internal shock jet model, and much higher than the optical fractional rms ($\lesssim 7 \%$). These results represent an excellent case of multi-wavelength jet spectral-timing and demonstrate how rich, multi-wavelength time-resolved data of X-ray binaries over accretion state transitions can help refining models of the disk-jet connection and jet launching in these systems.

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.

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 present Atacama Large Millimetre/Sub-Millimetre Array (ALMA) observations of IRAS 19132+1035, a candidate jet-ISM interaction zone near the black hole X-ray binary (BHXB) GRS 1915+105. With these ALMA observations (combining data from the 12 m array and the Atacama Compact Array), we map the molecular line emission across the IRAS 19132+1035 region. We detect emission from the $^{12}$CO [$J=2-1$], $^{13}$CO [$\nu=0$, $J=2-1$], C$^{18}$O [$J=2-1$], ${\rm H}_{2}{\rm CO}$ [$J=3_{0,3}-2_{0,2}$], ${\rm H}_{2}{\rm CO}$ [$J=3_{2,2}-2_{2,1}$], ${\rm H}_{2}{\rm CO}$ [$J=3_{2,1}-2_{2,0}$], SiO [$\nu=0$, $J=5-4$], CH$_3$OH [$J=4_{2,2}-3_{1,2}$], and CS [$\nu=0$, $J=5-4$] transitions. Given the morphological, spectral, and kinematic properties of this molecular emission, we present several lines of evidence that support the presence of a jet-ISM interaction at this site, including a jet-blown cavity in the molecular gas. This compelling new evidence identifies this site as a jet-ISM interaction zone, making GRS 1915$+$105 the third Galactic BHXB with at least one conclusive jet-ISM interaction zone. However, we find that this interaction occurs on much smaller scales than was postulated by previous work, where the BHXB jet does not appear to be dominantly powering the entire IRAS 19132+1035 region. Using estimates of the ISM conditions in the region, we utilize the detected cavity as a calorimeter to estimate the time-averaged power carried in the GRS 1915+105 jets of $(8.4^{+7.7}_{-8.1})\times10^{32}\,{\rm erg\,s}^{-1}$. Overall, our analysis demonstrates that molecular lines are excellent diagnostic tools to identify and probe jet-ISM interaction zones near Galactic BHXBs.

This review arose from the European Radio Astronomy Technical Forum (ERATec) meeting held in Firenze, October 2015, and aims to highlight the breadth and depth of the high-impact science that will be aided and assisted by the use of simultaneous mm-wavelength receivers. Recent results and opportunities are presented and discussed from the fields of: continuum VLBI (observations of weak sources, astrometry, observations of AGN cores in spectral index and Faraday rotation), spectral line VLBI (observations of evolved stars and massive star-forming regions) and time domain observations of the flux variations arising in the compact jets of X-ray binaries. Our survey brings together a large range of important science applications, which will greatly benefit from simultaneous observing at mm-wavelengths. Such facilities are essential to allow these applications to become more efficient, more sensitive and more scientifically robust. In some cases without simultaneous receivers the science goals are simply unachievable. Similar benefits would exist in many other high frequency astronomical fields of research.

We present simultaneous radio through sub-mm observations of the black hole X-ray binary (BHXB) V404 Cygni during the most active phase of its June 2015 outburst. Our $4$ hour long set of overlapping observations with the Very Large Array, the Sub-millimeter Array, and the James Clerk Maxwell Telescope (SCUBA-2), covers 8 different frequency bands (including the first detection of a BHXB jet at $666 \,{\rm GHz}/450\mu m$), providing an unprecedented multi-frequency view of the extraordinary flaring activity seen during this period of the outburst. In particular, we detect multiple rapidly evolving flares, which reach Jy-level fluxes across all of our frequency bands. With this rich data set we performed detailed MCMC modeling of the repeated flaring events. Our custom model adapts the van der Laan synchrotron bubble model to include twin bi-polar ejections, propagating away from the black hole at bulk relativistic velocities, along a jet axis that is inclined to the line of sight. The emission predicted by our model accounts for projection effects, relativistic beaming, and the geometric time delay between the approaching and receding ejecta in each ejection event. We find that a total of 8 bi-polar, discrete jet ejection events can reproduce the emission that we observe in all of our frequency bands remarkably well. With our best fit model, we provide detailed probes of jet speed, structure, energetics, and geometry. Our analysis demonstrates the paramount importance of the mm/sub-mm bands, which offer a unique, more detailed view of the jet than can be provided by radio frequencies alone.

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.

XTE J1908+094 is an X-ray transient black hole candidate in the Galactic plane that was observed in outburst in 2002 and 2013. Here we present multi-frequency radio and X-ray data, including radio polarimetry, spanning the entire period of the 2013 outburst. We find that the X-ray behaviour of XTE J1908+094 traces the standard black hole hardness-intensity path, evolving from a hard state, through a soft state, before returning to a hard state and quiescence. Its radio behaviour is typical of a compact jet that becomes quenched before discrete ejecta are launched during the late stages of X-ray softening. The radio and X-ray fluxes, as well as the light curve morphologies, are consistent with those observed during the 2002 outburst of this source. The polarisation angle during the rise of the outburst infers a jet orientation in agreement with resolved observations but also displays a gradual drift, which we associate with observed changes in the structure of the discrete ejecta. We also observe an unexpected 90deg rotation of the polarisation angle associated with a second component.

We present the results of our observations of the early stages of the 2012--2013 outburst of the transient black hole X-ray binary (BHXRB), Swift J1745$-$26, with the VLA, SMA, and JCMT (SCUBA--2). Our data mark the first multiple-band mm & sub-mm observations of a BHXRB. During our observations the system was in the hard accretion state producing a steady, compact jet. The unique combination of radio and mm/sub-mm data allows us to directly measure the spectral indices in and between the radio and mm/sub-mm regimes, including the first mm/sub-mm spectral index measured for a BHXRB. Spectral fitting revealed that both the mm (230 GHz) and sub-mm (350 GHz) measurements are consistent with extrapolations of an inverted power-law from contemporaneous radio data (1--30 GHz). This indicates that, as standard jet models predict, a power-law extending up to mm/sub-mm frequencies can adequately describe the spectrum, and suggests that the mechanism driving spectral inversion could be responsible for the high mm/sub-mm fluxes (compared to radio fluxes) observed in outbursting BHXRBs. While this power-law is also consistent with contemporaneous optical data, the optical data could arise from either jet emission with a jet spectral break frequency of $\nu_{{\rm break}}\gtrsim1\times10^{14}\,{\rm Hz}$ or the combination of jet emission with a lower jet spectral break frequency of $\nu_{{\rm break}}\gtrsim2\times10^{11}\,{\rm Hz}$ and accretion disc emission. Our analysis solidifies the importance of the mm/sub-mm regime in bridging the crucial gap between radio and IR frequencies in the jet spectrum, and justifies the need to explore this regime further.

MAXI J1836-194 is a Galactic black hole candidate X-ray binary that was discovered in 2011 when it went into outburst. In this paper, we present the full radio monitoring of this system during its `failed' outburst, in which the source did not complete a full set of state changes, only transitioning as far as the hard intermediate state. Observations with the Karl G. Jansky Very Large Array (VLA) and Australia Telescope Compact Array (ATCA) show that the jet properties changed significantly during the outburst. The VLA observations detected linearly polarised emission at a level of ~1% early in the outburst, increasing to ~3% as the outburst peaked. High-resolution images with the Very Long Baseline Array (VLBA) show a ~15 mas jet along the position angle $-21 \pm 2^\circ$, in agreement with the electric vector position angle found from our polarisation results ($-21 \pm 4^\circ$), implying that the magnetic field is perpendicular to the jet. Astrometric observations suggest that the system required an asymmetric natal kick to explain its observed space velocity. Comparing quasi-simultaneous X-ray monitoring with the 5 GHz VLA observations from the 2011 outburst shows an unusually steep hard-state radio/X-ray correlation of $L_{\rm R} \propto L_{\rm X}^{1.8\pm0.2}$, where $L_{\rm R}$ and $L_{\rm X}$ denote the radio and X-ray luminosities, respectively. With ATCA and Swift monitoring of the source during a period of re-brightening in 2012, we show that the system lay on the same steep correlation. Due to the low inclination of this system, we then investigate the possibility that the observed correlation may have been steepened by variable Doppler boosting.