skip to main content


Title: Relativistic Bondi accretion for stiff equations of state
ABSTRACT We revisit Bondi accretion – steady-state, adiabatic, spherical gas flow on to a Schwarzschild black hole at rest in an asymptotically homogeneous medium – for stiff polytropic equations of state (EOSs) with adiabatic indices Γ > 5/3. A general relativistic treatment is required to determine their accretion rates, for which we provide exact expressions. We discuss several qualitative differences between results for soft and stiff EOSs – including the appearance of a minimum steady-state accretion rate for EOSs with Γ ≥ 5/3 – and explore limiting cases in order to examine these differences. As an example, we highlight results for Γ = 2, which is often used in numerical simulations to model the EOS of neutron stars. We also discuss a special case with this index, the ultrarelativistic ‘causal’ EOS, P = ρ. The latter serves as a useful limit for the still undetermined neutron star EOS above nuclear density. The results are useful, for example, to estimate the accretion rate on to a mini-black hole residing at the centre of a neutron star.  more » « less
Award ID(s):
2010394 2006066 1707526 1662211
NSF-PAR ID:
10232903
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
502
Issue:
2
ISSN:
0035-8711
Page Range / eLocation ID:
3003 to 3011
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT

    GW190425 was the second gravitational wave (GW) signal compatible with a binary neutron star (BNS) merger detected by the Advanced LIGO and Advanced Virgo detectors. Since no electromagnetic counterpart was identified, whether the associated kilonova was too dim or the localization area too broad is still an open question. We simulate 28 BNS mergers with the chirp mass of GW190425 and mass ratio 1 ≤ q ≤ 1.67, using numerical-relativity simulations with finite-temperature, composition dependent equations of state (EOS) and neutrino radiation. The energy emitted in GWs is $\lesssim 0.083\mathrm{\, M_\odot }c^2$ with peak luminosity of 1.1–$2.4\times ~10^{58}/(1+q)^2\, {\rm {erg \, s^{-1}}}$. Dynamical ejecta and disc mass range between 5 × 10−6–10−3 and 10−5–$0.1 \mathrm{\, M_\odot }$, respectively. Asymmetric mergers, especially with stiff EOSs, unbind more matter and form heavier discs compared to equal mass binaries. The angular momentum of the disc is 8–$10\mathrm{\, M_\odot }~GM_{\rm {disc}}/c$ over three orders of magnitude in Mdisc. While the nucleosynthesis shows no peculiarity, the simulated kilonovae are relatively dim compared with GW170817. For distances compatible with GW190425, AB magnitudes are always dimmer than ∼20 mag for the B, r, and K bands, with brighter kilonovae associated to more asymmetric binaries and stiffer EOSs. We suggest that, even assuming a good coverage of GW190425’s sky location, the kilonova could hardly have been detected by present wide-field surveys and no firm constraints on the binary parameters or EOS can be argued from the lack of the detection.

     
    more » « less
  2. 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
  3. ABSTRACT

    The merging of a binary system involving two neutron stars (NSs), or a black hole (BH) and an NS, often results in the emission of an electromagnetic (EM) transient. One component of this EM transient is the epic explosion known as a kilonova (KN). The characteristics of the KN emission can be used to probe the equation of state (EoS) of NS matter responsible for its formation. We predict KN light curves from computationally simulated BH–NS mergers, by using the 3D radiative transfer code possis. We investigate two EoSs spanning most of the allowed range of the mass–radius diagram. We also consider a soft EoS compatible with the observational data within the so-called 2-families scenario in which hadronic stars co-exist with strange stars. Computed results show that the 2-families scenario, characterized by a soft EoS, should not produce a KN unless the mass of the binary components are small (MBH ≤ 6 M⊙ and MNS ≤ 1.4 M⊙) and the BH is rapidly spinning (χBH ≥ 0.3). In contrast, a strong KN signal potentially observable from future surveys (e.g. the Vera Rubin Observatory) is produced in the 1-family scenario for a wider region of the parameter space, and even for non-rotating BHs (χBH = 0) when MBH = 4 M⊙ and MNS = 1.2 M⊙. We also provide a fit that allows for the calculation of the unbound mass from the observed KN magnitude, without running timely and costly radiative transfer simulations. Findings presented in this paper will be used to interpret light curves anticipated during the fourth observing run (O4), of the advanced LIGO, advanced Virgo, and KAGRA interferometers and thus to constrain the EoS of NS matter.

     
    more » « less
  4. ABSTRACT

    Stellar-mass binary black holes (BBHs) embedded in active galactic nucleus (AGN) discs offer a promising dynamical channel to produce black hole mergers that are detectable by LIGO/Virgo. Modelling the interactions between the disc gas and the embedded BBHs is crucial to understand their orbital evolution. Using a suite of 2D high-resolution simulations of prograde equal-mass circular binaries in local disc models, we systematically study how their hydrodynamical evolution depends on the equation of state (EOS; including the γ-law and isothermal EOS) and on the binary mass and separation scales (relative to the supermassive black hole mass and the Hill radius, respectively). We find that binaries accrete slower and contract in orbit if the EOS is far from isothermal such that the surrounding gas is diffuse, hot, and turbulent. The typical orbital decay rate is of the order of a few times the mass doubling rate. For a fixed EOS, the accretion flows are denser, hotter, and more turbulent around more massive or tighter binaries. The torque associated with accretion is often comparable to the gravitational torque, so both torques are essential in determining the long-term binary orbital evolution. We carry out additional simulations with non-accreting binaries and find that their orbital evolution can be stochastic and is sensitive to the gravitational softening length, and the secular orbital evolution can be very different from those of accreting binaries. Our results indicate that stellar-mass BBHs may be hardened efficiently under ideal conditions, namely less massive and wider binaries embedded in discs with a non-isothermal EOS.

     
    more » « less
  5. ABSTRACT

    Parametric equations of state (EoSs) provide an important tool for systematically studying EoS effects in neutron star merger simulations. In this work, we perform a numerical validation of the M*-framework for parametrically calculating finite-temperature EoS tables. The framework, introduced by Raithel et al., provides a model for generically extending any cold, β-equilibrium EoS to finite temperatures and arbitrary electron fractions. In this work, we perform numerical evolutions of a binary neutron star merger with the SFHo finite-temperature EoS, as well as with the M*-approximation of this same EoS, where the approximation uses the zero-temperature, β-equilibrium slice of SFHo and replaces the finite-temperature and composition-dependent parts with the M*-model. We find that the approximate version of the EoS is able to accurately recreate the temperature and thermal pressure profiles of the binary neutron star remnant, when compared to the results found using the full version of SFHo. We additionally find that the merger dynamics and gravitational wave signals agree well between both cases, with differences of $\lesssim 1\!-\!2\,{\textrm{per cent}}$ introduced into the post-merger gravitational wave peak frequencies by the approximations of the EoS. We conclude that the M*-framework can be reliably used to probe neutron star merger properties in numerical simulations.

     
    more » « less