The Asteroid Terrestrial impact Last Alert System (ATLAS) system consists of two 0.5m Schmidt telescopes with cameras covering 29 square degrees at plate scale of 1.86 arcsec per pixel. Working in tandem, the telescopes routinely survey the whole sky visible from Hawaii (above $\delta > -50^{\circ}$) every two nights, exposing four times per night, typically reaching $o < 19$ magnitude per exposure when the moon is illuminated and $c < 19.5$ per exposure in dark skies. Construction is underway of two further units to be sited in Chile and South Africa which will result in an all-sky daily cadence from 2021. Initially designed for detecting potentially hazardous near earth objects, the ATLAS data enable a range of astrophysical time domain science. To extract transients from the data stream requires a computing system to process the data, assimilate detections in time and space and associate them with known astrophysical sources. Here we describe the hardware and software infrastructure to produce a stream of clean, real, astrophysical transients in real time. This involves machine learning and boosted decision tree algorithms to identify extragalactic and Galactic transients. Typically we detect 10-15 supernova candidates per night which we immediately announce publicly. The ATLAS discoveries not only enable rapid follow-up of interesting sources but will provide complete statistical samples within the local volume of 100 Mpc. A simple comparison of the detected supernova rate within 100 Mpc, with no corrections for completeness, is already significantly higher (factor 1.5 to 2) than the current accepted rates.

We present our findings based on pre- and post-explosion data of the type II-Plateau SN 2018aoq that exploded in NGC 4151. As distance estimates to NGC 4151 vary by an order of magnitude, we utilised the well-known correlation between ejecta velocity and plateau brightness, i.e. the standard candle method, to obtain a distance of 18.2$\pm$1.2 Mpc, which is in very good agreement with measurements based on geometric methods. The above distance implies a mid-plateau absolute magnitude of $M_{V}^{50}=-15.76\pm$0.14 suggesting that it is of intermediate brightness when compared to IIP SNe such as SN 2005cs at the faint end, and more typical events such as SN 1999em. This is further supported by relatively low expansion velocities (Fe II $\lambda$5169 $\sim$3000 km s$^{-1}$ at +42 d). Using archival HST/WFC3 imaging data, we find a point source coincident with the supernova position in the F350LP, F555W, F814W, and F160W filters. This source shows no significant variability over the $\sim$2 month time span of the data. From fits to the spectral energy distribution of the candidate progenitor, we find $\log\left(L/L_\odot\right)\sim 4.7$ and $T_{\mathrm{eff}}\sim 3.5$ kK, implying an M-type red supergiant progenitor. From comparisons to single and binary star models, we find that both favour the explosion of a star with a zero-age main sequence mass of $\sim$$10 M_\odot$.

Gravitational waves were discovered with the detection of binary black hole mergers and they should also be detectable from lower mass neutron star mergers. These are predicted to eject material rich in heavy radioactive isotopes that can power an electromagnetic signal called a kilonova. The gravitational wave source GW170817 arose from a binary neutron star merger in the nearby Universe with a relatively well confined sky position and distance estimate. Here we report observations and physical modelling of a rapidly fading electromagnetic transient in the galaxy NGC4993, which is spatially coincident with GW170817 and a weak short gamma-ray burst. The transient has physical parameters broadly matching the theoretical predictions of blue kilonovae from neutron star mergers. The emitted electromagnetic radiation can be explained with an ejected mass of 0.04 +/- 0.01 Msol, with an opacity of kappa <= 0.5 cm2/gm at a velocity of 0.2 +/- 0.1c. The power source is constrained to have a power law slope of beta = -1.2 +/- 0.3, consistent with radioactive powering from r-process nuclides. We identify line features in the spectra that are consistent with light r-process elements (90 < A < 140). As it fades, the transient rapidly becomes red, and emission may have contribution by a higher opacity, lanthanide-rich ejecta component. This indicates that neutron star mergers produce gravitational waves, radioactively powered kilonovae, and are a nucleosynthetic source of the r-process elements.