Discovery of a Relativistic Stripped Envelope Type Ic-BL Supernova at z = 2.83 with JWST

Abstract

We present JWST NIRCam and NIRSpec observations of a Type Ic supernova (SN Ic) and its host galaxy (JADES-GS+53.13533-27.81457) at $z=2.83$. This SN (named SN 2023adta) was identified in deep James Webb Space Telescope (JWST)/NIRCam imaging from the JWST Advanced Deep Extragalactic Survey (JADES) Program. Follow-up observations with JWST/NIRSpec provided a spectroscopic redshift of $z=2.83$ and the classification as a SN Ic-BL. The light curve of SN 2023adta matches well with other stripped envelope supernovae and we find a high peak luminosity, $M_{V}=-19.0\pm 0.2$ mag, based on the distribution of best-fit SNe. The broad absorption features in its spectrum are consistent with other SNe Ic-BL 1-3 weeks after peak brightness. We measure a Ca II NIR triplet expansion velocity of $29{,}000\pm 2{,}000$  km s^-1. The host galaxy of SN 2023adta is irregular, and modeling of its spectral energy distribution (SED) indicates a metallicity of $Z=0.35^{+0.16}_{-0.08}Z_{\odot}$. This environment is consistent with the population of low-$z$ SNe Ic-BL which prefer lower metallicities relative to other stripped envelope supernovae, and track long duration $\gamma$-ray burst (LGRB) environments. We do not identify any GRBs that are coincident with SN 2023adta. Given the rarity of SNe Ic-BL in the local universe, the detection of a SN Ic-BL at $z=2.83$ could indicate that their rates are enhanced at high redshift.

1 Introduction

Stripped-envelope supernovae (SESNe) result from the core-collapse of massive stars that have experienced some form of mass-loss causing them to lose their outer layers of hydrogen and sometimes helium (Clocchiatti et al., 1996). A Type Ic supernova (SN Ic) is a SESN characterized by the lack of both H and He features in its spectra (Filippenko et al., 1995; Matheson et al., 2001) and requires a progenitor star that has experienced a large amount of stripping (see, Woosley et al., 1994; Dessart et al., 2012; Gal-Yam, 2017; Dessart et al., 2020). Two primary channels with different mass-loss mechanisms are thought to be able to produce a SN Ic. First, the collapse of a high-mass Wolf-Rayet (WR) star (Woosley & Weaver, 1995; Georgy et al., 2009), whose strong metal line blanketing-driven winds have removed its outer layers. Progenitor angular momentum and metallicity likely play an important role in these progenitors because they result in stronger stellar winds. Second, young massive stars in close binary systems whose mass-loss occurs through Roche-Lobe overflow or common envelope evolution (Podsiadlowski et al., 1992; Yoon et al., 2010; Sana et al., 2012; Lyman et al., 2016). It is important to note that these channels are not mutually-exclusive, and a combination of both mass-loss mechanisms is likely needed to explain the diversity of SESNe (Smith et al., 2015).

There have been limited searches to directly detect companion stars. To date, there are only two SNe Ic nearby enough to detect any potential companion. Deep post-explosion upper limits of the Type Ic SN 1994I, which had earlier data restricting any single star scenario, favored companion scenarios with non-conservative mass transfer with intermediate initial orbital periods and mass ratios (Van Dyk et al., 2016). Deep, post-explosion observations of the fully stripped Type Ib/c SN 2013ge resulted in a direct detection of a surviving companion consistent with a slightly reddened post-main sequence 12 M_⊙ star (Fox et al., 2022). Aside from these two post-explosion observations, pre-explosion spectral energy distributions (SEDs) and detailed modeling suggest potential companions in the Type Ib iPTF13bvn (e.g., Cao et al., 2013; Eldridge & Maund, 2016; Folatelli et al., 2016), Ib 2019yvr (Kilpatrick et al., 2021; Sun et al., 2022), Ic 2020oi (Gagliano et al., 2022), although none of these have been confirmed with post-explosion imaging.

Broad-line SNe Ic (SNe Ic-BL) are a rare subclass of the Type Ic events, exhibiting broad spectral features that indicate unusually high ejecta velocities of $\sim 20{,}000-30{,}000$  km s^-1($\sim 0.1c$, Galama et al. 1998; Modjaz et al. 2006, 2016; Sahu et al. 2018; Taddia et al. 2019). Additionally, the luminosities of these SNe are typically higher than those of other core-collapse SNe (Drout et al., 2011; Cano, 2013; Taddia et al., 2019). Models suggest that the kinetic energy of these explosions can be as high as 10⁵² ergs (Mazzali et al., 2002; Maeda et al., 2003; Janka et al., 2016; Prentice et al., 2016), which is a factor of 10 larger than other typical SNe. Furthermore, SNe Ic-BL are the only SN type to be associated with long-duration $\gamma$-ray bursts (LGRBs), whose gamma-ray emission lasts longer than $2$ seconds (Galama et al., 1998; Woosley & Bloom, 2006a; Modjaz et al., 2016, 2019). Iwamoto et al. (1999) showed these SNe can also be accompanied by $X-ray$ flashes (XRF), which further suggests that only SNe Ic-BL are associated with LGRBs.

LGRBs are considered to be a natural result of the collapsar model for SESNe (Woosley & MacFadyen, 1999). In this progenitor scenario, the core of a massive star collapses to create a rapidly rotating compact object, then accretion onto this compact object utilizes the rotational energy of the star via magnetic coupling and launches collimated jets that power the explosion and produce the GRB (Woosley, 1993; Woosley & MacFadyen, 1999). The absence of LGRBs and/or their associated afterglows in normal SNe Ib/c disfavors the collapsar theory for explaining these SNe, and challenges the connection between normal SNe Ic and SNe Ic-BL (see, Corsi et al., 2012; Barnes et al., 2018; Modjaz et al., 2019; Zenati et al., 2020). Many SNe Ic-BL, however, are observed without an associated LGRB, and the reason for this is not fully understood.

Other progenitor channels for SNe Ic-BL may be possible. Zapartas et al. (2017) provide a detailed analysis of deep, post-explosion upper-limits of the Type Ic-BL SN 2002ap. The SN progenitor has a low metallicity and high ejecta mass, ruling out most single star models. Instead, the deep upper limits, when combined with the explosion parameters of the primary star, suggest a more likely scenario consisting of a low-mass binary system in an eccentric orbit undergoing non-conservative mass transfer. Another more exotic, but tempting, possibility for SN 2002ap is a reverse merger of the companion star with the compact remnant from the primary explosion.

The results for SN 2002ap suggest a possible bimodality in the progenitor population of SNe Ic-BL, but one important avenue for understanding the progenitors of SESNe and the connection between SNe Ic-BL and LGRBs is through detailed analysis of their environments. Theory predicts that metallicity should play an important role in the production of LGRBs. Specifically, Woosley & Bloom (2006b) suggested that angular momentum loss due to the metal line-driven winds in a WR star could prohibit the formation of a compact object that was rotating fast enough to launch a jet. Therefore, they proposed a metallicity threshold of $Z<0.3Z_{\odot}$ for collapsars.

Several studies have confirmed that SNe Ic-BL prefer lower-metallicity environments relative to other core-collapse SNe (Sanders et al., 2012; Kelly & Kirshner, 2012; Arcavi, 2018; Modjaz et al., 2020). In particular, Modjaz et al. (2020) found that SNe Ic-BL with and without associated GRBs prefer statistically similar environments with low-metallicity and high-specific star formation rates (sSFRs). Modjaz et al. (2020) conclude that SNe Ic-BL without GRBs in their sample either produced jets that were choked within the star (Milisavljevic et al., 2015; Modjaz et al., 2016), or produced off-axis GRBs. Nonetheless, it is still debated whether metallicity is required to explain the presence of jet (Mannucci et al., 2011).

The clear connection between SESNe type and their host environment suggests that metallicity could be a key factor in discerning their progenitor scenarios. Observing SESNe at high-$z$ presents an opportunity to study their properties in a diversity of environments. At high-$z$ there is more access to very low metallicity environments ($Z<0.1Z_{\odot}$) like those present in local metal-poor dwarf galaxies. So far only a small sample of low-$z$ core collapse SNe have been studied in extreme environments like these (Arcavi et al., 2010; Taddia et al., 2016; Anderson et al., 2018; Gutiérrez et al., 2018; Tucker et al., 2024). Given the rarity of SNe Ic-BL, we do not yet know if a population of these SNe exists in very-low metallicity environments. If the environments of SNe Ic-BL deviate from the LGRB population at high-$z$, it may indicate that multiple progenitor scenarios are required to explain these SNe. Furthermore, we can constrain how their observational properties and rates may evolve with redshift.

In this work, we present James Webb Space Telescope (JWST) NIRCam and NIRSpec observations of SN 2023adta. This SN was identified in the JWST Advanced Deep Extragalactic Survey (JADES, Eisenstein et al. 2023). A JWST Director’s Discretionary Time (DDT) program was accepted to follow-up the most interesting transients ($1<z<4$) in this field with additional NIRCam photometry and NIRSpec spectroscopy (DeCoursey et al. 2024, hereafter D24). Details on the sample of high-$z$ SNe and the design of the JADES Transient Survey + DDT observations are presented in this companion paper.

In combination with the resulting light curve from these data, we used the NIRSpec Prism spectrum of SN 2023adta (covering $0.6\mu$m $-$ $5.3\mu$m at $R\sim 100$) to confirm its redshift ($z=2.83$), and classify it as a SN Ic-BL. To date, this is the highest-redshift SESN ever discovered. We also determine that it came from a low-metallicity environment which is consistent with low-$z$ analogs.

In Section 2, we summarize our observations, in Section 3 we describe the observational properties of SN 2023adta (and its host galaxy), we discuss its classification as a SN Ic-BL, and put it into the context of other local SESNe, and in Section 4 we summarize our findings, their implications, and discuss future observations. In all subsequent analysis we assume a standard $\Lambda$CDM cosmology with $H_{0}=70$km s^-1 Mpc^-1, $\Omega_{m}=0.315$.

2 Observations & Data Reduction

A detailed description of the JADES transient observing program methods and data reduction are presented in D24. In short, JADES observations were taken acquired over a 1 year time baseline from September 2022 to January 2024, allowing us to conduct a transient search of unprecedented $5\sigma$ m${}_{AB}\sim 30$ depth. The first observing window $2022$ September $29$-October $5$, and the second epoch took place between $2023$ September $29$-October $3$ with an overlap of $25\hbox{${}^{\prime}$}^{2}$ in the NIRCam F090W, F115W, F150W, F200W, F277W, F335M, F356W, F410M, F444W filters. Additional visits on $2023$ November $15$ and $2024$ January $1$ were performed due to failed observations.

A JWST DDT program was approved to follow and classify the most interesting transients with two additional NIRCam visits on $2023$ November $28$ and $2024$ January $1$. The second visit including the NIRSpec multi-object spectroscopy (MOS) mode using the micro-shutter assembly (MSA) and Prism (R$\sim 100$) grating. The MSA provided SN spectra for $\sim 10$ transients, some of which are described in companion papers (e.g., Pierel et al. 2024, Coulter et al., Egami et al. in preparation), as well as a variety of galaxy spectra. In this section, we describe the data products and reduction methods for SN 2023adta.

2.1 JWST NIRCam

Our data reduction methods are described in detail by D24, but we summarize the process here. We adopt the point-spread function (PSF) fitting method developed in Pierel et al. (2024) for measuring photometry on Level 3 (drizzled and resampled 2D data) JWST images. Unlike their scenario though, we have a template image for all epochs of SN 2023adta from the $2022$ JADES observations. We therefore first align the Level 2 (CAL) NIRCam images containing SN 2023adta to the template images (in each filter) using the JWST/HST Alignment Tool (JHAT; Rest et al., 2023) software and then produce aligned Level 3 images with the JWST pipeline (v$1.12.5$; Bushouse et al., 2024, and see Figure 1). We obtain difference images in all filters using the High Order Transform of PSF and Template Subtraction (HOTPANTS; Becker, 2015) code, and then implement the Level 3 PSF fitting routine from Pierel et al. (2024) using $5\times 5$ pixel cutouts and Level 2 PSF models from webbpsf, which are temporally and spatially dependent and include a correction to the infinite aperture flux. The Level 2 PSF models are drizzled to create a Level 3 PSF model consistent with the observations. The total measured fluxes, which are in units of MJy/sr, are converted to AB magnitudes using the native pixel scale of each image ($0.03\hbox{${}^{\prime\prime}$}/$pix for SW, $0.06\hbox{${}^{\prime\prime}$}/$pix for LW). Measured photometry is given in Table 1.

Figure 1: False color images of SN 2023adta and its host galaxy taken in 2022 (left) and 2023 (right). The white rectangles show the orientation of the MSA 3-shutter slitlet on sky. Two additional nodded exposures were acquired north and south of this central location. The SN (indicated by the red arrow) is located near the southern edge of the central slit.

Table 1: Photometry for SN 2023adta measured in Section 2.1.

PID	Instrument	MJD	Filter/Disperser	mAB
$1180$	NIRCam	$60220$	F090W	$>30.2$
$1180$	NIRCam	$60220$	F115W	$29.65\pm 0.11$
$1180$	NIRCam	$60220$	F150W	$29.11\pm 0.09$
$1180$	NIRCam	$60220$	F200W	$28.54\pm 0.08$
$1180$	NIRCam	$60220$	F277W	$28.91\pm 0.10$
$1180$	NIRCam	$60220$	F335M	$28.60\pm 0.13$
$1180$	NIRCam	$60220$	F356W	$28.69\pm 0.10$
$1180$	NIRCam	$60220$	F410M	$29.19\pm 0.22$
$1180$	NIRCam	$60220$	F444W	$29.31\pm 0.22$
$1180$	NIRCam	$60264$	F090W	$>29.9$
$1180$	NIRCam	$60264$	F115W	$29.73\pm 0.19$
$1180$	NIRCam	$60264$	F150W	$27.45\pm 0.06$
$1180$	NIRCam	$60264$	F200W	$26.53\pm 0.04$
$1180$	NIRCam	$60264$	F277W	$26.64\pm 0.04$
$1180$	NIRCam	$60264$	F335M	$26.52\pm 0.04$
$1180$	NIRCam	$60264$	F356W	$26.74\pm 0.04$
$1180$	NIRCam	$60264$	F410M	$26.86\pm 0.06$
$1180$	NIRCam	$60264$	F444W	$27.07\pm 0.06$
$6541$	NIRCam	$60276$	F115W	$>28.9$
$6541$	NIRCam	$60276$	F150W	$27.63\pm 0.09$
$6541$	NIRCam	$60276$	F200W	$26.55\pm 0.05$
$6541$	NIRCam	$60276$	F277W	$26.57\pm 0.05$
$6541$	NIRCam	$60276$	F356W	$26.59\pm 0.06$
$6541$	NIRCam	$60276$	F444W	$26.92\pm 0.08$
$6541$	NIRCam	$60310$	F150W	$28.44\pm 0.14$
$6541$	NIRCam	$60310$	F200W	$26.98\pm 0.05$
$6541$	NIRCam	$60310$	F277W	$26.88\pm 0.07$
$6541$	NIRCam	$60310$	F356W	$26.69\pm 0.06$
$6541$	NIRCam	$60310$	F444W	$26.94\pm 0.09$
$6541$	NIRSpec	$60310$	Prism	–
$1180$	NIRCam	$60311$	F090W	$>29.4$
$1180$	NIRCam	$60311$	F115W	$>29.9$
$1180$	NIRCam	$60311$	F150W	$28.85\pm 0.10$
$1180$	NIRCam	$60311$	F200W	$27.01\pm 0.04$
$1180$	NIRCam	$60311$	F277W	$26.79\pm 0.04$
$1180$	NIRCam	$60311$	F335M	$26.42\pm 0.04$
$1180$	NIRCam	$60311$	F356W	$26.63\pm 0.04$
$1180$	NIRCam	$60311$	F410M	$26.57\pm 0.05$
$1180$	NIRCam	$60311$	F444W	$26.79\pm 0.05$

2.2 JWST Spectroscopy

Figure 2: (top): The 2D, master-background subtracted, NIRSpec Prism spectrum of SN 2023adta and its host galaxy in MJy/steradian. As shown in Figure 1, our observations capture both the SN and part of its host galaxy. The central rows for the extraction of SN 2023adta and its host galaxy are noted by the dashed-blue and pink horizontal lines, respectively. (middle): The resulting 1D extractions of SN 2023adta (blue) and its host galaxy (pink) in flux converted to $f_{\lambda}$. The spectrum of SN 2023adta was scaled to its photometry at the same epoch, and the spectrum of the host galaxy was scaled to the 2022 photometry at the position of the SN with an identical $0.1^{\prime\prime}$ aperture. Notable narrow emission lines from the host galaxy are marked with the dashed-vertical lines. (bottom): The difference between the flux-calibrated SN and host galaxy spectra from the middle panel. Assuming the host galaxy emission underlying the SN is similar to the galaxy emission at our extracted location, this spectrum should be mostly isolated SN flux. However, the presence of narrow emission lines could indicate that this subtraction is imperfect.

Table 2: JWST SN 2023adta NIRSpec Observation Details

Instrument	NIRSpec
Mode	MOS
Wavelength Range	0.6 $-$ 5.3$\mu$m
Slit	3 Shutter ($0.46^{\prime\prime}\times 0.2^{\prime\prime}$ each)
Grating/Filter	Prism/CLEAR
R $=\lambda/\Delta\lambda$	$\sim 30-300$
Readout Pattern	NRSIRS2
Groups per Integration	19
Integrations per Exposure	2
Exposures/Nods	3
Total Exposure Time	16,631s

SN 2023adta (R.A. $=3$h$32$m$32.4657$s decl. $=-27$d$48$m$52.2371$s) was selected as one of the highest-priority targets for spectroscopic follow-up observations because of its brightening in Nov 28, 2023 images (F200W $=26.5$ mag) relatively to the discovery image taken on Oct 1, 2022 (F200W $=28.5$ mag). We observed SN 2023adta on 2024 January 1 with the JWST micro-shutter assembly (MSA, Ferruit et al. 2022) and NIRSpec Prism. Three nodded exposures were acquired using 3-shutter slitlets centered on each target in the MSA design. This observing pattern results in five unique shutter locations (each 0.46”$\times$ 0.20”) open to sky around each target. We show the orientation of the central 3-shutter slitlet for SN 2023adta in Figure 1. This orientation allowed for the independent extractions of SN 2023adta and part of its extended host galaxy.

We reduced the JWST data using the “jwst” pipeline (version 1.14.0; Bushouse et al., 2024) routines for bias and dark subtraction, background subtraction, flat-field correction, wavelength calibration, flux calibration, rectification, outlier detection, and resampling. Given that the host galaxy is extended and its flux falls within multiple shutters, the default nodded point source background subtraction method is not sufficient for this target. Therefore we re-ran Stage 2 of the pipeline (calwebb_spec2) using the master background subtraction strategy which makes use of designated background shutters. We performed manual interactive extractions of both SN 2023adta and its host galaxy from the 2D Stage 2 data product using the Specviz module of Jdaviz. Since the resolution of the Prism is non-linear with wavelength ($R\sim 30-300$, Jakobsen et al. 2022), we use an irregularly spaced wavelength grid with $\Delta\lambda$ ranging $34-201$ Å.

Some details of our data reduction choices and results are shown in Figure 2. In the top panel, we show the 2D, master-background subtracted NIRspec Prism spectrum of SN 2023adta and its host galaxy. The colormap has units of MJy/steradian. The trace of SN 2023adta is shown by the dashed-blue horizontal line, and the trace of the host galaxy trace is indicated by the dashed-pink horizontal line and is located 0.4” south of SN 2023adta. Using the Specviz module of JDaviz (Jdadf Developers et al., 2023), we perform a boxcar extraction with 3-pixel width on both of these traces.

We flux calibrate the host galaxy spectrum in two different ways. First, we scale the the spectrum to the 2022 host galaxy photometry (well before the SN explosion; Eisenstein et al., 2023). We use this spectrum as input for the SED fitting code BAGPIPES (Carnall et al., 2018, 2019) to derive an initial redshift, and the host galaxy parameters presented in Section 3.4. Second, we scale the original extracted host spectrum to the 2022 host galaxy photometry measured at the SN position in order to get an estimate of the underlying host galaxy flux. We do similarly for the extracted spectrum of SN 2023adta but instead scale to the NIRCam photometry observed at the same epoch. The flux-calibrated spectra of SN 2023adta and its host-galaxy (SN-position estimate) are shown in the middle panel of Figure 2 (blue and pink curves, respectively), and the difference between them is shown in the bottom panel (black).

Clear narrow host-galaxy emission from [O II], [O III], H$\alpha$, and [S III], are present in the spectrum of SN 2023adta (marked with dashed-vertical lines). The host galaxy exhibits a strong Balmer break (1.5$\mu$m observer frame) and has notably weaker emission lines relative to the continuum than what is present in the spectrum of SN 2023adta. We expect that our measured host galaxy continuum flux (determined from the extended emission in the southern-most shutter) should be a reasonable estimate for what is underlying the SN. However, we note that the SN is $0.1$” from a region of star formation whose flux also contributes to our SN spectrum (see Figure 1). This could be the reason for the imperfect subtraction of the host galaxy emission lines.

The difference spectrum should contain mostly isolated SN flux. In this final spectrum of SN 2023adta, we see broad absorption features at $3.0\mu$m and $3.7\mu$m. We analyze these features in more detail in Section 3.2.

2.3 Redshift and Host Galaxy

SN 2023adta was identified in host galaxy JADES-GS+53.13533-27.81457. We fit the NIRspec Prism data of the host galaxy using the SED fitting code BAGPIPES (Carnall et al., 2018, 2019) to determine an initial redshift estimate ($z=2.83$). We then further refine this estimate by fitting individual emission lines present in the spectrum following the methods of Bunker et al. (2023). At this redshift, several galaxy emission lines blended together. We model H$\alpha+$N II, and [S II] $\lambda\lambda$ $6716,6731$ as individual Gaussian components. We fit H$\alpha+$N II and [S II] $\lambda\lambda$ $6716,6731$, simultaneously with their centroids fixed relative to one another, and do similarly for [O III] $\lambda\lambda$ $4959,5007$. Using these gas-phase nebular emission lines we measure a final redshift of $z=2.830\pm 0.001$. This uncertainty is consistent with those measured for Prism spectra in Bunker et al. (2023) ($\Delta z=0.001-0.01$ at $z=2.83$). For a detailed analysis of the properties of the host galaxy of SN 2023adta see Section 3.4

3 Analysis

3.1 Light Curve

Table 3: Best-fit light curve parameters and properties for SN 2023adta

Parameter	Bounds	SN Ib/c	SN II	SN Ia
$z$	Fixed	$z=2.83$	$z=2.83$	$z=2.83$
$t_{pk}$	[60120,60411]	$60270.09\pm 0.32$	$60255.75\pm 0.51$	$60278.45\pm 0.41$
$A$	[0,$6\times 10^{-18}$]	$(4.86\pm.05)\times 10^{-19}$	$(2.73\pm.03)\times 10^{-19}$	…
$M_{V}$	–	$19.0\pm 0.2$†	$19.4\pm 0.3$†	$19.03\pm 0.03$
Host $E(B-V)$	[0,3]	$(1.17\pm 4)\times 10^{-3}$	$0.43\pm 0.0006$	…
Host $R_{V}$	[2,4]	$2.92\pm 0.08$	$2.03\pm 0.05$	…
$x_{0}$	[0,$3\times 10^{-6}$]	…	…	$(5.92\pm 0.05)\times 10^{-8}$
$x_{1}$	[-3,3]	…	…	$-0.29\pm 0.08$
$c$	[-1.5,1.5]	…	…	$0.63\pm 0.02$
$\chi^{2}/\nu$	…	1.00 ($1.6\pm 0.4$)†	1.56 ($2.4\pm 0.5$)†	1.75

Note. — Upper limits are 5$\sigma$.

Figure 3: (left): The multiband light curve of SN 2023adta (colored points), plotted with the offsets shown for visual clarity, and its best-matched SN Ib/c light curve (dashed-curves). The epoch of our NIRSpec Prism spectrum is marked with the dashed-vertical line, and indicates an approximate phase of +11 days relative to peak brightness in F150W. (middle): Similar to the left panel but showing the best-fit SN II. (right): Similar to the other panels but showing the best-fit SN Ia SALT3-NIR model (Pierel et al., 2022), with shape and color parameters of $x_{1}=-0.29$ and $c=0.63$, respectively.

We examine the light curve of SN 2023adta (originally presented in D24) in Figure 3. The data shown (colored points) are NIRCam wideband photometry. We fit the observed photometry using all existing core-collapse, SN light curve evolution models with rest-frame optical to near-IR (to observer-frame $\sim 4\mu$m) wavelength coverage (Pierel et al., 2018), and the SALT3-NIR SN Ia light curve model (Pierel et al., 2022). Prior to fitting, we correct for Galactic extinction ($E(B-V)=0.01$ mag with $R_{V}=3.1$). In Figure 3, we show the best-fit SN Ib/c model, SN II model, and SALT3-NIR SN Ia model from left to right. The epoch of our NIRSpec DDT observation is represented by the vertical-dashed line. For the SN Ib/c model, this occurs at a phase of $+11$ days relative to peak-brightness in F150W ($+14$d and $+8$d for the SN II and SN Ia model, respectively). For the SN Ia model, we find best-fit shape and color parameters of $x_{1}=-0.29$ and $c=0.63$, respectively. The core-collapse light curve models do not have uncertainties. Therefore, to determine relative goodness of fit, we add a constant uncertainty to each model light curve such that $\chi^{2}/\nu=1.0$ for the best-fit (SN Ib/c). We add identical uncertainty to the SN II and SALT3-NIR SN Ia models so that the $\chi^{2}/\nu$ values are directly comparable. This uncertainty is shown by the shaded regions in Figure 3. The model parameters for each fit are shown in Table 3.

We find that SESNe reproduce the multiband light curve evolution of SN 2023adta. Specifically, from each of the top 10 core-collapse SN Ib/c and SN II fits, we determine $\chi^{2}/\nu=1.6\pm 0.4$ and $2.4\pm 0.5$, for each SN type, respectively. The best-fit SN Ia model has $\chi^{2}/\nu=1.75$. While all fits perform well in the rest-frame optical (F150W, F200W, F277W, and F356W), the best-fit SN Ia is unable to reproduce the flux at rest-frame NUV wavelengths (F115W). The best-fit SN Ib/c and SN II are also better matched to the last epoch in F444W. Due to limitations of the data, we caution that the light curves models are most uncertain at the reddest wavelengths. If we assume $H_{0}=70$ km s^-1 Mpc^-1 (Riess et al., 2016, 2018), we find a high peak luminosity of $M_{V}$ = $-19.0\pm 0.2$ mag ($M_{R}$ = $-19.1\pm 0.2$ mag) for the SN Ib/c model. We note that the best-fit SN II models similarly require high-peak brightnesses ($M_{V}$ = $-19.4\pm 0.3$ mag) that are difficult to reconcile with the observed luminosity functions of normal-SN II (Li et al., 2011; Valenti et al., 2016). Similar to $\chi^{2}/\nu$, we derive these uncertainties from the distribution of the 10 best-fit SN Ib/c models.

Some studies have shown that SN Ic-BL light curves are very similar to those of other SNe Ib/c but tend toward higher absolute luminosities (Drout et al., 2011; Taddia et al., 2015; Prentice et al., 2016; Lyman et al., 2016). In particular, Drout et al. (2011) found a mean $M_{R,peak}=-19.0\pm 1.1$ mag for SNe Ic BL ($M_{R,peak}=-17.9\pm 0.9$ mag and $-18.3\pm 0.6$ mag for SNe Ib and SNe Ic, respectively). In Figure 4, we show a comparison of the SN 2023adta peak brightness ($M_{R}$) and light curve $\Delta m_{15,R}$ shape, to other SESNe from Drout et al. (2011) and SNe Ic-BL from Taddia et al. (2019). These measurements were derived from $R$- and $r$-band photometry, respectively. Given its sparsely-sampled light-curve, $\Delta m_{15,R}$ is uncertain but consistent with other SESNe. Its peak brightness is similar to those of other SN Ic-BL, in particular, those with associated GRBs (blue diamonds with red edges). Using the relationship derived in Drout et al. (2011):

$$log(M_{Ni})\approx-0.41M_{R}-8.3\\$$ (1)

we estimate that $M_{Ni}=0.3\pm 0.1$ $M_{\odot}$ for SN 2023adta.

Figure 4: Comparison of the light curve peak brightness ($M_{R,r}$) and shape ($\Delta m_{15,R,r}$) characteristics of SN 2023adta to other SESNe from Drout et al. (2011) and SNe Ic-BL from Taddia et al. (2019). Measurements were derived from $R$- and $r$-band photometry, respectively.

3.2 Spectroscopic Classification & Comparisons

Figure 5: Rest-frame NIRSpec Prism spectrum of SN 2023adta (black curves) compared with spectra of a variety of other SN types at a similar phases (colored curves). We mask negative values that result from imperfect host-galaxy subtraction. The dotted-gray curves show clipped host galaxy emission lines. The color of each other curve corresponds to the SN type. From top to bottom we show: five SNe Ic-BL (blue), one SN Ic (light-blue), two SNe Ib (green), one SN Ia (yellow), one Ca-rich SN (orange), and one SN IIb (dark red). A “*” next to the phase of a comparison spectrum indicates that it was a matched spectrum using the SN classification code NGSF (Goldwasser et al., 2022).

Table 4: Best-fit NGSF parameters for SN 2023adta . Below the double-line we show the highest rankings for specific SN types.

Rank	SN	Type	Phase (days)	$A_{V}$ (mag)	$\chi^{2}/\nu$
1	SN 2006aj	Ic-BL	+16	0.45	3.19
2	SN 2002ap	Ic-BL	+8	1.25	3.42
3	SN 1998bw	Ic-BL	+20	1.90	3.52
4	PTF10vgv	Ic-BL	+2	2.0	3.93
5	SN 2014G	II	+29	1.2	4.12
8	SN 1994L	Ic	+2	2.0	4.37
10	SN 2013df	IIb	+14	1.0	4.45
11	SN 2015N	Ia	+25	0.0	4.57
21	SN 2006ep	Ib	+11	0.0	5.18

Figure 6: Similar to Figure 5 but zoomed in on the O I and Ca II NIR triplet region of each SN Ic-BL. The velocity width of the primary absorption feature suggests ejecta moving at $>20{,}000$  km s^-1. The velocity we measure from the absorption minimum is $29{,}000\pm 2{,}000$  km s^-1.

The low resolution of the NIRSpec Prism and uncertainty in underlying contamination from host galaxy light provide additional challenges in the classification of SN 2023adta. The presence of strong-host galaxy $H{\alpha}$ emission with FWHM = $4{,}000$  km s^-1could bias classification against SNe II, especially those with relatively narrow emission features like SNe IIn.

The spectroscopic classification of SN 2023adta was performed with Next Generation SuperFit (Goldwasser et al., 2022) based on the IDL code Superfit (Howell et al., 2005). This code allows for the simultaneous fitting of the host galaxy light contamination and SN flux. We restrict the fit to use only post-maximum-light template SNe and the “SB1” host galaxy template which best-matches the continuum of our host galaxy spectrum (see Figure 2). We also restrict the extinction parameter to $0<A_{V}<2$ mag. We de-redshift the SN spectrum, clip prominent host-galaxy emission lines, and then perform NGSF fitting at a resolution of $30$ Å. This is similar to the native wavelength binning output by the JWST Stage 3 pipeline which ranges from $25-52$ Å across our de-redshifted spectrum of SN 2023adta.

A summary of the spectroscopic fitting results is presented in Table 4. Ordered by $\chi^{2}/\nu$, the four best-fit SN spectra are from SNe Ic-BL with phases ranging from $+2$ to $+20$ days, in good agreement with the results of our light curve fit ($+11$d). Interestingly, two of these SNe (SN 1998bw and SN 2006aj) have associated GRBs. In each of these fits, the SN contributions are $>98\%$ indicating little host contamination. For SNe Ic-BL, NGSF prefers a high amount of reddening ($A_{V}>0.45$ mag), which is consistent with our derived host galaxy parameters ($A_{V}=1.8$ mag, see Section 3.4), but not our best-fit light curve ($A_{V}=0$ mag). However, we caution that imperfect subtraction of the host galaxy light could strongly influence this parameter. Similar to the light curve fitting, the second best fit SN type is SN II, however, its phase ($+29$d) is significantly greater than what we would expect based on best-fit SN II light curve ($+13$d).

Given the limited SN template bank available in NGSF, additional comparisons are needed to provide a secure classification for SN 2023adta. We show spectroscopic matches to SN 2023adta from NGSF along with other useful comparison spectra in Figure 5. For log$(f_{\lambda})$ comparisons, we mask negative values in SN 2023adta that result from imperfect host-galaxy subtraction.

From top to bottom, we show our NIRspec Prism spectrum of SN 2023adta in comparison to five SNe Ic-BL (SN 2014ad, Kwok et al. 2022; Shahbandeh et al. 2022; SN 2002ap, Mazzali et al. 2002; SN 2013dk, Elias-Rosa et al. 2013; Shahbandeh et al. 2022; SN 1998bw, Galama et al. 1998; and SN 2006aj, Pian et al. 2006), one SN Ic (SN 2013ge, Drout et al. 2016; Shahbandeh et al. 2022), two SNe Ib (SN 2006ep, Smartt et al. 2015; and SN 2013ek, Smartt et al. 2015; Shahbandeh et al. 2022), one SN Ia (SN 2015N, Stahl et al. 2020), one Ca-rich SN (SN 2018byg, De et al. 2019), and one SN IIb (SN 2013ak, Smartt et al. 2015; Shahbandeh et al. 2022). Four of these SNe (SN 2014ad, SN 2013ge, SN 2013ek, and SN 2013ak) have optical spectra that were observed within one day of a NIR spectrum presented in Shahbandeh et al. (2022). In these cases, we have combined these data into one full optical-to-NIR spectrum. All comparison spectra have been re-binned to a have a similar wavelength spacing to the NIRSpec spectrum of SN 2023adta. We mark spectroscopic fits from NGSF with a “*”.

There are two primary broad absorption features in the spectrum of SN 2023adta that agree relatively well with other SNe Ic-BL, one at $\sim 8{,}000$ Å and another at $\sim 10{,}000$ Å. The former is likely produced by the Ca II NIR triplet, and the latter is likely produced from blended contributions from C I, He I, and Mg II (Shahbandeh et al., 2022). Absorption from O I and Ca II can appear as separate features in normal SNe Ic and SNe Ib due to their lower velocities (Matheson et al., 2001). However, in SNe Ic-BL, these features can blend together due to the larger ejecta velocities (Foley et al., 2003). SN Ic-BL comparison spectra in Figure 5 have been selected to showcase the diversity in this wavelength region. Additionally, the spectrum of SN 2013dk at +6 days (+12 in the NIR) exhibits some interesting similarity to SN 2023adta at $\sim 9{,}000$ Å. Shahbandeh et al. (2022) identified this absorption feature as C I in SN 2013dk. This SN also matches the fall-off in flux near $4{,}000$ Å better than most of the other comparison SNe.

We also make spectroscopic comparisons to normal SESNe and a SN Ia. Generally, these SNe are unable to reproduce the large blueshifts seen in SN 2023adta and the line strengths relative the continuum of these SNe are not well matched. The presence of He I and/or Si II in the spectrum of SN 2023adta is ambiguous. SN 2023adta has additional shallow absorption features near $6{,}000$ Å that could result from these elements but the resolution of the data, and uncertainty on phase make a definitive line identification in this wavelength region very difficult. Furthermore, He I $\lambda 10830$ Å could contribute to the broad absorption at $10{,}000$ Å, however, there is no way to differentiate from potential contributions from C I $\lambda 10693$ Å  and/or Mg II $\lambda 10927$ Å (Shahbandeh et al., 2022). Furthermore, For the SN Ia (yellow curve, SN 2015N), NGSF prefers a later phase of $+25$ days. This could be due to the color of SN 2023adta at the epoch of our spectrum. A SN Ia at an earlier phase as suggested by the light curve fit, would likely be too blue to reproduce the observed continuum.

Some peculiar transients like Ca-rich SNe are red at early times due strong line-blanketing (Shen et al., 2018; Polin et al., 2019; Jacobson-Galán et al., 2020). In Figure 5, we show one such event, SN 2018byg (orange curve, De et al. 2019), at $-5$ days relative to maximum-light. This spectrum reproduces the continuum of SN 2023adta well and also has a high Ca II velocity of $-22{,}000$ km s^-1 (De et al., 2019; Jacobson-Galán et al., 2020), however, the strength of this absorption does not match as well as the Ca II features in our SNe Ic-BL comparison spectra. Finally, we compare to a SN IIb and a SN II (dark red curves, SN 2013ak and SN 2014G). These SNe have low ejecta velocities that are difficult to reconcile with SN 2023adta. We measure a H$\alpha$ FWHM$=4{,}000$  km s^-1which is consistent with the resolution of the NIRSpec Prism at this wavelength indicating that this emission is from the host galaxy. However, a clear continuum is difficult to define in this region and broader H$\alpha$ may be present.

Even though there are significant caveats related to both the photometric and spectroscopic classification methods, the observational characteristics of SN 2023adta outlined below provide convincing evidence that it can be classified as a SN Ic-BL at z=2.83:

• Light curve fit strongly prefers SN Ib/c which reproduces the F115W flux at early epochs.

• NGSF fitting of the NIRSpec spectrum strongly suggests the subclass of SN Ic-BL at a phase consistent with the light curve.

• Broad, high-velocity absorption features are inconsistent with those observed in normal SNe Ia, Ib, and Ic.

• Lack of obvious H or He absorption/emission features.

• Red continuum shape and flux drop-off at rest-frame $4000$ Å.

3.3 Ejecta Velocity and Physical Parameters

Table 5: Estimated explosion parameters for SN 2023adta .

$v_{ph}$ ( km s-1)	$M_{Ni}$ $(M_{\odot})$	$M_{ej}$ $(M_{\odot})$	$E_{KE}$ ($10^{51}$ erg)
$29{,}000\pm 2{,}000$	$0.3\pm 0.1$	$5^{+4}_{-2}$	$24^{+29}_{-14}$
$20{,}000$	$0.3\pm 0.1$	$3^{+3}_{-2}$	$8^{+11}_{-5}$
$10{,}000$	$0.3\pm 0.1$	$1.6^{+1.8}_{-0.9}$	$1.0^{+2.0}_{-0.7}$

Using a combination of both photometric and spectroscopic properties of SN 2023adta, we can derive some physical parameters of the explosion. Drout et al. (2011) showed how peak luminosity ($M_{R}$), light curve shape ($\Delta m_{15,R}$), and photospheric velocity ($v_{ph}$) are related to the nickel mass ($M_{Ni}$), ejecta mass ($M_{ej}$), and kinetic energy of the explosion ($E_{KE}$). To get an estimate of $v_{ph}$ we examine the Ca II and O I absorption features of SNe Ic-BL in more detail in Figure 6. The spectrum of SN 2014ad reproduces the shape of the primary absorption features well but at larger blueshifts. Kwok et al. (2022) attribute the absorption at $\sim 7{,}000$ Å primarily to O I and estimate a photospheric velocity of $22{,}000$ km s^-1 via a best-fit TARDIS (Kerzendorf et al., 2014) model.

Assuming the absorption minimum in the spectrum of SN 2023adta is primarily caused by Ca II, using an intensity-weighted rest wavelength of 8567 Å for the Ca II NIR triplet, we estimate an ejecta velocity of $-29{,}000\pm 2{,}000$ km s^-1 which is consistent with other SN Ic-BL (Modjaz et al., 2016; Taddia et al., 2019). We find that SN 2023adta may have another weaker absorption feature blueward of Ca II, which if caused by O I, also indicates ejecta with velocities ranging from $20{,}000$ to $40{,}000$  km s^-1. An estimate of the Fe II $5169$ Å line velocity would be a better indicator of the photospheric velocity since this line does not saturate (Branch et al., 2002; Modjaz et al., 2016), however, the strong residual [O III] emission from the host galaxy prohibits our ability to make this measurement.

Using the properties of the light curve that we estimated in Section 3.1, we derive physical explosion parameters for SN 2023adta using the equations 1 and 2 and methods described in Drout et al. (2011). We determine these parameters using our measured Ca II velocity as the photospheric velocity, but also calculate results for fiducial velocities of $20{,}000$  km s^-1and $10{,}000$  km s^-1for comparison. These parameters are shown in Table 5. We estimate a large ejecta mass and kinetic energy of $5^{+4}_{-2}$ $M_{\odot}$ and $24^{+29}_{-14}\times 10^{51}$ ergs, respectively. The large uncertainties primarily arise from the poor constraints on the light curve shape and strong dependence ejecta velocity. Nevertheless, these estimates are consistent with what we expect for engine driven explosions (Drout et al., 2011; Cano, 2013; Taddia et al., 2019).

3.4 Host Galaxy Properties

Figure 7: A comparison of the host galaxy properties of SN 2023adta to the mass-metallity relationship observed in SDSS field galaxies (gray contours, $z\sim 0$, Eisenstein et al. 2011), core-collapse SN host galaxies from Kelly & Kirshner (2012) and additional SN Ic/Ic-BL host galaxies from Modjaz et al. (2020) (colored points), and high-$z$ galaxies observed with JWST (pink-yellow hexagons, Morishita et al. 2024). The dashed-yellow line is the empirical, $z=3.3$, mass-metallicity relationship derived in Sanders et al. (2021).

The presence of a SN Ic-BL at $z=2.83$ provides a unique opportunity to study the environment of this rare subclass of core-collapse supernova in the early Universe. Its host galaxy is irregular with several bright knots of star-formation (see Figure 1). We derive galaxy properties by performing a BAGPIPES (Carnall et al., 2018, 2019) fit to the host galaxy spectrum that was flux-calibrated to the total host flux. From this fit we find mass log${}_{10}(M_{*}/M_{\odot})=9.78^{+0.09}_{-0.11}$, metallicity $Z=0.35^{+0.16}_{-0.08}Z_{\odot}$, $A_{V}=1.832^{+0.10}_{-0.14}$, timescale $\tau=2.10^{+0.06}_{-0.05}$ Gyr (SFR$\propto e^{t/\tau}$), and redshift $z=2.833^{+0.002}_{-0.003}$. Assuming a solar metallicity of log(O/H) + 12 = 8.69 (Asplund et al., 2009), we convert the best-fit $Z$ from BAGPIPES to oxygen abundance and find log(O/H) + 12 = $8.23^{+0.16}_{-0.12}$.

In Figure 7, we compare the mass and metallicity of the host galaxy of SN 2023adta to the mass-metallicity relation of low-$z$ galaxies (SDSS DR8, Eisenstein et al. 2011), high-$z$ galaxies ($z\sim 3-9$, Curti et al. 2023; Nakajima et al. 2023; Morishita et al. 2024), core-collapse supernova host galaxies (Kelly & Kirshner, 2012; Modjaz et al., 2020). These data are presented as the gray contours, pink-yellow hexagons, and colored points, respectively. The SDSS and SN host galaxy metallicities were derived via the PP04 O3N2 calibration (Pettini & Pagel, 2004), and metallicities of the high-$z$ galaxies were derived using the direct electron-temperature ($T_{e}$) method. We also compare to the mass-metallicity relationship derived for galaxies at $z\sim 3.3$ in the MOSDEF survey (dashed-yellow line, Sanders et al. 2021).

We find that the metallicity of the host galaxy of SN 2023adta is in good agreement with those of other SN Ic-BL host galaxies (dark blue diamonds). This corroborates previous work that has determined that these high-energy events tend to occur in lower-metallicity environments than those of other SE and normal core-collapse SNe (Sanders et al., 2012; Kelly & Kirshner, 2012; Modjaz et al., 2020). In particular, our measured metallicity for SN 2023adta, is very similar to the mean SN Ic-BL PP04 O3N2 metallicity, log(O/H) + 12 = $8.31\pm 0.04$, and mean SN-GRB PP04 O3N2 metallicity log(O/H) + 12 = $8.20\pm 0.06$, reported in Modjaz et al. (2020). Furthermore, both the mass and metallicity of the host galaxy of SN 2023adta are consistent with what we expect for galaxies at this redshift. Given their rarity ($\sim 4$% the rate of core-collapse SNe by volume, Shivvers et al. 2017), the likelihood of observing a SN Ic-BL in the JADES DDT sample is small. However, since the metallicities of high-$z$ galaxies are lower, it is possible that SNe Ic-BL occur at an enhanced rate at higher redshifts. We note that the mean galaxy mass of the SN Ic-BL host galaxies displayed is log${}_{10}(M_{*}/M_{\odot})=9.1$ with a standard deviation of $0.8$. SN 2023adta has relatively high a host galaxy log${}_{10}(M_{*}/M_{\odot})=9.78^{+0.09}_{-0.11}$. However, SN 2023adta was not discovered in an untargeted survey like the majority of SNe Ic-BL displayed. Since the presence of a clear host galaxy was used to motivate its spectroscopic follow-up, there could be a significant selection effect in the JADES SN sample for SNe in more massive host galaxies.

3.5 Searching for Coincident GRBs with SN 2023adta

The connection between SNe Ic-BL and LGRBs would be interesting to constrain at high-$z$. We used the online database GRBWeb to search for spatially and temporally coincident GRBs that could be associated with SN 2023adta. We find that the best candidate for an associated GRB is GRB230911A. The MJD of this event was 60198, 66 days prior to our first detection of SN 2023adta. Based on our best-fit light curve, this GRB occurred at a phase of $-18$d. The location of this GRB (RA = 59.8, Dec = -34.4, GCN Circular 34652) is 4.3 degrees from SN 2023adta with an uncertainty of 4.1 degrees. However, GRB230911A also has an associated optical afterglow GOTO23akf/AT2023shv (Belkin et al., 2024) at $=3$h$50$m$0.51$s decl. $=-29$d$49$m$30.66$s with 0.3 arcsecond uncertainty. With a separation of 4.3 degrees from SN 2023adta, we conclude GRB230911A is not associated with this SN. Furthermore, a Swift-XRT afterglow was also detected (GCN Circular 34702 ) at R.A. $=3$h$50$m$0.65$s decl. $=-29$d$49$m$33.2$s with 3.4 arcsecond uncertainty, providing additional support for conclusion. We do not find any other GRBs that are likely to be spatially and temporally coincident with SN 2023adta.

4 Discussion & Conclusions

As part of a DDT follow-up program for the JADES Transient Survey (D24), we used JWST NIRCam and NIRSpec observations to confirm the presence of a SN Ic at $z=2.83$ (SN 2023adta). To date, this is the most distant SESN ever discovered, occurring when the age of the universe was approximately 2.2 Gyr. Analysis of the light curve and spectrum of SN 2023adta, provide strong evidence that this SN fits within the SN Ic-BL subclass of SESNe. These SNe are characterized by their broad absorption features indicating typical ejecta velocities of $\sim 20{,}000-30{,}000$ km s^-1, high explosion energies, lack of H and He features, and have been associated with LGRBs. This subclass also has a preference for low-metallicity environments relative to other SESNe (Sanders et al., 2012; Kelly & Kirshner, 2012; Modjaz et al., 2020).

SN 2023adta satisfies many of these criteria. First, we fit the observed light curve and found a best-fit photometric classification of SN Ib/c. In particular, this model reproduces the flux in F115W better than any other SN type. Our best fit indicates that the phase of NIRSpec observation is at $\sim+11$ days. While it difficult to distinguish SN Ib/c/c-BL using solely the light curve shape, we find that SN 2023adta had a peak luminosity in $M_{V}$ = $-19.0\pm 0.2$. This is consistent with the high-luminosities observed in SNe Ic-BL relative to other SESNe (Drout et al., 2011; Cano, 2013; Taddia et al., 2019), and the high peak luminosity of the best-fit SN II models, $M_{V}$ = $-19.4\pm 0.3$, is inconsistent with the luminosity distribution of normal-SN II. Furthermore, engine-driven explosions, the subset of SN Ic-BL characterized by their associations with LGRBs, tend to have the highest luminosities within the SN Ic-BL subclass (Drout et al., 2011; Taddia et al., 2019).

Second, we performed spectroscopic classification of SN 2023adta on our NIRSpec Prism spectrum using the fitting code NGSF based on the IDL code Superfit (Howell et al., 2005). This code models contributions from the SN and host galaxy and found results consistent with our photometric classification. SN 2023adta has broad absorption features in the optical and NIR that are consistent with those observed in SNe Ic-BL. The absorption feature at a rest wavelength of $\sim 8{,}000$ Å is likely produced by the Ca II NIR triplet, and implies ejecta velocities of $\sim 29{,}000$  km s^-1, in agreement with the large explosion energy that we infer from the light curve. Additionally, we do not find strong evidence for H or He features that would be typical for normal SNe Ib or SNe II. In combination with the light curve properties of SN 2023adta, we estimate nickel mass, ejecta mass and kinetic energy ($0.3\pm 0.1$ $M_{\odot}$, $5^{+4}_{-2}$ $M_{\odot}$ and $24^{+29}_{-14}\times 10^{51}$), that are consistent with typical values for SNe Ic-BL (Taddia et al., 2019).

Third, we estimated properties of the host galaxy of SN 2023adta and found that its metallicity is low relative to other core-collapse SNe observed in the local Universe. This is in agreement with previous analyses of SNe Ic-BL host galaxies and what is expected for the galaxy mass-metallicity relationship at $z=2.83$. Lastly, we do not identify any GRBs that are coincident with SN 2023adta.

SN Ic-BL are rare in the local universe ($\sim 4$% of the core-collapse SN rate, Shivvers et al. 2017), therefore, detection of a SN Ic-BL at $z=2.83$ within a small sample suggests that the rate of these events may be enhanced at this distance. Future observations are needed to constrain the rates of these SNe at high-$z$ and to understand how their host populations may or may not evolve with redshift.

JWST is a remarkable tool for this task and in general, the characterization of the high-$z$ transient Universe. The NIRSpec MOS mode is capable of observing several candidate SNe at once, making it very efficient for transient follow-up and classification. The work presented here highlights some additional challenges that the community will face as a sample of high-$z$ SNe continues to be constructed. Specifically, at these distances it can become very difficult to separate light from a SN and its host galaxy. Host galaxy continuum contamination can make SN features appear weak and difficult to distinguish, and at low-resolution (like the NIRSpec Prism), galaxy emission line contamination can preclude the identification of key SN features (e.g., broad H$\alpha$). Future high-$z$ transient survey design should consider these limitations and take steps to mitigate their impact. Furthermore, the Nancy Grace Roman Space Telescope will also include only a low-resolution dispersive element and SN classification will be similarly difficult. More work is needed to improve classification at low-resolution in order to get the most out of future high-$z$ SNe observed with these telescopes.

Acknowledgements

This work is based on observations made with the NASA/ESA/CSA JWST as part of programs 1180 and 6541. This paper is based on observations with the NASA/ESA Hubble Space Telescope and James Webb Space Telescope obtained from the Mikulski Archive for Space Telescopes at STScI. We thank the DDT and JWST/HST scheduling teams and instrument experts at STScI for extraordinary effort in getting the DDT observations used here scheduled quickly.

MRS is supported by an STScI Postdoctoral Fellowship. JDRP is supported by NASA through a Einstein Fellowship grant No. HF2-51541.001 awarded by the Space Telescope Science Institute (STScI), which is operated by the Association of Universities for Research in Astronomy, Inc., for NASA, under contract NAS5-26555. This research is based (in part) on observations made with the NASA/ESA Hubble Space Telescope obtained from the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5–26555. Part of the JWST data used in this paper can be found in MAST: https://dx.doi.org/10.17909/8tdj-8n28 (catalog 10.17909/8tdj-8n28) (JADES DR1).

Additionally, this work made use of the lux supercomputer at UC Santa Cruz which is funded by NSF MRI grant AST 1828315, as well as the High Performance Computing (HPC) resources at the University of Arizona which is funded by the Office of Research Discovery and Innovation (ORDI), Chief Information Officer (CIO), and University Information Technology Services (UITS). AJB acknowledges funding from the “FirstGalaxies” Advanced Grant from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation program (Grant agreement No. 789056).

PAC, EE, DJE, BDJ, are supported by JWST/NIRCam contract to the University of Arizona, NAS5-02015. DJE is also supported as a Simons Investigator. RM acknowledges support by the Science and Technology Facilities Council (STFC), by the ERC through Advanced Grant 695671 “QUENCH”, and by the UKRI Frontier Research grant RISEandFALL. RM also acknowledges funding from a research professorship from the Royal Society.

BER acknowledges support from the NIRCam Science Team contract to the University of Arizona, NAS5-02015, and JWST Program 3215. JDRP is supported by NASA through a Einstein Fellowship grant No. HF2-51541.001 awarded by the Space Telescope Science Institute (STScI), which is operated by the Association of Universities for Research in Astronomy, Inc., for NASA, under contract NAS5-26555. YZ thanks Matteo Cantiello for helpful discussions.