Engineering topics

Vartanyan, David (ORCID:0000000319389282)

Publications and source records attributed to Vartanyan, David (ORCID:0000000319389282).

9 records

DOE OSTI2024

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

Abstract 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.

Choi, Lyla (ORCID:0000000153860133)↗

DOE OSTI2024

Physical Correlations and Predictions Emerging from Modern Core-collapse Supernova Theory

Abstract In this paper, we derive correlations between core-collapse supernova observables and progenitor core structures that emerge from our suite of 20 state-of-the-art 3D core-collapse supernova simulations carried to late times. This is the largest such collection of 3D supernova models ever generated and allows one to witness and derive testable patterns that might otherwise be obscured when studying one or a few models in isolation. From this panoramic perspective, we have discovered correlations between explosion energy, neutron star gravitational birth masses, 56 Ni and α -rich freezeout yields, and pulsar kicks and theoretically important correlations with the compactness parameter of progenitor structure. We find a correlation between explosion energy and progenitor mantle binding energy, suggesting that such explosions are self-regulating. We also find a testable correlation between explosion energy and measures of explosion asymmetry, such as the ejecta energy and mass dipoles. While the correlations between two observables are roughly independent of the progenitor zero-age main-sequence (ZAMS) mass, the many correlations we derive with compactness cannot unambiguously be tied to a particular progenitor ZAMS mass. This relationship depends on the compactness/ZAMS mass mapping associated with the massive star progenitor models employed. Therefore, our derived correlations between compactness and observables may be more robust than with ZAMS mass but can nevertheless be used in the future once massive star modeling has converged.

Astronomy & Astrophysics↗

DOE OSTI2024

A Theory for Neutron Star and Black Hole Kicks and Induced Spins

Abstract Using 20 long-term 3D core-collapse supernova simulations, we find that lower compactness progenitors that explode quasi-spherically due to the short delay to explosion experience smaller neutron star recoil kicks in the ∼100−200 km s −1 range, while higher compactness progenitors that explode later and more aspherically leave neutron stars with kicks in the ∼300−1000 km s −1 range. In addition, we find that these two classes are correlated with the gravitational mass of the neutron star. This correlation suggests that the survival of binary neutron star systems may in part be due to their lower kick speeds. We also find a correlation between the kick and both the mass dipole of the ejecta and the explosion energy. Furthermore, one channel of black hole birth leaves masses of ∼10 M ⊙ , is not accompanied by a neutrino-driven explosion, and experiences small kicks. A second channel is through a vigorous explosion that leaves behind a black hole with a mass of ∼3.0 M ⊙ kicked to high speeds. We find that the induced spins of nascent neutron stars range from seconds to ∼10 ms, but do not yet see a significant spin/kick correlation for pulsars. We suggest that if an initial spin biases the explosion direction, a spin/kick correlation would be a common byproduct of the neutrino mechanism of core-collapse supernovae. Finally, the induced spin in explosive black hole formation is likely large and in the collapsar range. This new 3D model suite provides a greatly expanded perspective and appears to explain some observed pulsar properties by default.

Astronomy & Astrophysics↗

DOE OSTI2023

Characterizing the Directionality of Gravitational Wave Emission from Matter Motions within Core-collapse Supernovae

Abstract We analyze the directional dependence of the gravitational wave (GW) emission from 15 3D neutrino radiation hydrodynamic simulations of core-collapse supernovae (CCSNe). Using spin weighted spherical harmonics, we develop a new analytic technique to quantify the evolution of the distribution of GW emission over all angles. We construct a physics-informed toy model that can be used to approximate GW distributions for general ellipsoid-like systems, and use it to provide closed form expressions for the distribution of GWs for different CCSN phases. Using these toy models, we approximate the protoneutron star (PNS) dynamics during multiple CCSN stages and obtain similar GW distributions to simulation outputs. When considering all viewing angles, we apply this new technique to quantify the evolution of preferred directions of GW emission. For nonrotating cases, this dominant viewing angle drifts isotropically throughout the supernova, set by the dynamical timescale of the PNS. For rotating cases, during core bounce and the following tens of milliseconds, the strongest GW signal is observed along the equator. During the accretion phase, comparable—if not stronger—GW amplitudes are generated along the axis of rotation, which can be enhanced by the low T /∣ W ∣ instability. We show two dominant factors influencing the directionality of GW emission are the degree of initial rotation and explosion morphology. Lastly, looking forward, we note the sensitive interplay between GW detector site and supernova orientation, along with its effect on detecting individual polarization modes.

Astronomy & Astrophysics↗

DOE OSTI2023

Black Hole Formation Accompanied by the Supernova Explosion of a 40 M ⊙ Progenitor Star

Abstract We have simulated the collapse and evolution of the core of a solar-metallicity 40 M ⊙ star and find that it explodes vigorously by the neutrino mechanism, despite its very high “compactness.” Within ∼1.5 s of explosion, a black hole forms. The explosion is very asymmetrical and has a total explosion energy of ∼1.6 × 10 51 erg. At black hole formation, its baryon mass is ∼2.434 M ⊙ and gravitational mass is 2.286 M ⊙ . Seven seconds after black hole formation, an additional ∼0.2 M ⊙ is accreted, leaving a black hole baryon mass of ∼2.63 M ⊙ . A disk forms around the proto−neutron star, from which a pair of neutrino-driven jets emanates. These jets accelerate some of the matter up to speeds of ∼45,000 km s −1 and contain matter with entropies of ∼50. The large spatial asymmetry in the explosion results in a residual black hole recoil speed of ∼1000 km s −1 . This novel black hole formation channel now joins the other black hole formation channel between ∼12 and ∼15 M ⊙ discovered previously and implies that the black hole/neutron star birth ratio for solar-metallicity stars could be ∼20%. However, one channel leaves black holes in perhaps the ∼5–15 M ⊙ range with low kick speeds, while the other leaves black holes in perhaps the ∼2.5–3.0 M ⊙ mass range with high kick speeds. However, even ∼8.8 s after core bounce the newly formed black hole is still accreting at a rate of ∼2 × 10 −2 M ⊙ s −1 , and whether the black hole eventually achieves a significantly larger mass over time is yet to be determined.

Astronomy & Astrophysics↗

DOE OSTI2023

Neutrino signatures of 100 2D Axisymmetric Core-Collapse Supernova Simulations

ABSTRACT We present in this paper a public data release of an unprecedentedly large set of core-collapse supernova (CCSN) neutrino emission models, comprising 100 detailed 2D axisymmetric radiation-hydrodynamic simulations evolved out to as late as ∼5 s post-bounce and spanning an extensive range of massive-star progenitors. The motivation for this paper is to provide a physically and numerically uniform benchmark data set to the broader neutrino detection community to help it characterize and optimize subsurface facilities for what is likely to be a once-in-a-lifetime galactic supernova burst event. With this release, we hope to (1) help the international experiment and modelling communities more efficiently optimize the retrieval of physical information about the next galactic CCSN, (2) facilitate the better understanding of core-collapse theory and modelling among interested experimentalists, and (3) help further integrate the broader supernova neutrino community.

Astronomy & Astrophysics↗

DOE OSTI2022

The essential character of the neutrino mechanism of core-collapse supernova explosions

Abstract Calibrating with detailed 2D core-collapse supernova (CCSN) simulations, we derive a simple CCSN explosion condition based solely upon the terminal density profiles of state-of-the-art stellar evolution calculations of the progenitor massive stars. This condition captures the vast majority of the behaviour of the one hundred 2D state-of-the-art models we performed to gauge its usefulness. The goal is to predict, without resort to detailed simulation, the explodability of a given massive star. We find that the simple maximum fractional ram pressure jump discriminant we define works well ∼90 per cent of the time and we speculate on the origin of the few false positives and false negatives we witness. The maximum ram pressure jump generally occurs at the time of accretion of the silicon/oxygen interface, but not always. Our results depend upon the fidelity with which the current implementation of our code F ornax adheres to Nature and issues concerning the neutrino–matter interaction, the nuclear equation of state, the possible effects of neutrino oscillations, grid resolution, the possible role of rotation and magnetic fields, and the accuracy of the numerical algorithms employed remain to be resolved. Nevertheless, the explodability condition we obtain is simple to implement, shows promise that it might be further generalized while still employing data from only the unstable Chandrasekhar progenitors, and is a more credible and robust simple explosion predictor than can currently be found in the literature.

Astronomy & Astrophysics↗

DOE OSTI2022

On the Origin of Pulsar and Magnetar Magnetic Fields

Abstract In order to address the generation of neutron star magnetic fields, with particular focus on the dichotomy between magnetars and radio pulsars, we consider the properties of dynamos as inferred from other astrophysical systems. With sufficiently low (modified) Rossby number, convective dynamos are known to produce dipole-dominated fields whose strength scales with convective flux, and we argue that these expectations should apply to the convective protoneutron stars (PNSs) at the centers of core-collapse supernovae. We analyze a suite of three-dimensional simulations of core collapse, featuring a realistic equation of state and full neutrino transport, in this context. All our progenitor models, ranging from 9 M ⊙ to 25 M ⊙ , including one with initial rotation, have sufficiently vigorous PNS convection to generate dipole fields of order ∼10 15 Gauss, if the modified Rossby number resides in the critical range. Thus, the magnetar/radio pulsar dichotomy may arise naturally in part from the distribution of core rotation rates in massive stars.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

The collapse and three-dimensional explosion of three-dimensional massive-star supernova progenitor models

ABSTRACT The explosion outcome and diagnostics of core-collapse supernovae depend sensitively on the nature of the stellar progenitor, but most studies to date have focused exclusively on one-dimensional, spherically symmetric massive star progenitors. We present some of the first core-collapse supernovae simulations of three-dimensional massive star supernovae progenitors, a 12.5- and a 15-M⊙ model, evolved in three dimensions from collapse to bounce through explosion with the radiation-hydrodynamic code fornax. We compare the results using those starting from three-dimensional progenitors to three-dimensional simulations of spherically symmetric, one-dimensional progenitors of the same mass. We find that the models evolved in three dimensions during the final stages of massive star evolution are more prone to explosion. The turbulence arising in these multidimensional initial models serves as seed turbulence that promotes shock revival. Detection of gravitational waves and neutrinos signals could reveal signatures of pre-bounce turbulence.

79 ASTRONOMY AND ASTROPHYSICS↗