Core-collapse Supernova Models Research Articles

Nucleosynthesis in the Innermost Ejecta of Magnetorotational Supernova Explosions in Three Dimensions

Core-collapse supernova (CCSN) explosions powered by rotation and magnetic fields present an interesting astrophysical site for nucleosynthesis that potentially contributes to the production of r-process elements. Here we present yields of the innermost ejecta in 3D magnetorotational CCSN models simulated using the CoCoNuT-FMT code. Strong magnetic fields tap the rotational energy of the proto−neutron star and lead to earlier and more energetic (∼3 × 1051 erg) explosions than typical neutrino-driven CCSNe. Compared to a reference nonmagnetic model, the ejecta in the magnetorotational models have much more neutron-rich components with Y e down to ∼0.25. Our post-processing calculations with the reaction network SkyNet show significant production of weak r-process elements up to mass number ∼130. We find negligible differences in the synthesis of heavy elements between two magnetorotational models with different initial field strengths of 1010 and 1012 G, in accord with their similar explosion dynamics. The magnetorotational models produce about ∼0.19 and 0.14 M ☉ of radioactive 56Ni, on the low end of inferred hypernova nickel masses. The yields are publicly available at Zenodo (doi: 10.5281/zenodo.10578981) for comparison with stellar abundance patterns, inclusion in modeling galactic chemical evolution, and comparison with other yield calculations. Our results add to the yet-restricted corpus of nucleosynthesis yields from 3D magnetorotational supernova simulations and will help quantify yield uncertainties.

Nucleosynthesis in the Innermost Ejecta of Magnetorotational Supernova Explosions in Three Dimensions

Core-collapse supernova (CCSN) explosions powered by rotation and magnetic fields present an interesting astrophysical site for nucleosynthesis that potentially contributes to the production of r-process elements. Here we present yields of the innermost ejecta in 3D magnetorotational CCSN models simulated using the CoCoNuT-FMT code. Strong magnetic fields tap the rotational energy of the proto−neutron star and lead to earlier and more energetic (∼3 × 1051 erg) explosions than typical neutrino-driven CCSNe. Compared to a reference nonmagnetic model, the ejecta in the magnetorotational models have much more neutron-rich components with Y e down to ∼0.25. Our post-processing calculations with the reaction network SkyNet show significant production of weak r-process elements up to mass number ∼130. We find negligible differences in the synthesis of heavy elements between two magnetorotational models with different initial field strengths of 1010 and 1012 G, in accord with their similar explosion dynamics. The magnetorotational models produce about ∼0.19 and 0.14 M ☉ of radioactive 56Ni, on the low end of inferred hypernova nickel masses. The yields are publicly available at Zenodo (doi: 10.5281/zenodo.10578981) for comparison with stellar abundance patterns, inclusion in modeling galactic chemical evolution, and comparison with other yield calculations. Our results add to the yet-restricted corpus of nucleosynthesis yields from 3D magnetorotational supernova simulations and will help quantify yield uncertainties.

Typing supernova remnant G352.7−0.1 using XMM–Newton X-ray observations

ABSTRACT G352.7−0.1 is a mixed-morphology (MM) supernova remnant (SNR) with multiple radio arcs and has a disputed supernova origin. We conducted a spatially resolved spectroscopic study of the remnant with XMM–Newton X-ray data to investigate its explosion mechanism and explain its morphology. The global X-ray spectra of the SNR can be adequately reproduced using a metal-rich thermal plasma model with a temperature of ∼2 keV and ionization time-scale of ∼3 × 1010 cm−3 s. Through a comparison with various supernova nucleosynthesis models, we found that observed metal properties from Mg to Fe can be better described using core-collapse supernova models, while thermonuclear models fail to explain the observed high Mg/Si ratio. The best-fit supernova model suggests a ∼13 M⊙ progenitor star, consistent with previous estimates using the wind bubble size. We also discussed the possible mechanisms that may lead to SNR G352.7−0.1 being an MMSNR. By dividing the SNR into several regions, we found that the temperature and abundance do not significantly vary with regions, except for a decreased temperature and abundance in a region interacting with molecular clouds. The brightest X-ray emission of the SNR spatially matches with the inner radio structure, suggesting that the centrally filled X-ray morphology results from a projection effect.

On the core-collapse supernova explanation for LAMOST J1010 + 2358

ABSTRACT Low-metallicity very massive stars with an initial mass of ∼140–$260\, \mathrm{M}_\odot$ are expected to end their lives as pair-instability supernovae (PISNe). The abundance pattern resulting from a PISN differs drastically from regular core-collapse supernova (CCSN) models and is expected to be seen in very metal-poor (VMP) stars of [Fe/H] ≲ −2. Despite the routine discovery of many VMP stars, the unique abundance pattern expected from PISNe has not been unambiguously detected. The recently discovered VMP star LAMOST J1010 + 2358, however, shows a peculiar abundance pattern that is remarkably well fit by a PISN, indicating the potential first discovery of a bonafide star born from gas polluted by a PISN. In this paper, we study the detailed nucleosynthesis in a large set of models of CCSN of Pop III and Pop II star of metallicity [Fe/H] = −3 with masses ranging from 12 to $30\, \mathrm{M}_\odot$. We find that the observed abundance pattern in LAMOST J1010 + 2358 can be fit at least equally well by CCSN models of ∼12–$14\, \mathrm{M}_\odot$ that undergo negligible fallback following the explosion. The best-fitting CCSN models provide a fit that is even marginally better than the best-fitting PISN model. We conclude the measured abundance pattern in LAMOST J1010 + 2358 could have originated from a CCSN and therefore cannot be unambiguously identified with a PISN given the set of elements measured in it to date. We identify key elements that need to be measured in future detections in stars like LAMOST J1010 + 2358 that can differentiate between CCSN and PISN origin.

Chemical enrichment of ICM within the Centaurus cluster – I. Radial profiles

ABSTRACT We examine deep XMM–Newton European Photon Imaging Camera pn observations of the Centaurus cluster to study the hot intracluster medium (ICM) and radial metal distributions within such an environment. We found that the best-fitting spectral model corresponds to a lognormal temperature distribution, with discontinuities around ∼10, ∼50, and ∼100 kpc, also observed in the abundance distributions. We measured the radial profiles of O, Si, S, Ar, Ca, and Fe. These profiles reveal prominent negative gradients for distances <90 kpc, which then transition to flatter profiles. We modelled X/Fe ratio profiles with a linear combination of core-collapse supernova and Type Ia supernova (SNIa) models. The best-fitting model suggests a uniform SNIa percentage contribution to the total cluster enrichment, thus supporting an early enrichment of the ICM, with most of the metals present being produced before clustering.

The γ-process nucleosynthesis in core-collapse supernovae

Context. The γ-process nucleosynthesis in core-collapse supernovae is generally accepted as a feasible process for the synthesis of neutron-deficient isotopes beyond iron. However, crucial discrepancies between theory and observations still exist: the average yields of γ-process nucleosynthesis from massive stars are still insufficient to reproduce the solar distribution in galactic chemical evolution calculations, and the yields of the Mo and Ru isotopes are a factor of ten lower than the yields of the other γ-process nuclei. Aims. We investigate the γ-process in five sets of core-collapse supernova models published in the literature with initial masses of 15, 20, and 25 M⊙ at solar metallicity. Methods. We compared the γ-process overproduction factors from the different models. To highlight the possible effect of nuclear physics input, we also considered 23 ratios of two isotopes close to each other in mass relative to their solar values. Further, we investigated the contribution of C–O shell mergers in the supernova progenitors as an additional site of the γ-process. Results. Our analysis shows that a large scatter among the different models exists for both the γ-process integrated yields and the isotopic ratios. We find only ten ratios that agree with their solar values, all the others differ by at least a factor of three from the solar values in all the considered sets of models. The γ-process within C–O shell mergers mostly influences the isotopic ratios that involve intermediate and heavy proton-rich isotopes with A > 100. Conclusions. We conclude that there are large discrepancies both among the different data sets and between the model predictions and the solar abundance distribution. More calculations are needed; particularly updating the nuclear network, because the majority of the models considered in this work do not use the latest reaction rates for the γ-process nucleosynthesis. Moreover, the role of C–O shell mergers requires further investigation.

Simple method for determining asymptotic states of fast neutrino-flavor conversion

Neutrino-neutrino forward scatterings potentially induce collective neutrino oscillation in dense neutrino gases in astrophysical sites such as core-collapse supernovae (CCSN) and binary neutron star mergers (BNSM). In this paper, we present a detailed study of fast neutrino-flavor conversion (FFC), paying special attention to asymptotic states, by linear stability analysis and local simulations with a periodic boundary condition. We find that asymptotic states can be characterized by two key properties of FFC: (1) the conservation of lepton number for each flavor of neutrinos and (2) the disappearance of ELN(electron neutrino-lepton number)-XLN(heavy-leptonic one) angular crossings in the spatial- or time-averaged distributions. The system which initially has the positive (negative) ELN-XLN density reaches a flavor equipartition in the negative (positive) ELN-XLN angular directions, and the other part compensates it to preserve the conservation laws. These properties of FFCs offer an approximate scheme determining the survival probability of neutrinos in asymptotic states without solving quantum kinetic equations. We also demonstrate that the total amount of flavor conversion can vary with species-dependent neutrino distributions for identical ELN-XLN ones. Our results suggest that even shallow or narrow ELN angular crossings have the ability to drive large flavor conversion, exhibiting the need for including the effects of FFCs in the modeling of CCSN and BNSM.

Applications of machine learning to detecting fast neutrino flavor instabilities in core-collapse supernova and neutron star merger models

Neutrinos propagating in a dense neutrino gas, such as those expected in core-collapse supernovae (CCSNe) and neutron star mergers (NSMs), can experience fast flavor conversions on relatively short scales. This can happen if the neutrino electron lepton number ($\nu$ELN) angular distribution crosses zero in a certain direction. Despite this, most of the state-of-the-art CCSN and NSM simulations do not provide such detailed angular information and instead, supply only a few moments of the neutrino angular distributions. In this study we employ, for the \emph{first} time, a machine learning (ML) approach to this problem and show that it can be extremely successful in detecting $\nu$ELN crossings on the basis of its zeroth and first moments. We observe that an accuracy of $\sim95\%$ can be achieved by the ML algorithms, which almost corresponds to the Bayes error rate of our problem. Considering its remarkable efficiency and agility, the ML approach provides one with an unprecedented opportunity to evaluate the occurrence of FFCs in CCSN and NSM simulations \emph{on the fly}. We also provide our ML methodologies on GitHub.

Nucleosynthesis of Binary-stripped Stars

The cosmic origin of the elements, the fundamental chemical building blocks of the universe, is still uncertain. Binary interactions play a key role in the evolution of many massive stars, yet their impact on chemical yields is poorly understood. Using the MESA stellar evolution code, we predict the chemical yields ejected in wind mass loss and the supernovae of single and binary-stripped stars. We do this with a large 162-isotope nuclear network at solar metallicity. We find that binary-stripped stars are more effective producers of the elements than single stars, due to their increased mass loss and an increased chance to eject their envelopes during a supernova. This increased production by binaries varies across the periodic table, with F and K being more significantly produced by binary-stripped stars than single stars. We find that the 12C/13C could be used as an indicator of the conservativeness of mass transfer, as 13C is preferentially ejected during mass transfer while 12C is preferentially ejected during wind mass loss. We identify a number of gamma-ray-emitting radioactive isotopes that may be used to help constrain progenitor and explosion models of core-collapse supernovae with next-generation gamma-ray detectors. For single stars we find that 44V and 52Mn are strong probes of the explosion model, while for binary-stripped stars it is 48Cr. Our findings highlight that binary-stripped stars are not equivalent to two single stars and that detailed stellar modeling is needed to predict their final nucleosynthetic yields.

Core collapse supernova gravitational wave emission for progenitors of 9.6, 15, and 25M⊙

We present gravitational wave emission predictions based on three core collapse supernova simulations corresponding to three different progenitor masses. The masses span a large range, between 9.6 and $25{M}_{\ensuremath{\bigodot}}$, are all initially nonrotating, and are of two metallicities: zero and solar. We compute both the temporal evolution of the gravitational wave strains for both the plus and the cross polarizations, as well as their spectral decomposition and characteristic strains. The temporal evolution of our zero metallicity $9.6{M}_{\ensuremath{\bigodot}}$ progenitor model is distinct from the temporal evolution of our solar metallicity $15{M}_{\ensuremath{\bigodot}}$ progenitor model and our zero metallicity $25{M}_{\ensuremath{\bigodot}}$ progenitor model. In the former case, the high-frequency gravitational wave emission is largely confined to a brief time period $\ensuremath{\sim}75\text{ }\text{ }\mathrm{ms}$ after bounce, whereas in the latter two cases high-frequency emission does not commence until $\ensuremath{\sim}125\text{ }\text{ }\mathrm{ms}$ after bounce or later. The excitation mechanisms of the high-frequency emission in all three cases correspond to proto-neutron star convection and accretion onto the proto-neutron star from the convective gain layer above it, with the former playing the dominant role for most of the evolution. The low-frequency emission in all three models exhibits very similar behavior. At frequencies below $\ensuremath{\sim}250\text{ }\text{ }\mathrm{Hz}$, gravitational waves are emitted by neutrino-driven convection and the standing accretion shock instability (SASI). This emission extends throughout the simulations when a gain region is present. In all three models, explosion is observed at $\ensuremath{\sim}125$, $\ensuremath{\sim}500$, and $\ensuremath{\sim}250\text{ }\text{ }\mathrm{ms}$ after bounce in the 9.6, 15, and $25{M}_{\ensuremath{\bigodot}}$ progenitor models, respectively. At these times, the low-frequency gravitational wave emission is joined by very low-frequency emission, below $\ensuremath{\sim}10\text{ }\text{ }\mathrm{Hz}$. These very low-frequency episodes are the result of explosion and begin at the above designated explosion times in each of our models. Our characteristic strains tell us that, in principle, all three gravitational wave signals would be detectable by current-generation detectors for a supernova at a distance of 10 kpc. However, our $9.6{M}_{\ensuremath{\bigodot}}$ progenitor model is a significantly weaker source of gravitational waves, with strain amplitudes approximately 5--10 times less than in our other two models. The characteristic strain for this model tells us that such a supernova would be detectable only within a much more narrow frequency range around the maximum sensitivity of today's detectors. Finally, in our $9.6{M}_{\ensuremath{\bigodot}}$ progenitor model, we see very high-frequency gravitational radiation, extending up to $\ensuremath{\sim}2000\text{ }\text{ }\mathrm{Hz}$. This feature results from the interaction of shock- and deleptonization-induced convection with perturbations introduced in the progenitor by nuclear burning during core collapse. While unique to the $9.6{M}_{\ensuremath{\bigodot}}$ progenitor model analyzed here, this very high-frequency emission may, in fact, be a generic feature of the predictions for the gravitational wave emission from all core collapse supernova models when simulations are performed with three-dimensional progenitors.

New Constraints for Supernova Models from Presolar Silicon Carbide X Grains with Very High 26Al/27Al Ratios

We report C, N, Mg-Al, Si, and S isotope data of six 1–3 μm-sized SiC grains of Type X from the Murchison CM2 chondrite, believed to have formed in the ejecta of core-collapse supernova (CCSN) explosions. Their C, N, and Si isotopic compositions are fully compatible with previously studied X grains. Magnesium is essentially monoisotopic 26Mg which gives clear evidence for the decay of radioactive 26Al. Inferred initial 26Al/27Al ratios are between 0.6 and 0.78 which is at the upper end of previously observed ratios of X grains. Contamination with terrestrial or solar system Al apparently is low or absent, which makes the X grains from this study particularly interesting and useful for a quantitative comparison of Al isotope data with predictions from supernova models. The consistently high 26Al/27Al ratios observed here may suggest that the lower 26Al/27Al ratios of many X grains from the literature are the result of significant Al contamination and in part also of an improper quantification of 26Al. The real dispersion of 26Al/27Al ratios in X grains needs to be explored by future studies. The high observed 26Al/27Al ratios in this work provide a crucial constraint for the production of 26Al in CCSN models. We explored different CCSN models, including both “classical” and H ingestion CCSN models. It is found that the classical models cannot account for the high 26Al/27Al ratios observed here; in contrast, H ingestion models are able to reproduce the 26Al/27Al ratios along with C, N, and Si isotopic ratios reasonably well.

Explosions in the cosmos: core-collapse supernovae

Explosions in the cosmos: core-collapse supernovae This is a project supported by OTKA and by the ERC Consolidator Grant project RADIOSTAR. The aim is to shed light on the production of proton-rich heavy nuclei that are measured in the solar system abundance pattern, using different sets of core-collapse supernova models and the latest nuclear reaction rates.

MRI-driven αΩ dynamos in protoneutron stars

Context. Magnetars are highly magnetized neutron stars that can produce a wide diversity of X-ray and soft gamma-ray emissions that are powered by magnetic dissipation. Their magnetic dipole is constrained in the range of 1014–1015 G by the measurement of their spin-down. In addition to fast rotation, these strong fields are also invoked to explain extreme stellar explosions, such as hypernovae, which are associated with long gamma-ray bursts and superluminous supernovae. A promising mechanism for explaining magnetar formation is the amplification of the magnetic field by the magnetorotational instability (MRI) in fast-rotating protoneutron stars (PNS). This scenario is supported by recent global incompressible models, which showed that a dipole field with magnetar-like intensity can be generated from small-scale turbulence. However, the impact of important physical ingredients, such as buoyancy and density stratification, on the efficiency of the MRI in generating a dipole field is still unknown. Aims. We assess the impact of the density and entropy profiles on the MRI dynamo in a global model of a fast-rotating PNS. The model focuses on the outer stratified region of the PNS that is stable to convection. Methods. Using the pseudo-spectral code MagIC, we performed 3D Boussinesq and anelastic magnetohydrodynamics simulations in spherical geometry with explicit diffusivities and with differential rotation forced at the outer boundary. The thermodynamic background of the anelastic models was retrieved from the data of 1D core-collapse supernova simulations from the Garching group. We performed a parameter study in which we investigated the influence of different approximations and the effect of the thermal diffusion through the Prandtl number. Results. We obtain a self-sustained turbulent MRI-driven dynamo. This confirms most of our previous incompressible results when they are rescaled for density. The MRI generates a strong turbulent magnetic field and a nondominant equatorial dipole, which represents about 4.3% of the averaged magnetic field strength. Interestingly, an axisymmetric magnetic field at large scales is observed to oscillate with time, which can be described as a mean-field αΩ dynamo. By comparing these results with models without buoyancy or density stratification, we find that the key ingredient explaining the appearance of this mean-field behavior is the density gradient. Buoyancy due to the entropy gradient damps turbulence in the equatorial plane, but it has a relatively weak influence in the low Prandtl number regime overall, as expected from neutrino diffusion. However, the buoyancy starts to strongly impact the MRI dynamo for Prandtl numbers close to unity. Conclusions. Our results support the hypothesis that the MRI is able to generate magnetar-like large-scale magnetic fields. The results furthermore predict the presence of a αΩ dynamo in the protoneutron star, which could be important to model in-situ magnetic field amplification in global models of core-collapse supernovae or binary neutron star mergers.

Prospects for reconstructing the gravitational-wave signals from core-collapse supernovae with Advanced LIGO-Virgo and the BayesWave algorithm

Our current understanding of the core-collapse supernova explosion mechanism is incomplete, with multiple viable models for how the initial shock wave might be energized enough to lead to a successful explosion. Detection of a gravitational-wave signal emitted in the initial few seconds after stellar core-collapse would provide unique and crucial insight into this process. With the Advanced LIGO and Advanced Virgo detectors expected to approach their design sensitivities soon, we could potentially detect this signal from a supernova within our galaxy. In anticipation of such a scenario, we study how well the BayesWave algorithm can recover the gravitational-wave signal from core-collapse supernova models in simulated advanced detector noise, and optimize its ability to accurately reconstruct the signal waveforms. We find that BayesWave can confidently reconstruct the signal from a range of supernova explosion models in Advanced LIGO-Virgo for network signal-to-noise ratios $\gtrsim 30$, reaching maximum reconstruction accuracies of $\sim 90\%$ at SNR $\sim 100$. For low SNR signals that are not confidently recovered, our optimization efforts result in gains in reconstruction accuracy of up to $20-40\%$, with typical gains of $\sim 10\%$.

Magnetic support for neutrino-driven explosion of 3D non-rotating core-collapse supernova models

ABSTRACT The impact of the magnetic field on post-bounce supernova dynamics of non-rotating stellar cores is studied by performing 3D magnetohydrodynamics simulations with spectral neutrino transport. The explodability of strongly and weakly magnetized models of 20 and 27 M⊙ pre-supernova progenitors are compared. We find that although the efficiency for the conversion of the neutrino heating into turbulent energy including magnetic fields in the gain region is not significantly different between the strong and weak field models, the amplified magnetic field due to the neutrino-driven convection on large hot bubbles just behind stalled shock results in a faster and more energetic explosion in the strongly magnetized models. In addition, by comparing the difference between the 2nd- and 5th-order spatial accuracy of the simulation in the strong field model for 27 M⊙ progenitor, we also find that the higher order accuracy in space is beneficial to the explosion because it enhances the growth of neutrino-driven convection in the gain region. Based on our results of core-collapse supernova simulations for the non-rotating model, a new possibility for the origin of the magnetic field of the protoneutron star (PNS) is proposed. The magnetic field is accumulated and amplified to magnetar level, that is, $\mathcal {O}(10^{14})$ G, in the convectively stable shell near the PNS surface.

Utilizing Gaussian mixture models in all-sky searches for short-duration gravitational wave bursts

Coherent WaveBurst is a generic, multidetector gravitational wave burst search based on the excess power approach. The coherent WaveBurst algorithm currently employed in the all-sky short-duration gravitational wave burst search uses a conditional approach on selected attributes in the multidimensional event attribute space to distinguish between noisy events from that of astrophysical origin. We have been developing a supervised machine learning approach based on the Gaussian mixture modeling to model the attribute space for signals as well as noise events to enhance the probability of burst detection [Gayathri et al.Phys. Rev. D 102, 104023 (2020)]. We further extend the GMM approach to the all-sky short-duration coherent WaveBurst search as a postprocessing step on events from the first half of the third observing run (O3a). We show an improvement in sensitivity to generic gravitational wave burst signal morphologies as well as the astrophysical source such as core-collapse supernova models due to the application of our Gaussian mixture model approach to coherent WaveBurst triggers. The Gaussian mixture model method recovers the gravitational wave signals from massive compact binary coalescences identified by coherent WaveBurst targeted for binary black holes in GWTC-2, with better significance than the all-sky coherent WaveBurst search. No additional significant gravitational wave bursts are observed.

Spatially Resolved X-Ray Study of Supernova Remnant G306.3–0.9 with Unusually High Calcium Abundance

Abstract G306.3–0.9 is an asymmetric Galactic supernova remnant (SNR), whose progenitor has been thought to be a Type Ia supernova (SN), but its high Ca abundance appears inconsistent with the Type Ia origin. Hoping to uncover the reason for its asymmetry and the origin of this SNR, we performed a spatially resolved X-ray spectroscopic analysis of XMM-Newton and Chandra observation data. We divided the SNR into 13 regions and analyzed the spectra using two-temperature models (0.2 keV + 1 keV). Compared to the southwestern regions, the northeastern regions have higher metal abundances and a lower gas density. This suggests that the asymmetric morphology results from the nonuniform ambient environment. We found that neither Type Ia nor core-collapse SN models can account for the abnormally high abundance ratios of Ar/Si, Ca/Si, or the shape of the abundance curve. A comparison with the Ca-rich transient models suggests that G306.3–0.9 is likely to be the first identified Galactic Ca-rich transient remnant, although the theoretical production of element S is lower. We also note that the conclusion for the SNR’s origin relies on the measured abundance ratios and existing nucleosynthesis models. Between two groups of Ca-rich transient explosion models, we prefer the He shell detonation for an accreting white dwarf, rather than the merger of a white dwarf and a neutron star.

Gravitational wave signals from 2D core–collapse supernova models with rotation and magnetic fields

ABSTRACT We investigate the impact of rotation and magnetic fields on the dynamics and gravitational wave emission in 2D core–collapse supernova simulations with neutrino transport. We simulate 17 different models of $15\, {\rm M}_\odot$ and $39\, {\rm M}_\odot$ progenitor stars with various initial rotation profiles and initial magnetic fields strengths up to $10^{12}\, \mathrm{G}$, assuming a dipolar field geometry in the progenitor. Strong magnetic fields generally prove conducive to shock revival, though this trend is not without exceptions. The impact of rotation on the post-bounce dynamics is more variegated, in line with previous studies. A significant impact on the time-frequency structure of the gravitational wave signal is found only for rapid rotation or strong initial fields. For rapid rotation, the angular momentum gradient at the proto-neutron star surface can appreciably affect the frequency of the dominant mode, so that known analytic relations for the high-frequency emission band no longer hold. In case of two magnetorotational explosion models, the deviation from these analytic relations is even more pronounced. One of the magnetorotational explosions has been evolved to more than half a second after the onset of the explosion and shows a subsidence of high-frequency emission at late times. Its most conspicuous gravitational wave signature is a high-amplitude tail signal. We also estimate the maximum detection distances for our waveforms. The magnetorotational models do not stick out for higher detectability during the post-bounce and explosion phase.

All-sky search for short gravitational-wave bursts in the third Advanced LIGO and Advanced Virgo run

This paper presents the results of a search for generic short-duration gravitational-wave transients in data from the third observing run of Advanced LIGO and Advanced Virgo. Transients with durations of milliseconds to a few seconds in the 24–4096 Hz frequency band are targeted by the search, with no assumptions made regarding the incoming signal direction, polarization, or morphology. Gravitational waves from compact binary coalescences that have been identified by other targeted analyses are detected, but no statistically significant evidence for other gravitational wave bursts is found. Sensitivities to a variety of signals are presented. These include updated upper limits on the source rate density as a function of the characteristic frequency of the signal, which are roughly an order of magnitude better than previous upper limits. This search is sensitive to sources radiating as little as ∼10−10 M⊙c2 in gravitational waves at ∼70 Hz from a distance of 10 kpc, with 50% detection efficiency at a false alarm rate of one per century. The sensitivity of this search to two plausible astrophysical sources is estimated: neutron star f modes, which may be excited by pulsar glitches, as well as selected core-collapse supernova models.1 MoreReceived 8 July 2021Accepted 12 October 2021DOI:https://doi.org/10.1103/PhysRevD.104.122004© 2021 American Physical SocietyPhysics Subject Headings (PhySH)Research AreasGravitational wave sourcesGravitational wavesGravitation, Cosmology & Astrophysics

Core-collapse Supernovae: From Neutrino-driven 1D Explosions to Light Curves and Spectra

Abstract We present bolometric and broadband light curves and spectra for a suite of core-collapse supernova models exploded self-consistently in spherical symmetry within the PUSH framework. We analyze broad trends in these light curves and categorize them based on morphology. We find that these morphological categories relate simply to the progenitor radius and mass of the hydrogen envelope. We present a proof-of-concept sensitive-variable analysis, indicating that an important determining factor in the properties of a light curve within a given category is 56Ni mass. We follow spectra from the photospheric to the nebular phase. These spectra show characteristic iron-line blanketing at short wavelengths and Doppler-shifted Fe ii and Ti ii absorption lines. To enable this analysis, we develop a first-of-its-kind pipeline from a massive progenitor model, through a self-consistent explosion in spherical symmetry, to electromagnetic counterparts. This opens the door to more detailed analyses of the collective properties of these observables. We provide a machine-readable database of our light curves and spectra online at go.ncsu.edu/astrodata.