The Disappearance of the Blue and Luminous Progenitor of the Type IIn SN 2010jl

Type IIn supernovae (SNe) exhibit narrow hydrogen lines that arise from the strong interaction between ejecta and circumstellar material. It remains poorly understood, however, what progenitor stars give rise to these explosions. In this work, we perform a detailed analysis of the progenitor and environment of the nearby Type IIn SN 2010jl. With newer images taken by the Hubble Space Telescope, we confirm that the previously reported progenitor candidate is a blend of the progenitor itself and a field star cluster in its close vicinity. SN 2010jl has now become much fainter than the progenitor. The progenitor is very blue and luminous with an effective temperature of log T eff/K = 4.26 −0.09+0.11 and a luminosity of log L/L ⊙ = 6.52 −0.16+0.20 . It is located in a very young star-forming region, but its luminosity is much higher than that expected from the environmental stellar populations. We suggest that the progenitor was in outburst when observed. Its nature and evolutionary history remain to be investigated.


INTRODUCTION
Type IIn supernovae (hereafter SNe IIn) is an interesting and important subclass of core-collapse SNe.They are characterized by strong and narrow (a few hundreds to about 2000 km s −1 ) emission lines, sometimes with broad bases, in the spectra (Schlegel 1990).These features indicate dense circumstellar material (CSM) surrounding the SN, which requires enhanced mass loss (10 −3 -10 0 M ⊙ /years ) of the progenitor stars within a short timescale prior to the deaths.Therefore, SNe IIn are an effective tool to investigate the mass loss of massive stars at their final evolutionary stage.
The progenitor of SNe IIn is a long-standing question.The connection between SNe IIn and luminous blue variables (LBVs) is widely speculated and supported by the direct progenitor detections, such as SN 2005gl, SN 2009ip, SN 2010bt, and SN 2015bh (Elias-Rosa et al. 2016, 2018;Gal-Yam & Leonard 2009;Smith et al. 2010b).The detected progenitors have prominent pre-explosion variabilities and extremely high luminosities that even exceed the Humphrey-Davidson instability limit (Humphreys & Davidson 1994)
In search of SN progenitors, late-time imaging is crucial to confirm whether the detected progenitor candidate has vanished.This is important to distinguish a genuine progenitor from other scenarios, such as a star cluster, blend of multiple objects, or chance alignment.For SNe IIn, this is also important to confirm whether the explosion is indeed the terminal explosion of the progenitor or a non-terminal eruption that is often as extreme as a SN explosion (i.e.SN impostor; Filippenko et al. 1995;Van Dyk et al. 2000, 2002;Kochanek et al. 2011;Van Dyk & Matheson 2012).SN 2010jl is a Type IIn SN in the galaxy UGC 5189A.It attracted considerable attention for its high luminosity, reaching M V ≈ −19.9 mag and M I ≈ −20.5 mag at peak (Stoll et al. 2011), which is even close to the superluminous SNe (M R ≤ −21 mag, Gal-Yam 2012).Smith et al. (2011) identified a luminous blue object at the site of SN 2010jl from pre-explosion images taken by the Hubble Space Telescope (HST) Wide Field Planetary Camera 2 (WFPC2).They proposed several possible scenarios for the object, including a host stellar cluster, an extremely luminous LBV-like star in quiescence, an LBV-like star in eruption, or a combination of the above.They claimed that any of these scenarios would indicate a very massive progenitor with M ini > 30 M ⊙ for SN 2010jl.Later, Fox et al. (2017) re-examined the accurate position of SN 2010jl using post-explosion images taken in 2014, 2015 and 2016 by HST Wide Field Camera 3 (WFC3), which have a much higher spatial resolution than WFPC2.They found that SN 2010jl actually has an offset from the center of the previously identified object on the pre-explosion WFPC2 images.They assumed the SN progenitor was below the detection limit of the pre-explosion WFPC2 images, which indicates a faint star below 30 M ⊙ for the progenitor.
After Fox et al. (2017), SN 2010jl was again observed by HST WFC3 in 2018 and 2020, almost 8 and 10 years after its explosion, respectively.This allows us to further check whether it has faded sufficiently and investigate the nature of the luminous blue object near the SN position on the pre-explosion WFPC2 images.As we shall show in this letter, the pre-explosion object is a blend of the SN progenitor and a close source in its vicinity, and the SN brightness is now fainter than pre-explosion level, confirming the candidate is the genuine progenitor that has now disappeared.We further carry out a detailed analysis of the progenitor and its environment in order to understand its properties.
Throughout this paper, we adopt the same distance modulus (m − M = 33.45mag) and Milky Way extinction [E(B − V ) = 0.027 mag] as in Smith et al. (2011), as well as a standard extinction law with R V = 3.1 (Cardelli et al. 1989), as suggested by Li et al. (2022).All magnitudes are reported in the Vega system.

DATA
The site of SN 2010jl was observed before explosion by HST WFPC2 in F300W/F814W in 2001, and imaged at late time by HST WFC3 in F275W/F336W in 2014and in F336W/F814W in 2015, 2016, 2018and 2020.This work utilizes data from these observations except the 2014 ones due to lack of the F814W band (Table 1).Images calibrated by the standard pipeline were retrieved from the Mikulski Archive for Space Telescopes (MAST1 ), and we re-drizzled the 2016 images with the astrodrizzle2 package with driz_cr_grow = 3 because they still suffered from high cosmic-ray contamination after the standard calibration.We performed point-spread-function (PSF) photometry using the dolphot package (Dolphin 2000) with the drizzled images ( * _drc.fits) as reference and photometric parameters recommended in the user's guide.In particular, we adopted img_RAper = 3, FitSky = 2, and Force1 = 1 for better a modelling of the crowded SN field.

SN PROGENITOR
Near the SN position (Fig. 1), we detect a single point source on the lower-resolution pre-explosion WFPC2 images (hereafter s0 ) and two spatially resolved sources on the higher-resolution late-time WFC3 images (hereafter s1 and s2 ).The measured magnitudes of s0, s1 and s2 are listed in Table 1 and shown in Fig. 2. While s1 continues to fade on the late-time images, s2 has almost constant magnitudes in both bands.Therefore, we identify s1 as SN 2010jl and s2 as a nearby field object.We note that we didn't run dolphot with the mode of warmstart, namely, all reported sources were detected and selected on its own.The source identification is consistent with Fox et al. (2017).
Image alignment was performed between the preexplosion WFPC2 images and the late-time WFC3 images with 9 common objects, reaching a differential astrometric uncertainty of 0.028".The position of s0 agrees with that of s2 within the astrometric uncertainty (also found by Fox et al. 2017).The SN position (i.e.s1 ) offsets from s0 by 0.08", which is larger than the 1σ astrometric uncertainty but still smaller than the 3σ astrometric uncertainty (0.084") and the pixel size of the WFPC2 images (0.1").Moreover, the derived source positions have their own uncertainties.Therefore, we suggest that the positions of s0 and s1 are not inconsistent with each other and that s0 is a blend of s2 and the SN progenitor at the position of s1.
Figure 2 shows the late-time brightness evolution of SN 2010jl and s2 in comparison with the pre-explosion brightness of s0.At these late-time epochs, SN 2010jl has faded significantly so that s2 can be clearly resolved with its magnitudes accurately measured.We find that the light of s2 is not sufficient to account for the brightness of s0, suggesting non-negligible light contribution from the blended SN progenitor.In particular, the composite brightness of s2 and SN 2010jl in 2018 and 2020 is significantly fainter than that of s0 (although we note the filter difference between F300W and F336W).It is possible that the SN progenitor has disappeared after explosion (see Section 6).
Next we try to decompose the spectral energy distribution (SED) of s0 to isolate the SN progenitor and derive its stellar parameters.We use a blackbody to simulate the SED of the SN progenitor; as we shall show in the next section, s2 is likely a star cluster so we adopt the BPASS v2.2 population models including interacting binaries (Eldridge et al. 2017) to model its SED.We assume that the SED of s0 is the sum of both components; we argue that this is a reasonable assumption since s2 and the SN progenitor are so close to each other while the other sources in this field are too far away to contaminate s0 significantly.We only consider Galactic extinction since the host-galaxy extinction is negligible (Smith et al. 2011).We use the pysynphot package (STScI Development Team 2013) to synthesize model magnitudes and then compare them with the observed ones of s0 and s2.The Markov Chain Monte Carlo (MCMC) technique is used (Foreman-Mackey et al. 2013) to search for the best-fitting parameters.We find that the progenitor has an effective temperature of log(T eff /K)=4.26+0.11  −0.09 and a bolometric luminosity log(L/L ⊙ )=6.52 +0.20  −0.16 ; s2 has an initial mass of 10 4 M ⊙ and a very young age of only 5 Myr.The best-fitting SEDs are displayed in Fig 3 .Figure 3 also shows the position of SN2010jl's progenitor on the Hertzsprung-Russell (HR) diagram.This analysis shows that the progenitor of SN 2010jl is a very blue and extremely luminous star that is even above the Humphreys-Davidson (HD) limit (Humphreys & Davidson 1994).

SN ENVIRONMENT
In this section we investigate the surrounding stellar populations in the environment of SN 2010jl.An environmental analysis can provide useful insights into the SN progenitor, since most massive stars form in groups and stars within each group have very similar ages and chemical abundances.We performed photometry with dolphot on all WFC3 images together, so that the derived photometry in each band is the inverse varianceweighted average of those on the individual frames, and then we selected point sources from the combined photometry catalog with the following criteria: (1) object type = 1; (2) signal-to-noise ratio ≥ 3; (3) −0.5 ≤ sharpness ≤ 0.5; (4) crowding ≤ 2. 42 sources are finally selected; for those detected in only one filter, we use artificial star tests at their positions to estimate the detection limits in the other filter.We apply a hierarchical Bayesian approach to fit model stellar populations to the data.Detailed description about this method can be found in Maund & Ramirez-Ruiz (2016) and Sun et al. (2021) (see also Maund 2017;Sun et al. 2020Sun et al. , 2022Sun et al. , 2023a,b),b).In brief, simulations of the stellar populations are based on the parsec v1.2S stellar isochrones (Bressan et al. 2012) with solar metallicity.The detected sources are assumed to be single stars or non-interacting binaries from model stellar populations with different mean log-ages and same extinction, and stars within each model population has a Gaussian log-age distribution with a small dispersion of 0.05 dex.
We find that stars in the SN environment can be fitted with two model stellar populations with mean logages of log(t/yr) = 6.66 and 6.88 dex (i.e.4.57 and 7.58 Myr, respectively).They have an extinction of E(B − V ) = 0.03 ± 0.01 , which is consistent with that for the SN (Smith et al. 2011).parsec isochrones of the corresponding ages are overplotted on the colormagnitude diagram (Fig. 4) for comparison.This suggests that SN 2010jl is located in a very young environment with very recent star formation; therefore, its progenitor is likely to be very young and massive.
It is worth noting, however, that the SN progenitor appears much brighter than the surrounding stellar populations.It is above the isochrone of the younger population by >1 mag in the F814W band, and this difference cannot be accounted by (non-interacting) binarity.It is possible that the SN progenitor is not a quiescent star and may become over luminous in outburst before the final explosion.
Additionally, a few field objects, including s2, are also much brighter than the fitted stellar isochrones; they might be star clusters, binary merger products, multiple stellar systems, and/or LBV-like stars in outburst.Assuming s2 to be a star cluster, the derived age in Sec- tion 3 is consistent with that for the younger resolved stellar population.

DISCUSSION
Using pre-explosion WFPC2 images, Smith et al. (2011) identified the luminous blue object at the site of SN 2010jl (i.e.s0 ) as its progenitor.Later, Fox et al. (2017) found the center of s0 is more consistent with that of the nearby field cluster s2 resolved on the higherresolution WFC3 images; they assumed s0 corresponds to s2 and the SN progenitor is below the detection limit of the pre-explosion images.Now that new WFC3 images have been observed at later epochs, and we are able to confirm that s0 is actually a blend of s2 and the SN progenitor.
We find the SN progenitor is an LBV-like star by decomposing the SED of s0 (Section 3).In this analysis, we only corrected for the Galactic extinction since the spectra of SN 2010jl exhibit very weak Na i lines from dust within the host galaxy (Smith et al. 2011).We note, however, several works report possible evidence for pre-existing dust in the CSM (e.g.Andrews et al. 2011;Fransson et al. 2014).The SN progenitor would become even bluer and more luminous, if non-zero circumstellar extinction is taken into account into the SED analysis.
In 2018 and 2020, SN 2010jl has faded sufficiently and the composite brightness of s1 and s2 is below the preexplosion brightness of s0 at the near-ultraviolet wavelengths (Figs. 2 and 3).This indicates that pre-explosion source at the SN position is not a stellar cluster or a field star but rather the SN progenitor.It is likely that the SN progenitor has now disappeared.We can't rule out the possibility that SN 2010ji is not a terminal explosion and the pre-explosion counterpart is still alive; in this case, the star was in outburst with higher luminosity when observed by HST in 2001 and has now returned to a fainter quiescent state at late times.It is also possible the brightness decline is (partly) due to the newly formed dust (Smith et al. 2012;Maeda et al. 2013;Gall et al. 2014).
In addition to SN 2010jl, there are another five SNe IIn with directly detected progenitors, including SN 2005gl (Gal-Yam et al. 2007), SN 2009ip (Smith et al. 2010b), SN 2010bt (Elias-Rosa et al. 2018), SN 2015bh (Elias-Rosa et al. 2016), and SN 2016jbu (Kilpatrick et al. 2018).Similar to SN 2010jl, all these SNe have late-time brightness below the progenitors, suggesting the progenitors have likely disappeared after the SN explosion (Brennan et al. 2022a;Elias-Rosa et al. 2018;Gal-Yam et al. 2007;Jencson et al. 2022;Smith et al. 2022).Obvious from the HR diagram (Fig. 3), the detected SNe IIn progenitors exhibit significant diversity in both effective temperature and luminosity.While SN 2016jbu's progenitor has a relatively low effective temperature and luminosity and is likely to be a 18-30-M ⊙ yellow hypergiant or a 17-M ⊙ star in a binary system (Kilpatrick et al. 2018;Brennan et al. 2022a,b), the other detected progenitors (including SN 2010jl) have extremely high luminosities above the HD limit (the empirical luminosity upper limit of massive stars in hydrostatic equilibrium; Humphreys & Davidson 1994) that are difficult to explain with standard models of stars in quiescence.
It has been noticed that the blue colors and high luminosities of the detected SNe IIn progenitors resemble some of the known LBVs in the Milky Way and other Local Group galaxies (e.g.η Car; Mehner et al. 2019), but none of the detected progenitors align well with the commonly adopted S Dor-like instability strip; Smith et al. 2004).The LBV giant eruptions are also similar to the SNe IIn precursor outbursts in terms of brightness and duration (Smith et al. 2010b), which can generate the dense CSM that later interact with the SN ejecta.SN 2005gl (Gal-Yam et al. 2007), SN 2009ip (Smith et al. 2010b), SN 2010bt (Elias-Rosa et al. 2018), SN 2015bh (Elias-Rosa et al. 2016), and SN 2016jbu (Kilpatrick et al. 2018).For the last two SNe (2015bh and 2016jbu), the progenitors are observed at multiple epochs and exhibit significant variability; we show their positions at all the epochs as reported in the references.For the progenitors of the first three SNe (2005gl, 2009ip and 2010bt), the effective temperatures are not well constrained by the detection (often in a single band), so we conservatively leave effective temperature as a free parameter and re-fit the reported photometry with a blackbody SED; the results are shown by the hatched regions.For comparison, We also plot single-star isochrones of Mini=30, 50, 80 M⊙ provided by MIST (black lines; Dotter 2016), standard locations of the hot/cool phases for LBVs (grey-shaded regions; Smith et al. 2004), the HD limit (grey broken line; Humphreys & Davidson 1994) and the position of the Galactic LBV η Car (black dotted line; Groh et al. 2012;Mehner et al. 2019).
However, LBVs were believed to be a transitional phase between the MS and the Wolf-Rayet stages instead of direct SN progenitors (Heger et al. 2003), although recent works suggest that LBVs could be very heterogeneous and have complicated origins and final fates (see Smith 2017;Smith et al. 2019 and references therein).Therefore, the relation between SNe IIn progenitors and LBVs remains still uncertain.
The high luminosities of SN IIn progenitors have been interpreted as them being very massive stars of several tens of and even more than a hundred solar masses (Gal-Yam et al. 2007).It has also been suggested that the high luminosities may arise from the progenitor being in outburst rather than in quiescence.This is supported by the prominent variability of several SN IIn progenitors (Pastorello et al. 2019).In this work, we find that the progenitor of SN 2010jl is significantly brighter than expected from the environmental populations, which argues against the progenitor being a quiescent star.This is also supported by the CSM of SN 2010jl, which was estimated to have a very high mass exceeding 10 M ⊙ (Ofek et al. 2014).For the typical mass-loss rate of SN IIn progenitors (10 −4 ∼ 10 0 M ⊙ /yr, Smith et al. 2014), the CSM mass corresponds to a timescale of 10 1 ∼ 10 5 years before explosion during which the SN progenitor was highly unstable and formed the dense CSM via mass ejections.Therefore, it is very likely that the SN progenitor, when observed in 2001 (or 9 years before explosion), was out of quiescence.

SUMMARY
In this paper, we perform a detailed analysis of progenitor of the Type IIn SN 2010jl based on HST images.With newer late-time WFC3 images where the SN has faded sufficiently, we are able to confirm that the pre-explosion source near the SN position on the lower-resolution WFPC2 images is actually a blend of the SN progenitor and a field star cluster.SN 2010jl is now much fainter than its progenitor; it is likely that the progenitor has now disappeared, although we can-not rule out the possibility that the progenitor is still alive and has return to quiescence or become obscured by newly formed dust.
SED analysis shows that the SN progenitor has a very high effective temperature of log(T eff /K)=4.24+0.10 −0.08 and an extreme luminosity of log(L/L ⊙ )=6.52 +0.20  −0.16 .The SN progenitor is located in a very young environment, where active star formation occurred only a few Myrs ago.However, the SN progenitor appears much more luminous than expected from the environmental stellar populations.We suggest that the SN progenitor, observed before the explosion, is most likely to be in outburst, giving rise to its extremely high luminosity.The nature of the progenitor and its detailed evolutionary history still remains to be investigated.
Figure 4 plots the images and the color-magnitude diagram (CMD) of the SN environment.

Figure 1 .Figure 2 .
Figure 1.HST images centered at SN 2010jl.All stamps have the same dimension of 0.8" × 0.8", and are aligned with North up and East to the left.The crosses mark the positions of sources near the SN position, i.e. the pre-explosion source s0 detected on the WFPC2 images (blue), and the post-explosion sources s1 and s2 detected on the WFC3 images (red).For comparison, we have astrometrically aligned the images and overplot the positions of all three sources on each image.

Figure 3 .
Figure3.Left: Photometry decomposition of the pre-explosion source s0, which is a blend of the progenitor of SN 2010jl (i.e.s1 ; blue) and the nearby field cluster s2 (orange).The black data points are based on HST observations as derived in work, while the grey data points are based on Spizter/IRAC observations as reported inFox et al. (2017).The progenitor is consistent with a blackbody of log(T eff /K)=4.26+0.11−0.09 and log(L/L⊙)=6.52+0.20 −0.16 .Right: Position of SN 2010jl's progenitor on the HR diagram (red star with error bars).The other directly detected Type IIn SN progenitors are alsodisplayed, including  SN 2005gl (Gal-Yam et al. 2007), SN 2009ip (Smith et al. 2010b), SN 2010bt (Elias-Rosa et al. 2018), SN 2015bh (Elias-Rosa  et al. 2016), and SN 2016jbu (Kilpatrick et al. 2018).For the last two SNe (2015bh and 2016jbu), the progenitors are observed at multiple epochs and exhibit significant variability; we show their positions at all the epochs as reported in the references.For the progenitors of the first threeSNe (2005gl, 2009ip and 2010bt), the effective temperatures are not well constrained by the detection (often in a single band), so we conservatively leave effective temperature as a free parameter and re-fit the reported photometry with a blackbody SED; the results are shown by the hatched regions.For comparison, We also plot single-star isochrones of Mini=30, 50, 80 M⊙ provided by MIST (black lines; Dotter 2016), standard locations of the hot/cool phases for LBVs (grey-shaded regions;Smith et al. 2004), the HD limit (grey broken line;Humphreys & Davidson 1994) and the position of the Galactic LBV η Car (black dotted line;Groh et al. 2012;Mehner et al. 2019).

Figure 4 .
Figure 4. Left and middle: WFC3 images of the SN environment with 4 different epochs combined together.Both images have the same dimension of 1.6" × 1.6" and are centered on the SN position with North up and East to the left.The small circles correspond to the detected point sources for each band within 150 pc of the SN position, and they are color-coded in the same way as in the right-hand panel.Right: color-magnitude diagram of the resolved stellar population in the SN vicinity (circles).parsec v1.2S single-star isochrones are over-plotted, which correspond to the two best-fitting model stellar populations (see text).Sources which can be well fitted by model stellar populations are shown in orange while the over-luminous ones are displayed in blue.In addition, we also show the SN progenitor (red star; based on synthetic magnitude from our SED decomposition in Section 3), the late-time behavior of SN 2010jl (green diamonds; which becomes fainter and redder with time), and the field cluster s2 (blue star).The error bars reflect photometric uncertainties.

Table 1 .
HST data and photometry.