Formation of Double Neutron Star Systems Tauris, T. M.; Kramer, M.; Freire, P. C. C. ...
Astrophysical journal/The Astrophysical journal,
09/2017, Letnik:
846, Številka:
2
Journal Article
Recenzirano
Odprti dostop
Double neutron star (DNS) systems represent extreme physical objects and the endpoint of an exotic journey of stellar evolution and binary interactions. Large numbers of DNS systems and their mergers ...are anticipated to be discovered using the Square Kilometre Array searching for radio pulsars, and the high-frequency gravitational wave detectors (LIGO/VIRGO), respectively. Here we discuss all key properties of DNS systems, as well as selection effects, and combine the latest observational data with new theoretical progress on various physical processes with the aim of advancing our knowledge on their formation. We examine key interactions of their progenitor systems and evaluate their accretion history during the high-mass X-ray binary stage, the common envelope phase, and the subsequent Case BB mass transfer, and argue that the first-formed NSs have accreted at most . We investigate DNS masses, spins, and velocities, and in particular correlations between spin period, orbital period, and eccentricity. Numerous Monte Carlo simulations of the second supernova (SN) events are performed to extrapolate pre-SN stellar properties and probe the explosions. All known close-orbit DNS systems are consistent with ultra-stripped exploding stars. Although their resulting NS kicks are often small, we demonstrate a large spread in kick magnitudes that may, in general, depend on the past interaction history of the exploding star and thus correlate with the NS mass. We analyze and discuss NS kick directions based on our SN simulations. Finally, we discuss the terminal evolution of close-orbit DNS systems until they merge and possibly produce a short γ-ray burst.
The mass function for black holes and neutron stars at birth is explored for mass-losing helium stars. These should resemble, more closely than similar studies of single hydrogen-rich stars, the ...results of evolution in close binary systems. The effects of varying the mass-loss rate and metallicity are calculated using a simple semi-analytic approach to stellar evolution that is tuned to reproduce detailed numerical calculations. Though the total fraction of black holes made in stellar collapse events varies considerably with metallicity, mass-loss rate, and mass cutoff, from 5% to 30%, the shapes of their birth functions are very similar for all reasonable variations in these quantities. Median neutron star masses are in the range 1.32-1.37 regardless of metallicity. The median black hole mass for solar metallicity is typically 8-9 if only initial helium cores below 40 (ZAMS mass less than 80 ) are counted, and 9-13 , in most cases, if helium cores with initial masses up to 150 (ZAMS mass less than 300 ) contribute. As long as the mass-loss rate as a function of mass exhibits no strong nonlinearities, the black hole birth function from 15 to 35 has a slope that depends mostly on the initial mass function for main-sequence stars. These findings imply the possibility of constraining the initial mass function and the properties of mass loss in close binaries using ongoing measurements of gravitational-wave radiation. The expected rotation rates of the black holes are briefly discussed.
ABSTRACT Nucleosynthesis, light curves, explosion energies, and remnant masses are calculated for a grid of supernovae (SNe) resulting from massive stars with solar metallicity and masses from 9.0 to ...120 M ⊙ . The full evolution is followed using an adaptive reaction network of up to 2000 nuclei. A novel aspect of the survey is the use of a one-dimensional neutrino transport model for the explosion. This explosion model has been calibrated to give the observed energy for SN 1987A, using five standard progenitors, and for the Crab SN using a 9.6 M ⊙ progenitor. As a result of using a calibrated central engine, the final kinetic energy of the SN is variable and sensitive to the structure of each pre-SN star. Many progenitors with extended core structures do not explode, but become black holes (BHs), and the masses of exploding stars do not form a simply connected set. The resulting nucleosynthesis agrees reasonably well with the Sun provided that a reasonable contribution from SNe Ia is also allowed, but with a deficiency of light s-process isotopes. The resulting neutron star initial mass function has a mean gravitational mass near 1.4 M ⊙ . The average BH mass is about 9 M ⊙ if only the helium core implodes, and 14 M ⊙ if the entire pre-SN star collapses. Only ∼10% of SNe come from stars over 20 M ⊙ , and some of these are Type Ib or Ic. Some useful systematics of Type IIp light curves are explored.
We investigate systematically the dynamical mass ejection, r-process nucleosynthesis, and properties of electromagnetic counterparts of neutron-star (NS) mergers in dependence on the uncertain ...properties of the nuclear equation of state (EOS) by employing 40 representative, microphysical high-density EOSs in relativistic, hydrodynamical simulations. The crucial parameter determining the ejecta mass is the radius R sub(1.35) of a 1.35 M sub(middot in circle) NS. NSs with smaller R sub(1.35) ("soft" EOS) eject systematically higher masses. The r-process nucleosynthesis exhibits a remarkable robustness independent of the EOS, producing a nearly solar abundance pattern above mass number 130. By the r-process content of the Galaxy and the average production per event the Galactic merger rate is limited to 4 x 10 super(-5) yr super(-1) (4 x 10 super(-4) yr super(-1)) for a soft (stiff) NS EOS, if NS mergers are the main source of heavy super(r)-nuclei.
Light curves, explosion energies, and remnant masses are calculated for a grid of supernovae resulting from massive helium stars that have been evolved including mass loss. These presupernova stars ...should approximate the results of binary evolution for stars in interacting systems that lose their envelopes close to the time of helium core ignition. Initial helium star masses are in the range 2.5-40 M , which corresponds to main-sequence masses of about 13-90 M . Common SNe Ib and Ic result from stars whose final masses are approximately 2.5-5.6 M . For heavier stars, a large fraction of collapses lead to black holes, though there is an island of explodability for presupernova masses near 10 M . The median neutron star mass in binaries is 1.35-1.38 M , and the median black hole mass is between 9 and 11 M . Even though black holes less massive than 5 M are rare, they are predicted down to the maximum neutron star mass. There is no empty "gap," only a less populated mass range. For standard assumptions regarding the explosions and nucleosynthesis, the models predict light curves that are fainter than the brighter common SNe Ib and Ic. Even with a very liberal but physically plausible increase in 56Ni production, the highest-energy models are fainter than 1042.6 erg s−1 at peak, and very few approach that limit. The median peak luminosity ranges from 1042.0 to 1042.3 erg s−1. Possible alternatives to the standard neutrino-powered and radioactive-illuminated models are explored. Magnetars are a promising alternative. Several other unusual varieties of SNe I at both high and low mass are explored.
We demonstrate by a large set of merger simulations for symmetric binary neutron stars (NSs) that there is a tight correlation between the frequency peak of the postmerger gravitational-wave (GW) ...emission and the physical properties of the nuclear equation of state (EoS), e.g., expressed by the radius of the maximum-mass Tolman-Oppenheimer-Volkhoff configuration. Therefore, a single measurement of the peak frequency of the postmerger GW signal will constrain the NS EoS significantly. For optimistic merger-rate estimates a corresponding detection with Advanced LIGO is expected to happen within an operation time of roughly a year.
We present the first comprehensive study of r-process element nucleosynthesis in the ejecta of compact binary mergers (CBMs) and their relic black hole (BH)–torus systems. The evolution of the ...BH–accretion tori is simulated for seconds with a Newtonian hydrodynamics code including viscosity effects, pseudo-Newtonian gravity for rotating BHs, and an energy-dependent two-moment closure scheme for the transport of electron neutrinos and antineutrinos. The investigated cases are guided by relativistic double neutron star (NS–NS) and NS–BH merger models, producing ∼3–6 M⊙ BHs with rotation parameters of A
BH ∼ 0.8 and tori of 0.03–0.3 M⊙. Our nucleosynthesis analysis includes the dynamical (prompt) ejecta expelled during the CBM phase and the neutrino and viscously driven outflows of the relic BH–torus systems. While typically ∼20–25 per cent of the initial accretion-torus mass are lost by viscously driven outflows, neutrino-powered winds contribute at most another ∼1 per cent, but neutrino heating enhances the viscous ejecta significantly. Since BH–torus ejecta possess a wide distribution of electron fractions (0.1–0.6) and entropies, they produce heavy elements from A ∼ 80 up to the actinides, with relative contributions of A ≳ 130 nuclei being subdominant and sensitively dependent on BH and torus masses and the exact treatment of shear viscosity. The combined ejecta of CBM and BH–torus phases can reproduce the solar abundances amazingly well for A ≳ 90. Varying contributions of the torus ejecta might account for observed variations of lighter elements with 40 ≤ Z ≤ 56 relative to heavier ones, and a considerable reduction of the prompt ejecta compared to the torus ejecta, e.g. in highly asymmetric NS–BH mergers, might explain the composition of heavy-element deficient stars.
We present the new code alcar developed to model multidimensional, multienergy-group neutrino transport in the context of supernovae and neutron-star mergers. The algorithm solves the evolution ...equations of the zeroth- and first-order angular moments of the specific intensity, supplemented by an algebraic relation for the second-moment tensor to close the system. The scheme takes into account frame-dependent effects of the order
$\mathcal {O}(v/c)$
as well as the most important types of neutrino interactions. The transport scheme is significantly more efficient than a multidimensional solver of the Boltzmann equation, while it is more accurate and consistent than the flux-limited diffusion method. The finite-volume discretization of the essentially hyperbolic system of moment equations employs methods well-known from hydrodynamics. For the time integration of the potentially stiff moment equations we employ a scheme in which only the local source terms are treated implicitly, while the advection terms are kept explicit, thereby allowing for an efficient computational parallelization of the algorithm. We investigate various problem set-ups in one and two dimensions to verify the implementation and to test the quality of the algebraic closure scheme. In our most detailed test, we compare a fully dynamic, one-dimensional core-collapse simulation with two published calculations performed with well-known Boltzmann-type neutrino-hydrodynamics codes and we find very satisfactory agreement.
ABSTRACT We present the first special relativistic, axisymmetric hydrodynamic simulations of black hole-torus systems (approximating general relativistic gravity) as remnants of binary-neutron star ...(NS-NS) and neutron star-black hole (NS-BH) mergers, in which the viscously driven evolution of the accretion torus is followed with self-consistent energy-dependent neutrino transport and the interaction with the cloud of dynamical ejecta expelled during the NS-NS merging is taken into account. The modeled torus masses, BH masses and spins, and the ejecta masses, velocities, and spatial distributions are adopted from relativistic merger simulations. We find that energy deposition by neutrino annihilation can accelerate outflows with initially high Lorentz factors along polar low-density funnels, but only in mergers with extremely low baryon pollution in the polar regions. NS-BH mergers, where polar mass ejection during the merging phase is absent, provide sufficiently baryon-poor environments to enable neutrino-powered, ultrarelativistic jets with terminal Lorentz factors above 100 and considerable dynamical collimation, favoring short gamma-ray bursts (sGRBs), although their typical energies and durations might be too small to explain the majority of events. In the case of NS-NS mergers, however, neutrino emission of the accreting and viscously spreading torus is too short and too weak to yield enough energy for the outflows to break out from the surrounding ejecta shell as highly relativistic jets. We conclude that neutrino annihilation alone cannot power sGRBs from NS-NS mergers.
Muons can be created in nascent neutron stars (NSs) due to the high electron chemical potentials and the high temperatures. Because of their relatively lower abundance compared to electrons, their ...role has so far been ignored in numerical simulations of stellar core collapse and NS formation. However, the appearance of muons softens the NS equation of state, triggers faster NS contraction, and thus leads to higher luminosities and mean energies of the emitted neutrinos. This strengthens the postshock heating by neutrinos and can facilitate explosions by the neutrino-driven mechanism.