ABSTRACT
The electromagnetic observations of GW170817 were able to dramatically increase our understanding of neutron star mergers beyond what we learned from gravitational waves alone. These ...observations provided insight on all aspects of the merger from the nature of the gamma-ray burst to the characteristics of the ejected material. The ejecta of neutron star mergers are expected to produce such electromagnetic transients, called kilonovae or macronovae. Characteristics of the ejecta include large velocity gradients, relative to supernovae, and the presence of heavy r-process elements, which pose significant challenges to the accurate calculation of radiative opacities and radiation transport. For example, these opacities include a dense forest of bound–bound features arising from near-neutral lanthanide and actinide elements. Here we investigate the use of fine-structure, line-binned opacities that preserve the integral of the opacity over frequency. Advantages of this area-preserving approach over the traditional expansion–opacity formalism include the ability to pre-calculate opacity tables that are independent of the type of hydrodynamic expansion and thus eliminate the computational expense of calculating opacities within radiation-transport simulations. Tabular opacities are generated for all 14 lanthanides as well as a representative actinide element, uranium. We demonstrate that spectral simulations produced with the line-binned opacities agree well with results produced with the more accurate continuous Monte Carlo Sobolev approach, as well as with the commonly used expansion–opacity formalism. The agreement between the line-binned and expansion–opacity results is explained as arising from the similarity in their opacities in the limit of low optical depth, where radiation transport is important in the ejecta. Additional investigations illustrate the convergence of opacity with respect to the number of included lines, and elucidate sensitivities to different atomic physics approximations, such as fully and semirelativistic approaches.
Abstract
Depending upon the properties of their compact remnants and the physics included in the models, simulations of neutron star mergers can produce a broad range of ejecta properties. The ...characteristics of this ejecta, in turn, define the kilonova emission. To explore the effect of ejecta properties, we present a grid of two-component 2D axisymmetric kilonova simulations that vary mass, velocity, morphology, and composition. The masses and velocities of each component vary, respectively, from 0.001 to 0.1
M
⊙
and 0.05 to 0.3
c
, covering much of the range of results from the neutron star merger literature. The set of 900 models is constrained to have a toroidal low electron fraction (
Y
e
) ejecta with a robust
r
-process composition and either a spherical or lobed high-
Y
e
ejecta with two possible compositions. We simulate these models with the Monte Carlo radiative transfer code
SuperNu
using a full suite of lanthanide and fourth-row element opacities. We examine the trends of these models with parameter variation, show how they can be used with statistical tools, and compare the model light curves and spectra to those of AT2017gfo, the electromagnetic counterpart of GW170817.
We present the first three-dimensional smoothed particle hydrodynamics simulations of the induced gravitational collapse scenario of long-duration gamma-ray bursts (GRBs) associated with supernovae ...(SNe). We simulate the SN explosion of a carbon-oxygen core (COcore) forming a binary system with a neutron star (NS) companion. We follow the evolution of the SN ejecta, including their morphological structure, subject to the gravitational field of both the new NS ( NS) formed at the center of the SN and the one of the NS companion. We compute the accretion rate of the SN ejecta onto the NS companion, as well as onto the NS from SN matter fallback. We determine the fate of the binary system for a wide parameter space including different COcore and NS companion masses, orbital periods, and SN explosion geometry and energies. We identify, for selected NS nuclear equations of state, the binary parameters leading the NS companion, by hypercritical accretion, either to the mass-shedding limit or to the secular axisymmetric instability for gravitational collapse to a black hole (BH), or to a more massive, fast-rotating, stable NS. We also assess whether the binary remains gravitationally bound after the SN explosion, hence exploring the space of binary and SN explosion parameters leading to NS-NS and NS-BH binaries. The consequences of our results for the modeling of long GRBs, i.e., X-ray flashes and binary-driven hypernovae, are discussed.
Gamma-ray bursts (GRBs) are flashes of high-energy radiation arising from energetic cosmic explosions. Bursts of long (greater than two seconds) duration are produced by the core-collapse of massive ...stars
, and those of short (less than two seconds) duration by the merger of compact objects, such as two neutron stars
. A third class of events with hybrid high-energy properties was identified
, but never conclusively linked to a stellar progenitor. The lack of bright supernovae rules out typical core-collapse explosions
, but their distance scales prevent sensitive searches for direct signatures of a progenitor system. Only tentative evidence for a kilonova has been presented
. Here we report observations of the exceptionally bright GRB 211211A, which classify it as a hybrid event and constrain its distance scale to only 346 megaparsecs. Our measurements indicate that its lower-energy (from ultraviolet to near-infrared) counterpart is powered by a luminous (approximately 10
erg per second) kilonova possibly formed in the ejecta of a compact object merger.
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 Following the induced gravitational collapse (IGC) paradigm of gamma-ray bursts (GRBs) associated with type Ib/c supernovae, we present numerical simulations of the explosion of a ...carbon-oxygen (CO) core in a binary system with a neutron-star (NS) companion. The supernova ejecta trigger a hypercritical accretion process onto the NS thanks to a copious neutrino emission and the trapping of photons within the accretion flow. We show that temperatures of 1-10 MeV develop near the NS surface, hence electron-positron annihilation into neutrinos becomes the main cooling channel leading to accretion rates of 10−9- s−1 and neutrino luminosities of 1043-1052 erg s−1 (the shorter the orbital period the higher the accretion rate). We estimate the maximum orbital period, as a function of the NS initial mass, up to which the NS companion can reach by hypercritical accretion the critical mass for gravitational collapse leading to black hole formation. We then estimate the effects of the accreting and orbiting NS companion onto a novel geometry of the supernova ejecta density profile. We present the results of a particle simulation which show that the NS induces accentuated asymmetries in the ejecta density around the orbital plane. We elaborate on the observables associated with the above features of the IGC process. We apply this framework to specific GRBs: we find that X-ray flashes (XRFs) and binary-driven hypernovae are produced in binaries with and respectively. We analyze in detail the case of XRF 060218.
We constrain the properties of massive binaries by comparing radial velocity data on early-type stars in Cygnus OB2 with the expectations of Monte Carlo models. Our comparisons test several popular ...prescriptions for massive binary parameters. We explore a range of true binary fraction, F, a range of power-law slopes, a, describing the distribution of companion masses between the limits q sub(low) and 1, and a range of power-law slopes, beta , describing the distribution of orbital separations between the limits r sub(in) and r sub(out). We also consider distributions of secondary masses described by a Miller-Scalo type IMF and by a two-component IMF that Includes a substantial "twin" population with M sub(2) unk M sub(1). Several seemingly disparate prescriptions for massive binary characteristics can be reconciled by adopting carefully chosen values for F, r sub(in), and r sub(out). We show that binary fractions F < 0.7 are less probable than F greater than or equal to 0.8 for reasonable choices of r sub(in) and r sub(out). Thus, the true binary fraction is high. For F = 1.0 and a distribution of orbital separations near the canonical "Opik's law distribution (i.e., flat; beta = 0), the power-law slope of the mass ratio distribution is alpha = -0.6 to 0.0. For F unk 0.8, alpha is somewhat larger, in the range -0.4 to 1.0. In any case, the secondary star mass function is inconsistent with a Miller-Scalo-like IMF unless the lower end is truncated below similar to 2-4 M unk. In other words, massive stars preferentially have massive companions. The best-fitting models are described by a Salpeter or Miller-Scalo IMF for 60% of secondary star masses with the other similar to 40% of secondaries having M sub(2) unk M sub(1), i.e., "twins." These model parameters simultaneously predict the fraction of Type Ib/c supernovae to be 30%-40% of all core-collapse supernovae, in agreement with recent observational estimates.
We perform a suite of smoothed particle hydrodynamics simulations to investigate in detail the results of a giant impact on the young Uranus. We study the internal structure, rotation rate, and ...atmospheric retention of the post-impact planet, as well as the composition of material ejected into orbit. Most of the material from the impactor's rocky core falls in to the core of the target. However, for higher angular momentum impacts, significant amounts become embedded anisotropically as lumps in the ice layer. Furthermore, most of the impactor's ice and energy is deposited in a hot, high-entropy shell at a radius of ∼3 R⊕. This could explain Uranus' observed lack of heat flow from the interior and be relevant for understanding its asymmetric magnetic field. We verify the results from the single previous study of lower resolution simulations that an impactor with a mass of at least 2 M⊕ can produce sufficiently rapid rotation in the post-impact Uranus for a range of angular momenta. At least 90% of the atmosphere remains bound to the final planet after the collision, but over half can be ejected beyond the Roche radius by a 2 or 3 M⊕ impactor. This atmospheric erosion peaks for intermediate impactor angular momenta (∼3 × 1036 kg m2 s−1). Rock is more efficiently placed into orbit and made available for satellite formation by 2 M⊕ impactors than 3 M⊕ ones, because it requires tidal disruption that is suppressed by the more massive impactors.
Context.
As a result of their formation via massive single and binary stellar evolution, the masses of stellar-remnant black holes (BH) are subjects of great interest in this era of ...gravitational-wave detection from binary black hole (BBH) and binary neutron star merger events.
Aims.
In this work, we present new developments in the stellar-remnant formation and related schemes of the current
N
-body evolution program
NBODY7
. We demonstrate that the newly implemented stellar-wind and remnant-formation schemes in the stellar-evolutionary sector or
BSE
of the
NBODY7
code, such as the “rapid” and the “delayed” supernova (SN) schemes along with an implementation of pulsational-pair-instability and pair-instability supernova (PPSN/PSN), now produce neutron star (NS) and BH masses that agree nearly perfectly, over large ranges of zero-age-main-sequence (ZAMS) mass and metallicity, with those from the widely recognised
StarTrack
population-synthesis program. We also demonstrate the new, recipe-based implementations of various widely debated mechanisms of natal kicks on NSs and BHs, such as “convection-asymmetry-driven”, “collapse-asymmetry-driven”, and “neutrino-emission-driven” kicks, in addition to a fully consistent implementation of the standard, fallback-dependent, momentum-conserving natal kick.
Methods.
All the above newly implemented schemes are also shared with the standalone versions of
SSE
and
BSE
. All these demonstrations are performed with both the updated standalone
BSE
and the updated
NBODY7
/
BSE
.
Results.
When convolved with stellar and primordial-binary populations as observed in young massive clusters, such remnant-formation and natal-kick mechanisms crucially determine the accumulated number, mass, and mass distribution of the BHs retained in young massive, open, and globular clusters (GCs); these BHs would eventually become available for long-term dynamical processing.
Conclusions.
Among other conclusions, we find that although the newer, delayed SN remnant formation model gives birth to the largest number (mass) of BHs, the older remnant-formation schemes cause the largest number (mass) of BHs to survive in clusters, when incorporating SN material fallback onto the BHs. The SN material fallback also causes the convection-asymmetry-driven SN kick to effectively retain similar numbers and masses of BHs in clusters as for the standard, momentum-conserving kick. The collapse-asymmetry-driven SN kick would cause nearly all BHs to be retained in clusters irrespective of their mass, remnant-formation model, and metallicity, whereas the inference of a large population of BHs in GCs would potentially rule out the neutrino-driven SN kick mechanism. Pre-SN mergers of massive primordial binaries would potentially cause BH masses to deviate from the theoretical, single-star ZAMS to mass-remnant mass relation unless a substantial of the total merging stellar mass of up to ≈40% is lost during a merger process. In particular, such mergers, at low metallicities, have the potential to produce low-spinning BHs within the PSN mass gap that can be retained in a stellar cluster and be available for subsequent dynamical interactions. As recent studies indicate, the new remnant-formation modelling reassures us that young massive and open clusters would potentially contribute to the dynamical BBH merger detection rate to a similar extent as their more massive GC counterparts.
PDF
