Three-dimensional simulations of neutrino-driven core-collapse supernovae from low-mass single and binary star progenitors
11 Apr 2019-Monthly Notices of the Royal Astronomical Society (Oxford Academic)-Vol. 484, Iss: 3, pp 3307-3324
TL;DR: In this paper, a suite of seven 3D supernova simulations of non-rotating low-mass progenitors using multi-group neutrino transport is presented, where the mass outflow rate already exceeds the accretion rate onto the proto-neutron star and the mass and angular momentum of the compact remnant have closely approached their final value, barring the possibility of later fallback.
Abstract: We present a suite of seven 3D supernova simulations of non-rotating low-mass progenitors using multi-group neutrino transport. Our simulations cover single star progenitors with zero-age main sequence masses between $9.6 M_\odot$ and $12.5 M_\odot$ and (ultra)stripped-envelope progenitors with initial helium core masses between $2.8 M_\odot$ and $3.5 M_\odot$. We find explosion energies between $0.1\,\mathrm{Bethe}$ and $0.4\,\mathrm{Bethe}$, which are still rising by the end of the simulations. Although less energetic than typical events, our models are compatible with observations of less energetic explosions of low-mass progenitors. In six of our models, the mass outflow rate already exceeds the accretion rate onto the proto-neutron star, and the mass and angular momentum of the compact remnant have closely approached their final value, barring the possibility of later fallback. While the proto-neutron star is still accelerated by the gravitational tug of the asymmetric ejecta, the acceleration can be extrapolated to obtain estimates for the final kick velocity. We obtain gravitational neutron star masses between $1.22 M_\odot$ and $1.44 M_\odot$, kick velocities between $11\, \mathrm{km}\, \mathrm{s}^{-1}$ and $695\, \mathrm{km}\, \mathrm{s}^{-1}$, and spin periods from $20\, \mathrm{ms}$ to $2.7\,\mathrm{s}$, which suggests that typical neutron star birth properties can be naturally obtained in the neutrino-driven paradigm. We find a loose correlation between the explosion energy and the kick velocity. There is no indication of spin-kick alignment, but a correlation between the kick velocity and the neutron star angular momentum, which needs to be investigated further as a potential point of tension between models and observations.
Citations
More filters
TL;DR: In 2019, the LIGO Livingston detector observed a compact binary coalescence with signal-to-noise ratio 12.9 and the Virgo detector was also taking data that did not contribute to detection due to a low SINR but were used for subsequent parameter estimation as discussed by the authors.
Abstract: On 2019 April 25, the LIGO Livingston detector observed a compact binary coalescence with signal-to-noise ratio 12.9. The Virgo detector was also taking data that did not contribute to detection due to a low signal-to-noise ratio, but were used for subsequent parameter estimation. The 90% credible intervals for the component masses range from to if we restrict the dimensionless component spin magnitudes to be smaller than 0.05). These mass parameters are consistent with the individual binary components being neutron stars. However, both the source-frame chirp mass and the total mass of this system are significantly larger than those of any other known binary neutron star (BNS) system. The possibility that one or both binary components of the system are black holes cannot be ruled out from gravitational-wave data. We discuss possible origins of the system based on its inconsistency with the known Galactic BNS population. Under the assumption that the signal was produced by a BNS coalescence, the local rate of neutron star mergers is updated to 250-2810.
1,189 citations
Cites background from "Three-dimensional simulations of ne..."
...It is also difficult to form light NS with masses below ∼ 1.2M in current supernova explosion simulations (Burrows et al. 2019; Müller et al. 2019)....
[...]
TL;DR: In this article, the authors present 39 candidate gravitational wave events from compact binary coalescences detected by Advanced LIGO and Advanced Virgo in the first half of the third observing run (O3a) between 1 April 2019 15:00 UTC and 1 October 2019 15.00.
Abstract: We report on gravitational wave discoveries from compact binary coalescences detected by Advanced LIGO and Advanced Virgo in the first half of the third observing run (O3a) between 1 April 2019 15:00 UTC and 1 October 2019 15:00. By imposing a false-alarm-rate threshold of two per year in each of the four search pipelines that constitute our search, we present 39 candidate gravitational wave events. At this threshold, we expect a contamination fraction of less than 10%. Of these, 26 candidate events were reported previously in near real-time through GCN Notices and Circulars; 13 are reported here for the first time. The catalog contains events whose sources are black hole binary mergers up to a redshift of ~0.8, as well as events whose components could not be unambiguously identified as black holes or neutron stars. For the latter group, we are unable to determine the nature based on estimates of the component masses and spins from gravitational wave data alone. The range of candidate events which are unambiguously identified as binary black holes (both objects $\geq 3~M_\odot$) is increased compared to GWTC-1, with total masses from $\sim 14~M_\odot$ for GW190924_021846 to $\sim 150~M_\odot$ for GW190521. For the first time, this catalog includes binary systems with significantly asymmetric mass ratios, which had not been observed in data taken before April 2019. We also find that 11 of the 39 events detected since April 2019 have positive effective inspiral spins under our default prior (at 90% credibility), while none exhibit negative effective inspiral spin. Given the increased sensitivity of Advanced LIGO and Advanced Virgo, the detection of 39 candidate events in ~26 weeks of data (~1.5 per week) is consistent with GWTC-1.
768 citations
TL;DR: In this article, the authors provide an answer to the question "How Were the Elements from Iron to Uranium Made?" (Abridged) by combining new results and important breakthroughs in the related nuclear, atomic and astronomical fields of science.
Abstract: The production of about half of the heavy elements found in nature is assigned to a specific astrophysical nucleosynthesis process: the rapid neutron capture process (r-process). Although this idea has been postulated more than six decades ago, the full understanding faces two types of uncertainties/open questions: (a) The nucleosynthesis path in the nuclear chart runs close to the neutron-drip line, where presently only limited experimental information is available, and one has to rely strongly on theoretical predictions for nuclear properties. (b) While for many years the occurrence of the r-process has been associated with supernovae, more recent studies have cast substantial doubts on this environment. Alternative scenarios include the mergers of neutron stars, neutron-star black hole mergers, but possibly also rare classes of supernovae as well as hypernovae/collapsars with polar jet ejecta and also accretion disk outflows related to the collapse of fast rotating massive stars with high magnetic fields. Stellar r-process abundance observations, have provided insights into, and constraints on the frequency of and conditions in the responsible stellar production sites. One of them, neutron star mergers, was just identified and related to the Gravitational Wave event GW170817. High resolution observations, increasingly more precise due to improved experimental atomic data, have been particularly important in defining the heavy element abundance patterns of the old halo stars, and thus determining the extent, and nature, of the earliest nucleosynthesis in our Galaxy. Combining new results and important breakthroughs in the related nuclear, atomic and astronomical fields of science, this review attempts to provide an answer to the question "How Were the Elements from Iron to Uranium Made?" (Abridged)
321 citations
TL;DR: In this article, the authors conducted 19 state-of-the-art 3D core-collapse supernova simulations spanning a broad range of progenitor masses and found that while the majority of these models explode, not all do, and that even models in the middle of the available proggenitor mass range may be less explodable.
Abstract: We have conducted nineteen state-of-the-art 3D core-collapse supernova simulations spanning a broad range of progenitor masses. This is the largest collection of sophisticated 3D supernova simulations ever performed. We have found that while the majority of these models explode, not all do, and that even models in the middle of the available progenitor mass range may be less explodable. This does not mean that those models for which we did not witness explosion would not explode in Nature, but that they are less prone to explosion than others. One consequence is that the "compactness" measure is not a metric for explodability. We find that lower-mass massive star progenitors likely experience lower-energy explosions, while the higher-mass massive stars likely experience higher-energy explosions. Moreover, most 3D explosions have a dominant dipole morphology, have a pinched, wasp-waist structure, and experience simultaneous accretion and explosion. We reproduce the general range of residual neutron-star masses inferred for the galactic neutron-star population. The most massive progenitor models, however, in particular vis a vis explosion energy, need to be continued for longer physical times to asymptote to their final states. We find that while the majority of the inner ejecta have Y$_e = 0.5$, there is a substantial proton-rich tail. This result has important implications for the nucleosynthetic yields as a function of progenitor. Finally, we find that the non-exploding models eventually evolve into compact inner configurations that experience a quasi-periodic spiral SASI mode. We otherwise see little evidence of the SASI in the exploding models.
175 citations
TL;DR: The delayed neutrino-heating mechanism is emerging as the key driver of supernova explosions, but there remain many issues to address, such as the chaos of the involved dynamics.
Abstract: Most supernova explosions accompany the death of a massive star. These explosions give birth to neutron stars and black holes, and eject solar masses of heavy elements. However, determining the mechanism of explosion has been a half-century journey of great numerical and physical complexity. Here we present the status of this theoretical quest and the physics and astrophysics upon which its resolution seems to depend. The delayed neutrino-heating mechanism is emerging as the key driver of supernova explosions, but there remain many issues to address, such as the chaos of the involved dynamics.
120 citations
References
More filters
TL;DR: This work recognizes the need for additional dissipation in any higher-order Godunov method of this type, and introduces it in such a way so as not to degrade the quality of the results.
3,892 citations
"Three-dimensional simulations of ne..." refers methods in this paper
...Recent updates to the code include 6th-order extremum preserving reconstruction (Colella & Sekora 2008; Sekora & Colella 2009) instead of the original 4th-order piecewise parabolic method of Colella & Woodward (1984) and an alternative scheme for filtering near the polar axis....
[...]
Book•
01 Jan 2001
3,513 citations
"Three-dimensional simulations of ne..." refers background in this paper
...Yet another alternative consists in filtering the solution in Fourier space, which has a long tradition in meteorology (Boyd 2001)....
[...]
TL;DR: Fromm's second-order scheme for integrating the linear convection equation is made monotonic through the inclusion of nonlinear feedback terms in this paper, where care is taken to keep the scheme in conservation form.
2,005 citations
"Three-dimensional simulations of ne..." refers methods in this paper
...In principle, higher-order conservative reconstruction can be used for this purpose, but we opt for piecewise linear reconstruction (typically using van Leer’s harmonic limiter; van Leer 1974)....
[...]
TL;DR: In this article, Buchdahl and Tolman showed that the moment of inertia and the binding energy of a neutron star are nearly universal functions of the star's compactness, which can be understood by considering two analytic, yet realistic, solutions of Einstein's equations.
Abstract: The structure of neutron stars is considered from theoretical and observational perspectives We demonstrate an important aspect of neutron star structure: the neutron star radius is primarily determined by the behavior of the pressure of matter in the vicinity of nuclear matter equilibrium density In the event that extreme softening does not occur at these densities, the radius is virtually independent of the mass and is determined by the magnitude of the pressure For equations of state with extreme softening or those that are self-bound, the radius is more sensitive to the mass Our results show that in the absence of extreme softening, a measurement of the radius of a neutron star more accurate than about 1 km will usefully constrain the equation of state We also show that the pressure near nuclear matter density is primarily a function of the density dependence of the nuclear symmetry energy, while the nuclear incompressibility and skewness parameters play secondary roles In addition, we show that the moment of inertia and the binding energy of neutron stars, for a large class of equations of state, are nearly universal functions of the star's compactness These features can be understood by considering two analytic, yet realistic, solutions of Einstein's equations, by, respectively, Buchdahl and Tolman We deduce useful approximations for the fraction of the moment of inertia residing in the crust, which is a function of the stellar compactness and, in addition, the pressure at the core-crust interface
1,354 citations
Additional excerpts
...Approximately correcting for the binding energy of cold neutron stars following Lattimer & Prakash (2001) to obtain gravitational masses Mgrav, Mgrav ≈ Mby − 0.084M ( Mgrav M )2 , (3) we find values between 1.22M for z9....
[...]
TL;DR: In this paper, an equation of state for hot, dense matter is presented in a form that is sufficiently rapid to use directly in hydrodynamical simulations, for example, in stellar collapse calculations.
1,188 citations
"Three-dimensional simulations of ne..." refers methods in this paper
...At high densities, we use the nuclear equation of state of Lattimer & Swesty (1991) with a bulk incompressibility of K = 220 MeV, supplemented by a low-density equation of state for nuclei, nucleons, and lepton and photon radiation....
[...]