All ten LIGO/Virgo binary black hole (BH-BH) coalescences reported following the O1/O2 runs have near-zero effective spins. There are only three potential explanations for this. If the BH spin ...magnitudes are large, then: (i) either both BH spin vectors must be nearly in the orbital plane or (ii) the spin angular momenta of the BHs must be oppositely directed and similar in magnitude. Then there is also the possibility that (iii) the BH spin magnitudes are small. We consider the third hypothesis within the framework of the classical isolated binary evolution scenario of the BH-BH merger formation. We test three models of angular momentum transport in massive stars: a mildly efficient transport by meridional currents (as employed in the Geneva code), an efficient transport by the Tayler-Spruit magnetic dynamo (as implemented in the MESA code), and a very-efficient transport (as proposed by Fuller et al.) to calculate natal BH spins. We allow for binary evolution to increase the BH spins through accretion and account for the potential spin-up of stars through tidal interactions. Additionally, we update the calculations of the stellar-origin BH masses, including revisions to the history of star formation and to the chemical evolution across cosmic time. We find that we can simultaneously match the observed BH-BH merger rate density and BH masses and BH-BH effective spins. Models with efficient angular momentum transport are favored. The updated stellar-mass weighted gas-phase metallicity evolution now used in our models appears to be key for obtaining an improved reproduction of the LIGO/Virgo merger rate estimate. Mass losses during the pair-instability pulsation supernova phase are likely to be overestimated if the merger GW170729 hosts a BH more massive than 50
M
⊙
. We also estimate rates of black hole-neutron star (BH-NS) mergers from recent LIGO/Virgo observations. If, in fact. angular momentum transport in massive stars is efficient, then any (electromagnetic or gravitational wave) observation of a rapidly spinning BH would indicate either a very effective tidal spin up of the progenitor star (homogeneous evolution, high-mass X-ray binary formation through case A mass transfer, or a spin- up of a Wolf-Rayet star in a close binary by a close companion), significant mass accretion by the hole, or a BH formation through the merger of two or more BHs (in a dense stellar cluster).
ABSTRACT
We present 3D full-sphere supernova simulations of non-rotating low-mass (∼9 M⊙) progenitors, covering the entire evolution from core collapse through bounce and shock revival, through shock ...breakout from the stellar surface, until fallback is completed several days later. We obtain low-energy explosions (∼0.5–1.0 × 1050 erg) of iron-core progenitors at the low-mass end of the core-collapse supernova (LMCCSN) domain and compare to a super-AGB (sAGB) progenitor with an oxygen–neon–magnesium core that collapses and explodes as electron-capture supernova (ECSN). The onset of the explosion in the LMCCSN models is modelled self-consistently using the vertex-prometheus code, whereas the ECSN explosion is modelled using parametric neutrino transport in the prometheus-HOTB code, choosing different explosion energies in the range of previous self-consistent models. The sAGB and LMCCSN progenitors that share structural similarities have almost spherical explosions with little metal mixing into the hydrogen envelope. A LMCCSN with less second dredge-up results in a highly asymmetric explosion. It shows efficient mixing and dramatic shock deceleration in the extended hydrogen envelope. Both properties allow fast nickel plumes to catch up with the shock, leading to extreme shock deformation and aspherical shock breakout. Fallback masses of $\mathord {\lesssim }\, 5\, \mathord {\times }\, 10^{-3}$ M⊙ have no significant effects on the neutron star (NS) masses and kicks. The anisotropic fallback carries considerable angular momentum, however, and determines the spin of the newly born NS. The LMCCSN model with less second dredge-up results in a hydrodynamic and neutrino-induced NS kick of >40 km s−1 and a NS spin period of ∼30 ms, both not largely different from those of the Crab pulsar at birth.
It is difficult to establish the properties of massive stars that explode as supernovae. The electromagnetic emission during the first minutes to hours after the emergence of the shock from the ...stellar surface conveys important information about the final evolution and structure of the exploding star. However, the unpredictable nature of supernova events hinders the detection of this brief initial phase. Here we report the serendipitous discovery of a newly born, normal type IIb supernova (SN 2016gkg), which reveals a rapid brightening at optical wavelengths of about 40 magnitudes per day. The very frequent sampling of the observations allowed us to study in detail the outermost structure of the progenitor of the supernova and the physics of the emergence of the shock. We develop hydrodynamical models of the explosion that naturally account for the complete evolution of the supernova over distinct phases regulated by different physical processes. This result suggests that it is appropriate to decouple the treatment of the shock propagation from the unknown mechanism that triggers the explosion.
H and He features in photospheric spectra have seldom been used to infer quantitatively the properties of Type IIb, Ib and Ic supernovae (SNe IIb, Ib and Ic) and their progenitor stars. Most ...radiative transfer models ignored non-local thermodynamic equilibrium (NLTE) effects, which are extremely strong especially in the He-dominated zones. In this paper, a comprehensive set of model atmospheres for low-mass SNe IIb/Ib/Ic is presented. Long-standing questions, such as how much He can be contained in SNe Ic, where He lines are not seen, can thus be addressed. The state of H and He is computed in full NLTE, including the effect of heating by fast electrons. The models are constructed to represent iso-energetic explosions of the same stellar core with differently massive H/He envelopes on top. The synthetic spectra suggest that 0.06-0.14 M⊙ of He and even smaller amounts of H suffice for optical lines to be present, unless ejecta asymmetries play a major role. This strongly supports the conjecture that low-mass SNe Ic originate from binaries where progenitor mass loss can be extremely efficient.
Type Ia supernovae form an observationally uniform class of stellar explosions, in that more luminous objects have smaller decline-rates. This one-parameter behaviour allows type Ia supernovae to be ...calibrated as cosmological 'standard candles', and led to the discovery of an accelerating Universe. Recent investigations, however, have revealed that the true nature of type Ia supernovae is more complicated. Theoretically, it has been suggested that the initial thermonuclear sparks are ignited at an offset from the centre of the white-dwarf progenitor, possibly as a result of convection before the explosion. Observationally, the diversity seen in the spectral evolution of type Ia supernovae beyond the luminosity-decline-rate relation is an unresolved issue. Here we report that the spectral diversity is a consequence of random directions from which an asymmetric explosion is viewed. Our findings suggest that the spectral evolution diversity is no longer a concern when using type Ia supernovae as cosmological standard candles. Furthermore, this indicates that ignition at an offset from the centre is a generic feature of type Ia supernovae.
ABSTRACT
We report observations and analysis of the nearby gamma-ray burst GRB 161219B (redshift z = 0.1475) and the associated Type Ic supernova (SN) 2016jca. GRB 161219B had an isotropic gamma-ray ...energy of ∼1.6 × 1050 erg. Its afterglow is likely refreshed at an epoch preceding the first photometric points (0.6 d), which slows down the decay rates. Combined analysis of the SN light curve and multiwavelength observations of the afterglow suggest that the GRB jet was broad during the afterglow phase (full opening angle ∼42° ± 3°). Our spectral series shows broad absorption lines typical of GRB supernovae (SNe), which testify to the presence of material with velocities up to ∼0.25c. The spectrum at 3.73 d allows for the very early identification of an SN associated with a GRB. Reproducing it requires a large photospheric velocity ($35\, 000 \pm 7000$ km s−1). The kinetic energy of the SN is estimated through models to be Ekin≈4 × 1052 erg in spherical symmetry. The ejected mass in the explosion was Mej≈6.5 ± 1.5 M⊙, much less than that of other GRB-SNe, demonstrating diversity among these events. The total amount of 56Ni in the explosion was 0.27 ± 0.05 M⊙. The observed spectra require the presence of freshly synthesized 56Ni at the highest velocities, at least three times more than a standard GRB-SN. We also find evidence for a decreasing 56Ni abundance as a function of decreasing velocity. This suggests that SN 2016jca was a highly aspherical explosion viewed close to on-axis, powered by a compact remnant. Applying a typical correction for asymmetry, the energy of SN 2016jca was ∼(1–3) × 1052 erg, confirming that most of the energy produced by GRB-SNe goes into the kinetic energy of the SN ejecta.
When carbon is ignited off-centre in a CO core of a super-asymptotic giant branch star, its burning in a convective shell tends to propagate to the centre. Whether the C flame will actually be able ...to reach the centre depends on the efficiency of extra mixing beneath the C convective shell. Whereas thermohaline mixing is too inefficient to interfere with the C-flame propagation, convective boundary mixing can prevent the C burning from reaching the centre. As a result, a C-O-Ne white dwarf (WD) is formed, after the star has lost its envelope. Such a 'hybrid' WD has a small CO core surrounded by a thick ONe zone. In our 1D stellar evolution computations, the hybrid WD is allowed to accrete C-rich material, as if it were in a close binary system and accreted H-rich material from its companion with a sufficiently high rate at which the accreted H would be processed into He under stationary conditions, assuming that He could then be transformed into C. When the mass of the accreting WD approaches the Chandrasekhar limit, we find a series of convective Urca shell flashes associated with high abundances of 23Na and 25Mg. They are followed by off-centre C ignition leading to convection that occupies almost the entire star. To model the Urca processes, we use the most recent well-resolved data for their reaction and neutrino-energy loss rates. Because of the emphasized uncertainty of the convective Urca process in our hybrid WD models of Type Ia supernova (SN Ia) progenitors, we consider a number of their potentially possible alternative instances for different mixing assumptions, all of which reach a phase of explosive C ignition, either off or in the centre. Our hybrid SN Ia progenitor models have much lower C-to-O abundance ratios at the moment of the explosive C ignition than their pure CO counterparts, which may explain the observed diversity of the SNe Ia.
Hubble Space Telescope spectroscopic observations of the nearby Type Ia supernova (SN Ia) SN 2011fe, taken on 10 epochs from −13.1 to +40.8 d relative to B-band maximum light, and spanning the ...far-ultraviolet (UV) to the near-infrared (IR) are presented. This spectroscopic coverage makes SN 2011fe the best-studied local SN Ia to date. SN 2011fe is a typical moderately luminous SN Ia with no evidence for dust extinction. Its near-UV spectral properties are representative of a larger sample of local events (Maguire et al.). The near-UV to optical spectra of SN 2011fe are modelled with a Monte Carlo radiative transfer code using the technique of 'abundance tomography', constraining the density structure and the abundance stratification in the SN ejecta. SN 2011fe was a relatively weak explosion, with moderate Fe-group yields. The density structures of the classical model W7 and of a delayed detonation model were tested. Both have shortcomings. An ad hoc density distribution was developed which yields improved fits and is characterized by a high-velocity tail, which is absent in W7. However, this tail contains less mass than delayed detonation models. This improved model has a lower energy than one-dimensional explosion models matching typical SNe Ia (e.g. W7, WDD1; Iwamoto et al.). The derived Fe abundance in the outermost layer is consistent with the metallicity at the SN explosion site in M101 (∼0.5 Z). The spectroscopic rise-time (∼19 d) is significantly longer than that measured from the early optical light curve, implying a 'dark phase' of ∼1 d. A longer rise-time has significant implications when deducing the properties of the white dwarf and binary system from the early photometric behaviour.