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.

We present a systematic search for optical counterparts to 13 gravitational wave (GW) triggers involving at least one neutron star during LIGO/Virgo's third observing run. We searched binary neutron star (BNS) and neutron star black hole (NSBH) merger localizations with the Zwicky Transient Facility (ZTF) and undertook follow-up with the Global Relay of Observatories Watching Transients Happen (GROWTH) collaboration. The GW triggers had a median localization of 4480 deg^2, median distance of 267 Mpc and false alarm rates ranging from 1.5 to 1e-25 per yr. The ZTF coverage had a median enclosed probability of 39%, median depth of 20.8mag, and median response time of 1.5 hr. The O3 follow-up by the GROWTH team comprised 340 UVOIR photometric points, 64 OIR spectra, and 3 radio. We find no promising kilonova (radioactivity-powered counterpart) and we convert the upper limits to constrain the underlying kilonova luminosity function. Assuming that all kilonovae are at least as luminous as GW170817 at discovery (-16.1mag), we calculate our joint probability of detecting zero kilonovae is only 4.2%. If we assume that all kilonovae are brighter than -16.6mag (extrapolated peak magnitude of GW170817) and fade at 1 mag/day (similar to GW170817), the joint probability of zero detections is 7%. If we separate the NSBH and BNS populations, the joint probability of zero detections, assuming all kilonovae are brighter than -16.6mag, is 9.7% for NSBH and 7.9% for BNS mergers. Moreover, <57% (<89%) of putative kilonovae could be brighter than -16.6mag assuming flat (fading) evolution, at 90% confidence. If we further account for the online terrestrial probability for each GW trigger, we find that <68% of putative kilonovae could be brighter than -16.6mag. Comparing to model grids, we find that some kilonovae must have Mej < 0.03 Msun or Xlan>1e-4 or phi>30deg to be consistent with our limits. (Abridged)

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.

In this paper we report the results of the first $\sim$four years of spectroscopic and photometric monitoring of the Type IIn supernova SN 2015da (also known as PSN J13522411+3941286, or iPTF16tu). The supernova exploded in the nearby spiral galaxy NGC 5337 in a relatively highly extinguished environment. The transient showed prominent narrow Balmer lines in emission at all times and a slow rise to maximum in all bands. In addition, early observations performed by amateur astronomers give a very well-constrained explosion epoch. The observables are consistent with continuous interaction between the supernova ejecta and a dense and extended H-rich circumstellar medium. The presence of such an extended and dense medium is difficult to reconcile with standard stellar evolution models, since the metallicity at the position of SN 2015da seems to be slightly subsolar. Interaction is likely the mechanism powering the light curve, as confirmed by the analysis of the pseudo bolometric light curve, which gives a total radiated energy $\gtrsim10^{51}\,\rm{erg}$. Modeling the light curve in the context of a supernova shock breakout through a dense circumstellar medium allowed us to infer the mass of the prexisting gas to be $\simeq8\,\rm{M_{\odot}}$, with an extreme mass-loss rate for the progenitor star $\simeq0.6\,\rm{M_{\odot}}\,\rm{yr^{-1}}$, suggesting that most of the circumstellar gas was produced during multiple eruptive events. Near- and mid-infrared observations reveal a flux excess in these domains, similar to those observed in SN 2010jl and other interacting transients, likely due to preexisting radiatively heated dust surrounding the supernova. By modeling the infrared excess, we infer a mass $\gtrsim0.4\times10^{-3}\,\rm{M_{\odot}}$ for the dust.

We report deep radio observations of nearby Type Ia Supernovae (SNe Ia) with the electronic Multi-Element Radio Linked Interferometer Net-work (e-MERLIN), and the Australia Telescope Compact Array (ATCA). No detections were made. With standard assumptions for the energy densities of relativistic electrons going into a power-law energy distribution, and the magnetic field strength (epsilon_e = epsilon_B = 0.1), we arrive at the upper limits on mass-loss rate for the progenitor system of SN 2013dy (2016coj, 2018gv, 2018pv, 2019np), to be less than 12 (2.8,1.3, 2.1, 1.7)EE(-8) solar masses per year (for a wind speed of 100 km/s). To SNe 2016coj, 2018gv, 2018pv and 2019np we add radio data for 17 other nearby SNe Ia, and model their non-detections. With the same model as described, all 21 SNe Ia have mass-loss rates less than 4EE(-8) solar masses per year (for a wind speed of 100 km/s). We compare those limits with the expected mass loss rates in different single-degenerate progenitor scenarios. We also discuss how information on epsilon_e and epsilon_B can be obtained from late observations of SNe Ia and the youngest SN Ia remnant detected in radio, G1.9+0.3, as well as stripped-envelope core-collapse SNe. We highlight SN 2011dh, and argue for epsilon_e approximately equal to 0.1 and epsilon_B approximately equal to 0.0033. Finally, we discuss strategies to observe at radio frequencies to maximize the chance of detection, given the time since explosion, the distance to the supernova and the telescope sensitivity.

On 2019 August 14, the Advanced LIGO and Virgo interferometers detected the high-significance gravitational wave (GW) signal S190814bv. The GW data indicated that the event resulted from a neutron star--black hole (NSBH) merger, or potentially a low-mass binary black hole merger. Due to the low false alarm rate and the precise localization (23 deg$^2$ at 90\%), S190814bv presented the community with the best opportunity yet to directly observe an optical/near-infrared counterpart to a NSBH merger. To search for potential counterparts, the GROWTH collaboration performed real-time image subtraction on 6 nights of public Dark Energy Camera (DECam) images acquired in the three weeks following the merger, covering $>$98\% of the localization probability. Using a worldwide network of follow-up facilities, we systematically undertook spectroscopy and imaging of optical counterpart candidates. Combining these data with a photometric redshift catalog, we ruled out each candidate as the counterpart to S190814bv and we placed deep, uniform limits on the optical emission associated with S190814bv. For the nearest consistent GW distance, radiative transfer simulations of NSBH mergers constrain the ejecta mass of S190814bv to be $M_\mathrm{ej} < 0.04$~$M_{\odot}$ at polar viewing angles, or $M_\mathrm{ej} < 0.03$~$M_{\odot}$ if the opacity is $\kappa < 2$~cm$^2$g$^{-1}$. Assuming a tidal deformability for the neutron star at the high end of the range compatible with GW170817 results, our limits would constrain the BH spin component aligned with the orbital momentum to be $ \chi < 0.7$ for mass ratios $Q < 6$, with weaker constraints for more compact neutron stars. We publicly release the photometry from this campaign at http://www.astro.caltech.edu/~danny/static/s190814bv.

The first two months of the third Advanced LIGO and Virgo observing run (2019 April-May) showed that distant gravitational wave (GW) events can now be readily detected. Three candidate mergers containing neutron stars (NS) were reported in a span of 15 days, all likely located more than 100 Mpc away. However, distant events such as the three new NS mergers are likely to be coarsely localized, which highlights the importance of facilities and scheduling systems that enable deep observations over hundreds to thousands of square degrees to detect the electromagnetic counterparts. On 2019-05-10 02:59:39.292 UT the GW candidate S190510g was discovered and initially classified as a BNS merger with 98% probability. The GW event was localized within an area of 3462 deg2, later refined to 1166 deg2 (90%) at a distance of 227 +- 92 Mpc. We triggered Target of Opportunity observations with the Dark Energy Camera (DECam), a wide-field optical imager mounted at the prime focus of the 4m Blanco Telescope at CTIO in Chile. This Letter describes our DECam observations and our real-time analysis results, focusing in particular on the design and implementation of the observing strategy. Within 24 hours of the merger time, we observed 65% of the total enclosed probability of the final skymap with an observing efficiency of 94%. We identified and publicly announced 13 candidate counterparts. S190510g was re-classified 1.7 days after the merger, after our observations were completed, with a "binary neutron star merger" probability reduced from 98% to 42% in favor of a "terrestrial" classification.

Only three extragalactic supernovae have been detected at late times at millimeter wavelengths: SN 1987A, SN 1978K, and SN 1996cr. SN 1978K is a remarkably luminous Type IIn supernova that remains bright at all wavelengths 40 years after its explosion. Here we present Atacama Large Millimeter/submillimeter Array (ALMA) observations taken in 2016 using Bands 3, 4, 6, and 7 that show a steepening in the spectrum. An absorbed single power law model broadly fits all the radio and millimeter observations, but would require significant chromatic variability. Alternatively, a broken power law fits the radio-millimeter spectrum: this can be explained using an ultra-relativistic spherical blast wave in a wind scaling with a cooling break, as in a gamma-ray burst afterglow. Using updated Australia Telescope Compact Array (ATCA) light curves, we show the non-thermal radio continuum continues to decay as $t^{-1.53}$; in the fireball model, this independently defines the power law indices found in the radio-millimeter spectrum. Supernovae such as SN 1978K might be important contributors to the Universal dust budget: only SN 1978K was detected in a search for warm dust in supernovae in the transitional phase (age 10-100 years). Using Spitzer Space Telescope observations, we show that at least some of this dust emission has been decaying rapidly as $t^{-2.45}$ over the past decade, suggesting it is being destroyed. Depending on the modeling of the synchrotron emission, the ALMA observations suggest there may be emission from a cold dust component.

Core collapse supernova (CCSN) rates suffer from large uncertainties as many CCSNe exploding in regions of bright background emission and significant dust extinction remain unobserved. Such a shortfall is particularly prominent in luminous infrared galaxies (LIRGs), which have high star formation (and thus CCSN) rates and host bright and crowded nuclear regions, where large extinctions and reduced search detection efficiency likely lead to a significant fraction of CCSNe remaining undiscovered. We present the first results of project SUNBIRD (Supernovae UNmasked By Infra-Red Detection), where we aim to uncover CCSNe that otherwise would remain hidden in the complex nuclear regions of LIRGs, and in this way improve the constraints on the fraction that is missed by optical seeing-limited surveys. We observe in the near-infrared 2.15 {\mu}m $K_s$-band, which is less affected by dust extinction compared to the optical, using the multi-conjugate adaptive optics imager GeMS/GSAOI on Gemini South, allowing us to achieve a spatial resolution that lets us probe close in to the nuclear regions. During our pilot program and subsequent first full year we have discovered three CCSNe and one candidate with projected nuclear offsets as small as 200 pc. When compared to the total sample of LIRG CCSNe discovered in the near-IR and optical, we show that our method is singularly effective in uncovering CCSNe in nuclear regions and we conclude that the majority of CCSNe exploding in LIRGs are not detected as a result of dust obscuration and poor spatial resolution.

Nuclear starbursts and AGN activity are the main heating processes in luminous infrared galaxies (LIRGs) and their relationship is fundamental to understand galaxy evolution. In this paper, we study the star-formation and AGN activity of a sample of 11 local LIRGs imaged with subarcsecond angular resolution at radio (8.4GHz) and near-infrared ($2.2\mu$m) wavelengths. This allows us to characterize the central kpc of these galaxies with a spatial resolution of $\simeq100$pc. In general, we find a good spatial correlation between the radio and the near-IR emission, although radio emission tends to be more concentrated in the nuclear regions. Additionally, we use an MCMC code to model their multi-wavelength spectral energy distribution (SED) using template libraries of starburst, AGN and spheroidal/cirrus models, determining the luminosity contribution of each component, and finding that all sources in our sample are starburst-dominated, except for NGC6926 with an AGN contribution of $\simeq64$\%. Our sources show high star formation rates ($40$ to $167M_\odot\mathrm{yr}^{-1}$), supernova rates (0.4 to $2.0\mathrm{SN}\mathrm{yr}^{-1}$), and similar starburst ages (13 to $29\mathrm{Myr}$), except for the young starburst (9Myr) in NGC6926. A comparison of our derived star-forming parameters with estimates obtained from different IR and radio tracers shows an overall consistency among the different star formation tracers. AGN tracers based on mid-IR, high-ionization line ratios also show an overall agreement with our SED model fit estimates for the AGN. Finally, we use our wide-band VLA observations to determine pixel-by-pixel radio spectral indices for all galaxies in our sample, finding a typical median value ($\alpha\simeq-0.8$) for synchrotron-powered LIRGs.

Supernova 1978K is one of the oldest-known examples of the class of Type IIn supernovae that show evidence for strong interaction between the blast wave and a dense, pre-existing circumstellar medium. Here we report detections of SN 1978K at both 34 GHz and 94 GHz, making it only the third extragalactic supernova (after SN 1987A and SN 1996cr) to be detected at late-times at these frequencies. We find SN 1978K to be >400 times more luminous than SN 1987A at millimetre wavelengths in spite of the roughly nine year difference in ages, highlighting the risk in adopting SN 1987A as a template for the evolution of core-collapse supernovae in general. Additionally, from new VLBI observations at 8.4 GHz, we measure a deconvolved diameter for SN 1978K of ~5 milli-arcsec, and a corresponding average expansion velocity of <1500 km/s. These observations provide independent evidence of an extremely dense circumstellar medium surrounding the progenitor star.

In this paper we present the results of the first systematic search for counterparts to nearby ultraluminous X-ray sources (ULXs) in the near-infrared (NIR). We observed 62 ULXs in 37 galaxies within 10 Mpc and discovered 17 candidate NIR counterparts. The detection of 17 out of 62 ULX candidates points to intrinsic differences between systems that show and those that do not show infrared emission. For six counterparts we conclude from the absolute magnitudes and - - in some cases - additional information such as morphology and previously reported photometric or spectroscopic observations, that they are likely background active galactic nuclei or ULXs residing in star clusters. Eleven counterparts have absolute magnitudes consistent with them being single red supergiant stars. Alternatively, these systems may have larger accretion discs that emit more NIR light than the systems that we do not detect. Other scenarios such as emission from a surrounding nebula or from a compact radio jet are also possible, although for Holmberg II X-1 the NIR luminosity far exceeds the expected jet contribution. The eleven possible red supergiant counterparts are excellent candidates for spectroscopic follow-up observations. This may enable us to measure the mass function in these systems if they are indeed red supergiant donor stars where we can observe absorption lines.