We present high cadence multi-frequency radio observations of the long Gamma-Ray Burst (GRB) 190829A, which was detected at photon energies above 100 GeV by the High Energy Stereoscopic System (H.E.S.S.). Observations with the Meer Karoo Array Telescope (MeerKAT, 1.3 GHz), and Arcminute Microkelvin Imager - Large Array (AMI-LA, 15.5 GHz) began one day post-burst and lasted nearly 200 days. We used complementary data from Swift X-Ray Telescope (XRT), which ran to 100 days post-burst. We detected a likely forward shock component with both MeerKAT and XRT up to over 100 days post-burst. Conversely, the AMI-LA light curve appears to be dominated by reverse shock emission until around 70 days post-burst when the afterglow flux drops below the level of the host galaxy. We also present previously unpublished observations of the other H.E.S.S.-detected GRB, GRB 180720B from AMI-LA, which shows likely forward shock emission that fades in less than 10 days. We present a comparison between the radio emission from the three GRBs with detected very high energy (VHE) gamma-ray emission and a sensitivity-limited radio afterglow sample. GRB 190829A has the lowest isotropic radio luminosity of any GRB in our sample, but the distribution of luminosities is otherwise consistent, as expected, with the VHE GRBs being drawn from the same parent distribution as the other radio-detected long GRBs.

We have found an unexpected class of astronomical objects which have not previously been reported, in the Evolutionary Map of the Universe Pilot survey, using the Australian Square Kilometre Array Pathfinder telescope. The objects appear in radio images as circular edge-brightened discs about one arcmin diameter, and do not seem to correspond to any known type of object. We speculate that they may represent a spherical shock wave from an extra-galactic transient event, or the outflow, or a remnant, from a radio galaxy viewed end-on.

Long-duration gamma-ray bursts (GRBs) originate from ultra-relativistic jets launched from the collapsing cores of dying massive stars. They are characterised by an initial phase of bright and highly variable radiation in the keV-MeV band that is likely produced within the jet and lasts from milliseconds to minutes, known as the prompt emission. Subsequently, the interaction of the jet with the external medium generates external shock waves, responsible for the afterglow emission, which lasts from days to months, and occurs over a broad energy range, from the radio to the GeV bands. The afterglow emission is generally well explained as synchrotron radiation by electrons accelerated at the external shock. Recently, an intense, long-lasting emission between 0.2 and 1 TeV was observed from the GRB 190114C. Here we present the results of our multi-frequency observational campaign of GRB~190114C, and study the evolution in time of the GRB emission across 17 orders of magnitude in energy, from $5\times10^{-6}$ up to $10^{12}$\,eV. We find that the broadband spectral energy distribution is double-peaked, with the TeV emission constituting a distinct spectral component that has power comparable to the synchrotron component. This component is associated with the afterglow, and is satisfactorily explained by inverse Compton upscattering of synchrotron photons by high-energy electrons. We find that the conditions required to account for the observed TeV component are not atypical, supporting the possibility that inverse Compton emission is commonly produced in GRBs.

We present the discovery with Keck of the extremely infrared (IR) luminous transient AT 2017gbl, coincident with the Northern nucleus of the luminous infrared galaxy (LIRG) IRAS 23436+5257. Our extensive multi-wavelength follow-up spans ~900 days, including photometry and spectroscopy in the optical and IR, and (very long baseline interferometry) radio and X-ray observations. Radiative transfer modelling of the host galaxy spectral energy distribution shows the presence of a hitherto undetected dust obscured active galactic nucleus (AGN). The optical and near-IR spectra show broad 2000 km/s hydrogen, He I and O I emission features that decrease in flux over time. Radio imaging shows a fast evolving compact source of synchrotron emission spatially coincident with AT 2017gbl. We infer a lower limit for the radiated energy of 7.3 x 10^50 erg from the IR photometry. An extremely energetic supernova would satisfy this budget, but is ruled out by the radio counterpart evolution. Instead, we propose AT 2017gbl is related to an accretion event by the central supermassive black hole, where the spectral signatures originate in the AGN broad line region and the IR photometry is consistent with re-radiation by polar dust. Given the fast evolution of AT 2017gbl, we deem a tidal disruption event (TDE) of a star a more plausible scenario than a dramatic change in the AGN accretion rate. This makes AT 2017gbl the third TDE candidate to be hosted by a LIRG, in contrast to the so far considered TDE population discovered at optical wavelengths and hosted preferably by post-starburst galaxies.

We present radio and optical afterglow observations of the TeV-bright long Gamma Ray Burst (GRB) 190114C, which was detected by the MAGIC telescope at a redshift of z=0.425. Our observations with ALMA, ATCA, and uGMRT were obtained by our low frequency observing campaign and range from ~1 to ~140 days after the burst and the optical observations were done with the 0.7-m GROWTH-India telescope up to ~25 days after the burst. Long term radio/mm observations reveal the complex nature of the afterglow, which does not conform to the predictions of the standard afterglow model. We find that the microphysical parameters of the external forward shock, representing the share of shock-created energy in the non-thermal electron population and magnetic field, are evolving with time. The kinetic energy in the blast-wave is almost an order of magnitude higher than that measured in the prompt emission.

The Murchison Widefield Array (MWA) is an open access telescope dedicated to studying the low frequency (80$-$300 MHz) southern sky. Since beginning operations in mid 2013, the MWA has opened a new observational window in the southern hemisphere enabling many science areas. The driving science objectives of the original design were to observe 21\,cm radiation from the Epoch of Reionisation (EoR), explore the radio time domain, perform Galactic and extragalactic surveys, and monitor solar, heliospheric, and ionospheric phenomena. All together 60$+$ programs recorded 20,000 hours producing 146 papers to date. In 2016 the telescope underwent a major upgrade resulting in alternating compact and extended configurations. Other upgrades, including digital back-ends and a rapid-response triggering system, have been developed since the original array was commissioned. In this paper we review the major results from the prior operation of the MWA, and then discuss the new science paths enabled by the improved capabilities. We group these science opportunities by the four original science themes, but also include ideas for directions outside these categories.

The Murchison Widefield Array (MWA) is an electronically steered low frequency ($<300$\,MHz) radio interferometer, with a `slew' time less than 8seconds. Low frequency ($\sim 100$\,MHz) radio telescopes are ideally suited for rapid-response follow-up of transients due to their large field of view, the inverted spectrum of coherent emission, and the fact that the dispersion delay between a 1GHz and 100MHz pulse is on the order of $1-10$\,min for dispersion measures of $100-2000$\,pc/cm$^3$. The MWA has previously been used to provide fast follow up for transient events including gamma-ray bursts, fast radio bursts, and gravitational waves, using systems that respond to gamma-ray coordinates network (GCN) packet-based notifications. We describe a system for automatically triggering MWA observations of such events, based on VOEvent triggers, which is more flexible, capable, and accurate than previous systems. The system can respond to external multi-messenger triggers, which makes it well-suited to searching for prompt coherent radio emission from gamma-ray bursts, the study of fast radio bursts and gravitational waves, single pulse studies of pulsars, and rapid follow-up of high-energy superflares from flare stars. The new triggering system has the capability to trigger observations in both the regular correlator mode (limited to $\geq 0.5$\,s integrations) or using the Voltage Capture System (VCS, $0.1$\,ms integration) of the MWA, and represents a new mode of operation for the MWA. The upgraded standard correlator triggering capability has been in use since MWA observing semester 2018B (July-Dec 2018), and the VCS and buffered mode triggers will become available for observing in a future semester.

We report the discovery of a radio transient in the spiral galaxy M81. The transient was detected in early 2015 as part of a two-year survey of M81 made up of 12 epochs using the Karl G. Jansky Very Large Array. While undetected on 2014 September 12, the source was first detected on 2015 January 2, from which point it remained visible at an approximately constant luminosity of $L_{R} = 1.5 \pm 0.1 \times 10^{24}$erg s$^{-1}$ Hz$^{-1}$ at the observing frequency of 6 GHz for at least 2 months. Assuming this is a synchrotron event with a rise-time between 2.6 and 112 days, the peak luminosity (at equipartition) corresponds to a minimum energy of $10^{44} \lesssim E_{min} \lesssim 10^{46}$erg and jet power of $P_{min} \sim 10^{39}$erg s$^{-1}$, which are higher than most known X-ray binaries. Given its longevity, lack of short-term radio variability, and the absence of any multi-wavelength counterpart (X-ray luminosity $L_{X} \lesssim 10^{36}$erg s$^{-1}$), it does not behave like known Galactic or extragalactic X-ray binaries. The M81 transient radio properties more closely resemble the unidentified radio transient 43.78+59.3 discovered in M82, which has been suggested to be a radio nebula associated with an accreting source similar to SS 433. One possibility is that both the new M81 transient and the M82 transient may be the birth of a short-lived radio bubble associated with a discrete accretion event similar to those observed from the ULX Holmberg II X-1. However, it is not possible to rule out other identifications including long-term supernova shockwave interactions with the surrounding medium from a faint supernova or a background active galaxy.

We examine how fast radio burst (FRB)-like signals predicted to be generated during the merger of a binary neutron star (BNS) may be detected in low-frequency radio observations triggered by the aLIGO/Virgo gravitational wave detectors. The rapidity, directional accuracy, and sensitivity of follow-up observations with the Murchison Widefield Array (MWA) are considered. We show that with current methodology, the rapidity criteria fails for triggered MWA observations above 136 MHz for BNS mergers within the aLIGO/Virgo horizon, for which little dispersive delay is expected. A calculation of the expected reduction in response time by triggering on `negative latency' alerts from aLIGO/Virgo observations of gravitational waves generated by the BNS inspiral is presented. This allows for observations up to 300 MHz where the radio signal is expected to be stronger. To compensate for the poor positional accuracy expected from these alerts, we propose a new MWA observational mode that is capable of viewing one quarter of the sky. We show the sensitivity of this mode is sufficient to detect an FRB-like burst from an event similar to GW170817 if it occurred during the ongoing aLIGO/Virgo third science run (O3).

The presence and detectability of coherent radio emission from compact binary mergers (containing at least one neutron star) remains poorly constrained due to large uncertainties in the models. These compact binary mergers may initially be detected as Short Gamma-ray Bursts (SGRBs) or via their gravitational wave emission. Several radio facilities have developed rapid response modes enabling them to trigger on these events and search for this emission. For this paper, we constrain this coherent radio emission using the deepest available constraints for GRB 150424A, which were obtained via a triggered observation with the Murchison Widefield Array. We then expand this analysis to determine the properties of magnetar merger remnants that may be formed via a general population of binary neutron star mergers. Our results demonstrate that many of the potential coherent emission mechanisms that have been proposed for such events can be detected or very tightly constrained by the complementary strategies used by the current generation of low-frequency radio telescopes.

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.

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.

We present the results of an intensive multi-epoch radio frequency campaign on the energetic and nearby GRB 171010A with the Karl G. Janksy Very Large Array and Arcminute Microkelvin Imager Large Array. We began observing GRB 171010A a day after its initial detection, and were able to monitor the temporal and spectral evolution of the source over the following weeks. The spectra and their evolution are compared to the canonical theories for broadband GRB afterglows, with which we find a general agreement. There are, however, a number of features that are challenging to explain with a simple forward shock model, and we discuss possible reasons for these discrepancies. This includes the consideration of the existence of a reverse shock component, potential microphysical parameter evolution and the effect of scintillation.

We investigate the afterglow of GRB 140713A, a gamma-ray burst (GRB) that was detected and relatively well-sampled at X-ray and radio wavelengths, but was not present at optical and near-infrared wavelengths, despite searches to deep limits. We present the emission spectrum of the likely host galaxy at $z = 0.935$ ruling out a high-redshift explanation for the absence of the optical flux detection. Modelling the GRB multi-wavelength afterglow using the radiative transfer hydrodynamics code BOXFIT provides constraints on physical parameters of the GRB jet and its environment, for instance a relatively wide jet opening angle and an electron energy distribution slope $p$ below 2. Most importantly, the model predicts an optical flux about two orders of magnitude above the observed limits. We calculated that the required host extinction to explain the observed limits in the $r$, $i$ and $z$ bands was $A^{\rm host}_{V} > 3.2$ mag, equivalent to $E(B-V)^{\rm host} > 1.0$ mag. From the X-ray absorption we derive that the GRB host extinction is $A^{\rm host}_{\rm V} = 11.6^{+7.5}_{-5.3}$ mag, equivalent to $E(B-V)^{\rm host} = 3.7^{+2.4}_{-1.7}$ mag, which is consistent with the extinction required from our BOXFIT derived fluxes. We conclude that the origin of the optical darkness is a high level of extinction in the line of sight to the GRB, most likely within the GRB host galaxy.

The Space Technology 7 Disturbance Reduction System (ST7-DRS) is a NASA technology demonstration payload that operated from January 2016 through July of 2017 on the European Space Agency's LISA Pathfinder spacecraft. The joint goal of the NASA and ESA missions was to validate key technologies for a future space-based gravitational wave observatory targeting the source-rich milliHertz band. The two primary components of ST7-DRS are a micropropulsion system based on colloidal micro-Newton thrusters (CMNTs) and a control system that simultaneously controls the attitude and position of the spacecraft and the two free-flying test masses (TMs). This paper presents our main experimental results and summarizes the overall the performance of the CMNTs and control laws. We find that the CMNT performance to be consistent with pre-flight predictions, with a measured system thrust noise on the order of $100\,\textrm{nN}/\sqrt{\textrm{Hz}}$ in the $1\,\textrm{mHz}\leq f \leq 30\,\textrm{mHz}$ band. The control system maintained the TM-spacecraft separation with an RMS error of less than 2$\,$nm and a noise spectral density of less than $3\,\textrm{nm}/\sqrt{\textrm{Hz}}$ in the same band. Thruster calibration measurements yield thrust values consistent with the performance model and ground-based thrust-stand measurements, to within a few percent. We also report a differential acceleration noise between the two test masses with a spectral density of roughly $3\,\textrm{fm}/\textrm{s}^2/\sqrt{\textrm{Hz}}$ in the $1\,\textrm{mHz}\leq f \leq 30\,\textrm{mHz}$ band, slightly less than twice as large as the best performance reported with the baseline LISA Pathfinder configuration and below the current requirements for the Laser Interferometer Space Antenna (LISA) mission.

We present the discovery of a powerful candidate microquasar in NGC 300, associated with the S 10 optical nebula (previously classified as a supernova remnant). Chandra images show four discrete X-ray knots aligned in the plane of the sky over a length of $\approx$150 pc. The X-ray emission from the knots is well fitted with a thermal plasma model at a temperature of $\approx$0.6 keV and a combined 0.3-8 keV luminosity of $\approx$10$^{37}$ erg s$^{-1}$. The X-ray core, if present at all, does not stand out above the thermal emission of the knots: this suggests that the accreting compact object is either currently in a dim state or occulted from our view. We interpret the emission from the knots as the result of shocks from the interaction of a jet with the interstellar medium (possibly over different epochs of enhanced activity). Cooler shock-heated gas is likely the origin also of the optical bubble and lobes near the X-ray structure, detected in images from the Hubble Space Telescope and the Very Large Telescope. In the radio bands, we observed the region with the Australia Telescope Compact Array, and discovered an elongated radio nebula (about 170 $\times$ 55 pc in size) with its major axis aligned with the chain of Chandra sources. The radio nebula has an integrated 5.5 GHz radio luminosity of $\approx$10$^{34}$ erg s$^{-1}$ for a distance of 1.88 Mpc. The morphology, size and luminosity of the extended X-ray, optical and radio structure suggest that NGC 300-S 10 belongs to the same class of powerful ($P_{\rm jet} > 10^{39}$ erg s$^{-1}$) microquasars as SS 433, Ho II X-1 and NGC 7793-S 26.

The technique of Very Long Baseline Interferometry (VLBI) can provide accurate localization and unique physical information about radio transients. However, it is still underutilized due to the inherent difficulties of VLBI data analysis and practical difficulties of organizing observations on short notice. We present a brief overview of the currently available VLBI arrays and observing strategies used to study long- and short-duration radio transients.

We present the second installment of the Massive Star-forming Regions (MSFRs) Omnibus X-ray Catalog (MOXC2), a compilation of X-ray point sources detected in Chandra/ACIS observations of 16 Galactic MSFRs and surrounding fields. MOXC2 includes 13 ACIS mosaics, three containing a pair of unrelated MSFRs at different distances, with a total catalog of 18,396 point sources. The MSFRs sampled range over distances of 1.3 kpc to 6 kpc and populations varying from single massive protostars to the most massive Young Massive Cluster known in the Galaxy. By carefully detecting and removing X-ray point sources down to the faintest statistically-significant limit, we facilitate the study of the remaining unresolved X-ray emission. Through comparison with mid-infrared images that trace photon-dominated regions and ionization fronts, we see that the unresolved X-ray emission is due primarily to hot plasmas threading these MSFRs, the result of feedback from the winds and supernovae of massive stars. The 16 MSFRs studied in MOXC2 more than double the MOXC1 sample, broadening the parameter space of ACIS MSFR explorations and expanding Chandra's substantial contribution to contemporary star formation science.

We present LOFAR high-band data over the frequency range 115-189 MHz for the X-ray binary SS 433, obtained in an observing campaign from 2013 February - 2014 May. Our results include a deep, wide-field map, allowing a detailed view of the surrounding supernova remnant W 50 at low radio frequencies, as well as a light curve for SS 433 determined from shorter monitoring runs. The complex morphology of W 50 is in excellent agreement with previously published higher-frequency maps; we find additional evidence for a spectral turnover in the eastern wing, potentially due to foreground free-free absorption. Furthermore, SS 433 is tentatively variable at 150 MHz, with both a debiased modulation index of 11 per cent and a $\chi^2$ probability of a flat light curve of $8.2 \times 10^{-3}$. By comparing the LOFAR flux densities with contemporaneous observations carried out at 4800 MHz with the RATAN-600 telescope, we suggest that an observed $\sim$0.5-1 Jy rise in the 150-MHz flux density may correspond to sustained flaring activity over a period of approximately six months at 4800 MHz. However, the increase is too large to be explained with a standard synchrotron bubble model. We also detect a wealth of structure along the nearby Galactic plane, including the most complete detection to date of the radio shell of the candidate supernova remnant G 38.7-1.4. This further demonstrates the potential of supernova remnant studies with the current generation of low-frequency radio telescopes.

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 the discovery and interpretation of ionized nebulae around two ultraluminous X-ray sources in M 51; both sources share the rare property of showing X-ray eclipses by their companion stars, and are therefore prime targets for follow-up studies. Using archival Hubble Space Telescope images, we found an elongated, 100-pc-long emission-line structure associated with one X-ray source (CXOM51 J132940.0+471237; ULX-1 for simplicity), and a more circular, ionized nebula at the location of the second source (CXOM51 J132939.5+471244; ULX-2 for simplicity). We observed both nebulae with the Large Binocular Telescope's Multi-Object Double Spectrograph. From our analysis of the optical spectra, we argue that the gas in the ULX-1 bubble is shock-ionized, consistent with the effect of a jet with a kinetic power of $\approx$2 $\times 10^{39}$ erg s$^{-1}$. Additional X-ray photo-ionization may also be present, to explain the strength of high-ionization lines such as He II $\lambda 4686$ and [Ne V] $\lambda 3426$. On the other hand, the emission lines from the ULX-2 bubble are typical for photoionization by normal O stars suggesting that the nebula is actually an H II region not physically related to the ULX but is simply a chance alignment. From archival Very Large Array data, we also detect spatially extended, steep-spectrum radio emission at the location of the ULX-1 bubble (consistent with its jet origin), but no radio counterpart for ULX-2 (consistent with the lack of shock-ionized gas around that source).

The discovery of the first electromagnetic counterpart to a gravitational wave signal has generated follow-up observations by over 50 facilities world-wide, ushering in the new era of multi-messenger astronomy. In this paper, we present follow-up observations of the gravitational wave event GW170817 and its electromagnetic counterpart SSS17a/DLT17ck (IAU label AT2017gfo) by 14 Australian telescopes and partner observatories as part of Australian-based and Australian-led research programs. We report early- to late-time multi-wavelength observations, including optical imaging and spectroscopy, mid-infrared imaging, radio imaging, and searches for fast radio bursts. Our optical spectra reveal that the transient source afterglow cooled from approximately 6400K to 2100K over a 7-day period and produced no significant optical emission lines. The spectral profiles, cooling rate, and photometric light curves are consistent with the expected outburst and subsequent processes of a binary neutron star merger. Star formation in the host galaxy probably ceased at least a Gyr ago, although there is evidence for a galaxy merger. Binary pulsars with short (100 Myr) decay times are therefore unlikely progenitors, but pulsars like PSR B1534+12 with its 2.7 Gyr coalescence time could produce such a merger. The displacement (about 2.2 kpc) of the binary star system from the centre of the main galaxy is not unusual for stars in the host galaxy or stars originating in the merging galaxy, and therefore any constraints on the kick velocity imparted to the progenitor are poor.

Whilst astronomy as a science is historically founded on observations at optical wavelengths, studying the Universe in other bands has yielded remarkable discoveries, from pulsars in the radio, signatures of the Big Bang at submm wavelengths, through to high energy emission from accreting, gravitationally-compact objects and the discovery of gamma-ray bursts. Unsurprisingly, the result of combining multiple wavebands leads to an enormous increase in diagnostic power, but powerful insights can be lost when the sources studied vary on timescales shorter than the temporal separation between observations in different bands. In July 2015, the workshop "Paving the way to simultaneous multi-wavelength astronomy" was held as a concerted effort to address this at the Lorentz Center, Leiden. It was attended by 50 astronomers from diverse fields as well as the directors and staff of observatories and spaced-based missions. This community white paper has been written with the goal of disseminating the findings of that workshop by providing a concise review of the field of multi-wavelength astronomy covering a wide range of important source classes, the problems associated with their study and the solutions we believe need to be implemented for the future of observational astronomy. We hope that this paper will both stimulate further discussion and raise overall awareness within the community of the issues faced in a developing, important field.

We present the Arcminute Microkelvin Imager (AMI) Large Array catalogue of 139 gamma-ray bursts (GRBs). AMI observes at a central frequency of 15.7 GHz and is equipped with a fully automated rapid-response mode, which enables the telescope to respond to high-energy transients detected by Swift. On receiving a transient alert, AMI can be on-target within two minutes, scheduling later start times if the source is below the horizon. Further AMI observations are manually scheduled for several days following the trigger. The AMI GRB programme probes the early-time (< 1 day) radio properties of GRBs, and has obtained some of the earliest radio detections (GRB 130427A at 0.36 and GRB 130907A at 0.51 days post-burst). As all Swift GRBs visible to AMI are observed, this catalogue provides the first representative sample of GRB radio properties, unbiased by multi-wavelength selection criteria. We report the detection of six GRB radio afterglows that were not previously detected by other radio telescopes, increasing the rate of radio detections by 50% over an 18-month period. The AMI catalogue implies a Swift GRB radio detection rate of >15%, down to ~0.2 mJy/beam. However, scaling this by the fraction of GRBs AMI would have detected in the Chandra & Frail sample (all radio-observed GRBs between 1997 - 2011), it is possible ~44 - 56% of Swift GRBs are radio-bright, down to ~0.1 - 0.15 mJy/beam. This increase from the Chandra & Frail rate (~30%) is likely due to the AMI rapid-response mode, which allows observations to begin while the reverse-shock is contributing to the radio afterglow.

We have detected four flares from UV Ceti at 154 MHz using the Murchison Widefield Array. The flares have flux densities between 10--65 mJy --- a factor of 100 fainter than most flares in the literature at these frequencies --- and are only detected in polarization. The circular polarized fractions are limited to $>27$% at 3$\sigma$ confidence and two of the flares exhibit polarity reversal. We suggest that these flares occur periodically on a time scale consistent with the rotational period of UV Ceti. During the brightest observed flare, we also detect significant linear polarization with polarization fraction $>18$%. Averaging the data in 6-minute, 10 MHz frequency bins we find that the flux density of these flares does not vary over the 30 MHz bandwidth of the Murchison Widefield Array, however we cannot rule out finer time-frequency structure. Using the measured flux densities for the flares, we estimate brightness temperatures between $(10^{13}-10^{14})\,$K, indicative of a coherent emission mechanism. The brightness temperature and polarization characteristics point to the electron cyclotron maser mechanism. We also calculate the flare rates given our four observed flares and compare them to flare rates for the set of M dwarf stars with known 100--200 MHz flares. Our measurement is the first for flares with intensities $<100$ mJy at 100-200 MHz.

We present a radio light curve of supernova (SN) 2014C taken with the Arcminute Microkelvin Imager (AMI) Large Array at 15.7 GHz. Optical observations presented by Milisavljevic et al. demonstrated that SN 2014C metamorphosed from a stripped-envelope Type Ib SN into a strongly interacting Type IIn SN within 1 year. The AMI light curve clearly shows two distinct radio peaks, the second being a factor of 4 times more luminous than the first peak. This double bump morphology indicates two distinct phases of mass-loss from the progenitor star with the transition between density regimes occurring at 100-200 days. This reinforces the interpretation that SN 2014C exploded in a low density region before encountering a dense Hydrogen-rich shell of circumstellar material that was likely ejected by the progenitor prior to the explosion. The AMI flux measurements of the first light curve bump are the only reported observations taken within ~50 to ~125 days post-explosion, before the blast-wave encountered the Hydrogen shell. Simplistic synchrotron self-absorption (SSA) and free-free absorption (FFA) modelling suggest that some physical properties of SN 2014C, such as the mass-loss rate, are consistent with the properties of other Type Ibc and IIn SNe. However, our single frequency data does not allow us to distinguish between these two models, which implies they are likely too simplistic to describe the complex environment surrounding this event. Lastly, we present the precise radio location of SN 2014C obtained with eMERLIN, which will be useful for future VLBI observations of the SN.

We report on radio and X-ray monitoring observations of the BHC Swift J1753.5-0127 taken over a ~10 year period. Presented are daily radio observations at 15 GHz with the AMI-LA and X-ray data from Swift XRT and BAT. Also presented is a deep 2hr JVLA observation taken in an unusually low-luminosity soft-state (with a low disk temperature). We show that although the source has remained relatively radio-quiet compared to XRBs with a similar X-ray luminosity in the hard-state, the power-law relationship scales as $\zeta=0.96\pm0.06$ i.e. slightly closer to what has been considered for radiatively inefficient accretion disks. We also place the most stringent limit to date on the radio-jet quenching in an XRB soft-state, showing the connection of the jet quenching to the X-ray power-law component; the radio flux in the soft-state was found to be $<21~\mu$Jy, which is a quenching factor of $\gtrsim25$

The tidal disruption of a star by a supermassive black hole leads to a short-lived thermal flare. Despite extensive searches, radio follow-up observations of known thermal stellar tidal disruption flares (TDFs) have not yet produced a conclusive detection. We present a detection of variable radio emission from a thermal TDF, which we interpret as originating from a newly-launched jet. The multi-wavelength properties of the source present a natural analogy with accretion state changes of stellar mass black holes, suggesting all TDFs could be accompanied by a jet. In the rest frame of the TDF, our radio observations are an order of magnitude more sensitive than nearly all previous upper limits, explaining how these jets, if common, could thus far have escaped detection.

The lack of unambiguous detections of atomic features in the X-ray spectra of ultraluminous X-ray sources (ULXs) has proven a hindrance in diagnosing the nature of the accretion flow. The possible association of spectral residuals at soft energies with atomic features seen in absorption and/or emission and potentially broadened by velocity dispersion could therefore hold the key to understanding much about these enigmatic sources. Here we show for the first time that such residuals are seen in several sources and appear extremely similar in shape, implying a common origin. Via simple arguments we assert that emission from extreme colliding winds, absorption in a shell of material associated with the ULX nebula and thermal plasma emission associated with star formation are all highly unlikely to provide an origin. Whilst CCD spectra lack the energy resolution necessary to directly determine the nature of the features (i.e. formed of a complex of narrow lines or intrinsically broad), studying the evolution of the residuals with underlying spectral shape allows for an important, indirect test for their origin. The ULX NGC 1313 X-1 provides the best opportunity to perform such a test due to the dynamic range in spectral hardness provided by archival observations. We show through highly simplified spectral modelling that the strength of the features (in either absorption or emission) appears to anti-correlate with spectral hardness, which would rule out an origin via reflection of a primary continuum and instead supports a picture of atomic transitions in a wind or nearby material associated with such an outflow.

In preparing the way for the Square Kilometre Array and its pathfinders, there is a pressing need to begin probing the transient sky in a fully robotic fashion using the current generation of radio telescopes. Effective exploitation of such surveys requires a largely automated data-reduction process. This paper introduces an end-to-end automated reduction pipeline, AMIsurvey, used for calibrating and imaging data from the Arcminute Microkelvin Imager Large Array. AMIsurvey makes use of several component libraries which have been packaged separately for open-source release. The most scientifically significant of these is chimenea, which implements a telescope-agnostic algorithm for automated imaging of pre-calibrated multi-epoch radio-synthesis data, of the sort typically acquired for transient surveys or follow-up. The algorithm aims to improve upon standard imaging pipelines by utilizing iterative RMS-estimation and automated source-detection to avoid so called `Clean-bias', and makes use of CASA subroutines for the underlying image-synthesis operations. At a lower level, AMIsurvey relies upon two libraries, drive-ami and drive-casa, built to allow use of mature radio-astronomy software packages from within Python scripts. While targeted at automated imaging, the drive-casa interface can also be used to automate interaction with any of the CASA subroutines from a generic Python process. Additionally, these packages may be of wider technical interest beyond radio-astronomy, since they demonstrate use of the Python library pexpect to emulate terminal interaction with an external process. This approach allows for rapid development of a Python interface to any legacy or externally-maintained pipeline which accepts command-line input, without requiring alterations to the original code.

The Fermi Gamma-ray Space Telescope has greatly expanded the number and energy window of observations of gamma-ray bursts (GRBs). However, the coarse localizations of tens to a hundred square degrees provided by the Fermi GRB Monitor instrument have posed a formidable obstacle to locating the bursts' host galaxies, measuring their redshifts, and tracking their panchromatic afterglows. We have built a target-of-opportunity mode for the intermediate Palomar Transient Factory in order to perform targeted searches for Fermi afterglows. Here, we present the results of one year of this program: 8 afterglow discoveries out of 35 searches. Two of the bursts with detected afterglows (GRBs 130702A and 140606B) were at low redshift (z=0.145 and 0.384 respectively) and had spectroscopically confirmed broad-line Type Ic supernovae. We present our broadband follow-up including spectroscopy as well as X-ray, UV, optical, millimeter, and radio observations. We study possible selection effects in the context of the total Fermi and Swift GRB samples. We identify one new outlier on the Amati relation. We find that two bursts are consistent with a mildly relativistic shock breaking out from the progenitor star, rather than the ultra-relativistic internal shock mechanism that powers standard cosmological bursts. Finally, in the context of the Zwicky Transient Facility, we discuss how we will continue to expand this effort to find optical counterparts of binary neutron star mergers that may soon be detected by Advanced LIGO and Virgo.

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.

In this short paper, we present constraints on the Galileon inflationary model from the CMB bispectrum. We employ a principal-component analysis of the independent degrees of freedom constrained by data and apply this to the WMAP 9-year data to constrain the free parameters of the model. A simple Bayesian comparison establishes that support for the Galileon model from bispectrum data is at best weak.

On 2014 April 23, the Swift satellite detected a gamma-ray superflare from the nearby star system DG CVn. This system comprises a M-dwarf binary with extreme properties: it is very young and at least one of the components is a very rapid rotator. The gamma-ray superflare is one of only a handful detected by Swift in a decade. As part of our AMI-LA Rapid Response Mode, ALARRM, we automatically slewed to this target, were taking data at 15 GHz within six minutes of the burst, and detected a bright (~100 mJy) radio flare. This is the earliest detection of bright, prompt, radio emission from a high energy transient ever made with a radio telescope, and is possibly the most luminous incoherent radio flare ever observed from a red dwarf star. An additional bright radio flare, peaking at around 90 mJy, occurred around one day later, and there may have been further events between 0.1-1 days when we had no radio coverage. The source subsequently returned to a quiescent level of 2-3 mJy on a timescale of about 4 days. Although radio emission is known to be associated with active stars, this is the first detection of large radio flares associated with a gamma ray superflare, and demonstrates both feasibility and scientific importance of rapid response modes on radio telescopes.

GRB 130427A was extremely bright as a result of occurring at low redshift whilst the energetics were more typical of high-redshift gamma-ray bursts (GRBs). We collected well-sampled light curves at 1.4 and 4.8 GHz of GRB 130427A with the Westerbork Synthesis Radio Telescope (WSRT); and we obtained its most accurate position with the European Very Long Baseline Interferometry Network (EVN). Our flux density measurements are combined with all the data available at radio, optical and X-ray frequencies to perform broadband modeling in the framework of a reverse-forward shock model and a two-component jet model, and we discuss the implications and limitations of both models. The low density inferred from the modeling implies that the GRB 130427A progenitor is either a very low-metallicity Wolf-Rayet star, or a rapidly rotating, low-metallicity O star. We also find that the fraction of the energy in electrons is evolving over time, and that the fraction of electrons participating in a relativistic power-law energy distribution is less than 15%. We observed intraday variability during the earliest WSRT observations, and the source sizes inferred from our modeling are consistent with this variability being due to interstellar scintillation effects. Finally, we present and discuss our limits on the linear and circular polarization, which are among the deepest limits of GRB radio polarization to date.

We use WMAP 9-year bispectrum data to constrain the free parameters of an 'effective field theory' describing fluctuations in single-field inflation. The Lagrangian of the theory contains a finite number of operators associated with unknown mass scales. Each operator produces a fixed bispectrum shape, which we decompose into partial waves in order to construct a likelihood function. Based on this likelihood we are able to constrain four linearly independent combinations of the mass scales. As an example of our framework we specialize our results to the case of 'Dirac-Born-Infeld' and 'ghost' inflation and obtain the posterior probability for each model, which in Bayesian schemes is a useful tool for model comparison. Our results suggest that DBI-like models with two or more free parameters are disfavoured by the data by comparison with single parameter models in the same class.

We present one of the best sampled early time light curves of a gamma-ray burst (GRB) at radio wavelengths. Using the Arcminute Mircrokelvin Imager (AMI) we observed GRB 130427A at the central frequency of 15.7 GHz between 0.36 and 59.32 days post-burst. These results yield one of the earliest radio detections of a GRB and demonstrate a clear rise in flux less than one day after the gamma-ray trigger followed by a rapid decline. This early time radio emission probably originates in the GRB reverse shock so our AMI light curve reveals the first ever confirmed detection of a reverse shock peak in the radio domain. At later times (about 3.2 days post-burst) the rate of decline decreases, indicating that the forward shock component has begun to dominate the light-curve. Comparisons of the AMI light curve with modelling conducted by Perley et al. show that the most likely explanation of the early time 15.7 GHz peak is caused by the self-absorption turn-over frequency, rather than the peak frequency, of the reverse shock moving through the observing bands.

We present the Chasing the Identification of ASCA Galactic Objects (ChIcAGO) survey, which is designed to identify the unknown X-ray sources discovered during the ASCA Galactic Plane Survey (AGPS). Little is known about most of the AGPS sources, especially those that emit primarily in hard X-rays (2-10 keV) within the F_x ~ 10^-13 to 10^-11 erg cm^-2 s^-1 X-ray flux range. In ChIcAGO, the subarcsecond localization capabilities of Chandra have been combined with a detailed multi-wavelength follow-up program, with the ultimate goal of classifying the >100 unidentified sources in the AGPS. Overall to date, 93 unidentified AGPS sources have been observed with Chandra as part of the ChIcAGO survey. A total of 253 X-ray point sources have been detected in these Chandra observations within 3' of the original ASCA positions. We have identified infrared and optical counterparts to the majority of these sources, using both new observations and catalogs from existing Galactic plane surveys. X-ray and infrared population statistics for the X-ray point sources detected in the Chandra observations reveal that the primary populations of Galactic plane X-ray sources that emit in the F_x ~ 10^-13 to 10^-11 erg cm^-2 s^-1 flux range are active stellar coronae, massive stars with strong stellar winds that are possibly in colliding-wind binaries, X-ray binaries, and magnetars. There is also a fifth population that is still unidentified but, based on its X-ray and infrared properties, likely comprise partly of Galactic sources and partly of active galactic nuclei.

We report the design and development of a self-contained multi-band receiver (MBR) system, intended for use with a single large aperture to facilitate sensitive & high time-resolution observations simultaneously in 10 discrete frequency bands sampling a wide spectral span (100-1500 MHz) in a nearly log-periodic fashion. The development of this system was primarily motivated by need for tomographic studies of pulsar polar emission regions. Although the system design is optimized for the primary goal, it is also suited for several other interesting astronomical investigations. The system consists of a dual-polarization multi-band feed (with discrete responses corresponding to the 10 bands pre-selected as relatively RFI-free), a common wide-band RF front-end, and independent back-end receiver chains for the 10 individual sub-bands. The raw voltage time-sequences corresponding to 16 MHz bandwidth each for the two linear polarization channels and the 10 bands, are recorded at the Nyquist rate simultaneously. We present the preliminary results from the tests and pulsar observations carried out with the Green Bank Telescope using this receiver. The system performance implied by these results, and possible improvements are also briefly discussed.

We use transport techniques to calculate the trispectrum produced in multiple-field inflationary models with canonical kinetic terms. Our method allows the time evolution of the local trispectrum parameters, tauNL and gNL, to be tracked throughout the inflationary phase. We illustrate our approach using examples. We give a simplified method to calculate the superhorizon part of the relation between field fluctuations on spatially flat hypersurfaces and the curvature perturbation on uniform density slices, and obtain its third-order part for the first time. We clarify how the 'backwards' formalism of Yokoyama et al. relates to our analysis and other recent work. We supply explicit formulae which enable each inflationary observable to be computed in any canonical model of interest, using a suitable first-order ODE solver.

We present multi-wavelength observations of the radio magnetar PSR J1622-4950 and its environment. Observations of PSR J1622-4950 with Chandra (in 2007 and 2009) and XMM (in 2011) show that the X-ray flux of PSR J1622-4950 has decreased by a factor of ~50 over 3.7 years, decaying exponentially with a characteristic time of 360 +/- 11 days. This behavior identifies PSR J1622-4950 as a possible addition to the small class of transient magnetars. The X-ray decay likely indicates that PSR J1622-4950 is recovering from an X-ray outburst that occurred earlier in 2007, before the 2007 Chandra observations. Observations with the Australia Telescope Compact Array show strong radio variability, including a possible radio flaring event at least one and a half years after the 2007 X-ray outburst that may be a direct result of this X-ray event. Radio observations with the Molonglo Observatory Synthesis Telescope reveal that PSR J1622-4950 is 8' southeast of a diffuse radio arc, G333.9+0.0, which appears non-thermal in nature and which could possibly be a previously undiscovered supernova remnant. If G333.9+0.0 is a supernova remnant then the estimates of its size and age, combined with the close proximity and reasonable implied velocity of PSR J1622-4950, suggests that these two objects could be physically associated.

We present X-ray, infrared, optical and radio observations of four previously unidentified Galactic plane X-ray sources, AX J163252-4746, AX J184738-0156, AX J144701-5919 and AX J144547-5931. Detection of each source with the Chandra X-ray Observatory has provided sub-arcsecond localizations, which we use to identify bright infrared counterparts to all four objects. Infrared and optical spectroscopy of these counterparts demonstrate that all four X-ray sources are extremely massive stars, with spectral classifications Ofpe/WN9 (AX J163252-4746), WN7 (AX J184738-0156 = WR121a), WN7-8h (AX J144701-5919) and OIf+ (AX J144547-5931). AX J163252-4746 and AX J184738-0156 are both luminous, hard, X-ray emitters with strong Fe XXV emission lines in their X-ray spectra at ~6.7 keV. The multi-wavelength properties of AX J163252-4746 and AX J184738-0156 are not consistent with isolated massive stars or accretion onto a compact companion; we conclude that their X-ray emission is most likely generated in a colliding-wind binary system. For both AX J144701-5919 and AX J144547-5931, the X-ray emission is an order of magnitude less luminous and with a softer spectrum. These properties are consistent with a colliding-wind binary interpretation for these two sources also, but other mechanisms for the generation of X-rays cannot be excluded. There are many other as yet unidentified X-ray sources in the Galactic plane, with X-ray properties similar to those seen for AX J163252-4746, AX J184738-0156, AX J144701-5919 and AX J144547-5931. This may indicate a substantial population of X-ray-emitting massive stars and colliding-wind binaries in the Milky Way.

We present detailed optical, near-infrared, and radio observations of the X-ray flash 050416a obtained with Palomar and Siding Springs Observatories as well as HST and the VLA, placing this event among the best-studied X-ray flashes to date. In addition, we present an optical spectrum from Keck LRIS from which we measure the redshift of the burst, z=0.6528. At this redshift the isotropic-equivalent prompt energy release was about 10^51 erg, and using a standard afterglow synchrotron model we find that the blastwave kinetic energy is a factor of 10 larger, E_K,iso ~ 10^52 erg. The lack of an observed jet break to t ~ 20 days indicates that the opening angle is larger than 7 deg and the total beaming-corrected relativistic energy is larger than 10^50 erg. We further show that the burst produced a strong radio flare at t ~ 40 days accompanied by an observed flattening in the X-ray band which we attribute to an abrupt circumburst density jump or an episode of energy injection (either from a refreshed shock or off-axis ejecta). Late-time observations with HST show evidence for an associated supernova with peak optical luminosity roughly comparable to that of SN 1998bw. Next, we show that the host galaxy of XRF 050416a is actively forming stars at a rate of at least 2 M_solar per year with a luminosity of L_B ~ 0.5L* and metallicity of Z ~ 0.2-0.8 Z_solar. Finally, we discuss the nature of XRF 050416a in the context of short-hard gamma-ray bursts and under the framework of off-axis and dirty fireball models for X-ray flashes.

We present optical, near-IR, and radio follow up of sixteen Swift bursts, including our discovery of nine afterglows and a redshift determination for three. These observations, supplemented by data from the literature, provide an afterglow recovery rate of 60% in the optical/near-IR, much higher than in previous missions (BeppoSAX, HETE-2, INTEGRAL, and IPN). The optical/near-IR afterglows of Swift events are on average 1.7 mag fainter at t=12 hr than those of previous missions. The X-ray afterglows are similarly fainter compared to those of pre-Swift bursts. In the radio the limiting factor is the VLA threshold and the detection rate for Swift bursts is similar to that for past missions. The redshift distribution of pre-Swift bursts peaked at z~1, whereas the five Swift bursts with measured redshifts are distributed evenly between 1.3 and 3.2. From these results we conclude that (i) the pre-Swift distributions were biased in favor of bright events and low redshift events, (ii) the higher sensitivity and accurate positions of Swift result in a better representation of the true burst redshift and brightness distributions (which are higher and dimmer, respectively), and (iii) as many as 1/3 of the bursts can be optically dark, as a result of a high redshift and/or dust extinction. We remark that the apparent lack of low redshift, low luminosity Swift bursts, and the lower event rate compared to pre-launch estimates (90 vs. 150 per year), are the result of a threshold that is similar to that of BATSE. In view of these inferences, afterglow observers may find it advisable to make significant changes in follow up strategies of Swift events. [abridged]

We propose a simple model in which the cosmological dark matter consists of particles whose mass increases with the scale factor of the universe. The particle mass is generated by the expectation value of a scalar field which does not have a stable vacuum state, but which is effectively stabilized by the rest energy of the ambient particles. As the universe expands, the density of particles decreases, leading to an increase in the vacuum expectation value of the scalar (and hence the mass of the particle). The energy density of the coupled system of variable-mass particles (``vamps'') redshifts more slowly than that of ordinary matter. Consequently, the age of the universe is larger than in conventional scenarios.

Fluctuations of scalar fields produced at the stage of preheating after inflation are so large that they can break supersymmetry much stronger than inflation itself. These fluctuations may lead to symmetry restoration along flat directions of the effective potential even in the theories where the usual high temperature corrections are exponentially suppressed. Our results show that nonthermal phase transitions after preheating may play a crucial role in the generation of the primordial baryon asymmetry by the Affleck-Dine mechanism. In particular, the baryon asymmetry may be generated at the very early stage of the evolution of the Universe, at the preheating era, and not when the Hubble parameter becomes of order the gravitino mass.