-
Estimating outflow masses and velocities in merger simulations: impact of r-process heating and neutrino cooling
Authors:
Francois Foucart,
Philipp Moesta,
Teresita Ramirez,
Alex James Wright,
Siva Darbha,
Daniel Kasen
Abstract:
The determination of the mass, composition, and geometry of matter outflows in black hole-neutron star and neutron star-neutron star binaries is crucial to current efforts to model kilonovae, and to understand the role of neutron star merger in r-process nucleosynthesis. In this manuscript, we review the simple criteria currently used in merger simulations to determine whether matter is unbound an…
▽ More
The determination of the mass, composition, and geometry of matter outflows in black hole-neutron star and neutron star-neutron star binaries is crucial to current efforts to model kilonovae, and to understand the role of neutron star merger in r-process nucleosynthesis. In this manuscript, we review the simple criteria currently used in merger simulations to determine whether matter is unbound and what the asymptotic velocity of ejected material will be. We then show that properly accounting for both heating and cooling during r-process nucleosynthesis is important to accurately predict the mass and kinetic energy of the outflows. These processes are also likely to be crucial to predict the fallback timescale of any bound ejecta. We derive a model for the asymptotic veloicity of unbound matter and binding energy of bound matter that accounts for both of these effects and that can easily be implemented in merger simulations. We show, however, that the detailed velocity distribution and geometry of the outflows can currently only be captured by full 3D fluid simulations of the outflows, as non-local effect ignored by the simple criteria used in merger simulations cannot be safely neglected when modeling these effects. Finally, we propose the introduction of simple source terms in the fluid equations to approximately account for heating/cooling from r-process nucleosynthesis in future seconds-long 3D simulations of merger remnants, without the explicit inclusion of out-of-nuclear statistical equilibrium reactions in the simulations.
△ Less
Submitted 17 November, 2021; v1 submitted 1 September, 2021;
originally announced September 2021.
-
Electromagnetic Signatures from the Tidal Tail of a Black Hole - Neutron Star Merger
Authors:
Siva Darbha,
Daniel Kasen,
Francois Foucart,
Daniel J. Price
Abstract:
Black hole - neutron star (BH-NS) mergers are a major target for ground-based gravitational wave (GW) observatories. A merger can also produce an electromagnetic counterpart (a kilonova) if it ejects neutron-rich matter that assembles into heavy elements through r-process nucleosynthesis. We study the kilonova signatures of the unbound dynamical ejecta of a BH-NS merger. We take as our initial sta…
▽ More
Black hole - neutron star (BH-NS) mergers are a major target for ground-based gravitational wave (GW) observatories. A merger can also produce an electromagnetic counterpart (a kilonova) if it ejects neutron-rich matter that assembles into heavy elements through r-process nucleosynthesis. We study the kilonova signatures of the unbound dynamical ejecta of a BH-NS merger. We take as our initial state the results from a numerical relativity simulation, and then use a general relativistic hydrodynamics code to study the evolution of the ejecta with parameterized r-process heating models. The unbound dynamical ejecta is initially a flattened, directed tidal tail largely confined to a plane. Heating from the r-process inflates the ejecta into a more spherical shape and smooths its small-scale structure, though the ejecta retains its bulk directed motion. We calculate the electromagnetic signatures using a 3D radiative transfer code and a parameterized opacity model for lanthanide-rich matter. The light curve varies with viewing angle due to two effects: asphericity results in brighter emission for orientations with larger projected areas, while Doppler boosting results in brighter emission for viewing angles more aligned with the direction of bulk motion. For typical r-process heating rates, the peak bolometric luminosity varies by a factor of $\sim 3$ with orientation while the peak in the optical bands varies by $\sim 3$ magnitudes. The spectrum is blue-shifted at viewing angles along the bulk motion, which increases the $V$-band peak magnitude to $\sim -14$ despite the lanthanide-rich composition.
△ Less
Submitted 24 August, 2021; v1 submitted 4 March, 2021;
originally announced March 2021.
-
Inclination Dependence of Kilonova Light Curves from Globally Aspherical Geometries
Authors:
Siva Darbha,
Daniel Kasen
Abstract:
The merger of two neutron stars (NSs) or a neutron star and a black hole (BH) produces a radioactively-powered transient known as a kilonova, first observed accompanying the gravitational wave event GW170817. While kilonovae are frequently modeled in spherical symmetry, the dynamical ejecta and disk outflows can be considerably asymmetric. We use Monte Carlo radiative transfer calculations to stud…
▽ More
The merger of two neutron stars (NSs) or a neutron star and a black hole (BH) produces a radioactively-powered transient known as a kilonova, first observed accompanying the gravitational wave event GW170817. While kilonovae are frequently modeled in spherical symmetry, the dynamical ejecta and disk outflows can be considerably asymmetric. We use Monte Carlo radiative transfer calculations to study the light curves of kilonovae with globally axisymmetric geometries (e.g. an ellipsoid and a torus). We find that the variation in luminosity in these models is most pronounced at early times, and decreases until the light curves become isotropic in the late optically thin phase. The light curve shape and peak time are not significantly modified by the global asymmetry. We show that the projected surface area along the line of sight captures the primary geometric effects, and use this fact to provide a simple analytic estimate of the direction-dependent light curves of the aspherical ejecta. For the kilonova accompanying GW170817, accounting for asymmetry with an oblate (prolate) ellipsoid of axial ratio $2$ ($1/2$) leads to a $\sim 40 \%$ decrease (increase) in the inferred ejecta mass compared to the spherical case. The pole-to-equator orientation effects are expected to be significantly larger (a factor of $\sim 5 - 10$) for the more extreme asymmetries expected for some NS-BH mergers.
△ Less
Submitted 15 March, 2021; v1 submitted 1 February, 2020;
originally announced February 2020.
-
Ultra-deep tidal disruption events: prompt self-intersections and observables
Authors:
Siva Darbha,
Eric R. Coughlin,
Daniel Kasen,
Chris Nixon
Abstract:
A star approaching a supermassive black hole (SMBH) can be torn apart in a tidal disruption event (TDE). We examine ultra-deep TDEs, a new regime in which the disrupted debris approaches close to the black hole's Schwarzschild radius, and the leading part intersects the trailing part at the first pericenter passage. We calculate the range of penetration factors $β$ vs SMBH masses $M$ that produce…
▽ More
A star approaching a supermassive black hole (SMBH) can be torn apart in a tidal disruption event (TDE). We examine ultra-deep TDEs, a new regime in which the disrupted debris approaches close to the black hole's Schwarzschild radius, and the leading part intersects the trailing part at the first pericenter passage. We calculate the range of penetration factors $β$ vs SMBH masses $M$ that produce these prompt self-intersections using a Newtonian analytic estimate and a general relativistic (GR) geodesic model. We find that significant self-intersection of Solar-type stars requires $β\sim 50 - 127$ for $M/M_\odot = 10^4$, down to $β\sim 5.6 - 5.9$ for $M/M_\odot = 10^6$. We run smoothed-particle hydrodynamic (SPH) simulations to corroborate our calculations and find close agreement, with a slightly shallower dependence on $M$. We predict that the shock from the collision emits an X-ray flare lasting $t \sim 2$ s with $L \sim 10^{47}$ ergs/s at $E \sim 2$ keV, and the debris has a prompt accretion episode lasting $t \sim$ several min. The events are rare and occur with a rate $\dot{N} \lesssim 10^{-7}$ Mpc$^{-3}$ yr$^{-1}$. Ultra-deep TDEs can probe the strong gravity and demographics of low-mass SMBHs.
△ Less
Submitted 16 May, 2019;
originally announced May 2019.
-
Gravitational interactions of stars with supermassive black hole binaries - II. Hypervelocity stars
Authors:
Siva Darbha,
Eric R. Coughlin,
Daniel Kasen,
Eliot Quataert
Abstract:
Supermassive black holes (SMBHs) in galactic nuclei can eject hypervelocity stars (HVSs). Using restricted three-body integrations, we study the properties of stars ejected by circular, binary SMBHs as a function of their mass ratios $q = M_2 / M_1$ and separations $a$, focusing on the stellar velocity and angular distributions. We find that the ejection probability is an increasing function of…
▽ More
Supermassive black holes (SMBHs) in galactic nuclei can eject hypervelocity stars (HVSs). Using restricted three-body integrations, we study the properties of stars ejected by circular, binary SMBHs as a function of their mass ratios $q = M_2 / M_1$ and separations $a$, focusing on the stellar velocity and angular distributions. We find that the ejection probability is an increasing function of $q$ and $a$, and that the mean ejected velocity scales with $q$ and $a$ similar to previous work but with modified scaling constants. Binary SMBHs tend to eject their fastest stars toward the binary orbital plane. We calculate the ejection rates as the binary SMBHs inspiral, and find that they eject stars with velocities $v_\infty > 1000$ km/s at rates of $\sim 4 \times 10^{-2} - 2 \times 10^{-1}$ yr$^{-1}$ for $q = 1$ ($\sim 10^{-4} - 10^{-3}$ yr$^{-1}$ for $q = 0.01$) over their lifetimes, and can emit a burst of HVSs with $v_\infty > 3000$ km/s as they coalesce. We integrate the stellar distributions over the binary SMBH inspiral and compare them to those produced by the "Hills mechanism" (in which a single SMBH ejects a star after tidally separating a binary star system), and find that $N \sim 100$ HVS velocity samples with $v_\infty \gtrsim 200$ km/s are needed to confidently distinguish between a binary and single SMBH origin.
△ Less
Submitted 11 October, 2018;
originally announced October 2018.
-
Gravitational interactions of stars with supermassive black hole binaries. I. Tidal disruption events
Authors:
Siva Darbha,
Eric R. Coughlin,
Daniel Kasen,
Eliot Quataert
Abstract:
Stars approaching supermassive black holes (SMBHs) in the centers of galaxies can be torn apart by strong tidal forces. We study the physics of tidal disruption by a binary SMBH as a function of the binary mass ratio $q = M_2 / M_1$ and separation $a$, exploring a large set of points in the parameter range $q \in [0.01, 1]$ and $a/r_{t1} \in [10, 1000]$. We simulate encounters in which field stars…
▽ More
Stars approaching supermassive black holes (SMBHs) in the centers of galaxies can be torn apart by strong tidal forces. We study the physics of tidal disruption by a binary SMBH as a function of the binary mass ratio $q = M_2 / M_1$ and separation $a$, exploring a large set of points in the parameter range $q \in [0.01, 1]$ and $a/r_{t1} \in [10, 1000]$. We simulate encounters in which field stars approach the binary from the loss cone on parabolic, low angular momentum orbits. We present the rate of disruption and the orbital properties of the disrupted stars, and examine the fallback dynamics of the post-disruption debris in the "frozen-in" approximation. We conclude by calculating the time-dependent disruption rate over the lifetime of the binary. Throughout, we use a primary mass $M_1 = 10^6 M_\odot$ as our central example. We find that the tidal disruption rate is a factor of $\sim 2 - 7$ times larger than the rate for an isolated BH, and is independent of $q$ for $q \gtrsim 0.2$. In the "frozen-in" model, disruptions from close, nearly equal mass binaries can produce intense tidal fallbacks: for binaries with $q \gtrsim 0.2$ and $a/r_{t1} \sim 100$, roughly $\sim 18 - 40 \%$ of disruptions will have short rise times ($t_\textrm{rise} \sim 1 - 10$ d) and highly super-Eddington peak return rates ($\dot{M}_\textrm{peak} / \dot{M}_\textrm{Edd} \sim 2 \times 10^2 - 3 \times 10^3$).
△ Less
Submitted 21 February, 2018;
originally announced February 2018.
-
Stellar binaries incident on supermassive black hole binaries: implications for double tidal disruption events, calcium-rich transients, and hypervelocity stars
Authors:
Eric R. Coughlin,
Siva Darbha,
Daniel Kasen,
Eliot Quataert
Abstract:
We analyze the outcome of the interaction between a stellar binary and a supermassive black hole binary (SMBHB) by performing a large number of gravitational scattering experiments. Most of the encounters result in either the ejection of an intact binary or the ejection of two individual stars following the tidal breakup of the binary. However, tidal disruption events (TDEs) and mergers constitute…
▽ More
We analyze the outcome of the interaction between a stellar binary and a supermassive black hole binary (SMBHB) by performing a large number of gravitational scattering experiments. Most of the encounters result in either the ejection of an intact binary or the ejection of two individual stars following the tidal breakup of the binary. However, tidal disruption events (TDEs) and mergers constitute a few percent of the outcomes, and double, temporally distinct TDEs (i.e., separated by at least one orbit of the supermassive black hole binary) occur at the percent level. We also demonstrate that the properties of the ejected binaries are significantly altered through the interaction with the SMBHB, and their large eccentricities increase the merger rate and could lead to gravitational-wave inspirals far from the nucleus of the host galaxy. We discuss our results in the context of observed tidal disruption events, hypervelocity stars, and remote supernovae, such as calcium-rich transients.
△ Less
Submitted 8 August, 2018; v1 submitted 20 February, 2018;
originally announced February 2018.
-
Nickel-Rich Outflows Produced by the Accretion-Induced Collapse of White Dwarfs: Lightcurves and Spectra
Authors:
S. Darbha,
B. D. Metzger,
E. Quataert,
D. Kasen,
P. Nugent,
R. Thomas
Abstract:
The accretion-induced collapse (AIC) of a white dwarf to form a neutron star can leave behind a rotationally supported disk with mass of up to ~ 0.1 M_sun. The disk is initially composed of free nucleons but as it accretes and spreads to larger radii, the free nucleons recombine to form helium, releasing sufficient energy to unbind the remaining disk. Most of the ejected mass fuses to form Ni56 an…
▽ More
The accretion-induced collapse (AIC) of a white dwarf to form a neutron star can leave behind a rotationally supported disk with mass of up to ~ 0.1 M_sun. The disk is initially composed of free nucleons but as it accretes and spreads to larger radii, the free nucleons recombine to form helium, releasing sufficient energy to unbind the remaining disk. Most of the ejected mass fuses to form Ni56 and other iron group elements. We present spherically symmetric radiative transfer calculations of the transient powered by the radioactive heating of this ejecta. For an ejecta mass of 1e-2 M_sun (3e-3 M_sun), the lightcurve peaks after <~ 1 day with a peak bolometric luminosity ~ 2e41 erg/s (~ 5e40 erg/s), i.e., a "kilonova"; the decay time is ~ 4 (2) days. Overall, the spectra redden with time reaching U-V ~ 4 after ~ 1 day; the optical colors (B-V) are, however, somewhat blue. Near the peak in the lightcurve, the spectra are dominated by Doppler broadened Nickel features, with no distinct spectral lines present. At ~ 3-5 days, strong Calcium lines are present in the infrared, although the Calcium mass fraction is only ~ 1e-4.5. If rotationally supported disks are a common byproduct of AIC, current and upcoming transient surveys such as the Palomar Transient Factory should detect a few AIC per year for an AIC rate of ~ 1e-2 of the Type Ia rate. We discuss ways of distinguishing AIC from other rapid, faint transients, including .Ia's and the ejecta from binary neutron star mergers.
△ Less
Submitted 30 May, 2011; v1 submitted 6 May, 2010;
originally announced May 2010.
-
Electromagnetic Counterparts of Compact Object Mergers Powered by the Radioactive Decay of R-process Nuclei
Authors:
B. D. Metzger,
G. Martinez-Pinedo,
S. Darbha,
E. Quataert,
A. Arcones,
D. Kasen,
R. Thomas,
P. Nugent,
I. V. Panov,
N. T. Zinner
Abstract:
The most promising astrophysical sources of kHz gravitational waves (GWs) are the inspiral and merger of binary neutron star(NS)/black hole systems. Maximizing the scientific return of a GW detection will require identifying a coincident electro-magnetic (EM) counterpart. One of the most likely sources of isotropic EM emission from compact object mergers is a supernova-like transient powered by…
▽ More
The most promising astrophysical sources of kHz gravitational waves (GWs) are the inspiral and merger of binary neutron star(NS)/black hole systems. Maximizing the scientific return of a GW detection will require identifying a coincident electro-magnetic (EM) counterpart. One of the most likely sources of isotropic EM emission from compact object mergers is a supernova-like transient powered by the radioactive decay of heavy elements synthesized in ejecta from the merger. We present the first calculations of the optical transients from compact object mergers that self-consistently determine the radioactive heating by means of a nuclear reaction network; using this heating rate, we model the light curve with a one dimensional Monte Carlo radiation transfer calculation. For an ejecta mass ~1e-2 M_sun[1e-3 M_sun] the resulting light curve peaks on a timescale ~ 1 day at a V-band luminosity nu L_nu ~ 3e41[1e41] ergs/s (M_V = -15[-14]); this corresponds to an effective "f" parameter ~3e-6 in the Li-Paczynski toy model. We argue that these results are relatively insensitive to uncertainties in the relevant nuclear physics and to the precise early-time dynamics and ejecta composition. Due to the rapid evolution and low luminosity of NS merger transients, EM counterpart searches triggered by GW detections will require close collaboration between the GW and astronomical communities. NS merger transients may also be detectable following a short-duration Gamma-Ray Burst or "blindly" with present or upcoming optical transient surveys. Because the emission produced by NS merger ejecta is powered by the formation of rare r-process elements, current optical transient surveys can directly constrain the unknown origin of the heaviest elements in the Universe.
△ Less
Submitted 23 April, 2010; v1 submitted 28 January, 2010;
originally announced January 2010.