Draft version August 25, 2021 Typeset using LATEX twocolumn style in AASTeX631 A VLA Survey of Late-time Radio Emission from Superluminous Supernovae and the Host Galaxies Bunyo Hatsukade ,1 Nozomu Tominaga ,2, 3, 4 Tomoki Morokuma ,1 Kana Morokuma-Matsui ,1 Yuichi Matsuda ,2, 5 Yoichi Tamura ,6 Kotaro Niinuma ,7 and Kazuhiro Motogi 7 1Institute of Astronomy, Graduate School of Science, The University of Tokyo, 2-21-1 Osawa, Mitaka, Tokyo 181-0015, Japan 2National Astronomical Observatory of Japan, 2-21-1 Osawa, Mitaka, Tokyo 181-8588, Japan 3Department of Physics, Faculty of Science and Engineering, Konan University, 8-9-1 Okamoto, Kobe, Hyogo 658-8501, Japan 4Kavli Institute for the Physics and Mathematics of the Universe (WPI), The University of Tokyo, 5-1-5 Kashiwanoha, Kashiwa, Chiba 277-8583, Japan 5Graduate University for Advanced Studies (SOKENDAI), Osawa 2-21-1, Mitaka, Tokyo 181-8588, Japan 6Department of Physics, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8602, Japan 7Graduate School of Sciences and Technology for Innovation, Yamaguchi University, Yoshida 1677-1, Yamaguchi, Yamaguchi 753-8512, Japan (Received June 25, 2021; Revised August 12, 2021; Accepted August 23, 2021) ABSTRACT We present the results of 3 GHz radio continuum observations of 23 superluminous supernovae (SLSNe) and their host galaxies by using the Karl G. Jansky Very Large Array conducted 5{21 years after the explosions. The sample consists of 15 Type I and 8 Type II SLSNe at z < 0:3, providing one of the largest sample of SLSNe with late-time radio data. We detected radio emission from one SLSN (PTF10hgi) and 5 hosts with a significance of >5σ. No time variability is found in late-time radio light curves of the radio-detected sources in a timescale of years except for PTF10hgi, whose variability is reported in a separate study. Comparison of star-formation rates (SFRs) derived from the 3 GHz flux densities with those derived from SED modeling based on UV{NIR data shows that four hosts have an excess of radio SFRs, suggesting obscured star formation. Upper limits for undetected hosts and stacked results show that the majority of the SLSN hosts do not have a significant obscured star formation. By using the 3 GHz upper limits, we constrain the parameters for afterglows arising from interaction between initially off-axis jets and circumstellar medium (CSM). We found that the models 53 −3 with higher energies (Eiso & several × 10 erg) and CSM densities (n & 0:01 cm ) are excluded, but lower energies or CSM densities are not excluded with the current data. We also constrained the models of pulsar wind nebulae powered by a newly born magnetar for a subsample of SLSNe with model predictions in the literature. Keywords: Radio transient sources (2008); Extragalactic radio sources (508); Supernovae (1668); Circumstellar matter (241); Radio continuum emission (1340); Very Large Array (1766); Radio Astronomy (1338) arXiv:2108.10445v1 [astro-ph.HE] 23 Aug 2021 1. INTRODUCTION 2012; Moriya et al. 2019), and therefore can be pow- Superluminous supernovae (SLSNe) are very bright erful probes of the distant universe. The power source explosions (>10{100 times brighter) and rare events and mechanism for large luminosities are still a matter (<100 times lower) compared to ordinary SNe (Gal-Yam of debate. SLSNe are classified into two types accord- 2012, 2019; Moriya et al. 2018, for reviews). SLSNe are ing to the spectra: hydrogen-poor Type I SLSNe (or detected at high redshifts up to z ∼ 4 (Cooke et al. SLSNe-I) and hydrogen-rich Type II SLSNe (or SLSNe- II). SLSNe-II can be explained as an interaction between the SN ejecta and dense circumstellar medium (CSM; [email protected] e.g., Woosley et al. 2007), while a number of models have been proposed for SLSNe-I: a large amount of 56Ni pro- 2 Hatsukade et al. duced by a pair-instability SN (e.g., Woosley et al. 2007), ble scenarios are a low-luminosity active galactic nucleus spin-down of a newborn strongly magnetized neutron (AGN) in the host galaxy or a magnetar wind nebula by star (magnetar; e.g., Kasen & Bildsten 2010; Woosley considering the radio light curve and spectrum. 2010), fallback accretion onto a compact remnant (e.g., In order to constrain the models, it is also impor- Dexter & Kasen 2013), and interaction with dense CSM tant to understand the properties of their host galax- (e.g., Smith & McCray 2007; Chevalier & Irwin 2011). ies. Previous studies have shown that SLSN-I hosts are Radio observations can probe synchrotron emission typically dwarf galaxies with low-luminosity, low stellar arising from shock interaction of SN ejecta or jet with mass, low star-formation rate (SFR), and high specific surrounding material, providing useful constraints on SFR (sSFR) compared to local star-forming galaxies and the models of SLSNe. It is expected that late-time the hosts of core-collapse SNe, while SLSN-II hosts show radio emission may be caused by initially off-axis jets a wider range of those parameters (e.g., Lunnan et al. that decelerate and spread into the line of sight. It is 2014; Leloudas et al. 2015a; Angus et al. 2016; Per- also predicted that quasistate synchrotron emission may ley et al. 2016; Chen et al. 2017; Schulze et al. 2018; arise from pulsar wind nebulae (PWNe) powered by a Taggart & Perley 2021). The observations of SLSN newly born magnetar (Metzger & Bower 2014; Murase hosts have been conducted mainly in the optical/near- et al. 2016; Kashiyama & Murase 2017; Metzger et al. infrared (NIR) wavelengths, which are subject to dust 2017). The radio emission is initially absorbed in the extinction in contrast to longer wavelengths, and it is PWNe and SN ejecta, but the system can be transparent possible that we are missing dust-obscured star forma- on timescales of decades. These central engine models tion. Radio observations are important to probe the are also a possible model for long-duration gamma-ray properties of star formation in the SLSN host without bursts (LGRBs) (e.g., Metzger et al. 2015; Margalit et al. the effect of dust extinction (Schulze et al. 2018; Hat- 2018b), suggesting a connection between SLSNe-I and sukade et al. 2018; Arabsalmani et al. 2019; Hatsukade LGRBs. It is interesting that the magnetar model is also et al. 2020; Eftekhari et al. 2021). Schulze et al.(2018) one of the plausible models for fast radio bursts (FRBs), searched radio emission for a sample of SLSN hosts which are mysterious radio transients with millisecond- by using the survey data of Faint Images of the Ra- scale bright flashes (Cordes & Chatterjee 2019, for a re- dio Sky at Twenty-Centimeters (FIRST; Becker et al. view). Magnetar models have been applied to the origin 1995), the NRAO VLA Sky Survey (NVSS; Condon of FRB 121102 (e.g., Murase et al. 2016; Kashiyama & et al. 1998), and the Sydney University Molonglo Sky Murase 2017; Metzger et al. 2017; Margalit et al. 2018a; Survey (SUMSS; Bock et al. 1999), and did not find Margalit & Metzger 2018). Recently, FRB 200428 was radio detection. They also obtained upper limits on identified as a Galactic magnetar, SGR 1935+2154 (An- three SLSN hosts based on deep 1.5 GHz observations dersen et al. 2020; Bochenek et al. 2020), showing a mag- with the Karl G. Jansky Very Large Array (VLA). Hat- netar origin of at least one FRB. sukade et al.(2018) performed 3 GHz radio observations Previous radio observations of SLSNe resulted in non- of the host galaxies of 8 SLSNe (5 Type I and 3 Type detections in most cases, providing constraints on phys- II) at 0:1 < z < 0:3 with VLA. They found the ex- ical properties (such as energies, mass-loss rates, and cess of radio-based SFRs over extinction-corrected op- CSM densities) and models for SLSNe-I (e.g., Nicholl tically based SFRs (UV-based SFRs from SED fitting et al. 2016, 2018; Margalit et al. 2018b; Bose et al. 2018; or Hα-based SFRs) in three hosts, suggesting the exis- Coppejans et al. 2018; Law et al. 2019; Eftekhari et al. tence of dust-obscured star formation that were missed 2021). Eftekhari et al.(2019) found an unresolved radio in previous observations. This indicates the necessity of source coincident with the position of SLSN-I PTF10hgi. longer-wavelength observations in order to understand They argued that the radio emission is consistent with true star-forming activity in SLSN hosts. They also an off-axis jet or wind nebula powered by a magnetar found that three hosts, which were located within the born in the SLSN, suggesting the presence of a central range of the main sequence based on previous optical ob- engine. Law et al.(2019) and Mondal et al.(2020) also servations, are actually above the main sequence in our supported the model of magnetar-powered SLSN based radio observations, suggesting that they have a starburst on the observations at 0.6{15 GHz. Hatsukade et al. nature. Eftekhari et al.(2021) observed 15 SLSNe-I at 6 (2021) found a variability in late-time radio emission in GHz with VLA and 29 SLSNe-I at 100 GHz with the At- PTF10hgi for the first time among SLSNe. They con- acama Large Millimeter/submillimeter Array (ALMA), strained both the rise and decay phases of the radio seven of which overlapped. They did not detect emission light curve over three years, peaking at approximately from any of the SLSNe except for a 6 GHz detection of 8{9 yr after the explosion. They concluded that plausi- PTF10hgi.