skip to main content


Title: Spherically symmetric accretion on to a compact object through a standing shock: the effects of general relativity in the Schwarzschild geometry
ABSTRACT

A core-collapse supernova is generated by the passage of a shock wave through the envelope of a massive star, where the shock wave is initially launched from the ‘bounce’ of the neutron star formed during the collapse of the stellar core. Instead of successfully exploding the star, however, numerical investigations of core-collapse supernovae find that this shock tends to ‘stall’ at small radii (≲10 neutron star radii), with stellar material accreting on to the central object through the standing shock. Here, we present time-steady, adiabatic solutions for the density, pressure, and velocity of the shocked fluid that accretes on to the compact object through the stalled shock, and we include the effects of general relativity in the Schwarzschild metric. Similar to previous works that were carried out in the Newtonian limit, we find that the gas ‘settles’ interior to the stalled shock; in the relativistic regime analysed here, the velocity asymptotically approaches zero near the Schwarzschild radius. These solutions can represent accretion on to a material surface if the radius of the compact object is outside of its event horizon, such as a neutron star; we also discuss the possibility that these solutions can approximately represent the accretion of gas on to a newly formed black hole following a core-collapse event. Our findings and solutions are particularly relevant in weak and failed supernovae, where the shock is pushed to small radii and relativistic effects are large.

 
more » « less
Award ID(s):
2006684
NSF-PAR ID:
10372055
Author(s) / Creator(s):
;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
516
Issue:
4
ISSN:
0035-8711
Page Range / eLocation ID:
p. 4814-4821
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Observations and theory suggest that core-collapse supernovae can span a range of explosion energies, and when sub-energetic the shockwave initiating the explosion can decelerate to speeds comparable to the escape speed of the progenitor. In these cases, gravity will complicate the explosion hydrodynamics and conceivably cause the shock to stall at large radii within the progenitor star. To understand these unique properties of weak explosions, we develop a perturbative approach for modeling the propagation of an initially strong shock into a time-steady, infalling medium in the gravitational field of a compact object. This method writes the shock position and the post-shock velocity, density, and pressure as series solutions in the (time-dependent) ratio of the freefall speed to the shock speed, and predicts that the shock stalls within the progenitor if the explosion energy is below a critical value. We show that our model agrees very well with hydrodynamic simulations, and accurately predicts (for example) the time-dependent shock position and velocity and the radius at which the shock stalls. Our results have implications for black hole formation and the newly detected class of fast X-ray transients (FXTs). In particular, we propose that a “phantom shock breakout”—where the outer edge of the star falls through a stalled shock—can yield a burst of X-rays without a subsequent optical/UV signature, similar to FXTs. This model predicts the rise time of the X-ray burst,td, and the mean photon energy,kT, are anticorrelated, approximately asTtd5/8.

     
    more » « less
  2. Abstract

    Some massive stars likely fail to produce core-collapse supernovae, but these failed supernovae (FSNe) can generate an electromagnetic outburst prior to the disappearance of the star, as the mass lost to neutrinos during the stellar core collapse results in the formation and breakout of a second shock. We show that, when the mass lost to neutrinos is sufficiently small, there are two self-similar solutions that describe the propagation of a weak shock into a hydrodynamically expanding envelope that simultaneously yield accretion onto the black hole. The larger Mach number solution is unstable and yields the minimum Mach number that a shock must have to strengthen into the energy-conserving regime. Above a critical mass loss, there are no weak-shock solutions, implying that there are only strong explosions if the neutrino mass loss is above a critical value, and this value is a few percent of the mass of the star (and is physically achievable) for typical parameters. Our results imply that the fate of the explosion from an FSN—weak with little to no mass ejection or strong with the expulsion of the majority of the envelope—is a sensitive function of the stellar properties and the neutrino mass loss. We also show that there is a second type of self-similar solution for the shock that results in thesettlingof the gas near the compact object, which may be applicable to nonterminal stellar eruptions and the response of a gaseous disk to gravitational-wave induced mass loss from a binary black hole merger.

     
    more » « less
  3. Abstract

    For the first ∼3 yrs after the binary neutron star merger event GW 170817, the radio and X-ray radiation has been dominated by emission from a structured relativistic off-axis jet propagating into a low-density medium withn< 0.01 cm−3. We report on observational evidence for an excess of X-ray emission atδt> 900 days after the merger. WithLx≈ 5 × 1038erg s−1at 1234 days, the recently detected X-ray emission represents a ≥3.2σ(Gaussian equivalent) deviation from the universal post-jet-break model that best fits the multiwavelength afterglow at earlier times. In the context ofJetFitafterglow models, current data represent a departure with statistical significance ≥3.1σ, depending on the fireball collimation, with the most realistic models showing excesses at the level of ≥3.7σ. A lack of detectable 3 GHz radio emission suggests a harder broadband spectrum than the jet afterglow. These properties are consistent with the emergence of a new emission component such as synchrotron radiation from a mildly relativistic shock generated by the expanding merger ejecta, i.e., a kilonova afterglow. In this context, we present a set of ab initio numerical relativity binary neutron star (BNS) merger simulations that show that an X-ray excess supports the presence of a high-velocity tail in the merger ejecta, and argues against the prompt collapse of the merger remnant into a black hole. Radiation from accretion processes on the compact-object remnant represents a viable alternative. Neither a kilonova afterglow nor accretion-powered emission have been observed before, as detections of BNS mergers at this phase of evolution are unprecedented.

     
    more » « less
  4. Abstract

    We perform the first 3D ab-initio general-relativistic neutrino-radiation hydrodynamics of a long-lived neutron star merger remnant spanning a fraction of its cooling timescale. We find that neutrino cooling becomes the dominant energy loss mechanism after the gravitational-wave dominated phase (∼20 ms postmerger). Electron flavor antineutrino luminosity dominates over electron flavor neutrino luminosity at early times, resulting in a secular increase of the electron fraction in the outer layers of the remnant. However, the two luminosities become comparable ∼20–40 ms postmerger. A dense gas of electron antineutrinos is formed in the outer core of the remnant at densities ∼1014.5g cm−3, corresponding to temperature hot spots. The neutrinos account for ∼10% of the lepton number in this region. Despite the negative radial temperature gradient, the radial entropy gradient remains positive, and the remnant is stably stratified according to the Ledoux criterion for convection. A massive accretion disk is formed from the material squeezed out of the collisional interface between the stars. The disk carries a large fraction of the angular momentum of the system, allowing the remnant massive neutron star to settle to a quasi-steady equilibrium within the region of possible, stable, rigidly rotating configurations. The remnant is differentially rotating, but it is stable against the magnetorotational instability. Other MHD mechanisms operating on longer timescales are likely responsible for the removal of the differential rotation. Our results indicate the remnant massive neutron star is thus qualitatively different from a protoneutron stars formed in core-collapse supernovae.

     
    more » « less
  5. Abstract

    Progenitor models for the “luminous” subclass of Fast Blue Optical Transients (LFBOTs; prototype: AT2018cow) are challenged to simultaneously explain all of their observed properties: fast optical rise times of days or less; peak luminosities ≳1044erg s−1; low yields ≲0.1Mof56Ni; aspherical ejecta with a wide velocity range (≲3000 km s−1to ≳0.1–0.5cwith increasing polar latitude); presence of hydrogen-depleted-but-not-free dense circumstellar material (CSM) on radial scales from ∼1014cm to ∼3 × 1016cm; embedded variable source of non-thermal X-ray/γ-rays, suggestive of a compact object. We show that all of these properties are consistent with the tidal disruption and hyper-accretion of a Wolf-Rayet (WR) star by a black hole or neutron star binary companion. In contrast with related previous models, the merger occurs with a long delay (≳100 yr) following the common envelope (CE) event responsible for birthing the binary, as a result of gradual angular momentum loss to a relic circumbinary disk. Disk-wind outflows from the merger-generated accretion flow generate the56Ni-poor aspherical ejecta with the requisite velocity range. The optical light curve is powered primarily by reprocessing X-rays from the inner accretion flow/jet, though CSM shock interaction also contributes. Primary CSM sources include WR mass loss from the earliest stages of the merger (≲1014cm) and the relic CE disk and its photoevaporation-driven wind (≳1016cm). Longer delayed mergers may instead give rise to supernovae Type Ibn/Icn (depending on the WR evolutionary state), connecting these transient classes with LFBOTs.

     
    more » « less