SN Explosion Research Articles

A targeted search for strongly lensed supernovae with the Las Cumbres Observatory

ABSTRACT Gravitationally lensed supernovae (glSNe) are of interest for time delay cosmology and SN physics. However, glSNe detections are rare, owing to the intrinsic rarity of SN explosions, the necessity of alignment with a foreground lens, and the relatively short window of detectability. We present the Las Cumbres Observatory Lensed Supernova Search, LCOLSS, a targeted survey designed for detecting glSNe in known strong lensing systems. Using cadenced $r^\prime$-band imaging, LCOLSS targeted 114 galaxy–galaxy lensing systems with high expected SN rates, based on estimated star formation rates. No plausible glSN was detected by LCOLSS during our two year observing program. We carry out an analysis here to measure a detection efficiency for these observations. We then perform Monte Carlo simulations using the predicted supernova rates to determine the expected number of glSN detections. The results of the simulation suggest an expected number of detections and 68 per cent Poisson confidence intervals, $N_{\mathrm{ SN}} = 0.20, [0,2.1]$, $N_{\mathrm{ Ia}} = 0.08, [0,2.0]$, $N_{\mathrm{ CC}} = 0.12, [0,2.0]$, for all SNe, Type Ia SNe, and core-collapse (CC) SNe, respectively. These results are broadly consistent with the absence of a detection in our survey. Analysis of the survey strategy can provide insights for future efforts to develop targeted glSN discovery programs, especially considering the large anticipated yields of upcoming surveys.

Evolution and final fate of massive post-common-envelope binaries

Mergers of neutron stars (NSs) and black holes (BHs) are nowadays observed routinely thanks to gravitational-wave (GW) astronomy. In the isolated binary-evolution channel, a common-envelope (CE) phase of a red supergiant (RSG) and a compact object is crucial to sufficiently shrink the orbit and thereby enable a merger via GW emission. Here, we use the outcomes of two three-dimensional (3D) magneto-hydrodynamic CE simulations of an initially 10.0 solar-mass RSG with a 5.0 solar-mass BH and a 1.4 solar-mass NS, respectively, to explore the further evolution and final fate of the remnant binaries (post-CE binaries). Notably, the 3D simulations reveal that the post-CE binaries are likely surrounded by circumbinary disks (CBDs), which contain substantial mass and angular momentum to influence the subsequent evolution. The binary systems in MESA modelling undergo another phase of mass transfer and we find that most donor stars do not explode in ultra-stripped supernovae (SNe), but rather in Type Ib/c SNe. Without NS kicks, the final orbital configurations of our models with the BH companion are too wide to allow for a compact object merger within a Hubble time. NS kicks are actually required to sufficiently perturb the orbit and thus facilitate a merger via GW emission. Moreover, we explore the influence of CBDs observed in 3D CE simulations on the evolution and final fate of the post-CE binaries. We find that mass accretion from the disk widens the binary orbit, while resonant interactions between the CBD and the binary can shrink the separation and increase the eccentricity of the binary depending on the disk mass and lifetime. Efficient resonant contractions may even enable a BH or NS to merge with the remnant He stars before a second SN explosion, which may be observed as gamma-ray burst-like transients, luminous fast blue optical transients, and Thorne-Żytkow objects. For the surviving post-CE binaries, the CBD-binary interactions may significantly increase the GW-induced double compact merger fraction. We conclude that accounting for CBD may be crucial to better understand observed GW mergers.

Forecast of strongly lensed supernovae rates in the China Space Station Telescope surveys

Strong gravitationally lensed supernovae (SNe) are a powerful probe for cosmology and stellar physics. The relative time delays between lensed SN images provide an independent way of measuring a fundamental cosmological parameter --- the Hubble constant ---, the value of which is currently under debate. The time delays also serve as a ``time machine'', offering a unique opportunity to capture the extremely early phase of the SN explosion, which can be used to constrain the SN progenitor and explosion mechanism. Although there are only a handful of strongly lensed SN discoveries so far, which greatly hinders scientific applications, the sample size is expected to grow substantially with next-generation surveys. In this work, we investigate the capability of detecting strongly lensed SNe with the China Space Station Telescope (CSST), a two-meter space telescope to be launched around 2026. Through Monte Carlo simulations, we predict that CSST can detect 1008.53 and 51.78 strongly lensed SNe from its Wide Field Survey (WFS, covering 17,500 deg$^2$) and Deep Field Survey (DFS, covering 400 deg$^2$) over the course of ten years. In both surveys, about 35<!PCT!> of the events involve Type Ia SNe as the background sources. Our results suggest that the WFS and DFS of CSST, although not designed or optimized for discovering transients, can still make a great contribution to the strongly lensed SNe studies.

Evidence for bipolar explosions in Type IIP supernovae

Aims. Recent observations of core-collapse supernovae (SNe) suggest aspherical explosions. Globally, aspherical structures in SN explosions are thought to encode information regarding the underlying explosion mechanism. However, the exact explosion geometries from the inner cores to the outer envelopes are poorly understood. Methods. Here, we present photometric, spectroscopic, and polarimetric observations of the Type IIP SN 2021yja and discuss its explosion geometry in comparison to those of other Type IIP SNe that show large-scale aspherical structures in their hydrogen envelopes (SNe 2012aw, 2013ej and 2017gmr). Results. During the plateau phase, SNe 2012aw and 2021yja exhibit high continuum polarization characterized by two components with perpendicular polarization angles. This behavior can be interpreted as being due to a bipolar explosion, where the SN ejecta is composed of a polar (energetic) component and an equatorial (bulk) component. In such a bipolar explosion, an aspherical axis created by the polar ejecta would dominate at early phases, while the perpendicular axis along the equatorial ejecta would emerge at late phases after the photosphere in the polar ejecta has receded. Our interpretation of the explosions in SNe 2012aw and 2021yja as bipolar is also supported by other observational properties, including the time evolution of the line velocities and the line shapes in the nebular spectra. The polarization of other Type IIP SNe that show large-scale aspherical structures in the hydrogen envelope (SNe 2013ej and 2017gmr) is also consistent with the bipolar-explosion scenario, although this is not conclusive.

Gravitational-wave Radiation from the Magnetar-driven Supernovae

Rapidly spinning magnetars are potential candidates for the energy source of supernovae (SNe) and gamma-ray bursts and the most promising sources for continuous gravitational waves (GWs) detected by ground-based GW detectors. Continuous GWs can be radiated from magnetars due to magnetic-induced deformation or fluid oscillations, compatible with magnetic dipole (MD) radiation for spin-down energy. In this paper, we investigate the diverse light curves of magnetar-driven SNe in the scenario that the spin-down is dominated by GW radiation and/or MD radiation. By simulating the light curves of SNe and employing the Markov Chain Monte Carlo method, we constrain the parameters of the magnetars and SN explosions and show that the signature of GW radiation may be indicated by the bolometric luminosity curves of SNe Ic-BL 2007ru and 2009bb. We find that the ellipticity of magnetars in the order of 10−3 can be induced by the magnetic field of ∼1016 G. If such continuous GWs associated with SNe can be detected in the future by the Advanced LIGO and Virgo detectors, this would be a smoking gun for a magnetar engine powering SNe.

A closer look at the host-galaxy environment of high-velocity Type Ia supernovae

ABSTRACT Recent studies suggested that the ejecta velocity of Type Ia supernova (SN Ia) is a promising indicator in distinguishing the progenitor systems and explosion mechanisms. By classifying the SNe Ia based on their ejecta velocities, studies found that SNe Ia with high $\mathrm{Si}\, {\rm{\small II}}\, \lambda 6355$ velocities (HV SNe Ia; $v \gtrsim 12\,000$ km s−1) tend to be physically different from their normal-velocity counterparts (NV SNe Ia). In this work, we revisit the low-z sample studied in previous work and closely look into the spatially resolved environment local to the site of SN explosion. Our results reveal a possible trend (at 2.4σ significance) that HV SNe Ia are likely associated with older stellar populations than NV SNe Ia. While the trend is inconclusive, the local host-galaxy sample studied in this work is likely skewed towards massive galaxies, limiting the parameter space that we would like to investigate from the original parent sample. Nevertheless, our results do not rule out the possibility that parameters other than the host-galaxy age (such as metallicity) could be the underlying factors driving the differences between HV and NV SNe Ia due to the limitation of our data set.

Supernovae in 2023 (review): possible breakthroughs by late observations

I present a review of how late observations of supernovae, of the nebular phase, and much later of supernova remnants (SNRs), and their analysis in 2023 made progress towards possible breakthroughs in supporting the jittering jets explosion mechanism (JJEM) for core-collapse supernovae (CCSNe) and in introducing the group of lonely white dwarf (WD) scenarios for type Ia supernovae (SNe Ia). The new analyses of CCSN remnants (CCSNRs) reveal point-symmetric morphologies in a way unnoticed before in several CCSNRs. Qualitative comparison to multipolar planetary nebulae that are shaped by jets suggests that jets exploded these CCSNe, as predicted by the JJEM, but incompatible with the prediction of the delayed neutrino explosion mechanism. The spherical morphology of the ejecta Pa 30 of the historical type Iax supernova (SN Iax) of 1181 AD, which studies in 2023 revealed, is mostly compatible with the explosion of a lonely WD. Namely, at the explosion time, there is only a WD, without any close companion, although the WD was formed via a close binary interaction, i.e., binary merger. Identifying point-symmetry in SNR G1.9+0.3, a normal SN Ia and the youngest SN in the Galaxy, suggests an SN explosion of a lonely WD inside a planetary nebula (an SNIP). The group of lonely WD scenarios includes the core degenerate scenario and the double degenerate scenario with a merger to explosion delay (MED) time. SN Ia explosions of lonely WDs are common and might actually account for most (or even all) normal SNe Ia.

Searching for precursor activity of Type IIn Supernovae

We conducted a search for luminous outbursts prior to the explosion of Type IIn supernovae (SNe IIn). We built a sample of 27 objects spectroscopically classified as SNe IIn and all located at $z<0.015$. Using deep archival SN fields images ---taken up to nearly 20 years prior to the SN explosions themselves--- from transient surveys (PTF, ZTF, DES and CHASE) and major astronomical observatories (European Southern Observatory, ESO and National Optical Astronomy Observatory, NOAO), we found at least one outburst years to months before the explosion of seven SNe IIn, with the earliest precursor being 10 years prior to the explosion of SN 2019bxq. The maximum absolute magnitudes of the outbursts range between $-$11.5 mag and $-$15 mag, and the eruptive phases last for a few weeks to a few years. The $g-r$ colour measured for three objects during their outburst is relatively red, with $g-r$ ranging between 0.5 and 1.0 mag. This is similar to the colour expected during the eruptions of luminous blue variables. We note that the light curves of SNe with pre-SN outbursts have faster decline rates than those of the SNe that do not show pre-SN outbursts. SN 2011fh is remarkable, as it is still visible 12 years after the luminous SN-like event, indicating that the progenitor possibly survived, or that the interaction is still ongoing. We detect precursor activity in 29<!PCT!> of bona fide SNe IIn in our sample. However, a quantitative assessment of the observational biases affecting the sample suggests that this fraction is an underestimation of the intrinsic precursor occurrence rate.

A Bias-corrected Luminosity Function for Red Supergiant Supernova Progenitor Stars

The apparent tension between the luminosity functions of red supergiant (RSG) stars and of RSG progenitors of Type II supernovae (SNe) is often referred to as the RSG problem and it motivated some to suggest that many RSGs end their life without an SN explosion. However, the luminosity functions of RSG SN progenitors presented so far were biased to high luminosities, because the sensitivity of the search was not considered. Here, we use limiting magnitudes to calculate a bias-corrected RSG progenitor luminosity function. We find that only (36 ± 11)% of all RSG progenitors are brighter than a bolometric magnitude of −7 mag, a significantly smaller fraction than (56 ± 5)% quoted by Davies & Beasor. The larger uncertainty is due to the relatively small progenitor sample, while uncertainties on measured quantities such as magnitudes, bolometric corrections, extinction, or SN distances, only have a minor impact, as long as they fluctuate randomly for different objects in the sample. The bias-corrected luminosity functions of RSG SN progenitors and Type M supergiants in the Large Magellanic Cloud are consistent with each other, as also found by Davies & Beasor for the uncorrected luminosity function. The RSG progenitor luminosity function, hence, does not imply the existence of failed SNe. The presented statistical method is not limited to progenitor searches, but applies to any situation in which a measurement is done for a sample of detected objects, but the probed quantity or property can only be determined for part of the sample.

Discovery of an old supernova remnant candidate through carbon monoxide line emission

Most old supernova remnants (SNRs) in the Milky Way have not yet been identified. Considering their high potential number and the sufficient momentum-energy transfer to the interstellar medium (ISM), they are a key part of our understanding of the overall role of SNRs in the ISM. Here, we report our discovery of an expanding molecular shell identified by CO line observations, namely It covers a known SNR, specifically G16.0$-$0.5, and is larger in size (i.e.,\ 0.56$^ over 0.20$^ Based on its spatial and kinematic structures, weak nonthermal radio-continuum emission, and derived physical properties, we suggest that it is an old SNR. At a systemic velocity of +41.3 the best estimated kinematic distance of is sim 3.2 kpc, implying its radius of about 15.6 pc. The age of is estimated to be greater than sim $ yr, and, in a dense molecular environment, it has formed dense and thin shell layers. The kinetic energy of the expanding molecular gas of is about $6.4 $ erg, accounting for approximately six percent of the initial SN explosion energy. Although old SNRs have essentially become cold and hard to detect, our discovery suggests that they can be found by searching for CO line emissions.

From total destruction to complete survival: dust processing at different evolutionary stages in the supernova remnant Cassiopeia A

ABSTRACT The expanding ejecta of supernova remnants (SNRs) are believed to form dust in dense clumps of gas. Before the dust can be expelled into the interstellar medium and contribute to the interstellar dust budget, it has to survive the reverse shock that is generated through the interaction of the preceding supernova blast wave with the surrounding medium. The conditions under which the reverse shock hits the clumps change with remnant age and define the dust survival rate. To study the dust destruction in the SNR Cassiopeia A, we conduct magnetohydrodynamical simulations of the evolution of a supernova blast wave and of the reverse shock. In a second step, we use these evolving conditions to model clumps that are disrupted by the reverse shock at different remnant ages. Finally, we compute the amount of dust that is destroyed by the impact of the reverse shock. We find that most of the dust in the SNR is hit by the reverse shock within the first 350 yr after the SN explosion. While the dust destruction in the first 200 yr is almost complete, we expect greater dust survival rates at later times and almost total survival for clumps that are first impacted at ages beyond 1000 yr. Integrated over the entire evolution of the SNR, the dust mass shows the lowest survival fraction (17 per cent) for the smallest grains (1 nm) and the highest survival fraction (28 per cent) for the largest grains (1000 nm).

Do Type Ia Supernovae Explode inside Planetary Nebulae?

The nature of Type Ia supernova (SN Ia) explosions remains an open issue, with several contending progenitor scenarios actively being considered. One such scenario involves an SN Ia explosion inside a planetary nebula (PN) in the aftermath of a stellar merger triggered by a common envelope (CE) episode. We examine this scenario using hydrodynamic and nonequilibrium ionization simulations of the interaction between the SN ejecta and the PN cocoon into the supernova remnant (SNR) phase, focusing on the impact of the delay between the CE episode and the SN explosion. We compare the bulk dynamics and X-ray spectra of our simulated SNRs to the observed properties of known Type Ia SNRs in the Milky Way and the Magellanic Clouds. We conclude that models where the SN explosion happens in the immediate aftermath of the CE episode (with a delay ≲1000 yr) are hard to reconcile with the observations, because the interaction with the dense PN cocoon results in ionization timescales much higher than those found in any known Type Ia SNR. Models with a longer delay between the CE episode and the SN explosion (∼10,000 yr) are closer to the observations, and may be able to explain the bulk properties of some Type Ia SNRs.

Detection of Metal Enrichment by SN 2011jm in NGC 4809

Cosmic metals are believed to originate from stellar and supernovae (SNe) nucleosynthesis, dispersed into the interstellar medium (ISM) through stellar winds and supernova explosions. In this paper, we present the clear evidence of metal enrichment by a Type Ic SN 2011jm in the galaxy NGC 4809, utilizing high spatial resolution integral field unit observations obtained from the Very Large Telescope/Multi Unit Spectroscopic Explorer. Despite SN 2011jm being surrounded by metal-deficient ISM (∼0.25 Z ⊙) at a scale about 100 pc, we clearly detect enriched oxygen abundance (∼0.35 Z ⊙) and a noteworthy nitrogen-to-oxygen ratio at the SN site. Remarkably, the metal pollution is confined to a smaller scale (≲13 pc). We posit that the enhanced ionized metal stems from stellar winds emitted by massive stars or previous SN explosions. This observation may represent the first direct detection of chemical pollution by stellar feedback in star-forming galaxies beyond the Local Volume.

A spectral data release for 104 type II supernovae from the Tsinghua Supernova group

ABSTRACT We present 206 unpublished optical spectra of 104 type II supernovae (SNe II) obtained by the Xinglong 2.16-m telescope and Lijiang 2.4-m telescope during the period from 2011 to 2018, spanning the phases from about 1 to 200 d after the SN explosion. The spectral line identifications, evolution of line velocities and pseudo-equivalent widths, as well as correlations between some important spectral parameters are presented. Our sample displays a large range in expansion velocities. For instance, the Fe ii 5169 velocities measured from spectra at t ∼ 50 d after the explosion vary from 2000 to 5500 km s−1, with an average value of 3872 ± 949 km s−1. Power-law functions can be used to fit the velocity evolution, with the power-law exponent quantifying the velocity decline rate. We found an anticorrelation existing between H β velocity at mid-plateau phase and its velocity decay exponent, SNe II with higher velocities tending to have smaller velocity decay rate. Moreover, we noticed that the velocity decay rate inferred from the Balmer lines (i.e. H α and H β) have moderate correlations with the ratio of absorption to emission for H α (a/e). In our sample, two objects show possibly flash-ionized features at early phases. Besides, we noticed that multiple high-velocity components may exist on the blue side of hydrogen lines of SN 2013ab, possibly suggesting that these features arise from complex line-forming region.

Spectropolarimetry of Type II supernovae

The explosion processes of supernovae (SNe) are imprinted in their explosion geometries. The recent discovery of several highly aspherical core-collapse SNe is significant, and studying these is regarded as being crucial in order to understand the underlying explosion mechanism. Here, we study the intrinsic polarization of 15 hydrogen-rich core-collapse SNe and explore the relation between polarization and the photometric and spectroscopic properties of these objects. Our sample shows diverse properties of the continuum polarization. Most SNe show a low degree of polarization at early phases but a sudden rise to ∼1% at certain points during the photospheric phase followed by a slow decline during the tail phase, with a constant polarization angle. The variation in the timing of peak polarization values implies diversity in the explosion geometry: some SNe have aspherical structures only in their helium cores, while in other SNe such structures reach out to a significant part of the outer hydrogen envelope with a common axis from the helium core to the hydrogen envelope. Other SNe show high polarization from early phases and a change in polarization angle around the middle of the photospheric phase. This implies that the ejecta are significantly aspherical out to the outermost layer and have multi-directional aspherical structures. Exceptionally, the Type IIL SN 2017ahn shows low polarization at both the photospheric and tail phases. Our results show that the timing of the polarization rise in Type IIP SNe is likely correlated with their brightness, velocity, and the amount of radioactive Ni produced: brighter SNe with faster ejecta velocity and a larger 56Ni mass have more extended aspherical explosion geometries. In particular, there is a clear correlation between the timing of the polarization rise and the explosion energy; that is, the explosion asphericity is proportional to the explosion energy. This implies that the development of a global aspherical structure, such as a jet, might be the key for the realisation of an energetic SN in the mechanism of SN explosions.

Spectropolarimetry of Type II supernovae

We investigate the polarization spectra of hydrogen-rich core-collapse supernovae (Type II SNe). The polarization signal from SNe contains two independent components: intrinsic SN polarization and interstellar polarization (ISP). From these components, we can study the SN explosion geometry and the dust properties in their host galaxies or in the Milky Way. In this first paper, we employ a newly improved method to investigate the properties of the ISP components of 11 well-observed Type II SNe. Our analyses revealed that 10 of these 11 SNe showed a steady ISP component with a polarization degree of ≲1.0%, while one SN was consistent with zero ISP. As for the wavelength dependence, SN 2001dh (and possibly SN 2012aw) showed a non-Milky-Way-like ISP likely originating from the interstellar dust in their respective host galaxies: their polarization maxima were located at short wavelengths (≲4000 Å). Similar results have been obtained previously for highly reddened SNe. The majority of the SNe in our sample had uncertainties in the wavelength dependence of their ISP components that were too large for further consideration. Our work demonstrates that further investigation of the ISP component of the SN polarization, by applying this method to a larger SN sample, can provide new opportunities to study interstellar dust properties in external galaxies.

Characterizing the true descendants of the first stars

ABSTRACT The metal-poor stars in the Galactic halo are thought to show the imprints of the first (Pop III) stars, and thus provide a glance at the first episodes of star formation. In this work, we aim at understanding whether all very metal-poor stars formed in environments polluted by Pop III supernovae (SNe) and at what level. With a general parametric model for early metal enrichment, we study the chemical abundances (from C to Zn) of an environment imprinted by a single Pop III SN. We investigate how these abundances depend on the initial mass and internal mixing of Pop III stars, as well as on their SN explosion energy. We then study how subsequent generations of normal (Pop II) SNe affect the Pop III chemical signatures. By comparing the observed chemical abundances with our model predictions, we show that stars with [C/Fe] > +2.5 form in environments polluted purely by low-energy Pop III SNe (ESN < 2 × 1051 erg). At lower [C/Fe], stars can be imprinted either by Pop III only, or also by normal Pop II SNe. The probability of being enriched by Pop II SNe increases as [C/Fe] decreases. When Pop II stars contribute more to the pollution, they wash out the diverse chemical peculiarities left by the different Pop III SNe, and the chemical dispersion between their descendants decreases. We conclude that C-normal stars ($\rm [C/Fe] \le +0.7$) have likely been enriched by Pop II SNe at a $\ge 50~{{\ \rm per\ cent}}$ level and we identify in the abundance scatter a key diagnostic to pinpoint the signature of Pop III SNe.

3D-Spatiotemporal forecasting the expansion of supernova shells using deep learning towards high-resolution galaxy simulations

ABSTRACT Supernova (SN) plays an important role in galaxy formation and evolution. In high-resolution galaxy simulations using massively parallel computing, short integration time-steps for SNe are serious bottlenecks. This is an urgent issue that needs to be resolved for future higher-resolution galaxy simulations. One possible solution would be to use the Hamiltonian splitting method, in which regions requiring short time-steps are integrated separately from the entire system. To apply this method to the particles affected by SNe in a smoothed particle hydrodynamics simulation, we need to detect the shape of the shell on and within which such SN-affected particles reside during the subsequent global step in advance. In this paper, we develop a deep learning model, 3D-Memory In Memory (3D-MIM), to predict a shell expansion after a SN explosion. Trained on turbulent cloud simulations with particle mass mgas = 1 M⊙, the model accurately reproduces the anisotropic shell shape, where densities decrease by over 10 per cent by the explosion. We also demonstrate that the model properly predicts the shell radius in the uniform medium beyond the training data set of inhomogeneous turbulent clouds. We conclude that our model enables the forecast of the shell and its interior where SN-affected particles will be present.

The Dusty Red Supergiant Progenitor and the Local Environment of the Type II SN 2023ixf in M101

As one of the closest supernovae (SNe) in the last decade, SN 2023ixf is an unprecedented target to investigate the progenitor star that exploded. However, there is still significant uncertainty in the reported progenitor properties. In this work, we present a detailed study of SN 2023ixf’s progenitor with two independent analyses. We first modeled its spectral energy distribution (SED) based on Hubble Space Telescope optical, Spitzer mid-infrared (IR), and ground-based near-IR data. We find that stellar pulsation and circumstellar extinction have great impacts on SED fitting, and the result suggests a relatively massive red supergiant surrounded by C-rich dust with an initial mass of 16.2–17.4 M ⊙. The corresponding rate of mass loss occurring at least 3 yr before the SN explosion is about 2 × 10−4 M ⊙ yr−1. We also derived the star formation history of the SN environment based on resolved stellar populations, and the most recent star-forming epoch corresponds to a progenitor initial mass of 17–19 M ⊙, in agreement with that from our SED fitting. Therefore, we conclude that the progenitor of SN 2023ixf is close to the high-mass end for Type II SN progenitors.

Star Formation in Self-gravitating Disks in Active Galactic Nuclei. III. Efficient Production of Iron and Infrared Spectral Energy Distributions

Strong iron lines are a common feature of the optical spectra of active galactic nuclei (AGNs) and quasars from z ∼ 6−7 to the local universe, and [Fe/Mg] ratios do not show cosmic evolution. During active episodes, accretion disks surrounding supermassive black holes (SMBHs) inevitably form stars in the self-gravitating part, and these stars accrete with high accretion rates. In this paper, we investigate the population evolution of accretion-modified stars (AMSs) to produce iron and magnesium in AGNs. The AMSs, as a new type of star, are allowed to have any metallicity but without significant loss from stellar winds, since the winds are choked by the dense medium of the disks and return to the core stars. Mass functions of the AMS population show a pile-up or cutoff pile-up shape in top-heavy or top-dominant forms if the stellar winds are strong, consistent with the narrow range of supernovae (SNe) explosions driven by the known pair-instability. This provides an efficient way to produce metals. Meanwhile, SN explosions support an inflated disk as a dusty torus. Furthermore, the evolving top-heavy initial mass functions lead to bright luminosity in infrared bands in dusty regions. This contributes a new component in infrared bands, which is independent of the emissions from the central part of accretion disks, appearing as a long-term trending of the NIR continuum compared to optical variations. Moreover, the model can be further tested through reverberation mapping of emission lines, including LIGO/LISA detections of gravitational waves and signatures from spatially resolved observations of GRAVITY+/VLTI.