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

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 more » 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.

« less
Publication Date:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range or eLocation-ID:
p. 4814-4821
Oxford University Press
Sponsoring Org:
National Science Foundation
More Like this
  1. 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 mergermore »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.

    « less
  2. Abstract The proposal that core collapse supernovae are neutrino driven is still the subject of active investigation more than 50 years after the seminal paper by Colgate and White. The modern version of this paradigm, which we owe to Wilson, proposes that the supernova shock wave is powered by neutrino heating, mediated by the absorption of electron-flavor neutrinos and antineutrinos emanating from the proto-neutron star surface, or neutrinosphere. Neutrino weak interactions with the stellar core fluid, the theory of which is still evolving, are flavor and energy dependent. The associated neutrino mean free paths extend over many orders of magnitude and are never always small relative to the stellar core radius. Thus, neutrinos are never always fluid like. Instead, a kinetic description of them in terms of distribution functions that determine the number density of neutrinos in the six-dimensional phase space of position, direction, and energy, for both neutrinos and antineutrinos of each flavor, or in terms of angular moments of these neutrino distributions that instead provide neutrino number densities in the four-dimensional phase-space subspace of position and energy, is needed. In turn, the computational challenge is twofold: (i) to map the kinetic equations governing the evolution of these distributionsmore »or moments onto discrete representations that are stable, accurate, and, perhaps most important, respect physical laws such as conservation of lepton number and energy and the Fermi–Dirac nature of neutrinos and (ii) to develop efficient, supercomputer-architecture-aware solution methods for the resultant nonlinear algebraic equations. In this review, we present the current state of the art in attempts to meet this challenge.« less
  3. Abstract

    Long-duration gamma-ray bursts (lGRBs) originate in relativistic collimated outflows—jets—that drill their way out of collapsing massive stars. Accurately modeling this process requires realistic stellar profiles for the jets to propagate through and break out of. Most previous studies have used simple power laws or pre-collapse models for massive stars. However, the relevant stellar profile for lGRB models is in fact that of a star after its core has collapsed to form a compact object. To self-consistently compute such a stellar profile, we use the open-source code GR1D to simulate the core-collapse process for a suite of low-metallicity rotating massive stellar progenitors that have undergone chemically homogeneous evolution. Our models span a range of zero-age main-sequence (ZAMS) masses:MZAMS= 13, 18, 21, 25, 35, 40, and 45M. All of these models, at the onset of core-collapse, feature steep density profiles,ρrα, withα≈ 2.5, which would result in jets that are inconsistent with lGRB observables. We follow the collapses of four of the seven models until they form black holes (BHs) and the other three models until they form proto-neutron stars (PNSs). We find, across all models, that the density profile outside the newly formed BH or PNS is well represented by amore »flatter power law withα≈ 1.35–1.55. Such flat density profiles are conducive to the successful formation and breakout of BH-powered jets and are, in fact, required to reproduce observable properties of lGRBs. Future models of lGRBs should be initialized with shallower post-collapse stellar profiles, like those presented here, instead of the much steeper pre-collapse profiles that are typically used.

    « less
  4. ABSTRACT The core collapse of massive, rapidly-rotating stars are thought to be the progenitors of long-duration gamma-ray bursts (GRB) and their associated hyperenergetic supernovae (SNe). At early times after the collapse, relatively low angular momentum material from the infalling stellar envelope will circularize into an accretion disc located just outside the black hole horizon, resulting in high accretion rates necessary to power a GRB jet. Temperatures in the disc mid-plane at these small radii are sufficiently high to dissociate nuclei, while outflows from the disc can be neutron-rich and may synthesize r-process nuclei. However, at later times, and for high progenitor angular momentum, the outer layers of the stellar envelope can circularize at larger radii ≳ 107 cm, where nuclear reactions can take place in the disc mid-plane (e.g. 4He + 16O → 20Ne + γ). Here we explore the effects of nuclear burning on collapsar accretion discs and their outflows by means of hydrodynamical α-viscosity torus simulations coupled to a 19-isotope nuclear reaction network, which are designed to mimic the late infall epochs in collapsar evolution when the viscous time of the torus has become comparable to the envelope fall-back time. Our results address several key questions, such as the conditions for quiescent burningmore »and accretion versus detonation and the generation of 56Ni in disc outflows, which we show could contribute significantly to powering GRB SNe. Being located in the slowest, innermost layers of the ejecta, the latter could provide the radioactive heating source necessary to make the spectral signatures of r-process elements visible in late-time GRB-SNe spectra.« less
  5. Abstract

    Recent studies have highlighted the sensitivity of core-collapse supernovae (CCSNe) models to electron-capture (EC) rates on neutron-rich nuclei near theN= 50 closed-shell region. In this work, we perform a large suite of one-dimensional CCSN simulations for 200 stellar progenitors using recently updated EC rates in this region. For comparison, we repeat the simulations using two previous implementations of EC rates: a microphysical library with parametrizedN= 50 rates (LMP), and an older independent-particle approximation (IPA). We follow the simulations through shock revival up to several seconds post-bounce, and show that the EC rates produce a consistent imprint on CCSN properties, often surpassing the role of the progenitor itself. Notable impacts include the timescale of core collapse, the electron fraction and mass of the inner core at bounce, the accretion rate through the shock, the success or failure of revival, and the properties of the central compact remnant. We also compare the observable neutrino signal of the neutronization burst in a DUNE-like detector, and find consistent impacts on the counts and mean energies. Overall, the updated rates result in properties that are intermediate between LMP and IPA, and yet slightly more favorable to explosion than both.