Supernova Models Research Articles

SN 1054 as a Pulsar-Driven Supernova: Implications for the Crab Pulsar and Remnant Evolution

Abstract One of the most studied objects in astronomy, the Crab Nebula, is the remnant of the historical supernova SN 1054. Historical observations of the supernova imply a typical supernova luminosity, but contemporary observations of the remnant imply a low explosion energy and low ejecta kinetic energy. These observations are incompatible with a standard 56Ni-powered supernova, hinting at an an alternate power source such as circumstellar interaction or a central engine. We examine SN 1054 using a pulsar-driven supernova model, similar to those used for superluminous supernovae. The model can reproduce the luminosity and velocity of SN 1054 for an initial spin period of ∼ 14 ms and an initial dipole magnetic field of 1014 − 15 G. We discuss the implications of these results, including the evolution of the Crab pulsar, the evolution of the remnant structure, formation of filaments, and limits on freely expanding ejecta. We discuss how our model could be tested further through potential light echo photometry and spectroscopy, as well as the modern analogues of SN 1054.

Enhanced production of 60Fe in massive stars.

Massive stars are a major source of chemical elements in the cosmos, ejecting freshly produced nuclei through winds and core-collapse supernova explosions into the interstellar medium. Among the material ejected, long-lived radioisotopes, such as 60Fe (iron) and 26Al (aluminum), offer unique signs of active nucleosynthesis in our galaxy. There is a long-standing discrepancy between the observed 60Fe/26Al ratio by γ-ray telescopes and predictions from supernova models. This discrepancy has been attributed to uncertainties in the nuclear reaction networks producing 60Fe, and one reaction in particular, the neutron-capture on 59Fe. Here we present experimental results that provide a strong constraint on this reaction. We use these results to show that the production of 60Fe in massive stars is higher than previously thought, further increasing the discrepancy between observed and predicted 60Fe/26Al ratios. The persisting discrepancy can therefore not be attributed to nuclear uncertainties, and points to issues in massive-star models.

A Possible Formation Scenario of the Gaia BH1: Inner Binary Merger in Triple Systems

Based on astrometric measurements and spectral analysis from Gaia DR3, two quiescent black hole (BH) binaries, Gaia BH1 and BH2, have been identified. Their origins remain controversial, particularly for Gaia BH1. By considering a rapidly rotating (ω/ω crit = 0.8) and strongly magnetized (B 0 = 5000 G) merger product, we find that, at typical Galactic metallicity, the merger product can undergo efficient chemically homogeneous evolution. This results in the merger product having a significantly smaller radius during its evolution compared to that of a normally evolving massive star. Under the condition that the initial triple stability is satisfied, we use the Multiple Stellar Evolution code and the MESA code to identify an initial hierarchical triple that can evolve into Gaia BH1. It initially consists of three stars with masses of 9.03 M ⊙, 3.12 M ⊙, and 1 M ⊙, with inner and outer orbital periods of 2.21 days and 121.92 days, and inner and outer eccentricities of 0.41 and 0.45, respectively. This triple initially experiences triple evolution dynamics instability (TEDI) followed by Roche lobe overflow (RLOF). During RLOF, the inner orbit shrinks, and tidal effects gradually suppress the TEDI. Eventually, the inner binary undergoes a merger through contact (or collision). Finally, using models of rapidly rotating and strongly magnetic stars, along with standard core-collapse supernova (SN) or failed supernova (FSN) models, we find that a postmerger binary (PMB) consisting of an 12.11 M ⊙ merger product and a 1 M ⊙ companion star (originally an outer tertiary) can avoid RLOF. After an SN or FSN with a low ejected mass of ∼0.22 M ⊙ and a low kick velocity ( 46−33+25kms−1 or 9−8+16kms−1 ), the PMB can form Gaia BH1 in the Galactic disk.

Gravitational-wave and Gravitational-wave Memory Signatures of Core-collapse Supernovae

In this paper, we calculate the energy, signal-to-noise ratio (SNR), detection range, and angular anisotropy of the matter, matter memory, and neutrino memory gravitational-wave (GW) signatures of 21 three-dimensional initially nonrotating core-collapse supernova (CCSN) models carried to late times. We find that inferred energy, SNR, and detection range are angle-dependent quantities, and that the spread of possible energy, signal to noise, and detection ranges across all viewing angles generally increases with progenitor mass. When examining the low-frequency matter memory and neutrino memory components of the signal, we find that the neutrino memory is the most detectable component of a CCSN GW signal, and that DECIGO is best equipped to detect both matter memory and neutrino memory. Moreover, we find that the polarization angle between the h + and h × strains serves as a unique identifier of matter and neutrino memory. Finally, we develop a Galactic density- and stellar mass-weighted formalism to calculate the rate at which we can expect to detect CCSN GW signals with the Advanced Laser Interferometer Gravitational-Wave Observatory (aLIGO). When considering only the matter component of the signal, the aLIGO detection rate is around 65% of the total Galactic supernova rate, but increases to 90% when incorporating the neutrino memory component. We find that all future detectors (Einstein Telescope, Cosmic Explorer, DECIGO) will be able to detect CCSN GW signals from the entire Galaxy, and for the higher-mass progenitors even into the Local Group of galaxies.

Correction to: Carbon monoxide cooling in radiative transfer modelling of supernovae

Correction to: Carbon monoxide cooling in radiative transfer modelling of supernovae

Hydrodynamics and Nucleosynthesis of Jet-driven Supernovae. II. Comparisons with Abundances of Extremely Metal-poor Galaxies and Constraints on Supernova Progenitors

The spectra of several galaxies, including extremely metal-poor galaxies from EMPRESS, have shown that the abundances of some Si-group elements differ from “spherical” explosion models of massive stars. This leads to the speculation that these galaxies have experienced supernova explosions with high asphericity, where mixing and fallback of the inner ejecta with the outer material lead to the distinctive chemical compositions. In this paper, we consider the jet-driven supernova models by direct 2D hydrodynamics simulations using progenitors of about 20–25 M ⊙ at zero metallicity. We investigate how the abundance patterns depend on the progenitor mass, mass cut, and asphericity of the explosion. We compare the observable with available supernova and galaxy catalogs based on 56Ni, ejecta mass, and individual element ratios. The proximity of our results with the observational data signifies the importance of aspherical supernova explosions in chemical evolution of these galaxies. Our models will provide the theoretical counterpart for understanding the chemical abundances of high-z galaxies measured by the James Webb Space Telescope.

Insights into the Production of 44Ti and Nickel Isotopes in Core-collapse Supernovae

We report nucleosynthetic results for both 44Ti and nickel isotopes for 18 three-dimensional (3D) core-collapse supernova (CCSN) simulations extended to ∼20 s after bounce. We find that many of our long-term models are able to achieve 44Ti/56Ni ratios similar to that observed in Cassiopeia A, and modern supernova models can synthesize up to 2 × 10−4 M ⊙ of 44Ti. Neutrino-driven winds and the fact that there can be simultaneous accretion and explosion in 3D models of CCSNe play central roles in its production. We conclude that the 44Ti underproduction problem in previous CCSN models is no longer an issue. In addition, we discuss the production of both 57Ni and stable nickel/iron ratios and compare our results to observations of SN 1987A and the Crab.

Circumstellar Interaction in the Ultraviolet Spectra of SN 2023ixf 14–66 Days After Explosion

SN 2023ixf was discovered in M101 within a day of the explosion and rapidly classified as a Type II supernova with flash features. Here we present ultraviolet (UV) spectra obtained with the Hubble Space Telescope 14, 19, 24, and 66 days after the explosion. Interaction between the supernova ejecta and circumstellar material (CSM) is seen in the UV throughout our observations in the flux of the first three epochs and asymmetric Mg ii emission on day 66. We compare our observations to CMFGEN supernova models that include CSM interaction ( Ṁ<10−3 M ⊙ yr−1) and find that the power from CSM interaction is decreasing with time, from L sh ≈ 5 × 1042 erg s−1 to L sh ≈ 1 × 1040 erg s−1 between days 14 and 66. We examine the contribution of individual atomic species to the spectra on days 14 and 19, showing that the majority of the features are dominated by iron, nickel, magnesium, and chromium absorption in the ejecta. The UV spectral energy distribution of SN 2023ixf sits between that of supernovae, which show no definitive signs of CSM interaction, and those with persistent signatures assuming the same progenitor radius and metallicity. Finally, we show that the evolution and asymmetric shape of the Mg ii λ λ 2796, 2802 emission are not unique to SN 2023ixf. These observations add to the early measurements of dense, confined CSM interaction, tracing the mass-loss history of SN 2023ixf to ∼33 yr prior to the explosion and the density profile to a radius of ∼5.7 × 1015 cm. They show the relatively short evolution from a quiescent red supergiant wind to high mass loss.

Viewing Explosion Models of Type Ia Supernovae through Insights from Terrestrial Cellular Detonation.

The cellular structure is considered to be a key as a criterion in initiation, propagation, and quenching of terrestrial detonation. While a few studies on type Ia supernovae, which are known to involve detonation, have addressed the importance of the cellular structure, further detailed treatment will benefit enhanced understanding of the explosion outcomes. In the present study, we bridge this gap in the astrophysics and engineering fields, focusing on the detonation in a helium-rich white dwarf envelope as the triggering process for the so-called double-detonation model. The cellular structures are quantified via high-resolution two-dimensional simulations. We demonstrate that widely accepted terrestrial-experimental criteria for quenching and initiation of detonation can indeed explain the results of previous hydrodynamic simulations very well. The present study highlights the potential of continuing to apply the insight from terrestrial detonation experiments to astrophysical problems, specifically the long unresolved problem of the explosion mechanism of type Ia supernovae.

Neutrino quantum kinetics in a core-collapse supernova

Our understanding of neutrino flavor conversion in the supernova core is still preliminary, despite its likely relevance to the neutrino-driven supernova mechanism. We present multi-angle and multi-energy numerical simulations of neutrino quantum kinetics within a spherically symmetric shell in the proximity of the region of neutrino decoupling. We rely on inputs from a one-dimensional core-collapse supernova model with a mass of 18.6 M ⊙ and find that, at early post-bounce times (t pb ≲ 0.5 s), no crossing is present in the angular distribution of the electron neutrino lepton number and flavor conversion is triggered by slow collective instabilities. Angular crossings appear for t pb ≳ 0.5 s and fast flavor conversion leads to flavor equipartition, with the spectral energy distribution of ν e (ν̅e) and ν x (ν̅ x ) becoming comparable. Notably, flavor equipartition is not a generic outcome of fast flavor conversion, rather it is a consequence of the relatively similar properties of neutrinos of different flavors characterizing the late accretion phase. Artificially tweaking the collision term to introduce an electron lepton number angular crossing for t pb ≲ 0.5 s, we observe that flavor equipartition is not achieved. While our findings are restricted to a specific supernova model, and they only take into account the feedback of the neutrino background on the flavor conversion, they suggest a rich phenomenology in the supernova core as a function of the post-bounce time which needs to be further explored to assess its impact on the explosion mechanism.

Muon-induced collisional flavor instability in core-collapse supernova

Neutrinos are known to undergo flavor conversion among their three flavors. In the theoretical modeling of core-collapse supernova, there has been a great deal of attention to recent discoveries of a new type of neutrino flavor conversions, namely collisional flavor instability (CFI), in which the instability is induced by the flavor-dependent decoherence due to the disparity of neutrino-matter interactions among flavors. In this paper, we study how the appearance of on-shell muons and associated neutrino-matter interactions can impact CFIs based on linear stability analysis of flavor conversions. Some striking results emerge from the present study. First, we analytically show that breaking beta and pair equilibrium is a necessary condition to trigger CFIs. This also indicates that CFIs with on-shell muons could appear in eτ and μτ neutrino mixing sectors in very high-density region (≳1013 g/cm3), exhibiting a possibility of large impacts of CFIs on core-collapse supernova. Second, resonancelike CFIs, having a much higher growth rate than normal CFIs, can be triggered by muons. The resonance point of CFIs is different between eτ and μτ sectors; the former (latter) occurs at μe(μ)=μn−μp, where μi denotes the chemical potential of i constitute (n and p represent neutrons and protons, respectively). Our result suggests that the nonlinear evolution of CFI with on-shell muons would induce flavor conversions with the complex interplay among all three different neutrino-mixing sectors. Published by the American Physical Society 2024

Chemo-dynamical Evolution of Simulated Satellites for a Milky Way–like Galaxy

The chemical abundances of Milky Way’s (MW's) satellites reflect their star formation histories (SFHs), yet, due to the difficulty of determining the ages of old stars, the SFHs of most satellites are poorly measured. Ongoing and upcoming surveys will obtain around 10 times more medium-resolution spectra for stars in satellites than are currently available. To correctly extract SFHs from large samples of chemical abundances, the relationship between chemical abundances and SFHs needs to be clarified. Here, we perform a high-resolution cosmological zoom-in simulation of a MW-like galaxy with detailed models of star formation, supernova (SN) feedback, and metal diffusion. We quantify SFHs, metallicity distribution functions, and the α-element (Mg, Ca, and Si) abundances in satellites of the host galaxy. We find that star formation in most simulated satellites is quenched before infalling to their host. Star formation episodes in simulated satellites are separated by a few hundred Myr owing to SN feedback; each star formation event produces groups of stars with similar [α/Fe] and [Fe/H]. We then perform a mock observation of the upcoming Subaru Prime Focus Spectrograph (PFS) observations. We find that Subaru PFS will be able to detect distinct groups of stars in [α/Fe] versus [Fe/H] space, produced by episodic star formation. This result means that episodic SFHs can be estimated from the chemical abundances of ≳1000 stars determined with medium-resolution spectroscopy.

Performance of SK-Gd’s Upgraded Real-time Supernova Monitoring System

Among multimessenger observations of the next Galactic core-collapse supernova, Super-Kamiokande (SK) plays a critical role in detecting the emitted supernova neutrinos, determining the direction to the supernova (SN), and notifying the astronomical community of these observations in advance of the optical signal. In 2022, SK has increased the gadolinium dissolved in its water target (SK-Gd) and has achieved a Gd concentration of 0.033%, resulting in enhanced neutron detection capability, which in turn enables more accurate determination of the supernova direction. Accordingly, SK-Gd’s real-time supernova monitoring system has been upgraded. SK_SN Notice, a warning system that works together with this monitoring system, was released on 2021 December 13, and is available through GCN Notices. When the monitoring system detects an SN-like burst of events, SK_SN Notice will automatically distribute an alarm with the reconstructed direction to the supernova candidate within a few minutes. In this paper, we present a systematic study of SK-Gd’s response to a simulated Galactic SN. Assuming a supernova situated at 10 kpc, neutrino fluxes from six supernova models are used to characterize SK-Gd’s pointing accuracy using the same tools as the online monitoring system. The pointing accuracy is found to vary from 3° to 7° depending on the models. However, if the supernova is closer than 10 kpc, SK_SN Notice can issue an alarm with three-degree accuracy, which will benefit follow-up observations by optical telescopes with large fields of view.

Primary and secondary source of energy in the superluminous supernova 2018ibb

We examined a possible pair-instability origin for the superluminous supernova 2018ibb. As the base model, we used a non-rotating stellar model with an initial mass of 250 at about 1/15 solar metallicity. We considered three versions of the model as input for radiative transfer simulations done with the stella and artis codes: with 25 of 56Ni, with 34 of 56Ni, and a chemically mixed case with 34 of 56Ni. We present light curves and spectra in comparison to the observed data of SN\,2018ibb, and conclude that the pair-instability supernova model with 34 of 56Ni explains the broadband light curves reasonably well between $-100$ and 250 days around the peak. Our synthetic spectra have many similarities with the observed spectra. The luminosity excess in the light curves and the blue-flux excess in the spectra can be explained by an additional energy source, which may be an interaction of the supernova ejecta with circumstellar matter. We discuss possible mechanisms that could have caused the circumstellar matter to be ejected in the decades before the pair-instability explosion.

Isotope studies of presolar silicon carbide grains from supernovae: new constraints for hydrogen-ingestion supernova models

ABSTRACT We report isotope data for C, N, Al, Si, and S of 33 presolar SiC and Si3N4 grains (0.3–1.6 $\mu$m) of Type X, C, D, and N from the Murchison CM2 meteorite of likely core-collapse supernova (CCSN) origin which we discuss together with data of six SiC X grains from an earlier study. The isotope data are discussed in the context of hydrogen ingestion supernova (SN) models. We have modified previously used ad-hoc mixing schemes in that we considered (i) heterogeneous H ingestion into the He shell of the pre-SN star, (ii) a variable C-N fractionation for the condensation of SiC grains in the SN ejecta, and (iii) smaller mass units for better fine-tuning. With our modified ad-hoc mixing approach over small scales (0.2–0.4 M⊙), with major contributions from the O-rich O/nova zone, we find remarkably good fits (within a few per cent) for 12C/13C, 26Al/27Al, and 29Si/28Si ratios. The 14N/15N ratio of SiC grains can be well matched if variable C-N fractionation is considered. However, the Si3N4 isotope data point to overproduction of 15N in hydrogen ingestion CCSN models and lower C-N fractionation during SiC condensation than applied here. Our ad-hoc mixing approach based on current CCSN models suggests that the O-rich O/nova zone, which uniquely combines explosive H- and He-burning signatures, is favourable for SiC and Si3N4 formation. The effective range of C/O abundance variations in the He shell triggered by H ingestion events in the massive star progenitor is currently not well constrained and needs further investigation.

Carbon monoxide cooling in radiative transfer modelling of supernovae

ABSTRACT Carbon monoxide (CO) emission has been observed in a number of core-collapse supernovae (SNe) and is known to be an important coolant at late times. We have implemented a chemical reaction network in the radiative-transfer code cmfgen to investigate the formation of CO and its impact on SN ejecta. We calculate two 1D SN models with and without CO: a BSG explosion model at one nebular epoch and a full time-sequence (50–300 d) for a red supergiant explosion. In both models, CO forms at nebular times in the dense, inner regions at velocities &amp;lt;2000 km s$^{-1}$ where line emission from CO can dominate the cooling and reduce the local temperature by as much as a factor of 2, weakening emission lines and causing the optical light curve to fade faster. That energy is instead emitted in CO bands, primarily the fundamental band at $\sim 4.5 \mathrm{\mu m}$, which accounts for up to 20 per cent of the total luminosity at late times. However, the non-monotonic nature of the CO cooling function can cause numerical difficulties and introduce multiple temperature solutions. This issue is compounded by the sensitivity of the CO abundance to a few reaction rates, many of which have large uncertainties or disparate values across literature sources. Our results also suggest that, in many SNe, CO-level populations are far from their local thermodynamic equilibrium values. Unfortunately, accurate collisional data, necessary to compute non-local thermodynamic equilibrium populations, are limited to a few transitions.

Numerical relativity simulations of black hole and relativistic jet formation

ABSTRACT We investigate impacts of stellar rotation and magnetic fields on black hole (BH) formation and its subsequent explosive activities, by conducting axisymmetric radiation-magnetohydrodynamics simulations of gravitational collapse of a 70 $\mathrm{M}_\odot$ star with two-moment multi energy neutrino transport in full general relativity for the first time. Due to its dense stellar structure, all models cannot avoid the eventual BH formation even though a strongly magnetized model experiences the so-called magnetorotational explosion prior to the BH formation. One intriguing phenomenon observed in the strongly magnetized model is the formation of a relativistic jet in the post-BH formation. The relativistic jet is the outcome of a combination of strong magnetic fields and low-density materials above the BH. The jet further enhances the explosion energy beyond $\sim 10^{52}$ erg, which is well exceeding the gravitational overburden ahead of the shock. Our self-consistent supernova models demonstrate that rotating magnetized massive stars at the high-mass end of supernova progenitors could be a potential candidate of hypernova and long gamma-ray burst progenitors.

Peculiarities of the chemical enrichment of metal-poor stars in the Milky Way Galaxy

The oldest stars in the Milky Way are metal-poor with Fe/H $&lt;$ -- 1.0, displaying peculiar elemental abundances compared to solar values. The relative variations in the chemical compositions among stars is also increasing with decreasing stellar metallicity, allowing for the pure signature of unique nucleosynthesis processes to be revealed. The study of the $r$-process is, for instance, one of the main goals of stellar archaeology and metal-poor stars exhibit an unexpected complexity in the stellar production of the $r$-process elements in the early Galaxy. In this work, we report the atmospheric parameters, main dynamic properties, and the abundances of four metal-poor stars: HE 1523-–0901, HD 6268, HD 121135, and HD 195636 (--1.5 $&gt;$ Fe/H The abundances were derived from spectra obtained with the HRS echelle spectrograph at the Southern African Large Telescope, using both local and non-local thermodynamic equilibrium (LTE and NLTE) approaches, with the average error between 0.10 and 0.20 dex. Based on their kinematical properties, we show that HE 1523--0901 and HD 195636 are halo stars with typical high velocities. In particular, HD 121135 displays a peculiar kinematical behaviour, making it unclear whether it is a halo or an accreted star. Furthermore, HD 6268 is possibly a rare prototype of very metal-poor thick disk stars. The abundances derived for our stars are compared with theoretical stellar models and with other stars with similar metallicity values from the literature. HD 121135 is Al-poor and Sc-poor, compared to stars observed in the same metallicity range (--1.62 $&gt;$ Fe/H $&gt;$--1.12). The most metal-poor stars in our sample, HE 1523 -- 0901, HD 6268, and HD 195636, exhibit anomalies that are better explained by supernova models from fast-rotating stellar progenitors for elements up to the Fe group. Compared to other stars in the same metallicity range, their common biggest anomaly is represented by the low Sc abundances. If we consider the elements beyond Zn, HE 1523--0901 can be classified as an r-II star, HD 6268 as an r-I candidate, and HD 195636 and HD 121135 exhibiting a borderline $r$-process enrichment between limited-r and r-I star. Significant relative differences are observed between the $r$-process signatures in these stars.

SN 2020zbf: A fast-rising hydrogen-poor superluminous supernova with strong carbon lines

SN 2020zbf is a hydrogen-poor superluminous supernova (SLSN) at z = 0.1947 that shows conspicuous C II features at early times, in contrast to the majority of H-poor SLSNe. Its peak magnitude is Mg = −21.2 mag and its rise time (≲26.4 days from first light) places SN 2020zbf among the fastest rising type I SLSNe. We used spectra taken from ultraviolet (UV) to near-infrared wavelengths to identify spectral features. We paid particular attention to the C II lines as they present distinctive characteristics when compared to other events. We also analyzed UV and optical photometric data and modeled the light curves considering three different powering mechanisms: radioactive decay of 56Ni, magnetar spin-down, and circumstellar medium (CSM) interaction. The spectra of SN 2020zbf match the model spectra of a C-rich low-mass magnetar-powered supernova model well. This is consistent with our light curve modeling, which supports a magnetar-powered event with an ejecta mass Mej = 1.5 M⊙. However, we cannot discard the CSM-interaction model as it may also reproduce the observed features. The interaction with H-poor, carbon-oxygen CSM near peak light could explain the presence of C II emission lines. A short plateau in the light curve around 35–45 days after peak, in combination with the presence of an emission line at 6580 Å, can also be interpreted as being due to a late interaction with an extended H-rich CSM. Both the magnetar and CSM-interaction models of SN 2020zbf indicate that the progenitor mass at the time of explosion is between 2 and 5 M⊙. Modeling the spectral energy distribution of the host galaxy reveals a host mass of 108.7 M⊙, a star formation rate of 0.24−0.12+0.41 M⊙ yr−1, and a metallicity of ∼0.4 Z⊙.

General-relativistic Radiation Transport Scheme in Gmunu. II. Implementation of Novel Microphysical Library for Neutrino Radiation—Weakhub

We introduce Weakhub, a novel neutrino microphysics library that provides opacities and kernels beyond conventional interactions used in the literature. This library includes neutrino–matter, neutrino–neutrino interactions and plasma process, along with corresponding weak and strong corrections. A full kinematics approach is adopted for the calculations of β-processes, incorporating various weak corrections and medium modifications due to the nuclear equation of state. Calculations of plasma processes, electron neutrino–antineutrino annihilation, and nuclear de-excitation are also included. We also present the detailed derivations of weak interactions and the coupling to the two-moment based general-relativistic multigroup radiation transport in the general-relativistic multigrid numerical (Gmunu) code. We compare the neutrino opacity spectra for all interactions and estimate their contributions at hydrodynamical points in core-collapse supernovae and binary neutron star (BNS) postmerger remnants, and predict the effects of improved opacities in comparison to conventional ones for a BNS postmerger at a specific hydrodynamical point. We test the implementation of the conventional set of interactions by comparing it to an open-source neutrino library NuLib in a core-collapse supernova simulation. We demonstrate good agreement with discrepancies of less than ∼10% in luminosity for all neutrino species, while also highlighting the reasons contributing to the differences. To compare the advanced interactions to the conventional set in core-collapse supernova modeling, we perform simulations to analyze their impacts on neutrino signatures, hydrodynamical behaviors, and shock dynamics, showing significant deviations.