How are accretion discs affected by their outflows? To address this question for white dwarfs accreting from cool giants, we performed optical, radio, X-ray, and ultraviolet observations of the outflow-driving symbiotic star MWC 560 (=V694 Mon) during its 2016 optical high state. We tracked multi-wavelength changes that signalled an abrupt increase in outflow power at the initiation of a months-long outflow fast state, just as the optical flux peaked: (1) an abrupt doubling of Balmer absorption velocities; (2) the onset of a $20$ $\mu$Jy/month increase in radio flux; and (3) an order-of-magnitude increase in soft X-ray flux. Juxtaposing to prior X-ray observations and their coeval optical spectra, we infer that both high-velocity and low-velocity optical outflow components must be simultaneously present to yield a large soft X-ray flux, which may originate in shocks where these fast and slow absorbers collide. Our optical and ultraviolet spectra indicate that the broad absorption-line gas was fast, stable, and dense ($\gtrsim10^{6.5}$ cm$^{-3}$) throughout the 2016 outflow fast state, steadily feeding a lower-density ($\lesssim10^{5.5}$ cm$^{-3}$) region of radio-emitting gas. Persistent optical and ultraviolet flickering indicate that the accretion disc remained intact. The stability of these properties in 2016 contrasts to their instability during MWC 560's 1990 outburst, even though the disc reached a similar accretion rate. We propose that the self-regulatory effect of a steady fast outflow from the disc in 2016 prevented a catastrophic ejection of the inner disc. This behaviour in a symbiotic binary resembles disc/outflow relationships governing accretion state changes in X-ray binaries.

We obtained radio observations of the symbiotic binary and known recurrent nova T Coronae Borealis following a period of increased activity in the optical and X-ray bands. A comparison of our observations with those made prior to 2015 indicates that the system is in a state of higher emission in the radio as well. The spectral energy distributions are consistent with optically thick thermal bremsstrahlung emission from a photoionized source. Our observations indicate that the system was in a state of increased ionization in the companion wind, possibly driven by an increase in accretion rate, with the radio photosphere located well outside the binary system.

It has recently been discovered that some, if not all, classical novae emit GeV gamma rays during outburst, but the mechanisms involved in the production of the gamma rays are still not well understood. We present here a comprehensive multi-wavelength dataset---from radio to X-rays---for the most gamma-ray luminous classical nova to-date, V1324 Sco. Using this dataset, we show that V1324 Sco is a canonical dusty Fe-II type nova, with a maximum ejecta velocity of 2600 km s$^{-1}$ and an ejecta mass of few $\times 10^{-5}$ M$_{\odot}$. There is also evidence for complex shock interactions, including a double-peaked radio light curve which shows high brightness temperatures at early times. To explore why V1324~Sco was so gamma-ray luminous, we present a model of the nova ejecta featuring strong internal shocks, and find that higher gamma-ray luminosities result from higher ejecta velocities and/or mass-loss rates. Comparison of V1324~Sco with other gamma-ray detected novae does not show clear signatures of either, and we conclude that a larger sample of similarly well-observed novae is needed to understand the origin and variation of gamma rays in novae.

We present multi-wavelength observations of the unusual nova V1535 Sco throughout its outburst in 2015. Early radio observations were consistent with synchrotron emission, and early X-ray observations revealed the presence of high-energy (>1 keV) photons. These indicated that strong shocks were present during the first ~2 weeks of the nova's evolution. The radio spectral energy distribution was consistent with thermal emission from week 2 to week 6. Starting in week 7, the radio emission again showed evidence of synchrotron emission and there was an increase in X-ray emission, indicating a second shock event. The optical spectra show evidence for at least two separate outflows, with the faster outflow possibly having a bipolar morphology. The optical and near infrared light curves and the X-ray measurements of the hydrogen column density indicated that the companion star is likely a K giant.

Since the Fermi discovery of $\gamma$-rays from novae, one of the biggest questions in the field has been how novae generate such high-energy emission. Shocks must be a fundamental ingredient. Six months of radio observations of the 2012 nova V5589 Sgr with the VLA and 15 weeks of X-ray observations with Swift/XRT show that the radio emission consisted of: 1) a shock-powered, non-thermal flare; and 2) weak thermal emission from $10^{-5}$ M$_\odot$ of freely expanding, photoionized ejecta. Absorption features in the optical spectrum and the peak optical brightness suggest that V5589 Sgr lies 4 kpc away (3.2-4.6 kpc). The shock-powered flare dominated the radio light curve at low frequencies before day 100. The spectral evolution of the radio flare, its high radio brightness temperature, the presence of unusually hard ($kT_x > 33$ keV) X-rays, and the ratio of radio to X-ray flux near radio maximum all support the conclusions that the flare was shock-powered and non-thermal. Unlike most other novae with strong shock-powered radio emission, V5589 Sgr is not embedded in the wind of a red-giant companion. Based on the similar inclinations and optical line profiles of V5589 Sgr and V959 Mon, we propose that shocks in V5589 Sgr formed from collisions between a slow flow with an equatorial density enhancement and a subsequent faster flow. We speculate that the relatively high speed and low mass of the ejecta led to the unusual radio emission from V5589 Sgr, and perhaps also to the non-detection of $\gamma$-rays.

The importance of shocks in nova explosions has been highlighted by Fermi's discovery of \gamma-ray producing novae. Over three years of multi-band VLA radio observations of the 2010 nova V1723 Aql show that shocks between fast and slow flows within the ejecta led to the acceleration of particles and the production of synchrotron radiation. Soon after the start of the eruption, shocks in the ejecta produced an unexpected radio flare, resulting in a multi-peaked radio light curve. The emission eventually became consistent with an expanding thermal remnant with mass $2 \times 10^{-4} M_\odot$ and temperature $10^4$ K. However, during the first two months, the $\gtrsim 10^6$ K brightness temperature at low frequencies was too high to be due to thermal emission from the small amount of X-ray producing shock-heated gas. Radio imaging showed structures with velocities of 400 km s$^{-1}$ (d/6 kpc) in the plane of the sky, perpendicular to a more elongated 1500 km s$^{-1}$ (d/6 kpc) flow. The morpho-kinematic structure of the ejecta from V1723 Aql appears similar to nova V959 Mon, where collisions between a slow torus and a faster flow collimated the fast flow and gave rise to \gamma -ray producing shocks. Optical spectroscopy and X-ray observations of V1723 Aql during the radio flare are consistent with this picture. Our observations support the idea that shocks in novae occur when a fast flow collides with a slow collimating torus. Such shocks could be responsible for hard X-ray emission, \gamma -ray production, and double-peaked radio light curves from some classical novae.

Determining reliable distances to classical novae is a challenging but crucial step in deriving their ejected masses and explosion energetics. Here we combine radio expansion measurements from the Karl G. Jansky Very Large Array with velocities derived from optical spectra to estimate an expansion parallax for nova V959 Mon, the first nova discovered through its gamma-ray emission. We spatially resolve the nova at frequencies of 4.5-36.5 GHz in nine different imaging epochs. The first five epochs cover the expansion of the ejecta from 2012 October to 2013 January, while the final four epochs span 2014 February to 2014 May. These observations correspond to days 126 through 199 and days 615 through 703 after the first detection of the nova. The images clearly show a non-spherical ejecta geometry. Utilizing ejecta velocities derived from 3D modelling of optical spectroscopy, the radio expansion implies a distance between 0.9 +/- 0.2 and 2.2 +/- 0.4 kpc, with a most probable distance of 1.4 +/- 0.4 kpc. This distance implies a gamma-ray luminosity much less than the prototype gamma-ray-detected nova, V407 Cyg, possibly due to the lack of a red giant companion in the V959 Mon system. V959 Mon also has a much lower gamma-ray luminosity than other classical novae detected in gamma-rays to date, indicating a range of at least a factor of 10 in the gamma-ray luminosities for these explosions.

Classical novae are the most common astrophysical thermonuclear explosions, occurring on the surfaces of white dwarf stars accreting gas from companions in binary star systems. Novae typically expel ~10^(-4) solar masses of material at velocities exceeding 1,000 kilometres per second. However, the mechanism of mass ejection in novae is poorly understood, and could be dominated by the impulsive flash of thermonuclear energy, prolonged optically thick winds, or binary interaction with the nova envelope. Classical novae are now routinely detected in gigaelectronvolt gamma-ray wavelengths, suggesting that relativistic particles are accelerated by strong shocks in the ejecta. Here we report high-resolution radio imaging of the gamma-ray-emitting nova V959 Mon. We find that its ejecta were shaped by the motion of the binary system: some gas was expelled rapidly along the poles as a wind from the white dwarf, while denser material drifted out along the equatorial plane, propelled by orbital motion. At the interface between the equatorial and polar regions, we observe synchrotron emission indicative of shocks and relativistic particle acceleration, thereby pinpointing the location of gamma-ray production. Binary shaping of the nova ejecta and associated internal shocks are expected to be widespread among novae, explaining why many novae are gamma-ray emitters.

The recurrent nova T Pyx underwent its sixth historical outburst in 2011, and became the subject of an intensive multi-wavelength observational campaign. We analyze data from the Swift and Suzaku satellites to produce a detailed X-ray light curve augmented by epochs of spectral information. X-ray observations yield mostly non-detections in the first four months of outburst, but both a super-soft and hard X-ray component rise rapidly after Day 115. The super-soft X-ray component, attributable to the photosphere of the nuclear-burning white dwarf, is relatively cool (~45 eV) and implies that the white dwarf in T Pyx is significantly below the Chandrasekhar mass (~1 M_sun). The late turn-on time of the super-soft component yields a large nova ejecta mass (>~10^-5 M_sun), consistent with estimates at other wavelengths. The hard X-ray component is well fit by a ~1 keV thermal plasma, and is attributed to shocks internal to the 2011 nova ejecta. The presence of a strong oxygen line in this thermal plasma on Day 194 requires a significantly super-solar abundance of oxygen and implies that the ejecta are polluted by white dwarf material. The X-ray light curve can be explained by a dual-phase ejection, with a significant delay between the first and second ejection phases, and the second ejection finally released two months after outburst. A delayed ejection is consistent with optical and radio observations of T Pyx, but the physical mechanism producing such a delay remains a mystery.

Despite being the prototype of its class, T Pyx is arguably the most unusual and poorly understood recurrent nova. Here, we use radio observations from the Karl G. Jansky Very Large Array to trace the evolution of the ejecta over the course of the 2011 outburst of T Pyx. The radio emission is broadly consistent with thermal emission from the nova ejecta. However, the radio flux began rising surprisingly late in the outburst, indicating that the bulk of the radio-emitting material was either very cold, or expanding very slowly, for the first ~50 days of the outburst. Considering a plausible range of volume filling factors and geometries for the ejecta, we find that the high peak flux densities of the radio emission require a massive ejection of 1-30 x 10^{-5} solar masses. This ejecta mass is much higher than the values normally associated with recurrent novae, and is more consistent with a nova on a white dwarf well below the Chandrasekhar limit.

The radio light curves of novae rise and fall over the course of months to years, allowing for detailed observations of the evolution of the nova shell. However, the main parameter determined by radio models of nova explosions - the mass of the ejecta - often seems to exceed theoretical expectations by an order of magnitude. With the recent technological improvements on the Karl G. Jansky Very Large Array (VLA), new observations can test the assumptions upon which ejecta mass estimates are based. Early observations of the classical nova V1723 Aql showed an unexpectedly rapid rise in radio flux density and a distinct bump in the radio light curve on the rise to radio maximum, which is inconsistent with the simple model of spherical ejecta expelled in a single discrete event. This initial bump appears to indicate the presence of shocked material in the outer region of the ejected shell, with the emission from the shocks fading over time. We explore possible origins for this emission and its relation to the mass loss history of the nova. The evolution of the radio spectrum also reveals the density profile, the mass of the ejected shell, and other properties of the ejecta. These observations comprise one of the most complete, longterm set of multi-wavelength radio observations for any classical nova to date.

Novae, which are the sudden visual brightening triggered by runaway thermonuclear burning on the surface of an accreting white dwarf, are fairly common and bright events. Despite their astronomical significance as nearby laboratories for the study of nuclear burning and accretion phenomena, many aspects of these common stellar explosions are observationally not well-constrained and remain poorly understood. Radio observations, modeling and interpretation can potentially play a crucial role in addressing some of these puzzling issues. In this review on radio studies of novae, we focus on the possibility of testing and improving the nova models with radio observations, and present a current status report on the progress in both the observational front and theoretical developments. We specifically address the issues of accurate estimation of ejecta mass, multi-phase and complex ejection phenomena, and the effect of a dense environment around novae. With highlights of new observational results, we illustrate how radio observations can shed light on some of these long-standing puzzles.

We present multi-frequency radio observations of the 2010 nova event in the symbiotic binary V407 Cygni, obtained with the Karl G. Jansky Very Large Array and spanning 1-45 GHz and 17-770 days following discovery. This nova---the first ever detected in gamma rays---shows a radio light curve dominated by the wind of the Mira giant companion, rather than the nova ejecta themselves. The radio luminosity grew as the wind became increasingly ionized by the nova outburst, and faded as the wind was violently heated from within by the nova shock. This study marks the first time that this physical mechanism has been shown to dominate the radio light curve of an astrophysical transient. We do not observe a thermal signature from the nova ejecta or synchrotron emission from the shock, due to the fact that these components were hidden behind the absorbing screen of the Mira wind. We estimate a mass loss rate for the Mira wind of Mdot_w ~ 10^-6 M_sun/yr. We also present the only radio detection of V407 Cyg before the 2010 nova, gleaned from unpublished 1993 archival VLA data, which shows that the radio luminosity of the Mira wind varies by a factor of >~20 even in quiescence. Although V407 Cyg likely hosts a massive accreting white dwarf, making it a candidate progenitor system for a Type Ia supernova, the dense and radially continuous circumbinary material surrounding V407 Cyg is inconsistent with observational constraints on the environments of most Type Ia supernovae.