Shock Breakout Research Articles

Probing for the host galaxies of the fast X-ray transients XRT 000519 and XRT 110103

ABSTRACT Over the past few years, ∼30 extragalactic fast X-ray transients (FXRTs) have been discovered, mainly in Chandra and XMM-Newton data. Their nature remains unclear, with proposed origins, including a double neutron star merger, a tidal disruption event involving an intermediate-mass black hole and a white dwarf, or a supernova shock breakout. A decisive differentiation between these three promising mechanisms for their origin requires an understanding of the FXRT energetics, environments, and/or host properties. We present optical observations obtained with the Very Large Telescope for the FXRTs XRT 000519 and XRT 110103 and Gran Telescopio Canarias observations for XRT 000519 designed to search for host galaxies of these FXRTs. In the gs, rs, and R-band images, we detect an extended source on the north-west side of the $\sim \, 1^{\prime \prime }$ (68 per cent confidence) error circle of the X-ray position of XRT 000519 with a Kron magnitude of gs = 26.29 ± 0.09 (AB magnitude). We discuss the XRT 000519 association with the probable host candidate for various possible distances, and we conclude that if XRT 000519 is associated with the host candidate a supernova shock breakout scenario is likely excluded. No host galaxy is found near XRT 110103 down to a limiting magnitude of R > 25.8.

Supernova Shock Breakout/Emergence Detection Predictions for a Wide-field X-Ray Survey

There are currently many large-field surveys that are operational and are being planned including the powerful Vera C. Rubin Observatory Legacy Survey of Space and Time. These surveys will increase the number and diversity of transients dramatically. However, for some transients, like supernovae (SNe), we can gain more understanding by directed observations (e.g., shock breakout and γ-ray detections) than by simply increasing the sample size. For example, the initial emission from these transients can be a powerful probe of these explosions. Upcoming ground-based detectors are not ideally suited to observing the initial emission (shock emergence) of these transients. These observations require a large field-of-view X-ray mission with a UV follow-up within the first hour of shock breakout. The emission in the first 1 hr to even 1 day provides strong constraints on the stellar radius and asymmetries in the outer layers of stars, the properties of the circumstellar medium (e.g., inhomogeneities in the wind for core-collapse SNe and accreting companions in thermonuclear SNe), and the transition region between these two areas. This paper describes a simulation for the number of SNe that could be seen by a large field-of-view lobster-eye X-ray and UV observatory.

Numerical Simulations of Convective Three-dimensional Red Supergiant Envelopes

We explore the three-dimensional properties of convective, luminous (L ≈ 104.5–105 L ⊙), hydrogen-rich envelopes of red supergiants (RSGs) based on radiation hydrodynamic simulations in spherical geometry using Athena++. These computations comprise ≈30% of the stellar volume, include gas and radiation pressure, and self-consistently track the gravitational potential for the outer ≈3M ⊙ of the simulated M ≈ 15M ⊙ stars. This work reveals a radius, R corr, around which the nature of the convection changes. For r > R corr, though still optically thick, diffusion of photons dominates the energy transport. Such a regime is well studied in less luminous stars, but in RSGs, the near- (or above-)Eddington luminosity (due to opacity enhancements at ionization transitions) leads to the unusual outcome of denser regions moving outward rather than inward. This region of the star also has a large amount of turbulent pressure, yielding a density structure much more extended than 1D stellar evolution predicts. This “halo” of material will impact predictions for both shock breakout and early lightcurves of Type IIP supernovae. Inside of R corr, we find a nearly flat entropy profile as expected in the efficient regime of mixing-length theory (MLT). Radiation pressure provides ≈1/3 of the support against gravity in this region. Our comparisons to MLT suggest a mixing length of α = 3–4, consistent with the sizes of convective plumes seen in the simulations. The temporal variability of these 3D models is mostly on the timescale of the convective plume lifetimes (≈300 days), with amplitudes consistent with those observed photometrically.

Luminosity Function and Event Rate Density of XMM-Newton-selected Supernova Shock Breakout Candidates

Abstract A dozen X-ray supernova shock breakout (SN SBO) candidates were reported recently based on XMM-Newton archival data, which increased the X-ray-selected SN SBO sample by an order of magnitude. Assuming that they are genuine SN SBOs, we study the luminosity function (LF) by improving on the method used in our previous work. The light curves and the spectra of the candidates were used to derive the maximum volume within which these objects could be detected with XMM-Newton by simulation. The results show that the SN SBO LF can be described by either a broken power law (BPL) with indices (at the 68% confidence level) of 0.48 ± 0.28 and 2.11 ± 1.27 before and after the break luminosity at log ( L b / erg s − 1 ) = 45.32 ± 0.55 or a single power law (SPL) with an index of 0.80 ± 0.16. The local event rate densities of SN SBOs above 5 × 1042 erg s−1 are consistent for two models, i.e., 4.6 − 1.3 + 1.7 × 10 4 Gpc−3 yr−1 and 4.9 − 1.4 + 1.9 × 10 4 Gpc−3 yr−1 for BPL and SPL models, respectively. The number of fast X-ray transients of SN SBO origin can be significantly increased by wide-field X-ray telescopes such as the Einstein Probe.

Finding Fast Transients in Real Time Using a Novel Light-curve Analysis Algorithm

Abstract The current data acquisition rate of astronomical transient surveys and the promise for significantly higher rates in the next decade necessitate the development of novel approaches to analyze astronomical data sets and promptly detect objects of interest. The Deeper, Wider, Faster (DWF) program is a survey focused on the identification of fast-evolving transients, such as fast radio bursts, gamma-ray bursts, and supernova shock breakouts. It employs multifrequency simultaneous coverage of the same part of the sky over several orders of magnitude. Using the Dark Energy Camera mounted on the 4 m Blanco telescope, DWF captures a 20 s g-band exposure every minute, at a typical seeing of ∼1″ and an air mass of ∼1.5. These optical data are collected simultaneously with observations conducted over the entire electromagnetic spectrum—from radio to γ-rays—as well as cosmic-ray observations. In this paper, we present a novel real-time light-curve analysis algorithm, designed to detect transients in the DWF optical data; this algorithm functions independently from, or in conjunction with, image subtraction. We present a sample of fast transients detected by our algorithm, as well as a false-positive analysis. Our algorithm is customizable and can be tuned to be sensitive to transients evolving over different timescales and flux ranges.

AT 2018lqh and the Nature of the Emerging Population of Day-scale Duration Optical Transients

Abstract We report on the discovery of AT 2018lqh (ZTF 18abfzgpl)—a rapidly evolving extragalactic transient in a star-forming host at 242 Mpc. The transient g-band light curve’s duration above a half-maximum light is about 2.1 days, where 0.4/1.7 days are spent on the rise/decay, respectively. The estimated bolometric light curve of this object peaked at about 7 × 1042erg s−1—roughly 7 times brighter than the neutron star (NS)–NS merger event AT 2017gfo. We show that this event can be explained by an explosion with a fast (v ∼ 0.08 c) low-mass (≈0.07 M ⊙) ejecta, composed mostly of radioactive elements. For example, ejecta dominated by 56Ni with a timescale of t 0 ≅ 1.6 days for the ejecta to become optically thin for γ-rays fits the data well. Such a scenario requires burning at densities that are typically found in the envelopes of neutron stars or the cores of white dwarfs. A combination of circumstellar material (CSM) interaction power at early times and shock cooling at late times is consistent with the photometric observations, but the observed spectrum of the event may pose some challenges for this scenario. We argue that the observations are not consistent with a shock breakout from a stellar envelope, while a model involving a low-mass ejecta ramming into low-mass CSM cannot explain both the early- and late-time observations.

Wave-driven Mass Loss of Stripped Envelope Massive Stars: Progenitor-dependence, Mass Ejection, and Supernovae

Abstract The discovery of rapidly rising and fading supernovae powered by circumstellar interaction has suggested the pre-supernova mass eruption phase as a critical phenomenon in massive star evolution. It is important to understand the mass and radial extent of the circumstellar medium (CSM) from theoretically predicted mass ejection mechanisms. In this work, we study the wave heating process in massive hydrogen-poor stars, running a suite of stellar models in order to predict the wave energy and pre-explosion timescale of surface energy deposition. We survey stellar models with main-sequence progenitor masses from 20–70 M ⊙ and metallicity from 0.002–0.02. Most of these models predict that less than ∼1047 erg is deposited in the envelope, with the majority of the energy deposited in the last week of stellar evolution. This translates to CSM masses less than ∼10−2 M ⊙ that extend to less than ∼1014 cm, too small to greatly impact the light curves or spectra of the subsequent supernovae, except perhaps during the shock breakout phase. However, a few models predict somewhat higher wave energy fluxes, for which we perform hydrodynamical simulations of the mass ejection process. Radiative transfer simulations of the subsequent supernovae predict a bright but brief shock-cooling phase that could be detected in some Type Ib/c supernovae if they are discovered within a couple days of explosion.

Electromagnetic Precursors of Short Gamma-Ray Bursts as Counterparts of Gravitational Waves

Precursor emissions are found in some short gamma-ray bursts (SGRBs). In this paper, we review the theories and observations of the SGRB precursor and discuss its prospect as an electromagnetic counterpart of the gravitational wave event produced by neutron star (NS) mergers. The observed luminosity, spectrum, and duration of precursors are explained by the magnetospheric interaction model during the inspiral or the cocoon/jet shock breakout model during the jet propagation. In general, these two models predict that the precursor will be weaker than the main GRB, but will be of a larger opening angle, which makes it an advantageous gamma-ray counterpart for NS mergers in the local Universe, especially for NS - black hole mergers with very low mass ratios, in which the main GRBs are not expected. The joint observation of the precursor, SGRB, and gravitational wave will help to reveal the jet launch mechanism and post-merger remnant.

Three-dimensional Core-collapse Supernova Simulations with 160 Isotopic Species Evolved to Shock Breakout

Abstract We present three-dimensional simulations of core-collapse supernovae using the FLASH code that follow the progression of the explosion to the stellar surface, starting from neutrino radiation hydrodynamic simulations of the neutrino-driven phase performed with the Chimera code. We consider a 9.6 M ☉ zero-metallicity progenitor starting from both 2D and 3D Chimera models and a 10 M ☉ solar-metallicity progenitor starting from a 2D Chimera model, all simulated until shock breakout in 3D while tracking 160 nuclear species. The relative velocity difference between the supernova shock and the metal-rich Rayleigh–Taylor (R-T) “bullets” determines how the metal-rich ejecta evolves as it propagates through the density profile of the progenitor and dictates the final morphology of the explosion. We find maximum 56Ni velocities of ∼1950 and ∼1750 km s−1 at shock breakout from 2D and 3D 9.6 M ☉ Chimera models, respectively, due to the bullets’ ability to penetrate the He/H shell. When mapping from 2D, we find that the development of higher-velocity structures is suppressed when the 2D Chimera model and 3D FLASH model meshes are aligned. The development of faster-growing spherical-bubble structures, as opposed to the slower-growing toroidal structure imposed by axisymmetry, allows for interaction of the bullets with the shock and seeds further R-T instabilities at the He/H interface. We see similar effects in the 10 M ☉ model, which achieves maximum 56Ni velocities of ∼2500 km s−1 at shock breakout.

Direct Numerical Simulations of Cosmic-ray Acceleration at Dense Circumstellar Medium: Magnetic-field Amplification and Maximum Energy

Abstract Galactic cosmic rays are believed to be accelerated at supernova remnants. However, whether supernova remnants can be PeV is still very unclear. In this work we argue that PeV cosmic rays can be accelerated during the early phase of a supernova blast-wave expansion in dense red supergiant winds. We solve in spherical geometry a system combining a diffusive–convection equation that treats cosmic-ray dynamics coupled to magnetohydrodynamics to follow gas dynamics. A fast shock expanding in a dense ionized wind is able to trigger fast, non-resonant streaming instability over day timescales and energizes cosmic rays even under the effect of p–p losses. We find that such environments produce PeV blast waves, although the maximum energy depends on various parameters such as the injection rate and mass-loss rate of the winds. Multi-PeV energies can be reached if the progenitor mass-loss rates are of the order of 10−3 M ⊙ yr−1. It has been recently proposed that, prior to the explosion, hydrogen-rich massive stars can produce enhanced mass-loss rates. These enhanced rates would then favor the production of a PeV phase in early times after shock breakout.

Bolometric light curves of aspherical shock breakout

ABSTRACT The shock breakout emission is the first light that emerges from a supernova. In the spherical case, it is characterized by a brief UV flash. In an axisymmetric, non-spherical prolate explosion, the shock first breaches the surface along the symmetry axis, then peels around to larger angles, producing a breakout light curve which may differ substantially from the spherically symmetric case. We study the emergence of a non-relativistic, bipolar shock from a spherical star, and estimate the basic properties of the associated bolometric shock breakout signal. We identify four possible classes of breakout light curves, depending on the degree of asphericity. Compared to spherical breakouts, we find that the main distinguishing features of significantly aspherical breakouts are (1) a longer and fainter initial breakout flash and (2) an extended phase of slowly declining, or even rising, emission which is produced as ejecta flung out by the oblique breakout expand and cool. We find that the breakout flash has a maximum duration of roughly ∼R*/vbo, where R* is the stellar radius and vbo is the velocity of the fastest moving ejecta. For a standard Wolf−Rayet progenitor, the duration of the X-ray flash seen in SN 2008D exceeds this limit, and the same holds true for the prompt X-ray emission of low-luminosity GRBs such as GRB 060218. This suggests that these events cannot be explained by an aspherical explosion within a typical Wolf−Rayet star, implying that they originate from non-standard progenitors with larger breakout radii.

A Global Numerical Model of the Prompt Emission in Short Gamma-ray Bursts

Abstract We present the first global model of prompt emission from a short gamma-ray burst (GRB) that consistently describes the evolution of the central black hole (BH) torus system, the propagation of the jet through multicomponent merger ejecta, the transition into free expansion, and the photospheric emission from the relativistic jet. To this end, we perform a special relativistic neutrino-hydrodynamics simulation of a viscous BH-torus system, which is formed about 500 ms after the merger and is surrounded by dynamical ejecta as well as neutron star winds, along with a jet that is injected in the vicinity of the central BH. In a postprocessing step, we compute the photospheric emission using a relativistic Monte Carlo radiative transfer code. It is found that the wind from the torus leaves a strong imprint on the jet as well as on the emission, causing narrow collimation and rapid time variability. The dependence of the emission on viewing angle gives rise to correlations among the spectral peak energy, E p , isotropic energy, E iso, and peak luminosity, L p , which may provide natural explanations for the Amati and Yonetoku relations. We also find that the degree of polarization is small for emission from the jet core (≲2%), while it tends to increase with viewing angle outside the core and can become as high as ∼10%–40% for energies larger than the peak energy. Finally, the comparison of our model with GRB 170817A strongly disfavors the photospheric emission scenario and therefore supports alternative scenarios, such as cocoon shock breakout.

SN2017jgh: a high-cadence complete shock cooling light curve of a SN IIb with the Kepler telescope

ABSTRACT SN 2017jgh is a type IIb supernova discovered by Pan-STARRS during the C16/C17 campaigns of the Kepler/K2 mission. Here, we present the Kepler/K2 and ground based observations of SN 2017jgh, which captured the shock cooling of the progenitor shock breakout with an unprecedented cadence. This event presents a unique opportunity to investigate the progenitors of stripped envelope supernovae. By fitting analytical models to the SN 2017jgh light curve, we find that the progenitor of SN 2017jgh was likely a yellow supergiant with an envelope radius of $\sim 50{\!-\!}290\, {\rm R}_{\odot }$, and an envelope mass of $\sim 0{\!-\!}1.7\, {\rm M}_{\odot }$. SN 2017jgh likely had a shock velocity of ∼7500−10 300 km s−1. Additionally, we use the light curve of SN 2017jgh to investigate how early observations of the rise contribute to constraints on progenitor models. Fitting just the ground based observations, we find an envelope radius of $\sim 50{\!-\!}330\, {\rm R}_{\odot }$, an envelope mass of $\sim 0.3{\!-\!}1.7\, {\rm M}_{\odot }$ and a shock velocity of ∼9000−15 000 km s−1. Without the rise, the explosion time cannot be well constrained that leads to a systematic offset in the velocity parameter and larger uncertainties in the mass and radius. Therefore, it is likely that progenitor property estimates through these models may have larger systematic uncertainties than previously calculated.

Shock breakout in winds of red supergiants: Type IIP supernovae

AbstractWhen a supernova shockwave launched deep inside the star exits the surface, it probes the circumstellar medium established by prior mass loss from the pre supernova star. The bright electromagnetic display accompanying the shock breakout is influenced by the properties of the star and scripts the history of the stellar mass loss. We investigate with MESA and STELLA codes the radiative display resulting from a set of progenitors that we evolved to core collapse. We simulate with different internal convective overshoot and compositional mixing and two sets of mass loss schema, one the standard “Dutch” scheme and another, an enhanced, episodic mass loss at a late stage. Shock breakout from the star shows double peaked bolometric light curves for the Dutch wind, as well as high velocity ejecta accelerated during shock breakout. We contrast the breakout flash from an optically thick CSM with that of the rarified medium.

The Infrared Echo of SN2010jl and Its Implications for Shock Breakout Characteristics

Abstract SN 2010jl is a Type IIn core-collapse supernova whose radiative output is powered by the interaction of the supernova (SN) shock wave with its surrounding dense circumstellar medium (CSM). After day ∼60, its light curve developed a near-infrared (NIR) excess emission from dust. This excess could be a thermal IR echo from preexisting CSM dust, or emission from newly formed dust either in the cooling post-shock region of the CSM, or in the cooling SN ejecta. Recent analysis has shown that dust formation in the CSM can commence only after day ∼380, and has also ruled out newly formed ejecta dust as the source of the NIR emission. The early (<380 days) NIR emission can therefore only be attributed to an IR echo. The H–K color temperature of the echo is about 1250 K. The best-fitting model requires the presence of about 1.6 × 10−4 M ⊙ of amorphous carbon dust at a distance of 2.2 × 1016 cm from the explosion. The CSM-powered luminosity is preceded by an intense burst of hard radiation generated by the breakout of the SN shock through the stellar surface. The peak burst luminosity seen by the CSM dust is significantly reduced by Thomson scattering in the CSM, but still has the potential of evaporating the dust needed to produce the echo. We show that the survival of the echo-producing dust provides important constraints on the intensity, effective temperature, and duration of the burst.

Supernova explosions in active galactic nuclear discs

ABSTRACT Active galactic nuclei (AGNs) are prominent environments for stellar capture, growth, and formation. These environments may catalyse stellar mergers and explosive transients, such as thermonuclear and core-collapse supernovae (SNe). SN explosions in AGN discs generate strong shocks, leading to unique observable signatures. We develop an analytical model that follows the evolution of the shock propagating in the disc until it eventually breaks out. We derive the peak luminosity, bolometric light curve, and breakout time. The peak luminosities may exceed 1045 erg s−1 and last from hours to days. The brightest explosions occur in regions of reduced density: either off-plane, or in discs around low-mass central black holes (${\sim} 10^6\ \rm {M}_\odot$), or in starved subluminous AGNs. Explosions in the latter two sites are easier to observe due to a reduced AGN background luminosity. We perform suites of 1D Lagrangian radiative hydrodynamics snec code simulations to validate our results and obtain the luminosity in different bands, and 2D axisymmetric Eulerian hydrodynamics code hormone simulations to study the morphology of the ejecta and its deviation from spherical symmetry. The observed signature is expected to be a bright blue, UV or X-ray flare on top of the AGN luminosity from the initial shock breakout, while the subsequent red part of the light curve will largely be unobservable. We estimate the upper limit for the total event rate to be $\mathcal {R}\lesssim 100\ \rm yr^{-1}\ Gpc^{-3}$ for optimal conditions and discuss the large uncertainties in this estimate. Future high-cadence transient searches may reveal these events. Some existing tidal disruption event candidates may originate from AGN SNe.

The Role of Jet–Cocoon Mixing, Magnetization, and Shock Breakout in Neutrino and Cosmic-Ray Emission from Short Gamma-Ray Bursts

Abstract We perform general relativistic magnetohydrodynamic and relativistic magnetohydrodynamic simulations of weakly and highly magnetized gamma-ray burst (GRB) jets propagating in binary neutron star (BNS) merger ejecta. Using the simulations, we first find that mixing between the jet and cocoon, which is present in all types of jets, inhibits the formation of subphotospheric collisionless shocks. However, we show that a mild magnetization may lead to the formation of collisionless subshocks, which allow efficient proton acceleration. We consider shear acceleration and diffusive shock acceleration at collimation shocks, internal shocks, shock breakout, and external shocks to provide the first estimate for neutrino and cosmic-ray (CR) signals from self-consistent simulations of GRBs in BNS mergers. We find that short GRBs do not produce detectable neutrino signals with current-day facilities. Shock breakout yields ∼10 PeV neutrinos at viewing angles ∼20°, independent of the jet magnetization. However, a neutrino signal from shock breakout is well below the detection limits of current detectors. Such a signal would allow a coincident neutrino–γ-ray detection, providing a testable prediction for shock breakout as a neutrino production site. Using the numerical modeling that fits GW 170817 afterglow emission, we find that blast waves in BNS mergers can account for 5%–10% of the Galactic CR luminosity in the PeV–EeV energy range. Based on these estimates, the observed level of CR anisotropy places a constraint on the distance of the latest Galactic BNS merger to ≲3 kpc.

Free-surface velocity measurements of opaque materials in laser-driven shock-wave experiments using photonic Doppler velocimetry

We present a novel photonic Doppler velocimetry (PDV) design for laser-driven shock-wave experiments. This PDV design is intended to provide the capability of measuring the free-surface velocity of shocked opaque materials in the terapascal range. We present measurements of the free-surface velocity of gold for as long as ∼2 ns from the shock breakout, at pressures of up to ∼7 Mbar and a free-surface velocity of 7.3 km/s with an error of ∼1.5%. Such laboratory pressure conditions are achieved predominantly at high-intensity laser facilities where the only velocity diagnostic is usually line-imaging velocity interferometry for any reflector. However, that diagnostic is limited by the lower dynamic range of the streak camera (at a temporal resolution relevant to laser shock experiments) to measure the free-surface velocity of opaque materials up to pressures of only ∼1 Mbar. We expect the proposed PDV design to allow the free-surface velocity of opaque materials to be measured at much higher pressures.

Supernova 1987A: 3D Mixing and Light Curves for Explosion Models Based on Binary-merger Progenitors

Abstract Six binary-merger progenitors of supernova 1987A (SN 1987A) with properties close to those of the blue supergiant Sanduleak −69°202 are exploded by neutrino heating and evolved until long after shock breakout in 3D and continued for light-curve calculations in spherical symmetry. Our results confirm previous findings for single-star progenitors: (1) 3D neutrino-driven explosions with SN 1987A-like energies synthesize 56Ni masses consistent with the radioactive light-curve tail; (2) hydrodynamic models mix hydrogen inward to minimum velocities below 40 km s−1 compatible with spectral observations of SN 1987A; and (3) for given explosion energy the efficiency of outward radioactive 56Ni mixing depends mainly on high growth factors of Rayleigh–Taylor instabilities at the (C+O)/He and He/H composition interfaces and a weak interaction of fast plumes with the reverse shock occurring below the He/H interface. All binary-merger models possess presupernova radii matching the photometric radius of Sanduleak −69°202 and a structure of the outer layers allowing them to reproduce the observed initial luminosity peak in the first ∼7 days. Models that mix about 0.5 M ⊙ of hydrogen into the He-shell and exhibit strong outward mixing of 56Ni with maximum velocities exceeding the 3000 km s−1 observed for the bulk of ejected 56Ni have light-curve shapes in good agreement with the dome of the SN 1987A light curve. A comparative analysis of the best representatives of our 3D neutrino-driven explosion models of SN 1987A based on single-star and binary-merger progenitors reveals that only one binary model fulfills all observational constraints, except one.

Determining the Structure of Rotating Massive Stellar Cores with Gravitational Waves

Abstract The gravitational wave (GW) signal resulting from stellar core collapse encodes a wealth of information about the physical parameters of the progenitor star and the resulting core-collapse supernova (CCSN). We present a novel approach to constrain CCSN progenitor properties at collapse using two of the most detectable parts of the GW signal: the core-bounce signal and evolution of the dominant frequency mode from the protoneutron star. We focus on the period after core bounce but before explosion and investigate the predictive power of GWs from rotating CCSNe to constrain properties of the progenitor star. We analyze 34 2D and four 3D neutrino-radiation-hydrodynamic simulations of stellar core collapse in progenitors of varied initial mass and rotation rate. Extending previous work, we verify the compactness of the progenitor at collapse to correlate with the early ramp-up slope, and in rotating cases, also with the core angular momentum. Combining this information with the bounce signal, we present a new analysis method to constrain the pre-collapse core compactness of the progenitor. Because these GW features occur less than a second after core bounce, this analysis could allow astronomers to predict electromagnetic properties of a resulting CCSN even before shock breakout.