On 2019 August 14, the LIGO and Virgo interferometers detected a high-significance event labelled S190814bv. Preliminary analysis of the GW data suggests that the event was likely due to the merger of a compact binary system formed by a BH and a NS. ElectromagNetic counterparts of GRAvitational wave sources at the VEry Large Telescope (ENGRAVE) collaboration members carried out an intensive multi-epoch, multi-instrument observational campaign to identify the possible optical/near infrared counterpart of the event. In addition, the ATLAS, GOTO, GRAWITA-VST, Pan-STARRS and VINROUGE projects also carried out a search on this event. Our observations allow us to place limits on the presence of any counterpart and discuss the implications for the kilonova (KN) possibly generated by this NS-BH merger, and for the strategy of future searches. Altogether, our observations allow us to exclude a KN with large ejecta mass $M\gtrsim 0.1\,\mathrm{M_\odot}$ to a high ($>90\%$) confidence, and we can exclude much smaller masses in a subsample of our observations. This disfavours the tidal disruption of the neutron star during the merger. Despite the sensitive instruments involved in the campaign, given the distance of S190814bv we could not reach sufficiently deep limits to constrain a KN comparable in luminosity to AT 2017gfo on a large fraction of the localisation probability. This suggests that future (likely common) events at a few hundreds Mpc will be detected only by large facilities with both high sensitivity and large field of view. Galaxy-targeted observations can reach the needed depth over a relevant portion of the localisation probability with a smaller investment of resources, but the number of galaxies to be targeted in order to get a fairly complete coverage is large, even in the case of a localisation as good as that of this event.

The interaction of a supernova with a circumstellar medium (CSM) can dramatically increase the emitted luminosity by converting kinetic energy to thermal energy. In 'superluminous' supernovae (SLSNe) of Type IIn -- named for narrow hydrogen lines in their spectra -- the integrated emission can reach $\sim 10^{51}$ erg, attainable by thermalising most of the kinetic energy of a conventional SN. A few transients in the centres of active galaxies have shown similar spectra and even larger energies, but are difficult to distinguish from accretion onto the supermassive black hole. Here we present a new event, SN2016aps, offset from the centre of a low-mass galaxy, that radiated $\gtrsim 5 \times 10^{51}$ erg, necessitating a hyper-energetic supernova explosion. We find a total (SN ejecta $+$ CSM) mass likely exceeding 50-100 M$_\odot$, with energy $\gtrsim 10^{52}$ erg, consistent with some models of pair-instability supernovae (PISNe) or pulsational PISNe -- theoretically-predicted thermonuclear explosions from helium cores $>50$ M$_\odot$. Independent of the explosion mechanism, this event demonstrates the existence of extremely energetic stellar explosions, detectable at very high redshifts, and provides insight into dense CSM formation in the most massive stars.

We present the first observations of a Type I superluminous supernova (SLSN) at $\gtrsim 1000$ days after maximum light. We observed SN 2015bn using the Hubble Space Telescope Advanced Camera for Surveys in the F475W, F625W and F775W filters at 721 days and 1068 days. SN 2015bn is clearly detected and resolved from its compact host, allowing reliable photometry. A galaxy template constructed from these data further enables us to isolate the SLSN flux in deep ground-based imaging. We measure a light curve decline rate at $>700$ days of $0.19 \pm 0.03$ mag (100 d)$^{-1}$, much shallower than the earlier evolution, and slower than previous SLSNe (at any phase) or the decay rate of $^{56}$Co. Neither additional radioactive isotopes nor a light echo can consistently account for the slow decline. A spectrum at 1083 days shows the same [O I] and [Ca II] lines as seen at $\sim300-400$ days, with no new features to indicate strong circumstellar interaction. Radio limits with the Very Large Array rule out an extended wind for mass-loss rates $10^{-2.7} \lesssim \dot{M}/v_{10} \lesssim 10^{-1.1}$ M$_\odot$ yr$^{-1}$ (where $v_{10}$ is the wind velocity in units of 10 km s$^{-1}$). The optical light curve is consistent with $L \propto t^{-4}$, which we show is expected for magnetar spin-down with inefficient trapping; furthermore, the evolution matches predictions from earlier magnetar model fits. The opacity to magnetar radiation is constrained at $\sim 0.01$ cm$^2$ g$^{-1}$, consistent with photon-matter pair-production over a broad $\sim$GeV-TeV range. This suggests the magnetar spectral energy distribution, and hence the 'missing energy' leaking from the ejecta, may peak in this range.

SN 2015bn is a nearby hydrogen-poor superluminous supernova (SLSN-I) that has been intensively observed in X-rays with the goal to detect the spin-down powered emission from a magnetar engine. The early-time UV/optical/infrared (UVOIR) data fit well to the magnetar model, but require leakage of energy at late times of $\lesssim 10^{43}$ erg s$^{-1}$, which is expected to be partially emitted in X-rays. Deep X-ray limits until $\sim$300 days after explosion revealed no X-ray emission. Here, we present the latest deep 0.3--10 keV X-ray limit at 805 days obtained with \textit{XMM-Newton}. We find $L_X < 10^{41}$ erg s$^{-1}$, with no direct evidence for central-engine powered emission. While the late-time optical data still follow the prediction of the magnetar model, the best-fit model to the bolometric light curve predicts that $\sim$97\% of the total input luminosity of the magnetar is escaping outside of the UVOIR bandpass at the time of observation. Our X-ray upper limit is $<$1.5\% of the input luminosity, strongly constraining the high-energy leakage, unless non-radiative losses are important. These deep X-ray observations identify a missing energy problem in SLSNe-I and we suggest future observations in hard X-rays and $\gamma$-rays for better constraints. Also, independent of the optical data, we constrain the parameter spaces of various X-ray emission scenarios, including ionization breakout by magnetar spin-down, shock interaction between the ejecta and external circumstellar medium, off-axis $\gamma$-ray burst afterglow, and black hole fallback accretion.

SN 2008es is one of the rare cases of a Type II superluminous supernova (SLSN) showing no relatively narrow features in its early-time spectra, and therefore its powering mechanism is under debate between circumstellar interaction (CSI) and magnetar spin-down. Late-time data are required for better constraints. We present optical and near-infrared (NIR) photometry obtained from Gemini, Keck, and Palomar Observatories from 192 to 554 days after explosion. Only broad H$\alpha$ emission is detected in a Gemini spectrum at 288 days. The line profile exhibits red-wing attenuation relative to the early-time spectrum. In addition to the cooling SN photosphere, a NIR excess with blackbody temperature $\sim1500$ K and radius $\sim10^{16}$ cm is observed. This evidence supports dust condensation in the cool dense shell being responsible for the spectral evolution and NIR excess. We favour CSI, with $\sim2$--3 $\textrm{M}_\odot$ of circumstellar material (CSM) and $\sim$10--20 $\textrm{M}_\odot$ of ejecta, as the powering mechanism, which still dominates at our late-time epochs. Both models of uniform density and steady wind fit the data equally well, with an effective CSM radius $\sim 10^{15}$ cm, supporting the efficient conversion of shock energy to radiation by CSI. A low amount ($\lesssim 0.4$ $\textrm{M}_\odot$) of $^{56}$Ni is possible but cannot be verified yet, since the light curve is dominated by CSI. The magnetar spin-down powering mechanism cannot be ruled out, but is less favoured because it overpredicts the late-time fluxes and may be inconsistent with the presence of dust.

[Abridged] We present UV/optical observations of PS16aqv (SN 2016ard), a Type I superluminous supernova (SLSN-I) classified as part of our search for low-$z$ SLSNe. PS16aqv is a fast evolving SLSNe-I that reached a peak absolute magnitude of $M_{r} \approx -22.1$. The lightcurves exhibit a significant undulation at 30 rest-frame days after peak, with a behavior similar to undulations seen in the slowly fading SLSN-I SN 2015bn. This similarity strengthens the case that fast and slow SLSNe-I form a continuum with a common origin. At $\approx\!80$ days after peak, the lightcurves exhibit a transition to a slow decline, followed by significant subsequent steepening, indicative of a plateau phase or a second significant undulation. Deep limits at $\approx280$ days after peak imply a tight constraint on the nickel mass, $M_{\rm Ni} \lesssim 0.35$ M$_{\odot}$ (lower than for previous SLSNe-I), and indicate that some SLSNe-I do not produce significantly more nickel than normal Type Ic SNe. Using MOSFiT, we model the lightcurve with a magnetar central engine model and find $P_{\rm spin} \approx 0.9$ ms, $B \approx 1.5 \times 10^{14}$ G, and $M_{\rm ej} \approx 16$ M$_{\odot}$. The implied rapid spin-down time and large reservoir of available energy coupled with the high ejecta mass may account for the fast evolving lightcurve and slow spectroscopic evolution. We also study the location of PS16aqv in its host galaxy and find that it occurred at an offset of $2.46 \pm 0.21$ kpc from the central star-forming region. We find the host galaxy exhibits low metallicity and spatially varying extinction and star formation rate, with the explosion site exhibiting lower values than the central region. The complexity seen in the lightcurves of PS16aqv and other events highlights the importance of obtaining well-sampled lightcurves for exploring deviations from a uniform decline.