We used adaptive optics imaging and integral field spectroscopy from the Very Large Telescope, together with images from the Hubble Space Telescope, to study the near-infrared (NIR) evolution of the equatorial ring (ER) of SN 1987A. We studied the NIR flux and morphology over time in order to lay the groundwork for James Webb Space Telescope observations of the system. We also studied the differences in the interacting ring structure and flux between optical, NIR, and other wavelengths, and between line and continuum emission, to constrain the underlying physical processes. For the most part, the evolution is similar in the NIR and optical. The morphology of the ER has been skewed toward the west side (with roughly two-thirds of the NIR emission originating there) since around 2010. A steady decline in the ER flux, broadly similar to the mid-infrared and the optical, has been ongoing since roughly this time as well. The expansion velocity of the ER hotspots in the NIR is fully consistent with the optical. However, continuum emission forms roughly 70% of the NIR luminosity, and has been stronger outside the hotspot-defined extent of the ER (relative to the hotspots themselves) than the optical emission or the NIR line emission since 2012-2013, suggesting a faster-expanding continuum component. We find that this outer NIR emission can have a significant synchrotron contribution. Even if emission from hot dust (2000 K) is dominant within the ER, the mass of this dust must be vanishingly small (a few times 10(-12) M-circle dot) compared to the total dust mass in the ER (greater than or similar to 10(-5) M-circle dot) to account for the observed HKs flux. The NIR continuum emission, however, expands more slowly than the more diffuse 180-K dust emission that dominates in the MIR, indicating a different source, and the same hot dust component cannot account for the J-band emission.

Context. SN 2020qlb (ZTF20abobpcb) is a hydrogen-poor superluminous supernova (SLSN-I) that is among the most luminous (maximum M-g = -22.25 mag) and that has one of the longest rise times (77 days from explosion to maximum). We estimate the total radiated energy to be > 2.1 x 10(51) erg. SN 2020qlb has a well-sampled light curve that exhibits clear near and post peak undulations, a phenomenon seen in other SLSNe, whose physical origin is still unknown. Aims. We discuss the potential power source of this immense explosion as well as the mechanisms behind its observed light curve undulations. Methods. We analyze photospheric spectra and compare them to other SLSNe-I. We constructed the bolometric light curve using photometry from a large data set of observations from the Zwicky Transient Facility (ZTF), Liverpool Telescope (LT), and Neil Gehrels Swift Observatory and compare it with radioactive, circumstellar interaction and magnetar models. Model residuals and light curve polynomial fit residuals are analyzed to estimate the undulation timescale and amplitude. We also determine host galaxy properties based on imaging and spectroscopy data, including a detection of the [O III]lambda 4363, auroral line, allowing for a direct metallicity measurement. Results. We rule out the Arnett Ni-56 decay model for SN 2020qlb's light curve due to unphysical parameter results. Our most favored power source is the magnetic dipole spin-down energy deposition of a magnetar. Two to three near peak oscillations, intriguingly similar to those of SN 2015bn, were found in the magnetar model residuals with a timescale of 32 +/- 6 days and an amplitude of 6% of peak luminosity. We rule out centrally located undulation sources due to timescale considerations; and we favor the result of ejecta interactions with circumstellar material (CSM) density fluctuations as the source of the undulations.

During the Zwicky Transient Facility (ZTF) Phase I operations, 78 hydrogen-poor superluminous supernovae (SLSNe-I) were discovered in less than 3 yr, constituting the largest sample from a single survey. This paper (Paper I) presents the data, including the optical/UV light curves and classification spectra, while Paper II in this series will focus on the detailed analysis of the light curves and modeling. Our photometry is primarily taken by ZTF in the g, r, and i bands, and with additional data from other ground-based facilities and Swift. The events of our sample cover a redshift range of z = 0.06 - 0.67, with a median and 1 sigma error (16% and 84% percentiles) of zmed=0.265-0.135+0.143 M ( g,peak) <= -19.8 mag, with a median value of -21.48-0.61+1.13 t (rise) = 41.9 +/- 17.8 days. The luminosity and timescale distributions suggest that low-luminosity SLSNe-I with a peak luminosity & SIM;-20 mag or extremely fast-rising events (< 10 days) exist, but are rare. We confirm previous findings that slowly rising SLSNe-I also tend to fade slowly. The rest-frame color and temperature evolution show large scatters, suggesting that the SLSN-I population may have diverse spectral energy distributions. The peak rest-frame color shows a moderate correlation with the peak absolute magnitude, i.e., brighter SLSNe-I tend to have bluer colors. With optical and UV photometry, we construct the bolometric luminosity and derive a bolometric correction relation that is generally applicable for converting g, r-band photometry to the bolometric luminosity for SLSNe-I.

We present analysis of the light curves (LCs) of 77 hydrogen-poor superluminous supernovae (SLSNe I) discovered during the Zwicky Transient Facility Phase I operation. We find that the majority (67%) of the sample can be fit equally well by both magnetar and ejecta-circumstellar medium (CSM) interaction plus 56Ni decay models. This implies that LCs alone cannot unambiguously constrain the physical power sources for an SLSN I. However, 23% of the sample show inverted V-shape, steep-declining LCs or features of long rise and fast post-peak decay, which are better described by the CSM+Ni model. The remaining 10% of the sample favors the magnetar model. Moreover, our analysis shows that the LC undulations are quite common, with a fraction of 18%-44% in our gold sample. Among those strongly undulating events, about 62% of them are found to be CSM-favored, implying that the undulations tend to occur in the CSM-favored events. Undulations show a wide range in energy and duration, with median values (and 1σ errors) being as 1.7 % − 0.7 % + 1.5 % E rad , total and 28.8 − 9.1 + 14.4 days, respectively. Our analysis of the undulation timescales suggests that intrinsic temporal variations of the central engine can explain half of the undulating events, while CSM interaction (CSI) can account for the majority of the sample. Finally, all of the well-observed He-rich SLSNe Ib either have strongly undulating LCs or the LCs are much better fit by the CSM+Ni model. These observations imply that their progenitor stars have not had enough time to lose all of the He-envelopes before supernova explosions, and H-poor CSM are likely to present in these events.

We present spectroscopy of the ejecta of SN 1987A in 2017 and 2018 from the Hubble Space Telescope and the Very Large Telescope, covering the wavelength range between 1150 and 10000 angstrom. At 31 yr, this is the first epoch with coverage over the ultraviolet-to-near-infrared range since 1995. We create velocity maps of the ejecta in the H alpha, Mg II lambda lambda 2796, 2804 and [OI] lambda lambda 6302, 6366 (vacuum) emission lines and study their morphology. All three lines have a similar morphology, but Mg II is blueshifted by similar to 1000 km s(-1) relative to the others and stronger in the north-west. We also study the evolution of the line fluxes, finding a brightening by a factor of similar to 9 since 1999 in Mg II, while the other line fluxes are similar in 1999 and 2018. We discuss implications for the power sources of emission lines at late times: thermal excitation due to heating by the X-rays from the ejecta-ring interaction is found to dominate the ultraviolet Mg II lines, while the infrared Mg II doublet is powered mainly by Ly alpha fluorescence. The X-ray deposition is calculated based on merger models of SN 1987A. Far-ultraviolet emission lines of H-2 are not detected. Finally, we examine the combined spectrum of recently discovered hotspots outside the equatorial ring. Their unresolved Balmer emission lines close to zero velocity are consistent with the interaction of fast ejecta and a clumpy, slowly moving outflow. A clump of emission in this spectrum, south of the equatorial ring at similar to 1500 km s(-1), is likely associated with the reverse shock.

We present a sample of 14 hydrogen-rich superluminous supernovae (SLSNe II) from the Zwicky Transient Facility (ZTF) between 2018 and 2020. We include all classified SLSNe with peaks M-g < -20 mag with observed broad but not narrow Balmer emission, corresponding to roughly 20 per cent of all hydrogen-rich SLSNe in ZTF phase I. We examine the light curves and spectra of SLSNe II and attempt to constrain their power source using light-curve models. The brightest events are photometrically and spectroscopically similar to the prototypical SN 2008es, while others are found spectroscopically more reminiscent of non-superluminous SNe II, especially SNe II-L. Ni-56 decay as the primary power source is ruled out. Light-curve models generally cannot distinguish between circumstellar interaction (CSI) and a magnetar central engine, but an excess of ultraviolet (UV) emission signifying CSI is seen in most of the SNe with UV data, at a wide range of photometric properties. Simultaneously, the broad H alpha profiles of the brightest SLSNe II can be explained through electron scattering in a symmetric circumstellar medium (CSM). In other SLSNe II without narrow lines, the CSM may be confined and wholly overrun by the ejecta. CSI, possibly involving mass lost in recent eruptions, is implied to be the dominant power source in most SLSNe II, and the diversity in properties is likely the result of different mass loss histories. Based on their radiated energy, an additional power source may be required for the brightest SLSNe II, however - possibly a central engine combined with CSI.

We present multiwavelength modeling of the afterglow from the long gamma-ray burst (GRB) 160625B using Markov Chain Monte Carlo techniques of the afterglowpy Python package. GRB 160625B is an extremely bright burst with a rich set of observations spanning from radio to gamma-ray frequencies. These observations range from similar to 0.1 days to >1000 days, thus making this event extremely well suited to such modeling. In this work we compare top-hat and Gaussian jet structure types in order to find best-fit values for the GRB jet collimation angle, viewing angle, and other physical parameters. We find that a Gaussian-shaped jet is preferred (2.7 sigma-5.3 sigma) over the traditional top-hat model. Our estimate for the opening angle of the burst ranges from 126 to 390, depending on jet-shape model. We also discuss the implications that assumptions on jet shape, viewing angle, and particularly the participation a fraction of electrons have on the final estimation of GRB intrinsic energy release and the resulting energy budget of the relativistic outflow. Most notably, allowing the participation fraction to vary results in an estimated total relativistic energy of similar to 10(53) erg. This is two orders of magnitude higher than when the total fraction is assumed to be unity; thus, this parameter has strong relevance for placing constraints on long GRB central engines, details of the circumburst media, and host environment.